Modelling of Structures under Seismic, Impact, and Shock Vibrations
View this Special IssueResearch Article  Open Access
Peizhen Li, Tangzhenhao Li, Zheng Lu, Jin Li, "Parametric Study on Dynamic Response of FRP Masonry Structures under the Impacts of Debris Flow", Shock and Vibration, vol. 2018, Article ID 4527571, 20 pages, 2018. https://doi.org/10.1155/2018/4527571
Parametric Study on Dynamic Response of FRP Masonry Structures under the Impacts of Debris Flow
Abstract
The aim of this study was to investigate the influences of different parameters on the performance of fiber reinforced polymer (FRP) masonry structures under debris flow using finite element models that were established using the software LSDYNA. The overall structural responses under the impacts of viscous debris flows were analyzed based on an indepth parametric study of some key factors (fiber types, relative impact positions, etc.). The results show that the diagonal and intersecting parallel types of FRP arrangements elicit better performances than horizontal types. Use of wider fiber cloths leads to the minimization of the structural response after its impact by debris flow. In addition, glass fiber reinforced polymer (GFRP) yields the best results among all studied materials in reducing local damage, while carbon fiber reinforced polymer (CFRP) yields a better overall structural response. Impact positions at the center of the wall are more unfavorable than those at the corners.
1. Introduction
Debris flow, usually caused by blizzards or heavy rains, is a solidliquid, twophase fluid, containing a large number of sand and gravel particles [1]. Mountain areas, where debris flow is generated from frequent eruptions [2], account for twothirds of China’s land area [3]. In these areas, masonry structures constitute one of the main structural forms used in inhabited areas. According to published statistical data, masonry bricks are used to occupy 70% of construction materials in China. Besides, 80% of houses in mountainous areas have been built using bricks [4]. However, the resistance of such masonry structures to geological disasters is generally quite weak, and these parts of China have not been well developed economically [5]. Thus, debris flow has posed a great threat to the economic development and personnel safety in mountain areas. As a line of defense against debris flow, structures play an important role in ensuring the safety of residents’ lives and their properties. Additionally, according to the “National Comprehensive Disaster Prevention and Reduction Plan” promoted by the Chinese government, disaster prevention and reduction is one of the most important national strategies [6]. Consequently, taking effective measures to strengthen the structures in areas threatened by debris flow is of significant importance.
In order to decrease losses caused by debris flow, many scholars all over the world have been expending efforts to study debris flow, as well as propose various reinforcing methods, for example, energy dissipation methods [7–14]. Wu [15] used a cantilever beam, or a simply supported beam, to simulate the dam and other structures and obtained an equation of impact force based on material mechanics. This method assumed that rocks and structures were elastic and ignored the plastic deformation which might have occurred under debris impact. Accordingly, this assumption resulted in an overestimation of the impact force of rocks. Based on such studies, scholars began to consider viscoelastic and elastoplastic behaviors of the material during collisions and developed other formulations and equations [16]. The contact stiffness method was based on a springmass system in which an effective contact stiffness was needed. This method was recommended in the Highway Bridge Design Specification published by the National Association of Australian State Road Authorities [17]. Some scholars established calculation methods on the debris flow force under rock impacts from the viewpoint of energy conversion. Hungr [18] analyzed the relationship between the impact energy of debris flow and the bending deformation energy of a beam, by simplifying the representation of the support scheme to a cantilever scheme. In addition, he introduced a structural stiffness coefficient and proposed a practical calculation formula for the debris flow force.
In regard to the reinforcing method, FRP plays an important role in repairing and reinforcing building structures because of its high strength, light weight, convenience in construction, and abilities of corrosion and aging resistance [19]. In a seismic test, the loadbearing capacity of a brick wall, reinforced by carbon fiber reinforced polymer (CFRP), was significantly improved in the elastic stage [20]. Seible [21] proved that the ductility of a brick wall reinforced by CFRP can be improved under earthquake actions. Based on the experiment of damaged masonry walls strengthened by CFRP under lowcyclic reversal tests, Wang [22] proved that the function of CFRP was similar to that of the rod in the truss model. CFRP delayed or prevented the development of cracks and improved the dynamic response of masonry walls.
Based on these studies, it can be concluded that FRP reinforcement can effectively increase the structural stiffness to resist seismic loads. However, the performances of conventional and FRPreinforced structures under debris flow impact have been rarely investigated [23]. In fact, in poor mountainous area, there are plenty of masonry structures which are vulnerable to debris flow loads [24]. According to the study of Guo et al. [25], FRP is favorable for masonry walls for its impact resistance. Hence, this FRP strengthening method may constitute a costeffective way to increase the structural resistance to debris flow. However, such a reinforcement method cannot be applied in construction without determining its parameters since they are key factors in influencing the structural performance. Consequently, this study mainly focuses on the dynamic response of FRPreinforced masonry structures under debris flow impacts, based on an extensive parametric study of some of key factors, including fiber reinforcement arrangement, width and types of the FRP material, and relative impact positions. Simulations are developed herein in accordance with the following steps. Select proper material models and element types based on relevant research experience. Establish structural models. To complete this study, thirteen twostory models are established based on four assumptions. Analyze structural response under dead load by dynamic relaxation method. This is a necessary step to consider the influence of the dead load and the slurry load. Investigate failure process. The parametric study is developed mainly at this step by focusing on the rock impact time history curves, damage of the wall, and the displacement of the structure.
2. Numerical Stimulation
2.1. Debris Flow Model
The impact of large rocks in debris flow is the main cause of structural damage [26]. The most popular method used to calculate the impact force of rocks in debris flow is the quasistatic method. However, this method does not take into account the time effect and complex contacts. By using the LSDYNA software, these defects can be solved [27]. LSDYNA is a generalpurpose, finite element program, which is suitable for solving nonlinear collision problems. The modeled and calculated numerical results have been shown to be in good agreement [28–30]. Wei [31] simulated clay brick unreinforced masonry walls through the use of LSDYNA software. It is pointed out that the numerical model can reasonably predict the damage levels by comparing them with field test data. According to Elsanadedy [32], FRP can be wellsimulated through LSDYNA by studying reinforced concrete oneway slabs upgraded with FRP composites.
Impact position and velocity are two core parameters of the rock in debris flow. According to Proske [33], the largest impact force of debris flow concentrates in regions corresponding to the range of 20–33% of flow depths. He [34] found that rocks concentrated at 1/2 of the maximum depth are more unfavorable to debris flow. Based on Zeng’s study [35], the typical velocity of boulders in debris flow ranges between 3.5 m/s and 5 m/s. Based on these studies, it is assumed herein that the velocity of the rock is 5 m/s, and the impact position—whose depth is equal to the height of the first floor of the structure—is at 1/2 of the depth of the debris flow. The density is selected as 2000 kg/m^{3}. Furthermore, the rock’s size is 0.5 m × 0.5 m × 0.5 m, and its density is defined as 2800 kg/m^{3}, resulting in the total mass of approximately 350 kg.
2.2. Masonry Structural Model
The masonry structural model comprises a typical twostory masonry structure designed based on Chinese codes [36] and consists of clay bricks with a thickness of 240 mm. The story height is 3.3 m and the layouts are shown in Figure 1. The compressive strengths of mortar and brick are 7.5 MPa and 15.0 MPa, respectively. The FRP is arranged at the wall of axis ③, and every piece of FRP comprises only one layer.
(a)
(b)
For FRP structures, including the column, ring beam, and masonry wall, the minimum element mesh size was 0.03 m, and the maximum element mesh size was 0.12 m. For slab and the roof, the element mesh size was 0.12 m. In addition, for steel, the minimum element mesh size was 0.1 m, and the maximum element mesh size was 0.2 m. Full 3D finite element model with the boundary conditions is shown in Figure 2. All freedom of the bottom nodes of the model is constrained.
The construction of the masonry model was based on the following assumptions.
(a) FRP behaves as an elastic material during the impact. (b) The masonry walls and FRP share the same nodes in the constructed model [37]. (c) The hydrostatic pressure is ignored in the simulation. (d) The slurry of the debris flow is assumed to behave as a onedimensional Newtonian fluid, and it impacts the wall perpendicularly.
2.3. Material Model
The anisotropic brittle damage model (MAT96) is selected to simulate masonry materials [38, 39]. This an anisotropic brittle damage model which admits progressive degradation of tensile and shear strengths across smeared cracks that are initiated under tensile loadings [40]. The continuous surface cap model (CSCM), which considers the hardening, damage, and strain rate effects of the material [41], is used to simulate concrete. According to the results of Xu and Fang [42], carbon fiber has a large elastic modulus, and the tensile properties of this material were not affected by the strain rate. Accordingly, carbon fiber could be treated as an ideal elastomeric material. The constitutive relationship of carbon fiber is modeled by the stress–strain relationship corresponding to a linearly elastic material [43]. Moreover, the stiffness of rock is much larger than that of the structure, and further considerations of this matter are beyond the main focus of this study. Therefore, the rock is assumed to be rigid.
3. Comparative Analyses of Reinforced and Unreinforced Models
Before conducting a parametric study on the dynamic response of FRPreinforced masonry structures under debris flow, it is necessary to evaluate the effectiveness of the FRPreinforced masonry structure. Hence, the corresponding comparative analyses of reinforced and unreinforced models are discussed in this section. As it is shown in Figure 3, the FRP is arranged in a diagonal arrangement in the reinforced model. The width of FRP of the model is 150 mm, and its thickness is 0.167 mm. The unreinforced model is regarded as the control model. The rock impact occurs at the center of the wall.
3.1. Dynamic Relaxation Analyses
Dynamic relaxation is a numerical method aiming to find a geometry where all forces are in equilibrium. It is a step employed to consider the influences of the dead load and the slurry load before the rock impact.
According to Tang et al.’s study [44], hydrodynamic pressure of the slurry can be calculated using the formula . The empirical coefficient depends on the type and stage of debris flow. is the density of slurry. is the density of slurry. The slurry load can be simplified as uniformly distributed load [39], which is loading on the impacted cross wall. In this article, the slurry load is considered only in stationary stage. Hence, the empirical coefficient equals to 0.5. The parameters of hydrodynamic pressure of the slurry are shown in Table 1.

