SSS-MHD: a one-dimensional magneto-hydrodynamics multi-physics simulation platform for magnetically-driven high-energy-density dynamics experiments
-
摘要: 超高压、超高密度物质状态生成和性质研究是当代极端物理学的重要前沿领域，电磁驱动的高能量密度物理实验对于该领域的意义尤为重要。这类实验虽然形式上多种多样，但在物理上有内在统一性，即均以力学守恒定律和宏观电磁理论为基本框架。为了建立统一数值模拟平台、依靠负载电流实验数据（或驱动电路真实数据）确定各种极端实验条件下负载构形的力学运动及其与各个物理场的耦合问题，将经受大量实际检验的冲击、爆轰动力学和激光效应计算的一维拉格朗日编码SSS，实质性扩展成为磁流体力学多物理场耦合编码SSS-MHD。对于具有典型意义的平面准等熵斜波压缩、高速平面固体飞片发射、固体套筒电磁内爆和炸药内爆磁通量压缩实验等各类高能量密度动力学实验案例的模拟计算结果表明，编码SSS-MHD计算与美国Z装置、中国CQ和CJ系列装置的实验及美国编码ALEGRA-1D和2D计算数据的相对偏差基本不超过5%。该数值模拟平台为极端材料动力学实验（包括气体、液体、化合物和金属）提供了有力的支撑，还将有助于多维磁流体力学多物理场编码的开发。
-
关键词:
- 磁流体力学多物理场耦合计算 /
- 磁驱动斜波压缩 /
- 电磁驱动内爆动力学 /
- 爆炸磁通量压缩技术
Abstract: The research on generation and properties of materials under ultra-high pressure and density constitutes an important part of the extreme physics and hence a field of modern frontier science, especially the magnetically-driven high-energy-density physics herein is meaningful and in great need by core technologies. High pulsed power devices with tens-megampere output current and thousands-Tesla magnetic field were developed in past decades, e.g., the Z machine capable of 30 MA and 100 TW on load at the Sandia National Laboratories, USA; also a record intense magnetic field of 2800 T achieved with a cascade magneto-cumulative generator of MC-1 type at VNIIEF, Russia. It is available now to compress heavy metals up to 1 TPa or to launch thin Al flyer plates to super high speed over 45 km/s using isentropic compression experiments on the Z machine. Although these experiments take various forms, they have intrinsic unity in physics, which is based on the conservation laws of mechanics and the macroscopic electromagnetic theory. Therefore, it is feasible and necessary to establish a unified numerical simulation platform and determine the mechanical motion of the load configuration and its coupling with various physical fields under extreme experimental conditions by relying on the load current data (or the real data of the drive circuit). The magneto-hydrodynamics multi-physics codes have been successfully developed in USA, e.g., the excellent performance codes—ALEGRA series at the Sandia National Laboratories. This paper substantively extends the one-dimensional Lagrangian code SSS, which has been extensively validated by shock, detonation and laser radiation effect simulations, into a magneto-hydrodynamics multi-physics one and now it is renamed as SSS-MHD. The simulation results of various high-energy-density dynamic experiments with typical significance, such as planar quasi-isentropic ramp wave compression, ultra-high speed solid flyer launch, solid liner implosion, and explosively-driven magnetic flux compression, indicate that their relative deviations of the SSS-MHD simulations from the experimental data of America’s Z machine, China’s CQ and CJ series devices, and ALEGRA-1D/2D calculations are generally less than 5%. The SSS-MHD code turns into a powerful platform to simulate experiments of extreme material dynamics (including gases, liquids, metals and compounds) and its practice could be helpful to develop advanced multi-dimensional MHD multi-physics codes. -
内爆炸毁伤是打击水面舰船重要的毁伤方式，战斗部内爆炸产生的压力载荷作用于舰船舱壁结构，产生侧舷、横舱壁等结构变形或破坏，导致船舱进水，严重时可造成整船沉没。由于内爆炸载荷是一种多波峰叠加的复杂压力波，瞬态压力峰值高，载荷作用时间长，对船体结构的毁伤破坏比自由环境下的爆炸毁伤更剧烈、更复杂，因此，内爆炸载荷作用下结构的失效破坏成为了国内外学者研究的热点。
对于爆炸载荷作用下板结构的变形失效，已开展了大量的实验研究和理论分析。Nurick等[1]、Wen[2]开展了周边固支低碳钢圆板受均匀载荷作用的变形实验研究。Jacob等[3]、Gupta等[4]开展了固支圆板在空爆载荷作用下的边界撕裂破坏实验研究。Geretto等[5]开展了均布载荷作用下固支方板的失效模式实验。Teeling-Smith等[6]利用有限元分析软件开展了均布脉冲载荷作用下软钢薄圆板的变形撕裂研究，得到了板厚度、不同边界支撑等特性对板变形的影响。朱锡等[7]、陈长海等[8]开展了爆炸载荷作用下固支方板的变形实验研究。崔高领等[9]进行了均布冲击载荷对金属方板的仿真计算研究。这些研究均在空爆或半密闭环境下开展，而对于密闭环境下板结构变形与载荷关系的研究较少，如姚熊亮等[10]开展的舰船舱内爆炸载荷特征与板架毁伤规律仿真研究，侯海量等[11]开展的舱内爆炸载荷及舱室板架结构的失效模式研究。与空爆压力载荷相比，密闭环境下的压力载荷是一种多波峰叠加、准静态压力作用时间长的复杂压力波，受压后的板结构在产生塑性变形的过程中还存在极限变形[12]现象，即结构变形达到一定值后不会持续增大，后续阶段的压力载荷对结构变形不产生影响。这种产生板结构极限变形的压力载荷称为有效压力载荷或饱和压力载荷，研究有效压力载荷随时间的变化规律，掌握其与板结构的耦合作用关系，成为研究内爆炸结构变形破坏的关键。郑成等[13]、孔祥韶等[14]通过仿真计算方法，研究了密闭方形结构内方薄板的变形特性，并提出了方板变形预报模型和内爆炸压力载荷等效模型。而通过内爆炸实验方法研究板结构的变形特征与压力载荷关系，目前鲜有报道。
本文中基于双圆筒爆炸装置，开展TNT和某温压装药对铝质、钢质薄圆板的内爆炸实验，获取薄圆板的变形和压力载荷数据，分析薄圆板的破坏模式和压力载荷特点，并基于相同变形下作用载荷相等原理，确定薄圆板极限变形的有效比冲量和作用时间，提出圆薄板变形的预报模型并进行验证，以期为内爆炸结构毁伤研究提供技术参考。
1. 实验方法
1.1 实验装置及样品
实验双圆筒装置如图1所示，由两段圆筒、薄平板、前/后端盖、密封圈组成，在圆筒端部设有法兰盘。薄平板放置于两个圆筒端部的法兰盘之间，由螺钉将两个圆筒、薄平板密封连接，形成两个密闭圆筒腔室。圆筒长1 670 mm，内径800 mm，壁厚12 mm。为达到薄平板与法兰盘连接接近固支条件，在法兰盘的端面上铣出几组细浅同心圆，并对圆筒内端面角作倒圆角处理，减小其对薄板变形的影响。薄平板为2 mm厚的铝板和1 mm 厚的Q235钢板，在靠近边缘的同心圆周上，设置与圆筒法兰盘孔对应的圆孔。
图 1 双圆筒爆炸装置及测点布设位置Figure 1. Double cylinder explosive device and locations of measuring points实验装药为柱形TNT和某温压装药，长径比1∶1，其中TNT密度为1.58 g/cm3，质量分别为20和40 g；某温压装药（主要组成为HMX和Al粉）密度为1.84 g/cm3，质量分别为20、30、40、50、60、70和80 g。实验时，将药柱垂直悬挂于圆筒的几何中心，由电雷管起爆1 g C4扩爆药柱，引爆装药。
1.2 测点布局及测试系统
在爆炸圆筒一端圆形盖板的中心轴线上，按照从圆心到边缘的顺序均布4个压力测量点P1、P2、P3、P4，如图1(c)所示。测点布设的压力传感器为PCB公司113B24型壁面压力传感器，量程为0～6.9 MPa。传感器敏感面距盖板内表面约0.8 mm。安装完成后，在传感器头部敏感面与安装孔形成的凹槽内涂满油脂，使油脂表面与盖板表面平齐，其目的是延长爆炸热传递到传感器敏感面的时间，减小热冲击对压力载荷测量的影响。
基于爆心对称的两个刚性壁面上，冲击波在静止刚性壁面的反射压力与运动刚性壁面的反射压力差异不大[15]，可近似认为实验获取的盖板测点压力载荷等于圆薄板对称位置的压力载荷。
2. 实验结果及分析
2.1 圆板变形结果及分析
本文中共进行了10个工况实验，其中工况1～3为铝质薄圆板的变形破坏实验，工况4是为了补充获取工况2开展的40 g TNT装药在爆炸室的压力载荷数据实验，工况5～10为钢质薄圆板的变形破坏实验。典型铝质、钢质薄圆板实验后的状态如图2所示，各实验工况薄圆板的变形挠度及破坏情况如表1所示，其中WY代表某温压装药。
表 1 各实验工况圆板的变形结果Table 1. Deformation results of circular plate under various experimental conditions实验工况 平板材质 平板厚度/mm 装药类型 装药质量/g 圆板挠度/mm 变形破坏结果 1 铝 2 TNT 20 71.5 夹持边、中心塑性大变形 2 铝 2 TNT 40 夹持边被切断 3 铝 2 WY 20 93.0 夹持边间断性裂缝 4 钢 1 TNT 40 85.5 夹持边、中心塑性大变形 5 钢 1 WY 30 86.0 夹持边、中心塑性大变形 6 钢 1 WY 40 100.5 夹持边、中心塑性大变形 7 钢 1 WY 50 111.0 夹持边、中心塑性大变形 8 钢 1 WY 60 122.0 夹持边、中心塑性大变形 9 钢 1 WY 70 134.0 夹持边、中心塑性大变形 10 钢 1 WY 80 147.0 夹持边、中心塑性大变形 
图 2 典型铝质、钢质薄板实验后状态Figure 2. Pictures of typical aluminum and steel thin circular plates after experiments由图2可以看出，工况2的铝质薄圆板从夹持边界剪切断裂，切下圆板飞向邻室，落于后端盖底部，扭缩成团，如图2(b)所示，其剪切断面边缘有微小弯曲和大小不均匀的小尖刺。除工况2外，其他工况的铝质、钢质薄板均产生了四周顶起的锥形隆起结构，在薄板的夹持边界和中心区域产生了显著塑性大变形，在中心区域形成了比周边结构更陡峭的锥尖，且锥角随装药质量的增大而减小，更尖锐。在工况3的铝板夹持边界还产生了两段不连续的断裂裂缝，如图2(c)中红色圈内标注线及左侧放大图所示，从裂缝断面可以看出，夹持边界有显著的板厚拉伸减薄现象，断面边缘不规则，有长度不等的细长尖刺。工况5～10的装药为某温压药，钢质薄板产生了大变形，变形挠度随装药质量的增加而增大。
从上述圆板变形结果可以看出，在圆板夹持边界和中心产生了塑性大变形，在夹持边界还出现了局部撕裂和剪切断裂现象，主要原因是在圆板夹持边界和中心产生了应力极值点。参考方形薄板的应力响应分析[16]，在爆炸压力载荷作用下，固支夹持边界处形成弯曲转角，产生拉伸应变，当弯曲应变小于边界极限应变时，夹持边界产生大弯角，形成的运动塑性铰环从夹持边界向圆板中心移动收缩。在塑性铰环内部，圆板作刚体平行移动[17]，塑性铰环通过的外部则形成了曲面板，并按照锥形分布速度运动。当后续压力载荷作用下，圆板变形不持续增大时，塑性铰环收缩到圆板中心，形成锥尖；当弯曲应变接近板的极限应变时，在夹持边界产生局部的撕裂；当超过板的极限应变后，夹持边界产生横向剪切断裂。
圆板的实验结果表明，在内爆炸载荷作用下，薄圆板产生了3种破坏模式：塑性大变形、拉伸撕裂、横向剪切。该结果与文献[14]所得到的空爆压力载荷作用下梁、圆板的变形破坏模式一致。
2.2 圆板上的压力载荷及分析
工况7的4个测点测量的压力载荷曲线和积分处理后的比冲量曲线如图3所示。从图3(a)的压力载荷曲线可以看出，4个压力测点的压力载荷到达时间差小于0.12 μs，P1、P2、P3测点最大压力峰值基本一致，角隅测点P4的最大峰值约为其他测点最大峰值的2倍，4个测点的压力载荷随时间的变化趋势基本一致。从图3(b)的比冲量曲线可以看出：在3 ms内，4个测点的比冲量大小与比冲量增长量基本一致，虽然P4测点的最大压力峰值最高，但是比冲量与其他测点基本相当；3 ms后，除P3测点外，其他测点的比冲量增长量基本一致。观察P3测点的压力载荷曲线发现，3 ms后，其压力值比其他测点略小，原因可能是受到爆炸热冲击对压力传感器性能的影响，由于压电型传感器具有热释电效应，当传感器的隔热防护作用不够时，爆炸热通过热传导作用到传感器的敏感元件上，输出热冲击干扰信号，并且传感器的预紧结构受热膨胀，预紧力减小，导致传感器的灵敏度变小，输出的压电信号也减小。
图 3 工况7的4个测点压力载荷曲线和比冲量曲线Figure 3. Pressure load curves and specific impulse curves of four measuring points under condition 7其他实验工况下4个测点的载荷特性与工况7基本类似，除角隅位置压力载荷的最大峰值是其他测点的2倍外，其他时刻4个测点的压力载荷大致相等，比冲量也基本一致。因此可由一个测点的压力载荷反映整个圆板面上的压力载荷。各工况圆板中心测点P1的压力载荷曲线如图4所示。
图 4 各实验工况圆板中心的压力载荷曲线Figure 4. Pressure load curves of circular plates under various experimental conditions从图4所示的压力载荷曲线可以看出，薄圆板受到多波峰叠加的复杂压力载荷作用。在初始阶段，压力波以显著的一大一小波峰组合周期性衰减出现，变化周期为2.8～3.5 ms（药量越大周期越短），每一周期内的波峰值均比前一周期对应的波峰值小，3个周期后，两个波峰逐渐融合，形成了压力波动逐渐减小的准静态压力。工况2的压力曲线在7.9 ms内，波形特征与其他工况的变化规律类似，7.9 ms后压力逐渐减小，形成约12 ms的周期性振荡波。工况5～10的压力载荷最大峰值和匀化形成的准静态压力均随着装药质量的增大而增大，但振荡周期随着装药质量的增大而缩短，从3.5 ms逐渐缩短到2.8 ms。
将上述各工况圆板中心测点P1的压力载荷积分处理，得到比冲量时间曲线，如图5所示。
图 5 各实验工况圆板中心比冲量曲线Figure 5. Specific impulse curves of circular plates under various experimental conditions除工况2外，其他工况的比冲量时间曲线变化规律类似，在约0.35 ms内，比冲量曲线急剧上升，随后呈波浪式增长，并且随着时间的延长，波浪振幅逐渐减小，逐渐趋于线性增长。这主要是由于爆炸室内压力波趋于匀化，形成了压力波动较小的准静态压力所致。比冲量波动周期等于压力载荷波动周期，并且随着装药质量的增加，比冲量也随之增大。
工况2与工况4的装药类型和装药质量相等，比冲量值在7.9 ms内基本一致，7.9 ms后，工况2的比冲量增长率相对减缓，形成了周期约12 ms的新波动。结合工况2薄圆板的破坏结果，可得到在7.9 ms时刻，铝板受到的比冲量为1.720 MPa·ms，铝板夹持边界的弯曲应变超过板的极限应变，产生了瞬时剪切断裂现象。工况3的铝板夹持边界产生了两小段撕裂破坏，其破坏时刻的弯曲应变接近铝板的极限应变，受到的比冲量应小于1.720 MPa·ms。工况5～10的钢质薄圆板，在压力载荷作用时间内，夹持边界产生的弯曲应变均小于材料的极限应变。
基于圆板上比冲量时间曲线呈线性增长的特点，其增长量由圆筒内形成的准静态压力产生，因此，除去比冲量曲线前期0.35 ms内的急剧上升段，选取0.35～20.0 ms的比冲量数据作线性拟合处理，得到比冲量平均增长率，该增长率近似等于圆筒内的准静态压力。由于急剧上升段的比冲量主要由压力载荷的最大峰值脉冲积分产生，脉冲持续时间短，可以近似认为圆板受到的比冲量瞬时增长到
Δi ，然后线性增长，因此，圆板上的比冲量i可表示为：i=pt+Δi (1) 式中：p为圆筒内准静态压力，
t 为压力载荷作用时间，Δi 为比冲量快速增长量。对各工况的比冲量数据进行处理，可得到各工况的准静态压力p和比冲量快速增长量
Δi 。表2为工况5～10的准静态压力p和Δi 以及3、5和10 ms时的比冲量实验值和依据式(1)得到的计算值。表 2 工况5～10的 p和Δi及典型时刻模型比冲量的计算值和实验值Table 2. p and Δiof conditions 5−10 and the calculated and experimental specific impulses at typical times 实验
工况压力/
MPa比冲量快速增长量/
（MPa·ms）3 ms 时的比冲量/（MPa·ms） 5 ms 时的比冲量/（MPa·ms） 10 ms 时的比冲量/（MPa·ms） 计算值 实验值 计算值 实验值 计算值 实验值 5 0.199 0.290 0.861 0.758 1.236 1.264 2.172 2.217 6 0.235 0.332 1.056 0.970 1.512 1.409 2.653 2.568 7 0.294 0.393 1.237 1.200 1.769 1.701 3.097 3.244 8 0.322 0.471 1.408 1.289 2.010 1.856 3.515 3.520 9 0.342 0.518 1.571 1.424 2.240 2.123 3.912 3.927 10 0.377 0.566 1.727 1.648 2.460 2.446 4.291 4.342 对比表2中的比冲量计算结果和实验结果可知，3、5和10 ms时段内的计算结果与实验结果的偏差分别为13.6%、8.3%和4.5%。可以看出：随着载荷作用时间的加长，模型计算结果与实验结果的偏差减小。这表明比冲量计算模型可以近似预估圆板受到的比冲量，且载荷作用时间越长，计算结果越接近实验结果。
Barker等[18]提出了一种内爆炸准静态压力p计算模型：
p=A(m/V)γ (2) 式中：
m 为炸药装药质量，V 为内爆炸空间体积，A 、γ 为相关因子。将上述工况5～10的准静态压力p的实验结果、实验装药、圆筒体积，按照式(2)的幂指数模型进行拟合，可以得到该装药的准静态压力计算模型，其中：系数
A =1.865，γ =0.67。依据该计算模型得到工况5～10的准静态压力p为0.200、0.243、0.282、0.319、0.353、0.386 MPa，与表2中压力的实验结果进行比对可知，偏差小于4.7%。上述快速增长比冲量
Δi ，主要由压力载荷曲线中起始段最大峰值的压力脉冲积分得到，该脉冲的最大峰值、脉宽与装药质量及爆心距相关，因此，按照比冲量与比对距离成幂指数关系：Δi=K(3√m/r)δ (3) 式中：K、δ为因数，r为内爆炸爆心距。将表2中工况5～10的
Δi 实验值、装药质量、爆心距按照式（3）的幂指数模型进行拟合处理，得到K=2.051，δ=2.11。依据该模型计算工况5～10的Δi 分别为0.299、0.372、0.441、0.506、0.569、0.630 MPa·ms，与表2中Δi 的实验结果比对可知，偏差小于3.3%，表明该模型能够描述冲量曲线中快速增长段的比冲量Δi 。将式(2)、式(3)代入式(1)，则式(1)转化为：
i=A(m/V)γt+K(3√m/r)δ (4) 2.3 钢质圆板变形的有效载荷分析
Zhao等[12]提出了平板变形的极限响应现象，即在平板变形达到一定值后不会持续增加，后续的作用载荷对平板的变形不产生影响，存在平板极限变形的有效载荷及作用时间。对于薄板结构，一般遵循冲量毁伤准则，即对薄板的破坏主要由冲量载荷产生。在空爆环境下，冲量值可通过冲击波压力正压时间内的比冲量与板面积的乘积得到，而在内爆炸环境下，由于压力载荷为多波峰叠加并逐渐匀化为准静态压力的复杂波，正压时间长，无法确定载荷作用时间，因此不能依据获取的压力载荷数据直接得到有效比冲量和作用时间。
上述内爆炸压力载荷未完全匀化前的压力波，可近似看作一种瞬态升压后又缓慢衰减的单波峰压力波（忽略内部压力波动），与外爆炸压力波波形特点相似。因此，可基于相同变形下冲量载荷相等原理，认为空爆条件下圆板受到的冲量载荷等于内爆炸条件下圆板极限变形所受到的有效冲量，再依据有效冲量，得到载荷作用时间。
在空爆条件下，Jacob等[19]提出了钢质圆板的变形挠度计算模型：
δ0H=0.425ϕc (5) 式中：δ0为圆板变形挠度；H为圆板厚度；
ϕc 为无量纲数，其表达式为：ϕc=IπRH2√ρσ0 (6) 式中：
I 为圆板受到的冲量，R 为圆板半径，ρ 为平板材料的密度，σ0 为平板材料的静态屈服应力。空爆作用下圆板上的冲量
I 等于内爆炸作用下圆板上的有效冲量I0 ，因此，式(5)、式(6)变形后，得到内爆炸圆板上的有效冲量：I0 =δ0H×πRH2√ρσ00.425=πδ0RH√ρσ00.425 (7) 从图3中圆板各测点的比冲量曲线可以得到，圆板面上各测点的比冲量基本相等，随时间加长比冲量的增长率也基本一致，表明圆板面上各位置的瞬态比冲量基本一致，随时间的增长比冲量增量也基本相当，特别是装药质量越大，一致性越好。因此，圆板上任意一点的有效冲量可通过测量点的有效比冲量与圆板面积之积得到：
I0=i0s (8) 式中：
i0 为圆板测量点的有效比冲量，s 为圆板面积。将工况5～10中钢质圆板的变形挠度代入式(7)，并结合式(8)，得到钢质圆板的有效冲量
I0 、有效比冲量i0 ，如表3所示，其中钢板密度为7 800 kg/m3，σ0 =235 MPa。表 3 工况5～10钢质圆板的有效冲量和有效比冲量Table 3. Effective impulses and effective specific impulses of steel circular plates under conditions 5−10实验工况 圆板挠度/mm 有效冲量/(N⋅s) 有效比冲量/(MPa⋅ms) 修正有效冲量/(N⋅s) 修正有效比冲量/(MPa⋅ms) 作用时间/ms 5 86.0 344 0.656 389 0.777 2.93 6 100.5 402 0.767 456 0.912 2.76 7 111.0 444 0.847 505 1.009 2.55 8 122.0 476 0.908 542 1.084 2.43 9 134.0 536 1.023 612 1.224 2.36 10 147.0 588 1.122 673 1.345 2.26 式(6)中忽略了爆炸动态加载过程中材料应变率的影响，在瞬态强压力载荷作用下，平板结构的动力响应与材料的应变率显著相关，因此，需要考虑材料应变率变化的影响。Cowper-Symonds提出了材料动态屈服应力与静态屈服应力的关系模型[19]：
σdσ0=1+(v0δ3√2˙ε0R2)1/q (9) 式中：
σd 为平板材料的动态屈服应力，v0 为平板的瞬态速度，˙ε0 、q 为与材料相关的常数。按照动量守恒，平板上的比冲量i 与v0 的关系为i=ρδv0 ，代入式(9)得到：σdσ0=1+(i3√2ρ˙ε0R2)1/q (10) 对于Q235低碳钢材料，
˙ε0 一般取值[19]为40.4 s−1，q 取值为5。将表3中的有效比冲量i0 代入式(10)，得到对应的σd ，再将σd 代入式(7)，得到修正后的有效冲量Im 和有效比冲量im 。依据修正后的有效比冲量
im ，在各工况比冲量时间曲线上，由im 得到对应有效比冲量作用时间t ，如表3所示。由表3可以看出：随着圆板变形挠度的增大，有效比冲量im 随之增大，有效载荷作用时间t 随着圆板变形挠度的增大而缩短。将有效比冲量
im 和作用时间t 进行数据拟合处理，得到有效比冲量im 与作用时间t 的关系式：t=75.64i−0.49m (11) 因此，当确定了圆板材料、尺度和极限变形量后，可以通过式(7)和式(10)得到有效比冲量，再依据式(4)和式(11)预估圆板极限变形下的装药质量。
依据上述方法，预估工况5～10钢质圆板极限变形下的装药量，分别为26、35、44、52、63和75 g，与实验装药量比较可知，偏差小于13.3%，表明该模型可以预估内爆炸作用下圆板的变形情况。
3. 结 论
通过双圆筒型装置内薄圆板的爆炸实验，研究了圆板的变形破坏模式和比冲量载荷变化规律，并基于相同变形下的冲量载荷相等原理，确定了内爆炸圆板极限变形下的有效比冲量和作用时间，提出了薄圆板变形的预报模型，并进行了验证，结果表明：
（1）在内爆炸压力载荷作用下，圆板的夹持边界和中心区域是应力集中区，圆板产生了塑性大变形、拉伸撕裂和横向剪切3种破坏模式。
（2）圆板上的内爆炸压力载荷为多波峰的复杂压力波，初始阶段比冲量陡峭增长，随后由波浪式增长逐渐匀化为线性增长，增长率近似等于装药在爆炸室内形成的准静态压力。
（3）在实验圆筒装置内，当1 mm厚钢质圆板产生86～147 mm的变形挠度时，有效比冲量在0.777～1.345 MPa·ms范围，有效载荷作用时间在2.26～2.93 ms范围。通过验证，由圆板变形预估模型得到的装药质量与实验装药质量偏差小于13.3%。
-
图 1 适合SSS-MHD编码模拟的各种类型电磁驱动高能量密度物理实验
Figure 1. Kinds of the magnetically-driven high-energy-density physics experiments suitably simulated with the SSS-MHD code
图 5 SSS-MHD编码中计算程序包MHDBLK的组成
Figure 5. Flow chart for the routine package MHDBLK in the SSS-MHD code
图 6 磁驱动准等熵压缩实验的典型负载构形
Figure 6. Typical loading configurations for magnetically-driven isentropic compression experiments (ICE)
图 8 对于流体物理研究所CQ装置磁驱动准等熵压缩实验的SSS-MHD模拟（算例2～4）
Figure 8. SSS-MHD simulations for magnetically-driven quasi-isentropic compression experiments at IFP (examples 2－4)
图 9 对于Sandia实验室Z装置铝飞片实验“11 mm-2 s”的SSS-MHD模拟（算例5）
Figure 9. SSS-MHD simulations for the Al flyer experiment “11 mm-2 s”on the Z machine at the Sandia National Laboratories, USA (example 5)
图 10 对于Sandia实验室Z装置上金属套筒电磁内爆实验的SSS-MHD模拟（算例6）
Figure 10. SSS-MHD simulation of magnetically-driven liner implosions on the Z-machine at the Sandia National Laboratories, USA (example 6)
图 11 CJ-100型MC-1发生器的示意图和实物照片
Figure 11. Schematics and picture of the CJ-100 type MC-1 generator
图 12 CJ-100型MC-1发生器磁通量压缩实验 及其性能的SSS-MHD模拟计算（算例7）
Figure 12. SSS-MHD simulations for the magnetic flux compression experiment and the performances of the CJ-100 type MC-1 generator (example 7)
表 1 磁驱动准等熵压缩实验算例2～4的主要参数[14, 16]
Table 1. Parameters of examples 2－4 for magnetically-driven isentropic compression experiments[14, 16]
算例 实验 充电电压/kV 电极板 台阶样品 光学窗口 速度峰值/(km·s−1) 压力峰值/GPa 材料 厚度/mm 材料 厚度/mm 实验 电流计算 电路计算 电流计算 电路计算 2 CQ4-13251[14, 16] 85 Cu 1.350 Ta 1.013 LiF 1.589 1.607 1.616 80.7 81.9 1.350 1.201 1.604 1.586 1.575 80.5 81.6 3 CQ4-12143[14, 16] 75 Al 2.780 0.774 0.770 0.761 6.30 6.08 3.010 0.779 0.763 0.757 6.31 6.09 4 CQ1.5-11202[14, 16] 60 Al 0.353 JO-9159 0.599 LiF 0.306 0.292 0.313 4.27 4.61 0.353 0.789 0.333 0.318 0.326 4.63 4.82 
下载: 导出CSV
-
[1] HATFIELD P W, GAFFNEY J A, ANDERSON G J, et al. The data-driven future of high-energy-density physics [J]. Nature, 2021, 593(7859): 351–361. DOI: 10.1038/s41586-021-03382-w. [2] DRAKE R P. 高能量密度物理: 基础、惯性约束聚变和实验天体物理学 [M]. 孙承纬, 译. 北京: 国防工业出版社, 2013. [3] 孙承纬. 电磁加载下的高能量密度物理问题研究 [J]. 高能量密度物理, 2007(1): 41–46. [4] 孙承纬, 赵剑衡, 王桂吉, 等. 磁驱动准等熵平面压缩和超高速飞片发射实验技术原理、装置及应用 [J]. 力学进展, 2012, 42(2): 206–219. DOI: 10.6052/1000-0992-2012-2-20120208.SUN C W, ZHAO J H, WANG G J, et al. Progress in magnetic loading techniques for isentropic compression experiments and ultra-high velocity flyer launching [J]. Advances in Mechanics, 2012, 42(2): 206–219. DOI: 10.6052/1000-0992-2012-2-20120208. [5] OLIPHANT T A, WITTE K H. RAVEN code: LA-10826 [R]. USA: Los Alamos National Laboratory, 1987. [6] ROBINSON A C, BRUNNER T A, CARROLL S, et al. ALEGRA: an arbitrary lagrangian-eulerian multimaterial, multiphysics code [C] // Proceedings of the 46th AIAA Aerospace Sciences Meeting and Exhibit. Reno, Nevada, USA, 2008: 1235. DOI: 10.2514/6.2008-1235. [7] 孙承纬. 一维冲击波和爆轰波计算程序SSS [J]. 计算物理, 1986, 3(2): 142–154. DOI: 10.19596/j.cnki.1001-246x.1986.02.002.SUN C W. SSS: a code for computing one dimensional shock and detonation wave propagation [J]. Chinese Journal of Computational Physics, 1986, 3(2): 142–154. DOI: 10.19596/j.cnki.1001-246x.1986.02.002. [8] GU Z W, SUN C W, ZHAO J H, et al. One-dimensional numerical simulation of laser-driven flyer plates [J]. Journal of Applied Physics, 2004, 96(6): 3486–3490. DOI: 10.1063/1.1781765. [9] YUAN H, TONG H F, LI M, et al. Computational study of nanosecond pulsed laser ablation and the application to momentum coupling [J]. Journal of Applied Physics, 2012, 112: 023105. DOI: 10.1063/1.4737188. [10] 刘启泰, 孙承纬, 胡熙静. SSS程序在电磁内爆中的数值模拟 [J]. 高压物理学报, 2004, 18(2): 183–187. DOI: 10.11858/gywlxb.2004.02.015.LIU Q T, SUN C W, HU X J. The numerical simulation of solid liner implosion in SSS code [J]. Chinese Journal of High Pressure Physics, 2004, 18(2): 183–187. DOI: 10.11858/gywlxb.2004.02.015. [11] 贺佳, 赵剑衡, 谭福利, 等. 电爆炸箔驱动绝缘飞片的一维数值模拟 [J]. 高压物理学报, 2009, 23(6): 476–480. DOI: 10.11858/gywlxb.2009.06.013.HE J, ZHAO J H, TAN F L, et al. One-dimensional numerical simulation of a flyer accelerated by electrically exploded metal foil [J]. Chinese Journal of High Pressure Physics, 2009, 23(6): 476–480. DOI: 10.11858/gywlxb.2009.06.013. [12] 孙承纬, 赵剑衡, 刘仓理, 等. 材料准等熵压缩实验研究进展 [M]. 北京: 中国原子能出版社, 2015. [13] 孙承纬, 赵剑衡, 王桂吉, 等. 磁驱动准等熵压缩和高速飞片的实验研究 [C] // 庆祝中国力学学会成立50周年暨中国力学学会学术大会. 北京, 2007: 17. [14] 赵继波, 孙承纬, 罗斌强, 等. 磁驱动等熵压缩实验构形的磁流体力学计算模拟 [J]. 力学学报, 2014, 46(5): 685–693. DOI: 10.6052/0459-1879-14-019.ZHAO J B, SUN C W, LUO B Q, et al. Magneto-hydrodynamics modeling of configurations in magnetically driven isentropic compression experiments [J]. Chinese Journal of Theoretical and Applied Mechanics, 2014, 46(5): 685–693. DOI: 10.6052/0459-1879-14-019. [15] 赵继波, 孙承纬, 谷卓伟, 等. 爆轰驱动固体套筒压缩磁场计算及准等熵过程分析 [J]. 物理学报, 2015, 64(8): 080701. DOI: 10.7498/aps.64.080701.ZHAO J B, SUN C W, GU Z W, et al. Magneto-hydrodynamic calculation of magnetic flux compression with explosion driven solid liners and analysis of quasi-isentropic process [J]. Acta Physica Sinica, 2015, 64(8): 080701. DOI: 10.7498/aps.64.080701. [16] ZHAO J B, SUN C W, LUO B Q, et al. Loading circuit coupled magnetohydrodynamic simulation of sample configurations in isentropic compression experiments [J]. IEEE Transactions on Plasma Science, 2015, 43(4): 1068–1076. DOI: 10.1109/TPS.2015.2405575. [17] 王桂吉, 罗斌强, 陈学秒, 等. 磁驱动平面准等熵加载装置、实验技术及应用研究新进展 [J]. 爆炸与冲击, 2021, 41(12): 121403. DOI: 10.11883/bzycj-2021-0119.WANG G J, LUO B Q, CHEN X M, et al. Recent progress on the experimental facilities, techniques and applications of magnetically driven quasi-isentropic compression [J]. Explosion and Shock Waves, 2021, 41(12): 121403. DOI: 10.11883/bzycj-2021-0119. [18] 罗斌强, 陈学秒, 王桂吉, 等. 磁驱动压-剪联合加载下材料动态强度的直接测量 [J]. 中国科学: 物理学 力学 天文学, 2016, 46(11): 114601. DOI: 10.1360/SSPMA2016-00184.LUO B Q, CHEN X M, WANG G J, et al. 磁驱动压-Direct measurement of material dynamic strength under high pressure using magnetically driven pressure-shear loading [J]. Scientia Sinica: Physica, Mechanica and Astronomica, 2016, 46(11): 114601. DOI: 10.1360/SSPMA2016-00184. [19] LUO B Q, SUN C W, ZHAO J H, et al. Unified numerical simulation of metallic foil electrical explosion and its applications [J]. IEEE Transactions on Plasma Science, 2013, 41(1): 49–57. DOI: 10.1109/TPS.2012.2227827. [20] LUO B Q, WANG G J, MO J J, et al. Verification of conventional equations of state for tantalum under quasi-isentropic compression [J]. Journal of Applied Physics, 2014, 116: 193506. DOI: 10.1063/1.4902064. [21] ZHANG X P, WANG G J, ZHAO J H, et al. High velocity flyer plates launched by magnetic pressure on pulsed power generator CQ-4 and applied in shock Hugoniot experiments [J]. Review of Scientific Instruments, 2014, 85: 055110. DOI: 10.1063/1.4875705. [22] ZHAO X M, SUN C W, SUN Q Z, et al. Simulation on the compressed field-reversed configuration with alpha particle self-heating [J]. Plasma Physics and Controlled Fusion, 2019, 61(7): 075015. DOI: 10.1088/1361-6587/ab1e84. [23] REINOVSKY R E, ATCHISON W L, DIMONTE G, et al. Pulsed-power hydrodynamics: an application of pulsed-power and high magnetic fields to the exploration of material properties and problems in experimental hydrodynamics [J]. IEEE Transactions on Plasma Science, 2008, 36(1): 112–124. DOI: 10.1109/TPS.2007.914708. [24] SINARS D B, SWEENEY M A, ALEXANDER C S, et al. Review of pulsed power-driven high energy density physics research on Z at Sandia [J]. Physics of Plasmas, 2020, 27(7): 070501. DOI: 10.1063/5.0007476. [25] SLUTZ S A, STYGAR W A, GOMEZ M R, et al. Scaling magnetized liner inertial fusion on Z and future pulsed-power accelerators [J]. Physics of Plasmas, 2016, 23(2): 3933. DOI: 10.1063/1.4941100. [26] 吴其芬, 李桦. 磁流体力学 [M]. 长沙: 国防科技大学出版社, 2007: 23–31. [27] JOHNSON J D. The Sesame database [C] // Proceedings of the 12th Symposium on Thermophysics Properties. Boulder, Colorado, USA, 1994: 1–12. [28] STEINBERG D J, COCHRAN S G, GUINAN M W. A constitutive model for metals applicable at high-strain rate [J]. Journal of Applied Physics, 1980, 51(3): 1498–1504. DOI: 10.1063/1.327799. [29] 赵继波. 一维磁流体动力学程序SSS-MHD研究和实验构形模拟计算 [D]. 四川绵阳:中国工程物理研究院, 2014: 78.ZHAO J B. Research on one-dimensional magneto-hydrodynamics code SSS-MHD and simulation on confingurations on experiments [D]. Mianyang, Sichuan, China: China Academy of Engineering Physics, 2014: 78. [30] 陆禹. 柱筒构型的电磁内爆、压缩实验与磁流体力学数值模拟 [D]. 合肥: 中国科学技术大学, 2022: 26–29.LU Y. Experiments and MHD simulations on magnetically diriven implosion and compression of cylindrical configurations [D]. Hefei, Anhui, China: China University of Science and Technology, 2022: 26–29. [31] 黄奕勇, 李星辰, 田野, 等. COMSOL多物理场仿真入门指南 [M]. 北京: 机械工业出版社, 2021. [32] LEMKE R W, KNUDSON M D, ROBINSON A C, et al. Self-consistent, two-dimensional, magnetohydrodynamic simulations of magnetically driven flyer plates [J]. Physics of Plasmas, 2003, 10(5): 1867–1874. DOI: 10.1063/1.1557530. [33] VOGLER T J, AO T, ASAY J R. High-pressure strength of aluminum under quasi-isentropic loading [J]. International Journal of Plasticity, 2009, 25(4): 671–694. DOI: 10.1016/j.ijplas.2008.12.003. [34] LEMKE R W, DOLAN D H, DALTON D G, et al. Probing off-Hugoniot states in Ta, Cu, and Al to 1000GPa compression with magnetically driven liner implosions [J]. Journal of Applied Physics, 2016, 119(1): 015904. DOI: 10.1063/1.4939675. [35] WANG G J, SUN C W, TAN F L, et al. The compact capacitor bank CQ-1.5 employed in magnetically driven isentropic compression and high velocity flyer plate experiments [J]. Review of Scientific Instruments, 2008, 79: 053904. DOI: 10.1063/1.2920200. [36] LEMKE R W, KNUDSON M D, DAVIS J P. Magnetically driven hyper-velocity launch capability at the Sandia Z accelerator [J]. International Journal of Impact Engineering, 2011, 38(6): 480–485. DOI: 10.1016/j.ijimpeng.2010.10.019. [37] ALTGILBERS L L. 磁通量压缩发生器 [M]. 孙承纬, 译. 北京: 国防工业出版社, 2008: 1–5. [38] 谷卓伟, 罗浩, 张恒第, 等. 炸药柱面内爆磁通量压缩实验技术研究 [J]. 物理学报, 2013, 62(17): 170701. DOI: 10.7498/aps.62.170701.GU Z W, LUO H, ZHANG H D, et al. Experimental research on the technique of magnetic flux compression by explosive cylindrical implosion [J]. Acta Physica Sinica, 2013, 62(17): 170701. DOI: 10.7498/aps.62.170701. [39] ZHOU Z Y, GU Z W, LUO H, et al. A compact explosive-driven flux compression generator for reproducibly generating multimegagauss fields [J]. IEEE Transactions on Plasma Science, 2018, 46(10): 3279–3283. DOI: 10.1109/TPS.2018.2794761. 期刊类型引用(4)
1. 张学瑞，周涛. 密闭空间中复合装药的能量释放特性. 爆炸与冲击. 2024(06): 41-49 . 
本站查看2. 赵豫熙，袁浩天，王须民，张之凡. 水下接触爆炸作用下金属/CFRP复合层合板的防护性能. 高压物理学报. 2024(06): 177-190 . 
百度学术3. 赵伟成，翟红波，毛伯永，杨峰. 两装药同时起爆时金属靶板的动态响应. 爆炸与冲击. 2023(12): 3-13 . 本站查看4. 李荣，万家平，蔡灵，谷家扬，刘涛. 基于CONWEP法空爆载荷下船体结构动态响应研究. 江苏船舶. 2022(04): 1-5 . 百度学术其他类型引用(2)
-
下载:
百度学术
计量
- 文章访问数: 322
- HTML全文浏览量: 65
- PDF下载量: 74
- 被引次数: 6