#### Abstract

At present, the research on the new ground anchor counterfort grille dam is not perfect. In order to analyze and design the new grille dam accurately and reasonably, it is necessary to carry out theoretical research on the internal force calculation of the overall coordination of the grille beam-column under the impact of debris flow. Under the condition of full reservoir of debris flow, the grille beam-column meets the working principle of elastic foundation beam; therefore, the equation of the overall cooperative work of the grille beam-column is derived based upon Winkler’s theory of elastic foundation beam. Finally, combined with the engineering example, this calculation method is used for analysis; the results show that the bending moment of the grille beam-column increases compared with that regardless of considering the overall coordination, and the increased value of the bending moment at the bottom of the grille column and the middle of the grille beam is the largest. Therefore, from the perspective of safety, the overall coordination cannot be ignored when designing a blocking dam.

#### 1. Introduction

China is a region with frequent debris flow outbreaks, most of which are in the central and western regions. Every year, the loss of life and property caused by debris flow is extremely serious [1]. As a main protective measure, control structure is always the research focus. With the continuous progress of research, various control structures with different forms and functions have emerged. Liu et al. [2] put forward a new suspended debris flow grille-dam structure; the advantages and disadvantages of this type of dam are obtained by expounding its application in debris flow control engineering of rare debris flow and high sediment flood of debris flow. Su et al. [3] proposed a flexible stone cage intercepting dam structure; numerical simulation is performed using the software ANSYS LS-DYNA and focuses on the dynamic response of the structure to the impact of the debris flows, and the stress of the dam varies with height under the static calculation condition is obtained. Wang and Zhang [4] proposed a new debris flow grille-dam with friction damper and used ABAQUS finite element software to simulate the impact of debris flow block stone on this new dam. Li et al. [5] proposed the Gabion-pile-slab wall debris flow composite retaining structure; through the experiment, the impact force of each layer of pressure box over time is obtained, which shows that the structure has good stability and anti-overturning ability, and the blocking effect is obvious. Ran et al. [6] proposed a concrete filled steel tubular grille-dam, the dynamic response of the structure was studied by solid impact test, and the test process was numerically simulated. Li et al. [7] proposed a new steel-concrete composite retaining dam and compared the dynamic response of the traditional gravity dam and the improved new retaining dam by using experimental research and numerical simulation. Xie et al. [8] proposed a herringbone water-sediment separation structure, simplified the calculation model of the rib beam, and established the strength calculation formula of the rib beam. Chen et al. [9] proposed the closure dam according to the geographical conditions of Sichuan Province, used the latest method to evaluate the debris flow, and proposed the control measures of the closure dam according to the results. Zhang et al. [10] proposed an archless debris flow dam based on the mechanical characteristics of the arch and solved the internal force of the arch dam by using the legislative principle in structural mechanics.

In summary, the researchers put forward many new prevention measures for debris flow prevention project. However, the research on the new of retaining dam focuses on numerical simulation and experimental research, which needs a lot of time to establish and analyze the model. Therefore, in this study, Wang et al. [11] proposed a new ground anchor counterfort grille-dam structure (hereinafter called as the new grille-dam). In order to calculate the maximum internal force of grille beam-column more practical, under the condition of full reservoir flow, the debris flow impact force is simplified by calculation as concentrated load reverse action on the grille beam-column node, the grille beam-column is regarded as a whole, the elastic frame model of the whole grille beam-column is established by using the elastic foundation beam model, and then, through the displacement coupling of the connection node, the equation of the overall cooperative work of the grille beam-column is derived, and the internal force is calculated, which lays the foundation for the further analysis of the structural performance of the new grille-dam.

#### 2. Stress Analysis of New Grille-dam

##### 2.1. A Brief Introduction of New Grille-dam

The new grille-dam is mainly composed of grille (composed of grille columns, beams, and waste section steel beams), foundation piles, counterfort walls, lassos, anchor pier, floor, and floor beams, as shown in Figure 1. The grille structure of reinforced concrete grille beam and grille column large frame embedded steel beam is adopted to achieve the goal of debris flow interception and drainage. The ground anchor lassos can ensure the overall stability of the new grille-dam [11].

**(a)**

**(b)**

**(c)**

##### 2.2. Load Calculation of Debris Flow

According to [12], the calculation formula of debris flow impact force iswhere is the density of debris flow, is dynamic correction factor, and is debris flow velocity.

In this study, the vertical variation of debris flow impact force is simplified as triangular linear distribution [13].

##### 2.3. Simplification of the Calculation Model

The load acting on the plane normal line of the axis of the grille beam-column belongs to the calculation category of the spatial structure. The grille beam-column not only has bending moment but also has torque. However, in the actual grille beam-column system, the torsional stiffness of the grille beam-column is smaller than the flexural stiffness, which has little effect on the overall structure and can be ignored. Therefore, the main work of this study is to calculate the bending moment of the grille beam column under the overall coordination.

