基于SPH方法的变压器网侧套管内部爆炸数值模拟及稳定性分析

赵春风; 周志航; 朱玉富

doi:10.11858/gywlxb.20220674

基于SPH方法的变压器网侧套管内部爆炸数值模拟及稳定性分析

doi: 10.11858/gywlxb.20220674

合肥工业大学土木与水利工程学院, 安徽合肥　230009

基金项目: 安徽能源互联网联合基金（2008085UD12）

详细信息

作者简介:
赵春风（1983-），男，博士，教授，主要从事工程结构抗震与减震、组合结构抗爆研究.E-mail：zhaowindy@hfut.edu.cn

中图分类号: O382
计量
- 文章访问数: 218
- HTML全文浏览量: 121
- PDF下载量: 32
出版历程
- 收稿日期: 2022-10-11
- 修回日期: 2022-11-04
- 网络出版日期: 2023-04-26
- 刊出日期: 2023-06-05

Numerical Simulation and Stability Analysis of Internal Explosion of Transformer Net Side Bushing by SPH Method

School of Civil Engineering, Hefei University of Technology, Hefei 230009, Anhui, China

摘要

摘要: 网侧套管是变压器的重要组成部分，在使用过程中因绝缘击穿现象可能引起套管内部冷却油爆炸，给变压器箱体造成很大的安全隐患，因此开展网侧套管爆炸事故的定量评估具有重要意义。通过非线性有限元软件ANSYS/LS-DYNA建立了二维变压器套管模型，采用光滑粒子流体动力学法对变压器套管在内部爆炸作用下的动态响应进行了模拟，分析了不同参数对套管破坏特征的影响规律。通过基于套管径向粒子速度曲线的稳定性判断方法，评估了各套管的失稳时间。结果表明：在内部爆炸作用下，套管管壁中部在内外壁拉压联合作用下率先产生破坏，在冲击波传播过程中套管整体损伤呈凸状变化趋势。过高的爆炸当量、冷却油的存在和初始裂纹缺陷对套管保持稳定有较大影响。爆炸当量的减小可以使套管破坏模式由双向剪切破坏向受拉破坏转变，整体稳定性也随之增强。当爆源位于引线外壁时，冷却油流体所辐射出的冲击波不仅使套管失稳时间有所提前，还会造成外壁膨胀破坏范围变大。应力集中现象和有效壁厚的减小使含初始裂纹缺陷套管剪切破坏失稳现象发展迅速。
- 变压器套管 /
- 光滑粒子流体动力学法 /
- 失稳时间 /
- 破坏模式
Abstract: The net side bushing is an important part of the transformer. During operation, the cooling oil inside the bushing will explode due to insulation breakdown, which will cause great potential safety hazard to the transformer. Therefore, it is of great significance to carry out quantitative evaluation of the net side bushing explosion accident. In this paper, a two-dimensional transformer bushing model was established by nonlinear finite element software ANSYS/LS-DYNA with smooth particle hydrodynamics (SPH) method. The dynamic response of transformer bushing under internal explosion was simulated. And the influence of different parameters on the failure characteristics of bushing was analyzed. The instability time of bushing was evaluated by the stability judgment method based on the curves of radial particle velocity. The results show that under the action of internal explosion, the middle part of the bushing is first damaged under the combined action of tension and compression of the inner and outer surfaces, and the overall damage of the bushing shows a convex trend during the propagation of the shock wave. The high explosion equivalent, the existence of cooling oil and initial crack defects have a greater influence on the stability of bushing. The decrease of explosion equivalent can change the failure mode of bushing from bidirectional shear failure to tensile failure, and the influence on overall stability is also reduced. When the explosion source is located on the outer wall of the conductive rod, the shock wave radiated by the cooling oil not only makes the bushing instability time advance, but also causes the outer wall expansion damage range to become wider. The stress concentration and the decrease of effective wall thickness make the shear failure instability of bushing with initial crack defects develop rapidly.
- transformer bushing /
- smooth particle hydrodynamics method /
- instability time /
- failure mode

HTML全文

