• ISSN 1001-1455  CN 51-1148/O3
  • EI Compendex、CA收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊

仿虾螯结构薄壁管设计及耐撞性分析

杨欣 范晓文 许述财 黄晗 霍鹏

引用本文:
Citation:

仿虾螯结构薄壁管设计及耐撞性分析

    作者简介: 杨 欣(1974- ),男,博士,教授,yangxin@hebau.edu.cn;
    通讯作者: 许述财, xushc@tsinghua.edu.cn
  • 中图分类号: O382;U463.9

Design and crashworthiness analysis of thin-walled tubes based on a shrimp chela structure

    Corresponding author: XU Shucai, xushc@tsinghua.edu.cn ;
  • CLC number: O382;U463.9

  • 摘要: 为提高薄壁管的耐撞性能,以虾螯为生物原型,通过结构仿生原理设计了仿虾螯结构多晶胞薄壁管。以晶胞数(2~6)和冲击角度(0°、10°、20°、30°)为试验因素,利用有限元法分析了仿虾螯结构多晶胞薄壁管在不同冲击角度下的耐撞性能,通过落锤试验验证了仿真结果的可靠性。结果表明:2晶胞仿生管在轴向和斜向载荷下的耐撞性最优。同工况条件下,减少晶胞数可降低仿生管峰值载荷。斜向冲击载荷下,仿生管保持稳定叠缩变形模式的时间随晶胞数的增加而缩短,其耐撞性能随晶胞数的增加而降低。虾螯结构特征与普通圆管的结合有效提高了仿虾螯结构多晶胞薄壁管的耐撞性能。
  • 图 1  虾螯截面结构

    Figure 1.  The structure of the cross section of the chela

    图 2  晶胞结构

    Figure 2.  The structure of a cell

    图 3  仿生管的结构及尺寸

    Figure 3.  The structure and size of a bionic tube

    图 4  具有不同晶胞数的仿生管(BT)的截面结构图

    Figure 4.  Cross-sectional structures of bionic tubes (BTs) with different cell numbers

    图 5  仿生管的有限元模型

    Figure 5.  The finite element model for the bionic tube

    图 6  BT有限元模型可靠性验证流程图

    Figure 6.  Reliability verification flow chart of the finite element model for BT

    图 7  BT-4样件

    Figure 7.  Samples of BT-4

    图 8  落锤冲击试验机

    Figure 8.  Drop hammer impact tester

    图 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/kJm/kgEa,s/(kJ·kg−1)α/(°)薄壁管Ea/kJm/kgEa,s/(kJ·kg−1)
    0BT-2 6.010.223 1026.9620BT-24.920.223 1022.06
    BT-3 7.110.247 1028.77BT-35.850.247 1023.66
    BT-4 8.280.271 2030.53BT-45.680.271 2020.93
    BT-5 8.070.295 2031.61BT-53.410.295 2011.56
    BT-610.290.319 2032.23BT-64.310.319 2013.51
    CT 1.500.094 9915.77CT1.020.094 9910.70
    10BT-2 5.420.223 1024.2830BT-22.010.223 10 9.00
    BT-3 6.500.247 1026.30BT-32.080.247 10 8.43
    BT-4 7.540.271 2027.82BT-42.050.271 20 7.54
    BT-5 8.580.295 2029.08BT-52.390.295 20 8.09
    BT-6 9.510.319 2029.78BT-62.530.319 20 7.92
    CT 1.330.094 9914.00CT0.630.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/kNFp/kNη/%α/(°)薄壁管Fm/kNFp/kNη/%
    0BT-2 59.90123.0748.6720BT-251.64 66.3077.88
    BT-3 69.39137.5250.46BT-360.37 84.9471.08
    BT-4 82.41152.3554.09BT-465.19 97.7266.71
    BT-5 95.09167.2656.85BT-553.92105.0451.33
    BT-6101.88181.8856.02BT-664.09121.0352.95
    CT 15.30 57.8226.47CT12.51 19.6463.68
    10BT-2 54.82 67.5781.1330BT-229.91 59.5250.25
    BT-3 54.89 76.8271.46BT-331.28 67.2846.48
    BT-4 77.15 92.9383.02BT-432.27 67.4947.82
    BT-5 85.19103.2382.52BT-538.38 85.0445.13
    BT-6 97.42114.1385.37BT-642.18 93.6445.04
    CT 14.55 21.8666.52CT 9.06 17.9850.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.
  • [1] 李松晏郑志军虞吉林 . 高速列车吸能结构设计和耐撞性分析. 爆炸与冲击, 2015, 35(2): 164-170. doi: 10.11883/1001-1455(2015)02-0164-07
    [2] 徐峰祥张锁武昆迎 . 厚度幂指数分布管状结构耐撞性设计准则与方法研究. 爆炸与冲击, 2019, 39(3): 035103-1-035103-9. doi: 10.11883/bzycj-2018-0013
    [3] 常白雪郑志军赵凯何思渊虞吉林 . 具有恒定冲击载荷的梯度泡沫金属材料设计. 爆炸与冲击, 2019, 39(4): 041101-1-041101-9. doi: 10.11883/bzycj-2018-0431
    [4] 殷之平李玉龙黄其青 . 含诱导缺陷薄壁圆管耐撞性优化设计. 爆炸与冲击, 2011, 31(4): 418-422. doi: 10.11883/1001-1455(2011)04-0418-05
    [5] 谭丽辉徐涛崔晓梅张炜赵世佳 . 带有圆弧形凹槽金属薄壁圆管抗撞性优化设计. 爆炸与冲击, 2014, 34(5): 547-553. doi: 10.11883/1001-1455(2014)05-0547-07
    [6] 任冀宾王斌王振刘军索涛李玉龙 . 某型飞机机翼前缘抗鸟撞结构设计与试验验证. 爆炸与冲击, 2019, 39(2): 025101-1-025101-9. doi: 10.11883/bzycj-2017-0407
    [7] 章冠人 . 落锤试验中炸药样品表面能的多重分形性质. 爆炸与冲击, 1996, 16(1): 1-4.
    [8] 王世英胡焕性 . B炸药装药发射安全性落锤模拟加载实验研究. 爆炸与冲击, 2003, 23(3): 275-278.
    [9] 刘龙飞周强 . 表面粗糙度对6061铝合金薄壁管冲击膨胀断裂性能的影响. 爆炸与冲击, 2018, 38(4): 749-758. doi: 10.11883/bzycj-2016-0389
    [10] 李永池朱林法胡秀章袁福平曹结东 . 粘塑性薄壁管中复合应力波的传播特性研究. 爆炸与冲击, 2003, 23(1): 1-5.
    [11] 叶想平李英雷张祖根 . 45钢薄壁圆柱管的冲击膨胀断裂机理研究. 爆炸与冲击, 2014, 34(3): 322-327. doi: 10.11883/1001-1455(2014)03-0322-06
    [12] 张涛吴英友朱显明刘土光 . 多边形截面薄壁管撕裂卷曲吸能研究. 爆炸与冲击, 2007, 27(3): 223-229. doi: 10.11883/1001-1455(2007)03-0223-07
    [13] 郝文乾卢进帅黄睿刘志芳王志华 . 轴向冲击载荷下薄壁折纹管的屈曲模态与吸能. 爆炸与冲击, 2015, 35(3): 380-385. doi: 10.11883/1001-1455-(2015)03-0380-06
    [14] 张振华钱海峰王媛欣牟金磊梅志远牛闯 . 球头落锤冲击下金字塔点阵夹芯板结构的动态响应实验. 爆炸与冲击, 2015, 35(6): 888-894. doi: 10.11883/1001-1455(2015)06-0888-07
    [15] 王丙斌王蕊 . 空心率对中空夹层钢管混凝土组合柱耐撞性能影响. 爆炸与冲击, 2018, 38(1): 204-211. doi: 10.11883/bzycj-2016-0143
    [16] 姚小虎韩强张晓晴赵隆茂 . 飞机圆弧风挡抗鸟撞试验研究. 爆炸与冲击, 2005, 25(5): 417-422. doi: 10.11883/1001-1455(2005)05-0417-06
    [17] 路国运秦斌张国权韩志军雷建平 . 冲击作用下夹层充液薄壁半球壳组合结构的动力响应. 爆炸与冲击, 2012, 32(6): 561-567. doi: 10.11883/1001-1455(2012)06-0561-07
    [18] 毛勇建李玉龙 . SHPB试验中试件的轴向应力均匀性. 爆炸与冲击, 2008, 28(5): 448-454. doi: 10.11883/1001-1455(2008)05-0448-07
    [19] 龚汉生陈远龄 . 对落锤仪中落锤运动的监测与分析. 爆炸与冲击, 1992, 12(4): 357-362.
    [20] 毋玲郭英男李玉龙 . 蜂窝夹芯雷达罩结构的鸟撞数值分析. 爆炸与冲击, 2009, 29(6): 642-647. doi: 10.11883/1001-1455(2009)06-0642-06
  • 加载中
