能量吸收装置耐撞性研究毕业论文

《能量吸收装置耐撞性研究毕业论文》由会员分享，可在线阅读，更多相关《能量吸收装置耐撞性研究毕业论文（54页珍藏版）》请在装配图网上搜索。

1、能量吸收装置耐撞性研究毕业论文毕业设计（论文）原创性声明和使用授权说明原创性声明本人郑重承诺：所呈交的毕业设计（论文），是我个人在指导教师的指导下进行的研究工作及取得的成果。尽我所知，除文中特别加以标注和致谢的地方外，不包含其他人或组织已经发表或公布过的研究成果，也不包含我为获得 及其它教育机构的学位或学历而使用过的材料。对本研究提供过帮助和做出过贡献的个人或集体，均已在文中作了明确的说明并表示了谢意。作 者 签 名： 日 期： 指导教师签名： 日期： 使用授权说明本人完全了解XX大学关于收集、保存、使用毕业设计（论文）的规定，即：按照学校要求提交毕业设计（论文）的印刷本和电子版本；学校有权保

2、存毕业设计（论文）的印刷本和电子版，并提供目录检索与阅览服务；学校可以采用影印、缩印、数字化或其它复制手段保存论文；在不以赢利为目的前提下，学校可以公布论文的部分或全部内容。作者签名： 日 期： 1 引言1.1 圆环耐撞性研究的意义近年来，各种车辆、船舶、飞行器的数量越来越多，速度越来越快，碰撞事故也随之日益增加，每年都要造成严重的生命和财产损失。解决这个问题的主要途径，是研究和提高各种车辆、船舶、飞行器结构的耐撞性(structural crashworthiness) 1。术语“耐撞性”指的是当车辆等卷入或经历碰撞时响应性质的优劣。众多学者已对金属吸能结构，如圆管的卷曲、方管的卷曲和渐进屈

3、曲、截壳体、圆管和圆环以及在其中填充泡沫塑料等进行了大量的研究。工程上特别关心如何评估结构受冲击后的整体性能、如何提出合理的设计方案及如何制定有效的防护措施等问题。目前，工程实际中比较切实有效的手段就是设计并采用性能优异的能量吸收装置(energy absorbing devices)，或者叫能量耗散装置(energy dissipating devices) 2，这项工作现己成为结构安全防护中重点课题。同时，在某些特定的情况(例如飞行器的紧急着陆，核电站和高速公路旁重要设施的防护等等)下，已有结构难以满足吸收全部碰撞能的要求，需要设计一些特殊的结构元件作为能量吸收装置，这种装置必须具有良好的

4、塑性变形性质，使得撞击动能通过结构的大塑性变形和破坏过程被耗散掉，从而大大降低冲击力的幅值，减缓冲击效应，使之保持在结构所能耐受的指标水平上。受横向载荷的圆环（图1.1）就是以塑性弯曲为其主要能量耗散机制的。 图1.1在一对刚性平板间受压的圆环 （a）初始构形；（b）变形后的形状在工作载荷下，传统结构（例如，用于图谋工程和机器中的结构）只发生很小的弹性变形。这些结构通常要求在规定的载荷下具有一定的强度和刚度，所以材料的选择和结构的设计主要基于结构必须承受的弹性应力或者应变。失效多是由于疲劳、腐蚀，或者是由于长时间使用引起的材料老化。能量吸收结构设计和分析与传统结构设计和分析大不相同。能量吸收结

5、构必须承受强碰撞载荷，所以它们的变形和失效涉及到几何大变形应变强化效应应变率效应以及不同变形模式（如弯曲和拉伸）之间的各种交互作用。因为这些原因，大多数能量吸收器用韧性金属制成，低碳钢和铝合金是使用最广泛的。目前在金属吸能装置领域的研究可分为两类：一类是金属吸能元件（包括与其它吸能材料组合而成的吸能装置）；另一类是利用金属塑性成形和金属切削原理的能量吸收装置。对于能量吸收装置的具体要求必须根据使用场合的不同以及具体使用工况的不同而定，不存在处处适用或绝对最优的能量吸收装置。由于碰撞工况的不同和安装部位的不同，对能量吸收装置的具休要求各不相同，但一般说来，以下几点要求是共同的：(1)碰撞动能应尽

6、可能不可逆地转换成变形能，也就是说，应该以塑性变形而不是弹性变形来贮存这种能量；(2) 在碰撞条件下，能量吸收装置的变形模式应当稳定，具有可重复性和可靠性, 且传递的反力波动小，保持常幅值；(3) 在吸收能量的过程中，应控制碰撞力和减速度，在大变形下应具有接近定常数的承载能力，以保护人员和重要结构。生物力学的研究指出人的颅脑系统的忍受度可用下式表达：其中GSI称为Godd指标，是加速度(或减速度)，以重力加速度g为单位，t是时间，单位是毫秒，T是经历加速度的总时间，1000是正常成人颅脑损伤的门槛值。从上式易证明，在给定初速和碰撞距离的条件下，使GSI取极小的减速过程，是为常数的匀减速过程。因