随着材料科学研究的深入，研究者们对材料的研究不再局限于常压范围，开始更多地关注高压下材料的结构和物理性质。在高压作用下，晶体的体积被压缩，原子间的间距变小，电子轨道发生变化，甚至原子的周期性排列也会改变，使晶体结构发生相变^[1-3]。晶体结构在高压环境下具有较为丰富的相变序列，许多常压下无法合成的材料可在高压环境下合成，因此利用高压技术，可以丰富人们对物质世界的认识^[4-5]。

高压下聚合氮以及氮化物的物理、化学性质受到较多关注。将氮作为媒介来合成不同结构的含氮材料是研究聚合氮及氮化物的基础^[6-8]，如利用氮单质在高温高压环境下合成立方氮^[9]，利用氮元素和其他元素合成具有超硬性质的化合物等^[10]。在氮的化合物中，金属氮化物的潜在应用价值巨大，如超导材料^[11]、特殊化学制品^[12]、热稳材料、高能量密度材料和超硬材料等^[13-14]。针对金属氮化物的研究也愈发广泛。

对于碱土金属氮化物的理论研究较多。Be₃N₂的晶体结构研究结果表明，R3m和$P\bar 3m1$两个晶体结构的维氏硬度分别为51 GPa和54 GPa，说明两种结构的Be₃N₂是潜在的超硬材料^[15]。在高压下发现Ca与N的多个配比结构，得出CaN₃的C₂/c结构中N原子形成孤立的六环；进一步对键长和键级进行计算，发现该结构的CaN₃是潜在的高能量密度材料^[16]。Mg₃N₂在高压下具有丰富的物理性质。研究发现，H₂分子很容易在τ-Mg₃N₂、α-Mg₃N₂和α′-Mg₃N₂相的表面分解为H原子并进入材料内部，说明Mg₃N₂是潜在的储氢材料^[17]。

由上述研究可知，碱土金属氮化物的物理性质丰富，潜在应用价值巨大。为此，对高压下Mg的多氮化合物的晶体结构和物理性质进行研究具有十分重要的意义。本研究结合晶体结构预测技术和基于密度泛函理论的第一性原理方法，在0～100 GPa的压强范围内，对MgN₈的晶体结构进行预测，并系统分析预测得到的α-MgN₈、β-MgN₈和γ-MgN₈相，所得结论对进一步研究碱土金属氮化物具有一定的参考价值。

1. 计算细节

采用基于粒子群优化算法的CALYPSO软件，对MgN₈进行结构搜索^[18]，压强范围为0～100 GPa，模拟晶胞采用2倍胞和4倍胞。对于搜索出的结构，使用VASP软件包进行结构优化^[19]，其中描述电子间的交换关联势能采用广义梯度近似（GGA）下的Perdew-Burke-Ernzerh（PBE）交换关联泛函^[19]，赝势采用全电子投影缀加平面波方法^[20]，Mg原子的价电子为2p⁶3s²，N原子的价电子为2s²2p³。为了保证计算精确，使能量收敛精度小于1 meV/atom，经过能量收敛测试，得到平面波的截断能为600 eV；第一布里渊区积分采用Monkhorst-Pack网格方法，网格的间距为0.2 nm^–1；优化过程中，在确保空间群不改变的情况下，选用10^–5 eV为自洽能量收敛的最小值，优化应力收敛设置为0.001 eV/Å。对于优化后的结构，使用超晶胞的方法，利用PHONOPY软件计算声子色散关系和声子态密度，选取不同的扩胞数和k点网格：α-MgN₈相为1 × 1 × 1和7 × 7 × 6；β-MgN₈相为2 × 2 × 4和6 × 6 × 9；γ-MgN₈相为3 × 1 × 1和9 × 5 × 4。使用VASP软件包计算电子局域函数。

2. 结果与讨论

2.1 晶体结构预测

