《落锤冲击性能的数值模拟》由会员分享，可在线阅读，更多相关《落锤冲击性能的数值模拟（18页珍藏版）》请在金锄头文库上搜索。
1、 落锤冲击下泡沫铝的应力数值模拟 作者：R. Rajendran *, A. Moorthi, S. BasuBARC 原子设备研究中心，卡尔帕尔姆泰米尔纳德邦，印度摘要：在自由落体的铁锤的冲击下，我们进行泡沫铝的数值模拟仿真实验。所进行的静态轴向压缩试验是基于三种不同密度的材料，在此之前我们要获取材料性能。实验结果表明，粗网格性能材料能够验证真确的预测结果。通过仿真实验，我们对不同密度的泡沫，不同的冲击速度下的泡沫铝参数进行了研究，最终得出了实验的结论。关键词：泡沫铝位移应力能量. 1.导论 实验的对象我们选择了泡沫铝材料，其主要原因在于木桶运动中的能量吸收原理。在自由落体下的木桶，动能的变
2、化短暂的使动态应力发生非常高的程度的变化，这些应力可以损害完整的木桶。国际原子能机构（原子能机构）安全标准系列规定，木桶在 9 米的高度进行自由落体，模拟最有害的下降速度达到了 48 公里每小时。自由落体冲击实验进行了模型缩比，从而仿真了木桶的模拟下降状态 4 6 。泡沫铝之间的模型和刚性表面大大降低了作用在模型上的力，使其所衍生的应力小于它的承受应力。泡沫铝的数值变型实验被不同的科学家所研究。马吉德【7】等人建立了准静态的三维非线性有限元模型，并且进行了破碎的填充泡沫铝箱实验，通过假设一个米塞斯型材料以及标准化了的硬化铝盒泡沫来建成了水平的模型，而且确保每一层都含有一个单层的固体元素。每一个
3、单层的节点都与其相邻的节点相连接。然后根据 LS DYNA，aktay 软件进行了准静态破碎挤塑聚苯乙烯泡沫薄壁铝管分析实验，实验中这些节点就会紧密的结合在一起，使其数值解具有有限元代码功能。铝管节点模型使用（belytschko-tsay-40）材料的薄壳单元，同时实验中的聚苯乙烯泡沫模型使用压碎泡沫固体模型。里佐夫 9 调查弹性塑料行为的闭孔泡沫细胞受到点和线负载的实验和数值模拟的影响。这次数值模拟进行了利用 ABA QUS 软件。泡沫的塑性模型使其被作为很好的压碎泡沫硬化材料模型。张、程 10进行了泡沫填充方列和多细胞方柱的能量吸收特性实验，并且对两者进行了比较，实验结果观察到了很好的仿
4、真模型和实验数据。masso-moreu 和米尔斯 12通过进行了数值评估压缩冲击响应，以及用于保护性的包装锥体聚苯乙烯泡沫实验，开发了一个金属的模型和节点的使用，他们的结果通过了物理的检测。本研究我们提出的有限元模拟闭孔泡沫铝模型，并且同时进行了轴向冲击气缸破碎的自由落锤数值模拟实验。本实验仿真验证与实验结果都收录在参考文献 1，泡沫铝冲击的参数研究是基于铁锤不同的冲击速度和泡沫铝密度的变化，通过两个参数的变化，更加准确的进行实验的研究。 2实验细节 对冲击试验详细的过程被记录在参考文献 1。一个直径 155mm 和320mm 高的钢锤通过加载电磁而后与辊和盘融合在一起。该锤重达 106 公
5、斤。图示为钢锤跌落实验显示图，图示说明：1 泡沫铝是保持轴向垂直的状态，摄影角度如实验显示图所示。2 跌落试验进行不同下降高度的冲击测量，并且利用不同密度的泡沫变形来测量轴向的数据。图 1 用于跌落试验铁锤。图 2 泡沫铝冲击试验的摄影拍摄点。3. 材料参数负载挠度试验是基于密度为 292，471 和 760 千克/立方米的泡沫铝上进行的，通过密度的变化，我们就能测量出其负载的准确真应力和应变。密度 292千克/立方米泡沫的作用是在应力和应变的条件下，来建立参考的变，这一公式为： -（1）公式的 m 有一定的范围。这种线性变化的趋势表明该地区的参数会随着区域的浓度而增加，另外一个参数会随着区域
6、的密度而增加。此属性被利用来产生材料模型。由于与密度相符合的只有一种，所以材料模型将会变得越来越不精确。这一计算公式表述为: -（2 ）-（3）-（4 ）用来模拟真实应力应变图的泡沫密度分别为 380，450，570 和 579 千克/ 立方米，这同样也能被数值模拟的图所示，如图 3-b。同时泡沫密度的300，400，500，600，700 和 800 千克/立方米也能被用来应力应变的模拟图，如图 3-c。图 3 应力应变图（一）实验测量和预测。（2）应力应变图及模型的有限元分析。（三）应力 应变图参数研究图 4。位移时间变化图4.有限元建模 有限元是用于数值分析的软件。当从泡沫铝从规定的高度
7、自由落体时，106公斤的铁锤产生的动能就通过其底座传递到 155mm 直径的泡沫铝上。因此，我们将锤辊和夹持盘作为一个长度为 720.2 毫米，直径相同的等效模型。在实验过程中，八节点六面体固体元素与双线性材料模型（表 1），也得到了应用。泡沫铝的尺寸我们设定为 100 毫米的长度和直径 80 毫米，通过赋予八节点六面体实体元素使其成为一个合理的实验模型，这样固支边界条件也同样适用于泡沫的底部，在这里，压碎泡沫材料模型也得到了应用。自动节点之间我们采用面接触算法，在泡沫顶部表面和锤之间，我们运用了自动节点以及面接触算法。初始速度计算从自由下落高度开始计算，并且一些科学家提出了弹性闭孔泡沫的密度
8、呈线性变化理论。然而，敏感性研究显示该弹性模量变化并没有改变泡沫的反应。在这里泊松比为 0，张力降低为 4 兆帕，粘性阻尼系数0.05。5.网格细化研究实验中网格的研究是基于边长为 2，4 ，5，10，15 和 20 毫米泡沫网。这个研究是关于泡沫顶端的节点瞬时错位效应和铁锤底端的受力反应的，根据双方的时间位移（如图 4 所示）和时间反应力（如图 5 所示）来建立坐标系，从图中可以看出他们几乎不受网格尺寸的影响。因此，我们将泡沫铝的边缘选择为 10 毫米长度。同时，我们对有限元模型的铁锤泡沫选定分析具体在图 6中有所显示。图 5。当泡沫铝为 570 公斤/立方米 ，铁锤自由落体冲击速度为 5.
9、42 米/秒，我们所得出的时间-位移的坐标图。图 6。所建立的泡沫铝冲击的有限元模型图 7.比较实验与数值预测。（1）泡沫密度=570千克 /立方米；下降高度1米冲击能量=1040 （焦耳）/秒（2）泡沫密度=579千克 /立方米；下降高度1.5 米；冲击能量=1557 （焦耳）/秒（3） 泡沫密度=380千克 /立方米；下降高度1.75米；冲击能量=1820（焦耳）/秒（4）泡沫密度=450千克 /立方米；下降高度2米；冲击能量=2076 （焦耳）/秒6。模型验证当模拟的泡沫密度分别为 570，579 ，380 和 450 公斤/ 立方米时，其冲击速度分别为 4.43， 5.42，5.86
10、和 6.26 米/秒。此时与 106 公斤的铁锤相比可发现轴向影响的不同，结果在参考文献有记录。如图 7，相应的冲击能量为1040，,1557，1820 和 2076（下降高度为 1，1.5，1.75 和 2 米）。显然实验和模拟的结果是一致的。在速度 5.42 米/秒和泡沫密度为 579 千克/立方米，位移的等值线的变化反应在图 8。泡沫顶部有一个一次性永久 35.6 毫米的位移，从顶部到底部，位置偏移成线性变化，图 9 所示为铁锤的冯米塞斯应力图。同时，锤子受到了微弱的压力，泡沫受到最大压力为 14.2 兆帕，持续时间为 12毫秒。图 8。方向位移等值线（当泡沫泡沫密度 579 公斤/ 立
11、方米和冲击速度为 5.42米/秒）。图 9。冯米塞斯应力（泡沫泡沫密度 579 公斤/立方米和冲击速度为 5.42 米/ 秒）7。参数研究随着泡沫铝密度，锤冲冲击击速度的变化，泡沫所产生位移也随之变化。如图 10 所示，当铁锤冲击速度从 1 到 5 米/秒，密度为 500kg/立方米的泡沫变形从 4 到 38.5 毫米的，而同样的应力加载在 800 千克/立方米泡沫变形为 120.3毫米。泡沫顶部的应力随着铁锤的冲击速度和泡沫的密度发生变化，如图 11显示。当铁锤冲击速度从 1 增加到 5 米/秒，对于泡沫的密度为 500kg/立方米，冯米塞斯应力从 3.22 增加到 11.36 兆帕，而同样
12、的应力对于泡沫密度为 800公斤/立方米怎为 6.8817.50 兆帕。铁锤底部的应力变化也受着泡沫密度和速度的影响，见图 12。冯米塞斯应力乘以泡沫顶部的面积等于最终的反应力。当铁锤冲击速度从 1 增加到 5 米/秒，导致反应力从 16.02 到 56.41，而相应的反应为对于 800 千克/立方米泡沫增加量为 34.0286.71 千牛 。当铁锤的冲击速度为 4 米/秒，泡沫与泡沫密度位移和时间变化图，如图 13 所示。当泡沫密度为 300，400，500，600，700 和 800 公斤/立方米，并且铁锤冲击速度为 4米/秒，最终的变形为 43.8，36.9 ，29.3，22.8，18.
13、2 和 14.6 毫米。相同的时间下，达到最终塑性变形分别为 23.4，20.6 ，17.8，16，12.9 和 11.3 毫秒。因此在不断的铁锤冲后，需要由泡沫的密度和时间来最终确定变形量。一个密度为 600kg/立方米的泡沫受到不同速度的冲击后表面力的反应也不同，变形时间也不同，如图 14 所示。当冲击速度 1，2，3 ，4，5 ，6 和 7 米/秒，最后的变形是 2.4，7.8，15.1，22.2 ，30.1， 37.9 和 45.1 毫米。而相应的时间达到最终塑性变形却需要 12.3，14，15.1，15.4，15.4，15.4 和 15.1 毫秒。 泡沫铝变形的能量和其变形时间和速度
14、的关系，在图 15 中表现出来。当锤的冲击速度 1，2 ，3，4，5，6 和 7 米/秒，动能分别为53，212，477，848 ，1325，1908 和 2597。对于 600 公斤/立方米密度泡沫相应的非弹性变形能量为 33.8，167 ，409，745，1174，1695 和 2313，对于700 千克/立方米密度泡沫非弹性变形能量为31.2，169，415，754，1188，1715 和 2355，对于 800 公斤/立方米密度泡沫的弹性变形能量为 28.58，172，417，759，1204，1740 和 2364。当比较总能量（弹性变形+弹性）与势能下降的关系时，我们发现所有的速度
15、都有紧密的联系。图 16 是总的变形能中的一部分分散的弹性势能，当冲击速度 1 米/ 秒时，泡沫附近的弹性能量约 36%。随着速度的不断增加，势能将下降到 8%左右。此外，除了 1m/s 的下降速率外，泡沫密度的增加也会对弹性能量有影响。这种因为 1 米/秒的冲击速度而造成的趋势主要是归结于力，位移，时间的数值偏差。对于泡沫密度 600 公斤/立方米，当受到的冲击速度变化从 1 到 7 米/ 秒，所受的应力和位移的变化如图 17 所示。该反应力的加载与速度的变化率有较大的关系，而力的卸载去弹性形变有较大的关系。随着变形程度和冲击速度的变化卸载的范围从 12.31 到 27.09 千牛/毫米并且时间从零到它的最大值的范围为4.9 到 3.3 毫秒。弹性能量耗散随着卸荷曲线的面积增加而增加并且受到了速度的影响。然而，能量耗散的一小部分也总是随冲击速度的变化而变化的。图 10。变化的塑性变形与泡沫的密度和速度的关系 。图 13。当速度为 4m/s 时，泡沫的位移随密度和时间的变化。图 11。冯米塞斯应力随泡沫密度和冲击速度而变化。图 14。密度为 600kg/立方米的泡沫随着冲击速度和时间的变化 。图 12。锤底变化的反应力随着泡沫密度和速度的影响 图 15。变化，变形能随着锤冲击速度和时间的变化而变化 （速度）图 16。如图表示的泡沫铝所受弹性能量和冲击速度
