Advances in Civil Engineering

Advances in Civil Engineering / 2020 / Article

Research Article | Open Access

Volume 2020 |Article ID 9032857 | https://doi.org/10.1155/2020/9032857

Iuliia Glushakova, Qihan Liu, Yu Zhang, Guangchun Zhou, "Conjugate Cellular Automata and Neural Network Approach: Failure Load Prediction of Masonry Panels", Advances in Civil Engineering, vol. 2020, Article ID 9032857, 12 pages, 2020. https://doi.org/10.1155/2020/9032857

Conjugate Cellular Automata and Neural Network Approach: Failure Load Prediction of Masonry Panels

Academic Editor: Filippo Ubertini
Received12 Jan 2020
Revised10 Jun 2020
Accepted17 Jun 2020
Published11 Jul 2020

Abstract

The intricate interplay between the microscopic constituents and their macroscopic properties for masonry structures complicates their failure analysis modelling. A composite strategy incorporating neural network (NN) and cellular automata (CA) is developed to predict the failure load for masonry panels with and without openings subjected to lateral loadings. The discretized panels are modelled by the CA methodology using nine neighbour cells, which derive their state values from geometric parameters and opening location placement for the panels. An identification coefficient dictated by these geometric parameters and experimental data is fed together as the input training data for the NN. The NN uses a backpropagation algorithm and two hidden layers with sigmoid activation functions to predict failure loads. This method achieves greater accuracy in prediction when compared with the yield line and finite elemental analysis (FEA) methods. The results attained elucidate the feasibility of the current methodology to complement conventional approaches such as FEA to provide additional insight into the failure mechanism of masonry panels under varied loading conditions.

1. Introduction

Masonry structures have a long history spanning centuries of usage due to simplicity in their construction and durability [1]. They are continuously subjected to dynamic loadings and intermittent abusive loadings such as earthquakes and floods. Generally, structural responses, e.g., displacements, velocities, and bending, are recorded and inversely solved for dynamic loads due to technological and economic constraints in direct measurement of external dynamic loads [2]. Various models have been proposed and can broadly be grouped under frequency-domain [3, 4] or time-domain [5, 6]. Innovative techniques such as the application of stochastic simulation of Griffith flaws [7] or a combination of topology optimization method and phase-field design [8] have been proposed to improve the modelling of the load and resistance of building materials. More attention is now paid on finding ways to simplify the modelling of the dynamic behaviour of structures [911].

Inevitable degradation suffered over years of use by masonry constructions necessitates the periodic assessment for structural integrity to ensure safety, and it can serve as an essential tool to validate the design of new constructions [12]. However, even nowadays, it is a rather difficult task to find a reliable method that would encompass a variety of masonry materials with accuracy proximal to the experimental data [13]. Physical experiments on masonry structures itself are prohibitively expensive which requires significant consumption of material and time. Thus, the amount of available experimental data of masonry structures is sparse and warrants the use of various computational analysis approaches.

The discontinuous and increasingly nonlinear nature of masonry under stress, especially after crack initiation, complicates structural analysis. Broadly, studies of masonry behaviour fall into two approaches, homogeneous and heterogeneous [14]. A homogeneous approach regards masonry as a composite material, and it is used in macromodelling to study the overall structural response [15]. It simplifies calculations, but macromodelling cannot properly predict the local behaviour. The heterogeneous approach is used in micromodelling where every masonry component, unit, and mortar are modelled individually. Hence, it is more detailed and realistic but complicated and time-consuming. A homogenization approach stands as an intermediate between macro- and micromodelling that considers masonry as a heterogeneous structure divided into periodic cells. It allows the user to extract a representative element to describe the whole structure [16, 17].

Initially, constitutive models with continuum representations were utilized for engineering modelling of masonry yielding acceptable estimations. However, the inhomogeneous composition is better simulated using discontinuous micromodelling approaches, and thus such techniques were frequently used in the past two decades [1820]. Even so, collating the multitude of compositional and structural variables inside the finite/discrete element method (F/DEM) to solve equations of motions for material deformation, contact point interactions, failure load, and crack propagation remains a daunting task. Many researchers turned to the use of artificial intelligence (AI) techniques to overcome these issues. Among such techniques, neural networks (NNs) have progressively gained popularity due to its ability to solve a wide variety of problems at lower computational costs and simplified approach. However, only a few studies incorporate the use of NNs for the approximation of masonry behaviour in general. Recently, NNs have been successfully used in civil engineering to solve a variety of problems [2130].

The analysis of panels under biaxial bending is amongst the earliest applications of NNs for the prediction of masonry behaviour [31], which showed the ability of NNs to solve complex nonlinear problems. A multilayer perceptron NN [32] and a radial basis function NN [33] were able to predict creep deformations in masonry structures showing a relatively small prediction error. Garzón-Roca et al. used NN and fuzzy logic to estimate the axial load and compressive strength of masonry made of clay bricks and cement mortar [34, 35]. Asteris and Plevris employed a NN to approximate failure for masonry under biaxial stress. They proposed a computational procedure for approximating failure curves in 3D achieving improved prediction performance and providing valuable information about the influence of different loading angles [36, 37]. Cascardi et al. proposed an artificial NN model to predict the in-plane shear strength of masonry panels strengthened by Fibre Reinforced Polymer systems based on a large experimental database [38, 39].