运用CALYPSO软件搜索，在0～100 GPa压强范围内，对预测出的MgN₈晶体结构进行晶格常数和原子位置的优化，计算时设定的温度为0 K，根据G = H – TS（G为吉布斯自由能，H为焓，T为温度，S为熵），可以用焓代替系统的自由能。通过计算优化后结构的焓值随压强的变化关系，可以得到晶体结构的热力学稳定区间。绘制出的焓差随压强变化的曲线如图1所示。由图1（a）可知：在常压下，空间群为P4/mbm的MgN₈晶体结构的焓值最低，记为α-MgN₈相；当压强达到24.3 GPa时，空间群为P4/mnc的晶体结构具有更低的焓值，记为β-MgN₈相；当压强达到68.3 GPa时，空间群为Cmcm的晶体结构能量更低，记为γ-MgN₈相。图1（b）给出了α-MgN₈、β-MgN₈和γ-MgN₈相的晶胞体积随压强的变化。分析可知：α-MgN₈相的晶胞体积随着压强的增大而减小，晶胞体积随压强连续变化；当压强达到24.3 GPa时，MgN₈发生体积坍缩，坍缩率为18.87%；在68.3 GPa时，发生结构相变，体积塌缩率为2.62%；两次相变均属于一级相变。

图 1 MgN₈的焓差曲线以及α-MgN₈、β-MgN₈、γ-MgN₈相体积随压强的变化关系

Figure 1. Calculated enthalpies per formula unit of pressure of MgN₈ and the calculated volume versus pressure of α-MgN₈, β-MgN₈ and γ-MgN₈

下载: 全尺寸图片幻灯片

α-MgN₈、β-MgN₈和γ-MgN₈相的晶体结构见图2，其平衡态晶格常数以及原子位置信息列于表1和表2。

表 1 α-MgN₈、β-MgN₈和γ-MgN₈相的平衡态晶格常数

Table 1. Lattice parameters of α-MgN₈, β-MgN₈ and γ-MgN₈ in the equilibrium state

Phase	Pressure/GPa	a/Å	b/Å	c/Å	α/(°)	β/(°)	γ/(°)
α-MgN₈ (P4/mbm)	0	5.913	5.913	6.572	90.0	90.0	90.0
β-MgN₈ (P4/mnc)	24.3	6.219	6.219	3.724	90.0	90.0	90.0
γ-MgN₈ (Cmcm)	68.3	4.167	4.167	8.680	90.0	90.0	51.1

下载: 导出CSV

| 显示表格

图 2 α-MgN₈、β-MgN₈和γ-MgN₈相的晶体结构

Figure 2. Crystal structures of α-MgN₈, β-MgN₈ and γ-MgN₈

下载: 全尺寸图片幻灯片

表 2 α-MgN₈、β-MgN₈和γ-MgN₈相的平衡态晶格常数和原子位置

Table 2. Lattice parameters and atomic coordinate of α-MgN₈, β-MgN₈ and γ-MgN₈ in the equilibrium state

Phase	Pressure/GPa	Atoms	Wyckoff position	x	y	z
α-MgN₈ (P4/mbm)	0	Mg	2b	0.000	0.000	0.500
		N1	2a	0.000	0.000	0.000
		N2	2c	0.000	0.500	0.500
		N3	4e	0.000	0.000	0.821
		N4	4f	0.000	0.500	0.915
		N5	4h	0.859	0.359	0.500
β-MgN₈ (P4/mnc)	24.3	Mg	2a	0.000	0.000	1.000
β-MgN₈ (P4/mnc)	24.3	N	16i	0.416	0.219	0.823
γ-MgN₈ (Cmcm)	68.3	Mg	8g	0.624	–0.376	1.250
		N1	16h	0.226	–0.135	0.617
		N2	16h	1.115	0.471	1.386

下载: 导出CSV

| 显示表格