7、此，良好的能量吸收装置，在大变形下应具有接近定常数的承载能力；(4) 装在飞行器和汽车上的能量吸收装置，应该自重较轻，具有良好的“比耗能”，即单位自重所能吸收的能量值要高；(5) 为了吸收更多的总动能，它应能提供足够长的变形行程，而且在变形后不占据过大的空间或造成次生破坏；(6) 由于能量吸收装置通常都是一次使用结构，所以应该成本低廉，便于制造和更换3 18。1.2 圆环耐撞性的研究现状由于金属圆环结构具有变形稳定、承载力稳定、变形行程长、初始冲击力小、取材方便、易于更换以及产生的比能耗高等优势，受到了工程界的高度重视，取得了良好的经济效益和社会效益4。例如，在美国很多大型结构和高速公路的防护

8、装置都安装此类抗冲击吸能装置，取得了令人满意的结果。迄今为止，结构塑性动力学研究领域对梁、板、壳、拱等结构5以及多种管状、环状和其它方式的金属吸能结构在冲击载荷作用下的动力响应特性研究已取得了很大的进展。但是至今仍有一些极具研究价值和应用潜力的碰撞缓冲吸能金属结构尚未得到足够重视。这类吸能结构的动力塑性压溃变形模式与相应的准静态压溃变形模式有很大的差别。其原因使由于变形过程中往往掺杂了应变率效应和惯性效应的影响。圆环结构就是这类吸能结构之一。尽管圆管和圆环列作为新型吸能器已经得到一定的工程应用,但是以往的研究大多关注圆环或圆环列系统的变形行为，而对工程实际中更为关心的吸能特性(尤其是吸能效率)

9、 和机理方面的研究却很不充分，特别是有针对性的对圆环列系统的吸能效率方面的实验研究方面的工作几乎没有。Hwang最早给出了对径受压的理想弹塑薄圆环的小变形分析。自70年代以来，许多学者对完整或含缺口圆环在不同支撑下的塑性动力响应特性进行了实验研究。通过这些研究，获得了对圆环结构吸能特性的一些认识。Calladine和England对低碳钢圆环等进行了冲击实验。他们的实验表明，圆环使一种对冲击速度不敏感的能量吸收元件，也就是说，在相同的冲击动能条件下，改变落锤的质量和速度，圆环的最终变形相差不大。DeRuntz和Hodge6指出，理想刚塑性薄圆环在两刚性平板对压下的大变形为一个四铰机构。Reid

10、等7人详尽地研究了强化效应对圆环承载能力的影响。他们指出，材料线性强化的结果不但使塑性极限弯矩提高，而且使塑性铰扩展为塑性区1。Reid等人对强化受压圆环的处理同余同希1964年对强化受拉圆环的处理13是完全一致的。Reid等最先对冲击载荷下金属圆环列(若干个圆环排成一列)的动态力学行为进行了实验研究，得到了大量的不同连接条件下的圆环列系统的实验数据，并在此基础上提出了结构波理论。Reid和Reddy首先报道了圆环串接系统端部受集中质量块撞击的实验结果7。发现圆环的变形过程乃至完全压扁是一个接一个在进行 ,在宏观上表现出间断冲击波的特性, 为此Reid，Bell和Barr建立了一维冲击波理论模

11、型,这个模型能够反映出圆环串接系统的的瞬态塑性变形机制,并与实验进行了比较15。Reid等人也研究过圆环在V形槽内受压，以及在有侧向约束的情况下受压的问题。对于此种情形，圆环的大变形机构会在两个受载点中间形成塑性铰，相应的分析也很复杂。Choi等14讨论了材料的应变率效应对圆环列系统的影响，指出:忽略材料应变率效应可能导致对结构应力水平的估计值过低。有关圆环及圆管受准静态横向载荷作用的塑性变形分析方面的文献可参考6-9,而动载作用下圆环及圆管的响应特性研究可参考10-12,其中包括了理论研究和实验研究。近年来，国内关于这方而的研究也取得了大量的研究成果，余同希1研究了利用金属塑性变形原理的碰撞

12、能量吸收装置，从比耗能、相对行程和承载力的稳定性等指标对各类能量吸收元件进行比较。余同希3 15还对圆环、圆管、方管等结构的耐撞性进行了深入研究。宋宏伟17对圆管等能量吸收结构撞击吸能特性做了研究。姜锡权对多排圆环系统进行了大量的实验研究，其结果表明：在准静态条件下，多排圆环系统的力学行为基本上与单排圆环列系统的变形行为相类似；多排圆环系统的变形过程较单排圆环系统的变化更为平缓。曾首义等的实验结果也表明：圆环层数对变形影响很大，多层圆环系统的总变形较大，但其反力波形平缓、没有继生峰值。而顾红军等进行了相关的实验研究，对以往的模态解公式进行了修正，使得多排圆环系统的理论分析结果更接近实际情况。同

