Research Article  Open Access
Changjiang Liu, Fan Wang, Jian Liu, Xiaowei Deng, Zuoliang Zhang, Haibing Xie, "Theoretical and Numerical Studies on Damped Nonlinear Vibration of Orthotropic Saddle Membrane Structures Excited by Hailstone Impact Load", Shock and Vibration, vol. 2019, Article ID 9234832, 21 pages, 2019. https://doi.org/10.1155/2019/9234832
Theoretical and Numerical Studies on Damped Nonlinear Vibration of Orthotropic Saddle Membrane Structures Excited by Hailstone Impact Load
Abstract
The orthotropic membrane structures have been popular in architectural structures. However, because of its lightweight and small stiffness, large nonlinear deflection vibration may occur under impact load, which leads to structural failure. In this paper, the governing equations of the large deflection nonlinear damped vibration of orthotropic saddle membrane structures excited by hailstone impact load are proposed according to the von KĆ”rmĆ”nās large deflection theory and solved by applying the BubnovāGalerkin method and the method of KBM perturbation. The approximate theoretical solution of the frequency function and displacement function of the large deflection nonlinear damped vibration of saddle membrane structures with four edges fixed excited by hailstone impact was obtained. The analytical examples proved that the mode shape function (equation (43)) can be applied to calculate the singleorder mode shapes and the total superposed mode shapes of the damped large nonlinear deflection vibration of orthotropic saddle membrane structures excited by hailstone impact load succinctly. In addition, we compare and analyze the results of vibration frequency, amplitude, time histories, and total displacement of membrane structures with different pretensions and archtospan ratios under the impact of differently sized hailstones. The correctness of the analytical theory is verified by comparing with the results of numerical simulation. According to the results of this paper, we put forward some suggestions for the vibration control and dynamic design of practical spatial membrane structures.
1. Introduction
The membrane structure which is made of orthotropic membrane materials has been popular in architectural structures, and the saddle membrane structure is one of the most common shapes of the doublecurved membrane [1, 2]. However, because of its lightweight, large flexibility, and small stiffness, it is very sensitive to impact load. Large nonlinear deflection vibration occurs under impact load, which may lead to structural failure [3, 4]. Thus, it is necessary to study damped large nonlinear deflection vibration of orthotropic saddle membrane structures excited by hailstone impact load.
In recent decades, more and more attention has been focused on the dynamic characteristics of the membrane. By applying the Hamilton principle and Galerkin method, Shin et al. [5, 6] obtained the natural frequencies and mode shapes of the free vibration for an axially moving membrane. The results showed that the translating speed, aspect ratio, and boundary conditions have significant effects on the inplane vibrations of the moving membrane. Pan and Gu [7] adopted DāAlembertās principle to deduce the free oscillating systemās equivalent fundamental frequency of the square tensioned membrane. The effects of prestrain, size, elastic ratio, density, relative amplitude, and dead load on the nonlinearity of square pretensioned membrane were studied. Zheng et al. [8], Liu et al. [9], and Li et al. [10] investigated the large deflection nonlinear free vibration of orthotropic rectangular membrane structure by applying von KĆ”rmĆ”nās large amplitude theory, DāAlembertās principle, BubnovāGalerkin approximate method, and LindstedtāPoincarĆ© perturbation method and obtained the approximate analytical solution in power series of the nonlinear vibration frequency function and the displacement function of the rectangular membrane with four edges fixed. However, they only studied the free vibration of planar rectangular membrane structures.
On the basis of researches of the nonlinear free vibration of membrane structures, the investigations of the nonlinear forced vibration of membrane structures were performed. GonĆ§alves et al. [11] derived the equations of motion of the prestretched hyperelastic isotropic membrane with finite deformations and lateral pressure and obtained analytically the functions of natural frequencies and mode shapes of the membrane. The results show the stretching ratio significantly affected the linear and nonlinear vibrations of the membrane. Zheng et al. [12] applied DāAlembertās principle and the momentum theorem to derive the fundamental equations of forced vibrations of orthotropic membranes, which were solved according to the LindstedtāPoincarĆ© perturbation method, and obtained the formula of impact load and nonlinear forced vibration deflection of rectangular membranes with four edges fixed. However, the damping was not considered. By applying the KrylovāBogolubovāMitropolsky (KBM) perturbation method to solve the governing equations of largeamplitude nonlinear vibration of rectangular orthotropic membranes with viscous damping, Liu et al. [13] obtained the asymptotic analytical solutions for the frequency and displacement function of the planar rectangular membrane structure. However, the external load was not taken into consideration. According to the FĆ¶ppl large deflection theory, Galerkin method, and multiple scale perturbation method, Zheng et al. [14] investigated the dynamic response of rectangular prestressed membrane subjected to concentrated impact load. Based on the von KĆ”rmĆ”n large deflection theory, DāAlembertās principle, the BubnovāGalerkin method, and perturbation method, Liu et al. [15ā17] studied nonlinear forced vibration of pretension rectangular orthotropic membrane structures with damping and without damping under single impact load. Based on the stochastic pulse theory and the random vibration theory, Zheng et al. [18] and Li et al. [19, 20] investigated the stochastic vibration problem of the orthotropic membrane subjected to random impact load through experimental and theoretical researches. Based on thinplate theory and the Galerkin method, Li et al. [21] investigated the dynamic response of rectangular prestressed membrane subjected to uniform impact load theoretically and experimentally. The aforementioned studies, however, are limited to planar rectangular membrane structures.
On the current status, the studies on saddle membrane structure are concentrated on wind or rain loads. Yang and Liu [22] and Li and Sun [23] studied the aerodynamic critical unstable wind velocity of saddle membrane structure by applying the nonmoment theory and the potential flow theory. The aerodynamic interaction equations of the membrane structure were obtained and simplified by applying the BubnovāGalerkin approximate method. By studying the free vibrations of membrane structures with the static and the dynamic effects of wind and snow, Lazzari et al. [24] analyzed the structural failure mechanisms of the roof of the Montreal Stadium membrane structure. Rizzo et al. [25] and Rizzo and Sepe [26] measured the pressure of the incoming wind on hyperbolic paraboloid roofs by conducting tests and finite element analyses. The possibility of defining equivalent static pressure fields able to reproduce the envelope of dynamic displacements of the cables net was explored. Wu et al. [27] studied the aeroelastic instability mechanism of a tensioned membrane structure. The response and wind velocities above two closedtype saddleshaped tensioned membrane structures, with the different pretension levels, were measured in uniform flow and analyzed. The results indicate that the aeroelastic instability is caused by vortexinduced resonance. Xu et al. [28, 29] and Liu et al. [30] studied nonlinear windinduced aerodynamic stability of orthotropic saddle membrane structures by establishing the interaction governing equations of windstructure coupling based on von KĆ”rmĆ”nās large amplitude theory and DāAlembertās principle. They determined the critical velocity of divergence instability, by judging the stability of the characteristic equation of the system. Cui et al. [31] applied the EulerianāEulerian model according to multiphase flow theory, and the saddle membrane structure response analysis under the simultaneous actions of wind and rain was conducted. The influences of changes in wind speed and rain intensity on the saddleshaped membrane structure response were compared. There are no researches about the problem of saddle orthotropic membranes under impact load.
In this paper, the approximate formulas of hailstone terminal velocity were substituted into the governing equations of the large deflection nonlinear damped vibration of orthotropic saddle membrane structures excited by impact load. And, solving the governing equations by applying the BubnovāGalerkin method and the method of KBM perturbation, the approximate theoretical solution of the frequency function and displacement function of the large deflection nonlinear damped vibration of saddle membrane structures with four edges fixed excited by hailstone impact was obtained. In analytical examples, the dynamic responses of saddle membrane structures with different pretension levels and archtospan ratios excited by the impact of different diametral hailstones were compared and analyzed separately. The correctness of the analytical theory is verified by comparing with the results of numerical simulation. In addition, the results of this paper can be applied in computation for the vibration control and dynamic design of practical spatial membrane structure under impact load.
2. Modeling of Saddle Membrane Structure
In this paper, we study a saddle, namely, hyperbolic paraboloid, membrane structure with four edges simply supported under an impact load. The theoretical model of saddle membrane structure is shown in Figure 1. The orthogonal axes x and y are the two different Youngās modulus fiber directions of orthotropic saddle membrane structures. a and b, respectively, are the spans in x and y axes. and , respectively, denote the pretension in x and y axes. and , respectively, are the midspan arch in x and y axes. The point O is the center of the plane xoy. The sphere H is a hailstone; denotes the velocity of the hailstone; (, ) is the impact point on membrane surface.
The saddle membrane model can be represented by [29]where z_{0} denotes the initial surface function of saddle membrane structure.
According to equation (1), the two initial principal curvatures in x and y directions are
With the action of the pretensions N_{0x} and N_{0y} [28], we can obtain
3. Dynamic Governing Equations
According to the von KĆ”rmĆ”nās large deflection theory and DāAlembertās principle, the compatible equation and dynamic motion equation of orthotropic saddle membrane structures are [28]where is the areal density of membrane; and , respectively, are the stress increments in x and y directions; is the shear stress; and , respectively, are the initial stress in x and y directions; is the lateral deflection; h is the thickness of membrane; and , respectively, are Youngās modulus in x and y directions; G is the shear modulus; and , respectively, are Poissonās ratio in x and y directions; and , respectively, are the two principal curvatures in x and y directions; and , respectively, are the two initial principal curvatures in x and y directions; c is the damping coefficient of structure; is the impact load in z direction.
According to the basic theory of plates and shells, the principal curvatures in x and y directions are [21]where and denote the principal curvature increments in x and y directions, respectively.
By means of Airyās stress function , we can obtain [30]
The effect of the shear stress is so small that we may assume that ā=ā0. Therefore, according to equation (6), we can obtain
By substituting equation (3) and equations (5)ā(7) into equation (4), we can obtain
Assume that the membrane does not bear an external load apart from the impact load; then, when the point of impact lies in the coordinates on () (as shown in Figure 1), the equation becomes [10]where is the impact force and is the impulse function.
Assume that the initial displacement of the membrane is zero before the membrane is impacted by hailstone; when the hailstone is just in contact with the membrane, at the moment of tā=ā0, they have the same velocity, and this velocity is the initial velocity of the membrane, so the initial conditions are
Impact loads are shortduration loads, so the maximum displacement amplitude depends principally upon the magnitude of the applied impulse and is not strongly influenced by the form of the load; thus, we can obtain the relational expression between and by the impulse theorem:where M is the mass of hailstone.
Derivation calculus to equation (12) yields
Equations (8)ā(13) are the fundamental equations applied in the analysis of saddle membrane excited by impact load.
According to the approximate formulas of hailstone terminal velocity [32], we can obtain the mass and velocity of hailstones (as shown in Table 1 and Figure 2).

