深圳市远华伟业科技有限公司

远华电气技术(深圳)有限公司

服务热线:0755-32858828
®
CMCE多电场补偿器系列产品请添加微信fangleiren咨询
深圳市远华伟业科技有限公司
我们提供专业的技术服务,以解决您的后顾之忧……

 


    技术资料​  Technology

基于FEKO的直升机雷电间接效应仿真分析
来源: | 作者:深圳远华伟业 | 发布时间 :2016-04-19 | 6166 次浏览: | 分享到:
本文通过对机载设备雷电间接效应进行研究,使用CAE软件——FEKO对其进行数值仿真预测,得到了影响直升机表面电流分布和电磁环境分布的主要原因。

1 引言

    闪电是一种高电压和大电流的自然放电现象,在放电过程中会产生上升时间极快、持续时间极短的大电流脉冲。一般而言,闪电持续最短的时间可能只有几百万分之一秒,但电流却可高达几万安培。这种高脉冲电流会向空间辐射很强的瞬态电磁场,对直升机的飞行构成严重的威胁。同时,现代直升机为了提高综合性能,越来越多的使用先进的复合材料,高灵敏的电子设备和系统,因而在直升机遭到雷击时,造成损失的几率也就越来越大。为了保证直升机的飞行安全,必须对直升机的雷电防护设计和试验进行研究。

    为了解决飞行器的飞行安全,欧美等国率先进行了雷电防护试验,并制定了一系列的雷电试验标准。雷电防护试验虽然可以完成飞行器雷电防护设计的要求,但由于试验成本高昂,试验条件苛刻,人员组织困难,电子设备易损等因素增加了飞行器雷电防护试验的难度和降低了试验执行的效率。因此,采用数值仿真方法对飞行器的雷电效应进行仿真分析是近几年研究的重点。文献对雷电电磁场的场线耦合问题上进行了理论计算、数值分析和试验测试验证。文献对飞行器系统雷电效应的试验仪器、试验方法和测试波形进行回顾,并对飞行器的雷电效应进行了实验和仿真分析。文献使用VAM-LIFE对C-27J运输机的雷电间接效应进行了仿真分析。文献使用TLM方法对空客A320的间接效应进行了仿真分析。由于国内雷电效应的研究起步比较晚,在飞机雷电防护方面缺乏仿真和试验数据,在飞机设计前期,完成对飞机雷电效应的数值仿真分析的对飞机的雷电防护设计具有重要的意义。

    本文通过使用三维电磁仿真软件FEKO对直升机雷电间接效应进行了仿真分析。根据美国工程师协会的SAE-ARP5416和ASE-ARP5412等标准相关试验方法的规定,本文在设定的雷电注入路径下,对直升机表面电流分布、舱外电磁环境及舱内线缆进行了仿真分析。通过分析找出影响直升机表面电流分布、电场分布和磁场分布的重要因素,为直升机的防护设计提供支撑。

2 仿真环境设置

2.1电磁场软件的选择

    全波算法是求解电磁场问题的精确方法,对于给定的计算机硬件资源,此类方法所能仿真的电尺寸有其上限。电磁仿真软件有一个共性,就是它们都与要仿真物体的电尺寸相关,电尺寸定义为被仿真物体的几何尺寸除以所涉及最高频率对应的波长,单位是波长数。雷电的最高频率达到50MHz,即使对于几十米长的现代大型航空客机,也属于中小电尺寸的仿真对象。因此常用的数字仿真全波算法如MoM(矩量法)、FDTD(时域有限积分法)、FEM(有限元法)都是比较有效的工具,可以用来评定大型航空飞机的表面雷电电流分布。

    FEKO是针对天线设计、天线布局与电磁兼容性分析而开发的专业电磁场分析软件,它的核心算法是MoM,原理上MoM可以解决任意复杂结构的电磁问题。航空器表面雷电电流分布的求解问题,属于低频、电磁兼容问题的建模和仿真,采用FEKO的MoM进行分析是非常合适的。

    时域有限差分(FDTD)求解器的加入,标志着FEKO在算法的完整性和解决各类问题的有效性方面进入一个新的里程碑。新增的FDTD求解器易于激活和使用,并可与其它求解器便捷切换。尽管该算法设定在时域,但通过傅立叶变化可以计算宽带频域数据。它采用一阶数值微分可达到二阶的精度。这种算法非常适合非均匀材料建模,也更加适用于多种并行技术。当硬件系统有可用的GPU资源时,FDTD求解器支持GPU计算以提升性能,获得显著加速效果。新的时域求解器的加入,为关注宽带天线和非均匀结构仿真的FEKO用户提供了更加有效的解决方案。