图2(a)所示的结构为0 GPa下优化后的α-MgN₈相结构。在α-MgN₈相结构中，每个Mg原子被6个N原子包围，构成八面体结构。其中，Mg-N的键长有2类：一类是中心Mg原子与八面体上下顶点的两个N相连，键长为2.280 Å；另一类为中心的Mg原子与八面体中其余4个N相连，键长为2.109 Å。α-MgN₈相的晶格结构常数：a = b = 5.913 Å，c = 6.572 Å，α = β = γ = 90.0°。其中，Mg原子的Wyckoff占位为2b (0.000, 0.000, 0.500)，N原子的Wyckoff占位分别为2a (0.000, 0.000, 0.000)、2c (0.000, 0.500, 0.500)、4e (0.000, 0.000, 0.821)、4f (0.000, 0.500, 0.915)和4h (0.859, 0.359, 0.500)。图2(b)所示的结构为24.3 GPa下优化后的β-MgN₈相结构。在该结构中，每个Mg原子被8个N原子包围，8支Mg—N键的键长相同，为2.184 Å，构成正六面体。该结构的晶体结构常数：a = b = 6.219 Å，c = 3.724 Å，α = β = γ = 90.0°。其中，Mg原子的Wyckoff占位为2a (0.000, 0.000, 1.000)，N原子Wyckoff占位为16i (0.416, 0.219, 0.823)。图2(c)所示的结构为68.3 GPa下优化后的γ-MgN₈相结构。在该结构中，每个Mg原子与8个N原子相连，构成六面体结构。六面体中的8支Mg-N键有2种键长：4支2.081 Å和4支2.069 Å。该结构的晶体结构常数：a = 4.167 Å，b = 4.167 Å，c = 8.680 Å，α = β = 90.0°，γ = 51.1°。其中，Mg原子的Wyckoff占位为8g (0.624, –0.376, 1.250)，N原子Wyckoff占位分别为16h (0.226, –0.135, 0.617)和16h (1.115, 0.471, 1.386)。

2.2 预测的MgN₈晶体结构的稳定性

通过计算声子色散关系，可以判断晶格结构是否具有动力学稳定性。晶格结构具有动力学稳定性的判断条件是，所有简正声子频率都为有限的实值^[21]，如果是虚值，就可以判断此结构出现了声子软化现象，晶格结构不稳定。计算了预测结构的声子谱和声子态密度，如图3所示。

图 3 α-MgN₈、β-MgN₈、γ-MgN₈相声子色散关系和声子态密度

Figure 3. Phonon-dispersion curves and the phonon density of states of α-MgN₈, β-MgN₈ and γ-MgN₈

下载: 全尺寸图片幻灯片

由图3(a)可知，0 GPa时α-MgN₈相在整个布里渊区未出现声子软化现象，说明该相具有动力学稳定性，最大光学支的频率为71.35 THz。由图3(b)和图3(c)可知，24.3 GPa和68.3 GPa时，β-MgN₈和γ-MgN₈相具有动力学稳定性，最大光学支频率分别为44.64 THz和45.63 THz。

2.3 预测的MgN₈晶体结构的电子性质

为了探究预测得到的MgN₈晶体结构的电子性质，分别研究了α-MgN₈相在0 GPa、β-MgN₈相在24.3 GPa和γ-MgN₈相在68.3 GPa时的能带结构和电子态密度，如图4所示。

图 4 α-MgN₈、β-MgN₈、γ-MgN₈相的能带结构和电子态密度

Figure 4. Band structure and partial density of states of α-MgN₈, β-MgN₈ and γ-MgN₈

下载: 全尺寸图片幻灯片

由能带图可知，压强为0 GPa时，α-MgN₈相的导带与价带之间有3.093 eV的带隙，表明该相结构具有非金属性；当压强为24.3 GPa和68.3 GPa时，β-MgN₈和γ-MgN₈相的导带与价带跨越费米面发生交叠，表明这两相结构具有金属性。

由电子态密度的计算结果可知，α-MgN₈、β-MgN₈和γ-MgN₈相费米能级处的电子态密度均主要由N原子的p轨道贡献，而N原子的s轨道和Mg原子的s、p轨道对电子态密度的贡献相对较少。N原子的p轨道与s轨道之间存在杂化，形成了N—N共价键。

为了研究α-MgN₈、β-MgN₈和γ-MgN₈相晶体结构的化学键，计算了这3个相的电子局域函数（Electron localization function, ELF）^[22]，如图5所示。

图 5 α-MgN₈、β-MgN₈、γ-MgN₈相的电子局域函数

Figure 5. Electron localization function of α-MgN₈, β-MgN₈ and γ-MgN₈

下载: 全尺寸图片幻灯片

