Oblique impact resistance of a bionic thin-walled tube based on antles osteon
-
摘要: 为提高薄壁管结构的耐撞性和吸能性，基于鹿角骨单位结构特征，结合结构仿生学原理设计出内径相同、外径等梯度逐层递减的仿生薄壁管。采用有限元法对75种仿生薄壁管结构进行10°、20°、30°等3种斜向冲击角度的吸能特性模拟；通过多项式回归元模型和多目标粒子群优化算法进行优化，以Pareto前沿最优原则得到各目标最优化的配置方案；采用最小距离选择法进行优化分析，得到各配置方案的最优结构设计参数。结果表明：仅考虑单一冲击角度时，在10°、20°、30°冲击角度下的仿生薄壁管耐撞性最优的仿生层数n均为6，最大壁厚与厚度梯度值参数组合tmax-a分别为2.84 mm-0.38 mm、2.89 mm-0.29 mm、2.91 mm-0.34 mm；综合考虑多种冲击角度权重因数不同配置方案时，仿生薄壁管耐撞性最优的仿生层数n均为6，最大壁厚与厚度梯度值参数组合tmax-a分别为2.95 mm-0.28 mm、2.92 mm-0.30 mm、2.85 mm-0.33 mm。Abstract: Some achievements have been made in the study on mechanical properties of antler, but they have not been applied in engineering practice, especially in the study of thin-walled tubes similar to antler crashworthiness. In order to improve the crashworthiness and energy absorption of the thin-walled tube structures, a bionic thin-walled tube with the same inner diameter and equal gradient of outer diameter was designed based on the structural characteristics of antler bone and the principle of structural bionics. The finite element method was used to simulate the energy absorption characteristics of 75 kinds of bionic thin-walled tube structures under the oblique impacts with the impact angles of 10°, 20° and 30°. The polynomial regression element model and multi-objective particle swarm optimization algorithm were used to optimize, and the Pareto front optimization principle was used to obtain the optimal allocation scheme of each target. The minimum distance selection method was used in optimization analysis to obtain the optimal structural design parameters of each scheme. The optimization method used in this study can provide reference for the follow-up research on the crashworthiness of thin-walled tubes, and the optimal structure of bionic thin-walled tubes can provide reference for practical engineering application. The results show that when only considering a single impact angle, the optimal number of biomimetic layers n is 6, and the parameter combination of maximum wall thickness and thickness gradient tmax-a is 2.84 mm-0.38 mm, 2.89 mm-0.29 mm, 2.91 mm-0.34 mm, respectively under 10°, 20° and 30° impact angles. Considering various impact angle weight factors and different configuration schemes, the optimal number of biomimetic layers n is 6, and the parameter combination of maximum wall thickness and thickness gradient tmax-a is 2.95 mm-0.28 mm, 2.92 mm-0.30 mm and 2.85 mm-0.33 mm, respectively.
-
Key words:
- thin-walled tube /
- structural bionics /
- oblique impact /
- crashworthiness /
- multi-objective optimization
-
众所周知，薄壁管质量轻、比强度高、比刚度大及抗冲击性能好，被广泛用于车辆与运载、航空航天等领域的安全吸能装置中[1]。针对薄壁管的吸能特性已进行了大量研究，研究结果表明，薄壁管在轴向冲击下发生的渐进式压溃变形能有效吸收冲击能量。徐峰祥等[2]研究了一种管壁厚按照幂指数形式连续分布的薄壁吸能结构，给出该管耐撞性设计准则，并根据此分布特点推导出等质量条件下其他管状结构（如均匀管、拼焊管和锥管等）的相关结构参数定量解析关系，得到了幂指数分布管的几何设计准则，评估了不同幂指数梯度对幂指数分布管耐撞性能的影响。郝文乾等[3]在方管的基础上引入折纹结构，利用几何关系建立折纹管的折角公式，并采用LS-DYNA研究了6种折纹管的屈曲模态和吸能性能，结果表明折纹的引入能有效提高方管的吸能性能。谭丽辉等[4]在金属薄壁圆管的基础上，引入圆弧形凹槽诱导结构，并进行多目标优化设计，得出使薄壁管吸能效果最优的凹槽数量和半径，得到了理想的诱导凹槽优化结构。张涛等[5]在圆管和方管管壁上增加诱导缺陷并提出新型组合截面研究高速冲击下薄壁组合结构吸能特性，得到抗冲击吸能性能最优的组合结构。殷之平等[6]运用最小二乘曲面拟合和遗传算法对含诱导缺陷的薄壁圆管进行耐撞性优化设计，结果表明优化后的结构具有良好的吸能效果，便于在结构设计中应用。然而，斜向冲击下的吸能特性会随冲击角度的变化而变化，且在实际应用中发生斜向冲击的情况较多[7]。Sun等[8]基于整数编码遗传算法针对多个斜向冲击角度下的多胞管进行拓扑优化研究，结果表明优化后的结构显著提高多胞管耐撞性能。Alkhatib等[9]通过数值模拟研究了斜向冲击下波纹锥形管的变形模式和吸能特性，结果表明增大冲击角度会导致吸能效果下降，冲击角度在0°～40°时比吸能下降约54%。Asanjarani等[10]研究了一种带凹槽锥形薄壁方管的单目标和多目标优化设计方法，并得到吸能效果最优的结构几何参数。亓昌等[11]研究了长径比、壁厚和锥度等结构参数对轴对称锥形多胞薄壁方管不同冲击角度下的失效模式和吸能特性的影响，并拟合出可用于斜向冲击下比吸能和峰值载荷预测的解析公式，为后续薄壁管斜向冲击耐撞性研究提供了理论基础。
在常规结构的薄壁管设计中，已经取得了不小的成绩，随着工程仿生学的发展[12]，研究人员注意到将生物结构耐撞性特征运用到薄壁管结构中，仿生薄壁管表现出了优异的耐撞性。Huang等[13]以雀尾螳螂虾螯结构为仿生原型，设计了一种仿生单元多胞薄壁管，研究了其轴向和斜向冲击下的耐撞性能，通过多目标粒子群优化算法对结构参数进行优化并得到最优解。许述财等[14]受竹结构启发设计出由仿生节、仿生单元和仿生内管组成的仿生薄壁管，并通过响应面法优化得到仿生薄壁管的最优结构。Zou[15]、Song等[16]基于竹结构特征设计仿生薄壁管，研究发现节可以引导变形，具有仿生节的薄壁管比吸能较普通管更高，通过多目标优化，获得吸能效果最优的仿生薄壁管结构。Li等[17]提出一种仿莲藕填充薄壁管结构，同时将最优拉丁超立方体抽样、响应面法和第二代非支配排序遗传算法集合到一起进行多目标设计优化，所提出的仿生薄壁管结构具有优异的吸能效果。Yin等[18]受马尾草结构特征启发设计泡沫填充仿生薄壁管，提出一种基于集合元模型、第二代非支配排序遗传算法、“3-sigma”稳健设计和蒙特卡罗模拟的多目标鲁棒优化方法，结果表明鲁棒优化结果比确定性优化方法更可靠，优化后的结构具有更优异的耐撞性。
鹿角是由松质骨、皮质层和管壁组成的致密骨管，具有优异机械性能[19]。Picavet等[20]从鹿的性别、活动的地理位置，鹿角形态学、水合状态和矿物组成等几个方面，比较已有研究数据，证实了组织学特性在鹿角韧性中起关键作用，研究发现鹿角因其优异的组织特性和机械特性，可以作为绿色、生态和可持续的生物材料来源。Fang等[21]对鹿角致密骨微观结构进行了研究，用材料试验机分别测试了鹿茸试件在轴向、纵向、横向3个不同方向的弯曲力学性能，然后用扫描电子显微镜观察3个不同方向的裂纹扩展路径和断口形貌，试验结果表明鹿角具有各向异性的力学性能，鹿角致密骨的横向断裂能远大于纵向和径向的断裂能。Currey等[22]对鹿角的湿标本和干标本以及鹿股骨的湿标本进行了力学测试，结果发现：与鹿股骨湿标本相比，鹿角湿标本的弹性模量和抗弯强度较低，但断裂功较高；与鹿股骨湿标本相比，鹿角干标本的杨氏模量略低，但抗弯强度和折断功要高得多，鹿角干标本的抗冲击吸能特性优于鹿股骨湿标本；此外，与普通骨头相比，鹿角也表现出优异的吸能性能。Davison等[23]综述了决定鹿角强度的成分，如矿化度、结晶度、胶原特性和骨细胞活性等材料特性和骨小梁的直径、厚度、皮质壳的孔隙率、骨小梁网络的连通性、各向异性、骨小梁应力水平、微裂纹等骨单位结构特性，但无法量化这些成分对鹿角强度的影响，通过相互关联的性质只能粗略地了解每种成分各自独立发挥的作用。
对鹿角的机械性能已有一定的研究成果，但在工程实际中应用较少。在进行轻量化薄壁管吸能研究中受鹿角骨单位启发，基于结构仿生学，结合鹿角骨单位结构特征及力学特性，利用结构仿生原理设计一种内径相同、外径等梯度逐层递减的仿生薄壁管，对不同层数的仿生薄壁管进行仿真分析，研究其斜向冲击载荷下的吸能特性，采用多项式回归元模型和多目标粒子群优化算法进行优化，获取仿生薄壁管最优结构几何参数。
1. 仿生设计
1.1 仿生原型分析
鹿角是鹿头上形成的骨质突起，外形是一种细长的管状结构，在大量的打斗与冲撞事例中，关于鹿角断裂的观察结果证明了鹿角不寻常的韧性[24]。研究结果发现鹿角由外层密质骨和内芯层松质骨（骨小梁或海绵状）组成，密质骨由围绕中央血管的同心薄壁骨单位组成，是长骨骨干起支撑作用的主要结构单位，在经受高冲击载荷和大弯曲力矩时不会断裂，图1是鹿角的分层结构。
鹿角骨单位与薄壁管在功能和载荷方面具有以下相似性。
（1）功能相似性：薄壁管通过自身屈曲变形吸收冲击产生的能量，从而避免整体结构的破坏，而鹿角的存在同样是承受来自外界施加的载荷，通过骨单位的微变形将载荷进行削弱和化解。
（2）载荷相似性：在常见的车辆碰撞过程中，为了保护整体结构，保险杠前加装薄壁管来抵消来自各个方向的部分载荷。同样，在鹿群之间的争斗中，鹿角承受来自各个方向的冲击。
1.2 仿生结构设计
基于以上分析，对骨单位进行仿生设计，将其应用到薄壁管设计中，以下简称仿生薄壁管。如图2所示，本研究参照骨单位特征，设置仿生薄壁管内径d相同，将层数n、最大壁厚tmax和由下到上逐层等值递减的厚度梯度值a作为设计变量，其中层数n等间距选择3个水平，最大壁厚tmax和厚度梯度值a等间距选择5个水平，探讨各变量对斜向冲击下的仿生薄壁管吸能特性影响规律。
2. 仿真模型
2.1 有限元模型
利用有限元软件HypermMesh进行建模，材料选择为AA6061铝合金，结构选择其中一组参数，层数n设为6层，最大壁厚tmax设为2.80 mm，厚度梯度值a设为0.40 mm，内径d设为44 mm，管高度h设为120 mm，网格采用shell单元。通过对比网格大小对仿真结果的影响，发现网格大小虽影响仿真结果，但网格太小会延长运算时间。四边形网格在1.5 mm和2.0 mm时，仿真结果差距不大，因此选择四边形网格大小为2.0 mm[26]。设置AA6061铝合金对应的MATL24材料属性，密度ρ为2.7×103 kg/m3，杨氏模量E为70 GPa，泊松比μ为0.3，屈服强度为0.25 GPa。刚性墙质量设置为300 kg，冲击速度设置为10 m/s，冲击方向为轴向冲击，考虑到刚性墙与管的接触以及管自身变形产生的接触，分析过程中采用自动单面接触算法（CONTACT_AUTOMATIC_SINGLE_SURFACE），设置面与面之间的摩擦因数为0.3。由于铝合金材料对应变率不敏感，建模时不考虑应变率对材料参数的影响。利用非线性有限元软件LS-DYNA对建立好的有限元模型进行仿真，结果如图3（a）所示。
图 3 仿生薄壁管在10 m/s轴向冲击下的变形模式Figure 3. Deformation patterns of a bionic tube under the axial impact of 10 m/s2.2 有限元模型验证
为验证有限元模型的准确性，采用机械加工的方法制造出样件，材料为铝合金AA6061，结构参数同2.1节中的有限元模型参数，并进行冲击试验，试验中冲击力由安装在锤头上的力传感器测量，管的变形与压缩位移由高速摄像提取。试验室落锤试验台的落锤质量为186 kg，且下落高度为10 m，其最大冲击动能远大于仿真模型所设置的刚性墙的冲击动能，可通过调节落锤下落高度控制其冲击速度以使落锤的冲击动能与仿真模型中刚性墙的冲击动能相同。经过计算得到落锤下落高度为8.23 m，接触时的初始冲击速度为12.70 m/s，可与仿真模型中刚性墙的冲击动能一致。薄壁管样件的变形模式如图3（b）所示，与图3（a）中的有限元模型变形模式基本吻合，均呈现出逐层渐进折叠的变形模式。在186 kg的落锤下落高度为10 m的轴向冲击载荷下仿生薄壁管的力-位移曲线如图4（a）所示，在个别直径突变处表现出较大的差距，原因是相邻直径的突然变化造成应力集中，而使力突然增大随后降低；仿生薄壁管的吸能量-位移曲线如图4（b）所示，仿真值与实验值基本相同。整体来看，有限元建模方法是准确的，可用有限元建模方法来研究仿生管的耐撞性能。
图 4 在186 kg的落锤下落高度为10 m的轴向冲击载荷下仿生薄壁管的力-位移曲线和吸能量-位移曲线Figure 4. Force-displacement and energy absorption-displacement curves of the bionic tube under the axial impact of a drop hammer of 186 kg droping from the heigh of 10 m3. 仿真试验
3.1 试验方案
仿生薄壁管的防撞性能与其本身结构参数息息相关，本研究对层数n、最大壁厚tmax和厚度梯度值a这3种结构参数以及冲击角度α这一加载状态量进行分析。设计一个试验方案，分别将仿生薄壁管最大壁厚tmax与厚度梯度值a分为5个等间距水平，如表1所示，分别定义tmax-a组合，如2.80 mm-0.20 mm、2.80 mm-0.25 mm、······为组别1、2、······，以此类推，共25个组别，考虑到每个组别都有4、5、6等3种层数，因此共有75组不同的结构参数组合，考察其在冲击角度
α =10°，20°，30°时的耐撞性能，共225组试验样本。仿生薄壁管有限元模型基本参数设置同上。表 1 试验组别Table 1. Group of test factors因素 组别 1 2 3 4 5 6 7 8 9 10 11 12 tmax/mm 2.80 2.80 2.80 2.80 2.80 2.85 2.85 2.85 2.85 2.85 2.9 2.90 a/mm 0.20 0.25 0.30 0.35 0.40 0.20 0.25 0.30 0.35 0.40 0.20 0.25 因素 组别 13 14 15 16 17 18 19 20 21 22 23 24 25 tmax/mm 2.90 2.90 2.90 2.95 2.95 2.95 2.95 2.95 3.00 3.00 3.00 3.00 3.00 a/mm 0.30 0.35 0.40 0.20 0.25 0.30 0.35 0.40 0.20 0.25 0.30 0.35 0.40 3.2 评价指标
选择合适的设计准则对于获得有效的优化设计参数是至关重要的，本研究选择学者们一致认可的设计指标[1]：比吸能（specific energy absorption，SEA，e）和峰值载荷（peak crush force，PCF，f）。
比吸能表示单位质量薄壁管吸收的能量：
e=E/m (1) 式中：E为薄壁管吸收的能量，m为薄壁管质量。
峰值载荷表示薄壁管在变形过程中出现的最大冲击载荷。
薄壁管的设计目标为轻量化和高效吸能，需要单位质量吸收最大能量，因此选择e为性能指标，并在优化研究时作为目标函数取最大值。另一方面，冲击过程中薄壁管不应产生较大的峰值载荷，以避免可能产生的更严重伤害或损失，因此选择f为性能指标，并在优化研究时作为目标函数取最小值。
3.3 试验结果
利用LS-DYNA软件对75种具有不同层数n、最大壁厚tmax和厚度梯度值a的仿生薄壁管进行了仿真计算。选择仿生薄壁管变形过程的前100 mm压溃位移进行计算，耐撞性评价指标选择e和f。以下是75种仿生薄壁管结构在冲击角度α = 10°，20°，30°时的性能指标，分别对其不同角度下的吸能特性进行分析，共计225种工况。
在10°冲击下的75组仿生薄壁管的e和f如图5所示。图5（a）中，6层的仿生薄壁管e较高，而4层的仿生薄壁管e较低。其中，层数为6，最大壁厚为2.85 mm，厚度梯度值为0.25 mm的仿生薄壁管e最高，为46.68 kJ/kg。图5（b）中，4层的仿生薄壁管f较低，而6层的仿生薄壁管f较高。其中，层数为4，最大壁厚为2.80 mm，厚度梯度值为0.40的仿生薄壁管f最低。
在20°冲击下的75组仿生薄壁管的e和f如图6所示。图6（a）中，6层的仿生薄壁管e较高，而4层的仿生薄壁管e较低。其中，层数为6，最大壁厚为2.95 mm，厚度梯度值为0.40 mm的仿生薄壁管e最高，为37.42 kJ/kg。图6（b）中，4层的仿生薄壁管f较低，而6层的仿生薄壁管f较高。其中，层数为5，最大壁厚为2.80 mm，厚度梯度值为0.30 mm的仿生薄壁管f最低。
在30°冲击下的75组仿生薄壁管的e和f如图7所示。图7（a）中，6层的仿生薄壁管e较高，而4层的仿生薄壁管e较低。其中，层数为6层，最大壁厚为2.95 mm，厚度梯度值为0.40 mm的仿生薄壁管e最高，为46.68 kJ/kg。图7（b）中，4层的仿生薄壁管f较低，而6层的仿生薄壁管f较高。其中，层数为4，最大壁厚为3.00 mm，厚度梯度值为0.20 mm的仿生薄壁管最低。
通过上述分析发现，不同的层数n、最大壁厚tmax、厚度梯度值a以及冲击角度α对仿生薄壁管的耐撞性有显著影响，但4个参数对仿生薄壁管性能指标的影响程度不同，导致各类仿生薄壁管耐撞性存在差异，仍不能确定具有哪种参数的仿生薄壁管耐撞性最好。因此，为了协调性能指标受参数的影响，定义了性能指标的加权组合，综合考虑4个参数对仿生薄壁管耐撞性的影响。
4. 多目标优化
4.1 优化模型
本文中测试样本点的方法是全因子试验，共生成225个试验样本，同时，采用三次多项式回归元模型去解决样本点的拟合问题，回归多项式如下：
f(tmax,a)=∑p=0,q=0npqt(p)maxa(q) (2) 式中：p、q表示多项式的次数；p = 0，1，2，3；q = 0，1，2，3。
对不同角度冲击下各结构耐撞性指标与结构参数进行拟合，表征拟合效果的系数R2如表2所示。
表 2 设计样本的拟合系数R2Table 2. Fitting coefficient R2 of design samples性能指标 层数 4 5 6 10° 20° 30° 10° 20° 30° 10° 20° 30° e 0.993 8 0.995 3 0.989 9 0.989 1 0.995 2 0.993 5 0.993 3 0.991 2 0.983 3 f 0.993 2 0.994 7 0.996 0 0.992 0 0.993 4 0.991 4 0.991 2 0.993 1 0.995 2 设计样本与多项式回归元模型的拟合度接近1，说明多项式回归元模型可以足够准确地用于e和f的预测以及多目标优化设计，响应面图以5层结构
20∘ 冲击情况为例示于图8。
图 8 20°冲击下的5层结构性能指标响应面Figure 8. Response surfaces of performance indexes of five-storeyed structures under 20° impact为了获得仿生薄壁管的最佳耐撞性能，对仿生薄壁管进行了结构参数的优化设计，需要研究不同结构、不同冲击角度下的综合耐撞性能。因此，需要进行加权平均法计算，公式表示如下：
ei=eemax,fi=ffmax (3) ew=∑ki=1wiei,fw=∑ki=1wifi,∑ki=1wi=1 (4) 式中：ei和fi分别为比吸能和峰值载荷在情况i = 1，2，3下的归一化值；i = 1，2，3分别表示受冲击角度α = 10°，20°，30°；ew和fw分别为比吸能和峰值载荷的加权值；k = 3表示受斜向冲击的3种角度；wi表示权重因子。
如前所述，e和f作为重要的耐撞性指标，在优化模型里作为优化目标，期望e最大化的同时f最小化，约束是最大厚度tmax和厚度梯度值a，数学模型表示如下：
min[fw(tmax,a)−ew(tmax,a)]2.80mm≤tmax≤3.00mm,0.20mm≤a≤0.40mm (5) 4.2 多目标优化方法
为确定仿生薄壁管的最优结构参数组合，本研究以多目标粒子群优化算法为基础，目标函数、约束条件和决策支持为基本要素，将比吸能最大化和峰值载荷最小化作为模型的优化目标，基于Pareto最优原则得到层数n、最大厚度tmax、厚度梯度值a结构参数最优化的配置方案。优化方法的流程图如图9所示。
4.3 优化分析
在多目标优化过程中，设置了3种不同冲击角度的权重因数wi（i = 1，2，3），w1、w2、w3分别为10°、20°、30°工况下的权重因数。表3所示是3种权重因数配置方案，前3种方案为单一冲击角度，后3种方案综合考虑了多角度的影响。
表 3 不同设计方案的权重因数Table 3. Weighting factors for different design cases配置方案 w1 w2 w3 Ⅰ 1 0 0 Ⅱ 0 1 0 Ⅲ 0 0 1 Ⅳ 1/6 1/3 1/2 Ⅴ 1/3 1/3 1/3 Ⅵ 1/2 1/3 1/6 对于仿生薄壁管层数的影响，本文通过加权计算的方法，利用下式：
ejw=∑6j=4μjew,fjw=∑6j=4μjfw (6) 计算得到各种冲击角度下4、5、6层的综合加权性能指标ejw和fjw，考虑到层数的影响程度相同，即配置每层情况的权重因数相同，本文中采用
μj 为1/3。基于多项式回归元模型，通过式（6）计算得到ejw和fjw，得到Pareto前沿如图10所示。图10（a）表示单一冲击角度配置方案的Pareto前沿，当冲击角度从10°增大到30°时，Pareto前沿逐步向右上方移动，意味着在产生相同的峰值载荷时，随着冲击角度的增大，仿生薄壁管的e减小。为了选择单一冲击角度下的最优值，本文中采用最小距离选择法求得极限最大值。首先确定最优指标对应的结构参数tmax-a组合，然后在层数之间通过fi−ei的差值进行组内排名，获得其对应的最优结构设计参数。如表4所示，10°斜向冲击时，n = 6，tmax = 2.84 mm，a = 0.38 mm为吸能效果最好的结构参数，同理，20°、30°斜向冲击时，吸能效果最好的结构参数分别为n = 6，tmax = 2.89 mm，a = 0.29 mm；n = 6，tmax = 2.91 mm，a = 0.34 mm。
表 4 不同设计方案的最优结构设计参数Table 4. Optimum structural design parameters in different design cases配置方案 最优指标 结构设计参数 性能指标 fi−ei 组内排名 ejw fjw tmax/mm a/mm α/(°) n e/(kJ·kg−1) f/kN Ⅰ −0.839 0.804 2.84 0.38 10 4 36.55 75.08 0.151 3 5 39.31 76.64 0.047 2 6 38.50 78.65 −0.064 1 Ⅱ −0.819 0.829 2.89 0.29 20 4 32.06 65.77 0.131 2 5 30.72 60.39 0.244 3 6 31.92 68.00 −0.216 1 Ⅲ −0.803 0.846 2.91 0.34 30 4 22.07 59.25 0.266 2 5 22.41 58.58 0.434 3 6 21.82 59.53 −0.447 1 Ⅳ −0.824 0.816 2.95 0.28 10 4 40.72 71.04 0.010 3 5 42.49 74.11 −0.025 2 6 42.44 75.26 0.024 4 20 4 32.40 66.88 0.138 5 5 30.80 60.95 0.242 6 6 32.07 68.79 −0.213 1 30 4 22.50 55.96 0.216 7 5 22.04 57.16 0.443 9 6 22.90 57.19 0.423 8 Ⅴ −0.817 0.826 2.92 0.30 10 4 46.01 65.32 −0.176 2 5 38.87 64.35 0.057 5 6 38.42 67.75 −0.060 4 20 4 31.47 61.47 0.090 6 5 31.00 63.17 0.237 8 6 36.12 58.42 −0.120 3 30 4 23.10 52.70 0.161 7 5 20.98 48.75 0.466 9 6 24.71 53.39 −0.381 1 Ⅵ −0.830 0.848 2.85 0.33 10 4 37.80 71.82 0.083 5 5 39.05 75.32 0.053 4 6 39.81 76.63 −0.036 3 20 4 31.60 72.34 0.224 7 5 32.81 62.38 0.196 6 6 30.23 61.30 −0.255 2 30 4 20.04 68.63 0.429 9 5 22.79 60.95 0.425 8 6 21.18 62.58 −0.462 1 图10（b）表示多种冲击角度配置方案的Pareto前沿，从配置方案Ⅳ到Ⅵ，Pareto前沿逐步从右向左下方移动，相同f时的e逐渐增大，参考单一冲击角度配置方案，出现这种趋势的原因是小角度的权重因数增加。因此，相同条件下配置方案Ⅵ下的仿生管相比配置方案Ⅳ和Ⅴ下的仿生管的吸能效果好。为了选择多种冲击角度下的最优值，本文中采用最小距离选择法求得极限最大值，首先确定最优指标对应的结构参数tmax-a组合，然后对层数和冲击角度通过fi−ei的差值进行组内排名，获得其对应的最优结构设计参数。最优结构设计参数列于表4。
综合考虑多种冲击角度的影响，在配置不同权重因数的方案中，方案Ⅳ的结构n = 6，tmax = 2.95 mm，a = 0.28 mm吸能效果最好；方案Ⅴ的结构n = 6，tmax = 2.92 mm，a = 0.30 mm吸能效果最好；方案Ⅵ的结构n = 6，tmax = 2.85 mm，a = 0.33 mm吸能效果最好。由此可见，多种冲击角度下的优化结果与权重因数分配方案有很大关系，不同冲击工况下权重因数的选择对仿生薄壁管结构优化设计具有重要意义。
5. 结 论
（1）基于结构仿生学，结合鹿角骨单位耐撞性结构特征，利用结构仿生原理设计了一种内径相同外径等梯度逐层递减的仿生薄壁管。通过有限元仿真和实物试验，探究斜向冲击时仿生薄壁管的吸能特性，发现仿生薄壁管变形模式为逐层渐进叠缩，并未出现欧拉弯曲变形模式，吸能效果良好。
（2）基于多目标粒子群优化算法，采用最小距离选择法对Pareto前沿进行优化分析，在方案Ⅰ、Ⅱ、Ⅲ中，考虑单一冲击角度下的仿生薄壁管耐撞性，得到最优结构参数分别为n = 6，tmax = 2.84 mm，a = 0.38 mm；n = 6，tmax = 2.89 mm，a = 0.29 mm；n = 6，tmax = 2.91 mm，a = 0.34 mm。在方案Ⅳ、Ⅴ、Ⅵ中，考虑多种角度冲击下的仿生薄壁管耐撞性，得到最优结构参数分别为n = 6，tmax = 2.95 mm，a = 0.28 mm；n = 6，tmax = 2.92 mm，a = 0.30 mm；n = 6，tmax = 2.85 mm，a = 0.33 mm。
（3）通过实物试验验证了仿真结果的可靠性，证明了将鹿角骨单位的结构特征应用到薄壁管的设计中可以表现出优异的耐撞性能，优化出的结构参数可为车辆吸能元件的设计提供参考。
-
图 3 仿生薄壁管在10 m/s轴向冲击下的变形模式
Figure 3. Deformation patterns of a bionic tube under the axial impact of 10 m/s
图 4 在186 kg的落锤下落高度为10 m的轴向冲击载荷下仿生薄壁管的力-位移曲线和吸能量-位移曲线
Figure 4. Force-displacement and energy absorption-displacement curves of the bionic tube under the axial impact of a drop hammer of 186 kg droping from the heigh of 10 m
图 8 20°冲击下的5层结构性能指标响应面
Figure 8. Response surfaces of performance indexes of five-storeyed structures under 20° impact
表 1 试验组别
Table 1. Group of test factors
因素 组别 1 2 3 4 5 6 7 8 9 10 11 12 tmax/mm 2.80 2.80 2.80 2.80 2.80 2.85 2.85 2.85 2.85 2.85 2.9 2.90 a/mm 0.20 0.25 0.30 0.35 0.40 0.20 0.25 0.30 0.35 0.40 0.20 0.25 因素 组别 13 14 15 16 17 18 19 20 21 22 23 24 25 tmax/mm 2.90 2.90 2.90 2.95 2.95 2.95 2.95 2.95 3.00 3.00 3.00 3.00 3.00 a/mm 0.30 0.35 0.40 0.20 0.25 0.30 0.35 0.40 0.20 0.25 0.30 0.35 0.40 
下载: 导出CSV
表 2 设计样本的拟合系数R2
Table 2. Fitting coefficient R2 of design samples
性能指标 层数 4 5 6 10° 20° 30° 10° 20° 30° 10° 20° 30° e 0.993 8 0.995 3 0.989 9 0.989 1 0.995 2 0.993 5 0.993 3 0.991 2 0.983 3 f 0.993 2 0.994 7 0.996 0 0.992 0 0.993 4 0.991 4 0.991 2 0.993 1 0.995 2
下载: 导出CSV
表 3 不同设计方案的权重因数
Table 3. Weighting factors for different design cases
配置方案 w1 w2 w3 Ⅰ 1 0 0 Ⅱ 0 1 0 Ⅲ 0 0 1 Ⅳ 1/6 1/3 1/2 Ⅴ 1/3 1/3 1/3 Ⅵ 1/2 1/3 1/6
下载: 导出CSV
表 4 不同设计方案的最优结构设计参数
Table 4. Optimum structural design parameters in different design cases
配置方案 最优指标 结构设计参数 性能指标 fi−ei 组内排名 ejw fjw tmax/mm a/mm α/(°) n e/(kJ·kg−1) f/kN Ⅰ −0.839 0.804 2.84 0.38 10 4 36.55 75.08 0.151 3 5 39.31 76.64 0.047 2 6 38.50 78.65 −0.064 1 Ⅱ −0.819 0.829 2.89 0.29 20 4 32.06 65.77 0.131 2 5 30.72 60.39 0.244 3 6 31.92 68.00 −0.216 1 Ⅲ −0.803 0.846 2.91 0.34 30 4 22.07 59.25 0.266 2 5 22.41 58.58 0.434 3 6 21.82 59.53 −0.447 1 Ⅳ −0.824 0.816 2.95 0.28 10 4 40.72 71.04 0.010 3 5 42.49 74.11 −0.025 2 6 42.44 75.26 0.024 4 20 4 32.40 66.88 0.138 5 5 30.80 60.95 0.242 6 6 32.07 68.79 −0.213 1 30 4 22.50 55.96 0.216 7 5 22.04 57.16 0.443 9 6 22.90 57.19 0.423 8 Ⅴ −0.817 0.826 2.92 0.30 10 4 46.01 65.32 −0.176 2 5 38.87 64.35 0.057 5 6 38.42 67.75 −0.060 4 20 4 31.47 61.47 0.090 6 5 31.00 63.17 0.237 8 6 36.12 58.42 −0.120 3 30 4 23.10 52.70 0.161 7 5 20.98 48.75 0.466 9 6 24.71 53.39 −0.381 1 Ⅵ −0.830 0.848 2.85 0.33 10 4 37.80 71.82 0.083 5 5 39.05 75.32 0.053 4 6 39.81 76.63 −0.036 3 20 4 31.60 72.34 0.224 7 5 32.81 62.38 0.196 6 6 30.23 61.30 −0.255 2 30 4 20.04 68.63 0.429 9 5 22.79 60.95 0.425 8 6 21.18 62.58 −0.462 1
下载: 导出CSV
-
[1] 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. [2] 徐峰祥, 张锁, 武昆迎. 厚度幂指数分布管状结构耐撞性设计准则与方法研究 [J]. 爆炸与冲击, 2019, 39(3): 035103. DOI: 10.11883/bzycj-2018-0013.XU F X, ZHANG S, WU K Y. Study on crashworthiness design criteria and method of tubular structures with power exponent distribution of thickness [J]. Explosion and Shock Waves, 2019, 39(3): 035103. DOI: 10.11883/bzycj-2018-0013. [3] 郝文乾, 卢进帅, 黄睿, 等. 轴向冲击载荷下薄壁折纹管的屈曲模态与吸能 [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. [4] 谭丽辉, 徐涛, 崔晓梅, 等. 带有圆弧形凹槽金属薄壁圆管抗撞性优化设计 [J]. 爆炸与冲击, 2014, 34(5): 547–553. DOI: 10.11883/1001-1455(2014)05-0547-07.TAN L H, XU T, CUI X M, et al. Design optimization for crashworthiness of metal thin-walled cylinders with circular arc indentations [J]. Explosion and Shock Waves, 2014, 34(5): 547–553. DOI: 10.11883/1001-1455(2014)05-0547-07. [5] 张涛, 刘土光, 肖汉林, 等. 高速冲击下薄壁组合结构吸能特性研究 [J]. 爆炸与冲击, 2006, 26(5): 395–403.ZHANG T, LIU T G, XIAO H L, et al. Energy absorption performance of thin-walled structures with triggering holes subjected to high-speed axial impact [J]. Explosion and Shock Waves, 2006, 26(5): 395–403. [6] 殷之平, 李玉龙, 黄其青. 含诱导缺陷薄壁圆管耐撞性优化设计 [J]. 爆炸与冲击, 2011, 31(4): 418–422. DOI: 10.11883/1001-1455(2011)04-0418-05.YIN Z P, LI Y L, HUANG Q Q. Optimal crashworthiness design of thin-walled circular tubes with triggering holes [J]. Explosion and Shock Waves, 2011, 31(4): 418–422. DOI: 10.11883/1001-1455(2011)04-0418-05. [7] 李松晏, 郑志军, 虞吉林. 高速列车吸能结构设计和耐撞性分析 [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. [8] SUN G Y, LIU T Y, FANG J G, et al. Configurational optimization of multi-cell topologies for multiple oblique loads [J]. Structural and Multidisciplinary Optimization, 2018, 57(2): 469–488. DOI: 10.1007/s00158-017-1839-5. [9] ALKHATIB S E, TARLOCHAN F, HASHEM A, et al. Collapse behavior of thin-walled corrugated tapered tubes under oblique impact [J]. Thin-Walled Structures, 2018, 122: 510–528. DOI: 10.1016/j.tws.2017.10.044. [10] ASANJARANI A, DIBAJIAN S H, MAHDIAN A. Multi-objective crashworthiness optimization of tapered thin-walled square tubes with indentations [J]. Thin-Walled Structures, 2017, 116: 26–36. DOI: 10.1016/j.tws.2017.03.015. [11] 亓昌, 董方亮, 杨姝, 等. 锥形多胞薄壁管斜向冲击吸能特性仿真研究 [J]. 振动与冲击, 2012, 31(24): 102–107. DOI: 10.13465/j.cnki.jvs.2012.24.009.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.13465/j.cnki.jvs.2012.24.009. [12] 任露泉, 梁云虹. 仿生学导论[M]. 北京: 科学出版社, 2016: 208-210. [13] HUANG H, XU S C. Crashworthiness analysis and bionic design of multi-cell tubes under axial and oblique impact loads [J]. Thin-Walled Structures, 2019, 144: 106333. DOI: 10.1016/j.tws.2019.106333. [14] 许述财, 邹猛, 魏灿刚, 等. 仿竹结构薄壁管的轴向耐撞性分析及优化 [J]. 清华大学学报(自然科学版), 2014, 54(3): 299–304. DOI: 10.16511/J.CNKI.QHDXXB.2014.03.007.XU S C, ZOU M, WEI C G, et al. Axial crashworthiness analysis and optimization of a bionic thin-walled tube based on bamboo structure [J]. Journal of Tsinghua University (Science and Technology), 2014, 54(3): 299–304. DOI: 10.16511/J.CNKI.QHDXXB.2014.03.007. [15] ZOU M, XU S C, WEI C G, et al. A bionic method for the crashworthiness design of thin-walled structures inspired by bamboo [J]. Thin-Walled Structures, 2016, 101: 222–230. DOI: 10.1016/j.tws.2015.12.023. [16] 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. [17] LI Z, DUAN L B, CHEN T, et al. Crashworthiness analysis and multi-objective design optimization of a novel lotus root filled tube (LFT) [J]. Structural and Multidisciplinary Optimization, 2018, 57(2): 865–875. DOI: 10.1007/s00158-017-1782-5. [18] YIN H F, XIAO Y Y, WEN G L, et al. Multi-objective robust optimization of foam-filled bionic thin-walled structures [J]. Thin-Walled Structures, 2016, 109: 332–343. DOI: 10.1016/j.tws.2016.10.011. [19] CURREY J D. Bones: structure and mechanics [M]. New Jersey: Princeton University Press, 2006. [20] PICAVET P P, BALLIGAND M. Organic and mechanical properties of Cervidae antlers: a review [J]. Veterinary Research Communications, 2016, 40(3-4): 141–147. DOI: 10.1007/s11259-016-9663-8. [21] FANG Z Q, CHEN B, LIN S Y, et al. Investigation of inner mechanism of anisotropic mechanical property of antler bone [J]. Journal of the Mechanical Behavior of Biomedical Materials, 2018, 88: 1–10. DOI: 10.1016/j.jmbbm.2018.07.043. [22] CURREY J D, LANDETE-CASTILLEJOS T, ESTEVEZ J, et al. The mechanical properties of red deer antler bone when used in fighting [J]. Journal of Experimental Biology, 2009, 212(24): 3985–3993. DOI: 10.1242/jeb.032292. [23] DAVISON K S, SIMINOSKI K, ADACHI J D, et al. Bone strength: the whole is greater than the sum of its parts [J]. Seminars in Arthritis and Rheumatism, 2006, 36(1): 22–31. DOI: 10.1016/j.semarthrit.2006.04.002. [24] HENSHAW J. Antlers: the unbrittle bones of contention [J]. Nature, 1971, 231(5303): 469. DOI: 10.1038/231469a0. [25] CHEN P Y, STOKES A G, MCKITTRICK J. Comparison of the structure and mechanical properties of bovine femur bone and antler of the North American elk (cervus elaphus canadensis) [J]. Acta Biomaterialia, 2009, 5(2): 693–706. DOI: 10.1016/j.actbio.2008.09.011. [26] 杨欣, 范晓文, 许述财, 等. 仿虾螯结构薄壁管设计及耐撞性分析 [J]. 爆炸与冲击, 2020, 40(4): 043301. DOI: 10.11883/bzycj-2019-0280.YANG X, FAN X W, XU S C, et al. Design and crashworthiness analysis of thin-walled tubes based on a shrimp chela structure [J]. Explosion and Shock Waves, 2020, 40(4): 043301. DOI: 10.11883/bzycj-2019-0280. 期刊类型引用(13)
1. 李世岳，胡洪林，陈哲伦. 薄壁结构在汽车吸能盒中的应用与展望. 汽车文摘. 2024(02): 42-49 . 
百度学术2. 王海文，徐丹，陈墨林，叶仁传，任鹏. 冰雹形状及重复加载对叶表皮细胞壁仿生结构撞击损伤影响分析. 振动与冲击. 2024(03): 227-232+240 . 百度学术3. 杨欣，王阳，宋家锋，朱勇，黄彬兵，许述财. 基于虾螯结构的仿生夹层板设计及数值模拟. 吉林大学学报(工学版). 2024(03): 842-851 . 百度学术4. 陈雅婷，洪熠豪，吴俊，吴菲，郑氏韦，赵顺秋. 仿生多胞薄壁结构冲击吸能特性研究. 机械强度. 2024(03): 750-755 . 百度学术5. 林启权，刘雯婷，王凯，董文正，王镇柱，贺礼财. 多胞薄壁结构优化及结构参数对耐撞性影响分析. 振动与冲击. 2024(16): 60-66+110 . 百度学术6. 李兴山，刘殿坤，吕玉山，李昌璐. 五轴数控机床床身的全植物叶片结构仿生设计. 机床与液压. 2023(16): 114-119 . 百度学术7. 王勇辉，何强，顾航. 仿荷叶脉络薄壁管耐撞性能研究. 江苏科技大学学报(自然科学版). 2022(01): 59-66 . 百度学术8. 何强，王勇辉，顾航. 爆炸载荷作用下仿玉莲脉络加筋板的动态响应分析. 江苏科技大学学报(自然科学版). 2022(02): 15-21 . 百度学术9. 秦晓宇，马其华，胡沛源，查一斌. 准静态轴压载荷下Al-CFRP梯度混合管的多目标优化设计. 塑性工程学报. 2022(09): 188-198 . 百度学术10. 罗耿，赵剑南，陈亮，柴成鹏，王童，陈轶嵩. 新型波纹结构吸能盒耐撞性研究. 汽车安全与节能学报. 2022(03): 453-462 . 百度学术11. 刘飞明，雷建银，乔力，刘志芳. 木贼属植物仿生薄壁结构的耐撞性优化. 高压物理学报. 2022(05): 141-151 . 百度学术12. 李响，蔡明杰，徐兴兴，周绍国，焦元辰. 形内自相似层级类蜂窝面外冲击特性研究. 河北科技大学学报. 2022(05): 481-494 . 百度学术13. 牛枞，黄晗，向枳昕，闫庆昊，陈金宝，许述财. 仿生多胞薄壁管耐撞性分析及优化. 爆炸与冲击. 2022(10): 169-182 . 
本站查看其他类型引用(19)
-
下载:
百度学术
计量
- 文章访问数: 1752
- HTML全文浏览量: 1042
- PDF下载量: 98
- 被引次数: 32