龚德志 王新栋 叶聪杰 / GONG Dezhi WANG Xindong YE Congjie
(上海飞机设计研究院,上海 201210)
多梁式中央翼盒下壁板压缩稳定性研究
龚德志 王新栋 叶聪杰 / GONG Dezhi WANG Xindong YE Congjie
(上海飞机设计研究院,上海 201210)
研究多梁式翼盒加筋壁板在压缩载荷作用下的稳定性。针对端部支持、侧边支持、本身曲率以及上述因素的综合作用对加筋壁板压缩失稳临界应力的影响进行分析,对目前文献资料中关于加筋壁板压缩稳定性临界应力计算公式中端部支持系数进行适当修正,以得到适合的壁板屈曲应力。研究发现,端部夹持、侧边梁支持和蒙皮自身曲率对加筋壁板的压缩稳定性有较大影响,对于蒙皮较厚的加筋曲板(如机翼壁板),建议的等效端部支持系数为1.5~2.0。
多梁式翼盒;加筋板;屈曲;端部支持系数
民用飞机的中央翼结构设计常采用多肋式或多梁式的结构布置形式。譬如空客系列飞机的中央翼盒常采用多肋(桁架肋)式布置,而波音系列飞机的中央翼盒则习惯采用多梁式布置。在工作状态下,中央翼盒通过上、下壁板的展向拉、压承受机翼弯曲引起的翼根弯矩。承压结构失稳而引起结构失效是制约民用飞机中央翼盒强度的一个主要因素。无论多肋式还是多梁式中央翼盒,上壁板均有较为密集的地板梁支持,具有较高的展向压缩设计值;下壁板则一般按承受展向拉伸载荷设计,但对下壁板的稳定性设计也是不容忽视的。特别是多梁式中央翼盒,下壁板从端肋平面到龙骨梁之间处于无支持状态,在某些机翼下弯工况下,其稳定性问题将变得十分突出。
国内外许多专家学者对组合式加筋壁板的压缩稳定性进行了大量的理论和实践研究,给出了较为完善的分析方法,形成了设计手册,进而指导飞机设计员进行设计工作。然而,目前的分析方法对多梁式翼盒加筋壁板的稳定性计算要么不适用,要么过于保守。在文献[1]第446-453页和文献[2]第88-99页中给出的图表要求加筋壁板的加筋比(筋条截面积与蒙皮截面积之比)为0~0.4,而文献[4]第556页则要求民用飞机壁板设计的加筋比不应小于0.5。同时,文献[1]第446页和文献[2]第137页中给出的方法均未考虑本身曲率对加筋壁板稳定性的影响。
本文对多梁式中央翼盒下壁板在压缩载荷作用下的稳定性进行研究,针对不同的边界条件以及结构形式的加筋壁板开展有限元计算,并对工程计算公式的参数进行修正,在安全的前提下,减少其保守量,为飞机壁板设计提供参考。
相对于多肋式中央翼盒,多梁式中央翼盒除了前梁和后梁外,其内部至少还布置一根沿展向的翼梁(图1(b)),且布置有一根翼肋或没有翼肋,如图1所示。此外,上壁板上表面沿航向一般布置有多根地板梁(图1(c)),下壁板下表面则沿航向布置有龙骨梁(图1(d)),由长桁和蒙皮组成的加筋壁板构成了翼盒的上、下壁板。
下壁板长桁端头通常斜削并加宽,进而插入翼盒端肋齿形缘条与对接带板之间,如图2所示,即下壁板端部被肋缘条和对接带板夹持,这样的连接结构具有重量轻、可靠性高和固有破损—安全特性等优点,在民用飞机结构上被广泛应用。
图2 下壁板结构及端部连接示意图
加筋壁板承受轴向压缩载荷作用示意图如图3所示。根据文献[1]第469页,加筋壁板柱压缩强度曲线分为短柱区、中长柱区和长柱区。加筋壁板从肋平面到龙骨梁位置较长,单根长桁与蒙皮组成的壁板长细比较大,加筋壁板有效长度L′与壁板回转半径ρ的比值通常落在长柱区,壁板的失稳形式为Euler柱失稳(弯曲失稳且临界应力依然处于材料线弹性范围内[3]),其临界应力由式(1)计算。
图3 壁板受压示意图
式中,Ec为材料压缩弹性模量,MPa;L′为壁板的有效长度,mm;ρ为壁板的横截面回转半径,mm。
式(1)中壁板的有效长度和横截面回转半径分别由式(2)和式(3)计算。
式中,L为壁板在两端支持结构之间的长度,mm;C为壁板的端部支持系数;I为壁板截面惯性矩,mm4;A为壁板截面积,mm2。
采用文献[1]第469页中的公式来计算加筋壁板的压缩失稳临界应力尽管简单方便,但该文献中并未给出壁板的端部支持系数C的取值。文献[4]第552页中建议C=1.5,然而通过计算及试验发现其结果较为保守,某试验中壁板由公式计算得到压缩失稳临界应力远低于有限元计算值,亦低于工作应力值,而壁板在试验载荷下并未发生失稳,如图4所示。另外,各方法均未给出加筋壁板端部支持系数C的计算方法,并且不能合理地考虑多梁式中央翼盒中展向梁对壁板的边界支持作用。
图4 试验中壁板压缩应力归一化图表
下文将通过有限元分析计算,分别研究端部约束、壁板曲率以及侧边支持三个因素对加筋壁板压缩失稳临界应力的影响。
文献[4]第556页中的计算方法要求壁板加筋比(长桁截面积与蒙皮截面积之比)不大于0.5。然而,为了得到较高结构效率,理想状态下承压加筋壁板的加筋比应为1.2左右,见文献[4]第553-554页。因此,本文折中考虑取壁板加筋比为0.8左右。
单根长桁与蒙皮组成的加筋壁板截面如图5所示。考虑到目前单通道机型中央翼壁板两根梁之间一般布置5~7根长桁,模型中长桁数目分别取5、6和7,如图5所示;为了便于建模并减少计算量,取加筋壁板截面参数t1=t2=t3=t4=t,b1= 20 mm、b2= 45 mm、b3= 60 mm、b4= 160 mm,其中,t分别取3 mm、4 mm、5 mm、6 mm和7 mm;壁板总长L0= 3 600 mm,龙骨梁宽度b0= 500 mm,端肋与龙骨梁侧边长度L= 1 550 mm;加筋曲板曲率半径R分别取5 000 mm、7 000 mm、10 000 mm和15 000 mm。加筋壁板长度定义见图6。本文计算的壁板结构几何参数见表2。壁板材料取疲劳性能较好的2000系列铝合金,螺栓材料为Ti-6Al-4V,材料属性见表1。
表1 金属材料性能数据
图5 单根长桁和蒙皮组成的加筋壁板截面图 图6 加筋壁板长度定义
模型编号t1mmb1mmt2mmb2mmt3mmb3mmt4mmb4mmAmm2加筋比Imm4ρmm133034536031608760.8345211016.402430445460416011640.8262357116.683530545560516014500.8180677116.984630645660616017340.81100249617.295730745760716020160.80121151317.61
加筋壁板的有限元建模及计算通过MSC.Patran/Nastran平台进行,有限元模型采用四边形CQUAD4单元以及极少量的三角形CTRIA单元模拟,连接螺栓采用CWELD单元模拟。图7为加筋壁板有限元模型,图8为两侧有梁支持的加筋壁板有限元模型示意图,图9为加筋壁板有限元模型端部夹持方式,图10为带曲率加筋壁板有限元模型示意图。
图7 加筋平板有限元模型
图8 侧边梁支持加筋平板有限元模型
图9 端部夹持加筋平板有限元模型
图10 加筋曲板有限元模型
载荷和边界条件为:对壁板龙骨梁站位约束蒙皮面外位移(Y向位移);加筋壁板一端约束X、Y和Z向位移, 另一端约束X和Y向位移并沿Z向施加0.36 mm的压缩位移载荷。
失稳临界应力:σcr=λσ0,λ为由有限元屈曲分析得到的屈曲特征值。
加筋壁板等效端部支持系数C通过Euler式(1)反推式(4)得到。
图11 加筋平板屈曲模态
图12 加筋平板压缩屈曲应力
图13 加筋壁板等效端部支持系数
分别对由5、6和7根长桁组成的加筋平板在仅有龙骨梁支持状态下受轴向压缩载荷时失稳临界应力进行计算。加筋壁板屈曲符合弯曲(Euler)柱强度屈曲的特点(垂直于板面弯曲),如图11所示。失稳临界应力见表2。根据式(2)计算得到壁板等效端部支持系数如图12和图13所示。通过图12和图13可看出加筋壁板其等效端部支持系数随长桁数增加而变小。
(1)蒙皮厚度较小的加筋壁板等效端部支持系数小于1,由于壁板加载偏心造成。而蒙皮与长桁之间有一定的高度,分别对长桁和蒙皮进行约束,限制加筋壁板的转动自由度,厚度较大加筋壁板等效端部支持系数大于1;
(2)加筋壁板等效端部支持系数在厚度等于6 mm处出现拐折,由于前面蒙皮先于长桁失稳,而后面则是长桁先于蒙皮失稳引起。
本小节所考察的加筋壁板中,将侧边长桁弹性支持改为梁弹性支持,梁为壁板两侧提供的面外支持刚度高于长桁,并有效增大壁板的扭转刚度,从而影响加筋壁板的压缩屈曲强度。
选取梁厚度为壁板蒙皮厚度的1.2倍,分析得到侧边梁支持加筋平板失稳临界应力如图14所示;其等效端部支持系数如图15所示。
图14 侧边展向梁支持加筋平板失稳临界应力
图15 侧边展向梁支持加筋平板等效端部支持系数
从上图可看出,翼梁对加筋壁板端部支持系数的影响较为明显,较图12和图13提高22%左右;壁板长桁数量越少,翼梁支持对壁板等效端部支持系数影响越大。
端部夹持(如图9所示)能够增加对壁板面外弯曲的约束。本小节所述加筋壁板取长桁数为7,夹持结构的厚度分别取为加筋壁板蒙皮厚度尺寸的1.0、1.2、1.4、1.6、1.8和2.0倍。
通过计算得到两端简支加筋平板的失稳临界应力,如图16所示;其等效端部支持系数如图17所示。
图16 端部夹持加筋平板屈曲应力
图17 端部夹持加筋平板等效端部支持系数
通过图16和图17可看出,夹持结构厚度由加筋壁板厚度的1.0倍逐渐增加到2.0倍的过程中,壁板端部支持系数增加呈现放缓趋势。与图11进行对比,端部夹持对提高加筋壁板屈曲应力作用非常明显。
飞机机翼壁板通常带有一定的曲率,曲率能够增大壁板整体的扭转及弯曲刚度。
分别对曲率半径R为5 000 mm、7 000 mm、10 000 mm和15 000 mm的含有7根长桁的加筋曲板进行屈曲分析。计算得到失稳临界应力如图18所示,等效端部支持系数如图19所示。
图18 加筋曲板屈曲应力
图19 加筋曲板等效端部支持系数
从图18和图19可看出,加筋壁板等效端部支持系数随着曲率增加而增大。曲率半径为15 000 mm的加筋壁板,其等效端部支持系数比平加筋壁板高10%以上。曲率对加筋壁板的失稳临界应力影响较为显著。
通过以上计算分析,端部夹持、两侧梁支持、曲率半径因素均对含有7根长桁的加筋壁板的压缩失稳临界应力产生影响。综合以上因素对加筋壁板进行屈曲应力计算分析。
根据参考文献[5]第229-232页,插入件与被插入件在挤压应力,拉伸应力和螺栓直径相同的情况下,插入件与被插入件载荷传递系数比值约为0.89/0.735=1.21。因此,结构设计要求被插入件厚度为插入件的至少1.2倍。因此,对于端部夹持件厚度取加筋壁板厚度的1.2倍。
计算得到失稳临界应力如图20所示,等效端部支持系数如图21所示。
图20 加筋壁板屈曲应力
图21 加筋壁板等效端部支持系数
各因素对加筋壁板的屈曲的综合影响较大,曲率半径为15 000 mm的加筋壁板等效端部支持系数C最小,为1.76,但相比两端简支的加筋平板有很大的提高。
分别对侧边支持、端部夹持、曲率三种因素单独进行分析,以及综合以上因素对加筋壁板稳定性的影响进行了分析研究,得到结论如下:
(1)端部夹持对加筋壁板压缩稳定性影响最大,端部夹持件厚度取壁板厚度的1.2~1.4倍较好;
(2)加筋壁板侧边展向梁弹性支持对加筋壁板压缩屈曲影响较大,为了提高加筋壁板压缩失稳临界应力,可适当增加展向梁数量;
(3)曲率半径对加筋壁板压缩屈曲影响明显,特别是当加筋壁板曲率半径较小时尤为明显;
(4)随着加筋壁板厚度的增加,加筋壁板端部支持系数可取1.5~2.0。
[1] 《飞机设计手册》总编委会.飞机设计手册第9册:载荷、强度和刚度[M].北京:航空工业出版社,2001.
[2] 崔德刚. 结构稳定性设计手册[M].北京:航空工业出版社,1996.
[3] 苟文选,金保森,卫丰. 材料力学[M]. 西安: 西北工业大学出版社,2001.
[4] 牛春匀. 实用飞机结构应力分析及尺寸设计[M]. 北京:航空工业出版社,2009.
[5] 郑晓玲,等.民用飞机耐久性与损伤容限设计手册(上册)[M].北京:航空工业出版社,2003.
Stability Study of Buckling on Stiffened Lower Panel for Center Wing Box with Multi-beam
(Shanghai Aircraft Design and Research Institute, Shanghai 201210, China)
This paper studies the buckling of stiffened lower panel of center wing box with multi-beam under compressive loading. For the effect of end support, lateral support and curvature for the buckling of panel, the end support fact of critical compressive stress formula of stiffened panel in current design manual was analyzed, and proper buckling stress of stiffened panel was obtained. The result shows that the end support , lateral support and curvature of skin make a significant effect on the stabilization of stiffened panel under compressive loading, and for the stiffened panel with thick shin, the end support fact is advised to be 1.5~2.0.
multi-beam wing box; stiffened panel; buckling; end fixed fact
V224
A
10.19416/j.cnki.1674-9804.2017.04.017
龚德志男,硕士,工程师。主要研究方向:飞机结构强度设计分析;Tel:20865865;E-mail:gongdezhi@comac.cc
王新栋男,硕士,工程师。主要研究方向:飞机结构强度设计分析;Tel:20865923;E-mail:wangxindong@comac.cc
叶聪杰男,硕士,工程师。主要研究方向:飞机结构强度设计分析;Tel:20865919;E-mail:yecongjie@comac.cc