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

厚度幂指数分布管状结构耐撞性设计准则与方法研究

徐峰祥 张锁 武昆迎

引用本文:
Citation:

厚度幂指数分布管状结构耐撞性设计准则与方法研究

    作者简介: 徐峰祥(1985- ),男,博士,副教授,xufx@whut.edu.cn;
  • 中图分类号: O347; TU311; P315.9

Study on crashworthiness design criteria and method of tubular structures with power exponent distribution of thickness

  • CLC number: O347; TU311; P315.9

  • 摘要: 厚度或质量连续分布技术对车身薄壁结构的轻量化和性能设计有着非常重要,甚至起到决定性的作用,从设计方法上研究连续变厚度结构在车身零部件中的耐撞性应用是安全性设计所需的主要工作。本文研究一种较新颖的薄壁吸能结构,其管壁厚度按照幂指数形式连续分布,根据此分布特点推导出了该薄壁结构在等质量条件下与其他管状结构(比如均匀管、拼焊管和锥管等)之间相关参数的定量解析关系,给出了前者的耐撞性设计准则,评估了不同梯度对幂指数管耐撞性能的影响。分析结果显示,该新颖管状结构比其他截面管具有更理想的耐撞特性。然后,在2个设计区间内对梯度指数分别采样并构造近似模型,采用遗传算法作为求解器得出了非劣解前沿,研究发现高阶响应面近似模型得到的设计结果不一定是最优的。
  • 图 1  厚度连续递增的薄壁结构轴向碰撞示意图

    Figure 1.  Schematic of thin-walled structures with thickness increasing under axial crashing

    图 2  厚度比沿着长度方向的变化情况

    Figure 2.  Thickness changes of thin-walled structures along the length direction

    图 3  幂指数管状结构示意图

    Figure 3.  Schematic showing thickness grading patterns in the axial direction

    图 4  材料应力应变曲线

    Figure 4.  Stress-strain curve of material

    图 5  幂指数分布管(FGT)的试验测试和数值结果比较

    Figure 5.  Comparisons of experimental and numerical results for tubes with power exponent (FGT)

    图 6  变厚度管和锥管的示意图

    Figure 6.  Schematics of variable thickness and tapered tubes

    图 7  三种管结构的有限元模型

    Figure 7.  Finite element model of three tubular structures

    图 8  三种管状结构的吸能特性对比

    Figure 8.  Comparisons of energy absorption capacity of three kinds of tubes

    图 9  幂指数分布管和锥管在不同情况下的碰撞性能对比

    Figure 9.  Crashing performance of tubular structures with power exponent and tapered angles

    图 10  不同响应面模型(RSM)的幂指数分布管(FGT)和均匀管(UT)的Pareto前沿解

    Figure 10.  Pareto solutions of tubular structures with power exponent (FGT) and uniform thickness (UT) by different response surface models (RSMs)

    表 1  三种管状结构在相同质量下的厚度分布

    Table 1.  Thickness distributions of three tubes with the same weight

    ntU/mmt1/mmt2/mmntU/mmt1/mmt2/mm
    0.22.0411.7592.11741.0940.8151.356
    0.41.9431.5172.04161.0170.8021.218
    0.61.6681.3451.97280.9740.8001.137
    0.81.5751.2191.908100.9470.8001.085
    21.2740.9051.621
    下载: 导出CSV

    表 2  比吸能响应面近似模型的精度评估(0≤n≤1)

    Table 2.  Accuracy estimate of response surface model of SAE (0≤n≤1)

    目标函数阶数R2δavgδmax
    SAEG10.952 960.151 600.750 46
    20.953 350.148 300.752 92
    30.953 940.149 120.745 72
    40.954 960.146 500.740 61
    SAEU10.993 740.065 280.151 81
    20.994 780.067 320.116 77
    30.995 120.064 120.118 13
    40.995 130.063 760.120 22
    下载: 导出CSV

    表 3  比吸能响应面近似模型的精度评估(1<n≤10)

    Table 3.  Accuracy estimate of response surface model of SAE (1<n≤10)

    目标函数阶数R2δavgδmax
    SAEG10.987 450.072 090.295 03
    20.991 130.062 970.252 97
    30.993 730.054 780.186 99
    40.995 250.059 800.190 98
    SAEU10.924 910.166 160.709 09
    20.925 660.167 220.724 94
    30.932 140.200 660.658 86
    40.933 110.163 120.667 10
    下载: 导出CSV
  • [1] 杜继涛, 甘屹, 齐从谦. TRB及其轧制应用关键技术 [J]. 汽车技术, 2007(7): 45–48 doi: 10.3969/j.issn.1000-3703.2007.07.013
    DU Jitao, GAN Yi, QI Congqian. Tailor rolled blanks and keys technologies in its rolling applications [J]. Automobile Technology, 2007(7): 45–48 doi: 10.3969/j.issn.1000-3703.2007.07.013
    [2] 何安瑞, 邵健, 孙文权, 等. 冷轧无取向硅钢横向厚差控制 [J]. 机械工程学报, 2011, 47(10): 25–30
    HE Anrui, SHAO Jian, SUN Wenquan, et al. Transverse thickness deviation control of non-oriented silicon steel during cold rolling [J]. Journal of Mechanical Engineering, 2011, 47(10): 25–30
    [3] 赵春璋, 余海东, 王皓, 等. 基于绝对节点坐标法的变截面梁动力学建模与运动变形分析 [J]. 机械工程学报, 2014, 50(17): 38–45
    ZHAO Chunzhang, YU Haidong, WANG Hao, et al. Dynamic modeling and kinematic behavior of variable cross-section beam based on the absolute nodal coordinate formulation [J]. Journal of Mechanical Engineering, 2014, 50(17): 38–45
    [4] 杜继涛, 齐从谦. 变截面薄板在车身制造中的应用研究 [J]. 锻压技术, 2005, 30(S1): 39–43 doi: 10.3969/j.issn.1000-3940.2005.z1.011
    DU Jitao, QI Congqian. Research on applications of variable section blanks in the auto-body manufacturing [J]. Forging & Tamping Technology, 2005, 30(S1): 39–43 doi: 10.3969/j.issn.1000-3940.2005.z1.011
    [5] 徐春国, 任广升, 袁相勇, 等. 轴类件的柔性辊轧成形工艺与实验研究 [J]. 塑性工程学报, 2008, 15(2): 100–104
    XU Chunguo, REN Guangsheng, YUAN Xiangyong, et al. The research on the process of flexible rolling for shaft forgings and experiment [J]. Journal of Plasticity Engineering, 2008, 15(2): 100–104
    [6] 袁相勇, 徐春国, 任广升. 轴向辊轧轧件与轧轮接触面理论分析 [J]. 锻压技术, 2005, 30(S1): 156–159 doi: 10.3969/j.issn.1000-3940.2005.z1.046
    YUAN Xiangyong, XU Chunguo, REN Guangsheng. Theoretic analysis of contact surfaces between workpiece and roller of axial feed bar rolling [J]. Forging & Stamping Technology,, 2005, 30(S1): 156–159 doi: 10.3969/j.issn.1000-3940.2005.z1.046
    [7] 袁海伦, 王泽武, 曾青, 等. 异形截面环件虚拟轧制及其工艺优化 [J]. 塑性工程学报, 2006, 13(6): 15–18 doi: 10.3969/j.issn.1007-2012.2006.06.004
    YUAN Hailun, WANG Zewu, ZENG Qing, et al. Virtual ring rolling and the process optimization of profile ring [J]. Journal of Plasticity Engineering, 2006, 13(6): 15–18 doi: 10.3969/j.issn.1007-2012.2006.06.004
    [8] 王宏雁, 君毅. 汽车车身轻量化结构与轻质材料 [M]. 北京: 北京大学出版社, 2009.
    WANG Hongyan, JUN Yi. Lightweight structures and mateiral of automotive body [M]. Beijing: Beijing Univeristy Press, 2009.
    [9] KOPP R, WIEDNER C, MEYER A. Flexibly rolled sheet metal and its use in sheet metal forming [J]. Sheet Metal, 2005, 6-8: 81–92.
    [10] KOPP R, WIEDNER C, MEYER A. Flexible rolling for load-adapted blanks [J]. International Sheet Metal Review, 2005, 4: 20–24.
    [11] 包向军. 变截面薄板弯曲成形回弹的试验研究和数值模拟: [D]. 上海: 上海交通大学, 2003.
    BAO Xiangjun. Experimental and numerical study on springback of tailor rolled blanks in rotary bending [D]. Shanghai: Shanghai jiaotong University, 2003.
    [12] 任灏宇, 杜继涛. 连续变截面板的轧制控制 [J]. 机械制造与自动化, 2008, 37(2): 143–144 doi: 10.3969/j.issn.1671-5276.2008.02.050
    REN Haoyu, DU Jitao. Roll control of tailor rolled Blanks [J]. Machine Building & Automation, 2008, 37(2): 143–144 doi: 10.3969/j.issn.1671-5276.2008.02.050
    [13] 姜银方, 方雷, 李志飞, 等. 连续变截面板及其应用中存在的关键问题 [J]. 制造技术与机床, 2011(1): 144–148 doi: 10.3969/j.issn.1005-2402.2011.01.043
    JIANG Yinfang, FANG Lei, LI Zhifei, et al. The key technologies in tailor rolled blanks and its application [J]. Manufacturing Technology & Machine Tool, 2011(1): 144–148 doi: 10.3969/j.issn.1005-2402.2011.01.043
    [14] CHUANG C, YANG R, LI G, et al. Multidisciplinary design optimization on vehicle tailor rolled blank design [J]. Structural and Multidisciplinary Optimization, 2008, 35(6): 551–560. doi: 10.1007/s00158-007-0152-0
    [15] YANG R J, FU Y, LI G. Application of tailor rolled blank in vehicle front end for frontal impact [J]. SAE Technical Paper, 2007
    [16] 兰凤崇, 李佳光, 马芳武, 等. 连续变截面板(TRB 板)在汽车前纵梁中的应用及优化分析 [J]. 机械设计与制造, 2014(1): 25–28, 32 doi: 10.3969/j.issn.1001-3997.2014.01.008
    LAN Fengchong, LI Jiaguang, MA Fangwu, et al. The application and optimization analysis of tailor rolled board on automotive front rail structure [J]. Machinery Design & Manufacture, 2014(1): 25–28, 32 doi: 10.3969/j.issn.1001-3997.2014.01.008
    [17] ZHANG X, ZHANG H. Optimal design of functionally graded foam material under impact loading [J]. International Journal of Mechanical Sciences, 2013, 68: 199–211. doi: 10.1016/j.ijmecsci.2013.01.016
    [18] YIN H, WEN G, HOU S, et al. Crushing analysis and multiobjective crashworthiness optimization of honeycomb-filled single and bitubular polygonal tubes [J]. Materials & Design, 2011, 32(8−9): 4449–4460.
  • [1] 李松晏郑志军虞吉林 . 高速列车吸能结构设计和耐撞性分析. 爆炸与冲击, 2015, 35(2): 164-170. doi: 10.11883/1001-1455(2015)02-0164-07
    [2] 常白雪郑志军赵凯何思渊虞吉林 . 具有恒定冲击载荷的梯度泡沫金属材料设计. 爆炸与冲击, 2019, 39(4): 041101-1-041101-9. doi: 10.11883/bzycj-2018-0431
    [3] 殷之平李玉龙黄其青 . 含诱导缺陷薄壁圆管耐撞性优化设计. 爆炸与冲击, 2011, 31(4): 418-422. doi: 10.11883/1001-1455(2011)04-0418-05
    [4] 丁兆洋郑志军虞吉林 . 列车分布式吸能系统的波传播特性和参数分析. 爆炸与冲击, 2019, 39(3): 035101-1-035101-10. doi: 10.11883/bzycj-2018-0053
    [5] 谭丽辉徐涛崔晓梅张炜赵世佳 . 带有圆弧形凹槽金属薄壁圆管抗撞性优化设计. 爆炸与冲击, 2014, 34(5): 547-553. doi: 10.11883/1001-1455(2014)05-0547-07
    [6] 史安顺张先锋乔良魏胜张将何源何勇 . 多功能含能结构材料冲击压缩特性的理论计算. 爆炸与冲击, 2013, 33(2): 148-155. doi: 10.11883/1001-1455(2013)02-0148-08
    [7] 李磊马宏昊沈兆武 . 基于正交设计方法的双锥罩结构优化设计. 爆炸与冲击, 2013, 33(6): 567-573. doi: 10.11883/1001-1455(2013)06-0567-07
    [8] 蔡煜东 . 岩体可爆性等级判别的遗传程序设计方法. 爆炸与冲击, 1995, 15(4): 329-334.
    [9] 任冀宾王斌王振刘军索涛李玉龙 . 某型飞机机翼前缘抗鸟撞结构设计与试验验证. 爆炸与冲击, 2019, 39(2): 025101-1-025101-9. doi: 10.11883/bzycj-2017-0407
    [10] 程帅张德志刘文祥殷文骏师莹菊陈博李焰 . 强脉冲载荷作用下端盖法兰结构的螺栓预紧力设计方法. 爆炸与冲击, 2019, 39(2): 024101-1-024101-7. doi: 10.11883/bzycj-2018-0007
    [11] 王代华刘殿书杜玉兰柳慧鹏 . 含泡沫吸能层防护结构爆炸能量分布的数值模拟研究. 爆炸与冲击, 2006, 26(6): 562-567. doi: 10.11883/1001-1455(2006)06-0562-06
    [12] 王丙斌王蕊 . 空心率对中空夹层钢管混凝土组合柱耐撞性能影响. 爆炸与冲击, 2018, 38(1): 204-211. doi: 10.11883/bzycj-2016-0143
    [13] 蒋宝权李玉龙刘元镛虞青俊 . SiC增强颗粒分布规律对梯度装甲板抗侵彻过程的影响. 爆炸与冲击, 2005, 25(6): 493-498. doi: 10.11883/1001-1455(2005)06-0493-06
    [14] 曾斐潘艺胡时胜 . 泡沫铝缓冲吸能评估及其特性. 爆炸与冲击, 2002, 22(4): 358-362.
    [15] 何源何勇王传婷潘绪超焦俊杰郭磊杨相礼李铨 . 含电子相影响的多功能含能结构材料冲击压缩理论计算. 爆炸与冲击, 2018, 38(1): 217-223. doi: 10.11883/bzycj-2017-0034
    [16] 金骥良王中黔 . 建筑物拆除爆破设计原理和方法. 爆炸与冲击, 1988, 8(2): 139-145.
    [17] 张继春肖清华夏真荣 . 隧道爆破设计智能系统的组成与结构研究. 爆炸与冲击, 2007, 27(5): 455-460. doi: 10.11883/1001-1455(2007)05-0455-06
    [18] 董海平蔡瑞娇穆慧娜曹建华 . 火工品可靠性的感度参数设计法. 爆炸与冲击, 2009, 29(6): 613-616. doi: 10.11883/1001-1455(2009)06-0613-04
    [19] 王丹于海江 . 试验设计在爆炸网络可靠性研究中的应用. 爆炸与冲击, 2015, 35(2): 184-190. doi: 10.11883/1001-1455(2015)02-0184-07
    [20] 唐德高周布奎周早生 . 侧向爆炸冲击波加载作用下钢管吸能特性的实验研究. 爆炸与冲击, 2002, 22(2): 179-183.
  • 加载中