图(14)表(2)
计量
  • 文章访问数:  1709
  • HTML全文浏览量:  1520
  • PDF下载量:  47
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-07-15
  • 录用日期:  2019-11-19
  • 网络出版日期:  2020-04-02
  • 刊出日期:  2020-04-01

仿虾螯结构薄壁管设计及耐撞性分析

    作者简介:杨 欣(1974- ),男,博士,教授,yangxin@hebau.edu.cn
    通讯作者: 许述财, xushc@tsinghua.edu.cn
  • 1. 河北农业大学机电工程学院,河北 保定 071000
  • 2. 清华大学汽车安全与节能国家重点实验室,北京 100084

摘要: 为提高薄壁管的耐撞性能,以虾螯为生物原型,通过结构仿生原理设计了仿虾螯结构多晶胞薄壁管。以晶胞数(2~6)和冲击角度(0°、10°、20°、30°)为试验因素,利用有限元法分析了仿虾螯结构多晶胞薄壁管在不同冲击角度下的耐撞性能,通过落锤试验验证了仿真结果的可靠性。结果表明:2晶胞仿生管在轴向和斜向载荷下的耐撞性最优。同工况条件下,减少晶胞数可降低仿生管峰值载荷。斜向冲击载荷下,仿生管保持稳定叠缩变形模式的时间随晶胞数的增加而缩短,其耐撞性能随晶胞数的增加而降低。虾螯结构特征与普通圆管的结合有效提高了仿虾螯结构多晶胞薄壁管的耐撞性能。

