[转帖]碰撞动力学模型综述

MVH

西北工业大学工程力学系 鲍四元(翻译并添加整理)

——科研中国http://www.SciEi.com 特邀论文

摘要：本文目的是展现撞击分析的总体回顾和此领域内的一些重要方法。

1 撞击理论的模型

含动能约束的多体系统的动态分析是已经完善的力学分支。为了建立数学模型，物体都被假设成为刚性，且铰接处认为不含间隙。

撞击问题吸引着从天体物理学到机器人学等不同学科领域学者的注意力。他们的共同目标是发展能够预测撞击物行为的理论。本文主要集中于与刚体有关的撞击模型。

撞击理论的演化主要含有四个方面：经典力学、弹性应力波传播、接触力学和塑性变形。不同的撞击理论适用于不同撞击特性(速度和材料性质)、假设和相关结论。

(1) 经典力学

包含应用基本力学定理来预测撞击后的速度。脉冲-动量定理构成这种方法的核心。Goldsmith在著作[1]中用了一章的篇幅介绍了这种方法在几个问题中的应用。Brach[2]在模拟几个具有实用价值的问题时一律采用了此法。这种方法具有简便和易于实现的特点。实际问题中的能量损失是通过恢复系数实现的。然而，此法不能预报物体之间的接触力和物体的应力。

(2)弹性应力波传播

撞击通过以撞击点为起点，应力波在撞击物之间的传播描述。总能量中的一部分转化为振动，这样，经典理论就无法验证这种理论。Goldsmith把这种方法应用于如下问题中：两杆的纵向碰撞、质点和杆碰撞、粘弹性对碰撞的影响等。Zukas等[3]也广泛地应用了这一方法。波传播法用来研究细长杆的纵向碰撞问题。近年文献[4,5]使用符合运算软件给出两类典型问题：质点杆撞击和杆撞击地面问题的符合表达式解。文献研究了[6]平面波在含空洞材料中的传播与考虑径向剪力和惯性力时波在圆柱形杆中传播具有模拟关系。文献[7]于不对称粘弹性杆在频域的波传播解，给出了理论和实验分析。

(3)接触力学

两个物体撞击产生的接触应力是碰撞研究中的另一个研究热点。常规接触力学主要与静态接触有关，尽管此法在涉及撞击时已经延伸至近似解。对于球形接触面，Hertz理论常被用于撞击关系的获得，从而计算撞击时间和最大变形。此方法还被用于含塑性变形的情况。通常假设材料有一个屈服点。当Hertz理论不适用时，也可使用屈服区模型。撞击力变形关系常通过增加一个阻尼项来反映接触区域的能量耗散，从而允许把接触区作为一个弹簧-阻尼系统的模型。

(4) 塑性变形

当塑性应变超过容许变形时，弹性波模型不再适用于分析撞击问题。这类问题属于高速撞击问题，如发生爆炸和侵彻时。Goldsmith[1]提供了2种方法：水动力学理论和塑性波传播理论。水动力学理论中，假设物体密度发生变化，材料的状态方程于密度、温度的变化相关，同时利用了能量、动量和质量守恒定理。而塑性波传播理论中，塑性区的材料认为是不可压缩的。同样，与应变、应力、应变率有关的状态方程假设与温度无关。Maugin[8]和Lubliner[9]假设了脆性材料，荷载的加载是一个长时间的过程。Zukas[3]提供了分别使用应变相关和应变独立理论的塑性波传播理论。文献[10]考虑了梁与梁碰撞的问题，采用了质量-弹簧模型。梁之间的能量能够很好地近似刚塑性解。

工程师常需要解答如下2个基本问题：(1)撞击前后速度变化的关系。(2)撞击点的碰撞力多少?

当恢复系数给定时，脉冲-动量定理 方法能够回答第一个问题。但前面已经提到，此法不能确定撞击力，即解决不了第二个问题。波传播理论可以得到撞击物内的应力，但动力分析中的积分比较复杂。接触力学方法把接触区域作为弹簧-阻尼系统，使撞击问题作为连续时间动力问题处理。塑性大应变理论在解决弹道学领域中的爆炸、侵彻时最有效。但本文不涉及这方面中高速碰撞问题。