The boundary conditions of the saddle membrane structure with four simply supported edges are
More concretely, the corresponding displacement and stress boundary conditions of every edges are
4. Solution of Fundamental Equations
The functions that satisfy the displacement conditions of every edges equation (15) are separated as follows [12, 13]:where is the given deformation function; is the function of time; m and n are the positive integer.
According to the basic vibration theory and boundary conditions, the displacement function is given bywhere m and n, respectively, denote the orders of vibration displacement in x and y directions.
We take one term of equation (18) for computation; i.e.,
Let and .
The substitution of equation (19) into equation (9) yields
The stress function should satisfy the stress boundary condition (16), on the basic of the differential equation theory. Therefore, the general solution of may be assumed as follows:
Let
The substitution of equation (21) into equation (16) yieldswhere are the arbitrary constant and can be set to zero to simplify the computation according to differential equations theory.
Then, the solution of is
The substitution of equations (24), (19), and (10) into equation (8) yieldswhere .
By applying the BubnovāGalerkin method [4], equation (25) can be transformed intowhere
Equation (26) can be simplified into a homogeneous differential equation as follows:where
The KBM perturbation method [13] is applied to solve equation (28). Assume the perturbation parameter is ; then, equation (28) can be simplified by letting as follows:where
According to the perturbation method of KBM, let and the solution of equation (30) is
In equation (32), and are determined bywhere
The substitution of equaiton (34) into equation (33) yieldswhere denotes the amplitude of vibration and denotes the initial phase of vibration. They can be determined by the initial conditions of membrane vibration.
The substitution of equation (35) into equation (32) yields
Let
Equation (36) is the frequency function of vibration. According to equation (36), we can conclude that the amplitude of vibration and damping coefficient of structure have effect on the frequency of the nonlinear damped forced vibration of orthotropic saddle membrane structure.
The expression of initial conditions of membrane vibration can be obtained according to the principle of conservation of momentum:where ; is the initial velocity of hailstone; is the initial velocity of the system that consisted of hailstone and the membrane impact point; is the mass of hailstone; is the density of hailstone.
The substitution of equation (38) into equation (36) yields
In equation (40), the amplitude of vibration is and there will undoubtedly be ; therefore, . Taking and substituting it into equation (36) yields
The substitution of the first derivative of equation (41) into equation (39) yields
We solved equation (42) by applying the root formula of simple cubic equation and obtained the real root expression of amplitude:where
When and we can only obtain the amplitude of vibration .
By substituting of equation (41) into equation (17), the displacement expression is obtained:where is determined by equation (43). According to equation (45), we can obtain the mode of displacement of membrane surface.
Superposition of the initial surface function of saddle membrane structure equation (1) and its displacement expression equation (45), we can obtain the mode shape of the saddle membrane structure excited by impact of hailstone.where denotes the mode shape of the saddle membrane structure.
According to equation (46), we can obtain the modes of vibration and displacement time histories of the saddle membrane structures.
5. Analytical Examples
We take the orthotropic membrane that is widely used in practical engineering application as the analytical example: E_{1}ā=ā1400āMPa, E_{2}ā=ā900āMPa [16], the areal density of membrane ā=ā1.7ākg/m^{2}, the thickness hā=ā1.0āmm, the length aā=ā1.0ām, the width bā=ā1.0ām, and the viscous damping cā=ā120āNs/m. We take the center point of membrane (i.e., , ) as the impact point that is excited by hailstone.
5.1. Displacement Time Histories of Single Order
The archtospan ratio f_{1}ā=āf_{2}ā=ā1/10, hailstone diameter dā=ā5.0ācm, and pretension Nā=ā1000āN/m. According to equation (2), we can obtain and . According to equation (45), the impact pointās first threeorder vibration time histories are shown in Figure 3. From Figure 3, we can conclude that(i)The maximum vibration amplitude decreases gradually with the increase of the vibration order (i.e., concluded from first order, second order to third order). The maximum vibration amplitude decreases the rate of first order to second order to 16.8% and second order to third order to 8.5%, and the decreases are nonlinear.(ii)With the increasing of time, the maximum vibration amplitude of single order decreases gradually until it reaches zero.
(a)
(b)
(c)
5.2. Computation of Amplitude of Impact Point
The archtospan ratio are 1/12 and 1/10, the pretension levels increase from 1000 to 4000āN/m, and the hailstone diameters increase from 1.0 to 6.0ācm. According to equation (45), the max. amplitude of the impact point when tā=ā0 is obtained. The results of the first three orders are presented in Table 2.