13、时，赵凯等4对圆环列系统在弹性范围内的应力波效应讲行了较为深入的研究。赵凯等4还通过对圆环列系统实验和数值研究，得出研究结果表明: (1) 圆环列系统是高效的能量吸收装置；(2) 考虑安装实用的要求，23个薄壁圆环系统是比较适合工程实际应用的方案。圆环列系统能够满足结构防护的工程要求，具有广阔应用前景和实际工程意义，值得进一步推广和应用。杨嘉陵和蔡序杰的研究侧重点在圆环系统的压扁失效的传播分析方面,包括有限元数值模拟和理论模型。理论模型将金属圆环等效为一弹塑性“弹簧”,并以此作为基本单元,研究了纵向串接的一系列圆环末端受撞击的弹塑性动力响应,进而预测各圆环不同时刻的变形情况。他们在对弹塑性圆环

14、串接系统受强动载荷作用的失效分析中，得出结论：(1) 对于串接的弹塑性圆环系统,提出了宏观等效的弹塑性质量- 弹簧系统,相应的材料与结构参数是基于单个圆环受强动载荷作用弹塑性动力分析的结果。(2) 通过定义塑性屈服点和圆环的失效准则,可以近似描述塑性波的传播过程,观察到圆环依次失效的发生、发展的传播过程。(3) 通过与有限元结果和实验结果的定量比较,说明这一理论模型能够很好的反映出圆环系统的弹塑性动力特性,对于进一步研究这种系统的能量吸收特性是有意义的。1.3 本文要研究的内容本文根据课题研究需要，在前人工作的基础上针对圆环结构的动力响应问题，从理论分析和数值模拟两个方面对圆环在压缩和冲击载荷

15、作用下的动力响应问题进行了研究，对金属圆环结构进行力学分析和研究，研究径向压缩和冲击中金属圆环结构的塑性变形行为以及能量吸收特性，推出金属圆环结构的载荷位移关系曲线。研究圆环在冲击过程中的变形形态及其对能量的吸收。在经典刚塑性理论的基础上，采用动力有限元程序LS-DYNA对金属圆环结构在冲击下的动态响应进行数值分析，同时对金属圆环结构进行压缩和冲击实验。统计整理实验数据、绘制金属圆环结构的载荷位移曲线，在此基础上初步探讨金属圆环结构的耐撞性问题。2 金属圆环结构的理论分析塑性弯曲作为一种能量耗散机制，在碰撞能量吸收装置中有着广泛的应用。其中，圆环就是一种研究得最充分应用得最广泛的能量吸收元件。

16、2.1 一对集中力作用下的受压圆环考虑一对方向相反的集中力作用下的刚塑性圆环 。它需要四个塑性铰以形成一个破损机构。利用能量法，我们可以得到其载荷位移曲线。另一方面这个问题也可以通过考虑一圆弧段平衡得到解答。当圆弧段的转角为时，压缩量为： （2.1）AB现在的长度为 （2.2）作用在这段圆弧上的力和力矩等效于两个数值相等（=P/2）作用方向相反的力。平衡条件要求这两个力必须沿同一条线作用，该线称为推力作用线。因此，力的平移必须等于AB/2，这里为圆环的塑性极限弯矩。将此引入方程（2.2）得出 （a）形成破损机构需要四个塑性铰 （b）作用在四分之一圆环上的力图2.1一对集中力作用下的受压圆环的破

17、损机构或者 （2.3）当时，得到初始破损载荷 （2.4）其中，D是圆环直径。联合式（2.1）、式（2.3）和式（2.4）给出载荷位移曲线（2.5）这个关系式表明，载荷随着挠度增加而减少。上述方法称为等效结构技术（equivalent structure technique ）（Mechant,1965;Reddy等,1987）。2.2 一对集中力作用下的受拉圆环当圆环受两个大小相等方向相反且向外作用的集中力作用，它的变形可以用和上面完全相同的方法分析。一个关键不同之处是，对应于最大弯矩的两个边上的塑性铰位置随挠度增加而移动图2.2（a）的B铰。（a）破损过程中的中间两个塑性铰（b）存在轴力时塑

18、性铰（c）作用在弧段上力的详图 分离成四个，如在B点 处的应力状态图2.2一对外向受拉圆环的破损机构未变形的弧段与变形了的直线部分总是在当前塑性铰B处相切。此外，圆环的总长度在变形过程中保持不变。根据这些考虑并利用上述等效结构技术，得到 （2.6）或者 （2.7）和 （2.8）上述分析中没有考虑到轴力对屈服的影响。当圆环比较厚以及接近1时，这种影响可能比较重要。将这种轴力效应近似考虑进来，每个塑性铰处的轴力是,抗弯能力由减少至于是，式（2.6）变成 （2.9）是轴向屈服力。联合式（2.9）和式（2.8），得到载荷位移曲线，当时最大载荷达到2。当截面关于点转过角度是时，中点的纤维被拉长了。这里是