图5(a)显示了α-MgN₈结构的三维电子局域函数。α-MgN₈相的N原子之间形成N2和N3，其中：N2的两个N原子之间存在电子局域，形成了共价键，键长为1.116 Å，两个N原子外围存在孤对电子；N3中间的N原子与相邻的两个N原子之间存在电子局域，形成两种N3共价键，其键长分别为1.178 Å和1.182 Å，N3两侧的N原子周围也存在孤对电子。这3种N-N键的键长介于N≡N三键(1.10 Å)和N=N双键(1.25 Å)之间。图5(b)显示了压强为24.3 GPa时β-MgN₈结构的三维电子局域函数。该相N原子之间形成链状结构，每个N原子与相邻的两个N原子之间均有电子局域，形成共价键，键长为1.308 Å和1.316 Å，这两种N-N键的键长介于N=N双键(1.25 Å)和N—N单键(1.45 Å)之间。图5(c)显示68.3 GPa下γ-MgN₈结构的三维电子局域函数。由图5(c)可知，N原子之间形成链状结构，每个N原子与相邻的两个N原子之间均有电子局域，形成共价键，键长为1.282、1.297和1.319 Å，3种N-N键的键长介于N=N双键（1.25 Å）和N—N单键（1.45 Å）之间。比较上述3相的N原子键长数据可以发现，β相和γ相的N-N键更容易断裂，释放出大量的能量，从而成为潜在的高能量密度材料。

图5(d)、图5(e)、图5(f)为MgN₈晶体3种相ELF的二维切面图，晶格切面分别选取(010)、(110)和(010)，ELF的数值范围为0～1。ELF = 1表示电子完全局域化，ELF = 0表示电子完全离域化或该处没有电子。从图5(d)、图5(e)、图5(f)中可以看出：在Mg原子与N原子之间存在电子局域，并且N原子附近对应的ELF值接近1，为高度局域的电子分布，说明Mg原子与N原子之间存在极性共价键；N原子与N原子之间对应的ELF值接近1，表明N原子与N原子之间存在共价键。

为了清晰地描述N原子与Mg原子之间的电子转移情况，计算了Bader电荷转移^[23-24]，计算结果见表3。N原子的电负性强于Mg原子，故N原子吸引电子的能力更强，Mg原子与N原子之间电荷转移的施主为Mg原子，受主为N原子。当压强为0 GPa时，α-MgN₈相每个Mg原子向周围N原子的电荷转移为1.74e。由于该结构中的N原子不仅与Mg原子相连，还形成了N2和N3结构，N2与N3结构的存在导致N原子既可以作为受主，接收来自Mg原子和其他N原子的电荷，又可以作为施主，将电荷转移给其他N原子；当压强为24.3 GPa时，β-MgN₈相中每个Mg原子向周围N原子的电荷转移为1.94e；当压强为68.3 GPa时，γ-MgN₈相中每个Mg原子向周围N原子的电荷转移为2.20e。由此可见，随着压强的增加，电荷转移逐渐增多。

表 3 α-MgN₈、β-MgN₈、γ-MgN₈相电荷转移

Table 3. Calculated Bader charges of α-MgN₈, β-MgN₈ and γ-MgN₈

Phase	Pressure/GPa	Atoms	Number	Charge value/e	Charge transfer/e
α-MgN₈ (P4/mbm)	0	Mg	1	6.26	1.74
		N1	1	4.84	0.16
		N2	1	4.86	0.13
		N3	2	5.04	–0.04
		N4	2	5.46	–0.46
		N5	2	5.52	–0.52
β-MgN₈ (P4/mnc)	24.3	Mg	1	6.06	1.94
β-MgN₈ (P4/mnc)	24.3	N	8	5.24	–0.24
γ-MgN₈ (Cmcm)	68.3	Mg	1	5.80	2.20
		N1	4	5.28	–0.28
		N2	4	5.27	–0.27

下载: 导出CSV

| 显示表格

3. 结　论

运用CALYPASO软件并结合密度泛函理论第一性原理计算方法，对MgN₈的晶体结构进行预测，同时对预测得出的晶体结构进行了结构优化。研究表明：在0～100 GPa范围内，3个相的相变序列为P4/mbm→P4/mnc→Cmcm；两次相变过程中，体积均发生坍塌，坍塌率分别为18.87%和2.62%。电子性质研究表明：α-MgN₈相结构具有非金属性，带隙为3.093 eV；β-MgN₈和γ-MgN₈相结构具有明显的金属性。电子局域函数显示，Mg原子与N原子之间形成极性共价键，N原子与N原子之间形成共价键。高压β-MgN₈和γ-MgN₈相更适合用作储能材料，有望成为高能量密度材料的候选。Bader电荷转移计算结果表明，N原子具有更强的电负性，Mg原子与N原子之间电荷转移的施主为Mg原子，受主为N原子。在α-MgN₈、β-MgN₈、γ-MgN₈相结构中，每个Mg原子向周围N原子的电荷转移分别为1.74e、1.94e和2.20e，即随着压强的增加，电荷转移增多。研究结果可以为碱土金属氮化物的进一步研究提供参考。

