飞机蒙皮点蚀数值模拟.docVIP

飞机蒙皮点蚀的数值模拟 　　摘 要：文章利用数值模拟方法对飞机蒙皮的点蚀现象进行模拟，根据瞬态图的变化规律，得出以下模拟结果：在瞬态图中模拟的飞机蒙皮点蚀形貌近似为半圆形，但也有其他的不规则图形出现；在模拟条件相同的情况下，随着时间的增加，蚀坑深度逐渐增加，蚀坑面积逐渐增大，腐蚀越来越严重；通过腐蚀瞬态图可以观察到飞机蒙皮表面蚀坑的面积和蚀坑的深度变化规律有相似之处，随时间的增长二者均有增加的趋势。腐蚀概率越大腐蚀越严重，但不同的是腐蚀浓度的影响自始至终都比较明显 关键词：元胞自动机；蚀坑；点蚀；数值模拟 飞机在大气环境中的飞行过程中，由于气流环境以及所受的阻力等因素，飞机表面的蒙皮会发生点蚀，严重影响飞机的安全运行。点蚀是一个非常复杂的过程，不仅应考虑物化和环境因素，还应考虑自然界的各种随机参数[1-3]。除了已有的各种确定性和经验方法外，有必要发展可以预测初期多点腐蚀损害和成长期形态及力学参数的计算机模拟模型。关于腐蚀前期工作大部分专注于化学过程、电流和电位，以及有限生长模型。虽然腐蚀形态学已经引起了国内外极大的关注，但是对二维生长模型的仿真进行的工作还很少，还没有基于局部规则的多点腐蚀生长演化模型的相关文献。采用元胞自动机（Cellular Automata）（简称CA）方法[4]构建一个离散动力系统模型模拟金属材料多点腐蚀的生长过程，可以在腐蚀过程的CA模型中采用确定性规则描述腐蚀发展的大趋势，引进概率系数对其随机性进行描述，以此来得到更符合实际情况的腐蚀仿真模型。因此，文章利用数值模拟方法对飞机蒙皮的点蚀现象进行模拟，得到飞机蒙皮点蚀失效过程的相关规律，对飞机在实际环境中的长期服役具有重要意义 1 物理模型 在文章模拟中，假设飞机蒙皮表面仅发生了点蚀，而金属在腐蚀过程中会发生一系列局部的电化学溶解或是钝化反应。若不区分阴极、阳极区，则在水溶液中，金属的溶解反应及钝化反应分别为[5]： 溶解反应：Me+H20-MeOHaq+（1/2）H2 钝化反应：Me+H20-MeOHsolid+（1/2）H2 假定上述反应均不会改变局部区域内溶液环境的酸碱性 2 数值模拟 数值模拟采用元胞自动机方法，进行二维元胞自动机设置规则。元胞所分布在空间上的网格点的集合就是元胞空间。理论上，元胞空间可以任意维数的欧几里德空间的规则划分。二维元胞自动机的元胞空间通常可按三角形、四方形或六边形三种网格排列，而且元胞空间可以在各维上无限延展，以利于理论研究和推理[6，7]。但是实际应用中，我们无法在计算机上实现这一理想条件。因为不可能处理无限的网格，元胞空间必须是有限的、有边界的。在该二维元胞自动机中，涉及的元胞有金属元胞、腐蚀元胞、非腐蚀元胞： 文章采用最简单的冯VonNeuman）型邻居类型[8]，如图1所示，一个元胞的上、下、左、右相邻的四个元胞为该元胞的邻居。黑色元胞表示中心元胞，灰色为其邻居元胞，以它们的状态一起来计算中心元胞在下一个时刻的状态。在二维元胞自动机的基础上对点蚀的纵向扩展进行的数值仿真模拟，故左右边界采用周期型边界，上下边界采用对称型边界，为了配合局部规则的进行，每个边界都包含两层元胞 2.1 腐蚀元胞浓度的影响 在文章中，通过查阅相关资料，腐蚀元胞浓度的大小模拟腐蚀环境的恶劣程度，腐蚀元胞浓度大则代表腐蚀环境恶劣，腐蚀元胞浓度小则代表腐蚀环境较好。对腐蚀模拟结果有着重要影响。通过腐蚀瞬态图可以观察到飞机蒙皮表面蚀坑的面积和蚀坑的深度变化规律有相似之处，随时间的增长二者均有增加的趋势。当腐蚀元胞浓度为0.1时，虽然蚀坑面积随时间增加，但其大小基本没有什么变化；当腐蚀元胞浓度为0.3时，随时间增加，蚀坑面积大小变化较大，说明腐蚀元胞浓度越高，蚀坑面积变化越大 2.2 腐蚀概率的影响 除了腐蚀元胞浓度之外，腐蚀概率也是影响腐蚀程度的重要因素之一。腐蚀概率的大小表征腐蚀因子对金属腐蚀能力的强弱，通过瞬态图可以观察到，腐蚀概率越大腐蚀越严重，但不同的是腐蚀浓度的影响自始至终都比较明显。当腐蚀概率为0.1时，蚀坑面积变化不不明显；当腐蚀概率为0.3时，蚀坑面积变化较大，说明腐蚀概率较大时对于蚀坑面积的影响较大 3 结束语 （1）在瞬态图中模拟的飞机蒙皮点蚀形貌近似为半圆形，但也有其他的不规则图形出现 （2）在模拟条件相同的情况下，通过腐蚀瞬态图可以观察到飞机蒙皮表面蚀坑的面积和蚀坑的深度变化规律有相似之处，随着时间的增加，蚀坑深度逐渐增加，蚀坑面积逐渐增大，腐蚀越来越严重，且腐蚀元胞浓度越高，蚀坑面积变化越大 （3）在模拟条件相同的情况下，通过瞬态图可以观察到，腐蚀概率越大腐蚀越严重，但不同的是腐蚀浓度的影响自始至终都比较明显，且腐蚀概率较大