Improved Damping Constant of Hertz-Damp Model for Pounding between Structures
Earthquake-induced pounding between structures has been widely studied by using different impact models. In this paper, we first briefly introduced the existing impact analysis models and discussed their shortcomings. For the Hertz-damp impact analysis model, the damping constants do not accord with the physical facts; then, based on the Hertz theory, the expression of the damping constant in the Hertz-damp impact analysis model is improved. The approximate formula of the damping constant is theoretically derived, and the accuracy and effectiveness of the calculation formula of the damping constant are verified by the selected collision analysis problems.
Structural pounding is a complicated nonlinear problem. In order to effectively study the response of the adjacent structures under the impact, the reasonable impact analysis model is very significant. Quite a lot of scholars make a series of studies on the impact models. Since the Kelvin model is frequently used to simulate the collision of adjacent structures, the expression of the damping coefficient in Kelvin model has been derived by Ye et al. . In order to preferably simulate the force and displacement relationship in the process of simulation, many scholars have used the Hertz theory to study the impact response. Jankowski  puts forward a nonlinear viscoelastic model based on the Hertz impact model and then derived the expression of the damping constant of the nonlinear viscoelastic impact model. The damping constant is related to the recovery coefficient, and the analytical formulation is verified by numerical simulations. Because the Hertz model does not consider the energy loss in the process of collision, Muthukumar and DesRoches  introduced the Hertz-damp model to be used in the impact analysis of structural engineering, which is compared with other impact models. And a lot of scholars [4, 5] have conducted in-depth research on the structural collision.
In this paper, based on the review of the main impact models and the Hertz theory, a new approximate formula of the damping constant of Hertz-damp model is derived. Meanwhile, the numerical simulation of the pounding analysis problems is conducted by MATLAB/Simulink software, the results of impact simulation show the correctness, accuracy, and effectiveness of the approximate formula of the damping constant of Hertz-damp model.
2. Review of Existing Impact Models
At present, the classical dynamic method and the contact element method are two main methods for the study of structural impact. As the classical dynamics method for the numerical analysis of the impact is not very convenient, generally rarely used, the second method is through establishment of the mathematical model of impact force for structure impact analysis, adding contact element at the impact site, which uses equivalent spring element and equivalent damper elements to, respectively, simulate interaction and energy dissipation during structural collision. Therefore, it is widely used. The following will be a brief review about the existing collision analysis models.
2.1. Linear Elastic Model
The linear elastic model is the simplest kind of collision analysis model; the mathematical expression of the impact force is where is the impact spring stiffness and stands for the relative displacement of two colliding objects. This model does not consider the energy loss in the process of pounding, and the impact force is proportional to relative displacement .
2.2. Kelvin Model
The Kelvin model is most widely used, and the linear viscous damping of Kelvin model is to represent the loss of energy. The impact force is expressed aswhere is the relative displacement of two colliding objects, expresses the relative deformation velocity of the colliding objects during the impact process, is the impact spring stiffness of Kelvin model, and stands for the damping coefficient, which is determined by the following formula:where and are the qualities of the colliding objects and stands for the damping constant, which is related to recovery coefficient , and its expression is given by the following formula:
The shortcoming of linear elastic impact models is that the viscous damping coefficient is kept constant during the whole collision process, which leads to a uniform change in energy loss at the impact start and rebound process; meanwhile, the Kelvin model will cause the impact force to show the pulling force at the impact beginning stage and emerge as leap at the impact rebound stage, and these phenomena are not consistent with the actual impact.
2.3. Hertz Model
In the impact process, the increase of the impact force simulated by Hertz model is nonlinear. In order to simulate the impact process more accurately, many scholars have studied the Hertz theory. The impact force based on Hertz theory is expressed aswhere stands for the relative displacement of two colliding objects and is the impact spring stiffness of Hertz model. However, Hertz model does not consider the energy loss in the impact process.
2.4. Hertz-Damp Model
The Hertz-damp model makes up the deficiency of the Hertz model without considering the energy loss and adds the nonlinear damping term. The relationship between the impact relative invasion displacement and the impact force is expressed aswhere is the relative displacement of two colliding objects, expresses the relative deformation velocity of the colliding objects during the impact process, is the Hertz-damp model impact spring stiffness, and stands for the damping coefficient, and its expression is determined by the following formula: The relationship between damping constant , nonlinear stiffness , coefficient of restitution , and the relative initial velocity of structure impact is expressed as
In the nonlinear viscoelastic model, indicates that the impact is completely elastic impact and indicates that the impact is completely plastic impact . Through (8), it can be found that and fits perfectly elastic impact; however and does not conform to the full plastic impact. Therefore, there are errors in the pounding simulation results simulated by the Hertz-damp model.
2.5. Nonlinear Viscoelastic Model
In the nonlinear viscoelastic model, Jankowski divided the collision process into the approach period and the restitution period. He thought that the cracking, crushing, and plastic deformation in the process of collision mainly occurred in the approach period; therefore, the energy loss mainly occurred in this stage. According to this, the Hertz-damp model was improved by Jankowski . In the nonlinear viscoelastic model, the formula for the pounding force is expressed as follows:where stands for the relative displacement of two colliding objects, is the relative deformation velocity of the colliding objects during the impact process, expresses the impact spring stiffness of nonlinear viscoelastic model, and is the damping coefficient, its expression is determined by the following formula: The relationship between damping constant and coefficient of restitution is expressed as
However, for the nonlinear viscoelastic model, the time history curve of the impact force simulated by the nonlinear viscoelastic model does not vary smoothly between the approach phase and restitution period of the collision , which is because there is a defect hypothesis in the process of velocity approximation formula derivation. All speed changes are added to the approach phase, which influences the coherence between the approach phase and restitution period of the collision. Therefore, the damping constant needs further optimization.
3. Correction of Damping Constant
In analysis model of structural pounding as shown in Figure 1, and are the qualities of two colliding spheres; and stand for the displacements of two colliding spheres; and express respective velocities of two colliding spheres at the initial moment of contact; and are the velocities of colliding spheres after impact; stands for the damping coefficient; is the initial clearance of the colliding spheres; expresses the stiffness coefficient of impact contact element.
According to classical mechanics, energy loss caused by the collision can be expressed as a function of coefficient of restitution and the relative velocity of two colliding objects :According to the classical theory of impact, coefficient of restitution can be calculated from the following equation:At the same time, the energy loss dissipated by the damper can be expressed as
The variation of interpenetration displacement during impact is illustrated in Figure 2, where , , and stand for the initial time of contact, the time which maximum interpenetration displacement occurs and the time of separation between the two colliding bodies.
The two degrees of freedom system shown in Figure 3 can be equivalent to a single degree of freedom system, which can be seen as a single mass, spring, and damper model, as shown in Figure 3. The definition of equivalent model is as follows: () the spring stiffness and damping are the same as that of two degrees of freedom system; () the mass of the equivalent model is equivalent masses of the two colliding bodies: ; () in Figure 3, A represents the initial point and the initial interpenetration displacement is , and the initial interpenetration velocity is , when the maximum displacement point is reached, , , among them, corresponds to the maximum interpenetration displacement in the collision process. The approach phase is represented from A to the maximum displacement point, the same way from the maximum displacement point to A indicates the restitution phase. In this single degree of freedom system, the equation of motion of the equivalent model is derived aswhere stands for the equivalent mass of the single degree of freedom system, when , it represents Hertz-damp model. When , the analytical solution of the equation of motion (15) will be very difficult to be found; however, it can be solved by a simple iterative method.
Using to replace , thenSubstituting (16) into (15), we can getEmploying the traditional separation of variables, ignoring the influence of the damping based on the Hertz contact theory, and integrating (17) yieldwhere is the interpenetration velocity, stands for the interpenetration velocity at the moment of collision separation, and expresses the interpenetration displacement.
After integral calculation, one obtainsTransformation form of (19) and then the displacement deformation can be expressed asMaximum interpenetration displacement can be expressed asBased upon (20) and (21), the relationship between and can be get:From (22), the interpenetration velocity at the restitution period can be obtained:where indicates the interpenetration velocity of two objects at the moment of separation.
Substituting (23) into (14), the energy loss caused by the damping force can be expressed asThe evaluation of the integral of (24) can be found elsewhere . Applying this method to solve the following equation,Equating (25) and (12) leads toSubstituting (21) into (26), one obtainsThen, the expression of improved damping constant can be obtained:
4. Numerical Verification
In order to verify the correctness of the theoretical derivation for the expression of damping constant , the following two selected collision analysis problems are selected for numerical verification. The procedures for the verification are as follows: first, selecting the collision analysis problem and determining the initial value of the coefficient of restitution , which varies from 0.1 to 1.0 with interval 0.1; then, obtaining the numerical solution of the coefficient of restitution by the collision analysis based on the impact models; and finally comparing relative error  between the initial value of the coefficient of restitution and the numerical solution of the coefficient of restitution .
4.1. The Collision between Sphere and Static Rigid Plane
The collision between the steel ball and the rigid surface by the free-falling  is shown in Figure 4. The following parameters have been used: , , and . The dynamic equation of motion for such a model can be obtained by analyzing the form as follows:where is the quality of falling ball, stands for the vertical acceleration of falling ball, expresses its vertical displacement, indicates the acceleration of gravity, is the drop height of falling ball, impact force is calculated according to (6) and (9), and impact deformation is counted as
The previously described nonlinear viscoelastic impact model together with the Hertz-damp model with the improved damping constant is used to simulate this collision analysis problem, which has been conducted for impacts analysis of a ball falling on a rigid surface applying various impact models . In this paper, a set of data, the maximum impact force of 80.7 N in the test, is selected for numerical simulation, and the appropriate values of impact parameters (, ) have been obtained numerically through iterative simulations in order to satisfy the maximum impact force. A time-stepping integration procedure written through MATLAB/Simulink software with the fixed time increment step s has been applied to solve (29), the impact force time history curves of the nonlinear viscoelastic impact and the Hertz-damp model with the improved damping constant are presented in Figure 5.
4.2. The Collision between the Pendulum and the Static Rigid Wall
In order to further verify the correctness of the approximate formula of damping constant , the collision between the pendulum and the static rigid wall  is as shown in Figure 6. The following parameters have been used: kg, m, and and the initial impact velocity is . The motion dynamic equation of this model can be expressed as follows:where is the mass of the pendulum striker, and , respectively, express the horizontal acceleration and horizontal displacement of the pendulum striker, stands for the length of the pendulum striker, expresses the acceleration of gravity, and stands for the impact force, which is calculated, respectively, according to (6) and (9), in which the deformation is expressed as
The numerical analysis has been conducted for impacts of the pendulum and the static rigid wall applying the nonlinear viscoelastic impact and the Hertz-damp model with the improved damping constant. In this paper, the maximum impact force of 102.5 kN in the test  is selected for numerical analysis, and the appropriate values of impact parameters (, ) have been determined through iterative procedures in order to attain the maximum impact force. A time-stepping integration procedure written by MATLAB/Simulink software with constant time step s has been used for solving (31) numerically and the time history curves of the impact force of different impact models are shown in Figure 7.
4.3. Comparison between Analytical Formulations and Numerical Results
The analytical formulation, which is related to the impact damping ratio and the coefficient of restitution , has been derived as shown in (28), the comparison with the numerical results has been carried out as shown in Figure 8 for verifying the accuracy and effectiveness of the expression. As shown above, the numerical studies have been conducted for pounding with nonlinear viscoelastic model and modified Hertz-damp model.
Simulation results of different impact models and the relative error percentages of recovery coefficient at different initial values are presented in Tables 1 and 2. As can be seen from the tables, the impact model of this paper has higher simulation precision than the nonlinear viscoelastic model and Hertz-damp model. At the same time, the relative error decreases with the increase of initial value , and when the initial value of the coefficient of restitution is greater than 0.5, the relative error is less than 1.5%.
Contrasting pounding force histories presented in Figures 5 and 7, it can be seen that the time history curves of the impact force simulated by the modified Hertz-damp model vary smoothly; however, impact force curves simulated by the nonlinear viscoelastic model do not vary consistently at a certain moment, which is the same as the simulation results by Jankowski ; at the same time, the simulation results obtained by modified Hertz-damp are closer to actual value than that simulated by nonlinear viscoelastic model at approach phase, and the simulation results by these two models are relatively close during the restitution period.
Through Figure 8, it can be seen that two analytical formulations can be used to express well the relationship between the impact damping ratio and the coefficient of restitution ; however, by carefully observing, we can find that analytical formulation 2 can quite accurately simulate the relation between the impact damping ratio and the coefficient of restitution between 0.5 and 0.75. As we all know, in the collision of the actual engineering structure, it has been assessed that the coefficient of restitution is generally used to simulate real collisions between structures and ranges usually from 0.5 to 0.75 ; thus it can be concluded that reliable results of impact simulation in structural engineering can be provided by using the improved damping constant formula.
This paper first reviews the existing structure impact analysis models and discusses the problems between them; then, the damping constant of the Hertz-damp model is theoretically rededuced according to Hertz theory, and the accuracy and effectiveness of the damping constant are verified by numerical simulation. In the practical structure engineering, more reliable results can be obtained by using the Hertz-damp model with the improved damping constant.
The authors declare no competing interests regarding the publication of this paper.
This research was supported by National Natural Science Foundation of China under Grant no. 51308487, no. 41402261, and no. 51408526, Heibei Provincial Natural Science Foundation of China under Grant no. E2014203055, and the Excellent Youth Foundation of Hebei Educational Committee under Grant no. YQ2013015. The support for this research is greatly appreciated.
P. Liu, K.-F. Zheng, W. He, and H. Liu, “Improved algorithm of contact element model for earthquake-induced pounding between structures,” Journal of Chongqing University, vol. 35, no. 2, pp. 135–142, 2012.View at: Google Scholar
C. H. Edwards and D. E. Penney, Calculus and Analytic Geometry, Prentice Hall, Englewood Cliffs, NJ, USA, 1st edition, 1982.
X. Xiang, L. Weiqing, and X. Xiuli, “Linear viscoelastic model of earthquake-induced structural pounding,” Engineering Mechanics, vol. 30, no. 2, pp. 278–284, 2013.View at: Google Scholar