单脉冲电火花放电气泡行为仿真建模

doi:10.3969/j.issn.1004-132X.2026.05.007

摘要/Abstract

摘要：

基于水下爆炸理论和理想气体状态方程，考虑了流体在加工间隙内流动时因壁面阻力导致的能量损失，结合实验观测的气泡最大膨胀体积和稳态体积，确定了极间间隙内单脉冲电火花放电产生气泡的初始温度、压力和体积的计算方法。进一步采用流体体积（VOF）法构建了单脉冲放电产生气泡行为的仿真模型，单脉冲放电气泡观测实验结果表明仿真与实验结果吻合良好。利用该模型分析了峰值电流和脉冲宽度对气泡演变行为的影响。最后基于高速摄像机观测实验数据拟合得出了不同峰值电流和脉冲宽度下气泡最大膨胀体积、膨胀时间及稳态体积的计算公式，且无需观测实验即可对单脉冲放电气泡行为进行仿真计算。

关键词: 电火花加工, 单脉冲, 气泡, 流场

Abstract:

Based on the underwater explosion theory and the ideal gas equation of state， considering the energy loss caused by the wall resistance when the fluid flows in the machining gaps， and combining the maximum expansion volume and the steady-state volume of the bubbles observed in the experiments， a calculation method was determined for the initial temperature， pressure and volume of the bubbles generated by a single-pulse EDM within the inter-electrode gaps. A simulation model for single-pulse discharge bubble behaviors was further developed by using the volume of fluid（VOF） method. The results of the single-pulse discharge bubble observation experiments show that the simulation results are in good agreement with the experimental ones.Using this model， the influences of peak current and pulse duration on bubble evolution behaviors were analyzed. Finally， based on the experimental data obtained from high-speed camera observations， the calculation formulas for the maximum expansion volume， expansion time and steady-state volume of the bubbles were derived under different peak currents and pulse durations. Moreover， the simulation calculation of the single-pulse discharge bubble behaviors may be carried out without observing the experiments.

Key words: electrical discharge machining（EDM）, single-pulse, bubble, flow field

中图分类号:

TG661

王津, 乔春凯, 贾志新, 何虎. 单脉冲电火花放电气泡行为仿真建模[J]. 中国机械工程, 2026, 37(5): 1072-1081.

WANG Jin, QIAO Chunkai, JIA Zhixin, HE Hu. Simulation Modeling of Single-pulse Electrical Discharge Bubble Behaviors[J]. China Mechanical Engineering, 2026, 37(5): 1072-1081.

导出引用管理器 EndNote|Ris|BibTeX

链接本文: https://www.cmemo.org.cn/CN/10.3969/j.issn.1004-132X.2026.05.007

https://www.cmemo.org.cn/CN/Y2026/V37/I5/1072

图/表 22

图1 实验装置示意图

Fig.1 Schematic diagram of the experimental setup

表1 ITO导电玻璃参数

Tab.1 Parameters of ITO conductive glass

型号	IT0导电玻璃
尺寸/（mm×mm×mm）	20×20×1.1
方阻/Ω	10
膜层厚度/nm	185
透光率/%	≥84

图2 气泡稳态观察实验图（脉冲宽度60 μs，峰值电流10 A，极间间隙40 μm）

Fig.2 Experimental diagram of bubble steady-state observation（pulse duration of 60 μs， peak current of 10 A， inter-electrode gap of 40 μm）

图3 平行板间电火花加工极间间隙几何模型

Fig.3 Geometric model of the inter-electrode gap in parallel plate electrical discharge machining

表2 求解器算法设置

Tab.2 Solver algorithm settings

算法	设定值
求解器	压力基、隐式、瞬态
多相流	相数：2；VOF策略：显式、激活体积力；层流
控制方程	开启能量方程
求解控制参数	PISO； Discretization； PRESTO！； QUICK； Geometric
收敛判据	残差：10 $- 6$
迭代计算参数	库朗数：2；各时间步最大迭代次数：40

表2 求解器算法设置

Tab.2 Solver algorithm settings