2 关于恢复系数的历史与现状

根据Kozlov[11]，关于撞击的首次研究可追溯道1668年，由Wallis, Wren和Huygens进行。Netwon后来于1687年在他的著作《Mathematical Foundations of Natural Philosophy》中参考了Wren的工作。Huygens的工作成果是推导出了动量守恒定理，从而成为撞击理论的基础。这个理论的主要假设是认为物体是刚性的，因此撞击持续时间为0。单独使用动量守恒定理不足以确定撞击后撞击物和靶体各自的速度。因此初等撞击理论考虑了两种极限情况：完全弹性碰撞和完全非弹性碰撞。完全弹性碰撞指碰撞前后系统的动能守恒。而完全非弹性碰撞指撞击后撞击物和靶体连为一体共同运动，从而组合体的速度可以通过定理守恒定理确定。然而，通常的撞击既不是完全弹性碰撞，也不是完全非弹性碰撞。初始动能的损失是通过恢复系数e的引入(Netwon提出这一观点)来实现的。

v1f -v2f=-e(v1i-v2i)


其中下标1和2分别表示撞击物和靶体，而i和f分别表示初始(initial)状态和最终(final)状态。e是个无量纲的系数，其值介于0和1之间，0对应于完全弹性状态，1对应于完全非弹性状态。恢复系数的一个对能量损失的综合概念，可包括不同的能量损失，如材料的粘弹性、接触面的塑性变形和两个物体之间的振动等。


恢复系数不是仅仅依赖于材料的一种固有属性，它取决于撞击物和靶体的材料、接触面的几何性质和撞击速度[1，p.262]。近年来，文献[12]使用能量法研究细长杆与光滑界面碰撞的恢复系数，提出了影响恢复系数的2个因素：碰撞倾斜解和反映杆几乎和材料性质的常数Hr。使用恢复系数的优点在于数学表达上的简洁性。姚文莉[13]使用波传播理论分别提出质点与杆和梁碰撞的恢复系数的求法。得到损失波动能量在质点-杆碰撞问题中所占比例的数学表达式。Brach在文献[2]中广泛使用了恢复系数来解决撞击问题。Brach还注意到恢复系数可取-1和0之间的数。这表明在撞击过程中损失了一些能量，但并不产生速度方向的改变。如侵彻物在穿过靶体时虽然降低了自身速度，但速度方向没有改变。


若干文献研究了撞击物初始速度与恢复系数的关系。靶体是粘弹性材料时，提出如下观点[14-16]：

e(v)=1-f(v1/5)

上式表明撞击速率越大，恢复系数就会变低。也即撞击物高速碰撞时，损失的能量更多。式(2)仅考虑粘弹性材料。现实中，还有其他的因素需要考虑。高速碰撞时，弹性波传播时的耗散及塑性变形消耗的能量需要考虑。而低速碰撞时的粘性力和重力显得更加重要。文献[17]中利用恢复系数讨论了粘弹性地基上的撞击响应问题。


3 接触力-变形模型


关于撞击力初级理论的上述综述基于完全刚体的简化假设。撞击物的实际情形是复杂的，并且撞击持续时间决对大于0。更为接近实际的模型是采用连续时间动态模型。这个方法的成功之处在于基于完善的数学模型。通常，接触力-变形关系如下：

F=Fc(δ)+Fv(δ, dδ/dt)+Fp(δ, dδ/dt)


Fc是接触力的弹性部分，Fv是粘弹性阻尼部分，Fp是由塑性变形导致的耗散部分。以下主要介绍接触力的弹性部分。其中1882年Hert关于半无限固体的弹性接触工作具有重大意义。Johnson[18]对此理论做了很好的介绍，并于附录中列举了相关公式。Hertz理论指出了应力在接触区的分布，也给出计算法向应力和剪切应力在撞击体内的分布。一个很常用的结论是球体-球体接触时的接触力-变形关系[18]：

F=Kδ3/2
　

