Al2O3陶瓷引弧微爆炸加工温度场模拟

张保国; 林克凌; 田欣利; 薛春芳; 李富强

doi:10.11883/1001-1455(2012)02-0129-07

Al2O3陶瓷引弧微爆炸加工温度场模拟

doi: 10.11883/1001-1455(2012)02-0129-07

1.
装甲兵工程学院再制造技术重点实验室,北京100072;
2.
装甲兵工程学院机械工程系,北京100072

基金项目:

国家自然科学基金项目(51075399)

详细信息

通讯作者:
张保国

计量
- 文章访问数: 2391
- HTML全文浏览量: 72
- PDF下载量: 191
- 被引次数: 3
出版历程
- 刊出日期: 2012-03-25

NumericalsimulationoftemperaturefieldforAl2O3ceramics duringmicro-detonationofstrikingarcmachining

1.
ScienceandTechnologyonRemanufacturingLaboratory,
2.
AcademyofArmoredForcesEngineering,Beijing 100072,China;

Funds:

NationalNaturalScienceFoundationofChina(51075399)

More Information

Corresponding author: ZHANG Bao-guo

摘要

摘要: 建立了Al2O3陶瓷引弧微爆炸加工(micro-detonationofstrikingarcmachining,MDSAM)过程的传热模型,基于有限元理论,利用ANSYS软件对加工过程中的温度场分布进行了模拟。结合材料性质,对模拟和实验得到的蚀坑尺寸进行了比较,并分析了加工参数对温度场的影响。模拟结果表明,Al2O3陶瓷引弧微爆炸加工时在给定的加工参数下的最高温度可达13435℃,且高温影响区范围很小,加工实验与模拟结果符合较好。随着脉冲宽度和工作电流的增加,加工区域的温度以及蚀坑的半径和深度增大;随着喷嘴半径的增大,加工区域的温度降低而蚀坑的径深比增大。模拟结果可为Al2O3 陶瓷引弧微爆炸加工过程中表面形貌的预测、材料去除机理的揭示以及加工参数的选择等提供参考。
- 爆炸力学 /
- 温度场 /
- 引弧微爆炸加工 /
- Al2O3 陶瓷
Abstract: AthermaltransfermodelforAl2O3ceramicsduringmicro-detonationofstrikingarcmachining( MDSAM)wasestablished.Basedonthefiniteelementtheory,thetemperaturefieldofAl2O3 ceramicsduringMDSAM wassimulatedwiththeaidoftheANSYSsoftware.Combinedwiththematerialproperties, thesimulatedcavitydimensionwascomparedwiththemachiningexperimentalresultandtheeffectsofthemachiningparametersonthetemperaturef eldwereanalyzed.ThesimulatedresultsshowthatthehighesttemperatureofAl2O3ceramicsduringMDSAMcanreach13435℃ withthegivenmachiningparameters.Theheat-affectedzoneisquitesmallandthemachiningtestis consistentwiththesimulatedresults.Withtheincreaseofthepulsewidthandtheworkingcurrent, thetemperatureoftheprocessingregionandtheradiusanddepthofthecavityincrease.Thetemperaturedecreaseswiththeincreaseofthenozzleradius, whereasthediameter-to-depthratioofthecavity increases.Thesimulatedresultscanprovideanimportantreferenceforthepredictionofthesurface topography,thedisclosureofthematerialremovalmechanism,andtheoptimizationofthemachining parametersforAl2O3ceramicsduringMDSAM.
- mechanicsofexplosion /
- temperaturefield /
- micro-detonationofstrikingarcmachining /
- Al2O3ceramics

HTML全文

航天器外壳受到空间碎片的超高速撞击产生的高速碎片云, 可能对航天器内部的仪器设备造成灾难性的破坏, 甚至导致航天器过早失效。研究碎片云的形状与运动特性, 对于定量预测航天器内部仪器设备的损伤程度, 并进行相应的防护设计, 具有十分重要的意义。因此, 这个问题受到高度重视, 各种不同的碎片云模型被提出或改进。以球形弹丸正撞击靶板为例, H.F.Swift^[1]给出的碎片云运动模型为一个膨胀的球壳, 碎片云质量均匀地分布在球壳上。碎片云的运动可分解为质心沿弹道方向的直线运动和以质心为中心的膨胀。W.P.Schonberg等^[2]给出的碎片云模型与Swift模型类似, 不同的是内部为弹丸碎片云, 外部是一个缓冲板碎片云形成的球壳。然而许多实验结果表明, 碎片云的外形显然不是一个球^[3-4], 而是更象一个旋转椭球面, 因此又为碎片云建立了旋转椭球模型。F.K.Schäfer^[4]提出了一个球形弹丸正撞击金属薄板的碎片云模型, 在忽略反溅碎片云的前提下, 碎片云模型分为3部分:(1)缓冲板材料碎片云, 其质量均匀分布在旋转椭球壳上; (2)弹丸材料碎片云, 分布在一个球壳上, 且头部内切于缓冲板材料形成的旋转椭球壳碎片云; (3)弹丸材料中心大碎片, 位于弹丸材料球壳碎片云头部。然而, 无论是球形碎片云模型, 还是旋转椭球壳碎片云模型, 都没有通过对实验图像进行边缘提取、并与相关模型进行拟合、以证实模型的正确性。