图(10)表(3)
计量
  • 文章访问数:  80
  • HTML全文浏览量:  237
  • PDF下载量:  6
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-01-08
  • 录用日期:  2018-01-24
  • 网络出版日期:  2019-03-25
  • 刊出日期:  2019-03-01

厚度幂指数分布管状结构耐撞性设计准则与方法研究

    作者简介:徐峰祥(1985- ),男,博士,副教授,xufx@whut.edu.cn
  • 1. 武汉理工大学现代汽车零部件技术湖北省重点实验室,湖北 武汉 430070
  • 2. 武汉理工大学汽车零部件技术湖北省协同创新中心,湖北 武汉 430070

摘要: 厚度或质量连续分布技术对车身薄壁结构的轻量化和性能设计有着非常重要,甚至起到决定性的作用,从设计方法上研究连续变厚度结构在车身零部件中的耐撞性应用是安全性设计所需的主要工作。本文研究一种较新颖的薄壁吸能结构,其管壁厚度按照幂指数形式连续分布,根据此分布特点推导出了该薄壁结构在等质量条件下与其他管状结构(比如均匀管、拼焊管和锥管等)之间相关参数的定量解析关系,给出了前者的耐撞性设计准则,评估了不同梯度对幂指数管耐撞性能的影响。分析结果显示,该新颖管状结构比其他截面管具有更理想的耐撞特性。然后,在2个设计区间内对梯度指数分别采样并构造近似模型,采用遗传算法作为求解器得出了非劣解前沿,研究发现高阶响应面近似模型得到的设计结果不一定是最优的。