其中 F 是撞击物和靶体之间压缩时的法向力


δ是两个球体之间的缩进，也即两个表面之间总的变形


K是取决于球体半径和材料弹性常数的常数。

MVH

　　

4 近年来的进展：

构方法研究了刚性小球和均匀柔性杆的纵向碰撞及和均匀柔性梁的横向碰撞问题，导出了用模态坐标表示的动力学方程。

(2) 直接模态叠加法研究弹性撞击问题 邢誉峰等利用DMSM策略，讨论了等截面杆、梁的碰撞问题[19-26]。文献[26]指出：这种方法可以得到结构弹性碰撞问题的解析解;这种方法不但可以用来分析平动结构的碰撞问题[19-25]，还可以用来分析机构的各种弹性锁定问题[22];不但可以用来分析结构的点碰撞问题[19-20]，对结构的线、面接触和碰撞等问题同样有效[23]。 对于梁碰撞问题，文献[24]进行了如下研究：考虑线弹性接触变形的前提下,分别对质点、杆与简支Euler-Bernoulli梁的垂直正撞问题进行了研究。文献[25]基于不同梁理论：Euler梁、Timoshenko梁和翘曲理论，比较了结构遭受冲击的动态响应。

　文献中，如果用一个假想的弹簧来模拟两个结构相碰处的接触刚度，并通过该弹簧把撞击体和靶体联系成一个组合振动体系，就可把结构碰撞分析转化为常规的结构振动响应分析问题，即是该组合振动体系在其撞击部分具有给定初始速度模式下的振动响应问题。因此可以方便地直接使用常规的振动模态叠加法或时间积分法来求解撞击问题。文献具体报道了利用解析模态和有限元离散模态求解质点-弹性杆的撞击力变化过程，并讨论了各种因素以及有限元建模对结果的影响。

　　(3)纤维复合板

　纤维复合板复合板受到低速撞击问题已被许多学者研究过。Sun和Chattopadhyay[27]研究了一个四边简支各向同性板受到中心撞击的情形，并考虑了横向剪切变形。Dobyns[28]研究了受均布矩形荷载时的撞击情形。A.Carvalho和C Guedes Soares[29]也研究了板的撞击响应，对位移、转角采用Fourier级数展开，数值积分用Nemark方法，并与拉普拉斯解进行了比较。

　(4)有限元方面的进展

　文献[30]较早使用有限元方法研究了接触/撞击问题。文献[31]使用辛方法研究了非线性撞击问题。Jerome M. Solberg, Panayiotis Papadopoulos [32]基于非线性力学有限元原理，使用数值方法研究了接触/撞击问题。对于无摩擦问题，建立数值微分方程。在接触面上损失了一部分能量，以稳定接触面的动能场。数值解采用了Nemark积分法，较好地模拟了接触/撞击过程。文献[33]依据波传播理论提出一种新的数值算法：含有模态综合的有限元计算法，并与柔性杆受轴向撞击的经典St. Venant解进行了比较。

　　台湾学者R.-F. FUNG AND J.-H. SUN和 J.-W. WU[34]研究了研究了滑动曲柄机构在撞击下的轨迹控制。Khulief and Shabana[35]通过GMB途径来研究多体系统的撞击问题，同时发展了CFM方法来研究多体系统撞击问题。

　　除了上述研究，近年来许多学者对不等截面杆及受载梁的自由振动进行了大量研究。Q.S. LI等对等截面杆、不等截面杆含有集中质量-弹簧耦合系统进行了大量研究[36-47]。

　M. Gürgöze针对两个固支-自由纵向振动杆，端部带有质量块，由两弹簧-质量系统耦合文献，还讨论了梁含有阻尼器的自由振动[48]。

　　5.本论文主要工作：

　　1给出细长圆锥形的截面杆受到质点纵向弹性碰撞时的精确解析解。使用了一种新方法用于分析质点-圆锥形杆碰撞，即由叠加法给出杆的响应。其结果可验证数值解和其他解析解。算例显示，所提出方法的优点之一是响应解的解析形式简洁。算例表明一些描述杆几何形状的变量在撞击分析中具有重要作用。