2.2模型建立

    本文以某直升机为例进行雷电间接效应仿真分析,该直升机模型如图1所示,其尺寸为22m×18.6m×7.8m,为了使模型更精确使用FEKO的挡风玻璃扩展功能能将驾驶窗前面设置为三块玻璃,厚度为3mm,如图1(a)。为了分析雷电对直升机表面天线影响,设置了如图1(a)的3个单极子天线;为了分析雷电对机载设备的影响,在机舱内设置了一根线缆,线缆两端分别端接50Ω负载,线缆的类型为RG58的同轴线,线缆分布如图1(b)所示。

图1 阿帕奇直升机模型 

图1 阿帕奇直升机模型

2.3激励源设置

    根据SAE-ARP5416-2005试验试验方法的规定,全机闪电间接效应试验需进行大电流脉冲注入试验和扫频试验,本文主要针对大电流脉冲试验进行仿真分析。FEKO中设置外加电流为飞机雷电击中点,用与注入电流反相表示电流流出点。ASE-ARP5412中对雷电试验波形进行了规定,间接试验需采用试验波形A,D和H进行脉冲注入试验,由于直升机属于低空飞行器,因此选择波形A进行电流脉冲试验的仿真分析。该波形是一个双指数形式,其数学表达式为:

公式1 波形表达式 

公式1 波形表达式


其时域波形如图2所示。由图可知雷电注入6.4μs时达到峰值电流200KA,在69μs时,波形A衰减到最大值的一半100KA。

图2 雷电波形A的时域曲线 

图2 雷电波形A的时域曲线

2.4仿真设置

    根据图2可知当波形A的电流幅值下降到峰值的50%时的时间为69us,而ASE-ARP5412中规定试验波形A的整体时间为500us以内。由于仿真时间的需要至少应包括雷电流流过整个机身的时间t1=73ns和波形A的半宽度时间t2=69us。为了能准确、快速的完成直升机间接雷电的仿真分析,本文仿真设置时间设为200us,仿真频率范围为1kHz-50MHz。由于本文主要分析直升机在飞行过程中,遭受间接雷电的影响,因此设置直升机的六个面均为开放边界条件。其中雷电注入路径如图3所示,电流注入点设在机头部位,从机尾流出。

图3 雷电注入路径 

图3 雷电注入路径

3 仿真结果分析

3.1机身表面电流分布

    根据图2可知在6.4μs的时候,波形A到达峰值电流200KA,这时注入直升机上的电流最大,机身表面电流密度也达到最大。在69μs时,波形A到达100KA,即最大电流的一半,这时注入直升机上的电流比6.4μs时有所减小,仿真结果如图4所示,直升机表面电流在雷电注入点和分离点处达到最大,其数值达到几百KA/m2。并且在路径中,机头、驾驶舱,发动机整流罩、主旋翼、尾梁、尾翼等位置处的表面电流也比较大,驾驶舱玻璃上边缘处有小电流流过。同时,雷电流流经的直升机表面棱边(如:驾驶舱、尾梁的边沿)的电流密度也比较大,其数值达到100KA/m。

图4 直升机表面电流分布 

图4 直升机表面电流分布

    为了更加清晰、准确的了解机身不同表面电流密度的分布情况,得到直升机表面三个天线的耦合电流分布如图5所示。其中antenna_1、2、3分别为机头、机中、机尾单极子天线。