Another (AI) technique that has gained popularity in resolving problems of masonry structures analysis is cellular automata (CA). CA method for masonry panels is first mentioned by Zhou G. C. in 2002. They applied CA techniques to improve the standard finite element (FE) method used to calculate failure load, proposing a concept of similar zones and strength/stiffness corrector [40]. Zhang et al. applied NN to predict the cracking patterns of masonry walls loaded vertically at different orientations using the CA model and experimental data of the recorded cracking patterns [41]. By harnessing the advantages of CA and finite element analysis (FEA), Huang et al. developed a method for predicting the failure load of masonry wall panels relying on generalized strain-energy density (GSED) extracted from the so-called “base” panel [42]. And later on, Huang et al. used the combination of CA and GSED to map cracking patterns of laterally loaded masonry wall panels with openings based on displacements of CA cells calculated from the FE method and maximum correlation coefficient [43]. In our previous work, a combined NN and CA approach also showed the capacity to effectively predict the crack propagation of masonry panels with openings based on panels’ configuration information [44].

This study is focused on the homogenization approach through the investigation of the synergistic effect of NN and a modified CA method to predict failure load of masonry panels, with and without openings, subjected to lateral loading relying on panels’ geometric configuration. An opening-centric modified CA method is combined with a NN that uses a backpropagation algorithm and two hidden layers for modelling. To the best of our knowledge, this is the first report on a combined CA and NN approach for predicting the failure load of masonry panels with openings. The proposed method is verified through comparison with FEA, yield line theory (YL), and GSED methods.

2. Materials and Methods

As the geometry variation has a significant effect on the unreinforced masonry resistance to lateral loading [45], masonry panels were modelled as a structure of nine cells defined by the length-height ratio of the cell relative to the openings using the CA theory. Then, the NN is trained on the obtained CA information and wallette strength (WS) parameter to predict the failure load for the panels.

2.1. Experimental Data

This study is built on the experimental data of 55 masonry wall panels [46]. Of which, 44 are single leaf solid panels and 11 have an opening (Table 1).


PanelLength (m)Height (m)Thickness (m)Wallet strength (N/mm2)PanelLength (m)Height (m)Thickness (m)Wallet strength (N/mm2)

11202.71.30.10252.4411505.52.60.2151.78
11352.72.60.10252.4411535.53.60.2151.78
11162.73.60.10252.4412375.54.50.2151.78
11263.62.70.10252.4411482.72.60.2151.15
11902.75.20.10252.4412112.73.60.2151.15
11084.51.30.10252.4412462.74.50.2151.15
11094.53.60.10252.4411704.52.60.2151.15
11874.55.20.10252.4411784.53.60.2151.15
10945.51.30.10252.4412244.54.50.2151.15
11105.51.30.10252.4411495.52.60.2151.15
10955.53.60.10252.4411625.53.60.2151.15
11715.55.20.10252.4412315.54.50.2151.15
11212.71.30.10251.19ART014.7152.4650.10253.2
11232.72.60.10251.19SB015.6152.4750.0651.7
11172.73.60.10251.19SB055.6152.4750.1021.7
12032.75.20.10251.19SB065.6152.4750.0651.7
11074.51.30.10251.19Panels with an opening
11114.53.60.10251.19Panel 14.842.4750.10251.2
12014.55.20.10251.19Panel 24.842.4750.10251.2
10965.51.30.10251.19Panel 34.842.4750.10251.2
10975.53.60.10251.19ART024.7152.4650.10253.2
11575.55.20.10251.19ART034.7152.4650.10253.2
11722.72.60.2151.78ART044.7152.4650.10253.2
11922.73.60.2151.78ART064.7152.4650.10253.2
12612.74.50.2151.78SB025.6152.4750.0651.7
11694.52.60.2151.78SB035.6152.4750.0651.7
11734.53.60.2151.78SB045.6152.4750.0651.7
12444.54.50.2151.78SB075.6152.4750.0651.7

The experimental data was randomly divided into three groups: training, testing, and validation by 70, 15, and 15 percent, respectively.

It should be noted that while the amount of the available training data was limited, it is comparable to experimental training data sets used by other researchers in the field [38, 47], which proved to be enough to showcase the method for the current problem.

2.2. Cellular Automata

CA can be described as a model of a spatially extended decentralized system made up of a number of individual cells. Each cell is in a specific state which changes over time depending on the state of its local neighbours. CA is a collection of cells that each adapts one of a finite number of states. Single cells change in states by following a local rule that depends on the environment of the cell [48].

Herein, a new generalized CA model was developed using the concept of eight-neighbourhood Moore CA. In this opening-centric model, the opening of a panel is imparted with the role of the central cell surrounded by discretized neighbouring cells (Figure 1).

In this method, the CA model for masonry panels is described as a pattern of 3 × 3 cells. State values for cells on that model are formed to consider the size effect and the effect of length-height ratio (Figure 2). It describes the percentage of the solid area of the panel that the cell occupies. The dimensions and size of a cell on the CA model relate to the position of an opening. Hence, for the state values for every cell of a solid masonry, the panel is considered equal as shown in equation (1) (Figure 2(a)). But state values for cells of masonry panels with an opening relate to dimensions of a cell and are calculated by equation (2). Openings and areas outside the CA model have state values equal to zero, such as in Figures 2(b) and 2(c) and in Figure 2(d).

where and are the panel’s length (m) and height (m), respectively, and and are the relative length (m) and height (m) of the CA cell , respectively.

The size effect for masonry panels is taken into consideration in this way. The information on CA model parameters for masonry panels with openings is given in Table 2.


Panel

Panel 11.8531.1341.8530.91.050.525
Panel 21.6151.611.6150.91.050.525
Panel 31.9071.0251.9072.050.450
ART021.5651.5851.5651.041.060.365
ART031.791.1351.790.591.510.365
ART041.1152.4851.1151.4150.6850.365
ART062.4651.5850.6651.041.060.365
SB021.6772.261.6780.91.1250.45
SB031.342.9351.341.50.5250.45
SB042.3520.912.3532.0250.450
SB0710.910.90.90.65