算法	设定值
求解器	压力基、隐式、瞬态
多相流	相数：2；VOF策略：显式、激活体积力；层流
控制方程	开启能量方程
求解控制参数	PISO； Discretization； PRESTO！； QUICK； Geometric
收敛判据	残差：10 $- 6$
迭代计算参数	库朗数：2；各时间步最大迭代次数：40

图4 平行板间电火花气泡仿真与实验结果对照图

Fig.4 Comparison chart of the simulation and experimental results of the electric spark bubble between parallel plates

图5 气泡体积随时间变化曲线

Fig.5 Curve of bubble volume varying with time

图6 峰值电流为5 A时气泡体积随时间变化曲线

Fig.6 Curve of bubble volume varying with time when the peak current is as 5 A

图7 峰值电流为10 A时气泡体积随时间变化曲线

Fig.7 Curve of bubble volume varying with time when the peak current is as 10 A

图8 峰值电流为15 A时气泡体积随时间变化曲线

Fig.8 Curve of bubble volume varying with time when the peak current is as 15 A

图9 气泡膨胀时间随峰值电流变化曲线图（脉冲宽度90 μs）

Fig.9 Curve graph of bubble expansion time varying with peak current（pulse duration of 90 μs）

图10 脉冲宽度为60 μs时气泡体积随时间变化曲线

Fig.10 Curve of bubble volume varying with time when the pulse duration is as 60 μs

图11 脉冲宽度为120 μs时气泡体积随时间变化曲线

Fig.11 Curve of bubble volume varying with time when the pulse duration is as 120 μs

图12 脉冲宽度为180 μs时气泡体积随时间变化曲线

Fig.12 Curve of bubble volume varying with time when the pulse duration is as 180 μs

图13 气泡膨胀时间随脉冲宽度变化曲线图（峰值电流10 A）

Fig.13 Curve graph of bubble expansion time varying with pulse duration（peak current of 10 A）

图14 不同峰值电流与极间间隙内电火花气泡体积关系曲线（脉冲宽度90 μs）

Fig.14 Relationship curve among different peak currents and the volume of electric spark bubbles in the inter-electrode gap（pulse duration of 90 μs）

图15 不同脉冲宽度与极间间隙内电火花气泡体积关系曲线（峰值电流10 A）

Fig.15 Relationship curve among different pulse durations and the volume of electric spark bubbles in the inter-electrode gap（peak current of 10 A）

图16 气泡体积与电参数拟合结果

Fig.16 Fitting results of bubble volume and electrical parameters

图17 气泡膨胀时间与电参数拟合结果

Fig.17 Fitting results of bubble expansion time and electrical parameters

表3 气泡膨胀时间验证实验结果与误差

Tab.3 Experimental results and errors of bubble expansion time verification

实验组号	实验条件（峰值电流，脉冲宽度）	实验测量结果/μs	公式计算结果/μs	误差/%
1	I=1 A， $t o n$ =5 μs	208	199.936	$-$ 4.84
2	I=2 A， $t o n$ =15 μs	230	219.578	$-$ 4.56
3	I=5 A， $t o n$ =60 μs	369	376	2.17
4	I=10 A， $t o n$ =110 μs	378	396	4.76
5	I=12 A， $t o n$ =170 μs	409	401	$-$ 2.20
6	I=15 A， $t o n$ =180 μs	414	403	$-$ 2.66

表3 气泡膨胀时间验证实验结果与误差

Tab.3 Experimental results and errors of bubble expansion time verification

实验组号	实验条件（峰值电流，脉冲宽度）	实验测量结果/μs	公式计算结果/μs	误差/%
1	I=1 A， $t o n$ =5 μs	208	199.936	$-$ 4.84
2	I=2 A， $t o n$ =15 μs	230	219.578	$-$ 4.56
3	I=5 A， $t o n$ =60 μs	369	376	2.17
4	I=10 A， $t o n$ =110 μs	378	396	4.76
5	I=12 A， $t o n$ =170 μs	409	401	$-$ 2.20
6	I=15 A， $t o n$ =180 μs	414	403	$-$ 2.66

表4 气泡最大膨胀体积验证实验结果与误差

Table.4 Experimental results and errors of the verification for the maximum expansion volume of bubbles