Figures 4ā6 show the results of Table 2. According to Table 2 and Figures 4ā6, we can come to the following conclusions:(i)When the archtospan ratio is 1/10 and the pretension level is 3000āN/m, the singleorder maximum amplitude of the impact point decreases with respect to increasing vibration order; the singleorder maximum amplitude of the impact point increases with respect to hailstone diameter increasing, and the increase is nonlinear (i.e., the increment of the maximum amplitude became smaller and smaller).(ii)When the archtospan ratio is 1/10 and the hailstone diameter is 5.0ācm, the singleorder maximum amplitude of impact point decreases with respect to the increasing pretensions and the decrement of the maximum amplitude became smaller and smaller. This reflects the nonlinearity of the of the membrane vibration.(iii)When the pretension level is 3000āN/m and hailstone diameter is 5.0ācm, the impact point maximum amplitude of single order decreases with respect to increasing archtospan ratio and the decrease is nonlinear.
5.3. Total Displacement Time Histories
The archtospan ratio is 1/10, and the pretension level is 1000āN/m. According to equation (45), the total displacement time history of impact point is shown in Figure 7, when hailstone diameter is 5.0ācm.
From Figure 7, we can observe that the impact point did weaken the vibration; i.e., the amplitude rapidly increases to its maximum when the membrane is excited by hailstone. Soon afterwards, with the increase of time, the maximum amplitude decreases gradually until it reaches zero.
5.4. Computation of Total Displacement of Impact Point
The archtospan ratios are 1/12 and 1/10, the pretension levels increase from 1000 to 4000āN/m, and the hailstone diameters increase from 1.0 to 6.0ācm. According to equation (45), the total displacement values of the impact point are calculated and listed in Table 3.