2.3. Neural Network Application

Artificial neural networks are often described as being a simplified model of a human brain. They consist of neuron layers where different neurons of different layers are connected through applied synaptic weights [49]. The structure of a simple neuron for a backpropagation NN is shown in Figure 3. There the input is transmitted through a connection where it multiplies its strength by the scalar weight . The bias of the —neuron is added to the weighted input by the summing function. Then, the activation (transfer) function φ uses this sum as its argument. Backpropagation NN usually uses a sigmoid function as the activation function.

The weights and biases of NN are adjustable scalar parameters. The central principle of NN is that they can reach some desired behaviour by adjusting weights and biases. And, the network can be trained to solve a particular problem. This ability to solve complicated nonlinear problems incorporating multiple parameters and variables makes NN well-suited for the highly anisotropic properties displayed by masonry structures. In the backpropagation training algorithm, the output values calculated by the activation function are compared with the desired output or acceptable convergence. This forms a predefined error function, and the error is then fed back through the network [50]. The algorithm adjusts the weights to reduce the value of the error function in response to this information [51]. This training process repeats until the network reaches a state where the error of the calculations is minimized.

The convergence process is expected to be relatively slow due to the nature of the stochastic gradient descent algorithm used in backpropagation [52]. Diverse techniques were proposed to improve this situation (solve this problem) [5355]. However, the most popular solutions are based on “minibatches,” where the network’s training set is divided into small subsets, and each one of them is iterated through an epoch [56]. Generally, the state-of-the-art software packages and implementations are optimized to use these options [57]. However, the choice of the training method is constrained by the nature of the problem that needs to be solved with available computational resources [58]. Nonetheless, the difference in speed for small-medium size problems, such as the case of this research, is often marginal and can be ignored with the development of computing technology.

The initialization of weights is a crucial step for backpropagation. The size of the gradient changes relative to the size of the weights. If weights are too small, the gradient becomes too small as well (vanishing gradient problem) and weights can never reach the optimal global minima. This significantly slows down the learning if not prevents the network from learning at all. On the other hand, if the weights are too large, then it leads to an exploding gradient problem when the network keeps learning on the large weights and gradient keeps getting larger, never reaching the convergence point. But if all weights are initialized to zeros, the neurons learn the same features during training and evolve symmetrically unable to learn different features [59]. There are diverse methods to set the initial weights, such as Xavier initialization and He initialization, that change the randomly calculated weights depending on the activation function [6062]. However, no determined rule would work for all problems. But, stochastic initialization of weights that follows the standard normal distribution is the most used approach. Therefore, both vanishing gradients and exploding gradients are rarely a problem for large networks and modern backpropagation techniques [63].

The NN in this study is trained on two nonlinear layers using feed-forward backpropagation and the efficient Levenberg–Marquardt training algorithm [6466]. The activation functions for the first hidden layer is the hyperbolic tangent sigmoid function, which returns a matrix of elements in the interval , shown in equation (3). The second hidden layer uses the logarithmic-sigmoid activation function (equation (4)). The result’s range of this function varies from 0 to 1. And, the linear activation function is used in the output layer. The number of learning cycles was 1,000:

The most efficient configuration in time/accuracy comparison for the problem was achieved when the number of neurons in the first hidden layer is 11 and neurons in the second hidden layer are 5 (Figure 4).

2.4. Input, Training, and Output Data

Eleven parameters are used as the training data for the NN. Firstly, the panel's configuration coefficient is considered to describe every calculated CA cell as the part of a masonry panel with a given geometry:where is the panel’s configuration coefficient, and are the panel’s length (m) and height (m), respectively, and is the panel’s thickness (m).

Next, the input parameters are wallet strength () and state values from equations (1)–(2). Both von Neumann and Moore neighbourhood discretization strategies were evaluated, and the latter was found to be more accurate for the given data set and all results presented hereafter are based on it.

The backpropagation algorithm requires the input data to be less than 1 which is suitable for the sigmoid activation function. There are different methods to normalize the input data for backpropagation neural networks [6769].

It was decided to normalize parameter through multiplication on using several optimization trial experiments.

The output data in the proposed method are vectors of failure load. Those vectors were normalized to match the interval using the min-max normalisation method:where is the output failure load, is the experimental failure load, and and are minimum and maximum failure load equal to 0 and 100, respectively.

The topology of NN is shown in Figure 4.

Table 3 contains the example of input data for panels with an opening SB02, Panel 3, and ART06.


InputPanel
SB02Panel 3ART06

P0.147460.198400.09456
WS0.170000.120000.32000
S 1,10.108610.323150.45734
S 1,20.146360.000000.29407
S 1,30.108670.323150.12338
S 2,10.135760.070940.46614
S 2,20.000000.038130.00000
S 2,30.135840.070940.12575
S 3,10.054300.000000.16051
S 3,20.073180.000000.10321
S 3,30.054330.000000.04330

3. Results and Discussion

Initially, a thousand training epochs were set for the NN to validate the rationality of the devised approach. However, the training was stopped after the best performance was reached at 363 epochs. The performance, in this case, was calculated by the mean square error for the training output compared with the experimental data.

The results of the training process achieved performance and training regression and are shown in Figure 5.

The predicted results of the NN were compared with the FEA and yield line theory (YL) to evaluate the performance of the proposed method. The YL prediction was built according to the British Standard Institution BS 5628 by Chong V. L. The FEA prediction was also made by Chong V.L. using 8 × 8 FE mesh. The comparison between predictions made by different methods and experimental failure loads are listed in Table 4. Moreover, the accuracy of the different methods can be seen in Table 5.


