Experimental Research Progress on Physical Properties and “Phase Transition” of Polymers under Impact Loading
-
摘要: 聚合物被广泛应用于国防和国民经济的各个领域,在服役过程中,不可避免地暴露于高温高压等极端环境,因此,有必要研究其在冲击荷载下的物性及“相变”问题。聚合物具有独特的分子链结构,表现出异于金属等大多数材料的性能。聚合物的雨贡纽曲线低压外推的截距明显高于零压体波声速,且低压下的波剖面呈现带弧形的波系结构。在20~30 GPa压力范围,其雨贡纽曲线存在明显的转折,说明材料在冲击下发生了“相变”。首先,对该“相变”反映的化学分解和晶格结构转变进行了解释,并对其中蕴含的“相变”动力学问题进行了研究。然后,简单介绍了基于化学分解的物态方程的建模方法。最后,对聚合物在冲击荷载下的物性和“相变”研究提出了展望。Abstract: Polymers are widely used in various fields of national defense and national economy. They are inevitably exposed to extreme conditions of high temperature and high pressure during applying. Thus, it is necessary to study their physical properties and “phase transition” under impact loading. Because of their characteristic molecular chain structure, polymers show different properties from most materials such as metals. The intercept extrapolated from Hugoniot curve at low pressure is obviously higher than their body sound velocity at atmospheric pressure. The wave profile at low pressure presents a structure with an arc shape. At 20–30 GPa, the Hugoniot line turns obviously, indicating that the material has undergone a “phase transition” under impact loading. The “phase transition” is explained as chemical decomposition or lattice structure transformation, and the kinetics of “phase transition” is studied. In addition, the modeling method of equation of state based on chemical decomposition is briefly introduced. Finally, the prospect is put forward according to the doubtful points in the study of physical properties and “phase transition” of polymers under impact loading.
-
Key words:
- polymer /
- Hugoniot /
- equation of state /
- chemical decomposition /
- lattice structure transformation
-
半穿甲导弹是在穿甲弹基础上发展起来的用于攻击舰船目标的武器,兼具穿甲和爆破毁伤能力,是舰船防护领域的重点设防对象[1]。前舱物作为半穿甲导弹的重要组成部分,在导弹的制导和控制方面起着至关重要的作用。半穿甲导弹侵彻靶板的过程中,前舱物先于战斗部接触靶板,对靶板产生影响。因此,研究前舱物对战斗部侵彻靶板的影响具有实际意义。迄今为止,对于半穿甲导弹侵彻舰船结构等问题,研究人员通常忽略前舱物的影响,即只考虑战斗部的动能穿甲作用,将其抽象为低速大质量弹体的穿甲问题。陈斌等[2]利用忽略前舱物的半穿甲模拟弹,开展了侵彻陶瓷/钢复合装甲试验,明确了弹着角对半穿甲弹侵彻效能有明显影响;朱锡等[3]采用与实际半穿甲导弹战斗部几何形状相似的弹丸,开展了侵彻多种装甲结构试验,得到了不同装甲结构形式抗半穿甲导弹动能穿甲的特性。目前,针对前舱物对半穿甲导弹侵彻靶板影响的研究相对较少,这是由于前舱物的组成部件复杂,无法完全模拟所有器件,因此,通常对前舱物进行力学等效。
对前舱物的等效主要有两种方法。第一种方法为采用质量块对前舱物进行等效。楼建锋等[4]通过有限元方法,分别使用金属材料和脆性材料质量块等效模拟前舱物,得到了高速(675 m/s)和低速(250 m/s)两种速度条件下前舱物对战斗部侵彻过程的影响。陈刚[5]利用质量块模拟前舱物,对比分析了前舱物对战斗部高速撞击靶板的影响,得到了前舱物改进战斗部弹体在斜穿甲过程中的荷载环境条件。Chen等[6]考虑前舱物作用导致的靶板结构预响应,将前舱物剩余部分等效为附着在战斗部前部的质量块,并对剪切冲塞模型进行修正,得到了与实验结果相吻合的终点弹道公式。该等效方法的关键是质量块材料的选取,只有选择与前舱物名义力学特性相近的材料,才能正确表征前舱物对导弹弹体动能穿甲的影响。第二种方法为采用薄壁圆筒结构对前舱物进行等效。徐钰巍等[7]采用硬铝圆筒结构对前舱物进行等效,建立了变截面刚度/含缺陷的自由梁模型,并结合模型试验得到了前舱物在弹体侵彻靶板初期对弹体姿态变化的影响规律。该等效方法忽略了前舱物内部仪器设备强度的影响,仅考虑前舱段外壳在侵彻过程中所起的作用,在研究半穿甲导弹正侵彻时,与实际情况的差异较大。对于正侵彻而言,将前舱物等效为质量块的方法是可行的,但现有的等效方法为了体现前舱物稀松多孔、中空易碎的结构特点,所选取的等效材料模型和前舱物在撞击过程中所表现出的力学特性存在不相符的情况。
本研究采用轻质泡沫铝材料模拟半穿甲导弹前舱物,建立变截面刚度弹体侵彻金属薄板的有限元模型,分析带前舱物平头弹侵彻金属薄板的动态力学过程,研究不同靶板厚度、弹速条件下前舱物对平头战斗部侵彻金属薄板的影响。
1. 计算模型与材料参数
1.1 前舱物特性
前舱物主要由雷达导引头、导航系统、电源系统、高度表等仪器设备组成,其特点是质量较大(约为战斗部总质量的1/4~1/3)且分布均匀,整体结构中空、易碎、延展性差,结构强度远低于战斗部。前舱物由高分子聚合物、铝合金、单晶硅等多种材料构成,平均密度低(小于1 g/cm3)。前舱物在撞击靶板时,会经历“接触—屈服—压实—破碎飞散”的破坏过程,其中屈服阶段持续时间长,名义应力基本保持不变,即存在较长的平台应力,随后前舱物被压实,名义应力大幅增加,前舱段整体长度大幅减小,这与泡沫铝材料在冲击荷载作用下的位移变化曲线类似[8-10],如图1所示。因此,本研究中前舱物材料模型采用低强度、低密度、有长应力平台的泡沫铝材料进行等效。
1.2 带前舱物战斗部等效计算模型
半穿甲反舰导弹通常由前舱物、战斗部和推进系统3部分组成。战斗部是弹体中强度最高、刚度最大的部分,有较大的药室,装填炸药量较多,大多为钝头[1]。本研究建立的前舱物战斗部等效计算模型如图2所示。战斗部壳体模型外径为428.40 mm,内径为384.79 mm,内置装药模型直径为384.79 mm,前舱物等效为与战斗部截面尺寸相等的圆柱体。
采用有限元软件LS-DYNA模拟不同工况下带前舱物战斗部侵彻靶板的过程。计算中,结构采用拉格朗日实体单元进行网格划分,时间积分采用显式中心差分法。靶板平面尺寸为3 000 mm × 3 000 mm,边界条件设置为固支约束。计算模型设置为沿中轴面对称,建立1/2模型,战斗部与前舱物采用六面体网格划分。靶板自中心采用过渡网格,中心区域网格尺寸为2.5 mm × 2.5 mm,距离撞击区域越远,网格尺寸越大。
1.3 本构方程与材料参数
数值仿真中主要涉及两种材料模型。战斗部壳体和金属靶板采用金属材料常用的Johnson-Cook材料模型,战斗部内置装药和前舱物等效材料采用弹塑性材料常用的弹塑性硬化模型。
Johnson-Cook材料模型考虑了塑性强化、应变率强化、温升软化效应,能够反映材料在高温、高应变率下力学特性的变化,其具体形式为
s=(A+Bˉepn)(1+Cln˙ε˙ε0)[1−(T−T0Tm−T0)m] (1) 式中:s为应力;A为静态屈服强度;B、n通过拟合等效应力-应变曲线数据获得;
ˉep 为有效塑性应变;C为应变率系数;˙ε 为有效塑性应变率;˙ε0 为参考塑性应变率,一般取˙ε0 = 1 s–1;T0为参考温度(取室温);Tm为材料熔点;m为热软化指数。战斗部壳体和金属靶板的材料参数根据文献[11-12]并结合工程实践调整得到,如表1所示,其中cp为定压比热容,ρ为密度,μ为泊松比。
表 1 弹体和靶板的J-C模型参表Table 1. J-C model parameters of projectile body and targetMaterial ρ/(g·cm–3) E/GPa μ cp/(J·kg–1·K–1) T0/K Tm/K Warhead shell 7.80 205 0.28 400.90 300 1 765 Target plate 7.80 200 0.30 452.00 298 1 881 Material ˙ε/s–1 A/MPa B/MPa n C m Warhead shell 1 760 500 0.53 0.014 1.13 Target plate 1 355 450 0.36 0.022 1.00 弹塑性硬化模型考虑了材料应变的等向或随动塑性强化。其应变率效应由Cowper-Symonds模型描述
σy=[1+(˙εc)1/p](σ0+βEEtE−Etεeff) 式中:σy为屈服强度,
˙ε 为等效塑性应变率;c、p、β为材料常数;σ0为静态屈服强度,E为杨氏模量,Et为硬化模量,εeff为有效塑性应变。内置装药和前舱物等效泡沫材料的模型参数如表2所示,其中fs为失效应变。考虑到不同导弹的前舱物内包含仪器设备的数量和种类不同,前舱物等效材料的屈服强度无法准确给出,本研究只给出取值范围。表 2 弹塑性硬化模型材料参数Table 2. Material parameters of the Elastoplastic hardening modelMaterial ρ/(kg·m–3) E/GPa μ σy/MPa Et/MPa β c/s–1 p fs Charge 1 750 5 0.3 20 300 0 0 0 5 Nose cabin 700 0.5 0.3 0.5–10.0 0 0 0 0 0.4 2. 结果分析
2.1 带前舱物战斗部侵彻过程分析
图3为带前舱物战斗部侵彻金属薄板的过程。图4为侵彻过程中典型时刻的靶板内应力分布。典型工况下,带前舱物战斗部与不计前舱物战斗部侵彻金属薄板的速度时程曲线如图5所示。根据金属薄板受力和变形的特点,可以将带前舱物战斗部侵彻金属薄板的过程分为3个阶段。
(1)金属薄板结构预变形阶段。从前舱物接触靶板开始,接触区域内最大压强远大于前舱物等效材料的极限强度和金属薄板的屈服强度。因此,前舱物发生变形破坏,金属薄板产生大范围碟形塑性变形,并获得一定的横向速度,战斗部与靶板之间的速度差减小,弹靶接触区边缘的应力集中程度减轻。由图4可知,在金属薄板结构预变形阶段,前舱物与靶板接触区域内正应力较大,且随着与中心距离的增大而降低;剪应力很小,且只存在于接触区域及其附近区域。在此阶段,战斗部的速度变化可以分为两个时期,如图5所示。AB时期,由于前舱物与金属薄板撞击产生的压缩应力波尚未传递到战斗部,战斗部以初速度向前运动;BC时期,随着前舱物与金属薄板撞击产生的压缩应力作用于战斗部,战斗部的速度开始下降,但由于作用于战斗部的压缩应力强度较小,战斗部速度下降较慢。
(2)冲塞破坏阶段。从战斗部接触金属薄板开始,接触区域由于速度差产生剪切变形,在拉应力与剪应力的联合作用下,接触区域周围产生裂纹,裂纹直径略大于战斗部直径。如图4所示,此阶段战斗部与靶板接触区域内正应力最大,且由于裂纹拓展释放压缩波,接触区域内出现拉伸应力与压缩应力混合的情况;接触区域内剪应力强度远大于结构预变形阶段,剪应力在接触边缘附近达到峰值,并随着与中心距离的增大迅速衰减。如图5所示,在此阶段,战斗部直接侵彻金属薄板,弹体受到的应力强度大,速度下降明显,可近似视为匀减速运动。
(3)弹体贯穿阶段。随着弹体的进一步运动,裂纹贯穿整个金属薄板厚度方向,形成与战斗部直径相似的帽形失效块,同时,金属薄板由于惯性作用继续发生碟形变形。如图4所示,该阶段由于战斗部与靶板接触的区域形成失效块与靶板脱离,接触区域内正应力与剪应力均为零。靶板动能转化为变形能,挠度增加,靶板仍有较大正应力,由于接触区域附近靶体材料存在速度差,因此靶板接触区域附近存在一定的剪应力。如图5所示,在该阶段,战斗部的速度变化呈现两个时期:DE阶段,由于冲塞块附着在战斗部头部,战斗部推动冲塞块前进,战斗部的一部分动能转化为冲塞块的动能,导致战斗部速度略有下降;EF阶段,冲塞块与战斗部速度达到一致,战斗部向前匀速运动。
由图5中带前舱物战斗部与不计前舱物战斗部的速度时程曲线对比可知,不计前舱物战斗部侵彻靶板的贯穿时间为90 μs,带前舱物战斗部的侵彻贯穿时间为2 180 μs,侵彻作用时间大幅增加,但增加的侵彻时间主要是前舱物被压溃破坏的时间。在金属薄板预变形阶段,金属薄板与前舱物接触区域获得的速度较小,与战斗部之间的速度差较大,因此,战斗部侵彻贯穿金属薄板所需要的时间并未明显增加。此外,相比于不计前舱物战斗部速度变化情况,前舱物的存在使得战斗部在侵彻过程中的受力更加复杂,战斗部速度变化的阶段性和多样性更加明显。
从塑性变形范围来看,通过有限元仿真软件的计算结果可知,不计前舱物战斗部侵彻靶板过程[13] 中金属薄板塑性变形的直径为883.8 mm,带前舱物战斗部侵彻靶板过程中,由于存在金属薄板预变形阶段,金属薄板迎弹面塑性变形范围大幅增加,直径达到2 629.7 mm。其原因在于:带前舱物战斗部侵彻金属薄板的过靶时间长,弹靶接触区域产生的塑性波传播范围远大于不计前舱物战斗部。因此,前舱物的“缓冲垫”作用明显,在大幅延长战斗部过靶时间的同时大幅扩大了金属薄板塑性变形范围。
2.2 前舱物等效模型屈服强度的影响
由于不同导弹的前舱物内包含仪器设备的种类和数量不同,因此,无法定量给出前舱物等效材料的屈服强度。本研究给定前舱物等效材料的屈服强度范围为0.5~10.0 MPa,在此屈服强度范围内研究前舱物对平头战斗部侵彻金属薄板的影响。
图6给出了初速度为340 m/s的不同前舱物等效屈服强度战斗部侵彻20 mm厚金属薄板的应变云图。由图6可知,前舱物等效屈服强度在0.5~10.0 MPa范围内时,金属薄板的塑性应变分布相似,战斗部贯穿金属薄板形成的帽形失效块直径相等。表3列出了不同前舱物等效屈服强度战斗部侵彻靶板时金属薄板的塑性变形范围和战斗部剩余速度。由表3可知,前舱段等效屈服强度增加时,金属薄板塑性变形范围并未出现明显变化,而是在一定范围内波动,最大波动幅度为0.994%。其原因在于,模型计算结果的输出时间间隔为20 μs,侵彻作用时间选取最接近战斗部贯穿靶板的时间,时间选取存在不准确的情况。随着前舱物等效屈服强度的增加,战斗部剩余速度增大,但增大幅度很小,仅为0.147%。因此,可以认为前舱物等效屈服强度对带前舱物平头弹侵彻金属薄板的影响很小,可以忽略。这是由于前舱物存在破坏强度低、延展性差的特点,变形破坏所耗散的变形能远小于带前舱物战斗部所具有的动能,因此当前舱物等效模型的屈服强度远小于靶板时,前舱物的惯性力起主要作用,结构强度的影响可以忽略不计。
表 3 不同等效屈服强度下侵彻结果的比较Table 3. Comparison of penetration results with different yield stressEquivalent yield strength/MPa Distortion range/mm Residual velocity/(m·s–1) 0.5 2 629.70 325.88 1.0 2 624.41 326.30 2.0 2 629.31 326.37 5.0 2 625.23 326.75 10.0 2 659.20 326.86 2.3 靶板厚度和侵彻速度的影响
为了分析不同靶板厚度时前舱物对战斗部侵彻的影响,分别模拟了靶板厚度在10~40 mm范围内的7种工况。不同靶板厚度条件下,战斗部剩余速度如图7所示。由图7可知:带前舱物战斗部剩余速度与不计前舱物战斗部的剩余速度均随靶板厚度的增大而减小,且变化规律一致,均近似呈线性变化;靶板厚度相同时,两者的剩余速度相差很小,最大偏差为0.67%。因此,靶板厚度在10~40 mm范围内时,前舱物对战斗部侵彻的影响很小。
常见反舰半穿甲导弹的飞行马赫数Ma 在 0.8~1.0之间[14],为了分析不同侵彻速度下前舱物对平头战斗部侵彻金属薄板的影响,分别模拟了侵彻速度在238~340 m/s范围内的7种情况。不同侵彻速度条件下,战斗部的剩余速度如图8所示。由图8可知,随着侵彻速度的增大,不计前舱物平头战斗部和带前舱物平头战斗部的剩余速度均近似呈线性增长,且带前舱物平头战斗部的剩余速度略大于不计前舱物平头战斗部(最大偏差为1.15%)。Chen等[6]提出的带前舱物钝头弹侵彻金属靶的计算公式可得到与本模型相似的结论:侵彻速度相同时,带前舱物平头战斗部的剩余速度略大于不计前舱物平头战斗部。
通过对不同靶板厚度和侵彻速度条件下战斗部的剩余速度进行分析可知,带前舱物战斗部的剩余速度均大于不计前舱物战斗部,原因在于:与不计前舱物战斗部相比,带前舱物战斗部的质量更大,侵彻速度相同时动能更大。前舱物变形破坏和金属靶板结构预变形消耗部分前舱物具有的动能,使带前舱物战斗部以更大动能侵彻靶板。带前舱物战斗部对靶板的侵彻破坏模式为冲塞破坏,与不带前舱物战斗部侵彻靶板的破坏模式相同,在侵彻过程中消耗的能量相差不大。因此,带前舱物战斗部的剩余速度略大于不计前舱物战斗部。
3. 结 论
本研究建立了变截面刚度弹体模型,对带前舱物战斗部侵彻金属薄板进行仿真分析。通过对比不同前舱物等效屈服强度战斗部的侵彻性能,以及不同靶板厚度与初速度条件下带前舱物战斗部与不计前舱物战斗部的剩余速度,得出如下主要结论。
(1)根据靶板的受力变形特征,带前舱物平头战斗部侵彻金属薄板可以分为3个阶段,即靶板预变形阶段、冲塞破坏阶段和弹体贯穿阶段。虽然带前舱物平头战斗部与不计前舱物平头战斗部侵彻靶板的过程存在明显差异,但靶板破坏模式均为剪切冲塞破坏。
(2)当前舱物等效材料的强度远小于靶板的屈服强度时,前舱物等效泡沫铝材料的屈服强度对战斗部侵彻性能的影响可以忽略不计。
(3)在导弹侵彻目标的过程中,一方面,前舱物的存在增大了弹体的总动能,且先于战斗部撞击目标,造成目标结构变形甚至破坏,对战斗部侵彻靶板起到了积极作用;另一方面,由于前舱物的强度远小于战斗部,前舱物的存在延长了弹体的过靶时间,降低了靶板过载,起到了“缓冲垫”的作用,不利于平头战斗部对目标的侵彻。在这两方面综合影响下,当飞行马赫数Ma 在 0.8~1.0之间,侵彻厚度为10~40 mm的靶板时,带前舱物平头战斗部的剩余速度比不计前舱物平头战斗部均有所增加,但增加幅度小于2%。
-
Material Before phase transition After phase transition C0a/(km·s−1) sa up,a/(km·s−1) C0b/(km·s−1) sb up,b/(km·s−1) High density polyethylene 2.86 1.57 0.7−3.2 3.27 1.43 3.4−5.3 Epoxy 2.69 1.51 0.4−2.8 2.88 1.35 3.6−5.2 Polycarbonate 2.33 1.57 0.4−2.6 2.06 1.39 3.6−5.2 Polyimide 2.66 1.48 0.6−2.2 0.93 1.64 3.4−4.6 Polysulfone 2.35 1.55 0.7−2.4 1.58 1.51 3.4−5.1 Material At zero pressure pt/GPa ρ0/(g·cm−3) Cl/(km·s−1) Ct/(km·s−1) Cb/(km·s−1) High density polyethylene 0.954 2.462 1.014 2.166 24.7 Epoxy 1.192 2.641 1.177 2.264 23.1 Polycarbonate 1.196 2.187 0.886 1.933 20.0 Polyimide 1.414 2.723 1.217 2.332 17.8 Polysulfone 1.235 2.249 0.930 1.976 18.5 Constituent rs/Å ϵk−1B/K α Constituent rs/Å ϵk−1B/K α H2 3.67 36.9 10.6 CO 3.55 155 12 H2O 3.22 195 13.4 HCOOH 4.096 335 13.781 H 1.55 110 10 C 3.7 71.4 13 CH4 4.45 152.1 11.47 O 2.41 700 11 C2H6 4.86 246.7 13 CO2 4.096 335 13.781 O2 4.11 75 13.117 -
[1] 国家自然科学基金委员会, 中国科学院. 中国学科发展战略: 软凝聚态物理学 [M]. 北京: 科学出版社, 2020: 1−25.National Natural Science Foundation of China, Chinese Academy of Sciences. Chinese discipline development strategy: soft condensed matter physics [M]. Beijing: Science Press, 2020: 1−25. [2] 朱诚身. 聚合物结构分析 [M]. 2版. 北京: 科学出版社, 2010: 1−5.ZHU C S. Analysis of polymer structure [M]. 2nd ed. Beijing: Science Press, 2010: 1−5. [3] MARSH S P. LASL shock Hugoniot data [M]. Berkeley: University of California Press, 1980. [4] CARTER W J, MARSH S P. Hugoniot equation of state of polymers: LA-13006-MS [R]. Los Alamos: Los Alamos National Laboratory, 1995. [5] BOURNE N K. On the shock response of polymers to extreme loading [J]. Journal of Dynamic Behavior of Materials, 2016, 2: 33–42. doi: 10.1007/s40870-016-0055-5 [6] MORRIS C E, FRITZ J N, MCQUEEN R G. The equation of state of polytetrafluoroethylene to 80 GPa [J]. The Journal of Chemical Physics, 1984, 80(10): 5203–5218. doi: 10.1063/1.446591 [7] HARTLEY N J, BROWN S, COWAN T E, et al. Evidence for crystalline structure in dynamically-compressed polyethylene up to 200 GPa [J]. Scientific Reports, 2019, 9(1): 4196. doi: 10.1038/s41598-019-40782-5 [8] HUBER R C, DATTELBAUM D M, LANG J M, et al. Polyimide dynamically compressed to decomposition pressures: two-wave structures captured by velocimetry and modeling [J]. Journal of Applied Physics, 2023, 133(3): 035106. doi: 10.1063/5.0128515 [9] JONES A H, ISBELL W M, MAIDEN C J. Measurement of the very-high-pressure properties of materials using a light-gas gun [J]. Journal of Applied Physics, 1966, 37(9): 3493–3499. doi: 10.1063/1.1708887 [10] 金柯, 习锋, 杨慕松, 等. 化爆加载装置系列化设计 [J]. 含能材料, 2003, 11(3): 113–115, 122. doi: 10.3969/j.issn.1006-9941.2003.03.001JIN K, XI F, YANG M S, et al. Design of serialization explosive-loading device [J]. Energetic Materials, 2003, 11(3): 113–115, 122. doi: 10.3969/j.issn.1006-9941.2003.03.001 [11] PAISLEY D L, LUO S N, GREENFIELD S R, et al. Laser-launched flyer plate and confined laser ablation for shock wave loading: validation and applications [J]. Review of Scientific Instruments, 2008, 79(2): 023902. doi: 10.1063/1.2839399 [12] 种涛, 张朝辉, 王贵林, 等. 35 GPa斜波加载下RDX单晶炸药的动力学行为 [J]. 高压物理学报, 2022, 36(3): 030103. doi: 10.11858/gywlxb.20210803CHONG T, ZHANG Z H, WANG G L, et al. Dynamic behaviors of RDX single crystal under ramp wave compression up to 35 GPa [J]. Chinese Journal of High Pressure Physics, 2022, 36(3): 030103. doi: 10.11858/gywlxb.20210803 [13] 王贵林. 磁驱动平面加载实验技术及其在高压物态方程研究中的应用 [D]. 合肥: 中国科学技术大学, 2014.WANG G L. Magnetic loading techniques and its applications in high-pressure EOS [D]. Hefei: University of Science and Technology of China, 2014. [14] ZHERNOKLETOV M V, GLUSHAK B L. Material properties under intensive dynamic loading [M]. Berlin: Springer, 2006: 74–91. [15] HAYES B. Particle-velocity gauge system for nanosecond sampling rate of shock and detonation waves [J]. Review of Scientific Instruments, 1981, 52(4): 594–603. doi: 10.1063/1.1136643 [16] BARKER L M. The development of the VISAR, and its use in shock compression science [J]. AIP Conference Proceedings, 2000, 505(1): 11–18. doi: 10.1063/1.1303413 [17] WENG J D, TAN H, WANG X, et al. Optical-fiber interferometer for velocity measurements with picosecond resolution [J]. Applied Physics Letters, 2006, 89(11): 111101. doi: 10.1063/1.2335948 [18] DOLAN D H. Accuracy and precision in photonic Doppler velocimetry [J]. Review of Scientific Instruments, 2010, 81(5): 053905. doi: 10.1063/1.3429257 [19] WU J, LI J B, LI J, et al. A sub-nanosecond pyrometer with broadband spectral channels for temperature measurement of dynamic compression experiments [J]. Measurement, 2022, 195: 111147. doi: 10.1016/j.measurement.2022.111147 [20] TAKAHARA A, HIGAKI Y, HIRAI T, et al. Application of synchrotron radiation X-ray scattering and spectroscopy to soft matter [J]. Polymers, 2020, 12(7): 1624. doi: 10.3390/polym12071624 [21] CHEN S, LI Y X, ZHANG N B, et al. Capture deformation twinning in Mg during shock compression with ultrafast synchrotron X-ray diffraction [J]. Physical Review Letters, 2019, 123(25): 255501. doi: 10.1103/PhysRevLett.123.255501 [22] WINEY J M, GUPTA Y M. Shock-induced chemical changes in neat nitromethane: use of time-resolved Raman spectroscopy [J]. The Journal of Physical Chemistry B, 1997, 101(50): 10733–10743. doi: 10.1021/jp972588a [23] RASTOGI V, CHAURASIA S, RAO U, et al. Time-resolved Raman spectroscopy of polystyrene under laser driven shock compression [J]. Journal of Raman Spectroscopy, 2017, 48(7): 1007–1012. doi: 10.1002/jrs.5166 [24] DATTELBAUM D M, COE J D. Shock-driven decomposition of polymers and polymeric foams [J]. Polymers, 2019, 11(3): 493. doi: 10.3390/polym11030493 [25] 谭华. 实验冲击波物理 [M]. 北京: 国防工业出版社, 2018: 1−26.TAN H. Experimental shock wave physics [M]. Beijing: National Defense Industry Press, 2018: 1−26. [26] MILLETT J C F, BOURNE N K, GRAY G T. The response of polyether ether ketone to one-dimensional shock loading [J]. Journal of Physics D: Applied Physics, 2004, 37(6): 942–947. doi: 10.1088/0022-3727/37/6/021 [27] ROBERTS A, APPLEBY-THOMAS G J, HAZELL P. Experimental determination of Grüneisen gamma for polyether ether ketone (PEEK) using the shock-reverberation technique [J]. AIP Conference Proceedings, 2012, 1426(1): 824–827. doi: 10.1063/1.3686405 [28] BIAN Y L, CHAI H W, YE S J, et al. Compression and spallation properties of polyethylene terephthalate under plate impact loading [J]. International Journal of Mechanical Sciences, 2021, 211: 106736. doi: 10.1016/j.ijmecsci.2021.106736 [29] BOURNE N K, MILLETT J C F. On the influence of chain morphology on the shock response of three thermoplastics [J]. Metallurgical and Materials Transactions A, 2008, 39(2): 266–271. doi: 10.1007/s11661-007-9371-7 [30] YE S J, CHAI H W, XIAO X H, et al. Spallation of polycarbonate under plate impact loading [J]. Journal of Applied Physics, 2019, 126(8): 085105. doi: 10.1063/1.5108965 [31] COE J D, BROWN E, CADY C M, et al. Equation of state and damage in polyethylene: LA-UR-17-29234 [R]. Los Alamos: Los Alamos National Laboratory, 2017. [32] DATTELBAUM D M, COE J D, RIGG P A, et al. Shockwave response of two carbon fiber-polymer composites to 50 GPa [J]. Journal of Applied Physics, 2014, 116(19): 194308. doi: 10.1063/1.4898313 [33] MORRIS C E, LOUGHRAN E D, MORTENSEN G F, et al. Shock induced dissociation of polyethylene: LA-UR-89-2864 [R]. Los Alamos: Los Alamos National Laboratory, 1989. [34] ABBOTT A, BRANCH B, BROWN E N, et al. The dynamic response of polymers interrogated by 3rd generation X-ray light source: LA-UR-19-29436 [R]. Los Alamos: Los Alamos National Laboratory, 2019. [35] KRAUS D, VORBERGER J, PAK A, et al. Formation of diamonds in laser-compressed hydrocarbons at planetary interior conditions [J]. Nature Astronomy, 2017, 1(9): 606–611. doi: 10.1038/s41550-017-0219-9 [36] KRAUS D, HARTLEY N J, FRYDRYCH S, et al. High-pressure chemistry of hydrocarbons relevant to planetary interiors and inertial confinement fusion [J]. Physics of Plasmas, 2018, 25(5): 056313. doi: 10.1063/1.5017908 [37] FONTANA L, VINH D Q, SANTORO M, et al. High-pressure crystalline polyethylene studied by X-ray diffraction and ab initio simulations [J]. Physical Review B, 2007, 75(17): 174112. doi: 10.1103/PhysRevB.75.174112 [38] FONTANA L, SANTORO M, BINI R, et al. High-pressure vibrational properties of polyethylene [J]. The Journal of Chemical Physics, 2010, 133(20): 204502. doi: 10.1063/1.3507251 [39] BROWN E N, TRUJILLO C P, GRAY G T, et al. Soft recovery of polytetrafluoroethylene shocked through the crystalline phase Ⅱ-Ⅲ transition [J]. Journal of Applied Physics, 2007, 101(2): 024916. doi: 10.1063/1.2424536 [40] HUBER R C, WATKINS E B, DATTELBAUM D M, et al. In situ X-ray diffraction of high density polyethylene during dynamic drive: polymer chain compression and decomposition [J]. Journal of Applied Physics, 2021, 130(17): 175901. doi: 10.1063/5.0057439 [41] HUBER R C, PETERSON J, COE J D, et al. Polysulfone shock compressed above the decomposition threshold: velocimetry and modeling of two-wave structures [J]. Journal of Applied Physics, 2020, 127(10): 105902. doi: 10.1063/1.5124252 [42] DATTELBAUM D M, COE J D, KIYANDA C B, et al. Reactive, anomalous compression in shocked polyurethane foams [J]. Journal of Applied Physics, 2014, 115(17): 174908. doi: 10.1063/1.4875478 [43] KIPP M E, CHHABILDAS L C, REINHART W D, et al. Polyurethane foam impact experiments and simulations [J]. AIP Conference Proceedings, 2000, 505(1): 313–316. doi: 10.1063/1.1303481 [44] COE J D, LENTZ M, VELIZHANIN K A, et al. The equation of state and shock-driven decomposition of polymethy-lmethacrylate (PMMA) [J]. Journal of Applied Physics, 2022, 131(12): 125108. doi: 10.1063/5.0080369 [45] MAERZKE K A, COE J D, TICKNOR C, et al. Equations of state for polyethylene and its shock-driven decomposition products [J]. Journal of Applied Physics, 2019, 126(4): 045902. doi: 10.1063/1.5099371 [46] BOCK N, COFFEY D, WALLACE D C. Nonadiabatic contributions to the free energy from the electron-phonon interaction in Na, K, Al, and Pb [J]. Physical Review B, 2005, 72(15): 155120. doi: 10.1103/PhysRevB.72.155120 [47] BOCK N, WALLACE D C, COFFEY D. Adiabatic and nonadiabatic contributions to the free energy from the electron-phonon interaction for Na, K, Al, and Pb [J]. Physical Review B, 2006, 73(7): 075114. doi: 10.1103/PhysRevB.73.075114 [48] BENNETT B I. Computationally efficient expression for the zero-temperature isotherm in equations of state: LA-8616-MS [R]. Los Alamos: Los Alamos National Laboratory, 1980. [49] COWAN R D, ASHKIN J. Extension of the Thomas-Fermi-Dirac statistical theory of the atom to finite temperatures [J]. Physical Review, 1957, 105(1): 144–157. doi: 10.1103/PhysRev.105.144 [50] WUNDERLICH B. Thermal analysis of polymeric materials [M]. Berlin: Springer, 2005: 121−131. [51] ROSS M. A high-density fluid-perturbation theory based on an inverse 12th-power hard-sphere reference system [J]. The Journal of Chemical Physics, 1979, 71(4): 1567–1571. doi: 10.1063/1.438501 [52] COE J D, GAMMEL J T. A new 5-phase equation of state for carbon: LA-UR-16-26877 [R]. Los Alamos: Los Alamos National Laboratory, 2016. [53] ZHANG H R, ZHANG X X, FU X L, et al. Decomposition mechanisms of insensitive 2D energetic polymer TAGP using ReaxFF molecular dynamics simulation combined with Pyro-GC/MS experiments [J]. Journal of Analytical and Applied Pyrolysis, 2022, 162: 105453. doi: 10.1016/j.jaap.2022.105453 [54] JONES K, LANE J M D, MOORE N W. A reactive molecular dynamics study of phenol and phenolic polymers in extreme environments [J]. AIP Conference Proceedings, 2020, 2272(1): 070018. doi: 10.1063/12.0001031 [55] ISLAM M M, STRACHAN A. Decomposition and reaction of polyvinyl nitrate under shock and thermal loading: a ReaxFF reactive molecular dynamics study [J]. The Journal of Physical Chemistry C, 2017, 121(40): 22452–22464. doi: 10.1021/acs.jpcc.7b06154 [56] MATTSSON T R, LANE J M D, COCHRANE K R, et al. First-principles and classical molecular dynamics simulation of shocked polymers [J]. Physical Review B, 2010, 81(5): 054103. doi: 10.1103/PhysRevB.81.054103 [57] VAN DUIN A C T, DASGUPTA S, LORANT F, et al. ReaxFF: a reactive force field for hydrocarbons [J]. The Journal of Physical Chemistry A, 2001, 105(41): 9396–9409. doi: 10.1021/jp004368u -