　　研究了含弹簧的锥形杆结构撞击问题的解析解。振动过程中把杆和质点作为整体考虑，采用无量纲变量，从而简化方程模型。算例说明了杆中波传播情况和撞击端的响应，并且讨论了质量比和锥形杆截面倾角对波传播的影响。解决了锥形杆结构的纵向撞击问题，并且与等截面杆的纵向撞击问题进行了比较。

　　2 把杆的质量函数和刚度函数作为2个独立的函数，对于质量函数和刚度函数的若干种形式，进行适当的函数变换后，基本方程转化为Bessel方程或具有常系数的常微分方程。得到满足正交条件的基本解，并且建立了一阶非均匀杆碰撞时的频率方程。

　　研究了如下受载梁的撞击问题：一个一端固支，一端自由的杆含有一个弹簧-质量耦合系统。使用DMSM方法，撞击问题转化为系统具有初速度的轴向振动系统。把撞击物和靶体作为一个整体振动，获得系统的微分方程，分离时间变量后，化为有初始边值的常微分方程问题。考虑系统的一些参数对系统频率的影响，并且给出此杆结构受撞击后的动态响应。

　　3 对于质点-梁撞击问题，把撞击物和被撞击物分开考虑，引入撞击力-时间模型，得到如下两种预报撞击力的方法：第一：基于位移协调方程的解法;第二：基于动力微分方程的数值法。与文[76]相比，方法一大幅简化了计算过程，得到近似解，且此法可以推广到四边简支板中去。

　　4对于复合材料梁端部受撞击问题，本文把质量块看成质点，使用模态叠加法提供了弯扭耦合作用下的分析方法。算例表明此方法是有效的。

　　5 提出一种确定恢复系数的方法：即首先使用DMSM方法得到撞击结束时间，再得到恢复系数的步骤。算例表明，本文方法能够从理论上得到弹性碰撞恢复系数的表达式，且结果是有效的。

MVH

参考文献：

1. Goldsmith, Werner, Impact, Edward Arnold Publishers, London, 1960.

2. Brach, Raymond M., Mechanical Impact Dynamics: Rigid Body Collisions. JohnWiley & Sons, New York, 1991.

3. Zukas, Jonas A.; Nicholas, T.; Swift, H. F.; Greszczuk, L. B.; Curran, D. R. Impact Dynamics, Krieger Publishing Company, Malabar, FL, 1992.

4 B. Hu, P. Eberhard and W. Schiehlen. Symbolic Impact Analysis for a Falling Conical Rod against the Rigid Ground [J].Journal of Sound and Vibration. 2001, 240(1): 41-57.

5 Bin Hu, Peter Eberhard. Symbolic computation of longitudinal impact wave [J]. Comput. Methods Appl. Mech. Engrg. 2001, 190: 4805-4815.

6. THOMAS W. WRIGHT. ElASTIC WAVE PROPAGATION THROUGH A MATERIAL WITH VOIDS. Journal of the Mechanics and Physics of Solids, 1998, 46(10): 2033-2047.

7. A. Benatara, D. Rittelb, A.L. Yarinb. Theoretical and experimental analysis of longitudinal wave propagation in cylindrical viscoelastic rods. Journal of the Mechanics and Physics of Solids 2003(51) 1413-1431.

8. Maugin, Gerard A., The Thermomechanics of Plasticity and Fracture, Cambridge University Press, Cambridge, 1992.

9. Lubliner, Jacob, Plasticity Theory, Macmillan Publishing Company, New York, 1990.

10. H. H. Ruan and T. X. Yu. Collision between mass–spring systems. International Journal of Impact Engineering, 2005, 31(3): 267-288.

11. Kozlov, Valerii V. and Treshchëv, Dmitrii V., Billiards: A Genetic Introduction to the Dynamics of Systems with Impacts, American Mathematical Society, 1991.

12. Yildirim Hurmuzlu. ASME Journal of Applied Mechanics, Vol. 65, No.4, An Energy Based Coefficient of Restitution for Planar Impacts of Slender Bars with Massive External Surfaces: 952-962.

