陈书剑,肖守讷,朱涛,缪得祥,杨冰,阳光武
(西南交通大学牵引动力国家重点实验室,四川成都,610031)
薄壁金属管塑性变形缓冲器因其结构简单,缓冲力平稳,环境温度适应性强等优点,被广泛应用于各种需要缓冲防护的场合。例如,神舟五号载人飞船返回舱座椅设有薄壁金属管缓冲结构[1]。在列车被动安全性防护中,起关键作用的吸能防爬装置也常采用该结构[2]。薄壁金属管塑性变形缓冲器依靠扩径变形过程中的弹塑性变形和摩擦发热来吸收冲击能量[3],许文斌等[4]将该冲击扩径过程视为塑性成形问题,提出了薄壁金属管在准静态加载下缓冲力的理论计算公式。胡念军等[5]根据塑性流动理论,推导出了圆管扩径过程厚向压变形与环向拉变形之比的计算方法。罗昌杰等[6−7]通过试验方法研究了薄壁铝管、泡沫金属塑性变形缓冲器在准静态作用下的缓冲特性,发现对薄壁金属管进行预处理可有效提高其理想吸能效率。吴鸿超等[8]采用试验和数值模拟的方法,从冲击速度和结构参数的角度出发,研究了薄壁管在中高速冲击下的缓冲特性,但未对薄壁管材料因素加以考虑,且目前国内外大部分研究仍集中在薄壁管准静态吸能特性上[9−10]。在大多数缓冲过程中,由于冲击速度较大,缓冲载荷表现出明显的动态行为,薄壁管塑性变形的过程中往往伴随着明显的应变率效应[11],而目前关于材料应变率效应对该类缓冲器吸能特性影响的研究甚少。孟广成等[12-13]对结构的碰撞进行仿真研究,发现考虑应变率效应的仿真结果与实验结果更符合。尤其是在冲击问题中,材料应变率效应对缓冲结构的影响不可忽视。曾必强等[14]的研究表明,材料应变率强化效应会提高结构的承载能力和结构的吸能能力,但具有应变率负敏感性的5083H111 铝合金会降低结构的吸能量[15]。这说明不同材料的应变率效应对结构缓冲特性的影响不同。本文作者以薄壁金属管塑性变形缓冲器为研究对象,对金属管进行不同应变率下的材料试验,获得5083P-O和6008-T6 铝合金在不同应变率下的材料本构关系;结合冲击物理试验与理论模型验证薄壁管缓冲器有限元模型的准确性,基于该模型研究不同材料的应变率效应对薄壁管在动态冲击条件下缓冲特性的影响,分析产生差异的原因,以期为工程实践中吸能结构设计提供参考。
试验材料为薄壁金属管缓冲器用5083P-O和6008-T6 铝合金,其主要化学成分分别如表1和表2所示,对于材料热处理状态,5083P-O 为非热处理铝合金,而6008-T6铝合金则进行了固溶处理加完全人工时效。
在各应变率下分别进行准静态拉伸和动态拉伸试验,准静态试验在RPL100电子式蠕变疲劳试验机上进行;0.5~100 s−1中应变率试验在ZWICK HTM5020 动态拉伸试验机上进行;高应变率试验则采用分离式霍普金森SHTB 试验系统。5083P-O铝合金和6008-T6 铝合金塑性段动态应力−应变曲线如图1和图2所示。
表1 5083P-O铝合金的化学成分(质量分数)Table 1 Chemical composition of aluminum alloy 5083P-O %
表2 6008-T6铝合金的化学成分(质量分数)Table 2 Chemical composition of aluminum alloy 6008-T6 %
图1 5083P-O铝合金准静态和动态拉伸试验结果Fig.1 Quasi-static and dynamic tensile test results of 5083P-O aluminum alloy
图2 6008-T6铝合金准静态和动态拉伸试验结果Fig.2 Quasi-static and dynamic tensile test results of 6008-T6 aluminum alloy
由图1可知:随着应变率的增加,在同一塑性应变ε下,5083P-O 铝合金的流动应力σ先显著降低,存在应变率弱化效应;但当应变率超过320 s−1后,5083P-O铝合金又表现出应变率强化效应,即该种铝合金存在中低应变率范围内先弱化,高应变率范围内再强化的应变率效应,这与已有的研究成果相似[16−18]。这种反向的应变率敏感性相变机制,可用位错钉扎机制模型来解释[19]。由图2 可知:6008-T6铝合金在同一塑性应变下的流动应力随着应变率的增加而增加,表现出明显的应变率强化效应,这是因为铝合金塑性变形机理主要由热激活过程控制[19]。图1和图2的试验结果可为后续有限元计算提供材料参数。
首先,采用文献[8]中试验结构尺寸、材料参数建立有限元模型,单元为均匀六面体网格。锥体以50 m/s的冲击速度轴向冲击薄壁金属管,有限元计算模型如图3所示。然后,分别从试验结果和理论解2个方面对比来验证本文有限元计算方法的可靠性。
图3 有限元模型Fig.3 Finite element model
吴鸿超[20]提出考虑应变率效应和惯性效应的薄壁金属管缓冲力模型,该模型适于计算锥体和金属管相互作用趋于稳定时的缓冲力。该模型采用Cowper-Symonds方程:
式中:σd为动态流动应力;e为应变率;σs为相应的静态应力;D和q为材料常数。
假设材料为理想刚塑性,视锥体为刚体,应用Mises屈服准则,缓冲力Fd可表示为如下形式[2]:
式中:R1和R2分别为金属管在扩径前、后半径;t为金属管径向厚度;α为金属管半锥角;f为金属管与锥台间的摩擦因数;l为金属管扩径段长度;v为冲击速度;ρ为材料密度。根据文献[8]中所示结构尺寸与材料参数得到文献[20]中理论模型结果,将该理论结果与试验、有限元结果中的加速度a进行对比,如图4所示。由图4可知:理论模型未考虑初期缓冲力的复杂性和瞬时性,考虑的是金属管稳定缓冲区的缓冲力,理论结果为定值[20]。有限元结果的平均值与该理论结果的相对误差很小,验证了有限元方法的准确性。
图4 锥体加速度曲线对比Fig.4 Comparison of cone acceleration curves
冲击锥体加速度的仿真结果与试验结果变化趋势整体一致,仿真结果的平均加速度高于试验值,相对误差为5.2%。试验结果与仿真结果在初期阶段变化趋势不同,主要原因是试验所采用的空气炮测试系统为逆弹道冲击系统,冲击体速度由零增加到所需速度,试验结果与仿真结果存在差异[8]。另一方面,文献[8]中的材料参数由试验数据拟合得到,与实际材料力学性能不完全一致,也会产生一定误差。但相对误差较小,说明本文有限元计算方法是可靠的,并且可以用于工程仿真研究。
参考常用薄壁金属管塑性变形缓冲器结构尺寸[6],建立有限元模型,如图5所示。单元相关设置与模型验证时的保持一致。为充分利用试验数据并保证材料参数的模拟精度,采用列表插值法[15]。边界条件设置如下:薄壁金属管底部固定不动,锥体分别以恒速5,20和40 m/s轴向冲击薄壁使其扩径变形。
将2种铝合金材料在20 m/s冲击速度下的缓冲力−位移曲线进行对比分析,如图6所示。由图6可知:无论是否考虑材料的应变率效应,同一材料缓冲器的缓冲力随锥体位移增加的变化形式基本一致,说明材料应变率效应对缓冲力历程无明显影响。考虑材料应变率效应会使瞬时缓冲力和平均缓冲力明显整体降低或升高。结合2种材料在不同应变率下的应力-应变曲线可知:由于在该冲击速度下,5083P-O铝合金的材料应变率敏感性明显比6008-T6铝合金的高,所以,对应结构的缓冲力变化更为剧烈。
图5 缓冲器有限元模型Fig.5 Finite element model of buffer
图6 冲击速度为20 m/s时缓冲力−位移曲线对比Fig.6 Comparison of load-displacement curve under 20 m/s of impact speed
缓冲器的吸能特性计算结果如表3所示。由表3可知:在考虑材料应变率效应时,各冲击速度下5083P-O铝合金材料均表现出应变率弱化效应,材料的流动应力降低,平均应变率增大,导致整体结构中每个单元产生相同应变的时间减少,从而减弱了薄壁管的抗冲击能力,所以,平均缓冲力相比不考虑应变率效应时明显降低;6008-T6铝合金材料则表现出应变率强化效应,使得平均缓冲力明显升高。
对于5083P-O铝合金制成的薄壁管,其平均缓冲力缩小因子随冲击速度的增大而逐渐降低,这是由于随着冲击过程应变率的增大,5083P-O材料的流动应力弱化行为逐渐减弱,冲击端部的抵抗能力增加,所以,与不考虑材料应变率效应时的平均缓冲力逐渐接近。对于6008-T6铝合金制成的薄壁管,其平均缓冲力放大因子随冲击速度的增大而增大,与不考虑材料应变率效应时的平均缓冲力逐渐远离。
结合2种材料的动态应力应变曲线可知:不同速度下结构的缓冲力缩放趋势与2种材料的流动应力随应变率的变化趋势一致。
总吸能量是薄壁管塑性变形缓冲器最主要的性能评价指标。由能量守恒定律可将缓冲器总吸能量表示为
表3 缓冲器吸能特性计算结果Table 3 Calculation results of energy absorption characteristics of buffer
式中:S为锥体位移;E为总吸能量。
总吸能量为塑性变形吸能量和摩擦耗能量之和。不同材料总吸能对比如图7所示。由图7 可知:缓冲器被压缩位移一定时,无论是否考虑应变率效应,总吸能量随着冲击速度的增大而增大。考虑应变率效应时,总吸能量随冲击速度增大的敏感性更高。
各部分能量的变化见表3。由表3 可知:无论是否考虑应变率效应,缓冲器的塑性变形吸能量均大于摩擦耗能,在吸能机制中起到最重要的作用。
图7 不同材料总吸能量对比Fig.7 Total energy absorption comparison of different materials
在考虑材料应变率效应时,一方面,5083P-O铝合金薄壁管的总吸能量显著降低,其中,塑性变形吸能量降低比摩擦耗能降低更明显;另一方面,6008-T6 铝合金薄壁管的总吸能量显著提高,塑性变形吸能量提高比摩擦耗能提高更明显。这说明材料应变率效应对缓冲器塑性变形吸能量的影响更大。
锥体扩径行为分为不定径和定径2个阶段,下面分析不定径段薄壁管的变形模式,此时薄壁管发生初始塑性变形,其不定径段变形模式的剖视图分别如图8和图9所示,最大塑性应变如表4所示。
由表4 可知:对于5083P-O 薄壁管,当冲击速度较低时,初期塑性应变增大;当冲击速度较高时,初期塑性应变反而减小;对于6008-T6 薄壁管,初期塑性应变随冲击速度的增大而单调减小。随着不定径段半径的减小,冲击锥体受到阻力增大,此时对应的缓冲力急剧增大,因此,较大的塑性应变集中于不定径段末端。
图8 5083P-O铝合金不定径段变形模式对比Fig.8 Comparison of deformation modes of tapered section for 5083P-O aluminium alloy
图9 6008-T6铝合金不定径段变形模式对比Fig.9 Comparison of deformation modes of tapered section for 6008-T6 aluminium alloy
表4 薄壁管不定径段最大塑性应变Table 4 Maximum plastic strain of tapered section of thin-walled tubes ×10−4
当低速冲击时,扩径变形应变率较小,5083PO铝合金呈应变率弱化效应,流动应力有一定程度减小,削弱了不定径段抵抗弯曲变形的能力,使结构整体鼓胀变形增大;但当冲击速度达40 m/s时,取母线上典型位置处单元应变率的统计平均值表征结构不定径段应变率,如图10所示。由图10 可知:扩径变形的应变率明显增大,超过了320 s−1,此时间段5083P-O 铝合金表现出应变率强化效应,材料流动应力增强,增大了不定径段抵抗弯曲变形的能力,使结构整体鼓胀变形减小。
在整个有效吸能过程中,对于6008-T6 铝合金,材料流动应力与应变率呈正相关,随着冲击速度增大,扩径应变率增大,6008-T6铝合金呈应变率强化效应,流动应力随应变率单调增加,提高了整体结构的抗冲击能力。
图10 5083P-O铝合金不定径段应变率Fig.10 Strain rate of tapered section of 5083P-O aluminum alloy
理想吸能效率是评价缓冲器吸能特性的重要指标。
式中:L为缓冲器的理想吸能效率;Fmax为整个缓冲过程中最大峰值力。吸能效果越好,理想吸能效率越接近100%。
根据缓冲器变形机理,为提高缓冲器理想吸能效率,可将金属管的不定径部分按照由厚到薄进行加工。在准静态加载条件下,该类缓冲器理想吸能效率得到有效提高[5]。为了研究在动态加载条件下材料应变率效应对该预处理结构的影响。参考文献[5]中的预处理方式对原有限元模型不定径段进行均匀减薄,如图11所示。计算的边界条件、单元相关设置与模型验证时的设置保持一致。
图11 有限元模型不定径段示意图Fig.11 Diagram of tapered section of FEM
不同冲击速度下未预处理与预处理缓冲器的缓冲力−位移曲线如图12所示。由图12 可知:在动态冲击作用下,经预处理的缓冲器载荷历程发生改变,虽然总吸能量有所降低,但初期峰值力显著降低,且缓冲力随位移的变化趋于平缓,达到稳态力前不会出现明显波动,这表明即使在动态加载条件下,该预处理方法也能有效降低冲击过程中的载荷波动。为进一步研究应变率效应对结构吸能特性的影响,统计各项评价指标,如表5所示。由表5 可知:在动态冲击载荷下,无论是否考虑应变率效应,冲击结束时预处理结构与原结构的最大峰值力在同一水平,但预处理结构的总吸能量降低,所以预处理结构的理想吸能效率在一定程度上低于原结构的理想吸能效率。
对于5083P-O 结构,考虑应变率效应时,缓冲器预处理前、后的理想吸能效率差值明显降低,这表明5083P-O 铝合金的材料应变率效应会减少因预处理而产生的吸能量损失;但对于6008-T6 结构,该差值相比不考虑应变率效应时有所增大,这表明6008-T6 铝合金的材料应变率效应会使因预处理而产生的吸能量损失。
不同材料的应变率效应对缓冲器缓冲特性的影响差别较大,在进行结构设计或期望通过预处理方式使结构变形更加可控时,需要考虑材料的应变率效应产生的影响。
图12 不同冲击速度下未预处理与预处理缓冲器的载荷−位移曲线Fig.12 Load−displacement curves of unpretreated and pretreated buffers under different impact velocities
表5 缓冲器吸能特性评估Table 5 Assessments of energy absorption of buffer
1)材料的应变率效应会显著改变缓冲器的缓冲特性,因此,在动态冲击条件下,必须考虑缓冲器受应变率效应的影响。
2)材料应变率效应对缓冲力历程无明显影响,但在一定冲击速度范围内,5083P-O铝合金材料的应变率弱化效应使缓冲力整体减小,而6008-T6铝合金材料的应变率强化效应则使缓冲力整体放大,不同冲击速度下缓冲力的缩放趋势与2种材料流动应力随应变率变化的趋势一致。
3)随着冲击速度的提高,5083P-O铝合金先表现出应变率弱化效应,使缓冲器整体鼓胀变形增大,而后又呈应变率强化效应,提高了缓冲器抵抗变形的能力。6008-T6铝合金材料的流动应力随应变率单调增加呈应变率强化效应,随着冲击速度的提高缓冲器的抗冲击能力逐渐增强。
4)在动态冲击载荷下,对缓冲器进行预处理可有效降低冲击过程中的载荷波动,但经预处理的缓冲器吸能量减少,其理想吸能效率降低。在一定冲击速度范围内,5083P-O铝合金的材料应变率效应会使缓冲器因预处理而产生的吸能量损失减少,而6008-T6铝合金的材料应变率效应则会使吸能量损失增加。对缓冲器进行结构优化设计时,需要考虑材料应变率效应的影响。