English Abstract

  • 连续变截面或者变厚度板技术是一种新兴的薄板轧制工艺[1-3],在满足车身结构性能的基本前提下通过该技术可以使车身零部件进一步减轻重量。由于变厚度技术减少了复杂的焊接工序,所以具有一定成本优势,为进一步推广应用提供了良好的基础,研究变厚度板技术对车身轻量化具有较大作用[4]

    与拼焊板结构相比,变厚度板结构的特点是厚度属性可以灵活变化,再加上厚度的连续特性,使其不会出现强度或者应力突变等缺陷;其次,变厚度板技术也不会遇到拼焊板所固有的焊接难题。变厚度板技术的基本原理是通过辊轧轮的上、下运动轧制而形成具有一定厚度分布的板材[5-7],这样采用较少重量的材料就可以制成车身结构件,促使其在汽车车身零部件结构中具有一定的应用前景[8]

    Kopp等[9-10]最早从成形角度研究变厚度板轧制控制以及回弹等问题中发现用较少的材料就可以加工成比较高刚度和强度的车身零部件。包向军[11]从试验角度和数值模拟方面研究了变厚度板的弯曲和回弹;任灏宇等[12]在控制模型及理论基础等方面对连续变厚度板进行了研究;姜银方等[13]通过有限元模拟研究了变厚度板塑性变形过程。

    但上述研究重点基本都是围绕其成形工艺特性,仅有少数学者在轻量化和耐撞性等方面进行了研究。从Yang等[14-15]的研究可以得出结论,变厚度结构具有较高的材料利用率和较大的设计灵活性;兰凤崇等[16]引入充分利用材料特性的连续变厚度板,实现前纵梁的耐撞性和轻量化双重目标,设计结果表明,变厚度结构的变形量和加速度峰值较小,质量也有一定的降低,降幅达到5.21%,并且并没有明显削弱耐撞特性。然而,少量的研究不足以对变厚度结构进行合理的设计。

    本文中,研究厚度连续分布的管状结构的耐撞性设计准则及方法,假定厚度的梯度变化形式是幂指数分布,推导等质量条件下其与其它管状结构(比如均匀管、拼焊管以及锥管等)之间的定量解析关系并对比其吸能特性,并对幂指数梯度变化形式的管状结构进行耐撞性优化设计。

    • 理论上本文封闭的管状结构应该是通过板结构的折弯、激光焊接等工艺制造而成,由于焊缝较窄,对试件性能的影响较小,尤其是对动态大变形耐撞性问题,焊缝带来的影响更小,所以在数值模拟中可以将焊缝视为共节点连接方式,即不考虑焊缝的材料和尺寸属性。图1是厚度连续递增的薄壁结构轴向碰撞示意图,上端部壁厚ttop和下端部厚度tbot分别是碰撞开始端和末端,中间过渡厚度由某种特定函数关系给定。

      图  1  厚度连续递增的薄壁结构轴向碰撞示意图

      Figure 1.  Schematic of thin-walled structures with thickness increasing under axial crashing

      从管状结构的轴向压溃过程可知,薄壁结构的碰撞端即开始端首先发生压溃变形,所以管状吸能结构的厚度分布应该是从最薄逐渐过渡到最厚,即碰撞端的厚度应该是最薄的(ttop=tmin),固定端则最厚(tbot=tmax),这样就可以使管状结构更容易发生渐进式(或手风琴式)折叠。与文献[17]中功能梯度泡沫材料的密度分布形式相类似,本文将特定管状结构的厚度假设为幂指数分布形式,也就是厚度沿着轴向逐渐增大的形式(图2),最小和最大厚度分别用t0tm表示,即:

      图  2  厚度比沿着长度方向的变化情况

      Figure 2.  Thickness changes of thin-walled structures along the length direction

      $ t(x) = {t_{\rm{m}}} + ({t_{\rm{m}}} - {t_0}){(x/L)^n} $

      式中:n为梯度指数,x为距离管状结构上端部的长度,L为薄壁管的总长度。

    • 为了方便讨论问题,考虑到目前制造工艺方面的不足,此处通过简单的线切割手段得到厚度线性分布(n=1.0)的圆管结构(图3)。碰撞端(上端)有最小的厚度值tmin,固定端(下端)有最大的厚度值tmax。本文试件的尺寸为:t0=tmin=1.2 mm,tm=tmax=1.5 mm,D=45 mm,L=120 mm。所采用的材料为铝合金(AA6061-T5),弹性模量为68.2 GPa,密度为2 700 kg/m3,泊松比为0.33,材料属性通过准静态拉伸试验(2 mm/min)测得,得到的应力应变曲线如图4所示。

      图  3  幂指数管状结构示意图

      Figure 3.  Schematic showing thickness grading patterns in the axial direction

      图  4  材料应力应变曲线

      Figure 4.  Stress-strain curve of material

      如果将变厚度管分割成若干层,那么均匀管的等效厚度(与幂指数分布管有相同质量)为:

      $ {t_{{\rm{avg}}}} =\left. \sum\limits_{i = 1}^{{N_{\rm{S}}}} {({t_{\rm{i}}}{L_{\rm{e}}})} \right/({N_{\rm{S}}}{L_{\rm{e}}}) = \left.\sum\limits_{i = 1}^{{N_{\rm{S}}}} {{t_{\rm{i}}}} \right/{N_{\rm{S}}} $

      式中:ti为第i层的厚度,NS 为总层数,Le为单层长度。通过不同层数构造的收敛性分析,此处设定31层已满足要求。

      图5所示,数值模拟和试验得到的峰值碰撞力分别为31.5和30.1 kN,相对误差控制在5%以内;另外,从变形模型和褶皱数量来看,结果同样具有良好的吻合程度,这证明了该模型可以进一步预测其他幂指数分布管的变形模式和行为。

      图  5  幂指数分布管(FGT)的试验测试和数值结果比较

      Figure 5.  Comparisons of experimental and numerical results for tubes with power exponent (FGT)

    • 均匀等直管和锥管结构(图6)都是常见的薄壁吸能结构。决定锥管碰撞性能的参数是壁厚和锥角;而决定变厚度管耐撞性的是梯度函数。采用等质量原则进行参数的等效处理,得到不同管状结构几何参数之间的关系,这样可以方便改变结构参数,分析在等质量条件下不同截面形状的碰撞特性。

      图  6  变厚度管和锥管的示意图

      Figure 6.  Schematics of variable thickness and tapered tubes

    • 在相同体积或者质量的情况下,等直管和幂指数分布管的几何参数之间的关系可以写成:

      $ \int_0^h {t(x){\rm d}x} = {t_{\rm{U}}}h $

      式中:h是管状结构的高度,tU是等直管的等效厚度。

      等直管的等效厚度tU和幂指数分布管的指数n能够显示表达为

      $ {t_{\rm{U}}} = {t_0} + ({t_{\rm{m}}} - {t_0})/(n + 1) $

      由式(4)可知,对于特定的厚度分布,等直管具有相对应的等效厚度。

    • 与前面相似,幂指数分布管和锥管相关的几何参数之间的关系可以表达为:

      $ 2a\int_0^h {t(x){\rm d}x} = h{t_{\rm{T}}}(2a + h/\tan \theta ) $

      式中:tTθ分别为锥管的厚度和倾斜角度,a为方管宽度的一半或圆管的半径。

      如果锥管的厚度tT等于幂指数管的最小厚度t0,那么决定性能的梯度指数n和倾斜角度θ之间存在一定的关系:

      $ \frac{{n + 1}}{{\tan \theta }} = 2\left( {\frac{a}{h}} \right)\left( {\frac{{{t_{\rm{m}}}}}{{{t_0}}} - 1} \right) $

      所以存在有特定倾斜角度θ的锥管,即:

      $ \theta = \left\{ \begin{aligned} &{\displaystyle\frac{\text{π}}{2}\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad n = 0\;}\\ &\arctan\left[\displaystyle\frac{n+a}{2}\cdot \frac{h}{a}\cdot \left(\displaystyle\frac{1}{t_{\rm m}/{t_0}-1}\right)\right]\quad\quad 0{\text{<}}\theta {\text{≤}}\text{π} /2,\quad n\text{>}0 \end{aligned} \right. $

      更进一步地,n>0时,有:

      $ \theta {\rm{ = arctan[}}{C_0}{\rm{(}}n{\rm{ + 1)]}} $

      式中:C0=h/[2a(tm/t0-1)]。

      若锥管壁厚tT与幂指数分布管的下端厚度tm 相同,那么将存在以下关系:

      $ \theta {\rm{ = }}\left\{ {\begin{aligned} &{\displaystyle\frac{\text{π}}{2}\;\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\;\; n = 0}\\ &{\text{π} {\rm{ - arctan}}\left| {{C_{\rm{m}}}\left( {\displaystyle\frac{{n{\rm{ + 1}}}}{n}} \right)} \right|\;\;\;\;\;\quad\quad\text{π} /2 {\text{<}} \theta {\text{<}} {\text{π}} ,\quad{\rm{ }}n \text{>} 0} \end{aligned}} \right. $

      式中:C0=h/[2a(t0/tm-1)]。

      注意到上述解析关系对截面是方形和圆形的管状结构都是一样的,为了讨论方便,在此仅对方形截面的管状结构进行研究。

    • 在实际的工程应用中,锥管的倾斜角度不应过大(假设该角度不小于π/3),因此:

      $ {C_0} \text{≥} \sqrt 3 /(n + 1) $

      根据式(10),有:

      $ \frac{{2\sqrt 3 }}{{n + 1}}\left(\displaystyle\frac{{{t_{\rm{m}}}}}{{{t_0}}} - 1\right) \text{≤} \frac{h}{a} $

      即:

      $ n + 1 \text{≥} \frac{{2\sqrt 3 a}}{h}\left(\displaystyle\frac{{{t_{\rm{m}}}}}{{{t_0}}} - 1\right) $

      对于tT=tm的情况,锥管上端尺寸大于下端,即倾斜角度为钝角(θ>π/2),这就要求2a>2h tan (θ-π/2),则有:

      $ \frac{a}{h} \text{>} \frac{n}{{\left| {{C_{\rm{m}}}} \right|(n + 1)}} $

      Cm代入至式(13)可以得到幂指数分布管的几何参数与梯度指数之间的表达式:

      $ \frac{n}{{n + 1}}\left(1 - \frac{{{t_0}}}{{{t_{\rm{m}}}}}\right) \text{<} \frac{1}{2} $

    • 首先对比幂指数分布管(FGT)、拼焊管(TWB)、均匀管(UT)之间的吸能特性,有限元模型如图7所示。幂指数分布管的上、下端厚度分别为0.8 mm和2.2 mm;采用厚度不同、材料相同的两拼组合,其厚度分别为幂指数分布管的上/下半部分的等效厚度,均匀管等效厚度由式(2)给出。不同指数n下等效拼焊管及均匀管的厚度见表1

      图  7  三种管结构的有限元模型

      Figure 7.  Finite element model of three tubular structures

      ntU/mmt1/mmt2/mmntU/mmt1/mmt2/mm
      0.22.0411.7592.11741.0940.8151.356
      0.41.9431.5172.04161.0170.8021.218
      0.61.6681.3451.97280.9740.8001.137
      0.81.5751.2191.908100.9470.8001.085
      21.2740.9051.621

      表 1  三种管状结构在相同质量下的厚度分布

      Table 1.  Thickness distributions of three tubes with the same weight

      图8(a)可以看出,当n<1.0时,拼焊管的峰值载荷位于均匀管和幂指数管之间,幂指数分布管所得到的峰值载荷最小,这说明其有效降低了峰值碰撞力,具有明显的安全优势。而均匀管的峰值碰撞力在不同的梯度指数下均是最大的。另外,幂指数分布管在比吸能指标上也具有明显的优势(图8(b))。综上所述,幂指数管是一种理想的吸能薄壁结构。

      图  8  三种管状结构的吸能特性对比

      Figure 8.  Comparisons of energy absorption capacity of three kinds of tubes

      图9是幂指数管与锥管在不同情况下的碰撞性能对比,A、B分别代表0<n≤1和1<n≤10,1~6代表不同的锥管厚度,从1变化到6,锥管厚度不断增大,角度也逐步倾斜。可以看出幂指数管的碰撞载荷随着梯度值的不断增大而减小。

      图  9  幂指数分布管和锥管在不同情况下的碰撞性能对比

      Figure 9.  Crashing performance of tubular structures with power exponent and tapered angles

      综上所述,幂指数分布管的比吸能和峰值碰撞载荷等碰撞指标优于等直管和锥管,主要的原因在于幂指数管的厚度是按照一定梯度形式变化的,可使幂指数管的耐撞性效果能够进行合理预测和评估,使得压溃变形程度更平缓并可控。

    • 对于不同的梯度指数,幂指数分布管的厚度分布是不同的,所以在此将梯度指数作为设计变量,并将比吸能和峰值碰撞力设定为耐撞性指标,这样幂指数分布管状结构的多目标设计问题就可以表达函数:

      $ f = \max [e(n), - {F_{\max }}(n)] $

      幂指数管状结构的碰撞过程是非线性行为,在数学上不可能获取目标函数的解析函数,所以通过近似模型技术(如响应面等)来拟合耐撞性响应[18]。需在两个设计空间0≤n≤1和1≤n<10内随机选取测试点来构建SEA和Fmax的响应面模型。反求出响应面近似模型的参数之后,接下来需检验模型的精度。常用的误差估计方式有决定系数R2,相对平均绝对误差 δavg,相对最大绝对误差δmax

      $ {R^2} = \left. 1 - \sum\limits_{i = 1}^n {{{({y_{\rm{i}}} - {{\hat y}_{\rm{i}}})}^2}} \right/\sum\limits_{i = 1}^n {{{({y_{\rm{i}}} - \bar y)}^2}} $

      $ {\delta _{{\rm{avg}}}} = \left.\sum\limits_{i = 1}^n {|{y_{\rm{i}}} - {{\hat y}_{\rm{i}}}|} \right/\sum\limits_{i = 1}^n {|{y_{\rm{i}}} - \bar y|} $

      $ {\delta _{\max }} = \max \{ |{y_{\rm{i}}} - {\hat y_1}|, \ldots ,|{y_{\rm{i}}} - {\hat y_{\rm{n}}}|\} \left/\left( {\frac{1}{n}\sum\limits_{i = 1}^n {{y_{\rm{i}}} - \bar y} } \right)\right. $

      式中:yi${{\hat y}_{\rm{i}}}$$\bar y$分别为响应的实测值、预测值和预测值的平均值。表2表3是比吸能响应面近似模型的精度评估结果,其中上标G和U分别表示幂指数分布管和均匀管,可以看出阶数越高,模型精度越高。

      目标函数阶数R2δavgδmax
      SAEG10.952 960.151 600.750 46
      20.953 350.148 300.752 92
      30.953 940.149 120.745 72
      40.954 960.146 500.740 61
      SAEU10.993 740.065 280.151 81
      20.994 780.067 320.116 77
      30.995 120.064 120.118 13
      40.995 130.063 760.120 22

      表 2  比吸能响应面近似模型的精度评估(0≤n≤1)

      Table 2.  Accuracy estimate of response surface model of SAE (0≤n≤1)

      目标函数阶数R2δavgδmax
      SAEG10.987 450.072 090.295 03
      20.991 130.062 970.252 97
      30.993 730.054 780.186 99
      40.995 250.059 800.190 98
      SAEU10.924 910.166 160.709 09
      20.925 660.167 220.724 94
      30.932 140.200 660.658 86
      40.933 110.163 120.667 10

      表 3  比吸能响应面近似模型的精度评估(1<n≤10)

      Table 3.  Accuracy estimate of response surface model of SAE (1<n≤10)

      图10给出了通过响应面模型得到的幂指数分布管和均匀管的优化结果(Pareto前沿),幂指数分布管的Pareto前沿位于右上方,这表明其比均匀管具有更好的优势。通过不同阶数响应面模型得到的均匀管的设计结果基本一致,而对于幂指数分布管的Pareto前沿解来说具有不同的结果。比如,在图10(a)中点A的左边,具有更高的吸能能力(拥有阶数更高的响应面),在点B的右侧,通过二阶多项式近似模型可以得到更好的优化结果,说明较高精度的响应面近似模型并不一定可以得到更优的结果。

      图  10  不同响应面模型(RSM)的幂指数分布管(FGT)和均匀管(UT)的Pareto前沿解

      Figure 10.  Pareto solutions of tubular structures with power exponent (FGT) and uniform thickness (UT) by different response surface models (RSMs)

    • 本文研究了一种厚度幂指数分布管状薄壁吸能结构,分析了不同厚度梯度对其耐撞性能的影响,最后研究了其耐撞性优化设计问题,主要得到以下结论:

      (1) 根据质量等效原则,推导出了三种吸能结构(幂指数管、等直管和锥管等薄壁结构)的尺寸参数之间的关系,为对比分析其耐撞性能提供了良好的数据基础;

      (2) 得到了幂指数分布管的几何设计准则,即梯度指数应满足不等式条件n+1≥2$\sqrt 3 $a/h(tm/t0−1);

      (3) 将梯度指数作为优化设计变量,在两个设计区间内采样并分别构建合理的响应面近似模型,通过对比证实高阶响应面近似模型得到的设计结果不一定是最优的。

参考文献 (18)

目录

    /

    返回文章
    返回