冯振宇,周书婷,李恒晖,解江,李翰
(中国民航大学 天津市民用航空器适航与维修重点实验室,天津 300300)
近年来,随着恐怖袭击事件的不断增多,民航客机经常成为恐怖袭击的对象。而炸弹作为其常用的袭击破坏手段,一旦引爆将对飞机结构造成致命性的破坏,例如1988年的洛克比空难。2008年10月28日,美国联邦航空管理局(Federal Aviation Administration,简称FAA)颁布FAR 25-127修正案[1],要求制造商针对最大合格审定客座量大于60或起飞总重超过100 000 lb(45 359 kg)的飞机,必须设计一个“最小风险炸弹位置[2](Least Risk Bomb Location,简称LRBL)”,用于放置发现的可疑装置,通过结构和系统的综合设计,使得关系飞机安全的关键结构和系统在其一旦爆炸后能够得到最大程度的保护,从而有效提高飞机的安保水平。
针对此适航要求,FAA颁布了咨询通告AC 25.795-6[3],给出了推荐的LRBL符合性方法。该咨询通告指出,表明符合性时既要考虑LRBL周围结构和系统的设计与布局,即必须计算分析或/和试验评估机体结构的承载能力确保飞机结构及系统的完整性,将损失降到最小,也要评估爆炸发生后的次生影响以及应对措施,即在设计LRBL时,应当考虑到包括碎片被吸入发动机及乘员损伤在内的次生影响,并指出当飞行中发现可疑装置时,可以使用爆炸包容系统来降低爆炸的影响。
然而,爆炸力学在民航领域的应用较少,AC 25.795-6中未提及制定LRBL符合性方法的关键依据,且国外的相关技术信息敏感,不对外公开。另外,我国正在研制国产大型客机,飞机安保设计也在力求达到FAR 25-127修正案的水平。即在该修正案生效后申请取证的运输类飞机,都必须考虑这一适航要求。但是我国最新的CCAR 25-R4对LRBL适航要求暂无修订,更没有专门的针对LRBL的符合性研究。因此,适航当局与制造商在面对这一新要求时,迫切地需要开展相关的研究来支持适航审定与型号设计。
本文基于LRBL适航要求及相关的咨询通告,针对飞机结构的爆炸冲击响应、人体冲击损伤、爆炸碎片及抗爆容器等方面的研究进展进行综述,并介绍中国民航大学天津市民用航空器适航与维修重点实验室的相关技术研究,指出进一步需要研究的问题,旨在为LRBL的适航审定与型号设计提供有效的指导。
AC 25.795-6指出LRBL必须避开飞机的关键结构及系统,且需考虑该位置结构的破坏对飞机结构完整性的影响。由于爆炸物对飞机的破坏作用主要来自于炸药爆炸后产生的高速高压冲击波,特别是对完全增压的飞机,爆炸造成的损伤将越加显著。
预测爆炸对飞机结构的冲击响应的早期尝试是由R.C.DeHart等[4]进行的,R.C.DeHart等基于弹性动力分析理论,对飞机结构各构件的动态响应进行了理论推算,为机身的设计以及机身极限承载的确定提供了理论依据。到20世纪70年代,FAA针对特定机型(B727-100、DC-9、B747-100和L-1011)开展了最小风险炸弹位置的评估研究[5-8],说明了诸如厨房、应急舱门、后承压框、行李架等不同LRBL设计方案的设计原则、特征及优缺点;提出了LRBL的估测流程为:在选定的位置放置爆炸源,计算作用于机身内部的压力冲量,将动载简化等效为静载并与蒙皮壁板的承载能力进行比较,当等效静载超过极限载荷时,失效发生。这一方法流程是采用等效近似的方法进行的理论计算,尽管其为LRBL的评估提供了有效的指导,但已无法满足目前计算精度与安全性的需要。
洛克比空难后的机身损伤形貌图如图1所示[9],在1988年洛克比空难之后,美国和欧洲加大了航空安保领域技术研究的投入。开展了140次内爆载荷下的机身易损性试验(其中窄体客机98次,宽体客机42次)、200多次支持性试验(主要是研究乘员的手提行李以及货舱的托运行李对爆炸效应的抑制作用)以及104次减缓爆炸结构设计的验证试验(包括客舱抗爆衬垫、货舱衬垫、宽体客机强化的货舱集装箱和验证最小风险炸弹位置的程序);建立了简单的圆筒结构以及B-52、B707等全尺寸飞机模型,分析了冲击波和机身的相互作用,提出了商用飞机“易损性分布图”的概念,如图2~图3所示[10-12]。可以看出:确定位置下的造成飞机无法安全着陆的最小炸弹当量。
(a) 飞机蒙皮损伤
(b) 飞机机身损伤图1 洛克比空难后的机身损伤形貌图Fig.1 Aircraft fuselage damage morphology after theLockerbie bombing
图2 B-52机身易损性分布图Fig.2 Vulnerability map of B-52
图3 B-707机身易损性分布图Fig.3 Vulnerability map of B-707
20世纪末,英国民用航空管理局(Civil Aviation Authority,简称CAA)的Hardening项目[13]发展了数值仿真计算的能力,并在此研究基础上,提出了多种飞机结构强化技术,例如局部加强机身蒙皮或框段、强化货舱行李包容系统及衬垫、加强系统保护以及增加货舱行李包容系统与飞机蒙皮的距离等,以保障飞机的安全运营。
进入21世纪,欧洲的多个机构,包括高校、公司、科研机构,在欧盟第六框架项目的支持下发起了VULCAN项目,旨在发展新材料、新方法/策略,来防止因爆炸造成的灾难性事故。VULCAN项目采用“积木式方法”(如图4所示)进行了研究,从较小的结构元件(试样、加强板、装配体、加压蒙皮)到全机身结构,研究了不同材料(金属、复合材料以及GLARE)的结构在爆炸载荷作用下的结构响应。例如,C.Soutis等[14]对GLARE层板进行了爆炸仿真模拟,研究了GLARE在爆炸载荷作用下的非线性动态响应,并与G.S.Langdon等[15]的试验结果进行了对比,得到了较好的吻合度。T.Kotzakolios等[16]和A.Dacko等[17-18]建立了两个典型商用运输类飞机的机身结构有限元模型,分别采用铝和GLARE材料,针对不同的炸药位置进行了数值模拟,得到了结构不同损伤位置下的不同损伤程度。另外,VULCAN项目还在设计和材料应用方面提出了强化策略,例如,对于铝合金飞机结构,提出在蒙皮、桁条和框架之间增加螺栓连接,如图5(a)所示;对于复合材料飞机结构,提出将长桁与蒙皮先进行粘合,然后再在框架与蒙皮之间增加螺栓连接,如图5(b)所示。
图4 积木式研究图Fig.4 Building block research
(a) 铝合金飞机结构强化
(b) 复合材料飞机结构强化图5 强化策略Fig.5 Hardening strategy
AC 25.795-6指出在设计LRBL时,除了考虑结构本身的设计以外,还应考虑包括乘员损伤在内的次生影响。炸药在爆炸时对人体的作用及伤害是多重效应共同作用的结果,包括冲击损伤、破片损伤、加速-碰撞损伤、热损伤等。其中,冲击损伤的研究较为关键。经过几十年的发展,冲击损伤的研究经历了动物试验、损伤机制分析、数学建模以及计算机数值模拟研究几个阶段。
炸药爆炸产生的冲击波将造成人体器官的创伤,尤其是肺。其损伤程度取决于炸药的威力、受伤者离爆炸点的距离、爆炸发生的空间性质以及炸弹碎片的性质。针对冲击损伤研究,最为著名的是I.G.Bowen,J.H.Stuhmiller等和H.Axelsson等基于动物(羊)试验数据,提出的损伤模型。三种损伤模型的对比如表1所示。
表1 损伤模型的对比Table 1 Comparison of injury models
Bowen模型[19]集中在肺损伤,其根据爆炸试验建立了爆炸创伤病理学数据库,针对70 kg人体不同方向的受载绘制了著名的Bowen创伤曲线,如图6所示。Bowen创伤曲线的局限之处在于只针对理想的爆炸冲击波。
(a) 人体与爆炸波阵面平行
(b) 人体与爆炸波阵面垂直
(c) 人体靠近反射面图6 Bowen创伤曲线Fig.6 Bowen trauma curve
J.H.Stuhmiller等[20]于1984年在人体躯干模型周围布放了压力传感器研究自由场冲击波作用时人体表面的载荷变化过程。随后在1996年又建立了人体胸部简化的动力学模型[21],考虑了由于爆炸负荷 产生的荷载,由于肺的体积压缩产生的内部压力以及由于胸部运动产生的压缩波,以冲击波对肺组织所作的功为参考量,获得了肺组织各个创伤等级所占的比率,推导公式[21]如下:
p(t)≅p0+ρ0c0v
(1)
(2)
(3)
(4)
式中:p为压力;p0为正常肺的压力;ρ0为正常肺的密度;c0为正常肺的声速;v为活塞速度;γ为比热比;x为位移;m为质量;L为肺跟胸腔面积的比值;w*为正常功;w为总功。
2004年,J.H.Stuhmiller研制出了爆炸冲击波对肺创伤的评估软件Injury,此评估软件的开发对于人体爆炸冲击损伤的研究起到了巨大的推动作用,所需输入参数是爆炸载荷、物种、冲击次数、大气压力以及起始和结束时间,计算结果是肺损伤的概率,定义的5个伤害级别分别为严重、中度、轻微、留有痕迹和无痕迹。
Axelsson的损伤模型[22]包括对所有易受伤害的内部器官/系统的损伤。H.Axelsson基于J.T.Yelverton的试验研究结果,建立了由弹簧、质量、阻尼器系统构成的简化的人体胸部动力学模型,如图7所示。该模型将冲击波测量数据作为输入参数,通过计算胸壁运动速度及肺腔中的压力变化,评估了胸腔内各个器官的创伤程度。但是Axelsson方法只能了解冲击波对生物胸部内部的生物器官创伤程度的范围,并没有具体分析爆炸波对生物组织产生各个等级创伤所占比率。模型方程如下:
图7 Axelsson模型Fig.7 Model Axelsson
式中:A为有效面积;M为有效质量;V为肺的初始气体体积;X为位移;C为阻尼系数;K为弹簧系数;P0为环境压力;P(t)为即时超压;γ为肺中气体的多变指数。
除了上述三种著名的人体损伤模型研究,国内外其他学者的工作也在积极开展。周杰等[23]通过去除模型中的撞击系统及完善胸部左-右部分的数学模型,获得了爆炸作用下的人体胸部动力学模型。将试验数据加载到模型的前、后、左、右位置作为边界条件,计算分析胸部各个方向的运动速度、位移、肺腔内压及肺组织中的压力波变化规律,利用归一化功与肺组织创伤的关系,评估了肺组织的创伤。W.E.Baker等[24]以人垂直站在或水平躺在平整地面上时,冲击波传播方向与身高方向垂直,周围无障碍物(这是最典型的暴露情形)这一假定为基础,根据Lovelace基金会的研究成果,推导了关于肺伤害人员的死亡率曲线。
相比于动物试验的高成本和低精度,数值模拟研究逐渐体现了其优越性。周杰等[25-26]进一步采用数值模拟方法,分析了人体胸部在自由空间爆炸场中受冲击波作用的力学过程,参照Bowen损伤曲线评估出了肺的损伤状态,并改进了T.E.Lobdel[27]的人体胸部创伤模型,推导了模型的动力学方程,选择了合理的计算方法。康建毅等[28]通过CT扫描完成了人体胸部三维表面模型的建立,探讨了冲击波载荷下的人体胸部周围压力的分布及肺脏的应力分部规律。A.I.D’yachenko等[29]建立了一维有限元肺模型,并选用粘弹性材料模型研究了冲击波传播造成肺的损伤。
AC 25.795-6指出在设定LRBL时,除了结构本身的设计以外,还需要考虑包括碎片造成的飞机损伤在内的次生影响。因为抛射碎片具有非常大的速度和能量,如果击中飞机关键结构(例如机翼),或是被吸入发动机,将会造成非常严重的后果。
针对爆炸产生的碎片抛射危害的研究,目前多集中在碎片的产生、抛射、击中及对目标的损坏等方面,旨在明确爆炸碎片的抛射速度及飞行轨迹。U.Hauptmanns等[30-31]通过Monte-Carlo方法对碎片抛射的随机性进行了分析,求取了碎片对水平物体的破坏概率。R.Pula[32]和G.Gubinelli等[33]基于二维抛射的随机轨迹方程建立了各自的破坏概率模型,考虑了简化的目标物体形状对破坏概率的影响。二维抛射碎片飞行轨迹如图8所示[32]。Q.B.Nguyen等[34]基于碎片抛射轨迹和目标物体的三维方程建立了破坏概率模型,其认为当受撞击的目标物体的残余壁厚小于某个临界值时,结构就会发生破坏,三维抛射碎片轨迹抛射简图如图9所示。尽管起步较晚,国内学者也在积极开展研究,张新梅等[35]推导了爆炸碎片抛射速度的计算公式,构建了爆炸碎片的飞行加速度方程,确定了求解爆炸碎片飞行轨迹及抛射速度的计算方法。韩晨曦等[36]综述了过程设备中由于爆炸可能产生的碎片的抛射危害概率,对碎片的产生、抛射以及对目标碰撞过程中各不确定参数所符合的概率分布进行了整理分析。赵敏等[37]通过数值模拟的方法,得到了爆炸储罐在爆炸后所产生的碎片轨迹。朱东[38]研究了储罐在遭受撞击和爆炸荷载下的动力响应。
图8 二维抛射碎片轨迹示意图Fig.8 Schematic of trajectories of two-dimensionalprojectile fragments
图9 三维抛射碎片轨迹示意图Fig.9 Schematic of trajectories of two-dimensionalprojectile fragments
AC 25.795-6中说明可以采用炸弹包容系统(Containment System),即抗爆容器,来达到降低爆炸影响的目的。抗爆容器能限制系统内部炸药爆炸产生的冲击波、爆轰产物和破片的作用范围,有效保护人员、设备和环境的安全。最初,抗爆容器以传统的金属抗爆容器为主。随着新材料的开发及1946年纤维缠绕技术的出现,轻质高强的新型复合材料抗爆容器开始发展。相比传统的金属抗爆容器,复合材料容器质量轻、强度高、刚性好、抗爆性能佳并且设计灵活。经过几十年的发展,到20世纪90年代,经过各国学者的努力,抗爆容器在理论研究、试验研究及应用上都取得了较大的进展。我国对于抗爆容器的研究起步较晚,从20世纪80年代开始至今,包括中国工程物理研究院、北京理工大学、国防科技大学以及西北核技术所等都在开展相关方面的研究工作。
针对航空应用的抗爆容器的研究,始于20世纪90年代。当时,FAA评估认为当时的货舱行李集装箱的抗爆能力有限,发起了货舱集装箱的加强技术研究。CAA[13]开发了加强的行李集装箱来减弱飞行中爆炸产生的影响,如图10所示。但其重量大,成本高且不适用于大多数窄体客机,最终没有得到广泛的使用。
图10 强化的行李集装箱Fig.10 Hardened luggage container
21世纪,D.Zangani等[39]提出了采用新型纺织材料和复合材料组合使用的方法实现了行李集装箱的高柔,高强及灵活的可设计性,并且其质轻,可以同时在宽体和窄体客机中使用。2015年,科加雷姆航空公司的9268号航班发生炸弹爆炸,航空界开始更加关注增强飞机生存性的问题。对此,欧洲于2015年9月完成的一项名为“飞行包裹”的研究项目开发出了行李集装箱,如图11所示[40]。
图11 飞行包裹Fig.11 Flying bag
与20世纪90年代测试的重量大、成本高的爆炸硬化金属容器不同,“飞行包裹”采用了具有高强度、高抗冲击和防热特性的组合纤维织物,可以用在窄体和宽体飞机上。
可见,抗爆容器的设计特征为质轻、高强、能包容碎片,在内爆载荷下不发生撕裂。尽管发展迅速,可是目前抗爆容器的发展还存在尺寸效应(随着特征尺寸的增大,结构的强度和塑性变形能力下降,并有从韧性断裂转为脆性断裂的趋势,其存在会大大降低抗爆容器的极限承载能力)的问题,应用还不成熟,尚在研究之中。抗爆容器作为未来的设计发展趋势将会有很大的发展空间。
由于民航安全/反恐的迫切需要以及适航审定与符合性方法的技术基础需要,中国民航大学天津市民用航空器适航与维修重点实验室在民航局安全能力项目以及中央高校基金项目的支撑下,梳理了与最小风险炸弹位置相关的适航要求与典型案例;通过对典型航空铝合金进行动态力学性能试验与爆炸试验,获得其力学性能数据及破坏失效模式;通过对爆炸冲击载荷作用下的结构进行建模与数值计算,评估了材料的不同本构模型和失效参数的局限性与适用性;发展并建立了适用于LRBL适航要求的结构响应与失效分析方法。
为了确定航空铝合金材料的动态力学性能,采用霍普金森压杆(SHPB)试验测试铝合金2024-T3材料在不同应变率(10-3~103s-1)下的应力应变数据,基于获得的力学性能数据对Johnson-Cook本构方程中的参数进行拟合,从而建立材料的Johnson-Cook本构方程,为爆炸冲击载荷下结构的响应与失效分析提供可靠的材料模型。
为了明确航空铝合金薄板结构在爆炸冲击载荷下的动态响应及失效模式,针对固支铝合金圆板以药量及爆距为变量进行了爆炸冲击试验,试验布置如图12所示[41]。基于试验,利用LS-DYNA建立相应的数值仿真模型,分析了不同的ALE输运步算法、流固耦合方式、耦合点数量、平板网格尺寸、平板单元类型对于数值模拟结果的影响,以建立有效的数值模拟方法。最终通过经验证的数值模拟方法,在试验的基础上扩大分析矩阵,分析了近爆场下固支铝合金圆板在不同药量、不同爆距以及不同比例距离下的动态力学响应特性及变形失效行为,分析结果如图13所示。
图12 试验布置Fig.12 Test layout
(a) 不同药量对铝合金圆板变形失效模式的影响
(b) 不同爆距对铝合金圆板变形失效模式的影响
(c) 不同比例距离对铝合金圆板变形失效模式的影响图13 近爆场下固支铝合金圆板的动态响应Fig.13 Dynamic response of fully clamped circular aluminum alloy plate subjected to near field explosion
为了探索典型机身结构在内爆冲击载荷下的结构响应及失效模式,通过铝合金机身筒段模型(包含蒙皮、框、长桁、角片等典型航空结构件)试验分析了整体模型、连接件及连接方式对机身结构动态特性的影响。在试验的基础上,基于LS-DYNA,建立了机身筒段模型,研究了其在不同炸药量、不同爆距、以及不同炸药位置(机身中部或两框之间)下的结构动态响应及失效模式。机身筒段模型及分析结果如图14~图15所示[42]。
图14 机身筒段模型Fig.14 Aircraft fuselage substructure
图15 机身筒段模型分析结果Fig.15 Analysis results of aircraft fuselage substructure
针对目前国际反恐形势的日益严峻,美国以及欧洲发达国家在民用客机安保方面的研究已经积累了数十年经验,FAA颁布的针对民用客机安保措施的适航条款中提出了“最小风险炸弹位置(LRBL)”,给出了相应的咨询通告及修正案。并在过去的几十年中,开展了大量的相关技术研究。
考虑到目前我国正积极推进国产大型客机的研制,本文综述了飞机结构在内爆冲击载荷作用下的响应及失效模式的研究进展;总结了爆炸冲击波对人体的损伤问题,分析了人体器官在爆炸冲击波作用下的损伤机理和评估损伤模型;明确了爆炸碎片的破坏概率模型;针对国内外对抗爆容器的研究现状,总结了抗爆容器的设计特征及优缺点;简述了为获得典型机身结构在内爆冲击载荷作用下的动态特性开展的研究工作。
为了更全面、客观地研究飞机内爆问题、提高我国运输类飞机的安全技术水平。未来建议加强以下三个方面的研究:
(1) 加强客舱增压对爆炸冲击载荷下机身结构动态特性的影响研究
飞机在完全增压的情况下发生爆炸,损伤是最大的。在大多数情况下,当在飞机中发现可疑装置时,优先采取将客舱压差降低到0,作为即时的有效措施。客舱和外界空气之间的压差作为一个关键参数对于LRBL程序的有效性是非常显著的,因此建议针对客舱增压对机身结构动态特性的影响开展研究工作。
(2) 加强典型连接件在爆炸冲击载荷下的响应与失效模式研究
根据所开展的机身圆筒模型爆炸试验发现,连接件的失效是主要的形式之一,且与连接方式、工艺、连接件材料、规格等具有密切联系,是影响结构整体响应与失效的关键因素之一,在掌握平板爆炸响应与失效行为的基础上,有必要针对典型机身连接结构在爆炸冲击载荷作用下的响应与失效行为开展研究工作,进一步提高机身结构的抗爆能力。
(3) 加强包容系统的适航审定技术研究
包容系统是应急处置爆炸可疑装置的另一个有效方法。已有研究机构研发出具有抗爆功能的包容系统或类似设施设备。如何保证该类系统的有效性、应如何制定相关标准或规范以及制造方如何向局方表明符合性,是民航业应对炸弹袭击可能遇到的另一个问题。因此建议开展相关的工作,特别是包容系统抗爆符合性的适航审定技术研究工作。