In simulation process, loading curve was defined by key word DEFINECURVE. Then, the uniformly distributed load was applied on the surface elements of impacted wall by defining key word SETSEGM and LOADSEGMENTSET. Dynamic relaxation was completed by defining key word CONTROL_DYNAMIC_RELAXATION.
In this step, the maximum stress of the control model is 1.61 MPa, while the corresponding maximum stress value for the reinforced model is 1.45 MPa. For the entire model, the stress of the reinforced model is relatively smaller. Similar to the stress distribution, the displacement of the reinforced model is smaller than that of the control model. The maximum displacement of the control model is 2.52 mm, while that for the reinforced model is 1.39 mm.
3.2. Rock Impact Time History Curves
After the dynamic relaxation analysis, the rock began to move to the masonry wall at a speed of 5 m/s when = 0 s. As shown in Figure 4, the collision between the rock and the wall triggered the impulsive force.
(a)
(b)
It can be seen that the impact force reaches its maximum value within 10 ms. The peak forces for the control and reinforced models are 208 kN and 218 kN, respectively. The reason for this finding is that the stiffness of the reinforced model is larger than that of the control model, since the FRP contributes to the increase of the stiffness. Moreover, owing to the impedance of the FRP, there are two successive impulses in the reinforced model. However, a second impact on FRP is favorable since the rest of the structure can be made full use of. In other words, an increased proportion of the energy of the rock is absorbed by the reinforced model. Consequently, the damage in the reinforced house that is expected to be caused by the rock could be reduced. Figure 5 shows the condition during the impact at = 15 ms. It can be seen that the rock has crashed into the structure in the control model at this moment, while the rock is stopped by the FRP of the reinforced model. This result also indicates that the reinforced model performs better than the unreinforced one.
(a)
(b)
3.3. Displacement Response
Four typical nodes are selected to analyze the overall displacement response of the structure. The position and number of the nodes are shown in Figure 6. Figure 7 shows the displacement time history curves of these typical nodes. It can be seen that the fluctuation of the displacement of the reinforced model is much smaller than that of the control model. In addition, as it is shown in Table 2, the maximum displacement of the control model is a little larger than that of the reinforced model. This is because FRP absorbs part of the energy of the rock, and because the masonry walls of the reinforced model are fully used.