实验组号	实验条件（峰值电流，脉冲宽度）	实验测量结果/mm³	公式计算结果/mm³	误差/%
1	I=1 A， $t o n$ =5 μs	0.051 56	0.046 68	$-$ 9.46
2	I=2 A， $t o n$ =15 μs	0.096 68	0.091 95	$-$ 4.89
3	I=5 A， $t o n$ =60 μs	0.668 92	0.637 11	$-$ 4.76
4	I=10 A， $t o n$ =110 μs	1.863 86	2.016 76	8.20
5	I=12 A， $t o n$ =170 μs	2.623 85	2.754 27	4.97
6	I=15 A， $t o n$ =180 μs	4.216 26	3.915 90	$-$ 7.12

表4 气泡最大膨胀体积验证实验结果与误差

Table.4 Experimental results and errors of the verification for the maximum expansion volume of bubbles

实验组号	实验条件（峰值电流，脉冲宽度）	实验测量结果/mm³	公式计算结果/mm³	误差/%
1	I=1 A， $t o n$ =5 μs	0.051 56	0.046 68	$-$ 9.46
2	I=2 A， $t o n$ =15 μs	0.096 68	0.091 95	$-$ 4.89
3	I=5 A， $t o n$ =60 μs	0.668 92	0.637 11	$-$ 4.76
4	I=10 A， $t o n$ =110 μs	1.863 86	2.016 76	8.20
5	I=12 A， $t o n$ =170 μs	2.623 85	2.754 27	4.97
6	I=15 A， $t o n$ =180 μs	4.216 26	3.915 90	$-$ 7.12

表5 气泡稳态体积验证实验结果与误差

Table.5 Experimental results and errors of bubble steady-state volume verification

实验组号	实验条件（峰值电流，脉冲宽度）	实验测量结果/mm³	公式计算结果/mm³	误差/%
1	I=1 A， $t o n$ =5 μs	0.013 95	0.012 71	$-$ 8.89
2	I=2 A， $t o n$ =15 μs	0.039 56	0.037 06	$-$ 6.31
3	I=5 A， $t o n$ =60 μs	0.566 48	0.542 89	$-$ 4.16
4	I=10 A， $t o n$ =110 μs	1.483 70	1.575 41	6.18
5	I=12 A， $t o n$ =170 μs	1.836 26	1.938 09	5.55
6	I=15 A， $t o n$ =180 μs	2.369 64	2.297 42	$-$ 3.05

表5 气泡稳态体积验证实验结果与误差

Table.5 Experimental results and errors of bubble steady-state volume verification

实验组号	实验条件（峰值电流，脉冲宽度）	实验测量结果/mm³	公式计算结果/mm³	误差/%
1	I=1 A， $t o n$ =5 μs	0.013 95	0.012 71	$-$ 8.89
2	I=2 A， $t o n$ =15 μs	0.039 56	0.037 06	$-$ 6.31
3	I=5 A， $t o n$ =60 μs	0.566 48	0.542 89	$-$ 4.16
4	I=10 A， $t o n$ =110 μs	1.483 70	1.575 41	6.18
5	I=12 A， $t o n$ =170 μs	1.836 26	1.938 09	5.55
6	I=15 A， $t o n$ =180 μs	2.369 64	2.297 42	$-$ 3.05

参考文献 20

[1]	郑倩，张亚欧，高强，等.基于硼酸溶液的电火花加工机理分析及实验研究［J］.中国机械工程，2023， 34（24）： 2936-2945.
	ZHENG Qian， ZHANG Yaou， GAO Qiang， et al. Analysis and Experimental Research on the Electrospark Machining Mechanism Based on Boric Acid Solution［J］. China Mechanical Engineering， 2023， 34（24）： 2936-2945.