PanelExpt. (kN/m2)YL (kN/m2)FEA (kN/m2)NN (kN/m2)PanelExpt. (kN/m2)YL (kN/m2)FEA (kN/m2)NN (kN/m2)

112010.913.169.6510.975115010.711.058.8310.727
11356.59.876.853.71011536.79.067.045.973
11167.98.886.676.55112376.58.246.26.406
11266.48.886.676.551114815.822.0915.3716.057
11906.48.466.656.480121112.820.2514.912.900
11085.26.664.875.356124611.819.2915.0912.103
11093.23.82.773.07111708.69.967.538.938
11872.83.382.335.30311785.28.765.934.087
10943.65.044.084.24612244.88.115.497.872
111055.044.084.24611496.67.75.635.864
10952.92.771.953.31711624.36.374.54.504
11712.22.411.673.59312313.95.863.926.272
11215.86.884.695.424ART014.44.564.290
11234.45.353.041.887SB012.83.1642.4642.968
11174.34.983.024.030SB052.73.1592.4573.184
12033.44.742.972.580Panels with an opening
110753.432.315.530SB067.58.76.755.920
11114.32.061.172.528Panel 11.51.311.769
12012.11.871.013.044Panel 21.21.221.503
10962.81.651.892.227Panel 321.242.125
10971.81.50.913.067ART022.752.392.487
11571.41.350.731.840ART032.62.62.933
117220.931.3523.7520.372ART042.52.212.557
119218.828.0823.618.862ART062.52.692.437
12611527.2623.5915.238SB022.42.5921.8242.461
11699.514.411.6610.059SB032.32.4381.5872.478
1173812.319.495.274SB042.22.5741.8042.602
1244811.288.638.112SB075.56.9855.175.514


PanelYL error (%)FEA error (%)NN error (%)PanelYL error (%)FEA error (%)NN error (%)

112020.7311.470.6911503.2717.480.26
113551.855.3842.92115335.225.0710.86
111612.4115.5717.08123726.774.621.44
112638.754.222.36114839.812.721.63
119032.193.911.26121158.2016.410.78
110828.086.353.00124663.4727.882.57
110918.7513.444.03117015.8112.443.93
118720.7116.7989.39117868.4614.0421.40
109440.0013.3317.93122468.9614.3864.00
11100.8018.4015.09114916.6714.7011.15
10954.4832.7614.37116248.144.654.74
11719.5524.0963.32123150.260.5160.82
112118.6219.146.48ART013.642.50
112321.5930.9157.11SB011312.006.01
111715.8129.776.27SB05179.0017.93
120339.4112.6524.12SB061610.0021.06
110731.4053.8010.59Panels with an opening
111152.0972.7941.22Panel 112.6717.95
120110.9551.9044.98Panel 21.6725.29
109641.0732.5020.47Panel 338.006.26
109716.6749.4470.38ART0213.099.57
11573.5747.8631.44ART030.0012.79
117250.0013.642.53ART0411.602.29
119249.3625.530.33ART067.602.53
126181.7357.271.59SB02824.002.55
116951.5822.745.89SB03631.007.74
117353.8818.6334.07SB041718.0018.26
124441.007.881.40SB07276.000.25

As the micromacro scale variation in masonry structure from one panel to the next makes it a tremendously difficult task to model it analytically, there is a large discrepancy in NN, FEA, and YL predictions in comparison to experimental data. In this research, two comparative models that fundamentally approach analytical load modelling from different points of scale were used.

YL theory briefly approaches the panel (slab) load prediction from a macroscale level by assuming that the panel behaves like a perfectly plastic structure and develops positive and negative yield lines under an applied overload. It also permits the determination of the ultimate load for a defined collapse mechanism. While the results from the YL method provide reasonable experimental agreement, they tend to overestimate the failure strength. Another drawback of this method is its inability to precisely define the position of openings [46].

In comparison, FEA approaches the load prediction problem from the microscale level by allowing the division of panel into many discretized units (elements) and tuning their properties to illustrate the local behaviour. FEA requires the definition of each element’s material properties, loading, geometry, location, and relationship with neighbouring elements. However, for heterogeneous structures such as masonry panels, such exhaustive elemental definitions can lead to an exponential increase in the required computing power while offering little flexibility for adapting to different panels.

In this work, we have tried to bridge the gap between FEA and YL modelling techniques by combining the cellular automata discretization and the adaptable learning ability of ANNs. CA allowed the accurate definition of panel structures and opening locations, while the macroscopic panel properties such as wallet strength, length, height, and thickness, are fed as input training data for the NN.

Further evaluation of the method was performed by using different statistical parameters such as the root-mean-squared error (RMSE), the mean absolute percentage error (MAPE), the coefficient of determination (), and the integral absolute error (IAE). In theory, the closer the values of RMSE, MAPE, and IAE parameters to zero, the higher the accuracy of the proposed NN. Besides, the closer the values to 1, the greater are the similarities between predicted and experimental failure loads. Parameters RMSE, MAPE, , and IAE were calculated by equations (7)–(10). Table 6 shows the comparison between FEA, proposed NN, and YL methods based on these statistical parameters. The comparison is calculated for 47 panels for which the YL results are given in Tables 4 and 5 and for all 55 panels studied here for FEA and NN:where is the number of tested panels, is the experimental failure load, is the predicted failure load, and is the average experimental failure load.


MethodRMSEMAPEIAE (%)
IIIIIIIIIIII

YL3.701410.309800.2978037.44
FEA1.888561.750440.203620.190050.817190.8367318.8618.33
NN1.143361.060200.188860.175790.9330.9401112.0711.81