13 姚文莉. 考虑波动效应的碰撞恢复系数研究. 山东科技大学学报，2004, 23 (2)：83-86.

14. Ramírez, Rosa; Pöschel Thorsten; Brilliantov Nikolai V. and Schwager Thomas, Coefficient of restitution of colliding viscoelastic spheres, Physical Review E, 1999 60 (4), 4465-4472.

15. Falcon, E.; Laroche C.; Fauve S. and Coste C., Behavior of One Elastic Ball Bouncing Repeatedly off the Ground, The European Physical Journal B, Vol. 3, 1998, 45-57.

16. Kuwabara, G. and Kono K., Restitution Coefficient in a Collision Between Two Spheres, Jap. J. of Appl. Physics, 1987, 26(8): 1230-1233.

17. 尹邦信. 弹性板受撞击的动力响应分析. 应用数学和力学, 1996, 17 (7)：639-644.

18. Johnson K. L., Contact Mechanics, Cambridge University Press, 1985.

19 Xing YuFeng, Zhu DeChao, Analytical Solutions of Impact Problems of Rod Strucutures with Springs. Comput. Methods Appl. Mech. Engrg. , 1998,160: 315～323

20 诸德超，邢誉峰. 点弹性碰撞问题之解析解. 力学学报, 1996, 28( 1): 99～103.

21 Xing Yu-feng, Qiao Yuan-song, Zhu De-chao, Sun Guo-jiang. Elastic Impact on Finite Timoshenko Beam. ACTA Mechanica Sinica. 2002, 18(3): 252-263.

22邢誉峰，诸德超. 杆和板弹性正碰撞的瞬态响应. 航空学报, 1996, 17(7): S42～S46

23邢誉峰. 梁结构线弹性碰撞的解析解. 北京航空航天大学学报. 1998, 24 (6): 633-637.

24邢誉峰.有限长Timoshenko梁弹性碰撞接触瞬间的动态特性[J].力学学报,1999,31(1):67-74.

25邢誉峰，诸德超，乔元松. 复合材料叠层梁和金属梁的固有振动特性. 力学学报. 1998,30 (5): 628-634.

26 邢誉峰，诸德超. 杆和梁在锁定过程的响应. 计算力学学报, 1998, 15(2): 192～196.

27 C.T. Sun and S. Chattopadhyay, Dynamic response of anisotropic laminated plates under initial stress to impact of a mass. J. Appl. Mech. 42 (1975), p. 693

28. A.L. Dobyns, Analysis of simply-supported orthotropic plates subject to static and dynamic loads. AIAA J. 19 (1980), 642.

29 A. Carvalho and C. Gliedes Soares. Dynamic response of rectangular plates of composite materials subjected to impact loads. Composite Structures, 1996,34(1): 55-63.　

30. T.J.R. Hughes, R.L. Taylor, J.L. Sackman et al., A finite element method for a class of contact-impact problems, Comput. Methods. Appl. Mech. Engrg. 8 (1976) 249-276.

31. J.C. Simo, N. Tarnow, K.K. Wong, Exact energy-momentum conserving algorithms and symplectic schemes for nonlinear dynamics, Comput. Methods. Appl. Mech. Engrg. 100 (1992) 63-116.　

32. Jerome M. Solberg, Panayiotis Papadopoulos. A finite element method for contact/impact.Finite Elements in Analysis and Design. 1998 (30) 297-311.　

33 J.L. Escalona, J. Mayo, J. DomõÂ nguez. A new numerical method for the dynamic analysis of impact loads in flexible beams. Mechanism and Machine Theory. 1999(34) 765-780.　

34. R.-F. FUNG, J.-H. SUN, J.-W. WU. TRACKING CONTROL OF THE FLEXIBLE SLIDER CRANK MECHANISM SYSTEM UNDER IMPACT Journal of Sound and Vibration. 2002. 255(2), 337-355.　

35 Y. A. KHULIEF, A. A. SHABANA. Dynamic analysis of constrained system of rigid and flexible bodies with intermittent motion. American Society of Mechanical Engineers Journal of Mechanisms, Transactions, and Automation in Design 1986(108), 38–45.