The vertical distribution of debris flow impact force shows the maximum at the bottom and gradually decays to the fluid surface [12]. In order to consider the overall safety, the grille column and the counterfort walls are cast-in-situ as a whole. The counterfort walls section is designed to be the largest at the bottom and the smallest at the top, *T*-section with a linear change in the middle, as shown in Figure 2. In order to simplify the calculation and convenient engineering application, the stiffness of each section of the grille column is equivalently calculated according to the size of midspan *T*-sections, as shown in Figure 3. The grille beam is simplified by equal stiffness multispan continuous beam.

Cross-sectional areas are

The distance, respectively, from the centroid of each section to the edge is

The distance from the centroid of any T-section to the edge can be expressed as

Substituting (2) and (3) into (4), we can obtain

The moment of inertia of the T-section of the wall can be expressed aswhere , , , , , and represent the cross-sectional area and size of any section of the grille column and the counterfort walls, is the distance from the centroid of any section of the grille columns and the counterfort walls to the side of the grille column, and is the distance from centroid to edge of any T-section.

#### 3. Internal Force Calculation of Grille Beam-Column

At present, there are many studies using Winkel elastic foundation beams to calculate the internal forces of structures at home and abroad. Mudhaffar et al. [14] used the high-order integral shear deformation theory to study the behavior of ceramic-metal plates in viscoelastic foundations under the action of damp heat and thermal loads, the bending equation was investigated, and the effect of temperature and water concentration on the metal plate was investigated. Merazka et al. [15] studied the damp-heat-thermal bending response of a simply supported FG slab on a Winkler–Pasternak elastic foundation under the condition of considering the transverse shear strain and not applying any shear correction factor. Based on the above concepts, in the case of full reservoir flow, debris flow impacts the new grille-dam instantly and debris flow deposits produce compression deformation. It is assumed that the compression deformation of deposits is generated by the reaction of the impact force of debris flow on the grille beam-column joints. Combined with Winkler's theory of elastic foundation beam, considering the constraint of side span grille column, the grille beam is regarded as a short beam with fixed ends, as shown in Figure 4; the grille column is regarded as a semi-infinite beam with fixed end for calculation, as shown in Figure 5.

On the basis of Winkler’s assumption, the flexible characteristic values of grille beams and columns are as follows [16]:where and is the rigid and flexible characteristic value of the grille beam-column, is the foundation coefficient, also called the foundation resistance coefficient, and is the single span of grille beam-column in the and directions, is the elastic modulus of the grille beam-column, and and is the section moment of inertia of the grille beam-column.

Assuming that the grille beam is parallel to the *x*-axis direction and the grille column is parallel to the *y*-axis direction, the bottom of any grille column is taken as the origin of the coordinate axis, and the concentrated force at the intersection of the *i*th grille beam and the *j*th grille column is set to , as shown in Figure 6. When a grille beam is subjected to a force, Wenkel assumes that the deformation at any position is as follows [17]:where is the force of the debris flow impact force distributed on the grille beam at the intersection of the *i*th grille beam and the *j*th grille column, is the distance between the concentrated force on the grille beam and any calculated section position , and is the distance between the concentrated force on the grille beam and the end point of the grille beam, that is, the origin of the grille beam.

The grille beam and grille column are subjected to multiple concentrated forces at the same time; the deformation of the *i*th grille beam generated by the superposition method under the action of concentrated forces is as follows [17]:

Orderthanwhere is the deformation coefficient of the *j*th node position of the *i*th grille beam and is the deformation of the *j*th node position of the *i*th grille beam under the action of .

The deformation of all the nodes of the grille beam can be expressed as a matrix:

Likewise, as shown in Figure 7, when a grille column is subjected to a force, Winkle assumes that the deformation at any position is as follows:where is the force of the debris flow impact force distributed on the grille column at the intersection of the *i*th grille beam and the *j*th grille column, is the distance between the concentrated force on the grille column and any calculated section position , and is the distance between the concentrated force on the grille column and the end point of the grille column, that is, the origin of the grille column.

The deformation of the *i*th grille column by the superposition method under concentrated forces is as follows:

Orderthanwhere is the deformation coefficient of the *i*th node position of the *j*th grille column and is the deformation of the *i*th node position of the *j*th grille column under the action of .

The deformation of all the nodes of the grille column can be expressed as a matrix:

According to the coordination conditions at the grille beam-column joints,

Comparing the displacement matrix of the grille column and the grille beam node, it is found that the two are transpose matrices, namely, .

Thus, (17) can be converted to the following:

Combining (12) and (19), and can be obtained by software. Finally, and are substituted into (9) and (14) to calculate their deformation, and then, the bending moment equation of the grille beam-column at any cross-section position is obtained by quadratic derivation of (9) and (14).

#### 4. Calculation and Analysis of Engineering Example