Note: I, calculated for 47 panels with YL; II, calculated for all 55 panels without YL.

It was observed that for more than 75% of the panels the percentage error was smaller than 20%. Especially for the panels with an opening where all panels had an error deviation below 20%. The low accuracy of the prediction for a few panels can be explained by the fact that prediction in this case significantly depends on the WS parameter, limited training data, and possible variations in the experimental process. Overall, it is safe to assume that the optimized NN was able to outperform FEA and YL.

Figure 6 compares the results of NN and FEA prediction to the experimental data. While both methods can offer comparable accuracy for panels with low experimental failure loads, the error variation range for masonry panels with failure loads more than 10 kN/m2 is lower and significantly more accurate using the NN prediction.

In addition, the proposed method is compared with the GSED method proposed by Huang and based on CA and FEA [42]. Huang et al. used the strain-energy density of a known (base) panel and a new (to be predicted) panel and a criterion for matching zone similarities to calculate failure load for masonry panels.

Although the results from GSED and NN give a good prediction for failure load, a closer inspection reveals that in all cases except one, the percentage error is smaller for NN as compared with GSED results (Table 7).


PanelFailure load (kN/m2)Error (%)

Expt.GSED (1173 based)NNGSED (1173 based)NN
117220.920.0120.3724.262.53
11699.512.5310.05931.895.89
124486.428.11219.751.40
115010.710.3610.7272.90.26
11536.77.035.9734.9310.86
12376.55.36.40618.461.44

Expt.GSED (1172 based)NNGSED (1172 based)NN
11699.513.610.05943.165.89
124487.138.11210.881.40
115010.711.2310.7274.950.26
11536.77.735.97315.3710.86
12376.55.96.4069.231.44

An additional drawback of the GSED method is that its accuracy depends on the base panel used for prediction. Hence, it is crucial to have experimental information for a panel with a similar condition to the predicted panel. Besides, the NN prediction can be calculated regardless of the configuration of the new (predicted) panel. NN can be trained on available experimental data to calculate failure load for panels with any configuration.

4. Conclusions

In this study, the backpropagation NN with two hidden layers is developed to accurately predict the failure load of masonry panels subjected to lateral loading based on panels’ configuration and WS. An additional input parameter, panels’ configuration coefficient , calculated through length, height, and thickness of the panel, is introduced. This parameter together with WS and discretized cell locations is used as the input training data for the neural network. The modified CA technique is used to realize the state values of the cells where the central cell representing the opening in the panel is fixed. The most efficient configuration for the algorithm programmed for the NN was achieved when using the tangent-sigmoid activation function for the first hidden layer and logarithmic-sigmoid activation function for the second hidden layer with 11 and 5 neurons in the layers, respectively. It is shown that NNs can predict failure load for masonry panels based only on the panel’s configuration information. But, the accuracy and comprehensiveness of the network deeply depend on the training data. The results were compared with other well-established methods, and judging by RMSE, MAPE, , and IAE statistical parameters NN prediction excels the prediction accuracy of FEA and YL analysis by every parameter. While GSED and NN offered effective failure load prediction capability, the NN is found to be more accurate and general in its applicability.

Even though the proposed method has only been used to validate the experimental data used in this research, an increase in the training data for NN should allow for further refinement of its achieved results. Anyway, the proposed method can provide an additional tool for predicting failure loads and complement other methods like FEA and GSED for in-depth structural analysis. An inherent advantage of the neural modelling method over traditional numerical techniques is that it implicitly identifies and extracts the different characteristics of the panels, including their nonlinear material properties, without requiring explicit and rigorous mathematical expressions, enabling a simplification in numerical modelling and computational cost of masonry.

Data Availability

The code for neural network and relative weights and biases used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The authors would like to thank the Chinese Scholarship Council (CSC) for the PhD scholarship.