Figures 8ā10 show the results of Table 3. According to Table 3 and Figures 8ā10, we can come to the following conclusions:(i)The max. total displacement of the impact point increases with the increase of hailstone diameter. When the archtospan ratio is 1/10 and the pretension level is 3000āN/m and with the hailstone diameter increasing from 1.0 to 6.0ācm, the growth rates of maximum total displacement are 648%, 222%, 182%, 107%, and 83%, respectively. The growth rate of the max. total displacement became smaller and smaller.(ii)As the membrane surface is subjected to effect of stress stiffening, the max. total displacement of impact point decreases with respect to increasing pretension levels. When archtospan ratio is 1/10 and the hailstone diameter is 6.0ācm, and with the pretension levels increasing from 1000 to 4000āN/m, the decrease rates of maximum total displacement are 8.3%, 8.1%, and 7.5%, respectively. The decrease rate of the max. total displacement became smaller and smaller. This reflects the nonlinearity of the membrane vibration.(iii)The max. total displacement of impact point increases with the increase of archtospan ratio, and the increase is nonlinear.
5.5. Computation of Frequency
5.5.1. Hailstone Diameter
The archtospan ratio is 1/10, and the pretension level is 3000āN/m. The hailstone diameters increased from 1.0 to 6.0ācm. According to equation (37), the singleorder frequency is affected by the hailstone diameter and time. The frequencies of the first three orders with different hailstone diameters and time instants are calculated and listed in Table 4.