19、系数，是与之间的距离。容易得到，轴力为，。上面塑性铰的轴力为零，下面塑性铰的轴力为。作用在拉伸圆环总的外力图2.2（a）为 （2.10）等效结构中力作用线的位置图2.2（c）为 （2.11）和 （2.12）再次根据平衡的考虑，两个为零方向相反的力必须沿着同一作用线。根据刚发生破损时的几何形状，位移为零，则 （2.13）将式（2.11）和式（2.12）代入到式（2.13）中，我们得到正的值的表达式为 （2.14）考虑到两种情况，和，方程（2.13）分别给出和。初始破损载荷变为 （2.15）这个公式与DeRuntz和Hodge所给出的圆环在刚性平板压缩下的初始破损载荷是一样的，他们考虑了轴力对屈服

20、的影响。2.3 两平板对压下的圆环在两刚性干板时压下的圆环弹塑性圆环对径受压的极限载荷为 （2.16）其中为圆环的半径，是截面的塑性极限弯矩.对于理想刚塑性薄圆环在两刚性平板对压下的大变形，DeRuntz和Hodge6指出，其大变形机构是一个四铰机构，如图2.4(a) 、(b)所示。由此得出大变形下的载荷一位移关系为 （2.17）其中是两板相对位移之半.此式表明，在大变形过程中，环的承载能力将逐渐增加。实验测得的圆环大变形承载能力较上式为高。于是，从1978年起，Reid等人对强化效应对圆环承载能力的影响作了详尽的研究。他们指出，材料线性强化的结果不但使塑性极限弯矩提高，而且使图中的塑性铰扩展

21、为塑性区。塑性区的大小取决于无量纲参数，其中是材料的初始屈服应力，是材料的线性强化模量，是圆环壁厚。这一参数越大，圆环的承载能力增长越迅速。两平板对压下的圆环的破坏需要有四个塑性铰，两个常见模式如图所示（Burton和Craig,1963;de Runtz和Hodge，1963）第一个模式有四个不动的塑性铰，较适合于软刚，因为它有上屈服点和下屈服点。第二个模式中圆环在移动接触点处被展平。这两个模式的未变形段的力作用图是相同的，因此得到相同的力-挠度曲线。初始破坏载荷与集中力作用下的情形相同式2.4。由平衡条件有 （2.18）图2.3刚性平板对压下的圆环由几何关系有 （2.19）联立式（2.18

22、）和式（2.19），并注意到式式（2.4），有 （2.20）或者 （2.21）其中，L是圆环的宽度或者圆管的长度。这说明了载荷随挠度的增加而增加（图2.3）。值得注意的是，至今为止对圆环所提出的分析，也同样适用于在类似载荷作用的圆管，只要屈服应力选取适当的值。于是，长度不超过其壁厚几倍的短圆管可以看做是圆环，式（2.18）中Y等于简单拉伸试验得到的屈服应力。3 圆环的压缩和冲击试验研究3.1 圆环的压缩试验3.1.1 试验仪器名称：电子多功能试验机 型号：CMT 5105A 规格：100 KN 精度等级：0.5工作原理：试验机利用两个压头对试件以恒定的速度进行压缩，这个速度可以通过电脑设置；通

23、过力传感器和位移传感器，可以在电脑上显示其数据和图象，这样既可以知道压头的位移，同时可以看到实时的载荷，防止载荷太大产生过载。3.1.2 试验试件表3.1 圆环试件的尺寸(mm)圆环内径圆环外径圆环壁厚圆环高度5458210试件材料类型：低碳钢，型号：Q2353.1.3 试验内容试验前，首先把要进行试验的圆环准备好，圆环试件的尺寸如表3.1所示。接通电子多功能试验机的电源，打开试验机开关，检查试验机运行状况，确定运行正常。放置试件到压缩基座上，选用压头，使用电脑控制上压头，并使其下降，一直到与试件上段轻微接触的高度，停止上冲头的下降，把位移和轴力清零，设定下降速度、位移等；点击按钮“Run”，

24、开始进行静压试验，观察试件的变形情况，并用相机记录变形图片；进行静压直到冲头到合适位置，点击“Stop”按钮，停止下压，提升冲头，保存试验图象和数据；记录其最终变形情况。图3.1 中北大学万能试验机照片3.1.4 试验结果实验1： 图 3.2 试件1初始图图 3.3 试件1变形图一图 3.4 试件1变形图二 图 3.5 试件1变形图三图 3.6 试件1变形图四图 3.7 试件1变形图五图 3.8 试验测得的压缩圆环1的载荷位移曲线图 3.9 压缩圆环1的载荷位移曲线表3.2 试验得到的压缩圆环1的载荷位移数值载荷N位移mm载荷N位移mm载荷N位移mm1.10.01479.359.49636.5

25、27.0827.70.15481.559.99640.9527.4928.550.2485.9510.48649.827.9937.650.26490.411655.3528.5147.60.29495.9511.52671.9529.4964.20.35498.1512687.4530.4980.80.4502.612.5706.2531.595.20.45504.812.99730.632.52112.90.51509.213.49737.2533.02126.20.56511.4513.99748.3533.51137.250.59518.0514.49759.434152.750.65

26、520.315771.634.53167.150.7523.615.5783.7535181.550.75528.0515.99795.9535.49199.250.81532.4516.48808.135.99212.550.8653816.99822.536.51219.20.89541.317.52835.836.99234.70.94545.7517.99849.0537.5250.21549.0518.51865.6538352.051.5551.318.96878.9538.5396.32556.8619.46897.7539419.552.49560.1519.99917.739