##### 4.1. Project Profile

The main design parameters of a new grille-dam are as follows. Concrete grade adopts , stressed bars for the grille beam-column adopt , stirrups adopt , seven strands of steel strand adopt , the standard value of single ultimate strength is , and section area is ; section of grille column is , the spacing is ; the counterfort wall is the bottom section , the upper section is , and the middle section changes linearly; section of grille beam is and the spacing is . In the case of full reservoir flow, the maximum velocity is , and the grille beam-column facade is shown in Figure 8.

##### 4.2. Impact of Overall Collaboration on Internal Force Calculation

In order to analyze the influence of the overall coordination on the internal force of the grille beam-column, in the calculation process, the formula considering the overall coordination is first calculated; then, according to [18] regardless the overall coordination, the calculation results are shown in Figures 9 and 10. Because of the symmetry of the grille column, the bending moment only represents the side column, the middle column 1, and the middle column 2. In the figure, the bending moment is positive in the direction of debris flow inflow. The following conclusions can be obtained from Figures 9 and 10.(1)The bending moment near the inner support of the end span of the grille beam is relatively large. Compared with not considering the overall coordination, the bending moment value of the grille beam increases, the maximum increase is 112.1 kNm, the bending moment value of the grille beam 1 is relatively large, and the maximum value is 546.3 kNm at the midspan position. The cross sections with zero moment are almost on the same vertical line, which indicates that the calculation of grille beams as continuous multispan beams is reasonable, and beam 1 is used as the benchmark in the design of grille beams. Since the velocity of debris flow in the lateral direction is large in the middle and small nonlinear changes on both sides [19], the grille beam 1 with the largest increase in the span direction has an increase of 58.1%, 66.6%, 99.6%, 99.6%, and 66.6%, respectively. The maximum increase position is the midspan, so the design method of variable cross section can be considered when designing grille beams(2)Under the impact of debris flow, the bending moment at the bottom of the grille column is the largest, which is 698.5 kNm, the bending moment at the upper part is smaller, which is 176.5 kNm, the middle bending moment gradually decreases along the height direction of the grille column, and the decreasing range is 32.5%, 20.5%, 15.4%, 15.3%, and 15.0%, and the decreasing range gradually decreases, because with the increase of height, the stiffness of the grille column gradually decreases. It is small, so the distributed bending moment is reduced, and the impact force of the debris flow is the largest at the bottom and then gradually decays to the fluid surface. The bending moment of the side column is smaller than that of the central column, and after considering the coordination effect, the bending moment values of the central column 1 and 2 both increase, and the bending moment value of the side column shows a decreasing trend due to the constraint effect of the grille beam. It shows that the bearing capacity of the central column should be mainly strengthened in the design of the grille-dam structure.(3)When do not consider overall coordination, the bending moments of the side column and the center column of the grille-dam are exactly the same. Compared with the results considering the overall coordination, the grille column design without considering the overall coordination is likely to have bottom bending under the impact of debris flow. Shear failure: considering the overall coordination effect, since the stiffness of the grille column is much larger than that of the grille beam, the increase of the bending moment of the grille column is greater than that of the grille beam; when the model is simplified, the end is simplified as a fixed support, so there will be inflection points in the grille beams and columns at equal intervals in the *X* and *Y* directions.

**(a)**

**(b)**

**(a)**

**(b)**

#### 5. Conclusions

(1)In this study, the elastic foundation beam theory is used to establish the internal force equation of grille beams and grille columns, and the calculation equation of the overall coordinated internal force of the grille dam is established through equal displacements at the nodes of grille beams and columns, and the obtained results better reflect the grille beams and variation of the internal force of the column.(2)Engineering analysis and calculation show that, since the stiffness of the grille columns is greater than that of the grille beams, considering the overall coordination effect of the grille columns, the distributed bending moment value is larger, and the increase is more obvious. The increase of the bending moment in the bottom of the grille column and the middle of the grille beam is the largest; the increase in the bending moment of the central column is larger than that of the side column. From the perspective of durability and safety, the overall coordination must be considered when designing new grille dam.#### Data Availability

Due to privacy reasons, the data of this paper are temporarily not open to the public.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest or personal relationships that could have appeared to influence the work reported in this paper.

#### Authors’ Contributions

Yongsheng Wang conceptualized the study, developed methodology, collected resources, administrated the project, carried out funding acquisition, and supervised the study. Hao Zhu conceptualized the study, developed the methodology, curated the data, validated the study, and wrote the original draft. Bing Wang carried out proof reading.

#### Acknowledgments

This work was financially supported by National Natural Science Foundation of China (Grant no.51768039), Industrial Support Program of Higher Education of Gansu Province (2020C-40), and Hongliu Support Funds for Excellent Youth Talents of Lanzhou University of Technology (Grant nos. 04-062002) and are also highly appreciated.