Figures 11 and 12 show the results of Table 4. According to Table 4 and Figures 11 and 12, we can come to the following conclusions:(i)In the case of the hail diameter being constant, the vibration frequencies increase with respect to the increasing vibration order. The frequencies of each order are maximum at tā=ā0 and gradually decrease to with increase of time. This reflects the nonlinearity of the damped forced vibration of the membrane.(ii)When tā=ā0, according to equation (37), the vibration frequency value is dependent on and amplitude. Amplitudes of each order increase with the increase of hailstone diameter. This shows that the vibration frequency value of tā=ā0 decreases with respect to increasing hailstone diameter. But, the higherorder vibration (i.e., begin with the 3rd order) frequencies value of tā=ā0 increase when hailstone diameter increases from 5.0 to 6.0ācm. According to equation (43) and equation (38), in terms of the velocity and the mass of hailstone, the velocity has a higherorder (i.e., ) effect on vibration amplitude and frequency than the mass (i.e. m^{ā1}), which is more significant in higherorder vibration. When hailstone diameter increased from 5.0 to 6.0ācm, as shown in Table 1, the velocity squared difference increases rather than decreasing. Namely, the increase of velocity determines the increase of frequency as hailstone diameter increased from 5.0 to 6.0ācm. This is the reason why the higherorder vibration frequency value of tā=ā0 increases when hailstone diameter increased from 5.0 to 6.0ācm.(iii)According to equation (37), when , the vibration frequency value is only dependent on . The vibration frequency value of each order decreases with respect to increasing hailstone diameter.
5.5.2. Pretension Levels
The archtospan ratio is 1/10, and the hailstone diameter is 5.0ācm. The pretension levels increased from 1000 to 4000āN/m. Table 5 shows the first threeorder frequencies with different pretension levels and times, which are calculated according to equation (37).

Figure 13 shows the results of Table 5. According to Table 5 and Figure 13, we can conclude that the singleorder frequencies increase with the increasing pretension levels. Moreover, the frequency increment of each order is bigger and bigger, so the increase is nonlinear.
5.5.3. ArchtoSpan Ratio
The pretension level is 3000āN/m, and the hailstone diameter is 5.0ācm. The archtospan ratio increased from 1/12 to 1/10. According to equation (37), the first threeorder frequencies with different archtospan ratios and times are calculated and presented in Table 6.