27、.51430.63563.4520.52938.7540.01438.353.5256920.99955.3540.514454.01570.121.5198341.03448.354.49577.8522.021007.3541.5452.754.99577.8522.491031.742456.15.48581.1522.981058.342.49457.26584.523.521062.742.52459.46.5588.924.011087.0542.99463.857592.2524.51089.343.02467.157.49594.45251118.0543.48469.3585

28、98.925.491122.543.53473.88.5609.9525.961162.3544477.18.97624.3526.51205.544.51243.1545实验2：图 3.10 试件2的初始图 图 3.11 试件2变形图一图 3.12 试件2变形图二图 3.13 试件2变形图三 图 3.14 试件2变形图四图 3.15 试件2变形图五 图 3.16 压缩圆环2的载荷位移曲线 图 3.17 试验测得压缩圆环2的载荷位移曲线表 3.3 试验得到的压缩圆环2的载荷位移数值载荷N位移mm载荷N位移mm载荷N位移mm1.10612.155.01715.119.974.40.07616.5

29、56725.0520.999.950.13619.96.5729.521.4923.20.16623.27.5740.5522.5136.50.2632.058.52750.523.52520.24635.49.02753.8524.0173.050.3636.59.5760.524.4994.050.35642.059.99769.3525.0399.60.37640.9510.01772.6525.51115.10.41646.4510.5780.426131.70.45646.4511.04788.1526.5147.20.4965211.53795.927.01166.050.5465

30、4.212802.5527.5190.40.61657.5512.48810.328.01204.750.65658.6512.99818.0528.5224.70.7664.1513.4982828.99242.40.75666.413.99834.6529.48255.70.79667.514.05844.630275.60.85666.414.07853.4530.49295.550.9665.314.13862.3530.99313.250.96667.514.19872.331.49325.451668.614.21883.3532353.11.09667.514.24895.553

31、2.5375.251.16667.514.3908.833399.61.25669.714.3692133.48421.751.35670.814.45932.0533.99456.051.5670.814.59944.2534.49471.551.59671.914.79959.7535487.051.7674.1514.99977.4535.55071.85674.1515.041005.1536.5525.82.01674.1515.491040.5537.5561.22.5677.4515.51058.2538.02566.752.65680.816.021074.8538.5574.

32、52.8686.316.51094.838.97578.953687.416.981114.739.5588.93.24695.1517.511140.240597.753.7698.518.551151.2540.17609.954.51708.4519.47上面这些图片是圆环试件1和2在压力下的变形情况。在刚进入变形的阶段, 载荷和位移呈线性关系；当载荷继续增大，随位移的增加, 载荷非常缓慢地增大, 载荷位移曲线为一斜率很小的直线；载荷随位移的增加而快速增大。图3.2是圆环1的初始形状，图3.3中圆环顶部首先开始发生凹陷，此时圆环内已经形成两个塑性铰在圆环的顶部和底部，但是这尚不足使圆环成

33、为机构。随着压缩的继续，载荷继续增大，圆环的顶部和底部都向圆心方向产生凹陷，圆环的顶部和底部还有左右端都形成了四个塑性铰，圆环此时成为机构，如图3.4中所示；在整个圆环的压缩过程中，圆环处于一个对称的塑性变形阶段（如图3.5和3.6所示），圆环通过塑性变形吸收能量。从数据文件中可看到，在位移的增加过程中，有时在一定的范围内载荷会出现较小的减小，然后，再增大。当变形进行时，塑性弯矩抗力增加。而且塑性变形将发生在一个区域内，而不是限于一个局部化的塑性铰，即铰有一定长度，这个几何变化改变了力臂长度，导致载荷的改变。由于载荷作用点随位移的增加而外移，载荷对边铰的力臂减小，因而圆环的承载能力将随变形而增

34、大。下面是Reid于1983年给出的受压圆环大变形的载荷位移曲线：图3.18受压圆环大变形的载荷位移曲线的实验和理论比较（h/R=0.108,R=42.16mm,L=101.6mm）（Reid，1983）预测的力低于实验结果，缘于塑性弯矩抗力增加和力臂长度改变导致的载荷改变。在以上两个圆环的压缩试验中，得到的载荷位移曲线和经过数据处理后绘制的载荷位移曲线是一致的。并且与Reid给出的曲线及其相似。3.2 圆环的冲击试验冲击试验可以通过落锤，摆锤或倾斜的滑轨进行。当使用落锤时，有关装置将一个质量块（落锤）提升至一定高度后释放，从而使落锤对放置在试验台基座上结构的冲击。所能获得的下落最大速度是由落