References

  1. B. Haseltine, “The evolution of the design and construction of masonry buildings in the UK,” Gestão & Tecnologia de Projetos, vol. 7, no. 2, pp. 20–26, 2012. View at: Publisher Site | Google Scholar
  2. M. Xu, J. Huang, C. Wang, and Y. Li, “Fuzzy identification of dynamic loads in presence of structural epistemic uncertainties,” Computer Methods in Applied Mechanics and Engineering, vol. 360, Article ID 112718, 2020. View at: Google Scholar
  3. Y. Liu and W. S. Shepard, “Dynamic force identification based on enhanced least squares and total least-squares schemes in the frequency domain,” Journal of Sound and Vibration, vol. 282, no. 1-2, pp. 37–60, 2005. View at: Publisher Site | Google Scholar
  4. A. N. Thite and D. J. Thompson, “The quantification of structure-borne transmission paths by inverse methods. Part 2: use of regularization techniques,” Journal of Sound and Vibration, vol. 264, no. 2, pp. 433–451, 2003. View at: Publisher Site | Google Scholar
  5. G. R. Liu, W. B. Ma, and X. Han, “An inverse procedure for identification of loads on composite laminates,” Composites Part B: Engineering, vol. 33, no. 6, pp. 425–432, 2002. View at: Publisher Site | Google Scholar
  6. F. E. Gunawan, H. Homma, and Y. Kanto, “Two-step B-splines regularization method for solving an ill-posed problem of impact-force reconstruction,” Journal of Sound and Vibration, vol. 297, no. 1-2, pp. 200–214, 2006. View at: Publisher Site | Google Scholar
  7. M. Pathirana, N. Lam, S. Perera, L. Zhang, D. Ruan, and E. Gad, “Risks of failure of annealed glass panels subject to point contact actions,” International Journal of Solids and Structures, vol. 129, pp. 177–194, 2017. View at: Publisher Site | Google Scholar
  8. S. Jeong, H. K. Seong, C. W. Kim, and J. Yoo, “Structural design considering the uncertainty of load positions using the phase field design method,” Finite Elements in Analysis and Design, vol. 161, pp. 1–15, 2019. View at: Publisher Site | Google Scholar
  9. Y. Yang, N. T. K. Lam, and L. Zhang, “Evaluation of simplified methods of estimating beam responses to impact,” International Journal of Structural Stability and Dynamics, vol. 12, no. 3, Article ID 1250016, 2012. View at: Publisher Site | Google Scholar
  10. J. Sun, N. Lam, L. Zhang, D. Ruan, and E. Gad, “Contact forces generated by hailstone impact,” International Journal of Impact Engineering, vol. 84, pp. 145–158, 2015. View at: Publisher Site | Google Scholar
  11. J. Mazars and S. Grange, “Simplified strategies based on damage mechanics for concrete under dynamic loading,” Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, vol. 375, no. 2085, Article ID 20160170, 2017. View at: Publisher Site | Google Scholar
  12. P. Roca, M. Cervera, G. Gariup, and L. Pela’, “Structural analysis of masonry historical constructions. Classical and advanced approaches,” Archives of Computational Methods in Engineering, vol. 17, no. 3, pp. 299–325, 2010. View at: Publisher Site | Google Scholar
  13. F. Ferretti, B. Ferracuti, C. Mazzotti, and M. Savoia, “Destructive and minor destructive tests on masonry buildings: experimental results and comparison between shear failure criteria,” Construction and Building Materials, vol. 199, pp. 12–29, 2019. View at: Publisher Site | Google Scholar
  14. S. Komurcu and A. Gedikli, “Macro and micro modeling of the unreinforced masonry shear walls,” European Journal of Engineering and Natural Sciences, vol. 3, pp. 116–123, 2019. View at: Google Scholar
  15. H. O. Koksal, B. Doran, A. O. Kuruscu, and A. Kocak, “Elastoplastic Finite Element analysis of masonry shear walls,” KSCE Journal of Civil Engineering, vol. 20, no. 2, pp. 784–791, 2016. View at: Google Scholar
  16. S. Casolo and G. Milani, “A simplified homogenization-discrete element model for the non-linear static analysis of masonry walls out-of-plane loaded,” Engineering Structures, vol. 32, no. 8, pp. 2352–2366, 2010. View at: Publisher Site | Google Scholar
  17. G. Milani, “Simple homogenization model for the non-linear analysis of in-plane loaded masonry walls,” Computers & Structures, vol. 89, no. 17-18, pp. 1586–1601, 2011. View at: Google Scholar
  18. F. Portioli, C. Casapulla, M. Gilbert, and L. Cascini, “Limit analysis of 3D masonry block structures with non-associative frictional joints using cone programming,” Computers & Structures, vol. 143, pp. 108–121, 2014. View at: Google Scholar
  19. S. Calderón, “Cristián Sandoval and Oriol Arnau, “Shear response of partially-grouted reinforced masonry walls with a central opening: testing and detailed micro-modelling,” Materials & Design, vol. 118, pp. 122–137, 2017. View at: Google Scholar
  20. V. Koç and Y. Emiroğlu, “Examining the behaviour of a reinforced concrete system under horizontal loads and determination of preliminary design reinforcement ratios,” International Journal of Civil Engineering, vol. 15, no. 7, pp. 1035–1052, 2017. View at: Publisher Site | Google Scholar
  21. M. H. Baziar and Y. Jafarian, “Assessment of liquefaction triggering using strain energy concept and ANN model: capacity energy,” Soil Dynamics and Earthquake Engineering, vol. 27, no. 12, pp. 1056–1072, 2007. View at: Publisher Site | Google Scholar
  22. A. Panakkat and H. Adeli, “Recurrent neural network for approximate earthquake time and location prediction using multiple seismicity indicators,” Computer-Aided Civil and Infrastructure Engineering, vol. 24, no. 4, pp. 280–292, 2009. View at: Publisher Site | Google Scholar
  23. H. Erdem, “Prediction of the moment capacity of reinforced concrete slabs in fire using artificial neural networks,” Advances in Engineering Software, vol. 41, no. 2, pp. 270–276, 2010. View at: Publisher Site | Google Scholar
  24. T. Nishikawa, J. Yoshida, T. Sugiyama, and Y. Fujino, “Concrete crack detection by multiple sequential image filtering,” Computer-Aided Civil and Infrastructure Engineering, vol. 27, no. 1, pp. 29–47, 2011. View at: Google Scholar
  25. A. Joghataie and M. S. Dizaji, “Designing high-precision fast nonlinear dam neuro-modelers and comparison with finite-element analysis,” Journal of Engineering Mechanics, vol. 139, no. 10, pp. 1311–1324, 2013. View at: Publisher Site | Google Scholar
  26. Z. H. Duan, S. C. Kou, and C. S. Poon, “Using artificial neural networks for predicting the elastic modulus of recycled aggregate concrete,” Construction and Building Materials, vol. 44, pp. 524–532, 2013. View at: Publisher Site | Google Scholar
  27. M. Gul, F. N. Catbas, and H. Hattori, “Image-based monitoring of open gears of movable bridges for condition assessment and maintenance decision making,” Journal of Computing in Civil Engineering, vol. 29, no. 2, Article ID 4014034, 2015. View at: Google Scholar
  28. M. Nikoo, F. Torabian Moghadam, and Ł. Sadowski, “Prediction of concrete compressive strength by evolutionary artificial neural networks,” Advances in Materials Science and Engineering, vol. 2015, Article ID 849126, 8 pages, 2015. View at: Publisher Site | Google Scholar
  29. P. G. Asteris, K. G. Kolovos, M. G. Douvika, and K. Roinos, “Prediction of self-compacting concrete strength using artificial neural networks,” European Journal of Environmental and Civil Engineering, vol. 20, no. sup1, pp. s102–s122, 2016. View at: Publisher Site | Google Scholar
  30. N. Wang, Q. Zhao, S. Li, X. Zhao, and P. Zhao, “Damage classification for masonry historic structures using convolutional neural networks based on still images,” Computer-Aided Civil and Infrastructure Engineering, vol. 33, no. 12, pp. 1073–1089, 2018. View at: Publisher Site | Google Scholar
  31. A. Mathew, B. Kumar, B. P. Sinha, and R. F. Pedreschi, “Analysis of masonry panel under biaxial bending using ANNs and CBR,” Journal of Computing in Civil Engineering, vol. 13, no. 3, pp. 170–177, 1999. View at: Publisher Site | Google Scholar
  32. M. M. R. Taha, A. Noureldin, N. El-Sheimy, and N. G. Shrive, “Artificial neural networks for predicting creep with an example application to structural masonry,” Canadian Journal of Civil Engineering, vol. 30, no. 3, pp. 523–532, 2003. View at: Publisher Site | Google Scholar
  33. A. El-Shafie, T. Abdelazim, and A. Noureldin, “Neural network modeling of time-dependent creep deformations in masonry structures,” Neural Computing and Applications, vol. 19, no. 4, pp. 583–594, 2010. View at: Publisher Site | Google Scholar
  34. J. Garzón-Roca, J. M. Adam, C. Sandoval, and P. Roca, “Estimation of the axial behaviour of masonry walls based on artificial neural networks,” Computers & Structures, vol. 125, pp. 145–152, 2013. View at: Google Scholar
  35. J. Garzón-Roca, C. O. Marco, and J. M. Adam, “Compressive strength of masonry made of clay bricks and cement mortar: estimation based on neural networks and fuzzy logic,” Engineering Structures, vol. 48, pp. 21–27, 2013. View at: Google Scholar
  36. P. G. Asteris and V. Plevris, “Neural network approximation of the masonry failure under biaxial compressive stress,” in Proceedings of the SEECCM III 3rd South-East European Conference on Computational Mechanics, pp. 12–14, Kos Island, Greece, June 2013. View at: Google Scholar
  37. V. Plevris and P. G. Asteris, “Modeling of masonry failure surface under biaxial compressive stress using Neural Networks,” Construction and Building Materials, vol. 55, pp. 447–461, 2014. View at: Google Scholar
  38. A. Cascardi, F. Micelli, and M. A. Aiello, “Analytical model based on artificial neural network for masonry shear walls strengthened with FRM systems,” Composites Part B: Engineering, vol. 95, pp. 252–263, 2016. View at: Publisher Site | Google Scholar
  39. A. Cascardi, F. Micelli, and M. A. Aiello, “An artificial neural networks model for the prediction of the compressive strength of FRP-confined concrete circular columns,” Engineering Structures, vol. 140, pp. 199–208, 2017. View at: Google Scholar
  40. G. C. Zhou, M. Y. Rafiq, G. Bugmann, and D. J. Easterbrook, “Cellular automata model for predicting the failure pattern of laterally loaded masonry wall panels,” Journal of Computing in Civil Engineering, vol. 20, no. 6, pp. 400–409, 2006. View at: Publisher Site | Google Scholar
  41. Y. Zhang, G. C. Zhou, Y. Xiong, and M. Y. Rafiq, “Techniques for predicting cracking pattern of masonry wallet using artificial neural networks and cellular automata,” Journal of Computing in Civil Engineering, vol. 24, no. 2, pp. 161–172, 2010. View at: Google Scholar
  42. Y. Huang, Y. Zhang, M. Zhang, and G. Zhou, “Method for predicting the failure load of masonry wall panels based on generalized strain-energy density,” Journal of Engineering Mechanics, vol. 140, no. 8, Article ID 04014061, 2014. View at: Publisher Site | Google Scholar
  43. Y. Huang, Y. Zhang, M. Zhang, C. Lv, and G. Zhou, “Cellular automata method for mapping cracking patterns of laterally loaded wall panels with openings,” Engineering Review, vol. 35, no. 1, pp. 81–88, 2015. View at: Google Scholar
  44. I. Glushakova, Y. Zhang, and G. Zhou, “Application of deep-learning neural networks for mapping cracking patterns of masonry panels with openings,” in Proceedings of the 13th Canadian Masonry Symposium, Vancouver, Canada, June 2017. View at: Google Scholar
  45. J. Park, “Investigation of the geometric variation effect on seismic performance of low-rise unreinforced masonry structures through fragility analysis,” International Journal of Civil Engineering, vol. 16, no. 1, pp. 93–106, 2016. View at: Publisher Site | Google Scholar
  46. V. L. Chong, The Behaviour of Laterally Loaded Masonry Panels with Openings, University of Plymouth, 1993.
  47. Q. Zhou, F. Zhu, X. Yang, F. Wang, B. Chi, and Z. Zhang, “Shear capacity estimation of fully grouted reinforced concrete masonry walls using neural network and adaptive neuro-fuzzy inference system models,” Construction and Building Materials, vol. 153, pp. 937–947, 2017. View at: Publisher Site | Google Scholar
  48. F. Berto and J. Tagliabue, “Cellular automata,” in The Stanford Encyclopedia of Philosophy, Fall 2017., E. N. Zalta, Ed., Metaphysics Research Lab, Stanford University, Stanford, CA, USA, 2017. View at: Google Scholar
  49. D. J. Norris, “Machine learning: deep learning,” in Beginning Artificial Intelligence with the Raspberry Pi, pp. 211–247, Apress, New York, NY, USA, 2017. View at: Google Scholar
  50. V. Rashidian and M. Hassanlourad, “Application of an artificial neural network for modeling the mechanical behavior of carbonate soils,” International Journal of Geomechanics, vol. 14, no. 1, pp. 142–150, 2014. View at: Publisher Site | Google Scholar
  51. Q. Zhou, F. Wang, and F. Zhu, “Estimation of compressive strength of hollow concrete masonry prisms using artificial neural networks and adaptive neuro-fuzzy inference systems,” Construction and Building Materials, vol. 125, pp. 417–426, October 2016. View at: Publisher Site | Google Scholar
  52. N. M. Nawi, R. S. Ransing, M. N. M. Salleh, R. Ghazali, and N. A. Hamid, “An improved back propagation neural network algorithm on classification problems,” in Database Theory and Application, pp. 177–188, Bio-Science and Bio-Technology, Berlin, Germany, 2010. View at: Google Scholar
  53. N. M. Nawi, A. Khan, and M. Z. Rehman, “A new Levenberg marquardt based back propagation algorithm trained with cuckoo search,” Procedia Technology, vol. 11, pp. 18–23, 2013. View at: Publisher Site | Google Scholar
  54. M. Otair and A. S. Walid, “Speeding up back-propagation neural networks,” Proceedings of the 2005 Informing Science and IT Education Joint Conference, vol. 1, 2005. View at: Google Scholar
  55. D. Krotov and J. J. Hopfield, “Unsupervised learning by competing hidden units,” Proceedings of the National Academy of Sciences, vol. 116, no. 16, pp. 7723–7731, 2019. View at: Google Scholar
  56. D. Masters and C. Luschi, “Revisiting small batch training for deep neural networks,” 2018, https://arxiv.org/abs/1804.07612. View at: Google Scholar
  57. I. H. Witten, E. Frank, M. A. Hall, and C. J. Pal, “Deep learning,” in Data Mining, pp. 417–466, Elsevier, Amsterdam, Netherlands, 2017. View at: Google Scholar
  58. D. V. Raman, A. P. Rotondo, and T. O’Leary, “Fundamental bounds on learning performance in neural circuits,” Proceedings of the National Academy of Sciences, vol. 116, no. 21, pp. 10537–10546, 2019. View at: Publisher Site | Google Scholar
  59. X. Glorot and Y. Bengio, “Understanding the difficulty of training deep feedforward neural networks,” Proceedings of the Thirteenth International Conference on Artificial Intelligence and Statistics, vol. 9, pp. 249–256, 2010. View at: Google Scholar
  60. K. He, X. Zhang, S. Ren, and J. Sun, “Delving deep into rectifiers: surpassing human-level performance on ImageNet classification,” 2015, https://arxiv.org/abs/1502.01852. View at: Google Scholar
  61. P. Dolezel, P. Skrabanek, and L. Gago, “Weight initialization possibilities for feedforward neural network with linear saturated activation functions,” IFAC-PapersOnLine, vol. 49, no. 25, pp. 49–54, 2016. View at: Publisher Site | Google Scholar
  62. J. Qiao, S. Li, and W. Li, “Mutual information based weight initialization method for sigmoidal feedforward neural networks,” Neurocomputing, vol. 207, pp. 676–683, 2016. View at: Publisher Site | Google Scholar
  63. Y. LeCun, Y. Bengio, and G. Hinton, “Deep learning,” Nature, vol. 521, no. 7553, pp. 436–444, 2015. View at: Publisher Site | Google Scholar
  64. B. M. Wilamowski and H. Hao Yu, “Improved computation for levenberg-marquardt training,” IEEE Transactions on Neural Networks, vol. 21, no. 6, pp. 930–937, 2010. View at: Publisher Site | Google Scholar
  65. E. M. Golafshani, A. Rahai, M. H. Sebt, and H. Akbarpour, “Prediction of bond strength of spliced steel bars in concrete using artificial neural network and fuzzy logic,” Construction and Building Materials, vol. 36, pp. 411–418, 2012. View at: Publisher Site | Google Scholar
  66. T.-H. Tran and N.-D. Hoang, “Predicting colonization growth of algae on mortar surface with artificial neural network,” Journal of Computing in Civil Engineering, vol. 30, no. 6, Article ID 04016030, 2016. View at: Publisher Site | Google Scholar
  67. G. Zhou, D. Pan, X. Xu, and Y. M. Rafiq, “Innovative ANN technique for predicting failure/cracking load of masonry wall panel under lateral load,” Journal of Computing in Civil Engineering, vol. 24, no. 4, pp. 377–387, 2010. View at: Publisher Site | Google Scholar
  68. U. Naik and S. Kute, “Span-to-depth ratio effect on shear strength of steel fiber-reinforced high-strength concrete deep beams using ANN model,” International Journal of Advanced Structural Engineering, vol. 5, no. 1, p. 29, 2013. View at: Publisher Site | Google Scholar
  69. M. AlHamaydeh, I. Choudhary, and K. Assaleh, “Virtual testing of buckling-restrained braces via nonlinear autoregressive exogenous neural networks,” Journal of Computing in Civil Engineering, vol. 27, no. 6, pp. 755–768, 2013. View at: Publisher Site | Google Scholar

Copyright © 2020 Iuliia Glushakova et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder
Views255
Downloads267
Citations

Related articles

Article of the Year Award: Outstanding research contributions of 2020, as selected by our Chief Editors. Read the winning articles.