类向日葵夹芯圆柱壳径向冲击数值模拟

上一篇

下一篇

闫栋, 王根伟, 宋辉, 王彬. 类向日葵夹芯圆柱壳径向冲击数值模拟[J]. 高压物理学报, 2020, 34(5): 054201-1. doi: 10.11858/gywlxb.20190858
引用本文: 闫栋, 王根伟, 宋辉, 王彬. 类向日葵夹芯圆柱壳径向冲击数值模拟[J]. 高压物理学报, 2020, 34(5): 054201-1. doi: 10.11858/gywlxb.20190858
Dong YAN, Genwei WANG, Hui SONG, Bin WANG. Numerical Simulation of Radial Impact on Sunflower-Like Sandwich Cylindrical Shell[J]. gaoyawli, 2020, 34(5): 054201-1. doi: 10.11858/gywlxb.20190858
Citation: Dong YAN, Genwei WANG, Hui SONG, Bin WANG. Numerical Simulation of Radial Impact on Sunflower-Like Sandwich Cylindrical Shell[J]. gaoyawli, 2020, 34(5): 054201-1. doi: 10.11858/gywlxb.20190858

类向日葵夹芯圆柱壳径向冲击数值模拟

    作者简介: 闫 栋(1992-),男,硕士研究生,主要从事金属薄壁夹芯结构的耐撞性研究. E-mail: 1690617047@qq.com .
    通讯作者: 王根伟(1974-),男,博士,副教授,主要从事新能源汽车安全与轻量化研究. E-mail: gwang@tyut.edu.cn
  • 中图分类号: O347