图 1 变压器套管二维SPH模型

Figure 1. Two-dimensional SPH model of transformer bushing

下载: 全尺寸图片幻灯片

图 2 模拟验证中的SPH模型

Figure 2. SPH model for validation analysis

下载: 全尺寸图片幻灯片

图 3 柱壳外表面膨胀速度

Figure 3. Expansion velocity of outer surface of cylinder

下载: 全尺寸图片幻灯片

图 4 套管爆炸损伤过程

Figure 4. Damage process of bushing under explosion

下载: 全尺寸图片幻灯片

图 5 套管沿壁厚的塑性应变分布

Figure 5. Plastic strain distribution along the wall thickness of bushing

下载: 全尺寸图片幻灯片

图 6 粒子节点选取示意图

Figure 6. Schematic diagram of particle node selection

下载: 全尺寸图片幻灯片

图 7 套管失稳时间判断曲线

Figure 7. Judgement curves of bushing instability time

下载: 全尺寸图片幻灯片

图 8 套管损伤特征图像及粒子失稳时间判断曲线

Figure 8. Damage characteristics of bushing and judgment curves of particle instability time

下载: 全尺寸图片幻灯片

图 9 对比试件的二维SPH模型

Figure 9. Two-dimensional SPH model of comparison specimen

下载: 全尺寸图片幻灯片

图 10 对比试件的损伤特征

Figure 10. Damage characteristic of comparison specimen

下载: 全尺寸图片幻灯片

图 11 对比试件的套管失稳时间判断曲线

Figure 11. Judgment curves of bushing instability time of comparison specimen

下载: 全尺寸图片幻灯片

图 12 无冷却油套管的损伤特征

Figure 12. Damage characteristic of bushing without cooling oil

下载: 全尺寸图片幻灯片

图 13 无冷却油套管的失稳时间判断曲线

Figure 13. Judgment curves of instability time of bushing without cooling oil

下载: 全尺寸图片幻灯片

图 14 两类裂纹缺陷套管

Figure 14. Two types of bushing with different crack defects

下载: 全尺寸图片幻灯片

图 15 含裂纹缺陷套管的损伤特征

Figure 15. Damage characteristic of bushing with different crack defects

下载: 全尺寸图片幻灯片

图 16 含裂纹缺陷套管沿壁厚方向的塑性应变

Figure 16. Plastic strain of bushing with different crack defects along the wall thickness

下载: 全尺寸图片幻灯片

图 17 含裂纹缺陷套管的失稳时间判断曲线

Figure 17. Judgment curves of instability time of bushing with different crack defects

下载: 全尺寸图片幻灯片

表 1 套管的主要材料属性^[14]

Table 1. Main material properties of bushing^[14]

p_HEL/GPa	Density/(kg·m⁻³)	Shear modulus/GPa	Tensile strength/GPa	${\dot \varepsilon }{_0}$/s^–1
3.63	3 280	157	0.26	1.0

下载: 导出CSV

表 2 引线的主要材料属性^[15–16]

Table 2. Main material properties of wire^[15–16]

Density/(kg·m⁻³)	A/MPa	B/MPa	n	m
8 960	905	292	0.31	1.095

c_p/(J·kg⁻¹·K⁻¹)	C	c/(m·s⁻¹)	S₁	γ₀
383	0.025 3	4 578	1.33	1.67

下载: 导出CSV

表 3 炸药的主要材料属性^[18]

Table 3. Main material properties of explosive^[18]

Density/(kg·m⁻³)	A_e/GPa	B_e/GPa	R₁	R₂	ω	E₀/GPa
1 630	540	9.4	4.5	1.1	0.35	8

下载: 导出CSV

表 4 变压器油的主要材料属性^[15–16]

Table 4. Main material properties of transformer oil^[15–16]

Density/(kg·m^–3)	c/(m·s⁻¹)	S₁	γ₀
895	1 480	1.75	0.28