图5 3个单级子天线电流分布 

图5 3个单级子天线电流分布



从图中可以看出,在雷电注入点处的电流密度最大,如耦合到机头部天线的电流最大(蓝色);直升机机身中部天线耦合最小(红色);机身尾部天线耦合因靠近尾部雷电流出点,电流进行汇聚导致比中部耦合电流大。从上面的分析和电流分布图可以看出,除了在雷击的附着点和分离点处的电流比较大以外,在发动机整流罩,尾梁、窗户等位置的边沿的电流比较集中。因此,在进行直升机设计时,应该避免不连续棱边和不规则表面结构的设计,这些地方是进行雷电防护的重要位置。

3.2舱外电磁环境

    雷电是一种高能量的脉冲信号,它在空间中分布的电磁场对直升机电子系统造成了严重的威胁。在雷电击中飞机时,6.4μs(最大)和69μs(衰减至一半)时直升机表面附近的电场强度分布图如图6所示。

图6 直升机舱外电场分布图 

图6 直升机舱外电场分布图

    根据仿真结果可知,在此路径下直升机注入电流最大时表面附近的瞬时电场强度达到了上百MV/m。从图中可以看出,直升机附近的电场强度主要集中在机头、旋翼、尾翼等结构比较突出的地方。在雷击条件下,直升机附近的电场可以看成是机体表面电流产生的近场电磁环境,而近场电场的分布与直升机表面电荷的分布有着紧密的关系,而机身表面比较突出的结构,如机翼、尾翼、机头等尖端处容易集聚大量的电荷,从而造成尖端附近的电场环境比较恶劣。

    图7显示了6.4μs和69μs时机身表面电流和直升机外部磁场的分布图。从图中可以看出,机身表面电流密度越大,其附近的磁场强度也就越大,磁场强度的分布与电流密度的分布保持一致,这是因为电流是磁场产生的源。

图7 直升机空间磁场分布图 

图7 直升机空间磁场分布图

    磁场强度最大值达到了170KA/m,其中机头,发动机整流罩等关键部位附近的磁场环境比较严峻。机头、机尾是雷电的注入和输出点,因此这两处附近的磁场强度明显比较大。

    由机舱外部电磁环境的分析可知,直升机表面电场分布与机身不连续的表面结构有关;而磁场强度分布与表面电流的分布有关,而表面电流密度的分布不仅与雷电路径相关,同时和机身表面的结构有关系。因此直升机结构的对舱外电磁环境有着重要的影响,良好的直升机结构设计能最大限度的使雷电流均匀分布在直升机表面,从而使得舱外电磁环境更加良态。

3.3舱内线缆

    直升机由于装配、焊接等操作产生了的孔缝和窗、门等部件严重降低了直升机的屏蔽效能。预先设定的直升机内部线缆在100us以内耦合电流分布如图8所示。

图8 线缆耦合电流 

图8 线缆耦合电流

    由于模型的简化,只是将直升机驾驶舱设置为玻璃而机身采取精密的全封闭机身进行仿真,线缆耦合的电流比实际略小。这是因为实际情况中驾驶舱和机身的孔缝是不能避免的,对舱内线缆影响更加剧烈。

    因此,对于雷电间接效应的防护,应该尽量避免机身孔缝等结构出现,并对关键部位采取加入屏蔽隔板的措施增强电磁屏蔽效能。

4 结论

    本文根据SAE-ARP5416和ASE-ARP5412相关标准,使用FEKO的强大的时域分析功能对直升机机身表面电流分布、舱外电磁环境、舱内电磁环境和线缆感应电流进行了仿真分析。得到直升机表面电场分布与机身不连续的表面结构有关;而磁场强度分布与表面电流的分布有关。根据仿真结果给出了影响机身表面电流分布和电磁环境的一些主要的因素,仿真结果表明该方法能有效的模拟直升机在飞行过程中的雷电间接效应,为直升机的优化设计和雷电防护提供了依据。


双击此处添加文字
技术文章