Numerical Simulation of Radial Impact on Sunflower-Like Sandwich Cylindrical Shell

    Corresponding author: Genwei WANG, gwang@tyut.edu.cn
  • MSC: O347

  • 摘要: 金属薄壁结构由于其优异的轻质性和耐撞性,一直被广泛地应用在汽车、飞机和火车等交通工具的碰撞动能耗散系统中。以一种类向日葵薄壁夹芯吸能结构为研究对象,研究了在瓣尖压和瓣间压两种径向冲击载荷下,类向日葵薄壁夹芯结构的变形模式、能量吸收能力、比吸能和平均压缩力。结果表明,类向日葵薄壁夹芯结构的壁厚、花瓣数、加载速度以及加载方向都会对结构的耐撞性产生一定的影响。在质量恒定条件下,随着外壳厚度的增加,瓣尖压冲击方式下薄壁结构的吸能效率降低,瓣间压比瓣尖压的比吸能最高多出了44.6%。随着花瓣数的变化,金属薄壁结构的吸能效率存在一个最优值。
  • 加载中
  • 图 1  类向日葵夹芯结构有限元模型

    Figure 1.  Finite element model of sunflower sandwich structure

    图 2  两种加载方式下不同位移对应的变形云图

    Figure 2.  Deformation nephogram under two loading modes at different displacements

    图 3  不同冲击速度下的力-位移曲线

    Figure 3.  Load-displacement curves at different impact velocities

    图 4  不同冲击速度、不同外壳壁厚条件下类向日葵夹芯结构的SEA和MCF

    Figure 4.  SEA and MCF of sunflower sandwich structure with different wall thicknesses at different impact velocities

    图 5  不同冲击速度下的载荷-位移曲线

    Figure 5.  Load-displacement curves at different impact velocities

    图 6  不同冲击速度、不同花瓣数下类向日葵夹芯结构的SEA和MCF

    Figure 6.  SEA and MCF of sunflower sandwich structure with various numbers of petals at different impact velocities

    表 1  冲击速度v = 20 m/s下的SEA和MCF

    Table 1.  SEA and MCF at v = 20 m/s

    t/mmSEA/(J·g−1)MCF/NIncrease ratio/%
    TSPGSPTSPGSP
    0.82.894.4594.40145.3435.05
    0.92.684.0887.43133.2334.31
    1.02.664.4686.90 145.8440.36
    1.12.314.1775.33136.1844.60
    1.22.283.5174.42114.7935.04
    下载: 导出CSV

    表 2  冲击速度v = 30 m/s下的SEA和MCF

    Table 2.  SEA and MCF at v = 30 m/s

    t/mmSEA/(J·g−1)MCF/NIncrease ratio/%
    TSPGSPTSPGSP
    0.83.074.83100.26157.6636.44
    0.92.874.6493.87 151.7238.15
    1.02.844.8192.76 157.2540.96
    1.12.634.6785.89 152.4243.68
    1.22.463.7480.29 122.1434.22
    下载: 导出CSV

    表 3  冲击速度v = 50 m/s下的SEA和MCF

    Table 3.  SEA and MCF at v = 50 m/s

    t/mmSEA/(J·g−1)MCF/NIncrease ratio/%
    TSPGSPTSPGSP
    0.84.085.27133.24172.0622.58
    0.93.835.08125.23165.8524.61
    1.03.615.40 117.93176.2833.15
    1.13.434.79112.05156.4428.39
    1.23.394.67110.81152.4127.41
    下载: 导出CSV
  • [1] ALEXANDER J. An approximate analysis of the collapse of thin cylindrical shells under axial loading [J]. The Quarterly Journal of Mechanics and Applied Mathematics, 1960, 13(1): 10–15. doi: 10.1093/qjmam/13.1.10
    [2] WIERZBICKI T, ABRAMOWICZ W. On the crushing mechanics of thin-walled structures [J]. Journal of Applied Mechanics, 1983, 50(4a): 727–734. doi: 10.1115/1.3167137
    [3] GUPTA N, PRASAD G E, GUPTA S. Plastic collapse of metallic conical frusta of large semi-apical angles [J]. International Journal of Crash Worthiness, 1997, 2(4): 349–366. doi: 10.1533/cras.1997.0054
    [4] HANSSEN A G, LANGSETH M, HOPPERSTAD O S. Static crushing of square aluminium extrusions with aluminium foam filler [J]. International Journal of Mechanical Sciences, 1999, 41(8): 967–993. doi: 10.1016/S0020-7403(98)00064-2
    [5] HANSSEN A G, LANGSETH M, HOPPERSTAD O S. Static and dynamic crushing of circular aluminium extrusions with aluminium foam filler [J]. International Journal of Impact Engineering, 2000, 24(5): 475–507. doi: 10.1016/S0734-743X(99)00170-0
    [6] 赵凯, 卢国兴, 沈建虎, 等. 圆环列系统吸能特性研究 [J]. 北京大学学报, 2007, 65(3): 312–316. ZHAO K, LU G X, SHEN J H, et al. Energy absorption characteristics of ring train system [J]. Journal of Peking University, 2007, 65(3): 312–316.
    [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] 李志斌, 虞吉林, 郭刘伟. 具有诱导结构的铝合金薄壁方管轴向压缩吸能性能试验研究 [J]. 工程力学, 2012, 29(6): 346–352. doi: 10.6052/j.issn.1000-4750.2010.09.0663 LI Z B, YU J L, GUO L W. Experimental study on energy absorption performance of aluminum alloy thin-walled square tube with induced structure in axial compression [J]. Engineering Mechanics, 2012, 29(6): 346–352. doi: 10.6052/j.issn.1000-4750.2010.09.0663
    [9] 杨鹏飞. 波纹夹芯板和结构的压缩与冲击吸能特性研究 [D]. 哈尔滨: 哈尔滨工程大学, 2013. YANG P F. Study on compression and impact energy absorption characteristics of corrugated sandwich plates and structures [D]. Harbin: Harbin Engineering University, 2013.
    [10] 陈亦涛, 王帅, 刘凯欣. 径向冲击载荷下金属圆环系统的动态吸能特性 [C]//中国力学会第19届学术年会. 北京, 2013.
    [11] 于渤, 张钱城, 金峰, 等. 泡沫铝填充波纹板的动态压缩性能研究 [C]//中国力学大会-2013. 西安, 2013.
    [12] TARLOCHAN F, SAMER F, HAMOUDA A M S, et al. Design of thin wall structures for energy absorption applications: enhancement of crashworthiness due to axial and oblique impact forces [J]. Thin-Walled Structures, 2013, 71: 7–17. doi: 10.1016/j.tws.2013.04.003
    [13] 杨彬彬, 赵修平. 多胞金属管受径向冲击时的吸能特性 [J]. 海军航空工程学院学报, 2014, 29(6): 552–556. doi: 10.7682/j.issn.1673-1522.2014.06.010 YANG B B, ZHAO X P. Energy absorption characteristics of polycellular metal tube subjected to radial impact [J]. Journal of Naval Aeronautical Engineering College, 2014, 29(6): 552–556. doi: 10.7682/j.issn.1673-1522.2014.06.010
    [14] WANG J, ZHANG Y, HE N, et al. Crashworthiness behavior of Koch fractal structures [J]. Materials & Design, 2018, 144: 229–244.
    [15] 韩宾. 波纹强化复合型多孔材料的力学行为研究 [D]. 西安: 西安交通大学, 2018. HAN B. Study on mechanical behavior of corrugated reinforced composite porous materials [D]. Xi’an: Xi’an Jiaotong University, 2018.
    [16] NGOC S H, LU G X, XIANG X M. Energy absorption of a bio-inspired honeycomb sandwich panel [J]. Journal of Materials Science, 2019, 54(8): 6286–6300. doi: 10.1007/s10853-018-3163-x
  • 加载中