English Abstract

  • 提高吸能元件的耐撞性能可有效提高车辆在碰撞事故中的安全性,薄壁吸能结构在冲击载荷下通过折叠、弯曲等塑性变形方式可有效吸收和耗散冲击动能,其耐撞性受材料、结构、受载方式等因素影响[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种管的耐撞性能,发现仿生管的壁厚、内壁直径及晶胞数对其耐撞性能有显著影响。

    本文中,将以仿虾螯结构薄壁管为研究对象,对不同晶胞数的仿生管进行轴向和斜向载荷工况下的仿真分析,重点探讨其斜向冲击载荷下的变形模式和吸能特性。

    • 虾螯是雀尾螳螂虾的捕食器官,外形呈锤型,其锤部截面呈纺锤状,由前部类水晶物质层、后部平行纤维结构层和连接两部分的中部软组织层3部分组成[16],如图1所示。捕食时螳螂虾用前锤部以高达23 m/s的速度在任意方向敲击猎物外壳,瞬时冲击力超过700 N,纺锤形结构使虾螯有效地吸收来自各方向的冲击动能,具有良好的吸能特性,避免了结构的损坏[17]

      图  1  虾螯截面结构

      Figure 1.  The structure of the cross section of the chela

      分析虾螯的受载条件后发现,其与薄壁管有以下相似性:

      (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所示。

      图  2  晶胞结构

      Figure 2.  The structure of a cell

      图  3  仿生管的结构及尺寸

      Figure 3.  The structure and size of a bionic tube

      根据晶胞单元数目将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 numbers

    • 在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所示。

      图  5  仿生管的有限元模型

      Figure 5.  The finite element model for the bionic tube

    • 衡量吸能结构耐撞性的参数[19]主要包括:

      (1)比吸能(specific energy absorption), Ea,s

      Ea,s是单位质量的薄壁管吸收的能量,是衡量能量吸收能力的重要参数,其表达式为:

      ${\rm{ }}{E_{{\rm{a,s}}}} = {\rm{ }}{E_{\rm{a}}}/m$

      式中:m为薄壁管的总质量,Ea为冲击过程中薄壁管吸收的能量。

      (2)平均载荷Fm和峰值载荷Fp

      Fm是整个冲击过程中冲击载荷的平均值,其表达式为:

      ${F_{\rm{m}}}{\rm{ }} = {\rm{ }}{E_{\rm{a}}}/x$

      式中:x为冲击过程中薄壁管的有效变形位移。

      Fp是压溃过程中的最大载荷,反映薄壁管的最大过载情况。设计薄壁管时应尽量减小峰值载荷,避免过高的冲击加速度对乘员造成生命危险。

      (3)压溃力效率(crushing force efficiency), η

      η是平均载荷与峰值载荷之比,其数值越大,表明薄壁管载荷波动越小,吸能效率越高,其表达式为:

      $ \eta = {F_{\rm{m}}}/{F_{\rm{p}}} $

    • 由于仿真试验中,仿生管的数值计算在理想条件下进行,而实际工况中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

      图  7  BT-4样件

      Figure 7.  Samples of BT-4

      图  8  落锤冲击试验机

      Figure 8.  Drop hammer impact tester

      图9为BT-4的仿真试验与实物试验的载荷(F)-位移(x)曲线对比图。从图9可发现,通过仿真试验与实物试验得到的BT-4的2条F-x曲线的变化趋势基本一致,其波峰和波谷的数量、位置及所对应的载荷值相差较小。仿真试验和实物试验得到BT-4的FmFpη分别为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-4

    • 在压溃过程中,薄壁管通过发生塑性变形来吸收冲击动能,其结构的设计应保证其受冲击时发生稳定的渐进叠缩变形,且避免发生欧拉弯曲变形。随斜向冲击角度的增大,薄壁管总经历由渐进叠缩变形模式转变为欧拉弯曲变形模式,最终发生坍塌变形的过程,而转变起始角度越大,发生坍塌时间越晚,其斜向冲击承受范围越广[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 right

    • 表1列出了2~6晶胞BT和CT在不同角度(α)冲击载荷下的EamEa,s的仿真试验值。图11为2~6晶胞的BT和CT在不同冲击角度下的比吸能(Ea,s)-位移(x)曲线,x取100 mm。

      α/(°)薄壁管Ea/kJm/kgEa,s/(kJ·kg−1)α/(°)薄壁管Ea/kJm/kgEa,s/(kJ·kg−1)
      0BT-2 6.010.223 1026.9620BT-24.920.223 1022.06
      BT-3 7.110.247 1028.77BT-35.850.247 1023.66
      BT-4 8.280.271 2030.53BT-45.680.271 2020.93
      BT-5 8.070.295 2031.61BT-53.410.295 2011.56
      BT-610.290.319 2032.23BT-64.310.319 2013.51
      CT 1.500.094 9915.77CT1.020.094 9910.70
      10BT-2 5.420.223 1024.2830BT-22.010.223 10 9.00
      BT-3 6.500.247 1026.30BT-32.080.247 10 8.43
      BT-4 7.540.271 2027.82BT-42.050.271 20 7.54
      BT-5 8.580.295 2029.08BT-52.390.295 20 8.09
      BT-6 9.510.319 2029.78BT-62.530.319 20 7.92
      CT 1.330.094 9914.00CT0.630.094 99 6.67

      表 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

      图  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在不同角度冲击载荷下的FmFpη的仿真试验值。图12为2~6晶胞的BT和CT在不同冲击角度下的载荷(F)-位移(x)曲线,x取100 mm。

      α/(°)薄壁管Fm/kNFp/kNη/%α/(°)薄壁管Fm/kNFp/kNη/%
      0BT-2 59.90123.0748.6720BT-251.64 66.3077.88
      BT-3 69.39137.5250.46BT-360.37 84.9471.08
      BT-4 82.41152.3554.09BT-465.19 97.7266.71
      BT-5 95.09167.2656.85BT-553.92105.0451.33
      BT-6101.88181.8856.02BT-664.09121.0352.95
      CT 15.30 57.8226.47CT12.51 19.6463.68
      10BT-2 54.82 67.5781.1330BT-229.91 59.5250.25
      BT-3 54.89 76.8271.46BT-331.28 67.2846.48
      BT-4 77.15 92.9383.02BT-432.27 67.4947.82
      BT-5 85.19103.2382.52BT-538.38 85.0445.13
      BT-6 97.42114.1385.37BT-642.18 93.6445.04
      CT 14.55 21.8666.52CT 9.06 17.9850.37

      表 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

      图  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的FpFm均高于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的FpFm过高。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°冲击时,其载荷波动最大。

    • 在不同角度的冲击工况下,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的耐撞性指标各存在优势角度区和劣势角度区,因此将评估指标定义为FpEa,sη的加权组合,以全面考虑这3种指标对不同晶胞数BT耐撞性能的影响,加权公式为:

      $W = \frac{{{\delta _1}{E_{{\rm{a}},{\rm{s}}}}}} {{{E_{{\rm{a}},{\rm{s}},{\rm{max}}}}}} - \frac{{{\delta _2}{F_{\rm{p}}}}}{{{F_{{\rm{p,}}\max }}}} + \frac{{{\delta _3}\eta }}{{{\eta _{\max }}}}$

      式中:δ1 + δ2 + δ3 = 1。因实际应用中,${E_{{\rm{a}},{\rm{s}}}}$为衡量薄壁结构耐撞性优劣的主要指标,且Fpη为另外2个重要指标,因此将δ1 = 0.4,δ2 = δ3 = 0.3作为3个评价指标的权重因子[14]。由于各评价指标的单位不一致,因此进行无量纲化,分别将${E_{{\rm{a}},{\rm{s}}}}$/${E_{{\rm{a}},{\rm{s}},{\rm{max} }}}$Fp/Fp,maxη/ηmax作为${E_{{\rm{a}},{\rm{s}}}}$Fpη的无量纲形式。

      图14为所有BT在0°、10°、20°和30°角碰撞条件下的W值。由图14可知,轴向冲击时,BT-2的权值略低于其他BT,但在10°、20°及30°角冲击时,BT-2的权值均明显高于其他BT。

      图  14  不同冲击角度下各晶胞数BT的权值W

      Figure 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晶胞结构为仿虾螯薄壁管的最佳结构。

    • (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)通过实物试验验证了仿真结果的可靠性,仿真结果和实物试验结果证明了将虾螯的结构特征应用到吸能薄壁管的设计中,可有效地提高薄壁管的吸能特性和耐撞性能,为车辆吸能元件的设计提供了参考。

参考文献 (19)

目录

    /

    返回文章
    返回