35、锤的总高度决定的（通常利用滚柱轴承，尽量减少垂直导轨上的摩擦）。所使用的典型的仪器包括附着在冲击质量块上的加速度计，测量冲击前速度的仪器（通常是测量通过一段已知距离的时间间隔）记录冲击质量块运动的位移传感器以及放置于被测量结构下面的动态力传感器。图3.19 中北大学液晶全自动落锤冲击万能试验机照片3.2.1 试验仪器名称：液晶全自动落锤冲击试验机最大冲击能量：300 J最大冲击高度：2000 mm锤头重量范围：0.2516kg高度误差：10 mm3.2.2 试验原理1、在只考虑重力作用下的近似理想状态下，落锤式试验机的势能改变量等于其动能的改变量： 2、重物下落时的速度和位移间的关系：3.2.

36、3 试验试件 表3.4 圆环试件的尺寸(mm)圆环内径圆环外径圆环壁厚圆环高度5458210试件材料类型：低碳钢，型号：Q2353.2.4 试验内容试验之前，首先要准备好要进行冲击试验的圆环试件，并用相机记录它们的初始形状。然后，接通试验机的电源，打开试验机开关，检查试验机的运行状况是否正常，如正常，按照设计好的标准选取落锤，并安装落锤；放置试件到试样台上，把冲头放置到试件之上，调整使之在中心；关闭试验机的小门，点控制板上的“set”按钮，利用箭头设定下落高度为1.276m，点“set” 按钮确定；设定好高度后，点“功能选择”键，切换到“自动”状态，这时下面的一排按钮的功能就是键上方汉字标明的

37、，这时，按“定位”键，试验机就会自动定位，定位好以后，按“提锤”键，落锤的高度就会被提升到1.276m的高度，这是由试验机的位移传感器自动运行的，看液晶显示板上是否显示“可落锤”，如果显示，且无任何异常现象，则点击“落锤”按钮，锤就会自动落下，达到冲击试件的目的；待到试验机自动提锤结束后，点“功能选择”键，切换到“手动”状态，这时下面一排按钮的功能就是键下的字所标明的，然后，点“试样台降”键，使试样台下降一段距离，打开小门，拿出试件，观察其变形情况，并记录。表3.5 试件的冲击参量试件编号落锤质量(kg)落锤高度(m)冲击速度(m/s)冲击能量(J)111.276512.5221.276525

38、361.2765753.2.5 试验结果图 3.20 试件1变形图 图 3.21 试件2变形图 图 3.22 试件3变形图4 圆环冲击的有限元分析随着计算机技术的高速发展和以有限元法为基础的工程计算方法的发展已日趋成熟，数值方法已成为碰撞研究的最重要的手段之一。用于碰撞仿真研究的最常用的大型有限元软件是ANSYS/LS-DYNA。本章将重点讨论LS-DYNA软件的基本算法以及仿真研究中与之相关的主要问题的处理方法，如材料性质、单元选择、接触碰撞界面算法等，这些问题的理解和正确处理对于提高仿真结果的可靠性是非常重要的。4.1 LS-DYNA发展历程LS-DYNA程序是世界上最著名的通用显式动力分

39、析程序，最初称为DYNA程序，由J.O.Hallquist博士于1976年在美国Lawrence Livermore National Laboratory主持开发完成。能够模拟真实世界的各种复杂的问题，特别适合求解各种二维、三维非线性结构的高速碰撞、爆炸和金属成型等非线性动力冲击问题，同时可以求解传热、流体及流固耦合问题。其特点主要包括以下几个方面：1 强大的分析能力2 丰富的材料模型库3 易用的单元库4 充足的接触方式5 自适应网格剖分功能6 ALE和Euler列式7 SPH算法8 边界元法9 隐式求解10 热分析11 不可压缩流场分析12 跌落测试分析13 强大的软硬平台支持4.2 LS

40、-DYNA分析的一般流程与一般的CAE辅助分析软件操作过程类似，ANSYS/LS-DYNA分析过程也包括问题的规划、前处理以及后处理四个部分，如图4.1所示。图4.1 LS-DYNA分析流程在“分析问题的规划”中要综合考虑问题的特点、计算精度、时间；在“前处理”中要指定单元类型、实常数、材料模型；建立有限元模型；定义接触、边界条件，施加载荷；输出K文件；在“加载与求解”中要设置求解参数并求解。4.3 求解过程4.3.1 建立模型 依次建立圆环模型和上下实体模型。圆环尺寸为外直径58mm内直径54mm高度10mm。上下实体尺寸为长100mm宽70mm高70mm。本文采用国际单位制-mm-s-mN

41、-kPa，圆环的密度为7.8E-6，柏松比为0.27，杨氏模量为2.07E8，屈服极限为3E5，切变模量为1E7。上下实体的密度为7.8E-6，柏松比为0.3，杨氏模量为2.07E8。图4.2 实体模型4.3.2 所用单元简介在本文中对圆环模型进行划分时所用单元类型为SHELL163单元，SHELL163单元是用于3维显式结构实体单元，由4节点构成，图 4.3 描述了SHELL163几何特性、节点位置和坐标系。这个单元只用在动力显式分析，它支持所有的非线性特性。图4.3 SHELL163实体单元几何特性本文中对圆环进行冲击的实体模型进行划分时所用单元类型为3D-SOLID164单元，图 4.4

