Elastic compensation for dynamic rigid-plastic solutions of structures
-
摘要: 近年来,我国学者以膜力因子法和饱和分析方法相结合为理论工具,对梁、板等结构件在脉冲载荷作用下的塑性大变形行为作了全面深入的研究,为脉冲加载下结构的最终挠度提供了优于历史上各种刚塑性近似解的最佳刚塑性预测公式。然而,由于实际工程应用中金属结构弹塑性动力响应的复杂性和数值模拟的局限性,与考虑材料弹性效应的结果相比,刚塑性解对脉冲加载下结构所预测的最终挠度的误差有多大,是一个亟待解决的关键问题。对这个问题的首阶段研究成果厘清了材料弹性对脉冲加载下结构塑性动态大变形的影响,定量评估了由最佳刚塑性理论解与弹塑性数值模拟得到的最终挠度预测结果之间的差异。在此基础上,提出了补偿弹性效应的策略和方法,即:在已有的最佳刚塑性解预测的挠度基础上添加一个补偿项,将补偿项表达为脉冲载荷强度的效应与结构自身刚度的效应分离的变量函数,并尽量减少待定系数/指数的数量,以求表达式的简洁;根据这些原则在金属结构的主要工程应用领域内选定结构刚度和外载参数的变化范围,对固支梁和固支方板的案例实施拟合与补偿,最后得到了对梁和板增添补偿项后的简单而实用的最终挠度预测公式,其相对误差在3%的范围之内,很适合工程设计应用。文末列表给出了符号与公式的一览,并对梁和方板的结果作了综合和比较。Abstract: In recent years, by combining membrane factor method (MFM) and saturation analysis (SA) as a powerful theoretical tool, scholars in China have made comprehensive studies on the large dynamic plastic deformation of beams, plates and other structures under pulse loading, leading to the best rigid-plastic predictions for the final deflection of the pulse-loaded structures, which are superior to the previously proposed various approximate rigid-plastic solutions. However, due to the complexity of the dynamic elastic-plastic response of structures used in engineering and the limitations of numerical simulations, it is critical to clarify how large is the error generated from the rigid-plastic solutions in predicting the final deflection of pulse-loaded structures compared to the result that takes the elastic effect of material into consideration. Our preliminary study on this issue, which has been published in leading international journals, reveals the effect of material’s elasticity on the large dynamic plastic deformation of structures under pulse loading, and quantitatively evaluates the discrepancies between the final deflection predicted by the best theoretical rigid-plastic solutions and that extracted from elastic-plastic numerical simulations. On this basis, the present paper proposes a strategy to compensate for elastic effect; that is, (1) adding a compensation term to the final deflection predicted by the existing best rigid-plastic solution; (2) expressing the compensation term as an elemental function of variables separation to respectively represent the effects of the pulse intensity and the structural stiffness; and (3) adopting minimum number of undetermined coefficients (or power) in the fitting function to achieve concise formulae. Meanwhile, the variation ranges of structural stiffness and dimensionless load parameter are investigated with reference to metallic structures used in their main application fields. Finally, by implementing the fitting and compensation for the cases of fully-clamped beams and square plates, simple and engineer-friendly formulae for predicting the final deflection of beams and plates are eventually obtained. With the compensation terms being added, the rigid-plastic-solution-based predictions on the final deflection of beams and plates possess a relative error within the range of 3%, which are appropriate and suitable for applications in the engineering design stage. A table at the end of the paper summarizes the major notations and formulae, as well as the comparison between the results on beams and square plates.
-
高可靠性和高安全性是现代武器弹药的两项重要需求,传爆药在爆炸序列中起着能量放大和传递作用,对于武器弹药爆炸序列的可靠动作和异常环境下的弹药安全性至关重要[1]。兼具高能量和低感度特性的传爆药配方及其性能研究,成为了含能材料的研究热点[2-8]。
TATB/RDX传爆药不仅保留了RDX的高能量属性,又融入了TATB的力、热低感度特性。Qu等[9]的研究表明,TATB/RDX传爆药的撞击感度较RDX有大幅下降,同时热分解温度提升,安全性能大幅提升。以TATB/RDX传爆药为代表的美国PBXN-7传爆药和英国Rowanex3601传爆药在钝感弹药中得到广泛应用。TATB/RDX传爆药的低感度得益于其中的TATB成分。TATB为钝感炸药,感度较低,但其反应区较宽,易受环境温度的影响。已有研究结果表明,随着环境温度的降低,TATB基炸药的爆轰波波速提高;波阵面中心与边界之间的时间差增大,波阵面曲率半径减小[10-13];爆轰传播拐角过程中的“死区”区域明显变大[14-15],起爆难度增加[16-19]。PBXN-7传爆药的曲率效应实验结果[20]表明,TATB成分的加入使得PBXN-7传爆药的爆轰传播性能受温度的影响较大。宽温域环境在武器弹药全天候作战剖面中是难以避免的,对TATB/RDX传爆药在不同温度环境中的起爆及驱动性能开展研究,对武器弹药的动作可靠性评估、爆炸序列的优化设计及钝感复合炸药基础爆轰性能的认识有着重要的意义。研究武器爆炸序列中待测炸药的宽温域冲击起爆及驱动性能,对实验设计及测试技术均有着更高的要求,尤其是低温易凝霜环境,不利于传统光电测试技术的开展,用于研究炸药基础爆轰性能的实验技术也无法直接应用于爆炸序列的性能研究中。
本文中,采用不凝霜宽温域环境加载技术、1 550 nm激光多普勒测速(photon Doppler velocimetry,PDV)技术[21-22]和无损穿透瞬态太赫兹波多普勒干涉测速(terahertz-wave Doppler interferometric velocimetry,TDV)技术[23-27],对TATB/RDX传爆药在–45~70 ℃温度环境中的爆轰增长、传播及驱动性能开展实验研究,探讨宽温域环境对TATB/RDX传爆药的反应区宽度、到爆轰距离、爆轰波速及驱动性能的影响。
1. 实 验
1.1 实验装置
本文中包含3种实验装置和2类测试技术,分别如图1~3所示。图1为炸药/窗口界面粒子速度测量实验装置。采用雷管起爆传爆药,冲击波经惰性隔层起爆待测炸药。在待测炸药末端中心点设置测试窗口,采用PDV测试技术对炸药/窗口的界面粒子速度进行测量,用以研究待测炸药的爆轰反应区剖面[21-22]。窗口材料为LiF单晶,与炸药的接触面镀0.7 μm厚的反射铝膜。实验所用PDV测速探头的输出光斑直径小于0.3 mm,激光波长为1 550 nm。实验所用传爆药为JH-9005(RDX与黏结剂的质量比为97∶3),尺寸为
∅ 50 mm×12 mm,密度为1.644 g/cm3。主炸药(TATB、RDX与黏结剂的质量比为 60∶35∶5)的尺寸为∅ 150 mm×37 mm,密度为1.816 g/cm3。传爆药与待测炸药之间的惰性隔层材料为不锈钢,厚度为2 mm。图2~3实验装置中的装药序列同图1。图2为爆轰增长及传播测量实验装置,采用TDV测试技术[23]对炸药中心爆轰波波阵面的传播历程进行测量,用以研究待测炸药冲击转爆轰过程及爆轰传播速度。太赫兹波是一种振荡频率介于微波与远红外波之间的电磁波,具有对多数非极性物质(包括常用固体炸药)穿透性好的特点[23-25]。TDV技术采用的太赫兹波能够有效穿透炸药,并对其内部的爆轰波或冲击波面进行非侵入式测量,从而获取爆轰波或冲击波的速度变化历程[26-27]。实验所采用的太赫兹波频率为0.214 67 THz,太赫兹波光束入射至炸药内的光斑尺寸约为6 mm(以半高宽计)。
图3为炸药驱动大板实验装置,待测炸药驱动下端面尺寸为
∅ 150 mm×2 mm的紫铜板,采用多点PDV测试技术对大板飞片各位置的飞行速度进行测量,用以研究待测炸药的散心及滑移爆轰驱动性能。以紫铜板中心为起点,沿待测炸药径向,每间隔14 mm布置PDV探头,共计6个,最边缘测点所在位置距离紫铜板边缘为5 mm。1.2 环境条件
实验中采用分体式高低温风冷加载系统对炸药件进行环境温度的加载,加载系统如图4所示。图4(a)为冷源设备,负责冷热空气的产生及程序控温,图4(b)为温度实验箱,内含待测炸药件及测试所需窗口。为了保障炸药件所处环境温度的一致性及测试系统的正常工作,在实验过程中,温度实验箱的位置不再移动。在–45~70 ℃范围内,选取–45、20和70 ℃作为本实验的环境温度点。环境温度的加载曲线如图5所示,保温时长不低于180 min。
2. 结果与分析
2.1 环境温度对TATB/RDX传爆药起传爆性能的影响
采用TDV测试技术获取的TATB/RDX传爆药爆轰增长及传播的典型信号如图6所示。TDV信号的周期反映了所测波阵面的速度,周期越长,所测波阵面速度越低;信号幅值与所测波阵面反射太赫兹波的强度以及穿透深度有关。在炸药爆轰增长、传播及产物膨胀的不同阶段,所测冲击波/爆轰波波阵面的传播速度及其反射太赫兹波的强度均不同,因此在TDV信号中可明确分辨炸药内的爆轰增长(S1)、爆轰传播(S2)及爆轰产物的膨胀过程(S3)。以S2爆轰传播与S3爆轰产物膨胀过程的分界为起点,对TDV所测信号的周期进行统计,结合所用太赫兹波在所测炸药内传播的折射率系数[23],给出所测炸药的到爆轰距离(run distance to detonation,RDTD);对TDV所测信号的振荡频率进行分析[23],给出所测波阵面的传播速度。
采用TDV技术获取的–45~70 ℃温度环境中TATB/RDX传爆药起传爆过程中的速度如图7所示。冲击波在炸药内经过不同距离的增长过程后转为爆轰状态,爆轰波波速基本稳定。不同温度中TATB/RDX传爆药的到爆轰距离及爆轰波波速如图8~9所示。可见,随着环境温度的降低,TATB/RDX传爆药的到爆轰距离呈线性增加,爆轰波波速逐渐提高。这可能是因为低温环境造成炸药的点火阈值提高,起爆难度变大。环境温度降低引起炸药密度增大,从而使得炸药能量释放的时间变长,反应区宽度增大(图10),反应区内的平衡时间变长,到爆轰距离增大。且密度的增大提高了炸药的能量密度,反应区内驱动爆轰波波阵面传播的能量增多,爆轰波传播速度变快。该趋势与TATB基钝感炸药曲率效应实验结果[11, 20]一致,但RDX成分的加入使得TATB/RDX传爆药爆轰波波速随温度的降低呈现非线性增长趋势。在–45~70 ℃温度环境中,TATB/RDX传爆药在2 mm不锈钢隔层起爆条件下的到爆轰距离随环境温度的降低呈线性增大趋势,线性系数为0.015 mm/℃;TATB/RDX传爆药的爆轰反应区时间宽度随环境温度的降低而线性增大,线性系数为0.165 ns/℃,该趋势与TATB基钝感炸药PBX-9502反应区宽度的预测模型[13]相一致。
2.2 环境温度对TATB/RDX传爆药驱动性能的影响
采用PDV测试技术测量的TATB/RDX传爆药驱动2 mm紫铜飞片的典型速度曲线如图11所示。从图11可以看出,飞片在飞行过程中具有明显的多阶段特征,其中,阶段a为飞片外表面的层裂层在冲击波来回反射作用下自由飞行的过程;阶段b为飞片的主体部分追赶上表面层裂层的过程;阶段c为飞片主体部分与层裂层的融合过程;阶段d则为飞片主体与层裂层融合完成后,在爆轰产物驱动下整体飞行的过程。
在70 mm半径范围内,飞片表面各处速度曲线的对比如图12所示。从图12可以看出,不同半径位置的飞片速度曲线剖面具有相似性,但随着半径的增大,飞片在测量方向的速度存在下降趋势,层裂层的飞行持续时间缩短,在测点r=42 mm处的速度变化最明显,层裂层自由飞行的持续时间变化最突出(见图13),这可能是因为在r=42 mm附近,爆轰驱动模式由散心爆轰驱动转变成滑移爆轰驱动。以惰性隔层/待测炸药的接触面中心为原点,以起爆方向为基线,计算得到散心爆轰驱动向滑移爆轰驱动的转变角度(见图14)约为49°。
在–45~70 ℃温度环境中,TATB/RDX传爆药驱动大板飞片各位置处的飞行速度曲线的对比如图15所示。飞片各点的起跳速度及层裂层的飞行速度与温度有相同的变化趋势,即:常温速度最高,高温次之,低温最低。这可能是因为飞片主体/层裂层融合前(阶段a)的驱动性能受炸药爆轰反应区特性的影响较大。对本实验中(图1)获取的界面粒子速度(图16)进行处理[28],计算出不同环境温度中待测炸药的CJ点爆压pCJ如图17所示。可以看出,在–45~70 ℃温度环境中,TATB/RDX传爆药的常温pCJ最高,高温次之,低温最低。该结果支持了上述“飞片主体/层裂层融合前(阶段a)的驱动性能受炸药爆轰反应区特性的影响较大”的猜想。
在散心爆轰驱动区(r<42 mm),大板飞片层裂层的速度及其飞行持续时间随着环境温度的变化趋势一致,但在滑移爆轰驱动区(r>42 mm),环境温度对层裂层飞行持续时间的影响并不明显。这可能与2种爆轰驱动方式下飞片主体的速度随环境温度的变化特性不同有关。
大板飞片主体和层裂层融合完成后的加速飞行速度(阶段d)受低温环境的影响较大,除r=70 mm处,低温环境中飞片的d阶段飞行速度均高于常温和高温结果,而在常温和高温环境中结果基本相同。这可能是由于低温环境中爆轰产物的特性与高温及常温结果差异较大造成的。
3. 结 论
采用PDV及TDV测试技术对TATB/RDX传爆药在–45~70 ℃温度环境中的起传爆性能及驱动性能进行了实验研究,结果表明:TATB/RDX传爆药的到爆轰距离由高温70 ℃的1.3 mm增大到低温–45 ℃的3.1 mm,且近乎呈线性增大趋势,线性系数为0.015 mm/℃。TATB/RDX传爆药的爆轰传播速度随环境温度的降低呈非线性增大趋势;TATB/RDX传爆药的驱动性能随温度的变化趋势在飞片的主体-层裂层融合前后存在明显差异。
-
图 8 补偿前(实线)和按式(8)作补偿后(虚线)与弹塑性有限元模拟得到的固支梁最终挠度之间的相对偏差
Figure 8. The relative discrepancies in the final deflection of fully-clamped beams between rigid-plastic predictions (solid lines displaying before compensation, and broken lines displaying after compensation given by Eq.(8)) and elastic-plastic simulation results
表 1 低碳钢和铝合金的材料性质
Table 1. Material properties of mild steel and aluminum alloys
材料 杨氏模量E/GPa 泊松比μ 屈服应力Y/MPa 低碳钢Q235 210 0.3 235 铝合金6061 71 0.3 240 铝合金7075 71 0.3 505 表 2 金属板的半边长与板厚之比 的常用范围
Table 2. Commonly-used ranges of half length-to-plate thickness ratio for metallic plates
工程领域 a/h 船舶海洋结构 10~125 汽车与运载机械 150~900 航空航天结构 80~1500 表 3 主要结果整理及固支梁与固支方板的比较
Table 3. Summary of main results as well as comparison of fully-clamped beams and plates
项目 固支梁(矩形截面) 固支方板 几何参数 长度2L,宽度b,厚度h 边长2a,厚度h 材料参数 杨氏模量E,屈服应力Y 杨氏模量E,泊松比μ,屈服应力Y 塑性极限弯矩 Mp = Ybh2/4 M0 = Yh2/4 准静态坍塌压力 线载荷pY = 4Mp/(bL) = Yh2/L 面载荷pY = 12M0/a2 = 3Yh2/a2 无量纲压力 λ ≡ p0/pY 结构刚度ζ ζbeam≡hL√EY ζplate≡ha√EY(1−μ) 最佳刚塑性解预测的最终挠度 ηrpf≡wrp0fh=1.10λ−0.414 ηrpf≡wrp0fh=2.16λ−1.456 刚塑性解加上弹性补偿后预测的最终挠度 η*f=1.10λ−0.414+(−0.08λ+0.29)ζ−0.9 η*f=2.160λ−1.456+(0.03λ2−0.33λ+0.85)ζ−0.4 弹塑性模拟得到的最终挠度 ηepf≡wep0fh 上述二者之间的相对偏差 D*η=η*f−ηepfηepf≤3% -
[1] JOHNSON W. Impact strength of materials [M]. London: Edward Arnord, 1983. [2] JONES N. Structural impact [M]. 2nd ed. Cambridge: Cambridge University Press, 2011. [3] STRONGE W J, YU T X. Dynamic models for structural plasticity [M]. Berlin: Springer-Verlag, 2012. [4] 余同希, 邱信明. 冲击动力学 [M]. 北京: 清华大学出版社, 2011. [5] SYMONDS P S, FRYE C W G. On the relation between rigid-plastic and elastic-plastic predictions of response to pulse loading [J]. International Journal of Impact Engineering, 1988, 7(2): 139–149. DOI: 10.1016/0734-743X(88)90022-X. [6] YU T X. Elastic effect in the dynamic plastic response of structures [M]//JONES N, WIERZBICKI T. Structural Crashworthiness and Failure. London: Springer-Verlage, 1993: 341–384. DOI: 10.1201/9781482262544. [7] 余同希, 朱凌, 许骏. 结构冲击动力学进展(2010−2020) [J]. 爆炸与冲击, 2021, 41(12): 121401. DOI: 10.11883/bzycj-2021-0113.YU T X, ZHU L, XU J. Progress in structural impact dynamics during 2010−2020 [J]. Explosion and Shock Waves, 2021, 41(12): 121401. DOI: 10.11883/bzycj-2021-0113. [8] 余同希, 朱凌, 陈发良. 饱和冲量与膜力因子法: 强动载荷下结构塑性大变形的分析和预测方法 [M]//陈建康, 白树林. 材料的非线性力学性能研究进展. 北京: 机械工业出版社, 2021: 17–29. [9] 余同希, 田岚仁, 朱凌. 强脉冲载荷作用下结构塑性大变形的最大挠度直接预测 [J]. 力学学报, 2023, 55(5): 1113–1123. DOI: 10.6052/0459-1879-22-607.YU T X, TIAN L R, ZHU L. Direct prediction of maximum deflection for plastically deformed structures under intense dynamic pulse [J]. Chinese Journal of Theoretical and Applied Mechanics, 2023, 55(5): 1113–1123. DOI: 10.6052/0459-1879-22-607. [10] 朱凌, 田岚仁, 李德聪, 等. 饱和冲量及其等效方法在舱室内爆炸中的应用 [J]. 中国舰船研究, 2021, 16(2): 99–107. DOI: 10.19693/j.issn.1673-3185.01876.ZHU L, TIAN L R, LI D C, et al. Saturated impulse and application of saturation equivalent method in cabin explosion [J]. Chinese Journal of Ship Research, 2021, 16(2): 99–107. DOI: 10.19693/j.issn.1673-3185.01876. [11] TIAN L R, CHEN F L, ZHU L, et al. Saturated analysis of pulse-loaded beams based on Membrane Factor Method [J]. International Journal of Impact Engineering, 2019, 131: 17–26. DOI: 10.1016/j.ijimpeng.2019.04.021. [12] TIAN L R, CHEN F L, ZHU L, et al. Large deformation of square plates under pulse loading by the combination of Saturated Analysis and Membrane Factor Method [J]. International Journal of Impact Engineering, 2020, 140: 103546. DOI: 10.1016/j.ijimpeng.2020.103546. [13] HU Q J, ZHU L, YU T X. Elastic effects on the dynamic plastic deflection of pulse-loaded beams [J]. International Journal of Impact Engineering, 2023, 176: 104550. DOI: 10.1016/j.ijimpeng.2023.104550. [14] HU Q J, ZHU L, YU T X. Elastic effect on the final deflection of rigid-plastic square plates under pulse loading [J]. Thin-Walled Structures, 2023, 193: 111238. DOI: 10.1016/j.tws.2023.111238. [15] YOUNGDAHL C K. Correlation parameters for eliminating the effect of pulse shape on dynamic plastic deformation [J]. Journal of Applied Mechanics, 1970, 37(3): 744–752. DOI: 10.1115/1.3408605. [16] YOUNGDAHL C K. Influence of pulse shape on the final plastic deformation of a circular plate [J]. International Journal of Solids and Structures, 1971, 7(9): 1127–1142. DOI: 10.1016/0020-7683(71)90057-6. [17] ZHU L, TIAN L R, CHEN F L, et al. A new equivalent method for complex-shaped pulse loading based on saturation analysis and membrane factor method [J]. International Journal of Impact Engineering, 2021, 158: 104018. DOI: 10.1016/j.ijimpeng.2021.104018. [18] CHEN F L, TIAN L R, YU T X, et al. Complete solution of large plastic deformation of square plates under exponentially decaying pulse loading [J]. Acta Mechanica Solida Sinica, 2021, 34(6): 922–936. DOI: 10.1007/s10338-021-00280-6. 期刊类型引用(0)
其他类型引用(1)
-