图( 6) 表( 3)
计量
  • 文章访问数:  2413
  • HTML全文浏览数:  2413
  • PDF下载数:  19
  • 施引文献:  0
出版历程
  • 收稿日期:  2019-11-18
  • 刊出日期:  2020-10-05

类向日葵夹芯圆柱壳径向冲击数值模拟

    通讯作者: 王根伟(1974-),男,博士,副教授,主要从事新能源汽车安全与轻量化研究. E-mail: gwang@tyut.edu.cn
    作者简介: 闫 栋(1992-),男,硕士研究生,主要从事金属薄壁夹芯结构的耐撞性研究. E-mail: 1690617047@qq.com
  • 1. 太原理工大学应用力学研究所,山西 太原 030024
  • 2. 材料强度与结构冲击山西省重点实验室,山西 太原 030024
  • 3. 伦敦布鲁内尔大学机械与航空工程系,英国 伦敦 UB8 3PH

摘要: 金属薄壁结构由于其优异的轻质性和耐撞性,一直被广泛地应用在汽车、飞机和火车等交通工具的碰撞动能耗散系统中。以一种类向日葵薄壁夹芯吸能结构为研究对象,研究了在瓣尖压和瓣间压两种径向冲击载荷下,类向日葵薄壁夹芯结构的变形模式、能量吸收能力、比吸能和平均压缩力。结果表明,类向日葵薄壁夹芯结构的壁厚、花瓣数、加载速度以及加载方向都会对结构的耐撞性产生一定的影响。在质量恒定条件下,随着外壳厚度的增加,瓣尖压冲击方式下薄壁结构的吸能效率降低,瓣间压比瓣尖压的比吸能最高多出了44.6%。随着花瓣数的变化,金属薄壁结构的吸能效率存在一个最优值。

