CST飞行器表面沉积静电分布仿真研究

CST飞行器表面沉积静电分布仿真研究 飞机在高空中飞行时会因与空间粒子进行摩擦等因素在飞机表面形成静电积累。随着飞机表面的电荷的不断地积累飞机表面的电压不断升高。当飞机表面积累电荷所产生电场高于空气击穿场强时就会产生静电放电。静电放电所产生的电磁干扰会影响民用和军用飞机的安全运行。Beach R 对飞机飞行中表面电荷积累原理进行了研究。Joseph E 等研究了喷气式飞机的静电起电原因及放电电流得到了静电放电电流的一般规律。Terhaseborg J L 等研究了子弹在空气中高速飞行时的带电情况在一定范围内预测弹丸上的电荷数量和飞行弹丸周围的电场分布。Grosshans H 等详细阐述并验证了摩擦带电模型以预测直升机在尘土飞扬的空中盘旋时的电气化情况对电荷积累的位置进行详细分析。Lekas T I 提出了一种计算在粉尘云中飞行的飞机静电荷量的方法证明了由灰尘引起的飞机静电荷电与任何其他静电荷源同样重要该方法也可用于评价雨滴、雪花、冰晶等其它大气颗粒物引起的静电荷电。Revel I等提出了在飞机表面布置探测器定位飞行中的飞机沉积静电放电的方法。Andersen A 等研究了聚合物材料的静电放电分布改善了空间环境中的静电放电击穿场的估计。Yadav R 等人详细介绍航空中各种聚合物纳米复合材料的开发飞机抗静电应用的发展以及静电现象及其在航空电子系统中普遍使用条件。易鸣等研究了固定翼飞机静电分布特性及着陆时静电泄放用矩量法分别求取了飞机空中飞行和停靠地面时的电荷分布。杨真一等研究了放电刷对飞行器静电放电的抑制作用得到增加放电刷直径是排放机体电荷最为有效的方法。胡晓峰等人对飞机表面的摩擦起电规律进行了理论分析和实验研究得出了飞机的摩擦极性和电荷定律。张靖等研究了某型飞机静电放电刷的静电泄放特性和影响参数。刘浩等研究了某飞行器表面硅基热防护材料静电起电和泄漏特性及影响因素。左曦等人通过对飞机充氧阀头系统的两种材料进行试验探究了接触压力、环境相对湿度、分离速度等因素对不同材料接触静电积累的影响为进一步探索静电特性上的绝缘材料提供实验依据与研究方向。可以看出国内外对飞机表面的静电起电问题均开展了研究然而上述研究大多未考虑处于飞行状态中的飞机因此本文将针对飞行状态下的飞机研究其表面的沉积静电分布规律。本文研究了CST对飞行状态下塞斯纳飞机表面沉积静电电荷与电场分布。本文首先建立了某型飞机的 1:1 沉积静电仿真模型。然后基于该模型仿真了飞行状态下飞机的电容将仿真所得飞机电容与经验公式所得飞机电容进行对比验证了仿真的正确性。最后研究了模型结构、沉积电荷量对飞机表面电荷密度与电场分布的影响。得到了模型结构对飞机沉积静电分布的影响还得到了飞机不同放电刷处的沉积静电分布的规律。1. 飞机建模本文的研究对象为某型号塞斯纳飞机该飞机机身长度为 14.91 m翼展为 16.72 m高度为 4.57 m。在SolidWorks 中建立该飞机的 1:1 仿真模型如图 1 所示飞机主要由以下几个主要部分构成机身、机翼、引擎、驾驶舱、机舱舷窗和放电刷。其中放电刷结构一致由圆柱和锥台构成圆柱底部半径为 0.005 m长度为 0.15 m锥台底部半径为 0.005 m顶端半径为 0.001 m长度为 0.01 m。机舱舷窗边长为 0.40 m厚度为 0.02 m。飞机仿真模型如图 1 所示。图 1. 1:1 飞机仿真模型.2. 飞行状态下飞机电容计算与仿真由于飞机的几何形状复杂因此无法准确通过解析表达式计算其电容此处采用类比的方法求解飞机电容。根据文献[17]可知自由空间中金属球体的电容可以通过公式1求得图 2. 孤立球的电容仿真结果和相对误差.由图 2 可得当外层球体半径逐渐增大时球体电容逐渐下降最终趋于稳定。当外层球体半径为内层球体半径的 150 倍时仿真所得电容为 11.14 pF与式1所得理论值 11.13 pF相比二者的相对误差仅为0.09%。当外层球体半径大于内层球体半径的 20 倍时仿真所得电容与理论值的相对误差均在 7%以下。因此可得利用CST电磁工作室仿真所得电容与理论值具有很好的一致性因此下面将采用类比的方式计算飞行状态下的飞机电容。图 3. 飞行中的飞机电容计算.如图 3 所示在 CST 电磁工作室中导入 1:1 飞机仿真模型将内层球体替换为飞机进而计算飞行状态下飞机电容。根据图 2 所得结论当外层球体半径为内层球体半径的 20 倍以上时相对误差较小。飞机的半翼长为 8.36 m因此在此设置外层球体半径为 168 m。仿真得到飞行状态下飞机的电容约为 460.5 pF。为进一步验证本文仿真所得飞机电容的正确性根据文献[2]中飞机电容的经验公式可知3. 飞机模型对飞机静电分布的影响分析为了考察不同模型对飞机沉积静电电荷与电场分布规律的影响本文建立两种飞机模型分别为简化模型和细化模型如图 4 所示。简化模型在细化模型的基础上省略驾驶舱、机舱舷窗、引擎结构。为了对比不同飞机模型上同一位置处的电荷与电场在两模型上选取 11 处典型位置如图 4 所示。图 4. 飞机模型图 5. 不同模型的表面电荷密度分布图 6. 不同模型的表面电场强度分布.分别对如图4 所示的两个模型施加相同的电压100 kV仿真该电压作用下飞机表面的电荷密度及电场强度。两模型的机身材料均为PEC驾驶舱与机舱舷窗材料为铅玻璃背景材料为空气。仿真所得飞机表面电荷密度与表面电场强度如图5和图6所示。从图 5 和图 6 可知两种模型表面电荷密度与表面电场强度分布规律相近在飞机尖端部位表面电荷密度与表面电场强度较大如机身两侧机翼尖端尾翼尖端机头顶端在驾驶舱与机舱的舷窗处、机身、机身与机翼连接处的表面电荷密度与表面电场场强较小。为进一步证明不同飞机模型对飞机沉积静电分布的影响将图 4 中所标记位置处的表面电荷密度与表面电场强度的仿真结果列于表 1 和表 2 中其中百分比差异可以通过公式4计算得出其中f1、f2为计算差异二者的取值。由表 1 和表 2 可知细化模型与简化模型相比在机头尖端位置 A 的表面电荷密度与表面电场强度的变化量为 20%左右变化较明显。说明加入驾驶舱后使得电荷集中向机头尖端移动使机头尖端的表面电场强度增高。同时由表 1 和表 2 可知细化模型机头位置处的表面电场强度与表面电荷密最高而简化模型放电刷 2位置处的表面电场强度与表面电荷密度不是最高。这说明驾驶舱对机头尖端处的电荷分布规律有较大影响。表1. 不同位置处的表面电荷密度.表2. 不同位置处的表面电场强度.在机身与机翼连接位置 E 处细化模型与简化模型相比表面电荷密度与表面电场强度的变化量大于 12%。说明加入引擎后部分电荷被转移到了引擎表面这导致位置 E 处细化模型的表面电荷密度与表面电场强度减小。位置 E 处与其他位置处的表面电荷密度与表面电场强度相比其在简化模型与细化模型中均最低。虽然此处的表面电荷密度与表面电场强度变化量较大但对飞机静电放电影响较小。在飞机机身、机翼、尾翼、及放电刷等位置处两种模型的表面电荷密度与表面电场强度差别均在 6%以下这说明驾驶舱、机舱舷窗和引擎结构对飞机机身、机翼、尾翼、及放电刷处的表面电荷密度与表面电场强度影响较小。同时细化模型与简化模型相比其飞机机身、机翼、尾翼及放电刷等位置处的表面电荷密度与表面电场强度的分布规律保持不变。根据上述分析可得细化模型加入驾驶舱、机舱舷窗、引擎结构后对机头处的静电场有较大影响而对其他位置处的影响较小。因此若仅研究放电刷处的静电场分布时可采用简化模型以提高仿真效率。4. 沉积电荷量对飞机放电刷处静电分布的影响分析飞机飞行过程中其表面与空气中的粒子发生碰撞从而不断积累电荷因此其产生的电场强度不断地增加。当其电场强度高于周围空气的击穿场强时将会发生放电现象。为避免发生飞机静电放电放电刷将提前泄放飞机表面的沉积静电。因此放电刷处的沉积静电分布规律对飞机沉积静电放电的研究至关重要。图 7. 放电刷监测位置本节基于飞机细化模型研究了不同电压下放电刷处的飞机沉积静电分布。飞机模型左右两侧对称因此本节仅分析其中一侧放电刷处的电场强度。如图 7 所示飞机放电刷位于飞机机翼尖端部位 a 处、飞机水平尾翼b 处、飞机竖直尾翼 c 处、机尾尖端部位 d 处。其中放电刷 1 和放电刷 9 与机翼呈 45°角度放电刷 2~8、10~11垂直于机翼。在标准大气环境下均匀电场中空气击穿强度为 29 kV/cm[13]。飞机通常飞行在 10 000 m 以上的高空中例如当飞机飞行在 16 000 m 时此时空气临界击穿场强为 419 kV/m[18]。为了比较不同位置放电刷处电场强度的大小当飞机表面沉积静电产生电压为150 kV 时不同位置处放电刷处平均电场强度的仿真结果如图 8 所示。图 8. 150kV 下不同位置处放电刷的电场强度.图 9. 发生空气击穿时飞机沉积静电电压.由图 8 可知飞机机翼尖端部位 a 处放电刷 2 的电场强度最高为 482 kV/m飞机水平尾翼 b 处放电刷 7次之为 471 kV/m二者均超过空气临界击穿场强。飞机水平尾翼 b 处放电刷 8 的电场强度最低为 245 kV/m。放电刷 2 的电场强度比放电刷 8 高约 97%。按照电场强度由高到低依次为放电刷 2、7、3、10、11、4、5、6、12、9、1、8。这说明当发生飞机沉积静电放电时放电刷 2 将先于其他放电刷开始放电。通过仿真得到各放电刷处电场强度达到空气临界击穿场强时所需的飞机沉积静电产生电压如图 9 所示。由图 9 可知不同放电刷处电场强度到达空气临界场强时所对应的飞机沉积静电电压不同。当飞机沉积静电产生电压为 131 kV 时放电刷 2 的电场强度便达到空气临界击穿场强。当飞机沉积静电产生电压为 140 kV 时放电刷 7 的电场强度达到空气临界击穿场强。当飞机沉积静电产生电压为 256 kV 时放电刷 8 的电场强度才达到空气临界击穿场强。当发生空气击穿时放电刷 2 所需的静电产生电压比放电刷 8 低约 95%。这说明不同的飞机沉积静电电压下发生放电的放电刷数量不同。当飞机沉积静电电压为 131 kV 时放电刷 2 已开始放电随着飞机沉积静电电压的增大开始放电的放电刷数量逐步增加。为了分析不同飞机沉积静电电压下飞机机翼尖端部位 a 处、飞机尾翼 b、c 处不同放电刷的电场强度规律将放电刷 1~6 的电场强度绘于图 10将放电刷 7~11 的电场强度绘于图 11。图 10. 机翼尖端部位 a 处放电刷的电场强度.图 11. 尾翼 b、c 处放电刷的电场强度.由图 10 可知不同放电刷处的电场强度随着飞机沉积静电电压的增加而增加且呈线性变化。放电刷 2 的电场强度最高放电刷 1 的电场强度最低。放电刷 2 的电场强度随飞机沉积静电电压的变化率高于其他放电刷按照变化率由高到底依次为放电刷 2、3、4、5、6、1。放电刷 2 的电场强度随着静电电压的变化率比放电刷 1高约 68%。由图 10 可知放电刷 2 的电场强度高于放电刷 3~6放电刷 2 的电场强度随着静电电压的变化率比放电刷 6高约 51%。由图 11 可知放电刷 7 的电场强度高于放电刷 8。放电刷 7 的电场强度随着静电电压的变化率比放电刷 8 高约 82%。同时放电刷 10 的电场强度略高于放电刷 11。这说明垂直于机翼的放电刷中越靠近机翼尖端的放电刷电场强度越大。此外由图 10 可知放电刷 1 的电场强度低于放电刷 2~6放电刷 1 的电场强度随着静电电压的变化率比放电刷 6 低约 11%。同样由图 11 可知放电刷 9 的电场强度低于放电刷 10 与 11。放电刷 9 的电场强度随着静电电压的变化率比放电刷 11 低约 25%。因此可得与机翼呈 45°角度放电刷的电场强度低于同机翼的垂直放电刷电场强度。5. 结论本文建立了 1:1 的飞机飞行状态下的静电场仿真模型基于此模型仿真了飞机电容、飞机沉积静电电荷密度与电场分布得到以下结论通过类比球体电容的方法对飞机电容进行了计算得到飞行状态下飞机的电容约为 460.5 pF。模型对飞机机头处的静电场有影响增加驾驶舱、机舱舷窗后机头处的静电场变化了约 20%而对其他位置处的影响较小。若重点关注飞机放电刷处的静电场分布规律可使用简化模型以提高仿真效率。在不同沉积静电电荷量下飞机的各放电刷处的平均电场强度呈线性变化规律飞机各放电刷处静电场强度不同其中放电刷静电场强度最高处比最低处高约 97%。对于垂直于机翼的放电刷中越靠近机翼尖端的放电刷电场强度越大。垂直于机翼的放电刷处的电场强度高于同机翼的与机翼呈 45°角度的放电刷处的电场强度。资料来源达索官方【相关内容】超表面与天线在 CST 仿真中的根本界限如何使用CST2026版本的无参优化器进行拓扑优化EMC测试会为什么失败CST电磁仿真注意这几点