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
-
高压下材料屈服强度特性是当前爆炸力学、冲击动力学、高压物理和材料科学等学科领域重点关注的基础问题,研究成果在装甲防护、航空航天器防护等方面有重要的应用价值。自20世纪60年代以来,先后发展了静水压线比较法[1]、压剪方法[2]、横向应力计方法[3]和自洽方法[4](也称为双屈服面方法或AC方法)等高压强度研究方法。静水压线比较法受状态方程精度影响,在高压下精度并不高。压剪方法和横向应力计方法分别受到加载技术和测试技术的限制,往往局限于20 GPa以内。自洽方法原则上没有压力限制,但需要进行一组冲击加载-卸载实验和冲击加载-再加载实验才能获取一个强度数据。其中冲击加载-再加载实验相对复杂,需要采用由高、低阻抗材料构成的组合飞片进行冲击才能实现。由于组合飞片在气炮驱动的过程中容易发生分离,导致在待测材料中形成冲击加载-卸载-再冲击加载,而非预期的冲击加载-再加载,也就无法获得强度数据。在现有的高压强度研究方法中,自洽方法是目前应用最广泛的一种方法,提供了包括金属材料(如铝及合金[4-10]、铜[6]、钒[11]、钨[12-14]、钽[15]、铍[16]、钢[17-18]等)、陶瓷材料(B4C[19]、SiC[20])、金属玻璃(Zr51Ti5Ni10Cu25Al9)[21]和复合材料[22]等大量材料的高压强度数据。但正如前面所指出的,由于冲击加载-再加载实验难度大,目前自洽方法提供的高压强度数据大部分仅来自冲击加载-卸载实验,这类强度数据实际上并不完备。
针对上述自洽方法的问题,胡建波等[9-10]通过采用平面焊接技术加强组合飞片之间的结合力来避免其分离。为了简化波系作用,利于数据处理分析,组合飞片的第一层通常为待测样品材料。胡建波等[9-10]利用由此制作的TC4/LY12铝组合飞片成功地加速到3.2 km/s,并对LY12铝样品进行了对称碰撞,获得了34 GPa压力下(指第一次冲击压力,下同)LY12铝的冲击加载-再加载粒子速度剖面。在此之前,铝的冲击加载-再加载实验的最高压力仅为22 GPa[5]。类似地,M.D.Furnish等[18]则采用爆炸焊接技术制作了Ta/2169钢组合飞片。该飞片被成功地加速到2.0 km/s,并由此获得了48 GPa压力下2169钢的冲击加载-再加载粒子速度剖面。上述焊接技术的应用极大提升了冲击加载-再加载实验能力,拓展了自洽方法的实验压力范围。然而,上述焊接技术比较复杂,更重要的是焊接处理时会经历高温、高压作用过程,可能改变组合飞片中的待测样品材料初始状态,从而对高压强度结果造成影响。值得注意的是,大量的轻气炮实验数据表明:在轻气炮驱动飞片加速过程中,飞片会发生如弓形或马鞍形等复杂的弯曲变形,而且不同材料飞片的变形量及形状也不同[23-24]。事实上,上述因素极易造成组合飞片在气炮驱动的过程中发生分离。
基于上述分析,本文中提出一种实现冲击加载-再加载实验的简便方法—采用较高硬度材料为支撑,通过环氧树脂与待测样品粘接制成组合飞片,由此减小气炮加载下飞片的弯曲变形来避免分离,并通过铝、锡和锆基金属玻璃的冲击实验验证该方法的有效性。在此基础上,获得铝、锡和锆基金属玻璃再加载过程剪应力的变化情况,并初步分析其对强度结果的影响。
1. 实验方法
实验采用如图 1所示的反向碰撞方式,即由组合飞片直接撞击透明的单晶LiF窗口。组合飞片的第一层为待测材料(即样品),其后为用于支撑样品的具有较高硬度的材料,两层飞片之间采用环氧树脂进行粘接,环氧树脂厚度约10 μm。上述组合飞片的制作过程均在常温、常压下进行,不会改变组合飞片中的待测样品材料初始状态。LiF窗口尺寸为Ø28 mm×12 mm,碰撞面镀有1 μm厚的铝膜作为光学测试的反射面,同时为保护长历时测量过程中铝膜不受破坏,铝膜前粘接了8 μm铜箔。一种激光干涉测速技术—DISAR(displacement interferometer system for any reflector)技术[25]用于测量样品/LiF窗口界面粒子速度剖面,通过剖面可以直观判断样品中是否经历预期的冲击加载-再加载过程。组合飞片冲击LiF窗口的速度由磁测速技术进行测量。根据实测的飞片速度以及样品和窗口材料的Hugoniot参数,采用阻抗匹配法可以计算得到样品的冲击压力和冲击波速度,结合图 1所示的波系作用可以进一步计算得到再加载弹塑性波的声速cL。
组合飞片中支撑材料除具有较高硬度外,其阻抗应高于样品材料的阻抗,以实现对样品的冲击加载-再加载。当然,为确保再加载过程出现弹塑性特征信息,再加载压力幅度不能过高,以避免形成冲击波。根据上述条件,针对LY12铝样品采用了TC4钛合(成分为Ti90-Al6-V4)作为支撑材料,锡和锆基金属玻璃(Zr51Ti5Ni10Cu25Al9,详细参数参见文献[26])则采用45钢作为支撑材料。
2. 实验结果与初步分析
利用图 1所示的实验装置,在Ø30 mm二级轻气炮上进行了5发铝、锡和锆基金属玻璃为样品的冲击加载-再加载实验,冲击速度范围为2.49~4.39 km/s,在样品中产生的压力范围为28~48 GPa,样品和支撑材料的直径均为28 mm, 详细实验参数见表 1, 表中Hs为样品厚度,hs为支撑材料厚度,D为冲击波速度,p为冲击压力,τH和τc分别为Hugoniot状态剪应力和临界剪应力。为了尽量减小气炮加载下组合飞片的弯曲变形,支撑材料总厚度不小于3 mm。尤其是在最高冲击速度为4.39 km/s的实验中(实验2),除2 mm厚的TC4外另增加了3 mm厚的钽作为支撑材料。
表 1 平板冲击实验条件及结果Table 1. Experimental conditions and results for planar plate-impact experiments实验编号 样品材料 Hs/mm 支撑材料 hs/mm Ds/(km·s-1) p/GPa (τc-τH)/GPa 1 LY12铝 1.445 TC4 3.01 3.67 38.3 0.73 2 LY12铝 1.465 Ta/TC4 3.05/2.09 4.39 48.5 0.77 3 锡 2.013 45钢 4.50 2.49 28.7 0.07 4 锡 2.015 45钢 4.50 3.08 38.1 0.16 5 锆基金属玻璃 3.135 45钢 4.50 3.00 39.1 0.53 由DISAR测得的铝、锡和锆基金属玻璃为样品的冲击加载-再加载粒子速度剖面如图 2~4所示。从图中可以看到,5发实验均获得了预期的冲击加载-再加载波剖面,这表明采用较高硬度材料为支撑以减小组合飞片弯曲变形来避免分离的方法是简便、有效的。此外,铝、锡和锆基金属玻璃为样品的再加载波剖面均出现一定程度的弹塑性波剖面特征,表明相对应支撑材料的阻抗选取是合适的。图中uw, H和uw, 1分别为再加载弹性波和塑性波起始对应的粒子速度。
根据图 1所示的波传播特性,由图 2~4的粒子速度剖面可以计算得到沿着再加载过程的拉格朗日纵波声速:
cL=HsΔt−Hs/Ds (1) 式中:Ds为碰撞时样品中产生的冲击波的速度(由实测的飞片速度以及样品和窗口材料的Hugoniot参数,通过阻抗匹配法计算得到冲击波速度),Δt为如图 1所示的来自样品/支撑材料界面的再加载波到达样品/窗口界面时间与碰撞时间之差。
图 5给出了LY12铝再加载过程拉格朗日纵波声速随粒子速度u变化的典型结果(实验1)。其中,粒子速度u由样品/窗口界面粒子速度uw结合增量型阻抗匹配法[12, 27]计算得到,由此得到的粒子速度计及了加载波在样品/窗口界面反射造成的影响。uH和u1分别为Hugoniot状态对应的粒子速度和再加载进入塑性屈服时对应的粒子速度,由图 2对应的uw, H和uw, 1根据增量型阻抗匹配法计算得到。图 5中同时给出了由LY12铝冲击加载-卸载粒子速度剖面计算得到的卸载过程拉格朗日纵波声速的结果[28]。从图中可以看到,再加载过程也呈现与卸载过程相同的准弹性行为特征,即弹、塑性波速为光滑过渡,而没有发生突降。如图中虚线所示,将卸载过程塑性段声速线性外延可得到再加载过程的拉格朗日体波声速cb。
在获得纵波和体波的基础上,根据J.R.Asay等[4]提出的自洽强度方法,对再加载过程的声速进行如下计算可得到剪应力的变化:
τc−τH=34ρ0∫u1uHc2L−c2bcLdu (2) 式中:uH和u1分别为图 5所示的Hugoniot状态对应的粒子速度和再加载进入塑性屈服时对应的粒子速度,ρ0为材料的初始密度,cL和cb分别为图 5所示的再加载过程对应的拉格朗日纵波和体波声速。由此计算得到的LY12铝的τc-τH结果见表 1。结合胡建波等[29]和俞宇颖等[21]分别给出的锡和锆基金属玻璃卸载过程声速数据,采用上述相同方法计算得到了锡和锆基金属玻璃再加载过程的τc-τH,计算结果见表 1。
在获得τc-τH的基础上,如能获得相同冲击压力下卸载过程的声速则可计算得到τc+τH,进而可以确定在该冲击压力下材料的屈服强度Y=(τc-τH)+(τc+τH)。由于冲击加载-再加载实验的困难,先前大量文献中将屈服强度Y ≈τc+τH,即认为再加载过程的τc-τH可以忽略。就LY12铝而言,本文再加载实验得到的38.3和48.5 GPa冲击压力下的τc-τH分别为0.73和0.77 GPa;胡建波等[10]通过卸载实验得到的32.2和54.7 GPa冲击压力下的τc+τH分别为0.85和1.05 GPa。显然,在上述压力范围内仅由冲击加载-卸载实验得到的LY12铝强度Y≈τc+τH将比实际结果Y=2τc降低约45%。对于锡,本文再加载实验得到的28.7和38.1 GPa冲击压力下的τc-τH分别为0.08和0.19 GPa;依据胡建波等[29]的卸载实验剖面数据计算得到的28.3和39.2 GPa冲击压力下的τc+τH分别为0.07和0.17 GPa。仅由冲击加载-卸载实验得到的锡强度Y≈τc+τH将比实际结果Y=2τc降低约50%。同样地,就锆基金属玻璃而言,本文再加载实验得到的39.1 GPa冲击压力下的τc-τH为0.53 GPa,俞宇颖等[21]通过卸载实验得到的37.3 GPa冲击压力下的τc+τH为1.73 GPa,仅由冲击加载-卸载实验得到的强度Y≈τc+τH将比实际结果Y=2τc降低约20%。综上所述,在采用自洽方法计算高压强度时冲击加载-再加载数据不可或缺。
3. 结论
针对自洽高压强度方法存在的因组合飞片分离而难以实现冲击加载-再加载的难题,提出了一种简便方法—采用较高硬度材料为支撑,通过环氧树脂与待测样品粘接制成组合飞片,由此减小气炮加载下飞片的弯曲变形来避免分离。采用该方法,在二级轻气炮上进行了冲击速度2.49~4.39 km/s、冲击压力28~48 GPa范围内铝、锡和锆基金属玻璃为样品的验证实验,5发实验均获得了较理想的冲击加载-再加载粒子速度剖面,表明了该方法的有效性。由本文获得的冲击加载-再加载粒子速度剖面,根据自洽方法计算得到了铝、锡和锆基金属玻璃再加载过程剪应力变化数据。结合已有的冲击加载-卸载过程剪应力变化数据的分析表明,在本文涉及的压力范围内,仅由冲击加载-卸载实验得到的铝、锡和锆基金属玻璃屈服强度将比实际结果降低约20%~50%。因此,在采用自洽方法计算高压强度时冲击加载-再加载数据不可或缺。
-
表 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 表 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 表 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 表 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 -
[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. 期刊类型引用(0)
其他类型引用(1)
-