(a)
(b)
(c)
(d)
In summary, the reinforced model performs better than the unreinforced model. Therefore, a parametric study on the dynamic response of FRPreinforced masonry structures under debris flow impact is necessary.
4. Parametric Study on Dynamic Responses
4.1. Fiber Reinforcement Arrangement
4.1.1. Scheme of Reinforcement
Three finite element models are established. They are reinforced by FRP and classified in three types, namely, the diagonal type, the intersecting parallel type, and the horizontal type. The reinforcement scheme is shown in Table 3, and the corresponding finite element models are shown in Figure 8.

(a)
(b)
(c)
4.1.2. Dynamic Response
(a) Rock Impact Time History Curves. Figure 9 shows the impact time history curves for models 1 to 3. It can be seen that the peak values of the impact forces of models 1 to 3 are 209 kN, 212 kN, and 211 kN, respectively. In addition, based on the irregular pattern of the curves upon rock impact with the FRP, it can be noticed that the impact force is complex. This is because the vibration of the model caused by the impact triggers the interaction between the rock and the FRP. As it is shown in Figure 9, the arrangement of carbon fiber has influenced the magnitude of the impact force and the duration of the interaction. Among these arrangement types, the horizontal FRP arrangement reduces the force and prolongs the interaction between the rock and the fiber.
(a)
(b)
(c)
(b) Damage Response. Figure 10 shows the stress contours for models 1 to 3 at = 300 ms. It can be seen that local damage is the main failure mode. A certain degree of stress concentration occurs at the four corners of the rectangular hole and at the bottom of the wall. Maximum stress appears at the locations of the fixed supports that are connected to the ground. The maximum stress values of models 1 to 3 are 1.06 MPa, 0.94 MPa, and 1.15 MPa, respectively. Based on the stress contours, it can be deduced that the stress values of the diagonal and the intersecting parallel types are smaller, because the stress of the wall distributes more uniformly in these two types. Consequently, these two types of arrangements are more favorable for debris flow resistance. Moreover, it can also be noted that there is variation in the way local damage develops, in accordance with the arrangement type. It can be seen that the local damage develops along the direction of fiber arrangements. This damage results owing to the change of stiffness, which is affected by FRP. These bricks have a higher stiffness and tend to absorb more energy during the interaction, which results in the local damage.
(a)
(b)
(c)
(c) Displacement Response. Four typical nodes are selected to analyze the overall displacement response of the structure. The position and number of the nodes are the same as those shown in Figure 6. Figure 11 shows the displacement time history curves of these typical nodes. It can be seen that the displacement patterns of the nodes in these three models are basically the same. Model 3 yields the largest displacement amplitude, and the largest final residential strain, while model 1 yields the smallest. This is because the arrangement types of models 1 and 2 make the stress distribute uniformly on the wall. In addition, the wave of the time history curve of node 4 is denser than the others, indicating that node 4 has a higher vibration frequency. This is because node 4 is on the structural column, whose outofplan stiffness is larger than that of the masonry wall. A higher stiffness leads to a higher vibration frequency.
(a)
(b)
(c)
(d)
Table 4 lists the maximum displacement of each typical node. The results from these data are consistent with the elicited results shown in Figure 11. According to these analyses, the diagonal and the intersecting parallel types of FRP arrangements elicit better performances in resisting debris flow than the horizontal type.