从许多超高速撞击实验获得的碎片云图片可以看出, 碎片云形状与旋转椭球存在明显差异。图 1是4幅超高速撞击实验的碎片云图片^[5-8], 由此可以看出球形弹丸正撞击金属薄板产生碎片云的大致形状:尾部为旋转曲面壳体, 头部为前面呈草帽状的高密度区域^[9]。草帽形状的确定主要是“草帽顶”形状的确定。

图 1 超高速撞击实验碎片云图

Figure 1. Debris clouds from hypervelocity impact experiments

下载: 全尺寸图片幻灯片

本文中, 将通过图像处理方法提取碎片云“草帽顶”的边缘图像, 并以二次曲线的旋转曲面与此拟合, 以选择最优的拟合曲面。

1. 碎片云形状的拟合

由实验可以得到碎片云的三维图像, 但从文献中获取的往往是在平面上的投影。但对于正撞击的情形, 碎片云外形可以视为旋转曲面, 其形状由它在平面上的正投影图像(如不混淆, 以下简称碎片云图像)的形状唯一确定。对于碎片云图像, 通过图像处理方法, 例如图像锐化^[10], 可以提取出边缘曲线。假定(x_i, y_i)是碎片云的边缘曲线, i=1, 2, …, n, n为数据点个数。如果碎片云图像的对称轴平行于坐标轴, 与此拟合的二次函数的一般形式为:

${{y^2} + Ay + B{x^2} + Cx + D = 0}$

(1)

${y + B{x^2} + Cx + D = 0}$

(2)

式中:A、B、C、D为待定拟合系数。式(1)当B < 0时为双曲线, 当B > 0时为椭圆; 式(2)为抛物线。

由于图像本身以及获取过程的各种因素, 图像不可能严格遵循某一二次函数, 所谓其形状为某一二次函数, 只能是更加接近而已。于是就产生了“接近”程度的度量问题。

对于抛物线, 一种常用的度量指标为:

$M=\sum\limits_{i=1}^{n}\left(y_{i}-f\left(x_{i}\right)\right)^{2}$

(3)

式中:f(x)=-Bx²-Cx-D。

而对于双曲线和椭圆的情形, f(x)是由式(1)~(2)确定的y关于x的二值函数之一:f(x) 或

2. 碎片云图像函数拟合方法

对于经图像处理获取的碎片云的边缘曲线(x_i, y_i)的函数拟合, 是为了寻找拟合系数A、B、C、D, 使度量指标M达到最小值。

2.1 抛物线拟合

将f(x)=-Bx²-Cx-D代入式(3)中, 则目标函数为抛物线拟合转化为典型的最小二乘问题。记则度量指标达到最小的B、C、D, 可表示为(B, C, D)^T=(A^TA)^-1A^TY, 由此可得到拟合抛物线f(x)=-Bx²-Cx-D, 及其相应的拟合接近程度度量指标。

2.2 双曲线/椭圆拟合

为避免非线性最小二乘问题的求解, 拟合问题分为两步进行。首先, 利用最小二乘法, 通过与的拟合, 确定拟合系数A、B、C、D的初值, 并以指标M最小为条件选择或

其次, 以A、B、C、D的初值, 利用变半径随机搜索算法^[11]进一步寻求使度量指标M达到最小的拟合系数A、B、C、D。为此, 记则的最小二乘拟合方程为最小二乘解为以A、B、C、D的初值, 利用变半径随机搜索算法就可求出拟合系数使达到最小值, 其中或

3. 碎片云图像辨识

通过边缘提取可以得到碎片云的边缘曲线; 通过二次函数拟合可以得到两种拟合曲线及其与边缘曲线拟合程度的度量指标; 通过对度量指标的选择可以对碎片云图像进行辨识; 通过对图 1中不同来源碎片云图像的辨识, 可以对碎片云形状给出结论。

3.1 图像1

利用上述方法进行抛物线拟合, 拟合系数B、C、D分别为0.007 5、-5.066 7和883.394 7, 度量指标为37.870 5。

进行双曲线/椭圆拟合, 拟合系数初值A、B、C、D分别为260、-1.168、784.638 7和-1.392 7, 度量指标为129.014 2。利用变半径随机搜索算法进行搜索, 但度量指标始终不小于37.870 5, 因此辨识结果为抛物线。最优拟合结果如图 2所示。

图 2 图像1的抛物线拟合结果

Figure 2. Parabola fitting result of image 1

下载: 全尺寸图片幻灯片

3.2 图像2

进行抛物线拟合, 拟合系数B、C、D分别为0.004 3、-2.521 0和536.066 2, 度量指标为664.115 1。

进行双曲线/椭圆拟合, 拟合系数初值A、B、C、D分别为96.018 2、1.364 7、-793.928 6和127 314.910 6, 对应的度量指标为274.413 5, 因此辨识结果为双曲线/椭圆。利用变半径随机搜索算法进行搜索发现, 度量指标变大, 因此进一步的辨识结果是图像为椭圆。拟合系数初值对应的拟合结果如图 3所示。