Figure 14 shows the results of Table 6. According to Table 6 and Figure 14, we can conclude that the vibration frequencies of each order decrease with the decreasing archtospan ratio.
5.6. Mode Shape
The hailstone diameter is 6.0ācm, archtospan ratio is 1/10, pretension level is 1000āN/m, and time tā=ā0.002ās and 0.007ās. According to equation (46), the first threeorder mode shapes are presented in Figures 15ā20. By superimposing the first threeorder mode shapes, the superposed mode shapes are obtained and shown in Figures 21 and 22. The coordinate dimension of Figures 15ā22 is in meter.
By analyzing the mode shapes, we can conclude that the mode shape function equation (46) can be applied to calculate the singleorder mode shapes and the total superposed mode shapes of the damped large nonlinear deflection vibration of orthotropic saddle membrane structures excited by hailstone impact load succinctly. In practical engineering vibration, the loworder mode shapes have a much bigger affect than the higherorder mode shapes, so we take the first threeorder mode shapes into analysis.
5.7. Brief Summary
(1)The archtospan ratio is 1/10, and the pretension level is 3000āN/m. The hailstone diameters increased from 1.0 to 6.0ācm. We plot the first trough of time histories of the firstorder vibration in Figure 23(a). From Figure 23(a), we can conclude that the hailstone diameter increase has a significant effect on the increase of amplitude and small effect on the decrease of frequency and increase of period.(2)The hailstone diameter is 5.0ācm, and the archtospan ratios are 1/12 and 1/10. The pretension levels increased from 1000 to 4000āN/m. We plot the first trough of time histories of the firstorder vibration in Figure 23(b). From Figure 23(b), we can conclude that the archtospan ratio decrease has a more significant effect on the increase of amplitude and the increase of period than that of decreasing pretension.
(a)
(b)
6. Numerical Simulation
In this section, we apply the universal explicit dynamics finite element analysis software ANSYS/LSDYNA to simulate the process of hailstone impacting on the membrane based on the explicittoimplicit sequential solution method. In implicit computation, Shell181 element is applied for the membrane and Solid185 element is applied for the hailstone. In explicit dynamic analysis, the elements will be converted to Shell163 and Solid164 element accordingly. We adopt a mapping triangleshaped shell element to generate the mesh of the membrane surfaces and a mapping hexahedralshaped solid element to generate the mesh for the hailstone. In ANSYS/LSDyna, the simulation results can converge better by refining the finite element mesh [33]. When hailstone diameter is 6.0ācm and archtospan ratio is 1/10, 5408 shell elements are generated with element size 2ācm and 4000 solid elements are generated with element size 0.5ācm, which satisfy the convergence accuracy. The results of the meshing are shown in Figures 24 and 25. The pretension load 1000āN/m is applied to the shell elements in implicit computation and converted to explicit dynamic analysis. In explicit dynamic analysis, the velocity according to Table 1 is applied to hailstone, automatic surfacetosurface contact is defined, and the corresponding parameter of the contact surface is set to describe the complex interaction among membranes in the large deformation contact and dynamic impact of hailstone. The Rayleigh damping coefficient is set to 0.05 that is recommended in explicit dynamic analysis.
Dynamic response process of hailstone impacting on membrane is shown in Figure 26. Time histories of the vertical displacement of impact point are shown in Figure 27. The following conclusions are drawn based on Figures 26 and 27:(i)At the initial time (tā=ā0), the stress distribution in the membrane is more uniform. The maximum stress appears near the highest point and the lowest point and the minimum stress occurs near four corners. The mean stress in the membrane approximately is 1000āN/m.(ii)The membrane material is elastic; it can only dissipate energy by deforming after being subjected to impact load. Therefore, when the hailstone impacted the membrane surface, most of the energy of hailstone is converted into the membrane strain energy and the kinetic energy of the vibration [12], which results in the rapid increment of vertical displacement of the impact point. At the same time, the stress concentration occurs in the impact area and the membrane stress diffuses from the impact point to the two high points.(iii)When the displacement of impact point is maximized, because of the pretension, the hailstone is rebounded up from the membrane and the membrane begins to vibrate freely. Subsequently, the displacement of the membrane attenuates gradually, and the vibration energy wave diffuses from the center to all around but rebounds at the boundary of membrane. This is the reason why the membrane stress exhibited irregularity.
(a)
(b)
(c)
(d)
(e)
(f)
7. Comparison and Analysis of Numerical Simulation and Theoretical Solution
The comparison of numerical simulation results and the theoretical results shows that the simulation results tally with the actual theoretical situation.
Figure 28 shows the results of Table 7. According to Table 7 and Figure 28, we can come to the following conclusions:(i)The nonlinear dynamic response law of orthotropic saddle membrane structures excited by hailstone impact load which reflected in the numerical simulation results is consistent with the theoretical ones.(ii)We neglected the dead load of membrane in the theoretical calculation. As a result of this, the theoretical results are slightly smaller (the maximum difference is less than 2.63āmm) than those of numerical simulation. The dead load of membrane is much more than the impact load of a hailstone when the diameters of hailstone are low; this is the main reason why the relative errors of simulation and theory are high at low hailstone diameters.(iii)The theoretical calculation results basically fit the numerical simulation results, implying that our methods have been successful.

