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

弧形折纸模式薄壁结构的压缩变形与能量吸收

张天辉 邓健强 刘志芳 李世强

引用本文:
Citation:

弧形折纸模式薄壁结构的压缩变形与能量吸收

    作者简介: 张天辉(1992- ),男,硕士,zhangtianhui1234@163.com;
    通讯作者: 李世强, lishiqiang@tyut.edu.cn
  • 中图分类号: O383

Compression deformation and energy absorption of thin-walled structures with arc-shaped origami patterns

    Corresponding author: LI Shiqiang, lishiqiang@tyut.edu.cn
  • CLC number: O383

  • 摘要: 选用PolyMaxTM PLA为材料,利用3D打印技术制备了弧形折纸薄壁管件。基于准静态轴向压缩实验,运用ABAQUS软件对弧形折纸薄壁管件轴向准静态压缩和冲击行为进行了有限元计算,探讨了其变形模式和能量吸收特性,分析了预折角和薄壁单胞管件阵列数量对其压溃模式及能量吸收的影响。有限元计算结果与实验结果吻合较好。薄壁管件的变形过程可分为4个阶段:初始压溃阶段、预折角塑性旋转阶段、腹板塑性屈曲阶段和完全压溃密实化阶段。弧形折痕的引入能够有效地降低薄壁管件在压缩过程中的初始压溃载荷峰值,减小冲击载荷的振荡幅值。对比了高度相等、质量近似相等的方管与弧形折纸薄壁管在不同冲击速度下的压缩变形与能量吸收。在准静态压缩作用下,对于单胞模型,仅有折痕倾角为70°的模型的比吸能优于方管;对于多胞管件阵列模型,方管的比吸能均优于折纸管。折纸管的压缩力效率和比总体效率均优于方管,其中折痕倾角为50°的模型的压缩力效率和比总体效率最高。在动态冲击压缩下,阵列方管的比吸能均优于阵列折纸管。当冲击速度为10 mm/ms时,折纸管的压缩力效率和比总体效率均优于方管,其中折痕倾角为50°的模型的压缩力效率和比总体效率最高。当冲击速度为20 mm/ms时,仅有折痕倾角为50°的模型的压缩力效率和比总体效率优于方管。
  • 图 1  弧形折痕薄壁管建模过程

    Figure 1.  The modeling process of an origami thin-walled tube with the arc-shape pattern

    图 2  利用PolyMaxTM PLA 打印的哑铃形试件的拉伸真实应力应变关系及打印工艺

    Figure 2.  True stress-strain relation of a dumbbell-shaped sample printed with PolyMaxTM PLA and the corresponding printing process

    图 3  模拟与实验压缩载荷-位移曲线对比

    Figure 3.  Comparison of compressive load-displacement curves between simulation and experimentl

    图 4  弧形折痕薄壁管状结构单胞有限元模型

    Figure 4.  The finite element model of a curved origami thin-walled tube

    图 5  网格敏感性验证

    Figure 5.  Validation of mesh sensitivity

    图 6  有限元算法验证

    Figure 6.  Verification of the finite element algorithm

    图 7  不同单胞模型的应力-应变曲线

    Figure 7.  Stress-strain curves of different single-cell models

    图 8  面内方向阵列个数为2×2的阵列模型的应力-应变曲线

    Figure 8.  Stress-strain curves corresponding to the models with the number n of in-plane arrays equal to 2×2

    图 9  面内方向阵列个数不同的A60°模型的应力-应变曲线

    Figure 9.  Stress-strain curves of different in-plane arrays number of A60° models

    图 10  不同冲击速度下,不同阵列模型的应力-应变曲线

    Figure 10.  Stress-strain curves of different models under different impact velocities

    图 11  准静态压缩下不同的单胞薄壁管比吸能随位移的变化

    Figure 11.  Specific energy absorption varying with displacement for different single-cell thin-walled tubes under quasi-static compression

    图 12  不同冲击速度下面内方向阵列个数相同的不同薄壁管结构比吸能的对比

    Figure 12.  Comparison of specific energy absorption for different thin-walled tubes with the same number of in-plane arrays under different impact velocities

    图 13  不同薄壁管模型压缩力效率的对比

    Figure 13.  Comparison of crush force efficiencies for different thin-walled tube models

    图 14  面内方向阵列个数相同的不同薄壁管模型在相同冲击载荷下的压缩力效率随压缩位移的变化

    Figure 14.  Crush force efficiency varying with compressive displacement for different thin-walled tube models with the same number of in-plane arrays under the same impact velocity

    图 15  不同薄壁管模型比总体效率的对比

    Figure 15.  Comparison of specific total efficiencies for different thin-walled tube models

    图 16  比总体效率随压缩位移的变化

    Figure 16.  Specific total efficiency varying with compressive displacement

    表 1  模型参数

    Table 1.  The parameters of the models

    模型单胞质量/gh/mmw/mml/mmτ/mmc/mma/mmh0/mmξ /(°)α/(°)θ/(°)β/(°)
    A50°8.65304045.231.101014.2610.9260504261
    A60°8.66304041.571.201011.9810.3760605571
    A70°8.59304038.861.251010.7910.1460706778
    SQU8.66304041.891.47
    下载: 导出CSV

    表 2  不同模型的初始压溃载荷峰值

    Table 2.  Initial peak crushing loads of different models

    冲击速度/(mm·ms−1)n载荷/kN
    A50°A60°A70°SQU
    11×1 0.23 0.65 1.72 6.59
    2×2 1.72 3.91 8.4732.46
    102×2 4.92 9.0613.2242.84
    202×210.4719.7723.7246.14
    下载: 导出CSV
  • [1] LEE T U, YOU Z, GATTAS J M. Elastica surface generation of curved-crease origami [J]. International Journal of Solids and Structures, 2018, 136-137: 13–27. DOI: 10.1016/j.ijsolstr.2017.11.029.
    [2] SINGACE A A, EL-SOBKY H. Behaviour of axially crushed corrugated tubes [J]. International Journal of Mechanical Sciences, 1997, 39(3): 249–268. DOI: 10.1016/S0020-7403(96)00022-7.
    [3] ZHANG X, CHENG G D, YOU Z, et al. Energy absorption of axially compressed thin-walled square tubes with patterns [J]. Thin-Walled Structures, 2007, 45(9): 737–746. DOI: 10.1016/j.tws.2007.06.004.
    [4] ZHAO Z A, KUANG X, WU J T, et al. 3D printing of complex origami assemblages for reconfigurable structures [J]. Soft Matter, 2018, 14(39): 8051–8059. DOI: 10.1039/C8SM01341A.
    [5] MA J Y, YOU Z. Energy absorption of thin-walled square tubes with a prefolded origami pattern: Part I: geometry and numerical simulation [J]. Journal of Applied Mechanics, 2014, 81(1): 011003. DOI: 10.1115/1.4024405.
    [6] SONG J, CHEN Y, LU G X. Axial crushing of thin-walled structures with origami patterns [J]. Thin-Walled Structures, 2012, 54: 65–71. DOI: 10.1016/j.tws.2012.02.007.
    [7] ZHOU C H, ZHOU Y, WANG B. Crashworthiness design for trapezoid origami crash boxes [J]. Thin-Walled Structures, 2017, 117: 257–267. DOI: 10.1016/j.tws.2017.03.022.
    [8] 周昳鸣, 周才华, 王博. 预折纹吸能管的多样性可竞争优化设计 [J]. 振动与冲击, 2016, 35(19): 143–147. DOI: 10.13465/j.cnki.jvs.2016.19.024.
    ZHOU Y M, ZHOU C H, WANG B. Optimization design for pre-fold energy absorption tubes [J]. Journal of Vibration and Shock, 2016, 35(19): 143–147. DOI: 10.13465/j.cnki.jvs.2016.19.024.
    [9] 郝文乾, 卢进帅, 黄睿, 等. 轴向冲击载荷下薄壁折纹管的屈曲模态与吸能 [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.
    [10] HOU D G, CHEN Y, MA J Y, et al. Axial crushing of thin-walled tubes with kite-shape pattern [C] // Proceedings of ASME 2015 International Design Engineering Technical Conferences and Computers and Information in Engineering Conference. Boston, Massachusetts, USA: American Society of Mechanical Engineers, 2015: 5.
    [11] MA J Y, HOU D G, CHEN Y, et al. Quasi-static axial crushing of thin-walled tubes with a kite-shape rigid origami pattern: numerical simulation [J]. Thin-Walled Structures, 2016, 100: 38–47. DOI: 10.1016/j.tws.2015.11.023.
    [12] LI S Q, LI X, WANG Z H, et al. Sandwich panels with layered graded aluminum honeycomb cores under blast loading [J]. Composite Structures, 2017, 173: 242–254. DOI: 10.1016/j.compstruct.2017.04.037.
    [13] LI S Q, LI X, WANG Z H, et al. Finite element analysis of sandwich panels with stepwise graded aluminum honeycomb cores under blast loading [J]. Composites Part A: Applied Science and Manufacturing, 2016, 80: 1–12. DOI: 10.1016/j.compositesa.2015.09.025.
    [14] QIAO J X, CHEN C Q. In-plane crushing of a hierarchical honeycomb [J]. International Journal of Solids and Structures, 2016, 85-86: 57–66. DOI: 10.1016/j.ijsolstr.2016.02.003.
    [15] SANTOSA S P, WIERZBICKI T, HANSSEN A G, et al. Experimental and numerical studies of foam-filled sections [J]. International Journal of Impact Engineering, 2000, 24(5): 509–534. DOI: 10.1016/S0734-743X(99)00036-6.
  • [1] 郝文乾卢进帅黄睿刘志芳王志华 . 轴向冲击载荷下薄壁折纹管的屈曲模态与吸能. 爆炸与冲击, doi: 10.11883/1001-1455-(2015)03-0380-06
    [2] 林钦栋冯春唐德泓李世海殷凌云王然江 . 冲击载荷作用下路面结构的沉降及破坏特征. 爆炸与冲击, doi: 10.11883/bzycj-2018-0320
    [3] 陈晓坤李海涛王秋红金永飞邓军张嬿妮 . 瓦斯爆炸载荷作用下矿用柱壳救生舱抗爆性分析及结构优化. 爆炸与冲击, doi: 10.11883/bzycj-2016-0248
    [4] 王永虎石秀华 . 入水冲击问题研究的现状与进展. 爆炸与冲击, doi: 10.11883/1001-1455(2008)03-0276-07
    [5] 胡柳青李夕兵龚声武 . 冲击载荷作用下裂纹动态响应的数值模拟. 爆炸与冲击, doi: 10.11883/1001-1455(2006)03-0214-08
    [6] 赵铮李晓杰陶钢 . 冲击载荷下孔隙塌缩过程的数值模拟. 爆炸与冲击, doi: 10.11883/1001-1455(2009)03-0289-06
    [7] 周风华王永刚 . 影响冲击载荷下脆性材料碎片尺度的因素. 爆炸与冲击, doi: 10.11883/1001-1455(2008)04-0298-06
    [8] 邓树升耿继辉谭俊杰 . 高压储气罐泄漏喷射时冲击载荷数值计算. 爆炸与冲击, doi: 10.11883/1001-1455(2007)01-0001-06
    [9] 赵志衡汝楠马涌张超李春峰 . 强脉冲电磁力驱动的冲击载荷. 爆炸与冲击, doi: 10.11883/1001-1455(2016)05-0710-05
    [10] 曹源金先龙李政 . 冲击载荷下柔性储液罐动态响应数值模拟及规律分析. 爆炸与冲击, doi: 10.11883/1001-1455(2011)05-0469-06
    [11] 王俊奇王亮张杰 . 冲击载荷下岩石压胀对其性质的影响. 爆炸与冲击, doi: 10.11883/1001-1455(2012)03-0333-04
    [12] 张铱鈖赵隆茂 . 非均匀泡沫金属材料在冲击载荷下的变形模拟. 爆炸与冲击, doi: 10.11883/1001-1455(2006)01-0033-06
    [13] 刘红岩杨军陈鹏万 . 冲击载荷作用下岩体破坏规律的数值流形方法模拟研究. 爆炸与冲击, doi: 10.11883/1001-1455(2005)03-0255-05
    [14] 潘光杨悝 . 空投鱼雷入水载荷. 爆炸与冲击, doi: 10.11883/1001-1455(2014)05-0521-06
    [15] 敬霖王志华赵隆茂阎庆荣 . 撞击载荷下泡沫铝夹芯梁的塑性动力响应. 爆炸与冲击, doi: 10.11883/1001-1455(2010)06-0561-08
    [16] 温丽晶张春明郭超段璞张连生段卓平 . 飞机模型高速撞击钢筋混凝土载荷特性实验研究. 爆炸与冲击, doi: 10.11883/bzycj-2016-0337
    [17] 张先锋赵晓宁乔良 . 反应金属冲击反应过程的理论分析. 爆炸与冲击, doi: 10.11883/1001-1455(2010)02-0145-07
    [18] 武仁杰李海波 . SHPB冲击作用下层状千枚岩多尺度破坏机理研究. 爆炸与冲击, doi: 10.11883/bzycj-2019-0187
    [19] 钱莹杨军王兵臣 . 岩石巴西圆盘动态劈裂的流形元法模拟. 爆炸与冲击, doi: 10.11883/1001-1455(2009)01-0023-06
    [20] 周文海梁瑞陈金林朱冕陈鹏辉楼晓明王敦繁 . 时程稳定性系数确定的边坡逐孔起爆孔间微差降振时间. 爆炸与冲击, doi: 10.11883/bzycj-2018-0337
  • 加载中
图(16)表(2)
计量
  • 文章访问数:  1180
  • HTML全文浏览量:  395
  • PDF下载量:  1
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-09-15
  • 录用日期:  2019-12-19
  • 网络出版日期:  2020-06-25

弧形折纸模式薄壁结构的压缩变形与能量吸收

    作者简介:张天辉(1992- ),男,硕士,zhangtianhui1234@163.com
    通讯作者: 李世强, lishiqiang@tyut.edu.cn
  • 太原理工大学机械与运载工程学院应用力学研究所 山西 太原 030024

摘要: 选用PolyMaxTM PLA为材料,利用3D打印技术制备了弧形折纸薄壁管件。基于准静态轴向压缩实验,运用ABAQUS软件对弧形折纸薄壁管件轴向准静态压缩和冲击行为进行了有限元计算,探讨了其变形模式和能量吸收特性,分析了预折角和薄壁单胞管件阵列数量对其压溃模式及能量吸收的影响。有限元计算结果与实验结果吻合较好。薄壁管件的变形过程可分为4个阶段:初始压溃阶段、预折角塑性旋转阶段、腹板塑性屈曲阶段和完全压溃密实化阶段。弧形折痕的引入能够有效地降低薄壁管件在压缩过程中的初始压溃载荷峰值,减小冲击载荷的振荡幅值。对比了高度相等、质量近似相等的方管与弧形折纸薄壁管在不同冲击速度下的压缩变形与能量吸收。在准静态压缩作用下,对于单胞模型,仅有折痕倾角为70°的模型的比吸能优于方管;对于多胞管件阵列模型,方管的比吸能均优于折纸管。折纸管的压缩力效率和比总体效率均优于方管,其中折痕倾角为50°的模型的压缩力效率和比总体效率最高。在动态冲击压缩下,阵列方管的比吸能均优于阵列折纸管。当冲击速度为10 mm/ms时,折纸管的压缩力效率和比总体效率均优于方管,其中折痕倾角为50°的模型的压缩力效率和比总体效率最高。当冲击速度为20 mm/ms时,仅有折痕倾角为50°的模型的压缩力效率和比总体效率优于方管。

English Abstract

  • 薄壁结构作为吸能装置被广泛应用于吸能缓冲、抗冲击领域,近年来,折纸图案在管状薄壁结构上的应用越来越受到关注[1]。研究发现,如果在管件表面预制不同模式的折痕,可调控在轴向外力作用下结构的屈曲模式和能量吸收性能。在圆管表面引入波纹,可以产生较低的初始峰值力和更均匀的压缩过程[2]。Zhang等[3]通过数值模拟研究了两种具有金字塔形状的方形管,发现其能量吸收比传统方管分别提高了(15~33)%和(54~93)%。Zhao等[4]在折纸管中引入具有有限厚度和长度的铰链区域,并利用3D打印技术进行设计,这种设计充分展示了3D打印技术为增强折纸结构的承载能力和可控的变形模式提供了设计灵活性与制造效率。Ma等[5]提出了一种新的装置,该装置是根据一种特别设计的折纸图案,在方形管的表面预折而成。和传统方管相比,在该设计中,初始峰值力降低了20%以上,能量吸收增加了50%以上。Song等[6]应用了一种折纸图案做成方形管,使初始峰值力最小化。Zhou等[7]提出一种梯形折纸薄壁结构,其表面根据可展开的折纸图案预先折叠,并在方管表面引入一种梯形折痕以降低峰值力,将设计中最优的梯形折纸管与传统的等质量方管和八边形管比较,发现梯形折纸薄壁结构在吸能方面优于传统薄壁结构。周昳鸣等[8]针对这种折纹管进行了优化设计。郝文乾等[9]采用LS-DYNA软件研究了6种折纹管在轴向冲击载荷下的屈曲模式和能量吸收性能,发现折痕的引入有效地降低了初始峰值力,并且部分折纹管的压缩力效率和比总体效率高于方管的。Hou等[10]和Ma等[11]基于风筝形折纸模型,提出了一种截面为风筝形的管状结构,与传统方管相比,其峰值载荷降低了67.2%,能量吸收提升了36%。与传统的方形管相比,在保证能量吸收能力一定的情况下,具有初始折痕的薄壁管具有较低的初始峰值力和更均匀的压缩载荷。

    本文中,将折纸艺术与3D打印技术相结合,设计并制备了面外方向具有初始折叠角度的弧形折纸薄壁管件单胞及其阵列结构,基于准静态轴向压缩实验,利用ABAQUS对弧形折痕薄壁管件在轴向压缩下的力学行为和吸能特性进行了数值模拟分析,探讨了结构的变形模式,分析了预折角的大小、单胞、面内阵列个数以及不同的冲击速度对其压缩行为的影响。由于压缩行程、质量是吸能结构设计与评价的主要参量,所以本文中将高度相等、质量近似相等(结构最大质量差为0.07 g,最小质量差为0.01 g)的方管及其阵列结构进行对比,发现弧形折痕的引入能够有效地降低初始压溃载荷峰值并减小冲击载荷作用下冲击力位移曲线的振荡幅值,可以作为一种变形可控的吸能结构。

    • 本文中所研究的模型主要由一个矩形薄板沿着预定的折痕折叠而成,如图1(a)所示。其中,DE长度记为cAD=DG=GK= h0,折痕倾角∠ABE记为α,折痕BE=CF=HL=IM=a,可得:

      图  1  弧形折痕薄壁管建模过程

      Figure 1.  The modeling process of an origami thin-walled tube with the arc-shape pattern

      $ {h}_{0}=a\mathrm{sin}\alpha $

      将基本单元沿折痕BEEHHL折叠成图1(b)所示的几何形状,其中,倾斜翼缘BCFE和水平面的夹角β、腹板DEHG和腹板FEHI的夹角ξ以及BE与水平面的夹角θ由以下几何关系进行计算:

      $ \mathrm{sin}\beta =\frac{\mathrm{sin}\theta }{\mathrm{sin}\alpha } $

      $ \mathrm{cos}\frac{\xi }{2}=\frac{\mathrm{cos}\mathrm{\alpha }}{\mathrm{cos}\mathrm{\theta }} $

      因此,几何模型可根据caαξ来确定。为保证折纸管单胞的可周期阵列性,设计折痕时,令EF=2DEOF=4DE,腹板DEHG和腹板FEHI的夹角ξ=60°。将2个沿折痕折叠成的单元对称连接,得到弧形折纸薄壁管件,如图1(c)所示。长、宽、高分别用lwh表示(图1中弧形实折线(CFIM)为向外的峰线,弧形虚折线为向内的谷线),其中:

      $ \begin{split} l =& 2\left( {2c\sin \xi + a\sin \alpha \cos \beta } \right),\\w =& OF = 4c,\\h =& a\sin \alpha + 2a\sin \theta \end{split}$

      本文中所研究的模型均使用3D打印机打印制造,在建立3D打印模型时,连接面DEHG与其他竖直腹板的厚度一致。

    • 本研究中所有样品均采用Volunex Mars Plus+ 3D打印机打印。使用的材料为PolyMaxTM PLA,成型后材料密度为1 200 kg/m3,泊松比为0.35,弹性模量为1.97 GPa。利用SHIMADZU AG-25TB/SFL-5AG材料试验机,参照拉伸测试规范ASTM D638-10测试了拉伸试样在准静态(应变率为10−4 s−1)拉伸时的力学参数,初始屈服强度为40 MPa,断裂应变约为0.3。拉伸试样的几何参数及真应力应变曲线如图2(a)所示。

      图  2  利用PolyMaxTM PLA 打印的哑铃形试件的拉伸真实应力应变关系及打印工艺

      Figure 2.  True stress-strain relation of a dumbbell-shaped sample printed with PolyMaxTM PLA and the corresponding printing process

      3D打印试件力学性能与其打印工艺关系密切,对于熔融堆积法打印工艺,试件拉伸力学响应与打印角度相关性较大。本文研究的管状结构打印模式为沿着结构轴线向上铺层打印,为保证基体材料力学性能测试的可靠性,试件打印选取沿长度方向铺层打印,如图2(b)所示。

      根据前文所述几何模型,利用Volunex Mars Plus+3D打印机制造表观尺寸(l×w×h)为80 mm×80 mm×60 mm、预折角α=60°、壁厚τ=1 mm、在面内方向阵列个数n为2×2的薄壁管状结构,在ABAQUS中建立了相应的有限元模型。实验中,采用SHIMADZU AG-25TB/SFL-5AG压缩试验验机,压缩速度为1 mm/min;压头为钢制平板,与结构自由接触。结构在准静态压缩载荷下的实验和数值模拟对比结果如图3所示,两者吻合较好。数值模拟的模型尺寸与实验相同,有关单元类型、材料模型及接触关系的详细信息见第4节。由于数值模拟中未考虑摩擦效应的影响,结构的压缩力略低于实验值。

      图  3  模拟与实验压缩载荷-位移曲线对比

      Figure 3.  Comparison of compressive load-displacement curves between simulation and experimentl

      图3可以看出,结构变形可定性地分为4个阶段。初始压溃阶段:当结构受到初始压缩载荷时,压缩力迅速上升至峰值载荷,约为7.5 kN。之后,结构进入预折角塑性旋转阶段:预折角开始发生塑性旋转变形,同时伴随有局部破坏,此时结构压缩力呈非线性下降趋势。随着预折角的进一步旋转,结构压缩力逐渐上升,当腹板两侧倾斜部分旋转至水平位置时,压缩力主要有腹板承担,在力-位移曲线上出现第2个峰值载荷,约为10 kN,这个过程为腹板塑性屈曲阶段。当腹板受力发生屈曲后,压缩力下降,随后整个结构逐步进入压缩密实化阶段。

    • 针对不同几何参数的弧形折纸模式薄壁结构开展准静态和动态压缩数值模拟。弧形折纸薄壁管单胞根据其预折角α的不同取3种,分别记为A50°、A60°、A70°,与之对比的方管记为SQU。由于常用的多孔夹芯抗冲击防护结构芯层厚度在15~40 mm之间[12-13],因此本节中所涉及的弧形折纸模式单胞结构(实验组)和方管单胞模型(对照组)的相关参数如表1所示。

      模型单胞质量/gh/mmw/mml/mmτ/mmc/mma/mmh0/mmξ /(°)α/(°)θ/(°)β/(°)
      A50°8.65304045.231.101014.2610.9260504261
      A60°8.66304041.571.201011.9810.3760605571
      A70°8.59304038.861.251010.7910.1460706778
      SQU8.66304041.891.47

      表 1  模型参数

      Table 1.  The parameters of the models

    • 结构单胞有限元模型如图4所示,折纸管采用S4R壳单元,上下板为刚体。对于2×2阵列模型,在建模时直接把所有的单胞合并为一个整体。经过网格敏感性和有限元算法验证后(图56),数值模拟中的模型网格尺寸确定为1 mm,上下板和折纸管采用自接触,忽略摩擦的影响,上板采用1 mm/ms的压缩速度来模拟准静态加载过程(总动能必须比压缩过程的总内能小得多,压缩力位移响应必须与所施加的速度无关[14-15])。

      图  4  弧形折痕薄壁管状结构单胞有限元模型

      Figure 4.  The finite element model of a curved origami thin-walled tube

      图  5  网格敏感性验证

      Figure 5.  Validation of mesh sensitivity

      图  6  有限元算法验证

      Figure 6.  Verification of the finite element algorithm

      图5可以看出:折纸管在①处产生初始压溃应力峰值,随后进入预折角塑性旋转阶段(①~②),该阶段称为第1平台段;预折角塑性旋转阶段结束后,应力值开始上升,此时主要承载部分为折纸管的竖直腹板,在折纸管出现第2个压溃应力峰值(③)之后,折纸管进入腹板塑性屈曲阶段(③~④)该阶段称为第2平台段;当竖直腹板完全压溃之后(④之后),整个模型进入密实化阶段。

    • 本文中所涉及的应力应变曲线均为名义应力-名义应变曲线。名义应力σ为压板与薄壁管件接触面上的力与模型初始横截面面积之比,名义应变ε定义为模型沿压缩方向的整体压缩与原长度之比,即:

      ${\rm{\sigma }} = \frac{f}{{lw}},\quad\quad\quad \varepsilon = \frac{{\Delta l}}{h}$

      式中:f为压缩载荷,Δl为压缩位移。

      弧形折纸薄壁管和方管的单胞在准静态压缩作用下的应力-应变曲线如图7所示。由图7可知,随着预折角α的增大,折纸管在初始压溃应力峰值越大,折纸管的第一平台段应力应变曲线的变化趋势一致且应力值的大小基本保持恒定,在ε=0.62时,模型的主要承载部分由预折部分变为竖直的方管,折纸管产生第2个应力峰值,进入腹板塑性屈曲阶段。由于竖直管部分发生坍塌失效后,模型失去承载能力,折纸管的应力应变曲线会出现下降趋势。在此过程中,A50°单胞模型发生整体坍塌,模型承载能力急剧下降;A60°、A70°单胞模型产生渐进压溃模式,在压缩过程中,每个竖直板发生坍塌失效都会对应产生一个峰值应力,因此,A60°、A70°单胞模型的应力应变曲线会发生波动。方管的单胞模型在模拟中的初始压溃应力峰值最大。与只具有一个平台阶段的方管相比较,折纸管的第1平台段应力值与方管平台阶段的应力值相近,折纸管的第2平台段应力值高于方管平台阶段的应力值。

      图  7  不同单胞模型的应力-应变曲线

      Figure 7.  Stress-strain curves of different single-cell models

      弧形折纸薄壁管和方管面内方向阵列个数n为2×2模型的应力-应变曲线如图8所示(方管阵列采用密排结构,相邻单胞方管间紧密接触,接触位置的壁厚为2倍的单壁厚(τ),在有限元模型中设置该处壳单元壁厚为2τ),腹板塑性屈曲阶段,A50°、A60°阵列模型在出现第2个峰值应力之后发生整体坍塌,应力急剧下降,A70°阵列模型竖直管发生整体失稳,第2平台段应力值变小,由于方管阵列模型中间连接部分腹板壁厚为侧边壁厚的两倍,在压缩过程中,当中间腹板发生整体坍塌时,应力有明显增高,屈曲失稳后失去承载能力,应力迅速下降。所有的弧形折纸薄壁管的第2平台段应力值均大于方管。图9为A60°单胞、面内方向阵列个数n为2×2、4×4的阵列模型在准静态作用下的应力-应变曲线,由图可知,2×2、4×4阵列模型的应力应变曲线基本相同,增加面内方向阵列个数对折纸管的应力应变关系没有明显的影响。由于阵列后胞元之间存在一定的横向约束,多胞阵列模型的第一平台段应力明显高于单胞模型,且随着单胞阵列个数的增加,折纸管的第2平台段应力值明显提高。

      图  8  面内方向阵列个数为2×2的阵列模型的应力-应变曲线

      Figure 8.  Stress-strain curves corresponding to the models with the number n of in-plane arrays equal to 2×2

      图  9  面内方向阵列个数不同的A60°模型的应力-应变曲线

      Figure 9.  Stress-strain curves of different in-plane arrays number of A60° models

    • 图10为折纸管面内方向阵列个数n=2×2的阵列模型在冲击速度为10、20 mm/ms载荷作用下的应力-应变曲线。从图10可以看出,折纸管在冲击速度为10 mm/ms载荷作用下的初始压溃应力峰值均低于方管的,预折角的变化对折纸管的第1平台段应力-应变曲线变化趋势及其应力值影响较小。在腹板塑性屈曲阶段,所有折纸管在出现第2个峰值应力之后发生整体坍塌,应力急剧下降,但第2平台段应力值均大于方管的。当冲击速度为20 mm/ms时,A70°阵列模型在压缩过程中表现出逐层压溃模式,在ε=0.42时,折纸管上部预折部分首先达到密实化,产生第2个应力峰值,在ε=0.65时,折纸管预折部分全部密实化,产生第3个应力峰值,折纸管的第1、2平台段应力值均小于方管的。

      图  10  不同冲击速度下,不同阵列模型的应力-应变曲线

      Figure 10.  Stress-strain curves of different models under different impact velocities

    • 图1112为方管和弧形折纸薄壁管的单胞、n=2×2阵列模型在准静态、冲击速度为10、20 mm/ms的比吸能对比。从图1112可以看出:在准静态压缩作用下,对于单胞模型,A70°在压缩位移约22 mm之后,比吸能高于方管的,面内方向阵列个数n=2×2的阵列模型中,方管比吸能优于折纸管;在冲击速度为10、20 mm/ms的载荷作用下,n=2×2的阵列模型中,弧形折纸薄壁管中A70°模型具有较高的比吸能,但方管的比吸能均优于弧形折纸薄壁管。

      图  11  准静态压缩下不同的单胞薄壁管比吸能随位移的变化

      Figure 11.  Specific energy absorption varying with displacement for different single-cell thin-walled tubes under quasi-static compression

      图  12  不同冲击速度下面内方向阵列个数相同的不同薄壁管结构比吸能的对比

      Figure 12.  Comparison of specific energy absorption for different thin-walled tubes with the same number of in-plane arrays under different impact velocities

      对准静态,冲击速度为10、20 mm/ms的单胞、面内方向阵列个数n为2×2模型的数值模拟采用压缩力效率η和比总体效率δ作为薄壁管吸能特性的两种评价指标,压缩效率越大,模型的吸能效果越好。压缩力效率为平均载荷与初始压溃载荷峰值的比值,平均载荷是结构吸收的总能量与最后压缩量之比,比总体效率是单位长度下管的比吸能与初始压溃载荷峰值的比值,所有模型的初始压溃载荷峰值如表2所示。

      冲击速度/(mm·ms−1)n载荷/kN
      A50°A60°A70°SQU
      11×1 0.23 0.65 1.72 6.59
      2×2 1.72 3.91 8.4732.46
      102×2 4.92 9.0613.2242.84
      202×210.4719.7723.7246.14

      表 2  不同模型的初始压溃载荷峰值

      Table 2.  Initial peak crushing loads of different models

      图1314为弧形折纸薄壁管和方管的压缩力效率-压缩位移曲线,从图中可以看出,在准静态压缩作用下,折纸管单胞及面内方向阵列个数n=2×2阵列模型的压缩力效率优于方管。在冲击速度为10 mm/ms的载荷作用下,A50°阵列模型具有较高的压缩力效率,且所有的弧形折纸薄壁管的压缩力效率均高于方管,如图13所示。在冲击速度为20 mm/ms的载荷作用下,当压缩位移d<19.73 mm时,方管阵列模型的压缩力效率高于折纸管,当压缩位移d>19.73 mm时,A50°折纸管阵列模型的压缩力效率高于方管,如图14所示。

      图  13  不同薄壁管模型压缩力效率的对比

      Figure 13.  Comparison of crush force efficiencies for different thin-walled tube models

      图  14  面内方向阵列个数相同的不同薄壁管模型在相同冲击载荷下的压缩力效率随压缩位移的变化

      Figure 14.  Crush force efficiency varying with compressive displacement for different thin-walled tube models with the same number of in-plane arrays under the same impact velocity

      图1516为弧形折纸薄壁管和方管的比总体效率的比较。从图15可以看出,在准静态压缩及冲击速度为10 mm/ms的载荷作用下,折纸管单胞及面内方向阵列个数n=2×2阵列模型的比总体效率均优于方管。从图16可以看出,在冲击速度为20 mm/ms的载荷作用下,当压缩位移d<19.67 mm时,方管的阵列模型比总体效率高于折纸管,当压缩位移d>19.67 mm时,A50°折纸管阵列模型的比总体效率高于方管。

      图  15  不同薄壁管模型比总体效率的对比

      Figure 15.  Comparison of specific total efficiencies for different thin-walled tube models

      图  16  比总体效率随压缩位移的变化

      Figure 16.  Specific total efficiency varying with compressive displacement

    • 利用3D打印技术制备了弧形折纸薄壁管件,基于准静态轴向压缩实验,并利用ABAQUS有限元软件对弧形折纸薄壁管件在准静态,冲击速度10、20 mm/ms下的压缩行为进行了数值模拟,探讨了结构的变形模式,分析了不同预折角以及单胞的阵列个数对其变形过程及承载能力的影响,并与高度相等、质量近似相等的方管作对比,得出以下结论。

      (1)薄壁管件的变形过程可分为4个阶段,分别为初始压溃阶段、预折角塑性旋转阶段、腹板塑性屈曲阶段和完全压溃密实化阶段。

      (2)设计的弧形折痕能够有效降低初始压溃阶段的载荷峰值,减小压缩载荷-位移曲线在预折角塑性旋转阶段、腹板塑性屈曲阶段的振荡幅值,使结构压缩过程可控、按照预定的弧形折痕发展。

      (3)在准静态压缩下:对于单胞模型,A70°模型的比吸能优于方管的;对于多胞管件阵列模型,方管的比吸能均优于折纸管的;折纸管的压缩力效率和比总体效率均优于方管的,且A50°模型的压缩力效率和比总体效率最高。

      (4)在本文涉及的动态冲击载荷下,分别通过比吸能、压缩力效率和比总体效率3个评价指标对比了不同结构的吸能特性:阵列方管的比吸能均优于阵列折纸管的;当冲击速度为10 mm/ms时,折纸管的压缩力效率和比总体效率均优于方管的,其中A50°模型的压缩力效率和比总体效率最高;当冲击速度为20 mm/ms时,仅有A50°模型的压缩力效率和比总体效率优于方管的。

参考文献 (15)

目录

    /

    返回文章
    返回