高压物理学报

电爆炸金属桥箔早期过程中电磁-热-力多物理场耦合建模与分析

PDF下载

王刚华 , 谢龙 , 肖波 , 王强 , 唐久棚 , 欧海彬 , 阚明先 , 段书超

高压物理学报 | 高压理论与计算 2024,38(1): 012301-1-012301-9

收起

高压物理学报 |高压理论与计算 2024 , 38 (1) : 012301-1 -012301-9

电爆炸金属桥箔早期过程中电磁-热-力多物理场耦合建模与分析

全屏

[Author(id=1325132988484297200, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, orderNo=0, firstName=null, middleName=null, lastName=null, nameCn=null, orcid=null, stid=null, country=null, authorPic=null, dead=0, email=wanggh@caep.cn, emailSecond=null, emailThird=null, correspondingAuthor=0, authorType=1, ext={EN=AuthorExt(id=1325132988891144690, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, authorId=1325132988484297200, language=EN, stringName=Ganghua WANG, firstName=Ganghua, middleName=null, lastName=WANG, prefix=null, suffix=null, authorComment=null, nameInitials=null, affiliation=null, department=null, xref=null, address=Institute of Fluid Physics, CAEP, Mianyang 621999, Sichuan, China, bio=null, bioImg=null, bioContent=null, aboutCorrespAuthor=null), CN=AuthorExt(id=1325132989444792819, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, authorId=1325132988484297200, language=CN, stringName=王刚华, firstName=null, middleName=null, lastName=null, prefix=null, suffix=null, authorComment=null, nameInitials=null, affiliation=null, department=null, xref=null, address=中国工程物理研究院流体物理研究所, 四川绵阳　621999, bio={"content":"

王刚华（1976－），男，博士，副研究员，主要从事电磁驱动多物理场耦合理论与数值方法研究. E-mail：wanggh@caep.cn

"}, bioImg=null, bioContent=

王刚华（1976－），男，博士，副研究员，主要从事电磁驱动多物理场耦合理论与数值方法研究. E-mail：wanggh@caep.cn

, aboutCorrespAuthor=null)}, companyList=[AuthorCompany(id=1325132987817402860, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, xref=null, ext=[AuthorCompanyExt(id=1325132987855151597, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, companyId=1325132987817402860, language=EN, country=null, province=null, city=null, postcode=null, companyName=null, departmentName=null, remark=Institute of Fluid Physics, CAEP, Mianyang 621999, Sichuan, China), AuthorCompanyExt(id=1325132987980980718, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, companyId=1325132987817402860, language=CN, country=null, province=null, city=null, postcode=null, companyName=null, departmentName=null, remark=中国工程物理研究院流体物理研究所, 四川绵阳　621999)])]), Author(id=1325132989876806133, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, orderNo=1, firstName=null, middleName=null, lastName=null, nameCn=null, orcid=null, stid=null, country=null, authorPic=null, dead=0, email=null, emailSecond=null, emailThird=null, correspondingAuthor=0, authorType=1, ext={EN=AuthorExt(id=1325132990568866295, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, authorId=1325132989876806133, language=EN, stringName=Long XIE, firstName=Long, middleName=null, lastName=XIE, prefix=null, suffix=null, authorComment=null, nameInitials=null, affiliation=null, department=null, xref=null, address=Institute of Fluid Physics, CAEP, Mianyang 621999, Sichuan, China, bio=null, bioImg=null, bioContent=null, aboutCorrespAuthor=null), CN=AuthorExt(id=1325132992670212601, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, authorId=1325132989876806133, language=CN, stringName=谢龙, firstName=null, middleName=null, lastName=null, prefix=null, suffix=null, authorComment=null, nameInitials=null, affiliation=null, department=null, xref=null, address=中国工程物理研究院流体物理研究所, 四川绵阳　621999, bio=null, bioImg=null, bioContent=null, aboutCorrespAuthor=null)}, companyList=[AuthorCompany(id=1325132987817402860, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, xref=null, ext=[AuthorCompanyExt(id=1325132987855151597, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, companyId=1325132987817402860, language=EN, country=null, province=null, city=null, postcode=null, companyName=null, departmentName=null, remark=Institute of Fluid Physics, CAEP, Mianyang 621999, Sichuan, China), AuthorCompanyExt(id=1325132987980980718, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, companyId=1325132987817402860, language=CN, country=null, province=null, city=null, postcode=null, companyName=null, departmentName=null, remark=中国工程物理研究院流体物理研究所, 四川绵阳　621999)])]), Author(id=1325132993051894267, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, orderNo=2, firstName=null, middleName=null, lastName=null, nameCn=null, orcid=null, stid=null, country=null, authorPic=null, dead=0, email=null, emailSecond=null, emailThird=null, correspondingAuthor=0, authorType=1, ext={EN=AuthorExt(id=1325132993496490493, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, authorId=1325132993051894267, language=EN, stringName=Bo XIAO, firstName=Bo, middleName=null, lastName=XIAO, prefix=null, suffix=null, authorComment=null, nameInitials=null, affiliation=null, department=null, xref=null, address=Institute of Fluid Physics, CAEP, Mianyang 621999, Sichuan, China, bio=null, bioImg=null, bioContent=null, aboutCorrespAuthor=null), CN=AuthorExt(id=1325132994041750014, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, authorId=1325132993051894267, language=CN, stringName=肖波, firstName=null, middleName=null, lastName=null, prefix=null, suffix=null, authorComment=null, nameInitials=null, affiliation=null, department=null, xref=null, address=中国工程物理研究院流体物理研究所, 四川绵阳　621999, bio=null, bioImg=null, bioContent=null, aboutCorrespAuthor=null)}, companyList=[AuthorCompany(id=1325132987817402860, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, xref=null, ext=[AuthorCompanyExt(id=1325132987855151597, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, companyId=1325132987817402860, language=EN, country=null, province=null, city=null, postcode=null, companyName=null, departmentName=null, remark=Institute of Fluid Physics, CAEP, Mianyang 621999, Sichuan, China), AuthorCompanyExt(id=1325132987980980718, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, companyId=1325132987817402860, language=CN, country=null, province=null, city=null, postcode=null, companyName=null, departmentName=null, remark=中国工程物理研究院流体物理研究所, 四川绵阳　621999)])]), Author(id=1325132994595398144, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, orderNo=3, firstName=null, middleName=null, lastName=null, nameCn=null, orcid=null, stid=null, country=null, authorPic=null, dead=0, email=null, emailSecond=null, emailThird=null, correspondingAuthor=0, authorType=1, ext={EN=AuthorExt(id=1325132995258098178, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, authorId=1325132994595398144, language=EN, stringName=Qiang WANG, firstName=Qiang, middleName=null, lastName=WANG, prefix=null, suffix=null, authorComment=null, nameInitials=null, affiliation=null, department=null, xref=null, address=Institute of Fluid Physics, CAEP, Mianyang 621999, Sichuan, China, bio=null, bioImg=null, bioContent=null, aboutCorrespAuthor=null), CN=AuthorExt(id=1325132997338472963, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, authorId=1325132994595398144, language=CN, stringName=王强, firstName=null, middleName=null, lastName=null, prefix=null, suffix=null, authorComment=null, nameInitials=null, affiliation=null, department=null, xref=null, address=中国工程物理研究院流体物理研究所, 四川绵阳　621999, bio=null, bioImg=null, bioContent=null, aboutCorrespAuthor=null)}, companyList=[AuthorCompany(id=1325132987817402860, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, xref=null, ext=[AuthorCompanyExt(id=1325132987855151597, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, companyId=1325132987817402860, language=EN, country=null, province=null, city=null, postcode=null, companyName=null, departmentName=null, remark=Institute of Fluid Physics, CAEP, Mianyang 621999, Sichuan, China), AuthorCompanyExt(id=1325132987980980718, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, companyId=1325132987817402860, language=CN, country=null, province=null, city=null, postcode=null, companyName=null, departmentName=null, remark=中国工程物理研究院流体物理研究所, 四川绵阳　621999)])]), Author(id=1325132998013755910, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, orderNo=4, firstName=null, middleName=null, lastName=null, nameCn=null, orcid=null, stid=null, country=null, authorPic=null, dead=0, email=null, emailSecond=null, emailThird=null, correspondingAuthor=0, authorType=1, ext={EN=AuthorExt(id=1325132998575792648, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, authorId=1325132998013755910, language=EN, stringName=Jiupeng TANG, firstName=Jiupeng, middleName=null, lastName=TANG, prefix=null, suffix=null, authorComment=null, nameInitials=null, affiliation=null, department=null, xref=null, address=Institute of Fluid Physics, CAEP, Mianyang 621999, Sichuan, China, bio=null, bioImg=null, bioContent=null, aboutCorrespAuthor=null), CN=AuthorExt(id=1325132998986834441, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, authorId=1325132998013755910, language=CN, stringName=唐久棚, firstName=null, middleName=null, lastName=null, prefix=null, suffix=null, authorComment=null, nameInitials=null, affiliation=null, department=null, xref=null, address=中国工程物理研究院流体物理研究所, 四川绵阳　621999, bio=null, bioImg=null, bioContent=null, aboutCorrespAuthor=null)}, companyList=[AuthorCompany(id=1325132987817402860, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, xref=null, ext=[AuthorCompanyExt(id=1325132987855151597, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, companyId=1325132987817402860, language=EN, country=null, province=null, city=null, postcode=null, companyName=null, departmentName=null, remark=Institute of Fluid Physics, CAEP, Mianyang 621999, Sichuan, China), AuthorCompanyExt(id=1325132987980980718, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, companyId=1325132987817402860, language=CN, country=null, province=null, city=null, postcode=null, companyName=null, departmentName=null, remark=中国工程物理研究院流体物理研究所, 四川绵阳　621999)])]), Author(id=1325132999372710411, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, orderNo=5, firstName=null, middleName=null, lastName=null, nameCn=null, orcid=null, stid=null, country=null, authorPic=null, dead=0, email=null, emailSecond=null, emailThird=null, correspondingAuthor=0, authorType=1, ext={EN=AuthorExt(id=1325132999834083853, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, authorId=1325132999372710411, language=EN, stringName=Haibin OU, firstName=Haibin, middleName=null, lastName=OU, prefix=null, suffix=null, authorComment=null, nameInitials=null, affiliation=null, department=null, xref=null, address=Institute of Fluid Physics, CAEP, Mianyang 621999, Sichuan, China, bio=null, bioImg=null, bioContent=null, aboutCorrespAuthor=null), CN=AuthorExt(id=1325133001830572560, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, authorId=1325132999372710411, language=CN, stringName=欧海彬, firstName=null, middleName=null, lastName=null, prefix=null, suffix=null, authorComment=null, nameInitials=null, affiliation=null, department=null, xref=null, address=中国工程物理研究院流体物理研究所, 四川绵阳　621999, bio=null, bioImg=null, bioContent=null, aboutCorrespAuthor=null)}, companyList=[AuthorCompany(id=1325132987817402860, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, xref=null, ext=[AuthorCompanyExt(id=1325132987855151597, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, companyId=1325132987817402860, language=EN, country=null, province=null, city=null, postcode=null, companyName=null, departmentName=null, remark=Institute of Fluid Physics, CAEP, Mianyang 621999, Sichuan, China), AuthorCompanyExt(id=1325132987980980718, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, companyId=1325132987817402860, language=CN, country=null, province=null, city=null, postcode=null, companyName=null, departmentName=null, remark=中国工程物理研究院流体物理研究所, 四川绵阳　621999)])]), Author(id=1325133002279363090, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, orderNo=6, firstName=null, middleName=null, lastName=null, nameCn=null, orcid=null, stid=null, country=null, authorPic=null, dead=0, email=null, emailSecond=null, emailThird=null, correspondingAuthor=0, authorType=1, ext={EN=AuthorExt(id=1325133002744930836, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, authorId=1325133002279363090, language=EN, stringName=Mingxian KAN, firstName=Mingxian, middleName=null, lastName=KAN, prefix=null, suffix=null, authorComment=null, nameInitials=null, affiliation=null, department=null, xref=null, address=Institute of Fluid Physics, CAEP, Mianyang 621999, Sichuan, China, bio=null, bioImg=null, bioContent=null, aboutCorrespAuthor=null), CN=AuthorExt(id=1325133003248247318, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, authorId=1325133002279363090, language=CN, stringName=阚明先, firstName=null, middleName=null, lastName=null, prefix=null, suffix=null, authorComment=null, nameInitials=null, affiliation=null, department=null, xref=null, address=中国工程物理研究院流体物理研究所, 四川绵阳　621999, bio=null, bioImg=null, bioContent=null, aboutCorrespAuthor=null)}, companyList=[AuthorCompany(id=1325132987817402860, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, xref=null, ext=[AuthorCompanyExt(id=1325132987855151597, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, companyId=1325132987817402860, language=EN, country=null, province=null, city=null, postcode=null, companyName=null, departmentName=null, remark=Institute of Fluid Physics, CAEP, Mianyang 621999, Sichuan, China), AuthorCompanyExt(id=1325132987980980718, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, companyId=1325132987817402860, language=CN, country=null, province=null, city=null, postcode=null, companyName=null, departmentName=null, remark=中国工程物理研究院流体物理研究所, 四川绵阳　621999)])]), Author(id=1325133004162605592, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, orderNo=7, firstName=null, middleName=null, lastName=null, nameCn=null, orcid=null, stid=null, country=null, authorPic=null, dead=0, email=null, emailSecond=null, emailThird=null, correspondingAuthor=0, authorType=1, ext={EN=AuthorExt(id=1325133006196843034, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, authorId=1325133004162605592, language=EN, stringName=Shuchao DUAN, firstName=Shuchao, middleName=null, lastName=DUAN, prefix=null, suffix=null, authorComment=null, nameInitials=null, affiliation=null, department=null, xref=null, address=Institute of Fluid Physics, CAEP, Mianyang 621999, Sichuan, China, bio=null, bioImg=null, bioContent=null, aboutCorrespAuthor=null), CN=AuthorExt(id=1325133006851154459, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, authorId=1325133004162605592, language=CN, stringName=段书超, firstName=null, middleName=null, lastName=null, prefix=null, suffix=null, authorComment=null, nameInitials=null, affiliation=null, department=null, xref=null, address=中国工程物理研究院流体物理研究所, 四川绵阳　621999, bio=null, bioImg=null, bioContent=null, aboutCorrespAuthor=null)}, companyList=[AuthorCompany(id=1325132987817402860, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, xref=null, ext=[AuthorCompanyExt(id=1325132987855151597, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, companyId=1325132987817402860, language=EN, country=null, province=null, city=null, postcode=null, companyName=null, departmentName=null, remark=Institute of Fluid Physics, CAEP, Mianyang 621999, Sichuan, China), AuthorCompanyExt(id=1325132987980980718, tenantId=1146029695717560320, journalId=1255847867265597444, articleId=1325132959933669822, companyId=1325132987817402860, language=CN, country=null, province=null, city=null, postcode=null, companyName=null, departmentName=null, remark=中国工程物理研究院流体物理研究所, 四川绵阳　621999)])])]

王刚华, 谢龙, 肖波, 王强, 唐久棚, 欧海彬, 阚明先, 段书超

作者信息

中国工程物理研究院流体物理研究所, 四川绵阳　621999

作者简介:

王刚华（1976－），男，博士，副研究员，主要从事电磁驱动多物理场耦合理论与数值方法研究. E-mail：wanggh@caep.cn

Electric Explosion Early Process Analysis of Metal Bridge Foil Based on an Electromagnetic-Thermal-Mechanical Model

Ganghua WANG, Long XIE, Bo XIAO, Qiang WANG, Jiupeng TANG, Haibin OU, Mingxian KAN, Shuchao DUAN

Affiliations

Institute of Fluid Physics, CAEP, Mianyang 621999, Sichuan, China

出版时间: 2024-02-05 doi: 10.11858/gywlxb.20230711

文章导航

摘要

收起

金属桥箔的电爆炸过程对电炮、冲击片雷管等的性能影响极为关键。这一过程中物质性质变化、几何构型和动力学过程非常复杂，早期的大量理论模拟工作均采用极为简化的物理模型。为此，建立了描述爆炸箔早期行为的三维电磁-热-固体力学耦合求解全物理模型，模拟爆炸箔在大电流加载下的早期受热膨胀过程，并对早期膨胀过程中桥箔上的磁场、电流、温度等演化进行分析，观察到电流在构型和电阻率两种因素影响下的角扩散和线扩散现象，模拟得到的桥区温度场分布与参考文献中的实验图像及模拟结果定性符合。

关键词

爆炸箔 / 冲击片雷管 / 磁流体 / 多场耦合

Abstract

收起

The electrical explosion of metal bridge foil is a key physical process in many fields, e.g. electric guns and slapper detonator. Due to the complexity of material properties, dynamic processes, and effects of geometric configurations, a large amount of earlier theoretical simulations used simplified physical models. In this paper, a three-dimensional electromagnetic-thermal-mechanical model was established. which was a full physical model describing the early behavior of explosive foil. The early thermal expansion process of the explosive foil under large current loading was simulated, and the evolution of magnetic field, current, temperature, expansion speed of the foil was analyzed. The phenomena of angular and linear current diffusions induced by the configuration and the resistivity were observed in the simulation. The temperature distribution in the bridge region obtained from the simulation is qualitatively consistent with experiment and simulation results of literature.

Key words

explosive foil / slapper detonator / magnetic fluid / multi-field coupling

引用本文

王刚华, 谢龙, 肖波, 王强, 唐久棚, 欧海彬, 阚明先, 段书超. 电爆炸金属桥箔早期过程中电磁-热-力多物理场耦合建模与分析. 高压物理学报, 2024 , 38 (1) : 012301-1 -012301-9 . DOI: 10.11858/gywlxb.20230711

Ganghua WANG, Long XIE, Bo XIAO, Qiang WANG, Jiupeng TANG, Haibin OU, Mingxian KAN, Shuchao DUAN. Electric Explosion Early Process Analysis of Metal Bridge Foil Based on an Electromagnetic-Thermal-Mechanical Model[J]. Chinese Journal of High Pressure Physics, 2024 , 38 (1) : 012301-1 -012301-9 . DOI: 10.11858/gywlxb.20230711

正文

收起

利用大电流流经金属箔/金属丝时的欧姆加热作用，可以将金属箔/丝汽化电离，产生急速膨胀的高能量密度等离子体，进而可以作为高压驱动源。该类高压驱动源应用广泛，例如：基于该原理发展的电炮可以用于开展超高速飞片发射和材料动力学特性研究^[1–4]；利用金属箔电爆炸发射小尺寸飞片来制造高安全性、低感度的冲击起爆雷管^[5–10]等。

电爆炸过程中，导体被快速加热，短时间内经历固–液–气/等离子体多个状态，欧姆电阻可能增长几个数量级。因此，采用数值模拟方法研究这一过程并非易事。早期的研究主要基于电路方程和电导率模型以实现对实验电流的模拟和温度的估算^[11–12]。1981年，董玉斌等^[12]将金属箔电爆炸过程分为沸点前和沸点后两个阶段，分阶段计算了放电电流和桥箔电阻的变化以及金属箔电爆炸驱动飞片的加速过程，实现了一维流体动力学方程、放电回路方程以及沙哈方程的耦合求解。1989年，Osher等^[13]报道了用于计算金属箔电爆炸全过程的一维磁流体动力学程序，该程序与2010年Saxena等^[14]报道的基于Burgess电导率模型的一维磁流体程序SEEFA一样，均对金属箔区域做零维处理。后来，罗斌强^[15]、贺佳等^[16]和Luo等^[17]结合大量实验完善了电箔爆炸一维磁流体程序，该程序能较好地模拟爆炸箔的放电波形，同时给出与实验基本一致的飞片速度波形。

然而，对于冲击片雷管而言，提升雷管工作的可靠性需要对桥箔电爆炸驱动飞片的过程有更为深入的认识，一维简化的磁流体模型虽然能够预估飞片的速度，但不能给出飞片的形状、温度、速度的空间分布这些可靠起爆炸药的关键信息，也无法给出爆炸箔上电流密度、磁场的空间分布，因而无法为具体的改进设计提供参考。因此，发展更全面的物理模型和更高维度的模拟程序是深入理解冲击片雷管爆炸物理过程的基础。2017年，英国原子武器研究中心（Atomic Weapons Establishment，AWE）和美国圣地亚实验室（Sandia National Laboratories，SNL）报道了他们利用三维磁流体模拟程序ALEGRA MHD开展的爆炸箔起爆的数值仿真工作^[18]，首次获得了与实验照相结果基本相似的图像，是目前报道的最先进的冲击片雷管模拟能力。

近年来，我国的爆炸箔数值仿真工作取得了不少进展，在工程应用中发挥了较大作用^[19–23]。然而，目前国内相关数值模拟工作还存在物理模型不全、空间维度降阶等问题，对于最为关键的桥箔三维结构和电磁场分布造成的等离子体空间分布等问题缺乏研究模型和研究工具。例如，北京理工大学与中国工程物理研究院化工材料研究所合作开展了三维冲击片雷管的数值模拟研究，对桥箔进行了简化处理，使用的是等效电阻模型，无法分析电磁场和电流的三维空间分布。虽然基于磁流体的数值模型在物理上较为全面，但是现有的报道还是一维的，中国工程物理研究院流体物理研究所在二维磁流体方面的研究正在进行，发展二维、三维磁流体数值模拟能力尚需时日。

尽管桥箔电爆炸的全过程模拟极其困难，但是其后续复杂的物理过程是在早期能量沉积和分布的基础上演化而来的，因此，在大量等离子体形成和爆炸膨胀之前，弄清楚桥箔上的电流分布和温度分布演化过程极其重要，有利于开展桥箔性能改进设计。本研究综合考虑电磁场在三维结构上的扩散、欧姆加热效应、材料受热膨胀效应，拟建立三维电爆炸桥箔的电磁-热-固体力学耦合模型，开展桥箔早期电磁能量沉积过程模拟。这一模型充分考虑爆炸桥箔发生剧烈汽化膨胀前的各种物理效应，是一种宏观尺度上的全三维、全物理模拟方法，对清晰理解爆炸桥箔早期过程和指导工程设计具有重要意义。

1 电磁-热-固体力学耦合模型

收起

桥箔的早期能量沉积过程取决于电流密度和电导率，而电导率又依赖于温度场分布。由于早期温度相对较低，桥箔膨胀较小，可以使用固体力学来描述桥箔的膨胀过程。为此，建立了电磁-热-固体力学耦合的物理模型以模拟桥箔早期的能量沉积过程。采用电磁模块进行桥箔的电磁和电流演化分析，采用固体传热模块考虑桥箔的温度变化，采用固体力学模块模拟桥箔的早期受热膨胀情况。

电磁模块求解如下控制方程。

对于桥箔区域，使用A形式的矢量势公式（安培定律）描述

$ \nabla \times {\boldsymbol{H}} = {\boldsymbol{J}} $

$ {\boldsymbol{B}} = \nabla \times {\boldsymbol{A}} $

$ \boldsymbol{E}=-\frac{\partial\boldsymbol{A}}{\partial\mathit{{t}}} $

$ {\boldsymbol{J}} = \sigma {\boldsymbol{E}} $

对于空气区域，使用标量势公式（磁通量守恒）描述

$ \left\{ \begin{array}{l}\nabla {{{\cdot}} {\boldsymbol{B}}}=0 \\ \boldsymbol{B}=\mu_{\text{0}}\mu_{\text{r}}\boldsymbol{H} \\ \boldsymbol{H}=-\nabla V_{\mathrm{m}} \\ \end{array}\right. $

式中：H为磁场强度，A/m；B为磁通密度，T；E为电场，V/m；A为磁矢势，Wb/m；

$ \sigma $

为电导率，S/m；

$ {\mu _0} $

为真空磁导率，H/m；

$ {\mu _{\text{r}}} $

为相对磁导率；

$ {V_{\text{m}}} $

为磁标势，A。

固体传热模块求解如下控制方程

$ \rho c_p\frac{\partial T}{\partial t}+\rho c_p{\boldsymbol{u}}{\cdot}\nabla T=\nabla\mathbf{\cdot}\left(k\nabla T\right)+Q\mathrm{_e} $

式中：

$ \rho $

为密度，kg/m³；

$ c_p $

为比定压热容，J/(kg∙K)；

$ T $

为温度，K；

$ \boldsymbol{u} $

为速度，m/s；

$ k $

为热导率，W/(m∙K)；

$ {Q_{\text{e}}} $

为热源，W/m³，耦合计算中热源为欧姆加热，即

$ Q_{\text{e}}={{{\boldsymbol{J}}\cdot {\boldsymbol{E}}}} $

。

固体力学模块求解如下控制方程

$ \rho\frac{\partial^2s}{\partial t^2}=\nabla{\cdot{\boldsymbol{S}}}+F\mathrm{_v} $

式中：s为位移，m；

$ \boldsymbol{S} $

为应力张量，Pa；

$ F_{\mathrm{v}} $

为外力，N。

耦合计算中调用热膨胀模块，即考虑由热效应引起的热应变

$ \boldsymbol{\varepsilon}_{\text{th}}=\alpha\left(T\right)\left(T-T_{\text{ref}}\right) $

$ \boldsymbol{S} = \boldsymbol{C}:{\boldsymbol{\varepsilon} _{{\text{th}}}} $

式中：

$ \alpha \left( T \right) $

为随温度变化的热膨胀系数，K⁻¹；

$ \boldsymbol{\varepsilon}_{\text{th}} $

为弹性应变；

$ \boldsymbol{C} $

为弹性模量张量，Pa；

$ {T_{{\text{ref}}}} $

为参考温度，K。

计算初始时刻，电场、磁场强度、电流密度等均为零，温度为常温。电磁边界条件主要包括：对称面上指定磁标势

$ {V_0} $

，区域边界上法相磁场为零，空气与桥箔交界面上切向磁场连续。固体传热模块边界条件：将所有面均设定为绝热条件。力学边界条件主要包括：对称面上位移为零，其余面为自由面。计算中将图1所示的实验电流作为激励电流，对应的回路电阻为280

$ {\text{m}}\Omega $

，回路电感为310 nH，电容为0.25 μF，充电电压为5.1 kV。桥箔材料为铜，考虑了铜材料的热导率（

$ k $

）随温度（T）和密度（

$ \rho $

）的变化以及电导率随温度和密度的变化情况（如图2和图3所示）。

具体求解过程中，3个模块的耦合关系如下：电磁求解模块通过电流密度J将电流分布传递给热传导模块，传热模块通过温度T与固体力学模块进行耦合；同时，材料温度的改变通过电导率模型传递给电磁模块，用于计算实时的电磁场分布；固体力学计算产生的几何变形，也会耦合传递到电磁模块和传热模块。

为了减少计算量，采用如图4所示的1/2几何模型进行分析，并通过镜像还原整体情况（图4中未显示外围空气区）。由于物理模型只适用于小变形，因此计算时使用拉氏网格，网格划分采用自适应四面体网格，在所关心的桥箔区采用密网格，空气区使用较粗网格，总网格数约为2.084×10⁷，如图5所示。图4中，由于桥箔的厚度（z方向）很薄，仅为5 μm（为了便于显示，图4中将厚度放大了10倍），远小于x、y方向上的尺寸（毫米量级），为了减小计算量，划分网格时，铜箔区在z方向只划分了2层网格。为便于后续沿标记线进行物理量分布分析，在如图4所示的桥面上标记两条线L1、L2，其中：L1为桥面上对称轴线，L2为边缘轮廓线，弧长为l。

2 计算分析

收起

2.1 基于常数电导率的简要分析

采用常温电导率模型分析金属桥箔的磁场和电流密度分布，不考虑材料电导率、热导率随温度和密度的变化。由于常温电导率模型只能用于分析放电早期的电磁场，为了更接近实际情况，需要进一步考虑欧姆加热的影响，因此，只对常数电导率模型的计算结果作简要介绍。

鉴于材料性质不发生变化，可采用归一化电流进行分析。图6给出了桥箔区附近的磁场强度分布。受几何构型的影响，磁场强度在桥区和过渡区的交界处表现出明显的非均匀性，特别是边缘附近非均匀性更为显著。因此，边缘附近的电流密度也很大，如图7所示。

依据图7的电流密度分布，绘制电流密度的高度图，如图8所示。图9给出了沿L1、L2的电流密度随弧长l的分布，弧长起点为桥箔左端。从图8和图9可以非常直观地看出：随着过渡区变窄，电流密度逐步增加，在桥区和过渡区交界处的4个角点上电流最为集中；在桥区主体部分虽然还存在一定的边界效应，但桥面主体部分的电流分布较为均匀。

爆炸箔上的电流分布决定了各处产生的电功率密度P_e（图10），也随即决定了各处的温升和早期等离子体的分布，因为在电功率高的地方，材料被率先加热，进而汽化形成等离子体。

本节的所有分析是基于常数电导率模型开展的，即爆炸箔的电导率在空间上均匀分布，因此，电流和磁场只与爆炸箔的几何尺寸相关，各个时刻获得的磁场强度、电流密度的空间分布类似，不同时刻电流的大小仅影响其绝对值。然而，实际上，金属的电导率与其自身的温度和密度密切相关，需要考虑这种因素的影响。

2.2 基于Burgess电导率的计算分析

常数电导率模型只能分析放电早期的电磁场，在中后期，为了更接近实际情况，需要进一步考虑欧姆加热的影响。事实上，随着电流加热形成等离子体，材料的电参数会发生复杂的变化，为此，在电导率模型中引入随温度变化的Burgess电导率模型，并引入固体传热模块和固体力学模块，以分析桥箔放电一段时间内的电磁场分布情况。

图11给出了引入Burgess电导率模型、固体传热模块和固体力学模块后计算得到的不同时刻爆炸箔上的磁场强度分布。可以看出，与常数电导率模型类似，改进后的模型计算结果显示，磁场先集中分布在爆炸箔4个角点和爆炸箔边缘，这些地方的磁场梯度较大，磁场分布极不均匀。在桥区，由于宽度变窄，磁场相对较大，边缘效应更加明显，特别是桥区与过渡区衔接的4个角点的磁场强度很大。与常数电导率模型计算结果不同的是，随着时间的增长，磁场向爆炸箔中间扩散，中间区域的磁场逐步增大。磁场扩散呈现两种扩散特征：一是从4个角点出发的点扩散，二是从边缘向里的线扩散。

电流分布主要由两个因素决定：构型和电导率。前述的磁场角点和边缘效应是构型引起的。由于电流密度对应磁场强度的旋度，因此，磁场不均匀的地方电流密度更大。图12给出了由改进后的模型计算获得的不同时刻的电流密度分布。可以看出：电流一开始集中在边缘，然后向中间扩散，此时构型是电流分布的主导因素，在70 ns时，桥区主体的电流分布基本均匀；此后，电导率因素占据主导，边缘由于温升快，其电导率低于中心处的电导率，使得电流向中心集中；后续电流分布在两种因素的竞争中不断演化。图13给出了90 ns时刻沿线L1、L2的电流密度分布，可以看出，其大致特征与常数电导模型率给出的结果一致，只是桥区的电流分布更均匀。

爆炸箔的主要温升机制是欧姆加热，温度分布取决于电流密度、电导率和时间。图14给出了温度分布的演化。可以看出，早期边缘加热效应明显，4个角点的温升相比于其他地方更快，温度也更高。但是，随着磁场和电流的扩散作用，电流向温度较低的中间区域集中，桥区主体部分的温度分布逐渐变得匀滑。图15给出了沿线L1、L2的温度分布。两条曲线的平台区对应桥区，两条曲线的平台有错位，这是L1、L2的弧长不一致导致的。两条曲线的平台有高度差，反映了桥区边缘与桥区中心有温度差。92 ns时，桥区大部分区域的温度均在1245～1345 K之间，温度差在100 K以内，这意味着在该结构设计下爆炸箔的温度分布已十分均匀。通常情况下，在设计冲击片雷管时，将桥箔爆炸时刻设计在电流峰值附近，虽然计算时刻与电流峰值时刻相距时间较长，尚不能确定后续温度分布是否更加均匀，但是从目前的计算结果来看，温度均匀性已经较好。此外，爆炸箔紧邻桥区的过渡区中有相当多的区域温升较为明显，意味着该处沉积了较多能量。在冲击片雷管设计中，可以充分利用该处产生的等离子体能量来驱动飞片，使飞片获得更多动能。因此，从提高能量利用率的角度来讲，设计冲击片雷管时可考虑使炮筒口径略大于桥区。

图16给出了90 ns时爆炸箔沿z轴的膨胀速度分布，4个角点附近的膨胀速度相对较快，桥区边缘的膨胀速度比桥区中心慢，其原因是边缘处金属箔除了z向膨胀以外，还向宽度方向膨胀。

图17显示了采用改进模型计算的温度分布与文献[16]报道结果的对比，其中：右图为3个时刻的温度分布灰度图，左图为文献[16]中ALEGRA MHD3D的模拟结果，中间图为文献[16]给出的实验照片。可以看出，文献[16]的模拟结果没有抓拍到实验观测到的边缘加热效应，而本研究的模拟结果较为清晰地显示了前期的边缘加热效应和持续时间较长的角加热效应，4个角点附近的扩散现象也与实验结果定性相符。

以上分析只适应于电爆炸箔的初步设计及其工作过程的早期情况。在角点出现明显的受热膨胀以后，建立的模型受制于固体力学模型所能描述的范围，无法进行后续过程的模拟。事实上，随着电流加热形成等离子体，材料的电参数会发生复杂的变化，其动力学行为将是复杂的磁流体力学过程，需要开展更为深入的磁流体力学建模研究。

3 结　论

收起

建立了三维电磁-热-固体力学耦合求解模型，模拟了爆炸箔在大电流加载下的早期受热膨胀过程。基于常数电导率模型，给出了爆炸箔的磁场、电流密度、电功率密度的空间分布。基于Burgess电导率模型，分析了早期膨胀过程中桥箔的磁场、电流密度、温度的演化情况，指出影响其分布的决定因素为构型和桥箔电导率分布。模拟了爆炸箔的膨胀速度，计算结果较好地捕捉到了实验中观测到的边缘加热效应和角加热效应。

总体而言，建立的三维电磁-热-力耦合模型可有效地模拟爆炸箔早期行为。本研究初步分析了金属桥箔结构和电导率对电流、磁场和温度分布的影响，研究结果是不同构型、不同材料金属桥箔物理设计的基础。后续将会进一步总结电爆炸过程的物理规律。

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

GUENTHER A H, WUNSCH D C, SOAPES T D, et al. Acceleration of thin plates by exploding foil techniques [M]//CHACE W G, MOORE H K. Exploding Wires. Boston: Springer, 1962: 279–298.

KELLER D V, PENNING J R JR. Exploding foil-the production of plane shock waves and the acceleration of thin plates [M]//CHACE W G, MOORE H K. Exploding Wires. Boston: Springer, 1962: 263–277.

罗斌强, 张旭平, 郝龙, 等. 7 km/s以上超高速发射技术研究进展 [J]. 爆炸与冲击, 2021, 41(2): 021401.

LUO B Q, ZHANG X P, HAO L, et al. Advances on the techniques of ultrahigh-velocity launch above 7 km/s [J]. Explosion and Shock Waves, 2021, 41(2): 021401.

WANG G J, HE J, ZHAO J H, et al. The techniques of metallic foil electrically exploding driving hypervelocity flyer to more than 10 km/s for shock wave physics experiments [J]. Review of Scientific Instruments, 2011, 82(9): 095105.

CUTTING J L, LEE R S, HONGIN R I. Slapper detonator perfonnance: UCRL-IDl 17787 [R]. USA: Lawrence Livermore National Labratory, 1994.

王桂吉, 蒋吉昊, 邓向阳, 等. 电爆炸驱动小尺寸冲击片实验与数值计算研究 [J]. 兵工学报, 2008, 29(6): 657–661.

WANG G J, JIANG J H, DENG X Y, et al. Experiments and numerical simulation of small-scale slapper driven by electrical explosion [J]. Acta Armamentarii, 2008, 29(6): 657–661.

张玉若, 梁车平, 程涛, 等. 爆炸箔换能元失效模式分析 [J]. 火工品, 2018(4): 5–7.

ZHANG Y R, LIANG C P, CHENG T, et al. Experimental study on failure mode of exploding foil [J]. Initiators & Pyrotechnics, 2018(4): 5–7.

秦国圣, 张蕊, 王寅, 等. 多点阵列冲击片雷管仿真与试验研究 [J]. 兵工学报, 2016, 37(Suppl 2): 81–85.

QIN G S, ZHANG R, WANG Y, et al. Numerical simulation and experimental study of multi-point array EFI [J]. Acta Armamentarii, 2016, 37(Suppl 2): 81–85.

先明春, 谢浚尧, 王成玲, 等. 考虑压缩空气边界的爆炸箔起爆器飞片运动模型 [J]. 爆破器材, 2023, 52(2): 1–7.

XIAN M C, XIE J Y, WANG C L, et al. A Motion model of flyer in exploding foil initiator considering compressed air boundary [J]. Explosive Materials, 2023, 52(2): 1–7.

侯新朋, 彭志凌, 宋进宇, 等. 桥区结构对爆炸箔起爆器发火性能的影响 [J]. 兵器装备工程学报, 2022, 43(5): 103–107, 113.

HOU X P, PENG Z L, SONG J Y, et al. Influence of bridge structure on ignition performance of exploding foil detonator [J]. Journal of Ordnance Equipment Engineering, 2022, 43(5): 103–107, 113.

LOGAN J D, LEE R S, WEINGART R C, et al. Calculation of heating and burst phenomena in electrically exploded foils [J]. Journal of Applied Physics, 1977, 48(2): 621–628.

董玉斌, 李献文, 董维申. 电加热金属箔爆炸的有限元法计算 [J]. 爆炸与冲击, 1981, 1(1): 28–36.

DONG Y B, LI X W, DONG W S. A finite-element method calculation of electrically heated bursting metallic foils [J]. Explosion and Shock Waves, 1981, 1(1): 28–36.

OSHER J E, BARNES G, CHAU H H, et al. Operating characteristics and modeling of the LLNL 100-kV electric gun [J]. IEEE Transactions on Plasma Science, 1989, 17(3): 392–402.

SAXENA A K, KAUSHIK T C, GUPTA S C. Shock experiments and numerical simulations on low energy portable electrically exploding foil accelerators [J]. Review of Scientific Instruments, 2010, 81(3): 033508.

罗斌强. 金属箔电爆炸及其在冲击动力学中的应用 [D]. 合肥: 中国科学技术大学, 2012.

LUO B Q. Electrical explosion of metallic foils and its applications in dynamic mechanics [D]. Hefei: University of Science and Technology of China, 2012.

贺佳, 罗斌强, 庞树财, 等. 微型爆炸箔电爆炸过程的数值模拟 [J]. 高压物理学报, 2017, 31(1): 21–26.

HE J, LUO B Q, PANG S C, et al. Numerical simulation of electrical explosion of micro-exploding foil [J]. Chinese Journal of High Pressure Physics, 2017, 31(1): 21–26.

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.

NEAL W, GARASI C. High fidelity studies of exploding foil initiator bridges, part 3: ALEGRA MHD simulations [J]. AIP Conference Proceedings, 2017, 1793: 080008.

苏谦, 褚恩义, 解瑞珍, 等. 九点爆炸箔电爆炸性能试验与仿真模拟 [J]. 火工品, 2019(5): 15–18.

SU Q, CHU E Y, XIE R Z, et al. Experimental study and simulation on the electrical explosion performance of nine-point explosive foil [J]. Initiators & Pyrotechnics, 2019(5): 15–18.

钱石川, 甘强, 任志伟, 等. 爆炸箔起爆器发火阈值影响因素的数值模拟 [J]. 含能材料, 2018, 26(3): 248–254.

QIAN S C, GAN Q, REN Z W, et al. Numerical simulation of the factors affecting the ignition threshold of an exploding foil initiator [J]. Chinese Journal of Energetic Materials, 2018, 26(3): 248–254.

伍俊英, 于红新, 汪龙, 等. 金属桥箔水中电爆炸流场数值模拟研究 [J]. 兵工学报, 2016, 37(Suppl 1): 51–56.

WU J Y, YU H X, WANG L, et al. Numerical simulation of electric exploding of metal bridge foil in water [J]. Acta Armamentarii, 2016, 37(Suppl 1): 51–56.

王亮, 邹苑楠, 蒋小华, 等. 短脉冲电流作用下铜微桥箔的电热分析 [J]. 含能材料, 2013, 21(4): 500–505.

WANG L, ZOU Y N, JIANG X H, et al. Thermal-electric analysis of small-scale copper bridge foils excited by short pulse currents [J]. Chinese Journal of Energetic Materials, 2013, 21(4): 500–505.

曹始发, 黄寅生. 小飞片冲击起爆HNS-IV的临界速度的研究 [J]. 计算机仿真, 2016, 33(4): 23–26, 149.

CAO S F, HUANG Y S. Research on threshold velocity of shock initiation of HNS-IV by small flyer [J]. Computer Simulation, 2016, 33(4): 23–26, 149.

2024年第38卷第1期

PDF下载

引用本文

BibTeX

文章信息

doi: 10.11858/gywlxb.20230711

接收时间：2023-08-11
首发时间：2026-11-05
出版时间：2024-02-05

补充材料

相关文章

文章信息

作者

出版历史

收稿日期：2023-08-11
修回日期：2023-10-30

基金

作者信息

中国工程物理研究院流体物理研究所, 四川绵阳　621999

参考文献

分享链接

https://castjournals.cast.org.cn/joweb/gywlxb/CN/10.11858/gywlxb.20230711

分享至

全文二维码

扫描看全文

引用本文

BibTeX

本文的引用情况

2种不同金属材料的力学参数

科 Family	属数 Number of genus	种数 Number of species	占总种数比例 Percentage of total species (%)	属 Genus	种数 Number of species	占总种数比例 Percentage of total species (%)
鹅膏菌科Amanitaceae	2	11	5.26	鹅膏菌属 Amanita	10	4.78
小菇科 Mycenaceae	2	12	5.74	丝盖伞属 Inocybe	5	2.39
多孔菌科 Polyporaceae	8	14	6.70	蜡蘑属 Laccaria	5	2.39
红菇科 Russulaceae	3	23	11.00	小皮伞属 Marasmius	6	2.87
				小菇属 Mycena	11	5.26
				光柄菇属 Pluteus	5	2.39
				红菇属 Russula	17	8.13
				栓菌属 Trametes	5	2.39

关闭全屏

BibTeX
EndNote
RefWorks
TxT