4.2. Fiber Width
4.2.1. Scheme of Reinforcement
Three finite element models are established. The FRP of these models are all arranged in a diagonal arrangement, but with different widths. The scheme of reinforcement is shown in Table 5, and the corresponding finite element models are shown in Figure 12.

(a)
(b)
(c)
4.2.2. Dynamic Response
(a) Rock Impact Time History Curves. Figure 13 shows the impact time history curves from model 4 to model 6. It can be seen that the peak values of the impact forces of models 4 to 6 are 202 kN, 208 kN, and 209 kN, respectively. The differences of the force values result from the different wall stiffness values, which are affected by the width of the FRP. As it is mentioned in Section 4.1.2, irregular curve patterns result from the vibration of the model.
(a)
(b)
(c)
(b) Damage Response. Figure 14 shows the stress contours for models 4 to 6 at = 300 ms. It can be found that local damage is the main failure mode. A certain degree of stress concentration occurs at the four corners of the hole, and at the bottom of the wall. The maximum stress appears at the fixed supports connected to the ground.
(a)
(b)
(c)
By comparing the damage of the three models, model 4 yields the smallest final damage, while model 6 yields the largest. There are two main reasons for this. First, at increasing FRP widths, the duration of the interaction between the rock and FRP becomes longer. Second, the stiffness of the wall increases when the width of the carbon fiber cloth increases. According to (where is the impact force, is the mass, is the velocity, and is the stiffness), the impact force increases with increases in stiffness. Consequently, an increased energy is absorbed by the wall when the fiber cloth is wider, which results in an increased damage.
(c) Displacement Response. Figure 15 shows the displacement time history curves of selected typical nodes. From these curves, the displacement patterns of the nodes in the three models are generally the same. Model 4 has the largest displacement amplitude, and the largest final residential strain, while model 6 has the smallest. This is because the deformation of the FRP cloth contributes to energy absorption, which results in decreases of the wall displacement. In addition, the time history curve of node 4 is denser than the others, indicating that node 4 has a higher vibration frequency. As it was mentioned before, node 4 is on the structural column, whose outofplan stiffness is larger than that of the masonry wall. A higher stiffness leads to a higher vibration frequency.
(a)
(b)
(c)
(d)
Table 6 lists the maximum displacement of each typical node. The results from these data are consistent with the elicited results shown in Figure 15. In general, a wider fiber cloth minimizes the response of the structure after its impact by debris flow.

4.3. Fiber Types
To analyze the effect of different fiber types on the reinforced structure under debris flow, carbon fiber reinforced polymer (CFRP), glass fiber reinforced polymer (GFRP), and aramid fiber reinforced polymer (AFRP) were chosen. The main differences of these FRPs exist in their physical parameters, including their elastic moduli, tensile strengths, and ultimate elongation rates.
4.3.1. Scheme of Reinforcement
Some of the main physical parameters of the CFRP, GFRP, and AFRP are listed in Table 7. It can be found that CFRP has the largest tensile and elastic modulus values, whereas GFRP has the largest density and AFRP the largest ultimate elongation rate. The adopted scheme of reinforcement is shown in Table 8.