下载: 导出CSV

表 5 不同TNT当量下响应的数值模拟结果对比

Table 5. Comparison of numerical results of response under different TNT equivalents

W/(g·cm⁻¹)	Failure mode	Expansion velocity/(m·s⁻¹)	Instability time/μs
6.0	Bidirectional shear failure (three locations)	289.18	48.8
4.0	Bidirectional shear failure (one location)	196.09	60.1
3.2	Mixed tensile-shear failure	188.06	66.9
2.1	Tensile failure	164.97	82.6

下载: 导出CSV

参考文献(21)

[1]	谢强, 孙新豪, 赖炜煌. 变压器-套管体系抗震加固理论分析及振动台试验 [J]. 中国电机工程学报, 2020, 40(19): 6390–6399. XIE Q, SUN X H, LAI W H. Theoretical analysis and shaking table test on seismic reinforcement of transformer-bushing system [J]. Proceedings of the CSEE, 2020, 40(19): 6390–6399.
[2]	陈星, 谢强, 李晓璇, 等. 地震作用下变压器侧壁套管的理论建模及摆动效应分析 [J]. 电网技术, 2020, 44(1): 114–121. CHEN X, XIE Q, LI X X, et al. Seismic theoretical modelling and rocking effect analysis on transformer lateral bushing [J]. Power System Technology, 2020, 44(1): 114–121.
[3]	何畅, 谢强, 马国梁, 等. ±800 kV换流变压器-套管体系的抗震性能 [J]. 高电压技术, 2018, 44(6): 1878–1883. HE C, XIE Q, MA G L, et al. Seismic behavior of ±800 kV UHV converter transformer and bushing system [J]. High Voltage Engineering, 2018, 44(6): 1878–1883.
[4]	HENDERSON C N, DEFRANCE C S, PREDECKI P, et al. Ballistic fragmentation confinement of coated brittle transformer bushing models [J]. International Journal of Impact Engineering, 2018, 122: 363–373.
[5]	张少华, 蒋慧灵, 刘国平. 某特高压变压器爆炸事故数值模拟重构 [J]. 消防科学与技术, 2021, 40(3): 356–359. ZHANG S H, JIANG H L, LIU G P. Numerical simulation and reconstruction of explosion accident of an UHV transformer [J]. Fire Science and Technology, 2021, 40(3): 356–359.
[6]	余阳, 阮娅琳, 杨黎波, 等. 换流变压器油蒸汽爆炸的三维仿真研究 [J]. 电瓷避雷器, 2021(2): 155–161. YU Y, RUAN Y L, YANG L B, et al. 3D simulation of oil vapor explosion in converter transformer [J]. Insulators and Surge Arresters, 2021(2): 155–161.
[7]	赵春风, 何凯城, 卢欣, 等. 弧形双钢板混凝土组合板抗爆性能数值研究 [J]. 爆炸与冲击, 2022, 42(2): 025101. doi: 10.11883/bzycj-2021-0205 ZHAO C F, HE K C, LU X, et al. Numerical study of blast resistance of curved steel-concrete-steel composite slabs [J]. Explosion and Shock Waves, 2022, 42(2): 025101. doi: 10.11883/bzycj-2021-0205
[8]	赵铮, 李晓杰, 闫鸿浩, 等. 爆炸压实过程中颗粒碰撞问题的SPH法数值模拟 [J]. 高压物理学报, 2007, 21(4): 373–378. doi: 10.11858/gywlxb.2007.04.007 ZHAO Z, LI X J, YAN H H, et al. Numerical simulation of partides impact in explosive-diven compaction process using SPH method [J]. Chinese Journal of High Pressure Physics, 2007, 21(4): 373–378. doi: 10.11858/gywlxb.2007.04.007
[9]	赵春风, 何凯城, 卢欣, 等. 双钢板混凝土组合板抗爆性能分析 [J]. 爆炸与冲击, 2021, 41(9): 095102. doi: 10.11883/bzycj-2020-0291 ZHAO C F, HE K C, LU X, et al. Analysis on the blast resistance of steel concrete composite slab [J]. Explosion and Shock Waves, 2021, 41(9): 095102. doi: 10.11883/bzycj-2020-0291
[10]	LIU G R, LIU M B. Smoothed particle hydrodynamics: a meshfree particle method [M]. SG Singapore: World Scientific Publishing Company, 2003.
[11]	FRAGASSA C, TOPALOVIC M, PAVLOVIC A, et al. Dealing with the effect of air in fluid structure interaction by coupled SPH-FEM methods [J]. Materials, 2019, 12(7): 1162.
[12]	姚安林, 田晓建, 徐涛龙, 等. 管道爆炸对同沟邻管的冲击效应及防爆墙抗爆性能 [J]. 石油学报, 2020, 41(6): 753–761. doi: 10.7623/syxb202006010 YAO A L, TIAN X J, XU T L, et al. Impact effect of pipeline explosion to adjacent pipelines in the same trench and explosion-proof performance of explosion-proof wall [J]. Acta Petrolei Sinica, 2020, 41(6): 753–761. doi: 10.7623/syxb202006010
[13]	强洪夫, 孙新亚, 王广, 等. 钢箱内部爆炸破坏的SPH数值模拟 [J]. 爆炸与冲击, 2019, 39(5): 052201. doi: 10.11883/bzycj-2017-0439 QIANG H F, SUN X Y, WANG G, et al. Numerical simulation on steel box damage under internal explosion by smoothed particle hydrodynamics [J]. Explosion and Shock Waves, 2019, 39(5): 052201. doi: 10.11883/bzycj-2017-0439
[14]	CRONIN D S, BUI K, KAUFMANN C, et al. Implementation and validation of the Johnson-Holmquist ceramic material model in LS-DYNA [C]//4th European LS-DYNA Users Conference, 2003: 47-60.
[15]	HALLQUIST J O. LS-DYNA keyword user’s manual version 971 [M]. Livermore, USA: Livermore Software Technology Corporation, 2007.
[16]	辛春亮. 有限元分析常用材料参数手册 [M]. 北京: 机械工业出版社, 2019. XIN C L. Manual of material parameters for finite element analysis [M]. Beijing: Mechanical Industry Press, 2019.
[17]	李晓杰, 赵春风, 于娜, 等. TNT炸药和乳化炸药驱动飞板的通用状态方程特征线法研究 [J]. 高压物理学报, 2012, 26(4): 462–468. doi: 10.11858/gywlxb.2012.04.016 LI X J, ZHAO C F, YU N, et al. Characteristic curve method for movement of flyer driven by TNT and emulsion explosives with general equation of state [J]. Chinese Journal of High Pressure Physics, 2012, 26(4): 462–468. doi: 10.11858/gywlxb.2012.04.016
[18]	赵铮, 陶钢, 杜长星. 爆轰产物JWL状态方程应用研究 [J]. 高压物理学报, 2009, 23(4): 277–282. doi: 10.3969/j.issn.1000-5773.2009.04.007 ZHAO Z, TAO G, DU C X. Application research on JWL equation of state of det onation products [J]. Chinese Journal of High Pressure Physics, 2009, 23(4): 277–282. doi: 10.3969/j.issn.1000-5773.2009.04.007
[19]	吴思思, 董新龙, 俞鑫炉. 45钢柱壳爆炸膨胀断裂的SPH模拟分析 [J]. 爆炸与冲击, 2021, 41(10): 103101. doi: 10.11883/bzycj-2021-0172 WU S S, DONG X L, YU X L. An investigating on explosive expanding fracture of 45 steel cylinders by SPH method [J]. Explosion and Shock Waves, 2021, 41(10): 103101. doi: 10.11883/bzycj-2021-0172
[20]	谢富佩, 徐绯, 曾卓, 等. 复合圆柱壳冲击压缩数值模拟及稳定性研究 [J]. 爆炸与冲击, 2021, 41(11): 112201. doi: 10.11883/bzycj-2020-0431 XIE F P, XU F, ZENG Z, et al. Numerical simulation on stability of compositecylindrical shell under impact compression [J]. Explosion and Shock Waves, 2021, 41(11): 112201. doi: 10.11883/bzycj-2020-0431
[21]	MOTT N F. Fragmentation of shell cases [J]. Proceedings of the Royal Society of London, Series A, Mathematical and Physical Sciences, 1947, 189(1018): 300–308.