English Abstract

  • 金属薄壁结构的耐撞性研究在人身和设备安全防护方面方兴未艾。在汽车等交通工具发生的碰撞事件中,薄壁结构或构件对人或汽车实施保护的能力就是结构的耐撞性。在碰撞过后,如果金属薄壁结构及其人或汽车的损伤越小,则结构的耐撞性就越高。在实际的工程应用中,由于安全保护的需要,公路上的护栏及各种海陆空交通工具如汽车、飞机、轮船等,都必须满足一定的耐撞性要求,须依靠薄壁结构的塑性变形来吸收缓冲碰撞时的能量,确保人员安全和重要器材的完好。我国载人航天飞船返回着陆时人体防护研究中的一项重要内容就是改善宇航员座椅的冲击吸能性能。近年来由于交通事故所造成的人员伤亡不断增加,汽车在碰撞下的安全引起了全世界的高度关注,所以研究薄壁结构的耐撞性显得尤为重要。

    金属薄壁结构的耐撞性研究的核心问题在于能量的吸收。在结构的设计过程中一般需要遵循以下基本原则。首先,所设计的薄壁结构和所用的材料能够满足大部分冲击动能通过薄壁结构的塑性变形或其他方式的能量耗散转换为非弹性能,而不是以弹性能量的方式储存下来。其次,能量吸收器的初始峰值力应当低于一个临界值。在金属薄壁结构的大变形过程中,理想状况下的反作用力应当尽可能的保持恒定,是一个近似矩形的力-位移曲线特性。再者,力所作的功等于其大小乘以沿力作用线发生的位移,所以金属薄壁结构要吸收的能量应该足够大,行程足够长,同时冲程效率要足够高。这是因为当吸能结构受到冲击载荷开始变形时,一直到结构变实,这段距离是吸能效率最高的一段,当结构由实到完全密实化,有效冲程很短,触头受到的反作用力非常大,对能量的吸收是无效的,而且反复冲击会对触头产生很大的损害,影响触头的使用寿命。最后,为了应对实际工程中不稳定的工作载荷,所设计的金属薄壁结构的变形模式和能量吸收能力应当是稳定、可重复的,确保金属薄壁结构在实际复杂工程应用条件下的可靠性和实用性。

    薄壁构件作为汽车车身的主要吸能装置,学者们很早就开始了研究。1960年,Alexander[1]首先建立了金属薄壁模型在轴向加载下的塑性变形吸能理论,为金属薄壁结构的耐撞性研究做出了巨大的贡献。在Alexander研究的基础上,Wierzbicki等[2]创建了一个新的超折叠单元分析理论,他们还发现了金属薄壁结构在管壁发生大变形的情况下,折叠管壁是曲线而不是直线。研究者也对薄壁结构横截面形状进行了大量的研究。Gupta等[3]研究了长径比对金属薄壁结构吸能效率的影响,通过试验得到了多种破坏模式、载荷-变形行为以及初始峰值力和平均压缩载荷。Hanssen等[4-5]对泡沫填充管的吸能特性进行了理论研究,并做了大量准静态和动态轴向压缩实验,验证了泡沫填充管良好的吸能效率。赵凯等[6]对圆环在径向冲击下的吸能特性进行了模拟和实验研究,结果表明:圆环壁厚越薄,其变形越大,圆环的吸能效率越高;多圆环系统会使单个圆环的吸能效率降低,但是随着冲击速度的增大,多个圆环系统与单个圆环的吸能效率差值逐渐减小。Najafi等[7]利用有限元模拟和实验研究了不同多胞、多角构型的轴向冲击金属薄壁结构的吸能特性。李志斌等[8]对设置了诱导结构的铝合金管进行了轴向冲击实验,发现诱导结构可以降低结构响应的初始峰值力。杨鹏飞[9]进行了波纹板的压缩实验,通过变化波纹板的厚度、高度以及波间距,得到了承载能力最优的模型参数。陈亦涛等[10]对不同排列方式的圆环系统进行了实验研究,得出:在抗力允许的范围内,减小圆环的径厚比有助于提高多圆环系统的比吸能;随着圆环层数的增加,正方形排列的圆环系统的比吸能降低,而六边形排列的圆环系统的比吸能几乎没有变化,说明六边形排列的圆环系统比正方形排列的圆环系统更具有工程应用价值。于渤等[11]对泡沫铝填充波纹板的动态压缩性能做了数值模拟研究,在不同冲击速度下,泡沫铝填充波纹板展现出了不同的压缩性能,具体划分为准静态模式、中间模式和冲击波模式。在准静态模式下,由于填充泡沫所产生的约束,使波纹板的抗弯曲能力得到提高;在中间模式下,波纹板的波长由于惯性作用减小,泡沫的约束作用不明显;在冲击模式下,高速冲击使泡沫内部产生了冲击波,导致波纹板不发生屈曲。Tarlochan等[12]研究了圆管、方管以及多边形管的金属薄壁结构在动态冲击载荷下的耐撞性。随着社会的发展和科技的进步,由于实验和理论研究的局限性,学者利用有限元模拟做了大量的金属薄壁构件的耐撞性研究。杨彬彬等[13]研究了传统多胞金属管和夹层多胞金属管受径向冲击时的吸能特性,结果表明:在两种类型的多胞管中,夹层多胞管的吸能特性整体较优;与夹层多胞管相比,传统多胞管在缓冲后出现了特别明显的回弹现象,对触头产生很大的反作用力,影响触头的使用寿命。Wang等[14]研究了不同科赫分形阶数对吸能效率的影响,通过实验测试和数值模拟,发现科赫分形结构与外壳之间所产生的协同作用极大地提高了薄壁结构的吸能效率,并且得出二阶科赫分形设计具有更好的吸能效率。韩宾[15]对闭孔泡沫铝和不锈钢波纹复合夹芯板进行了面外压缩实验研究,发现泡沫填充使复合芯体的面外压缩强度以及能量吸收能力都得到了极大的提升,此外复合芯体中的蜂窝泡沫铝与波纹板之间的强耦合作用,使得复合芯体和蜂窝泡沫铝的承载能力和能量吸收效率都大大的提高。Ngoc等[16]根据啄木鸟的喙设计出了波纹夹芯结构,并且进行了轴向压缩研究,取得了很好的吸能效果。本研究将波纹板弯曲形成类向日葵结构,放入两个同心圆柱壳之间,组成双圆柱壳夹芯结构,即类向日葵夹芯薄壁圆柱壳结构,分析此结构在径向冲击载荷下的吸能特性。

  • 在本研究中,提出了一种类向日葵夹芯金属薄壁吸能结构,利用有限元分析软件ABAQUS对薄壁构件进行了模拟。有限元模型如图1所示。在薄壁构件的上下两端,有两个质量均为6 kg的刚性板,薄壁构件质量m = 2.45 g恒定,网格为实体单元。为了保证计算精度,最大单元网格控制为0.3 mm × 0.3 mm。在外壳中层直径为90 mm、内壳中层直径为30 mm的条件下,通过改变壁厚(t = 0.8、0.9、1.0、1.1、1.2 mm)和花瓣数(9、10、11、12、13、14、15、16瓣),采用图1中的瓣尖压(Tip side pressure,TSP)和瓣间压(Gap side pressure,GSP)两种方式进行径向冲击数值模拟。试件材料为铝合金6061-O,密度ρ = 2 700 kg/m3,泊松比μ = 0.3,杨氏模量E = 68.2 GPa。铝合金6061-O的应力-应变曲线来自文献[14]中的图4。

  • 尽管薄壁吸能结构各不相同,有圆管、多边形管、多胞管、填充泡沫管以及复合材料管等,研究人员也从不同的角度、不同内容,运用理论分析、数值模拟以及实验验证等不同方法研究薄壁结构的吸能特性,但最终都有同样的参考指标,对比不同薄壁结构的吸能特性,综合各方面的优势,选出最佳的薄壁吸能结构。一般用总吸能(EA)、平均压缩力(MCF)、比吸能(SEA)等评价薄壁结构的耐撞性能。

    EA是压缩过程中吸收的总能量,反映了在整个压缩过程中吸收能量的大小,表达式为

    式中:Et为总吸能;d为有效变形位移,在本研究中为75 mm;F(x)为瞬时压缩载荷。

    MCF代表薄壁结构的平均压缩强度,反映了在整个压缩过程中压缩载荷的大小,其表达式为

    式中:Fave为平均压缩力。

    SEA为单位质量所吸收的比能,反映了薄壁结构的吸能效率,是区分不同材料吸能能力的关键指标,其表达式为

    式中:e为比吸能,M为结构的总质量。

  • 在保持相同质量的前提下,改变外壳壁厚tt为0.8、0.9、1.0、1.1、1.2 mm),花瓣壁厚等厚,内壳壁厚随外壳壁厚的增加而减小。采用瓣尖压和瓣间压两种加载方式,以不同冲击速度(v为20、30、50 m/s)对薄壁构件进行径向压缩,得到了两种加载方式下不同位移对应的变形云图,如图2所示。

    图3给出了类向日葵夹芯结构(花瓣数相同、外壳壁厚不同)在不同冲击速度下的载荷-位移曲线。从图3可以看出,当刚性板接触到薄壁构件时,瓣尖压首先出现一个峰值,而瓣间压则逐渐增大,达到一个峰值。这是因为在瓣尖压加载方式下,刚性板首先接触到花瓣尖端,尖端为主要承载部分,而瓣间压是花瓣尖端之间承受冲击力,在压缩到两个尖端开始承受力时,出现一个峰值。在瓣尖压加载方式下,出现第二个峰值是因为下端的花瓣尖端开始成为主要承载部分,当最下端花瓣开始屈曲,曲率达到最大的时候,力减到最小,最后,吸能元件逐渐密实化。瓣间压在达到峰值后,力逐渐减小,是因为下端的花瓣也开始变形,最后力趋于一个稳定的波动范围,直到密实化。相对于瓣尖压,瓣间压的力的波动性更小,更符合理想化的“矩形”载荷-位移曲线。

    图4给出了类向日葵夹芯结构(花瓣数相同、外壳壁厚不同)在不同冲击速度下的SEA和MCF,可见,外壳厚度的变化对SEA和MCF有显著的影响。瓣尖压加载下,SEA随着外壳壁厚的增大而减小,这是因为外壳壁厚越小,外壳进入塑性的区域越大,外壳质量占整个构件质量的比例比内壳高,材料利用率也更高。瓣尖压加载下的初始峰值力大于瓣间压加载,所以瓣间压模式对车辆的保护更好。随着冲击速度的增大,SEA和MCF都增大,瓣间压加载下的SEA比瓣尖压加载下高,在t = 1.1 mm、v = 20 m/s时,瓣间压下的SEA比瓣尖压下高44.6%。表1表2表3列出了3种冲击速度下的SEA和MCF。

  • 考虑到花瓣个数可能会影响薄壁构件的吸能效率,本研究采用瓣尖压和瓣间压两种径向加载模式,以不同冲击速度(v = 20、30、50 m/s)对8种类向日葵夹芯结构,花瓣数分别为9、10、11、12、13、14、15、16,进行径向压缩模拟计算。

    图5给出了类向日葵夹芯结构(花瓣数不同)在不同冲击速度下的载荷-位移曲线。从图5可以看出,在瓣间压径向加载下,花瓣数越多,载荷越先到达峰值。这是因为花瓣数越多,花瓣之间的角度越小,相邻的两个花瓣越早开始承受载荷,最先达到峰值。随着冲击速度增大,初始峰值力增大,载荷-位移曲线的波动性更明显,并且瓣尖压径向加载下载荷-位移曲线的波动比瓣间压径向加载下更大。这是因为瓣尖压加载下,上下两个花瓣在受到冲击时,屈曲部分变形的抖动比较剧烈,而瓣间压加载下花瓣变形相对稳定。

    图6给出了类向日葵夹芯结构(花瓣数不同)在不同冲击速度下的SEA和MCF。从图6可以看出:在冲击速度v = 20 m/s的条件下,不同花瓣数的SEA和MCF之间的差距很大,花瓣数为13的SEA和MCF都是最小的;随着冲击速度的增大,花瓣数为15的SEA和MCF增加明显;在冲击速度v = 50 m/s下,瓣间压加载下的SEA和MCF都是最高的。同时从图6还可以看出,随着花瓣数的增加,瓣尖压和瓣间压径向加载下的吸能效率差值减小。这是因为随着花瓣数的增加,花瓣间的角度逐渐减小,瓣尖压和瓣间压两种加载方式对薄壁构件吸能效率的影响减小。

  • 本研究提出了一种新型的金属薄壁结构——类向日葵夹芯圆柱结构,探讨了具有不同外壳壁厚和花瓣数的类向日葵夹芯结构在瓣尖压和瓣间压两种径向冲击载荷下的吸能特性。

    (1)随着冲击速度的增大,初始峰值力增大,SEA和MCF都增大。

    (2)SEA的大小与外壳壁厚的变化有关,在质量一定的前提下,随着外壳壁厚的增加,瓣尖压加载下的SEA逐渐减小。

    (3)SEA和MCF的大小与类向日葵夹芯结构的花瓣数有关。在质量一定的前提下,随着花瓣数的增加,瓣尖压和瓣间压两种加载方式对薄壁构件吸能效率的影响减小;随着冲击速度的增大,不同花瓣数的薄壁构件的吸能效果逐渐呈现规律性,在某一花瓣数下,SEA和MCF达到最大。

参考文献 (16)

目录

/

返回文章
返回