42、 描述了SOLID164几何特性、节点位置和坐标系。图4.4 SOLID164实体单元几何特性依次选择材料模型,通过路径 LS-DYNANolineearInelasticIsotropic HardeningBilinear Isotropic，分别输入DENS密度、EX 杨氏模量、泊松比NUXY 、Yield Stress屈服极限、Tangent Moduls切变模量的数值，通过路径LS-DYNA LineearElasticIsotropic，分别输入上下实体的DENS密度、EX杨氏模量、NUXY泊松比。利用MeshTool工具分别对圆环、上下实体进行网格划分,得到有限元模型,如图3.6

43、所示。圆环被划分成1408个单元，590个节点。上下实体各被划分成138个单元，60个节点。图4.5 有限元模型选择Main MenuPreprocessorLS-DYNA OptionsParts Options命令,生成PART。选择Main MenuPreprocessorLS-DYNA OptionsInitial VelocityOn Part，输入速度-5000mm/s，将计算终止时间设为0.008s，将结果输出文件.RST的输出步数设为100步,时间历程文件.HIS的输出步数设为10。选择Main MenuSolutionWrite Jobname.k命令,生成k文件。打开k文件

44、,修改文件中的关键字和参数并保存。启动LS-DYNA求解器,选择修改后的k文件,调用LS-DYNA970求解器求解。4.3.3 后处理本文采用ANSYS/LS-DYNA内含的LS-PREPOST处理器进行圆环冲击的数值模拟的后处理分析。运行LS-PREPOST,选择FileOpenBinary Plot命令,读入d3plot结果数据文件。以下是圆环在冲击下的变形过程图：图4.6 圆环冲击的初始图片图4.7 圆环冲击变形的临界状态图片图4.8 圆环冲击变形图图4.9 圆环冲击变形图图4.10 圆环冲击变形图图4.11 圆环冲击变形图图4.12 圆环冲击变形最后的图片从以上的圆环冲击变形图中，可以

45、看到圆环的变形形态与压缩试验得到的圆环的变形形态是一致的。从而可以通过以下节点和单元在模拟冲击过程中的应力、位移、能量的各种曲线和云图，对圆环的冲击动态响应进行研究。图4.13 节点9和节点11在圆环中的位置节点9和节点11在圆环在冲击过程中的位移曲线图如下：图4.14 节点9和节点11在冲击过程中的X轴位移时间曲线图4.15 节点9和节点11在冲击过程中的Y轴位移时间曲线图4.16 节点9和节点11在冲击过程中的Z轴位移时间曲线图4.17 节点9和节点11在冲击过程中的总位移时间曲线单元613的应力时间曲线和应变时间曲线如下：图4.18 单元613在圆环上的位置图4.19 单元613在冲击过

46、程中的X轴应力时间曲线图4.20 单元613在冲击过程中的Y轴应力时间曲线图4.21 单元613在冲击过程中的Z轴应力时间曲线图4.22 单元613在冲击过程中的Von mises 有效应力时间曲线图4.23 单元613在冲击过程中的有效塑性应变时间曲线节点10在冲击过程中的位移时间曲线：图4.24 节点10在圆环上的位置图4.25 节点10在冲击过程中的X轴位移时间曲线图4.26 节点10在冲击过程中的Y轴位移时间曲线图4.27 节点10在冲击过程中的Z轴位移时间曲线图4.28 节点10在冲击过程中的总位移时间曲线圆环顶部单元1378 437 946 1370在冲击过程中的应力时间曲线和应变

47、时间曲线如下：图4.29 圆环顶部单元1378 437 946 1370的位置图4.30 圆环顶部单元1378 437 946 1370在冲击过程中的Y轴应力时间曲线图4.31 圆环顶部单元1378 437 946 1370在冲击过程中的Z轴应力时间曲线图4.32 圆环顶部单元1378 437 946 1370在冲击过程中的Von mises 有效应力时间曲线图4.33 圆环顶部单元1378 437 946 1370在冲击过程中的有效塑性应变时间曲线圆环在0.0010399s、0.0021599s、0.00616s、0.00725s的Von mises应力云图：图4.34 圆环冲击过程中不同时

48、间的Von mises stress云图圆环在0.0010399s、0.0021599s、0.00616s、0.00725s的Y轴应力云图图4.35 圆环冲击过程中不同时间的y轴方向应力云图圆环在冲击过程中的内能时间曲线如下：图4.36 冲击过程中圆环的内能时间曲线从圆环内能时间曲线可以看到圆环吸收了能量使内能增加。圆环在冲击过程中的动能时间曲线如下：图4.37 冲击过程中圆环的动能-时间曲线图4.38 冲击过程中圆环的总能量-时间曲线5 总结随着人类社会的进步及科学技术的不断发展，耐撞性的研究工作越来越引起人们的重视。结构耐撞性的研究，主要关心的是低、中速碰撞，但碰撞物所携带的动能仍然可以是

