Design and crashworthiness analysis of thin-walled tubes based on a shrimp chela structure
-
摘要: 为提高薄壁管的耐撞性能，以虾螯为生物原型，通过结构仿生原理设计了仿虾螯结构多晶胞薄壁管。以晶胞数（2～6）和冲击角度（0°、10°、20°、30°）为试验因素，利用有限元法分析了仿虾螯结构多晶胞薄壁管在不同冲击角度下的耐撞性能，通过落锤试验验证了仿真结果的可靠性。结果表明：2晶胞仿生管在轴向和斜向载荷下的耐撞性最优。同工况条件下，减少晶胞数可降低仿生管峰值载荷。斜向冲击载荷下，仿生管保持稳定叠缩变形模式的时间随晶胞数的增加而缩短，其耐撞性能随晶胞数的增加而降低。虾螯结构特征与普通圆管的结合有效提高了仿虾螯结构多晶胞薄壁管的耐撞性能。Abstract: In order to improve the crashworthiness of thin-walled tubes, the multi-cell bionic thin-walled tubes based on a shrimp chela structure were designed by the principle of structural bionics. By taking the cell number (2−6) and the impact angle (0°, 10°, 20°, 30°) as experimental factors, the finite element method was used to simulate the crashworthiness of the bionic tubes, the reliability of the results by the simulation test was verified by the drop-weight tests. The results show that the two-cell bionic tube has the best crashworthiness under axial and oblique loads. Under the same working conditions, the reduction of the number of unit cells can reduce the peak loads of the bionic tubes. Under the oblique impact load, the time for the bionic tubes to maintain the stable collapse deformation mode is shortened with the increase of the number of the cells, and the crashworthiness of the bionic tubes decreases with the increase of the number of the cells. The combination of a shrimp cheek structure and an ordinary circular tube effectively improves the crashworthiness of the designed structures. So it can provide a reference for the design of energy-absorbing components in vehicles.
-
Key words:
- thin-walled tube /
- structural bionic /
- crashworthiness /
- drop-hammer test /
- cell
-
提高吸能元件的耐撞性能可有效提高车辆在碰撞事故中的安全性，薄壁吸能结构在冲击载荷下通过折叠、弯曲等塑性变形方式可有效吸收和耗散冲击动能，其耐撞性受材料、结构、受载方式等因素影响[1]。除普通方管、圆管外，研究者们设计出了蜂窝铝、多晶胞薄壁管、功能梯度厚度管等具有优异耐撞性能的薄壁结构[2]。王博等[3]、郝文乾等[4]在普通方管管壁上引入折角，设计出预折纹管，有效地降低了管的峰值载荷，减小了冲击力的波动幅度，增大了管的有效塑性变形体积，提高了管的吸能效率。韩会龙等[5]设计了具有负泊松比特性的节点层级蜂窝材料，研究发现，较正方形蜂窝，其动态承载能力和能量吸收能力明显提高。李松晏等[6]研究了由泡沫铝填充的不同圆角半径和厚度的带圆角方管的耐撞性，发现填充了泡沫铝的大圆角半径的厚管缓冲吸能效果明显增强。Najafi等[7]研究了多种具有角元素的多胞元薄壁铝管在轴向冲击载荷下的压溃吸能特性，推导出了平均挤压力的计算公式，并分析了胞元结构对管冲击特性的影响。Fang等[8]向多单元管中引入功能梯度厚度（functionally graded thickness，FGT）的结构特征，发现FGT对多胞管的耐撞性有显著的影响。
实际碰撞事故中，吸能结构很少发生单一轴向或径向压溃，多为斜向压溃，斜向载荷会使管的压溃稳定性和吸能性能显著降低。因此，研究者们研发了多种抗斜向冲击薄壁结构。Li等[9]比较了斜向冲击载荷下FGT管、锥形均匀厚度(tapered uniform thickness，TUT)管和直均匀厚度(straight uniform thickness，SUT)管的吸能性能，发现梯度指数为2的FGT管在斜向冲击工况下的吸能性能最优，且随冲击载荷角度的增大，较TUT管和SUT管，FGT管仍保持良好的吸能性能。亓昌等[10]分析了轴对称锥形多胞薄壁方管的长径比、壁厚和锥度对其斜向冲击吸能特性的影响，并拟合出了可用于预测其在斜向冲击下比吸能和冲击峰值力的解析公式。高强等[11]设计了新型椭圆形泡沫填充管，采用有限元法分析了椭圆向心率、壁厚和泡沫铝密度等参数对其斜向冲击吸能特性的影响，发现减小向心率可降低其峰值冲击力并提高其比吸能，峰值冲击力随着泡沫密度的增大而增大。
近年来，随着工程仿生技术的发展[12]，通过结构仿生已设计了多种仿生薄壁吸能管。邹猛等[13]基于牛角结构特征，设计了仿牛角结构锥形薄壁管，仿真结果表明其比吸能比普通锥管高1.8倍。受竹材截面微观结构和竹节启发，Song等[14]设计并研究了不同的仿竹结构薄壁结构，发现竹截面结构可有效提高薄壁结构的耐撞性能。Yin等[15]等基于木贼类植物的茎设计了6种不同截面结构的仿生管，采用有限元法研究了在轴向及径向载荷下6种管的耐撞性能，发现仿生管的壁厚、内壁直径及晶胞数对其耐撞性能有显著影响。
本文中，将以仿虾螯结构薄壁管为研究对象，对不同晶胞数的仿生管进行轴向和斜向载荷工况下的仿真分析，重点探讨其斜向冲击载荷下的变形模式和吸能特性。
1. 仿虾螯薄壁管设计
1.1 仿生原型分析
虾螯是雀尾螳螂虾的捕食器官，外形呈锤型，其锤部截面呈纺锤状，由前部类水晶物质层、后部平行纤维结构层和连接两部分的中部软组织层3部分组成[16]，如图1所示。捕食时螳螂虾用前锤部以高达23 m/s的速度在任意方向敲击猎物外壳，瞬时冲击力超过700 N，纺锤形结构使虾螯有效地吸收来自各方向的冲击动能，具有良好的吸能特性，避免了结构的损坏[17]。
分析虾螯的受载条件后发现，其与薄壁管有以下相似性：
（1）载荷相似性。车辆受到外部撞击时，薄壁管可能受到任意方向的冲击载荷。同样，螳螂虾用虾螯击碎猎物躯壳时，虾螯也可能承受来自任意方向的冲击载荷。
（2）功能相似性。薄壁管的功能是承受外部的冲击载荷并吸收冲击动能，而虾螯的功能也是承受冲击猎物硬壳时的冲击载荷并吸收冲击动能。
因此，选取虾螯锤部作为仿生原型，对薄壁管的结构进行优化设计。
1.2 仿生管设计
提取虾螯截面结构外轮廓线后发现其外轮廓由前部椭圆弧与后部两条切线组成，考虑到轻量化设计和加工成本，将虾螯简化为前部椭圆弧和后部2条切线。如图2所示，将简化后的结构应用于普通圆管（circular tube，CT）的截面中，由于虾螯前部和后部由中部软组织层连接，因此增加中间管作为连接管，并增加内管连接后部，形成仿虾螯截面的晶胞结构。当普通圆管管长与管壁厚之比小于50且管长与管直径之比小于3时，其发生轴向压溃时可维持稳定的环形渐进叠缩变形模式[2]。考虑实际应用，将具有仿虾螯截面晶胞结构的仿生管（bionic tube, BT）的长度H、外管直径D1、壁厚t分别定为160、70、1.5 mm。作为中间连接管，外管直径D1与中间管直径D2之差和中间管直径D2与内管直径D3之差相等，合理设计后将中间管直径D2与内管直径D3分别定为44 mm和18 mm，如图3所示。
根据晶胞单元数目将BT分为2、3、4、5和6晶胞管，分别以BT-2、BT-3、BT-4、BT-5和BT-6表示，以长度、壁厚、外径相同的CT作为对照管。因下文所采用的耐撞性指标为比吸能，计算时管的质量被消掉，所以不再对不同晶胞数的BT进行等质量化处理。各BT的截面结构如图4所示。
图 4 具有不同晶胞数的仿生管(BT)的截面结构图Figure 4. Cross-sectional structures of bionic tubes (BTs) with different cell numbers2. 仿虾螯薄壁管仿真分析
2.1 有限元模型建立
在Hypermesh软件中建立2～6晶胞的BT及CT的有限元模型。管的材料选择为AA6061铝合金，其密度ρ为2.73 × 103 kg/m3，弹性模量E为71 GPa，泊松比μ为0.3，屈服强度σs为240 MPa，伸长率σb为19%。赋予模型对应AA6061铝合金的MATL24材料特性，模型有限元网格边长为2 mm[6-8]。在BT的自由端设置质量为200 kg、冲击速度为10 m/s的刚性墙（试验室落锤试验台的落锤质量为186 kg且下落高度为10 m，其最大冲击动能远大于所设置的刚性墙的冲击动能，可通过调节落锤下落高度控制其冲击速度以使落锤冲击动能与刚性墙冲击动能相同），在BT的固定端施加约束以固定BT。刚性墙冲击角度α设定为0°（轴向）、10°、20°和30°[18]。将导出的K文件导入LS-DYNA中进行仿真计算。进行全因素仿真试验的结果表明，较对BT其他位置进行轴向和斜向冲击，正对BT晶胞位置进行轴向和斜向冲击时，BT的耐撞性能最优（耐撞性评价指标及加权评价方法在下文给出）。因此，下文只描述了正对晶胞位置冲击的仿真试验结果分析。仿真试验中，刚性墙的冲击方向和冲击角度与对BT晶胞冲击的位置如图5所示。
2.2 耐撞性指标
衡量吸能结构耐撞性的参数[19]主要包括：
(1)比吸能(specific energy absorption), Ea,s
Ea,s是单位质量的薄壁管吸收的能量，是衡量能量吸收能力的重要参数，其表达式为：
Ea,s=Ea/m (1) 式中：m为薄壁管的总质量，Ea为冲击过程中薄壁管吸收的能量。
(2)平均载荷Fm和峰值载荷Fp
Fm是整个冲击过程中冲击载荷的平均值，其表达式为：
Fm=Ea/x (2) 式中：x为冲击过程中薄壁管的有效变形位移。
Fp是压溃过程中的最大载荷，反映薄壁管的最大过载情况。设计薄壁管时应尽量减小峰值载荷，避免过高的冲击加速度对乘员造成生命危险。
(3)压溃力效率(crushing force efficiency), η
η是平均载荷与峰值载荷之比，其数值越大，表明薄壁管载荷波动越小，吸能效率越高，其表达式为：
η=Fm/Fp (3) 2.3 有限元模型可靠性验证
由于仿真试验中，仿生管的数值计算在理想条件下进行，而实际工况中BT的耐撞性能是否与仿真结果一致，需通过实物试验进行验证。
图6为BT有限元模型可靠性验证的方案流程图。首先确定BT的结构和参数，选用BT-4进行实物验证试验，其长度为160 mm，壁厚为1.5 mm，外管、中间管和内管直径分别为70、44和18 mm。使用软件Hypermesh建立BT-4的有限元模型，网格大小为2 mm×2 mm。所制作的BT-4样件的材料为AA6061铝合金，其尺寸与有限元模型的尺寸完全一致，如图7所示。采用落锤式冲击试验机对BT-4样件进行冲击试验，如图8所示，落锤的质量为186 kg，最大掉落高度为10 m，最大冲击速度为14 m/s。仿真与实物试验中BT-4的工况均为轴向冲击，通过设定落锤与BT-4样件的垂直高度来确定落锤与BT-4样件接触时的初始冲击速度。为保证仿真试验中刚性墙与实物试验中落锤的冲击动能相同，通过计算得到落锤与BT-4样件的垂直高度应设定为5.48 m（仿真试验中刚性墙的质量设定为200 kg，冲击速度为10 m/s），即落锤与BT-4样件接触时的初始冲击速度为10.37 m/s。
图 6 BT有限元模型可靠性验证流程图Figure 6. Reliability verification flow chart of the finite element model for BT图9为BT-4的仿真试验与实物试验的载荷（F）-位移（x）曲线对比图。从图9可发现，通过仿真试验与实物试验得到的BT-4的2条F-x曲线的变化趋势基本一致，其波峰和波谷的数量、位置及所对应的载荷值相差较小。仿真试验和实物试验得到BT-4的Fm、Fp、η分别为82.67 kN、152.35 kN、54.26%和89.23 kN、168.66 kN、52.90%，误差分别为7.3%、9.6%、2.5%。BT-4的仿真试验与实物试验的对比结果验证了BT-4的仿真试验的可靠性。
图 9 BT-4的仿真与实物试验的载荷-位移曲线对比Figure 9. Comparison of load-displacement curves between simulation and physical test for BT-42.4 仿真结果分析
2.4.1 仿生管变形模式分析
在压溃过程中，薄壁管通过发生塑性变形来吸收冲击动能，其结构的设计应保证其受冲击时发生稳定的渐进叠缩变形，且避免发生欧拉弯曲变形。随斜向冲击角度的增大，薄壁管总经历由渐进叠缩变形模式转变为欧拉弯曲变形模式，最终发生坍塌变形的过程，而转变起始角度越大，发生坍塌时间越晚，其斜向冲击承受范围越广[18]。
图10为2～6晶胞的BT和CT在不同冲击角度下的变形模式。由图10可知，轴向冲击和10°角冲击时，BT截面的仿生结构使其受到的冲击载荷均匀分布，随着冲击位移的增大，折叠数量逐渐增多，各仿生管均发生稳定的渐进叠缩变形模式。但轴向冲击时，CT上端发生渐进叠缩变形的同时，管身发生倾斜，且底端也出现折叠变形，变形稳定性差。10°角冲击时，由于BT-3特殊的截面结构，在压溃中发生2次大幅度叠缩变形。20°角冲击时，BT-4、BT-5、BT-6均由渐进叠缩变形转变为欧拉弯曲变形，最终发生坍塌变形，导致其能量吸收量大幅降低，而BT-2、BT-3则呈现渐进叠缩变形模式。30°角冲击时，所有管均经历渐进叠缩变形、欧拉弯曲变形及坍塌变形的过程，且BT的晶胞越多，发生坍塌的时间越早。
图 10 不同冲击角度（由左至右分别为0°、10°、20°和30°）下2～6晶胞的BT和CT的变形模式Figure 10. Deformation modes of the BTs with 2−6 cells and CTs under impact angles of 0°, 10°, 20° and 30° from left to right2.4.2 仿生管吸能特性分析
表1列出了2～6晶胞BT和CT在不同角度(α)冲击载荷下的Ea、m及Ea,s的仿真试验值。图11为2～6晶胞的BT和CT在不同冲击角度下的比吸能（Ea,s）-位移（x）曲线，x取100 mm。
表 1 2～6晶胞BT和CT在不同角度冲击载荷下的Ea、m及Ea,s的仿真试验值Table 1. Simulation test values of Ea, m and Ea,s of BTs with 2−6 cells and CTs under different-angle impact loadsα/(°) 薄壁管 Ea/kJ m/kg Ea,s/(kJ·kg−1) α/(°) 薄壁管 Ea/kJ m/kg Ea,s/(kJ·kg−1) 0 BT-2 6.01 0.223 10 26.96 20 BT-2 4.92 0.223 10 22.06 BT-3 7.11 0.247 10 28.77 BT-3 5.85 0.247 10 23.66 BT-4 8.28 0.271 20 30.53 BT-4 5.68 0.271 20 20.93 BT-5 8.07 0.295 20 31.61 BT-5 3.41 0.295 20 11.56 BT-6 10.29 0.319 20 32.23 BT-6 4.31 0.319 20 13.51 CT 1.50 0.094 99 15.77 CT 1.02 0.094 99 10.70 10 BT-2 5.42 0.223 10 24.28 30 BT-2 2.01 0.223 10 9.00 BT-3 6.50 0.247 10 26.30 BT-3 2.08 0.247 10 8.43 BT-4 7.54 0.271 20 27.82 BT-4 2.05 0.271 20 7.54 BT-5 8.58 0.295 20 29.08 BT-5 2.39 0.295 20 8.09 BT-6 9.51 0.319 20 29.78 BT-6 2.53 0.319 20 7.92 CT 1.33 0.094 99 14.00 CT 0.63 0.094 99 6.67 
图 11 不同冲击角度下2～6晶胞的BT和CT的比吸能-位移曲线Figure 11. Specific energy absorption-displacement curves of the BTs with 2−6 cells and CTs under different impact angles由表1和图11可知，不同冲击工况下，各晶胞数的BT的Ea,s均高于CT的。轴向和10°角冲击时：2～6晶胞的BT的Ea,s较接近，但远高于CT的Ea,s；所有管中，BT-6的Ea,s最大，分别为32.23 kJ/kg和29.78 kJ/kg，较BT中Ea,s最小的BT-2分别提高19.55%和22.65%，但较CT分别提高了104.37%、112.71%。可见，受轴向和10°角冲击时，BT具有优异的吸能特性。20°角冲击时，BT-2、BT-3在有效变形范围内始终保持良好的变形稳定性，未发生欧拉弯曲变形，表现出良好的吸能能力，其Ea,s分别为22.06 kJ/kg和23.66 kJ/kg，较CT分别提高了106.17%和121.12%。而BT-4、BT-5和BT-6分别在压缩位移为54、38和40 mm时发生欧拉弯曲变形，其Ea,s-x曲线斜率明显减小，表明其吸能能力大幅降低。CT虽未发生欧拉弯曲变形，但其Ea,s仅为10.7 kJ/kg，吸能能力较弱，导致其Ea,s明显低于其他BT，表现出较差的吸能特性。30°角冲击时，各管均经历了渐进叠缩变形模式转变为欧拉弯曲变形模式的过程，吸能能力显著降低，Ea,s-x曲线的斜率均先增大后减小。此时各BT的Ea,s仍均高于CT的，其中BT-2的比吸能最高，为9 kJ/kg，较CT的提高了34.93%。
表2列出了2～6晶胞BT和CT在不同角度冲击载荷下的Fm、Fp及η的仿真试验值。图12为2～6晶胞的BT和CT在不同冲击角度下的载荷（F）-位移（x）曲线，x取100 mm。
表 2 2～6晶胞BT与CT在不同角度冲击载荷下的Fm、Fp及η的仿真试验值Table 2. Simulation test values of Fm, Fp and η of BTs with 2−6 cells and CTs under different-angle impact loadsα/(°) 薄壁管 Fm/kN Fp/kN η/% α/(°) 薄壁管 Fm/kN Fp/kN η/% 0 BT-2 59.90 123.07 48.67 20 BT-2 51.64 66.30 77.88 BT-3 69.39 137.52 50.46 BT-3 60.37 84.94 71.08 BT-4 82.41 152.35 54.09 BT-4 65.19 97.72 66.71 BT-5 95.09 167.26 56.85 BT-5 53.92 105.04 51.33 BT-6 101.88 181.88 56.02 BT-6 64.09 121.03 52.95 CT 15.30 57.82 26.47 CT 12.51 19.64 63.68 10 BT-2 54.82 67.57 81.13 30 BT-2 29.91 59.52 50.25 BT-3 54.89 76.82 71.46 BT-3 31.28 67.28 46.48 BT-4 77.15 92.93 83.02 BT-4 32.27 67.49 47.82 BT-5 85.19 103.23 82.52 BT-5 38.38 85.04 45.13 BT-6 97.42 114.13 85.37 BT-6 42.18 93.64 45.04 CT 14.55 21.86 66.52 CT 9.06 17.98 50.37 
图 12 不同冲击角度下2～6晶胞的BT和CT的载荷-位移曲线Figure 12. Force-displacement curves of the BTs with 2−6 cells and CTs under different impact angles由表2和图12可知，各冲击角度下，BT的Fp和Fm均高于CT的，且随晶胞数的增加而增大。这是由于BT的横截面具有多晶胞结构，压溃过程中参与塑性变形的结构增多，使其抵抗冲击变形的能力增强。0°角（轴向）冲击时BT-5和10°角冲击时BT-6的η最高，分别为56.85%和85.37%，分别较同工况下的CT、BT-2的高114.77%和28.34%、16.8%和5.2%，但BT-5、BT-6较BT-2的Fp和Fm过高。20°角冲击时，BT-2的η最高，为77.88%，比CT的高22.3%。30°角冲击时，CT的η在所有薄壁管中最高，为50.37%；而BT中BT-2的η最高，BT-6的最低，分别为50.25%、45.04%；但CT的η仅比BT-2和BT-6的η分别高0.2%和11.8%，BT依然呈现良好的耐撞性能。不同冲击角度时，BT-2的载荷波动最平稳。10°角冲击时，BT-3由于其特殊的截面结构，在50～70 mm与80～100 mm两段位移中发生大幅度折叠变形，导致载荷出现大的波动，且在20°冲击时，其载荷波动最大。
3. 仿生管最优结构选取
在不同角度的冲击工况下，5种BT的耐撞性能各存在优势角度区和劣势角度区，其耐撞性能受多个评价指标的综合影响，因此，以变形模式和多个评价指标的加权组合的权值来综合评价不同结构的BT的耐撞性能，获取最优结构的BT。
BT由渐进叠缩变形模式转变为欧拉弯曲模式的区间长度受冲击角度和转变时间的综合影响，BT所能承受的最大冲击角度越大、转变时间越长，其承受的斜向冲击角度区间越大。为探求不同晶胞数结构的BT由渐进叠缩变形转变为欧拉弯曲变形和坍塌变形的极限角度与转变区间，在(10～30)°斜向冲击角度内，以2°为间隔进一步开展仿真试验。
图13为各晶胞数BT在相同压溃位移内的最大比吸能（Ea,s,max）- 冲击角度（α）曲线。由图13可知，冲击角度小于20°时，BT-2的Ea,s,max-α的曲线斜率保持不变。这是由于冲击角度小于20°时，BT-2保持稳定的渐进叠缩变形模式，在相同压溃时间内，其Ea,s,max随冲击角度的增大而减小。冲击角度为(20～24)°时，曲线斜率迅速降低，这是由于BT-2由渐进叠缩变形模式逐渐转变为欧拉弯曲模式，且最终发生坍塌变形。欧拉弯曲变形模式下，BT的能量吸收量降低，较渐进叠缩变形模式时，其Ea,s,max大幅减小。冲击角度大于24°时，BT-2总逐步发生3种变形，Ea,s,max随冲击角度的增大而减小。而BT-3、BT-4、BT-5和BT-6的转变区间分别为(20～22)°、(18～22)°、(18～20)°和(18～20)°，因此，BT-2的斜向变形稳定性最优，且随晶胞的增加，BT斜向变形稳定性降低。
图 13 各晶胞数BT的最大比吸能-冲击角度曲线Figure 13. Maximum specific energy absorption-impact angle curves of BTs with different number of cells由于5种BT的耐撞性指标各存在优势角度区和劣势角度区，因此将评估指标定义为Fp、Ea,s和η的加权组合，以全面考虑这3种指标对不同晶胞数BT耐撞性能的影响，加权公式为：
W=δ1Ea,sEa,s,max−δ2FpFp,max+δ3ηηmax (4) 式中：δ1 + δ2 + δ3 = 1。因实际应用中，
Ea,s 为衡量薄壁结构耐撞性优劣的主要指标，且Fp和η为另外2个重要指标，因此将δ1 = 0.4，δ2 = δ3 = 0.3作为3个评价指标的权重因子[14]。由于各评价指标的单位不一致，因此进行无量纲化，分别将Ea,s /Ea,s,max 、Fp/Fp,max和η/ηmax作为Ea,s 、Fp和η的无量纲形式。图14为所有BT在0°、10°、20°和30°角碰撞条件下的W值。由图14可知，轴向冲击时，BT-2的权值略低于其他BT，但在10°、20°及30°角冲击时，BT-2的权值均明显高于其他BT。
图 14 不同冲击角度下各晶胞数BT的权值WFigure 14. Weights of BTs with different number of the cells at different impact angles综上所述，BT-2在轴向和斜向冲击载荷下的变形模式在所有BT中最稳定，当冲击角度为0°、10°和20°时，BT-2在压溃过程中始终保持稳定的渐进叠缩变形模式，未发生欧拉弯曲变形。当冲击角度为30°时，较其他晶胞数的BT，其发生欧拉弯曲变形的时间点最晚。观察BT-2的压溃过程发现，晶胞数量少可保证BT-2在受轴向及斜向冲击时有足够的空间进行稳定的渐进叠缩变形且发生变形干涉的管壁少于其他仿生管，这也大大降低了BT-2在压溃过程中的Fp和载荷波动值，进而提高了BT-2的η。同时，较其他BT，BT-2的横截面曲线结构少，使其在受斜向冲击时维持渐进叠缩变形模式的时间长于其他BT，并有效避免了由于上端无法继续变形而导致下端发生欧拉弯曲变形的情况，这使其受斜向冲击时仍保持良好的吸能能力。因此，综合评价各晶胞数BT在轴向及斜向冲击载荷下的变形模式、吸能特性及压溃载荷值，确定2晶胞结构为仿虾螯薄壁管的最佳结构。
4. 结 论
（1）轴向冲击时，仿生管呈现出稳定的渐进叠缩变形模式。斜向冲击时，随角度的增加，仿生管由渐进叠缩变形模式转变为欧拉弯曲变形模式，最终发生坍塌变形，转变起始角度随晶胞数的增加而减小，转变区间缩短。
（2）仿真试验结果表明，仿生管的耐撞性能明显优于普通圆管。2晶胞仿生管的耐撞性能最优，在轴向和斜向冲击载荷下，其变形模式最稳定，载荷效率最高，载荷波动最小，在受0°、10°、20°、30°角冲击时的比吸能分别为26.96 kJ/kg （0°）、24.28 kJ/kg （10°）、22.06 kJ/kg （20°）和9 kJ/kg （30°），较普通圆管分别提高了67.2%、65.6%、92.5%和35.3%
（3）通过实物试验验证了仿真结果的可靠性，仿真结果和实物试验结果证明了将虾螯的结构特征应用到吸能薄壁管的设计中，可有效地提高薄壁管的吸能特性和耐撞性能，为车辆吸能元件的设计提供了参考。
-
图 4 具有不同晶胞数的仿生管(BT)的截面结构图
Figure 4. Cross-sectional structures of bionic tubes (BTs) with different cell numbers
图 6 BT有限元模型可靠性验证流程图
Figure 6. Reliability verification flow chart of the finite element model for BT
图 9 BT-4的仿真与实物试验的载荷-位移曲线对比
Figure 9. Comparison of load-displacement curves between simulation and physical test for BT-4
图 10 不同冲击角度（由左至右分别为0°、10°、20°和30°）下2～6晶胞的BT和CT的变形模式
Figure 10. Deformation modes of the BTs with 2−6 cells and CTs under impact angles of 0°, 10°, 20° and 30° from left to right
图 11 不同冲击角度下2～6晶胞的BT和CT的比吸能-位移曲线
Figure 11. Specific energy absorption-displacement curves of the BTs with 2−6 cells and CTs under different impact angles
图 12 不同冲击角度下2～6晶胞的BT和CT的载荷-位移曲线
Figure 12. Force-displacement curves of the BTs with 2−6 cells and CTs under different impact angles
图 13 各晶胞数BT的最大比吸能-冲击角度曲线
Figure 13. Maximum specific energy absorption-impact angle curves of BTs with different number of cells
图 14 不同冲击角度下各晶胞数BT的权值W
Figure 14. Weights of BTs with different number of the cells at different impact angles
表 1 2～6晶胞BT和CT在不同角度冲击载荷下的Ea、m及Ea,s的仿真试验值
Table 1. Simulation test values of Ea, m and Ea,s of BTs with 2−6 cells and CTs under different-angle impact loads
α/(°) 薄壁管 Ea/kJ m/kg Ea,s/(kJ·kg−1) α/(°) 薄壁管 Ea/kJ m/kg Ea,s/(kJ·kg−1) 0 BT-2 6.01 0.223 10 26.96 20 BT-2 4.92 0.223 10 22.06 BT-3 7.11 0.247 10 28.77 BT-3 5.85 0.247 10 23.66 BT-4 8.28 0.271 20 30.53 BT-4 5.68 0.271 20 20.93 BT-5 8.07 0.295 20 31.61 BT-5 3.41 0.295 20 11.56 BT-6 10.29 0.319 20 32.23 BT-6 4.31 0.319 20 13.51 CT 1.50 0.094 99 15.77 CT 1.02 0.094 99 10.70 10 BT-2 5.42 0.223 10 24.28 30 BT-2 2.01 0.223 10 9.00 BT-3 6.50 0.247 10 26.30 BT-3 2.08 0.247 10 8.43 BT-4 7.54 0.271 20 27.82 BT-4 2.05 0.271 20 7.54 BT-5 8.58 0.295 20 29.08 BT-5 2.39 0.295 20 8.09 BT-6 9.51 0.319 20 29.78 BT-6 2.53 0.319 20 7.92 CT 1.33 0.094 99 14.00 CT 0.63 0.094 99 6.67 
下载: 导出CSV
表 2 2～6晶胞BT与CT在不同角度冲击载荷下的Fm、Fp及η的仿真试验值
Table 2. Simulation test values of Fm, Fp and η of BTs with 2−6 cells and CTs under different-angle impact loads
α/(°) 薄壁管 Fm/kN Fp/kN η/% α/(°) 薄壁管 Fm/kN Fp/kN η/% 0 BT-2 59.90 123.07 48.67 20 BT-2 51.64 66.30 77.88 BT-3 69.39 137.52 50.46 BT-3 60.37 84.94 71.08 BT-4 82.41 152.35 54.09 BT-4 65.19 97.72 66.71 BT-5 95.09 167.26 56.85 BT-5 53.92 105.04 51.33 BT-6 101.88 181.88 56.02 BT-6 64.09 121.03 52.95 CT 15.30 57.82 26.47 CT 12.51 19.64 63.68 10 BT-2 54.82 67.57 81.13 30 BT-2 29.91 59.52 50.25 BT-3 54.89 76.82 71.46 BT-3 31.28 67.28 46.48 BT-4 77.15 92.93 83.02 BT-4 32.27 67.49 47.82 BT-5 85.19 103.23 82.52 BT-5 38.38 85.04 45.13 BT-6 97.42 114.13 85.37 BT-6 42.18 93.64 45.04 CT 14.55 21.86 66.52 CT 9.06 17.98 50.37
下载: 导出CSV
-
[1] BARTCZAK B, GIERCZYCKA-ZBROZEK D, GRONOSTAJSKI Z, et al. The use of thin-walled sections for energy absorbing components: a review [J]. Archives of Civil and Mechanical Engineering, 2010, 10(4): 5–19. DOI: 10.1016/S1644-9665(12)60027-2. [2] BAROUTAJI A, SAJJIA M, OLABI A G. On the crashworthiness performance of thin-walled energy absorbers: recent advances and future developments [J]. Thin-Walled Structures, 2017, 118: 137–163. DOI: 10.1016/j.tws.2017.05.018. [3] 王博, 周才华, 由衷. 预折纹管在低速冲击载荷作用下的能量吸收 [J]. 爆炸与冲击, 2015, 35(4): 473–481. DOI: 10.11883/1001-1455(2015)04-0473-09.WANG B, ZHOU C H, YOU Z. Energy absorption of pre-folded origami under low speed impact [J]. Explosion and Shock Waves, 2015, 35(4): 473–481. DOI: 10.11883/1001-1455(2015)04-0473-09. [4] 郝文乾, 卢进帅, 黄睿, 等. 轴向冲击载荷下薄壁折纹管的屈曲模态与吸能 [J]. 爆炸与冲击, 2015, 35(3): 380–385. DOI: 10.11883/1001-1455-(2015)03-0380-06.HAO W Q, LU J S, HUANG R, et al. Buckling and energy absorption properties of thin-walled corrugated tubes under axial impacting [J]. Explosion and Shock Waves, 2015, 35(3): 380–385. DOI: 10.11883/1001-1455-(2015)03-0380-06. [5] 韩会龙, 张新春, 王鹏. 负泊松比蜂窝材料的动力学响应及能量吸收特性 [J]. 爆炸与冲击, 2019, 39(1): 013103. DOI: 10.11883/bzycj-2017-0281.HAN H L, ZHANG X C, WANG P. Dynamic responses and energy absorption properties of honeycombs with negative Poisson’s ratio [J]. Explosion and Shock Waves, 2019, 39(1): 013103. DOI: 10.11883/bzycj-2017-0281. [6] 李松晏, 郑志军, 虞吉林. 高速列车吸能结构设计和耐撞性分析 [J]. 爆炸与冲击, 2015, 35(2): 164–170. DOI: 10.11883/1001-1455(2015)02-0164-07.LI S Y, ZHENG Z J, YU J L. Energy-absorbing structure design and crashworthiness analysis of high-speed trains [J]. Explosion and Shock Waves, 2015, 35(2): 164–170. DOI: 10.11883/1001-1455(2015)02-0164-07. [7] NAJAFI A, RAIS-ROHANI M. Mechanics of axial plastic collapse in multi-cell, multi-corner crush tubes [J]. Thin-Walled Structures, 2011, 49(1): 1–12. DOI: 10.1016/j.tws.2010.07.002. [8] FANG J G, GAO Y K, SUN G Y, et al. Dynamic crashing behavior of new extrudable multi-cell tubes with a functionally graded thickness [J]. International Journal of Mechanical Sciences, 2015, 103: 63–73. DOI: 10.1016/j.ijmecsci.2015.08.029. [9] LI G Y, XU F X, SUN G Y, et al. A comparative study on thin-walled structures with functionally graded thickness (FGT) and tapered tubes withstanding oblique impact loading [J]. International Journal of Impact Engineering, 2015, 77: 68–83. DOI: 10.1016/j.ijimpeng.2014.11.003. [10] 亓昌, 董方亮, 杨姝, 等. 锥形多胞薄壁管斜向冲击吸能特性仿真研究 [J]. 振动与冲击, 2012, 31(24): 102–107. DOI: 10.3969/j.issn.1000-3835.2012.24.021.QI C, DONG F L, YANG S, et al. Energy-absorbing characteristics of a tapered multi-cell thin-walled tube under oblique impact [J]. Journal of Vibration and Shock, 2012, 31(24): 102–107. DOI: 10.3969/j.issn.1000-3835.2012.24.021. [11] 高强, 王良模, 王源隆, 等. 椭圆形泡沫填充薄壁管斜向冲击吸能特性仿真研究 [J]. 振动与冲击, 2017, 36(2): 201–206. DOI: 10.13465/j.cnki.jvs.2017.02.033.GAO Q, WANG L M, WANG Y L, et al. Energy-absorbing characteristics of foam-filled oval tubes under oblique impact [J]. Journal of Vibration and Shock, 2017, 36(2): 201–206. DOI: 10.13465/j.cnki.jvs.2017.02.033. [12] 任露泉, 梁云虹. 仿生学导论[M]. 北京: 科学出版社, 2016: 32−65. [13] 邹猛, 于用军, 张荣荣, 等. 仿牛角结构薄壁管吸能特性仿真分析 [J]. 吉林大学学报(工学版), 2015, 45(6): 1863–1868. DOI: 10.13229/j.cnki.jdxbgxb201506020.ZOU M, YU Y J, ZHANG R R, et al. Simulation analysis of energy-absorption properties of thin-walled tube based on horn structure [J]. Journal of Jilin University (Engineering and Technology Edition), 2015, 45(6): 1863–1868. DOI: 10.13229/j.cnki.jdxbgxb201506020. [14] SONG J F, XU S C, WANG H X, et al. Bionic design and multi-objective optimization for variable wall thickness tube inspired bamboo structures [J]. Thin-Walled Structures, 2018, 125: 76–88. DOI: 10.1016/j.tws.2018.01.010. [15] YIN H F, XIAO Y Y, WEN G L, et al. Crushing analysis and multi-objective optimization design for bionic thin-walled structure [J]. Materials and Design, 2015, 87: 825–834. DOI: 10.1016/j.matdes.2015.08.095. [16] WEAVER J C, MILLIRON G W, MISEREZ A, et al. The stomatopod dactyl club: a formidable damage-tolerant biological hammer [J]. Science, 2012, 336(6086): 1275–1280. DOI: 10.1126/science.1218764. [17] PATEK S N, CALDWELL R L. Extreme impact and cavitation forces of a biological hammer: strike forces of the peacock mantis shrimp odontodactylus scyllarus [J]. Journal of Experimental Biology, 2005, 208(19): 3655–3664. DOI: 10.1242/jeb.01831. [18] SADJAD P, MOHAMMAD-HOSSEIN E, SOBHAN E M. Crashworthiness of double-cell conical tubes with different cross sections subjected to dynamic axial and oblique loads [J]. Journal of Central South University, 2018, 25(3): 632–645. DOI: 10.1007/s11771-018-3766-z. [19] CHEN B C, ZOU M, LIU G M, et al. Experimental study on energy absorption of bionic tubes inspired by bamboo structures under axial crushing [J]. International Journal of Impact Engineering, 2018, 115: 48–57. DOI: 10.1016/j.ijimpeng.2018.01.005. 期刊类型引用(16)
1. 滕晨皓，刘坤澎，宋志博，陈家玉，周才华，涂昕昊，朱冬妹. 斜向低速冲击载荷下多边形旋转拼接折纸管耐撞性研究. 应用数学和力学. 2025(05): 633-649 . 
百度学术2. 杨欣，王瑞祥，许述财，秦豪毅，宋家锋. 仿枪虾虾螯纵向混合负泊松比结构设计与耐撞性研究. 重庆理工大学学报(自然科学). 2024(01): 77-86 . 百度学术3. 张旭，赵锐，石云波，祖凯旋，张洁宇，闫晓朋. 仿虾螯式抗冲击微加速度传感器设计. 机械强度. 2024(02): 410-415 . 百度学术4. 杨欣，王阳，宋家锋，朱勇，黄彬兵，许述财. 基于虾螯结构的仿生夹层板设计及数值模拟. 吉林大学学报(工学版). 2024(03): 842-851 . 百度学术5. 周志宏，马娇，彭金艳，康睿超，杨胜兵. 梯度吸能盒耐撞性研究. 机械设计. 2024(12): 122-128 . 百度学术6. 胡敬坤，徐鹏，范志强，李耀宙，谭晓丽. 仿平行叶脉次级肋板增强多胞管轴向压溃响应. 振动与冲击. 2023(01): 38-45 . 百度学术7. 覃凌云，杨书仪，戴巨川，汪国胜，李桂兵，凌启辉，陈哲吾. 爆炸冲击下乘员下肢保护装置仿生设计与优化. 兵器装备工程学报. 2023(03): 95-103 . 百度学术8. 邓敏杰，刘志芳. 仿马尾草薄壁结构的设计与耐撞性研究. 高压物理学报. 2022(03): 111-120 . 百度学术9. 何强，王勇辉，顾航. 爆炸载荷作用下仿玉莲脉络加筋板的动态响应分析. 江苏科技大学学报(自然科学版). 2022(02): 15-21 . 百度学术10. 贺璞，邓庆田，李新波. 层合多孔圆柱壳的轴向冲击吸能特性. 高压物理学报. 2022(04): 145-154 . 百度学术11. 罗耿，赵剑南，陈亮，柴成鹏，王童，陈轶嵩. 新型波纹结构吸能盒耐撞性研究. 汽车安全与节能学报. 2022(03): 453-462 . 百度学术12. 牛枞，黄晗，向枳昕，闫庆昊，陈金宝，许述财. 仿生多胞薄壁管耐撞性分析及优化. 爆炸与冲击. 2022(10): 169-182 . 
本站查看13. 黄晗，许述财，陈姮. 仿生波纹夹层结构耐撞性分析及优化. 爆炸与冲击. 2021(08): 36-46 . 本站查看14. 于鹏山，刘志芳，李世强. 新型仿竹薄壁圆管的设计与吸能特性分析. 高压物理学报. 2021(05): 104-114 . 百度学术15. 徐峰祥，邵天如，牛小强. 薄壁圆管轴向压溃特性仿真分析与理论预测. 塑性工程学报. 2021(12): 207-215 . 百度学术16. 霍鹏，许述财，范晓文，李建平，杨欣，黄晗. 鹿角骨单位仿生薄壁管斜向冲击耐撞性研究. 爆炸与冲击. 2020(11): 127-138 . 本站查看其他类型引用(17)
-
下载:
百度学术
计量
- 文章访问数: 5258
- HTML全文浏览量: 2482
- PDF下载量: 143
- 被引次数: 33