8. Conclusions
In this paper, the approximate formulas of hailstone terminal velocity were substituted into the governing equations of the large deflection nonlinear damped vibration of orthotropic saddle membrane structures excited by impact load. And, solving the governing equations by applying the BubnovāGalerkin method and the method of KBM perturbation, the approximate theoretical solution of the frequency function and displacement function of the large deflection nonlinear damped vibration of saddle membrane structures with four edges simply supported excited by hailstone impact was obtained.
The analytical examples proved that the mode shape function equation (46) can be applied to calculate the singleorder mode shapes and the total superposed mode shapes of the damped large nonlinear deflection vibration of orthotropic saddle membrane structures excited by hailstone impact load succinctly. In addition, we compare and analyze the results of vibration frequency, amplitude, time histories, and total displacement of membrane structures with different pretensions and archtospan ratios under the impact of different size hailstones, and the following conclusions can be drawn:(i)The increasing hailstone diameter has a significant effect on the increase of amplitude and small effect on the decrease of frequency and increase of period(ii)The decreasing archtospan ratio has a more significant effect on the increase of amplitude and the increase of period than that of decreasing pretension
The correctness of the analytical theory is verified by comparing with the results of numerical simulation. In addition, the results of this paper can be applied in computation for the vibration control and dynamic design of practical spatial membrane structure under impact load. Therefore, we put forward some suggestions for the vibration control and dynamic design of practical spatial membrane structures:(i)In the preliminary design phase, the membrane structure with high archtospan ratio should be adopted as far as possible for a strong resistance to external load and vibration(ii)After the archtospan ratio of membrane structure was determined, increasing the pretension helps resist external load and vibration control
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (Project Numbers 51608060, 51678168, and 51878586), the Natural Science Foundation of Guangdong Province (Project Number 2017A030313267), and the Science and Technology Plan of Guangzhou City (Project Number 201607010107).
References
 K. Ando, A. Ishii, T. Suzuki, K. Masuda, and Y. Saito, āDesign and construction of a double membrane airsupported structure,ā Engineering Structures, vol. 21, no. 8, pp. 786ā794, 1999. View at: Publisher Site  Google Scholar
 P. D. Gosling, B. N. Bridgens, A. Albrecht et al., āAnalysis and design of membrane structures: results of a round robin exercise,ā Engineering Structures, vol. 48, pp. 313ā328, 2013. View at: Publisher Site  Google Scholar
 J. Gade, R. Kemmler, M. Drass, and J. Schneider, āEnhancement of a mesoscale material model for nonlinear elastic finite element computations of plainwoven fabric membrane structures,ā Engineering Structures, vol. 177, pp. 668ā681, 2018. View at: Publisher Site  Google Scholar
 Y. Tang and T. Li, āEquivalentforce density method as a shapefinding tool for cablemembrane structures,ā Engineering Structures, vol. 151, pp. 11ā19, 2017. View at: Publisher Site  Google Scholar
 C. Shin, J. Chung, and W. Kim, āDynamic characteristics of the outofplane vibration for an axially moving membrane,ā Journal of Sound and Vibration, vol. 286, no. 45, pp. 1019ā1031, 2005. View at: Publisher Site  Google Scholar
 C. Shin, W. Kim, and J. Chung, āFree inplane vibration of an axially moving membrane,ā Journal of Sound and Vibration, vol. 272, no. 12, pp. 137ā154, 2004. View at: Publisher Site  Google Scholar
 J. Pan and M. Gu, āGeometric nonlinear effect to square tensioned membraneās free vibration,ā Journal of Tongji University Natural Science, vol. 35, pp. 1450ā1454, 2007. View at: Google Scholar
 Z. L. Zheng, C. J. Liu, X. T. He, and S. L. Chen, āFree vibration analysis of rectangular orthotropic membranes in large deflection,ā Mathematical Problems in Engineering, vol. 2009, Article ID 634362, 9 pages, 2009. View at: Publisher Site  Google Scholar
 C. J. Liu, Z. L. Zheng, X. T. He et al., āLP perturbation solution of nonlinear free vibration of prestressed orthotropic membrane in large amplitude,ā Mathematical Problems in Engineering, vol. 2010, Article ID 561364, 17 pages, 2010. View at: Publisher Site  Google Scholar
 D. Li, Z. L. Zheng, and M. D. Todd, āNonlinear vibration of orthotropic rectangular membrane structures including modal coupling,ā Journal of Applied MechanicsTransactions of the ASME, vol. 85, no. 6, Article ID 061004, 2018. View at: Publisher Site  Google Scholar
 P. B. GonĆ§alves, R. M. Soares, and D. Pamplona, āNonlinear vibrations of a radially stretched circular hyperelastic membrane,ā Journal of Sound and Vibration, vol. 327, no. 12, pp. 231ā248, 2009. View at: Publisher Site  Google Scholar
 Z. L. Zheng, W. J. Song, C. J. Liu, X. T. He, J. Y. Sun, and Y. P. Xu, āStudy on dynamic response of rectangular orthotropic membranes under impact loading,ā Journal of Adhesion Science and Technology, vol. 26, no. 1011, pp. 1467ā1479, 2012. View at: Google Scholar
 C. J. Liu, Z. L. Zheng, J. Long, J. J. Guo, and K. Wu, āDynamic analysis for nonlinear vibration of prestressed orthotropic membranes with viscous damping,ā International Journal of Structural Stability and Dynamics, vol. 13, no. 2, Article ID 1350018, 2013. View at: Publisher Site  Google Scholar
 Z. L. Zheng, C. Y. Liu, D. Li, and T. Zhang, āDynamic response of orthotropic membrane structure under impact load based on multiple scale perturbation method,ā Latin American Journal of Solids and Structures, vol. 14, no. 8, pp. 1490ā1505, 2017. View at: Publisher Site  Google Scholar
 Z. L. Zheng, C. J. Liu, W. C. Gong et al., āA new methodāejection method for nondestructive online monitoring of the pretension of building membrane structure,ā Structural Control and Health Monitoring, vol. 20, no. 4, pp. 445ā464, 2013. View at: Publisher Site  Google Scholar
 C. J. Liu, Z. L. Zheng, X. Y. Yang, and H. Zhao, āNonlinear damped vibration of prestressed orthotropic membrane structure under impact loading,ā International Journal of Structural Stability and Dynamics, vol. 14, no. 1, Article ID 1350055, 2014. View at: Publisher Site  Google Scholar
 C. Liu, Z. Zheng, and X. Yang, āAnalytical and numerical studies on the nonlinear dynamic response of orthotropic membranes under impact load,ā Earthquake Engineering and Engineering Vibration, vol. 15, no. 4, pp. 657ā672, 2016. View at: Publisher Site  Google Scholar
 Z. L. Zheng, F. M. Lu, X. T. He, J. Y. Sun, C. X. Xie, and C. He, āLarge displacement analysis of rectangular orthotropic membranes under stochastic impact loading,ā International Journal of Structural Stability and Dynamics, vol. 16, no. 1, Article ID 1640007, 2016. View at: Publisher Site  Google Scholar
 D. Li, Z. L. Zheng, R. Yang, and P. Zhang, āAnalytical solutions for stochastic vibration of orthotropic membrane under random impact load,ā Materials, vol. 11, no. 7, p. 1231, 2018. View at: Publisher Site  Google Scholar
 D. Li, Z. Zheng, Y. Tian, J. Sun, X. He, and Y. Lu, āStochastic nonlinear vibration and reliability of orthotropic membrane structure under impact load,ā ThinWalled Structures, vol. 119, pp. 247ā255, 2017. View at: Publisher Site  Google Scholar
 D. Li, Z.L. Zheng, C.Y. Liu et al., āDynamic response of rectangular prestressed membrane subjected to uniform impact load,ā Archives of Civil and Mechanical Engineering, vol. 17, no. 3, pp. 586ā598, 2017. View at: Publisher Site  Google Scholar
 Q. S. Yang and R. X. Liu, āOn aerodynamic stability of membrane structures,ā Engineering Mechanics, vol. 23, no. 9, pp. 18ā29, 2006. View at: Google Scholar
 Q. X. Li and B. N. Sun, āWindinduced aerodynamic instability for closed membrane roofs,ā Journal of Vibration Engineering, vol. 19, pp. 346ā353, 2006. View at: Google Scholar
 M. Lazzari, M. Majowiecki, R. V. Vitaliani, and A. V. Saetta, āNonlinear F.E. analysis of Montreal Olympic Stadium roof under natural loading conditions,ā Engineering Structures, vol. 31, no. 1, pp. 16ā31, 2009. View at: Publisher Site  Google Scholar
 F. Rizzo, P. DāAsdia, M. Lazzari, and L. Procino, āWind action evaluation on tension roofs of hyperbolic paraboloid shape,ā Engineering Structures, vol. 33, no. 2, pp. 445ā461, 2011. View at: Publisher Site  Google Scholar
 F. Rizzo and V. Sepe, āStatic loads to simulate dynamic effects of wind on hyperbolic paraboloid roofs with square plan,ā Journal of Wind Engineering and Industrial Aerodynamics, vol. 137, pp. 46ā57, 2015. View at: Publisher Site  Google Scholar
 Y. Wu, Z.Q. Chen, and X.Y. Sun, āResearch on the windinduced aeroelastic response of closedtype saddleshaped tensioned membrane models,ā Journal of Zhejiang UniversityScience A, vol. 16, no. 8, pp. 656ā668, 2015. View at: Publisher Site  Google Scholar
 Y. P. Xu, Z. L. Zheng, C. J. Liu, K. Wu, and W. J. Song, āAerodynamic stability analysis of geometrically nonlinear orthotropic membrane structure with hyperbolic paraboloid in sag direction,ā Wind and Structures, vol. 26, pp. 355ā367, 2018. View at: Google Scholar
 Y.P. Xu, Z.L. Zheng, C.J. Liu, W.J. Song, and J. Long, āAerodynamic stability analysis of geometrically nonlinear orthotropic membrane structure with hyperbolic paraboloid,ā Journal of Engineering Mechanics, vol. 137, no. 11, pp. 759ā768, 2011. View at: Publisher Site  Google Scholar
 C. Liu, X. Deng, and Z. Zheng, āNonlinear windinduced aerodynamic stability of orthotropic saddle membrane structures,ā Journal of Wind Engineering and Industrial Aerodynamics, vol. 164, pp. 119ā127, 2017. View at: Publisher Site  Google Scholar
 W. Cui, L. Zhang, and H. Zhang, āThe saddle membrane structure response analysis under the simultaneous actions of wind and rain based on multiphase flow theory,ā Industrial Construction, vol. 48, pp. 60ā64, 2018. View at: Google Scholar
 H. ChiaLiu, āApproximate formulas of hailstone drag coefficients and its terminal velocity,ā Journal of Lanzhou University, vol. 1, pp. 90ā103, 1978. View at: Google Scholar
 S. C. Her and Y. C. Liang, āThe finite element analysis of composite laminates and shell structures subjected to low velocity impact,ā Composite Structures, vol. 66, no. 1ā4, pp. 277ā285, 2004. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2019 Changjiang Liu 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.