4.3.2. Dynamic Response
(a) Rock Impact Time History Curves. Figure 16 shows the impact time history curves for models 7 to 9. It can be seen that the peak values of the impact force of models 7 to 9 are 221 kN, 219 kN, and 213 kN, respectively. The magnitude of the impact forces and the duration of the interactions exhibit minor differences in these models. This is primarily owing to the different elastic moduli of the three FRP materials. Furthermore, as previously mentioned, the generated irregular patterns of the curves after the rock penetrates the wall are complex, and this is an outcome of the induced vibration owing to the specific model.
(a)
(b)
(c)
(b) Damage Response. Figure 17 shows the stress contours of models 7 to 9 at = 300 ms. It can be seen that the local damage is the main failure mode of the walls. In addition, a certain degree of stress concentration occurs at the four corners of the hole and at the bottom of the wall. By comparing the damage of these three models, model 8 yields the smallest final damage, while model 7 yields the most severe final damage. This finding possibly relates to the elastic modulus of the FRP. As it is shown in Tables 7 and 8, model 7 is reinforced by GFRP whose elastic modulus is 72.5 GPa, and model 8 is reinforced by CFRP whose elastic modulus is 230 GPa.
(a)
(b)
(c)
(c) Stress–Strain Response. Figure 18 shows the stress fields for models 7 to 9 when the velocity of the rock is reduced to zero upon impact. Figure 19 shows the strain of each model at the same instant in time. At that moment, the fiber stress and strain reach their maximum values. It can be seen that the stress values of the elements around the impact position also reach the value of the largest stress. The maximum stresses of models 7 to 9 are 290 MPa, 180 MPa, and 230 MPa, respectively. These values are much smaller than the tensile strengths of the fibers. The magnitude distribution pattern of the strain is coincident with that of stress. The maximum strain values of models 7 to 9 are 0.005, 0.012, and 0.009, respectively, which are also smaller than the FRPs’ ultimate strain. The conclusion that can be inferred is that increased energy is absorbed by the fiber, and increased fiber stress is developed when the elastic modulus is increased.
(a)
(b)
(c)
(a)
(b)
(c)
(d) Displacement Response. Figure 20 shows the displacement time history curves of selected typical nodes. From these curves, the displacement patterns of the nodes in three models are generally the same. Model 8 yields the largest displacement amplitude, and the largest final residential strain, while model 7 yields the smallest displacement amplitude and final residential strain. Because model 8 is reinforced by CFRP with a larger elastic modulus, it can absorb more energy upon deformation, resulting in the reduction of the displacement of the entire model. In addition, the time history curve of node 4 is also denser than the others, which means that node 4 has a higher vibration frequency. As it was mentioned before, node 4 attains a higher stiffness than other nodes.
(a)
(b)
(c)
(d)
Table 9 lists the maximum displacement of each typical node. The results from these data are consistent with the results presented in Figure 20. In general, CFRP performs better under debris flow, because of its high elastic modulus.

4.4. Relative Impact Position
4.4.1. Scheme of Reinforcement
To analyze the effect of different impact positions, two models were established. As shown in Figure 3, a diagonal FRP arrangement is adopted. The widths of the FRP of the two models are 150 mm, and the thickness values are 0.167 mm. The reinforcement scheme is shown in Table 10. Figure 21 shows the relative impact positions of the two models.

(a)
(b)
4.4.2. Dynamic Response
(a) Rock Impact Time History Curves. Figure 22 shows the impact time history curves for models 10 to 11. It can be seen that the peak values of the impact force of models 10 and 11 are 218 kN and 209 kN, respectively. This difference results from the different stiffness. The stiffness of the wall is a little larger at the center, because the area of the FRP at the crossposition is larger than that of other parts. Furthermore, as mentioned above, irregular curve patterns are elicited after the rock penetrates the wall owing to the model vibrations.
(a)
(b)
(b) Displacement Response. Figure 23 shows the displacement time history curves of four typical nodes. As seen, the displacement patterns of the nodes in the two models are generally the same. However, the value of the displacement in model 11 is much smaller than that of model 10. This is because the impact position of model 11 is close to the foundation of the structure. In terms of mechanics, for the forces that have the same value, the force that is closer to the anchorage end causes a smaller displacement. Table 11 lists the maximum displacements for each of the studied nodes. The results from these data are consistent with the results plotted in Figure 23.