49、巨大的。因此，在碰撞中结构所吸收的总能量是这个问题中最主要的一个参数。在某些特定的情况，现有的结构无法满足吸收全部碰撞动能的要求，这是就需要采用一些经过特殊设计的结构元件作为能量吸收装置。设计良好的能量吸收装置可以使能量的耗散以一种有控制的方式进行，即不但能吸收设计要求的碰撞总能量，而且能按设计要求控制碰撞力和减速度的变化。除此以外，结构的耐撞性研究还涉及到分析结构在撞击下特性的细节，包括结构是如何控制碰撞减速度的，以及结构变形和破坏的具体模式。在本篇毕业论文中，我们学习、介绍了已有的研究成果，对已有的理论进行了深入的学习、理解。前面对耐撞性和能量吸收的概念进行了介绍，同时，介绍了耐撞性理论在

50、现实的生产和生活中的广泛应用，和其不可或缺的地位。由于现实中撞击的普遍性和缓冲吸能的重要性，就有必要对缓冲吸能材料进行深入的研究。材料的吸能有多种形式，有的在吸能过程中发生极小的变形，有的发生极大的变形，但是要得最好的吸能材料，就需要对其吸能性进行试验分析，以设计出更好的缓冲吸能元件。塑性弯曲作为一种能量耗散机制，在碰撞能量吸收装置中有着广泛的应用。由于金属圆环结构具有变形稳定、承载力稳定、变形行程长、初始冲击力小、取材方便、易于更换以及产生的比能耗高等优势，受到了工程界的高度重视，取得了良好的经济效益和社会效益。因此，圆环成为一种研究得最充分应用得最广泛的能量吸收元件。第二章中我们简单介绍了

51、一对集中力作用下的受压、受拉圆环以及在两平板对压下圆环的理论分析。第三章中对圆环进行了压缩和冲击试验，通过试验研究了圆环在压缩和冲击过程中的变形模态，从试验中得到的圆环在压缩过程中的载荷位移曲线。同时对试验得到的数据分析，绘制了圆环在压缩过程中的载荷位移曲线。第四章中，在经典刚塑性理论的基础上，采用动力有限元程序LS-DYNA对金属圆环结构在冲击下的动态响应进行数值分析，从计算机模拟得到的圆环冲击变形图中，可以看到低速冲击下的圆环的变形形态与压缩试验得到的圆环的变形形态是相似的。通过分析圆环在冲击过程中关键节点和单元的应力、位移、能量的各种曲线和云图，对圆环的冲击动态响应进行研究。研究结果表明

52、，圆环结构是高效的能量吸收装置。具有广泛的应用前景和实际工程意义，值得进一步推广和研究。参考文献1 余同希.利用金属塑性变形原理的碰撞能量吸收装置.力学进展,1986,16(1)：28392 Johnson W, Reid S R. Metallic energy dissipating systems.Applied Mechanic Review,1978 ,31 :2772883 余同希.结构的耐撞性和能量吸收装置.力学与实践,1985,3:294 赵凯,沈建虎,刘凯欣等. 圆环列系统吸能特性研究.北京大学学报.2007，43(3):3123165 杨嘉陵,余同希,王仁. 结构塑性动力响

53、应当前的研究进展和重点.力学进展, 1993,23(1):23336 Deruntz J A,Hodge P G. Crashing of a tube between rigid plate J.Appl Mech ASME 1963,30:3913957 Reid S R ,Reddy T Y. Effects of strain hardening on the lateral compression of tube between rigid plate J.Int J Solids Structures,1978,14 :2132258 Redwood R G. Discussion

54、 of Ref. 9 J.Appl Mech ,1964,31:3573589 Reid S R ,Reddy T Y. Phenomena associated with the crushing of metaltubes between rigid plateJ. Int J Solids Structures,1979,16:54556210 Shim V P W,Quah S E. Effects of geometrical and loading parameters on wave propagation in a circular ringJ.International Jo

55、urnal of Impact Engineering,2001,25:10322211 Shu D W,Lim B B. Metal tube under lateral compression,engineering plasticity and impact dynamicsZ. Editor by L Zhang,2001,759512 张铁光,郭晓刚. 受横向冲击圆环的粘塑性大变形分析J.爆炸与冲击,1990,10(1):556113 余同希.对径受拉圆环的塑性大变形.力学学报,1979,1:899114 Choi S T,Carney J F,Veillette J R. Ener

56、gy-dissipation characteristics of steel tubes.Experimental Mechanics 1986,6:30130515 杨嘉陵,蔡序杰.弹塑性圆环串接系统受强动载荷作用的失效分析.航空学报，2002,23(5):42142616 余同希,章亮炽.塑性弯曲理论及其应用.北京:科学出版社,199217 宋宏伟,赵桂范,杜星文.提高汽车耐撞性的能量吸收结构撞击吸能特性研究设计.计算.研究,2000,1:111418 顾力强.轿车保险杠和金属缓冲吸能结构的耐撞性研究.博士学位论文. 上海:上海交通大学,200019 余同希，卢国兴.材料与结构的能量吸收.北京：化学工业出版社.200620 余同希. 结构塑性动力响应的研究进展.冲击动力学进展. 合肥：中国科学技术大学出版社,1992致 谢本文是在孙华东老师的悉心指导、启发和热诚的帮助下完成的。孙老师严谨的治学风格、谦