36. Q. S. Li, L. F. Yang, Y. L. Zhao and G. Q. Li. Dynamic analysis of non-uniform beams and plates by finite elements with generalized degrees of freedom. International Journal of Mechanical Sciences, 2003, 45(5): 813-830.　

37. Ke Yang , Q. S. Li and Lixiang Zhang. Longitudinal vibration analysis of multi-span liquid-filled pipelines with rigid constraints. Journal of Sound and Vibration, 2004, 273, (1,2):125-147.

38. Q. S. Li. Analytical solutions for buckling of multi-step non-uniform columns with arbitrary distribution of flexural stiffness or axial distributed loading. International Journal of Mechanical Sciences, 2001, 43(2): 349-366.

39. Q. S. Li. Non-conservative stability of multi-step non-uniform columns. International Journal of Solids and Structures, 2002, 39(9): 2387-2399.

40. Q. S. Li. Buckling of multi-step non-uniform beams with elastically restrained boundary conditions. Journal of Constructional Steel Research, 2001, 57(7): 753-777.

41. Q. S. Li. Classes of exact solutions for buckling of multi-step non-uniform columns with an arbitrary number of cracks subjected to concentrated and distributed axial loads. International Journal of Engineering Science, March 2003, , 41(6): 569-586.

42. Q. S. Li. Torsional vibration of multi-step non-uniform rods with various concentrated elements. Journal of Sound and Vibration, 2003, 2607(4): 637-651.　

43. H. Qiao, Q. S. Li and G. Q. Li. Vibratory characteristics of flexural non-uniform Euler–Bernoulli beams carrying an arbitrary number of spring–mass systems. International Journal of Mechanical Sciences, 2002, 44(4): 725-743.

44. Q. S. Li, J. R. Wu and Jiayun Xu. Longitudinal vibration of multi-step non-uniform structures with lumped masses and spring supports. Applied Acoustics, 2002, 63(3):333-350.

45. Q. S. Li. Buckling analysis of non-uniform bars with rotational and translational springs. Engineering Structures, 2003, 25(10):1289-1299.

46. Q. S. Li. FREE VIBRATION ANALYSIS OF NON-UNIFORM BEAMS WITH AN ARBITRARY NUMBER OF CRACKS AND CONCENTRATED MASSES. Journal of Sound and Vibration, 2002, 252(3): 509-525.

47. Q.S. Li, G.Q. Li, D.K. Liu. International Journal of Mechanical Sciences. 2000, (42) 1135-1152 Exact solutions for longitudinal vibration of rods coupled by translational springs.

48 M. Gürgöze. ON THE EIGENVALUES OF VISCOUSLY DAMPED BEAMS, CARRYING HEAVY MASSES AND RESTRAINED BY LINEAR AND TORSIONAL SPRINGS. Journal of Sound and Vibration, 1997208(1): 153-158.　　

前半部分翻译自文献：

Salah Faik, Holly Witteman. Modeling of Impact Dynamics: A Literature Survey. International ADAMS User Conference. 2000. 1-11.

并结合自己的阅读写成。

yejet

这篇不错，值得大家仔细看看

Anonymous

不错，收藏！

ChinaMaker

大家可以到http://www.sciei.com/bbs/dispbbs ... p;ID=215&page=2 交流

frogfish

《柔性多体系统接触碰撞子结构动力学模型》
摘要 对柔性多体系统接触碰撞的子结构动力学模型进行了深入的研究. 将柔性杆
纵向碰撞的解析解与该对象的子结构动力学模型的数值解进行了对照, 揭示了子结构
动力学模型能描述在接触碰撞期间应力波在柔性杆中的传播现象以及波的传播速度与
杆的某些参数间的关系, 研究了影响该模型精确性的参数范围, 给出了该参数的一些
优化值.

zzzxxxx

不错，就是比较难 啊

zylzl

不错啊，要好好借鉴借鉴！

ansys10

hao a ,学习了

cadpair

zswx1866

好多东西，学习了。