[2]	WANG J， HAN F， CHENG G， et al. Debris and Bubble Movements during Electrical Discharge Machining［J］. International Journal of Machine Tools and Manufacture， 2012， 58： 11-18.
[3]	HAYAKAWA S， SASAKI Y， ITOIGAWA F， et al. Relationship between Occurrence of Material Removal and Bubble Expansion in Electrical Discharge Machining［J］. Procedia CIRP， 2013， 6： 174-179.
[4]	KITAMURA T， KUNIEDA M， ABE K. High-speed Imaging of EDM Gap Phenomena Using Transparent Electrodes［J］. Procedia CIRP， 2013， 6： 314-319.
[5]	KITAMURA T， KUNIEDA M， ABE K. Observation of Relationship between Bubbles and Discharge Locations in EDM Using Transparent Electrodes［J］. Precision Engineering， 2015， 40： 26-32.
[6]	LI G， NATSU W， YU Z. Elucidation of the Mechanism of the Deteriorating Interelectrode Environment in Micro EDM Drilling［J］. International Journal of Machine Tools and Manufacture， 2021， 167： 103747.
[7]	LI G， NATSU W， YANG J， et al. Bubble Flushing Effect in Micro EDM Drilling and Its Relation with Debris［J］. Journal of Materials Processing Technology， 2022， 305： 117590.
[8]	孟慧玲. 电火花小孔加工中极间气泡运动及分布研究［D］. 哈尔滨：哈尔滨工业大学， 2022.
	MENG Huiling. Study on Bubble Behavior and Distribution in Deep Hole Drilling with Electrical Discharge Machining［D］. Harbin： Harbin Institute of Technology， 2022.
[9]	WANG C， SASAKI T， HIRAO A. Bubble Behavior of Single-pulse Discharge in EDM［J］. Journal of Manufacturing Processes， 2024， 126： 60-73.
[10]	WANG J， HAN F. Simulation Model of Debris and Bubble Movement in Consecutive-pulse Discharge of Electrical Discharge Machining［J］. International Journal of Machine Tools and Manufacture， 2014， 77： 56-65.
[11]	黄青宇. 电火花加工蚀除物和气泡在间隙中分布状态仿真研究［D］. 哈尔滨：哈尔滨工业大学， 2017.
	HUANG Qingyu. The Simulation Study of the Distribution of Debris and Bubbles in the Gap during EDM［D］. Harbin： Harbin Institute of Technology， 2017.
[12]	GONG S W， OHL S W， KLASEBOER E， et al. Scaling Law for Bubbles Induced by Different External Sources： Theoretical and Experimental Study［J］. Physical Review E：Statistical， Nonlinear， and Soft Matter Physics， 2010， 81（5）： 056317.
[13]	LI T， WANG S， LI S， et al. Numerical Investigation of an Underwater Explosion Bubble Based on FVM and VOF［J］. Applied Ocean Research， 2018， 74： 49-58.
[14]	GONG S W， GOH B H T， OHL S W， et al. Interaction of a Spark-generated Bubble with a Rubber Beam： Numerical and Experimental Study［J］. Physical Review E：Statistical， Nonlinear， and Soft Matter Physics， 2012， 86（2）： 026307.
[15]	GOH B H T， GONG S W， OHL S W， et al. Spark-generated Bubble near an Elastic Sphere［J］. International Journal of Multiphase Flow， 2017， 90： 156-166.
[16]	ZHAO B， CHEN J， WU C， et al. Numerical Simulation of Bubble and Arc Dynamics during Underwater Wet Flux-cored Arc Welding［J］. Journal of Manufacturing Processes， 2020， 59： 167-185.
[17]	LEE M， KLASEBOER E， KHOO B C. On the Boundary Integral Method for the Rebounding Bubble［J］. Journal of Fluid Mechanics， 2007， 570： 407-429.
[18]	RAJESH J V， GIRIDHARAN A. Exploring the Influence of Powder Particles in Dielectric Liquid of EDM： Pulse Train， Discharge Energy， Material Removal Rate and Surface Topography［J］. Physica Scripta， 2023， 98（10）： 105001.
[19]	周佰正. 基于达西-魏斯巴赫方程的直升机主减速器滑油系统故障分析［J］. 航空维修与工程， 2018（7）： 91-92.
	ZHOU Baizheng. Failure Analysis of Helicopter Main Reducer Lubricating Oil System Based on Darcy-Weisbach［J］. Aviation Maintenance & Engineering， 2018（7）： 91-92.
[20]	KUNIEDA M， TAKEUCHI H， KUSAKABE Y， et al. Clarifying EDM Gap Phenomena by Gas Chromatography Analysis of Bubbles［J］. International Journal of Electrical Machining， 2011， 16： 15-19.