(a)
(b)
(c)
(d)
(c) FRP Response. Figure 24 shows the stress distribution of FRP for the two models when maximum stress occurs. It can be seen that the maximum stress concentrates at the impact position. The maximum stress of FRP in model 10 is 293 GPa when = 130 ms, while its corresponding value is 340 GPa in model 11 when = 150 ms. This is because the impact is burdened by two FRP strips in model 10, but by one FRP strip in model 11. Figure 25 shows the stress time history curves of the maximum stress element in the two models. Based on the curves, it can be concluded that FRP is in the elastic stage of its deformation during the impact. The strain is calculated to be more than 0.0002, which is far lower that the ultimate tensile strain. Accordingly, FRP is not exposed to tensile failure. In other words, the assumption that FRP is elastic is proven.
(a)
(b)
(a)
(b)
In general, it can be concluded that the impact at the center is more unfavorable than at the corner. Consequently, it is reasonable to simulate the constructed models using an impact position set at the center of the wall.
5. Application Proposals
Based on the detailed parametric study and the discussion on the corresponding mechanism introduced in Section 4, some suggestions of engineering applications in debris flow hazard areas can be proposed:
For FRPreinforced masonry structures, the dynamic performance under debris flow varies based on different FRP arrangements and FRP material types. The importance of these differences should be reflected in the design codes in the future. However, further studies should be carried out to determine the details of these parameters.
The width of FRP is an essential parameter for reinforcing masonry structures. If civil engineers need to control local damage, narrower widths of FRPs are favorable, but wider FRP widths are more suitable in controlling the overall displacement. However, the achievement of a consistent onetoone match between values of the FRP width and the degrees of local damage, or the overall displacements, requires further studies.
Although there is a reinforcement code in China, FRP reinforcement regulations for masonry structures under debris flow cannot be referred. These regulations, such as the ones pertaining to failure criteria, demand more indepth studies. Suggestions on the failure criteria can be determined by the destruction of walls, failure of FRP, or by surpassing the limitation of the overall displacement.
6. Conclusion
In this study, the influences of different parameters on the dynamic response of FRPreinforced masonry structures under viscous debris flow are discussed. Because the application of FRP in resisting debris flow has been seldomly studied before, this study is valuable and necessary for its further application in masonry structures in resisting debris flow. The following conclusions are obtained based on parametric studies by analyzing the numerical simulation results:
The reinforced model performed better than the unreinforced model under debris flow. The parametric study on the dynamic response of FRPreinforced masonry structures under debris flow impacts was necessary.
The arrangement of reinforced FRP was an important parameter that influenced responses of masonry structures under debris flow. Generally, the diagonal and the intersecting parallel types of FRP arrangements yielded better performances than the horizontal type. The carbon fibers arranged in diagonal and intersecting parallel arrangements induced a more uniform wall stress distribution. Additionally, different arrangements also led to different shapes of local damage of masonry walls.
The widths of the reinforced FRP materials had a profound influence on the structural behavior. A wider fiber cloth led to a minimized structural response impacted by debris flow, and resulted in more severe local damage.
The elastic modulus was the dominant physical parameter that influenced the dynamic response of the masonry structure reinforced by different types of FRP. GFRP yielded the smallest elastic modulus value among all the three types and had the most significant effect in reducing the local damage, while the CFRP yielded a better structural response under debris flow because of its high elastic modulus, which ultimately helped the fiber absorb more energy.
The impact at the center was more unfavorable than that at the corner. Consequently, it was reasonable to conduct a parametric study on the dynamic response of FRPreinforced masonry structures with the debris impact at the center of the wall.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The financial support that was received from the National Key Technology R&D Program through Grant 2014BAL05B01 is highly appreciated.
References
 X. Tang, Debris Flows in China, The Commercial Press, Beijing, China, 2000.
 B.H. Jun, “Numerical simulation of the topographical change in Korea mountain area by intense rainfall and consequential debris flow,” Advances in Meteorology, vol. 2016, Article ID 9363675, 11 pages, 2016. View at: Publisher Site  Google Scholar
 W. Zhang, Q. Wang, J. Chen, H. Li, J. Que, and Y. Kong, “Grainsize analysis of debris flow alluvial fans in Panxi Area along Jinsha River, China,” Sustainability, vol. 7, no. 11, pp. 15219–15242, 2015. View at: Publisher Site  Google Scholar
 F. Tang, “Finite element analysis of seismic resistance of masonry structure in rural China,” Advanced Materials Research, vol. 827, pp. 328–333, 2014. View at: Publisher Site  Google Scholar
 Z. Lu, Y. Yang, X. Lu, and C. Liu, “Preliminary study on the damping effect of a lateral damping buffer under a debris flow load,” Applied Sciences, vol. 7, no. 2, article no. 201, 2017. View at: Publisher Site  Google Scholar
 S. Q. Yang, “Introduction of National Comprehensive Disaster Prevention and Reduction Plan,” Disaster Reduction in China, no. 1, pp. 20–22, 2017. View at: Google Scholar
 Z. Lu, D. Wang, S. F. Masri, and X. Lu, “An experimental study of vibration control of windexcited highrise buildings using particle tuned mass dampers,” Smart Structures and Systems, vol. 18, no. 1, pp. 93–115, 2016. View at: Publisher Site  Google Scholar
 Z. Lu, X. Chen, D. Zhang, and K. Dai, “Experimental and analytical study on the performance of particle tuned mass dampers under seismic excitation,” Earthquake Engineering & Structural Dynamics, vol. 46, no. 5, pp. 697–714, 2017. View at: Publisher Site  Google Scholar
 Z. Lu, X. Lu, H. Jiang, and S. F. Masri, “Discrete element method simulation and experimental validation of particle damper system,” Engineering Computations, vol. 31, no. 4, pp. 810–823, 2014. View at: Publisher Site  Google Scholar
 K. Dai, J. Wang, R. Mao, Z. Lu, and S.E. Chen, “Experimental investigation on dynamic characterization and seismic control performance of a TLPD system,” The Structural Design of Tall and Special Buildings, vol. 26, no. 7, Article ID e1350, 2017. View at: Publisher Site  Google Scholar
 Z. Lu, X. Chen, X. Lu, and Z. Yang, “Shaking table test and numerical simulation of an RC framecore tube structure for earthquakeinduced collapse,” Earthquake Engineering & Structural Dynamics, vol. 45, no. 9, pp. 1537–1556, 2016. View at: Publisher Site  Google Scholar
 Z. Lu, Z. Wang, S. F. Masri, and X. Lu, “Particle impact dampers: Past, present, and future,” Structural Control and Health Monitoring, vol. 25, no. 1, article e2058, 2018. View at: Google Scholar
 Z. Lu, X. Chen, and Y. Zhou, “An equivalent method for optimization of particle tuned mass damper based on experimental parametric study,” Journal of Sound and Vibration, 2017. View at: Publisher Site  Google Scholar
 Z. Lu, B. Huang, and Y. Zhou, “Theoretical study and experimental validation on the energy dissipation mechanism of particle dampers,” Structural Control and Health Monitoring, 2017. View at: Publisher Site  Google Scholar
 J. Wu, DebrisFlow and Comprehensive Management, Science Press, Beijing, China, 1993.
 G. Kuwabara and K. Kono, “Restitution coefficient in a collision between two spheres,” Japanese Journal of Applied Physics, vol. 26, no. 8R, pp. 1230–1233, 1987. View at: Publisher Site  Google Scholar
 American Association of State Highways and Transportation Officials, “LRFD bridge design specifications,” Washington, DC, USA, 1998. View at: Google Scholar
 O. Hungr, G. C. Morgan, and R. Kellerhals, “Quantitative analysis of debris torrent hazards for design of remedial measures,” Canadian Geotechnical Journal, vol. 21, no. 4, pp. 663–677, 1984. View at: Publisher Site  Google Scholar
 J. Xie, Study on the new technique of repairing and strengthening masonry walls with continuous carbon fiber sheet, Doctor’s Thesis, Tianjin University, Tianji, China, 2005.
 G. Schwegler, “Masonry construction strengthened with fiber composites in seismically endangered zones,” in Proceedings of the 10th European Conference on Earthquake Engineering, vol. 3, pp. 2299–2303, Vienna, Austria, 1994. View at: Google Scholar
 F. Seible, “Repair and seismic retest of a fullscale reinforced masonry building,” in Proceedings of the 6th International Conference on Structural Faults and Repair, vol. 3, pp. 229–236, Edinburgh, UK, 1995. View at: Google Scholar
 X. Wang, Z. Lu, and X. Lu, “An investigation of the aseismatic behavior of masonry structures strengthened by carbon fiber reinforced plastics,” Industrial Construction, vol. 33, no. 9, pp. 11–13, 2003. View at: Google Scholar
 D. Laigle, P. Lachamp, and M. Naaim, “SPHbased numerical investigation of mudflow and other complex fluid flow interactions with structures,” Computational Geosciences, vol. 11, no. 4, pp. 297–306, 2007. View at: Publisher Site  Google Scholar
 C. Zeng, P. Cui, and Y. Ge, “Characteristic and mechanism of buildings damaged by debris flows on 11 July, 2013 in Qipangou of Wenchuan, Sichuan,” Jourmal Earth Science and Environment, vol. 36, no. 2, pp. 81–91, 2014. View at: Google Scholar
 Y. Guo, Z. Yu, and L. Guo, “Impacting testing methods for masonry wall based on the drop hammer,” Structural Engineering, vol. 28, no. 6, pp. 123–127, 2012. View at: Google Scholar
 P. Li, T. Li, Z. Lu, and J. Li, “Study on dynamic response of novel masonry structures impacted by debris flow,” Sustainability, vol. 9, no. 7, article no. 1122, 2017. View at: Publisher Site  Google Scholar
 T.L. Teng, C.C. Liang, and T.T. Tran, “Development and validation of a finite element model for road safety barrier impact tests,” SimulationTransactions of the Society for Modeling and Simulation International, vol. 92, no. 6, pp. 565–578, 2016. View at: Publisher Site  Google Scholar
 O. I. Abdelkarim and M. A. ElGawady, “Dynamic and static behavior of hollowcore FRPconcretesteel and reinforced concrete bridge columns under vehicle collision,” Polymer, vol. 8, no. 12, article no. 432, 2016. View at: Publisher Site  Google Scholar
 C. Reuter, K.H. Sauerland, and T. Tröster, “Experimental and numerical crushing analysis of circular CFRP tubes under axial impact loading,” Composite Structures, vol. 174, pp. 33–44, 2017. View at: Publisher Site  Google Scholar
 G. S. Chandekar and A. D. Kelkar, “Experimental and numerical investigations of textile hybrid composites subjected to low velocity impact loadings,” The Scientific World Journal, vol. 2014, Article ID 325783, 14 pages, 2014. View at: Publisher Site  Google Scholar
 X. Wei and M. G. Stewart, “Model validation and parametric study on the blast response of unreinforced brick masonry walls,” International Journal of Impact Engineering, vol. 37, no. 11, pp. 1150–1159, 2010. View at: Publisher Site  Google Scholar
 H. M. Elsanadedy, T. H. Almusallam, S. H. Alsayed, and Y. A. ALSalloum, “Experimental and FE study on RC oneway slabs upgraded with FRP composites,” KSCE Journal of Civil Engineering, vol. 19, no. 4, pp. 1024–1040, 2015. View at: Publisher Site  Google Scholar
 D. Proske, J. Suda, and J. Hubl, “Debris flow impact estimation for breakers,” Georisk Assessment & Management of Risk for Engineered Systems & Geohazards, vol. 5, no. 2, pp. 143–155, 2011. View at: Publisher Site  Google Scholar
 N. He, N. S. Chen, and C. Zeng, “Current situation and tendencies of debris flow initiation mechanism,” Journal of Catastrophology, vol. 28, no. 1, pp. 121–125, 2013. View at: Google Scholar
 C. Zeng, P. Cui, Z. Su, Y. Lei, and R. Chen, “Failure modes of reinforced concrete columns of buildings under debris flow impact,” Landslides , vol. 12, no. 3, pp. 561–571, 2015. View at: Publisher Site  Google Scholar
 Ministry of Housing and UrbanRural Construction of the People’s Republic of China. Code for design of masonry structures GB 50032011, China Architecture & Building Press, Beijing, China, 2011.
 S. Zhang, D. Yang, Y. Sheng, S. W. Garrity, and L. Xu, “Numerical modelling of FRPreinforced masonry walls under inplane seismic loading,” Construction and Building Materials, vol. 134, no. Supplement C, pp. 649–663, 2017. View at: Publisher Site  Google Scholar
 G. Wang, Debris flow impact pressure signal processing by the wavelet analysis, Doctoral Thesis, University of Science and Technology of China, Hefei, China, April 2017.
 P. Li, J. Li, and T. Li, “The failure mechanism of masonry structure under debris flow impact,” Journal of Huazhong University of Science and Technology (Natural Science Edition), vol. 45, no. 7, p. 15+29, 2017. View at: Google Scholar
 S. Govindjee, G. J. Kay, and J. C. Simo, “Anisotropic modelling and numerical simulation of brittle damage in concrete,” International Journal for Numerical Methods in Engineering, vol. 38, no. 21, pp. 3611–3633, 1995. View at: Publisher Site  Google Scholar
 “LSDYNA Theoretical Manual,” Livermore Software Technology Corporation, May 2007, http://ftp.lstc.com/anonymous/outgoing/jday/manuals/lsdyna_971_manual_k_rev1.pdf. View at: Google Scholar
 S.L. Xu and Q. Fang, “Numerical analysis on blastresistant capacity of masonry walls retrofitted with elastomeric polymer and CFRP,” Journal of PLA University of Science and Technology (Natural Science Edition), vol. 11, no. 3, pp. 306–311, 2010. View at: Google Scholar
 Y. R. Guo and H. Zheng, “Numerical simulation of CFRP strengthened perforated masonry walls under blast loading,” Advanced Materials Research, vol. 640, no. 1, pp. 727–731, 2012. View at: Google Scholar
 J. B. Tang, K. H. Hu, G. D. Zhou, H. Y. Chen, X. H. Zhu, and C. Ma, “Debris flow impact pressure signal processing by the wavelet analysis,” Journal of Sichuan University (Engineering Science Edition), vol. 45, pp. 8–13, 2013. View at: Google Scholar
Copyright
Copyright © 2018 Peizhen Li 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.