連接翼布局的非線性屈曲分析

竇怡彬,張曉宏

(上海機(jī)電工程研究所,上海 201109)

連接翼布局的非線性屈曲分析

竇怡彬,張曉宏

(上海機(jī)電工程研究所,上海 201109)

以連接翼結(jié)構(gòu)為研究對(duì)象，通過非線性有限元方法對(duì)其進(jìn)行非線性屈曲分析，研究連接翼幾何參數(shù)和載荷分布對(duì)結(jié)構(gòu)非線性屈曲特性的影響規(guī)律。連接翼采用基于Co-Rotational格式的三角形平板殼單元進(jìn)行建模，采用弧長(zhǎng)法計(jì)算非線性屈曲的位移突跳現(xiàn)象以及后臨界變形模式。研究結(jié)果表明連接翼典型的載荷-位移曲線模式包括剛度軟化、位移突跳不穩(wěn)定和剛度剛化三部分。連接翼夾角增大后非線性屈曲載荷因子逐步降低，但存在一個(gè)臨界夾角，當(dāng)角度超過臨界角后載荷-位移曲線模式中的位移突跳不穩(wěn)定階段就不會(huì)出現(xiàn)。增加連接翼高度可提高結(jié)構(gòu)的非線性屈曲載荷。在連接翼夾角不變條件下，減小下翼展長(zhǎng)可有效提高非線性屈曲載荷，同時(shí)其載荷-位移曲線模式中的位移突跳不穩(wěn)定階段也將不存在。調(diào)整連接翼的載荷分布，將載荷更多地施加于上翼的同時(shí)減少下翼的載荷可提高結(jié)構(gòu)的非線性屈曲載荷。

連接翼； 非線性屈曲； Co-rotational格式； 位移突跳； 弧長(zhǎng)法； 剛度軟化； 剛度剛化； 臨界夾角

0 引言

連接翼布局結(jié)構(gòu)飛行器，也稱為盒型翼或Prandtl機(jī)翼，其結(jié)構(gòu)布局為兩個(gè)水平放置的機(jī)翼通過翼稍的垂直部分相連接而成。與單翼布局相比，連接翼的誘導(dǎo)阻力更小。文獻(xiàn)[1]比較了Prandtl最佳機(jī)翼系統(tǒng)和單翼系統(tǒng)(具有相同展長(zhǎng)和總升力)誘導(dǎo)阻力比值隨翼尖高度/展長(zhǎng)比(h/b)的變化，發(fā)現(xiàn)采用Prandtl布局可使誘導(dǎo)阻力減小20%～30%。這意味采用連接翼的飛行器具更高的經(jīng)濟(jì)性、更長(zhǎng)的飛行距離和更多的留空時(shí)間。自從WOLKOVITCH首次提出連接翼概念以來，國(guó)外對(duì)連接翼進(jìn)行了大量研究。早期的研究著重于氣動(dòng)力和結(jié)構(gòu)優(yōu)化等[2-7]。文獻(xiàn)[8]總結(jié)了之前關(guān)于連接翼氣動(dòng)彈性的研究成果并預(yù)測(cè)了可能面臨的挑戰(zhàn)，指出由于下翼會(huì)向上翼傳遞壓縮載荷，使連接翼的氣動(dòng)彈性模式較傳統(tǒng)單翼更為復(fù)雜，分析時(shí)須考慮機(jī)翼屈曲與氣動(dòng)彈性的耦合，并認(rèn)為屈曲是比顫振更重要的失穩(wěn)模式。文獻(xiàn)[9]認(rèn)為研究連接翼結(jié)構(gòu)的非線性屈曲行為有助于設(shè)計(jì)氣動(dòng)彈性靜/動(dòng)穩(wěn)定的飛行器，文獻(xiàn)[10]則在非線性屈曲分析的基礎(chǔ)上用耦合渦格法定常氣動(dòng)力對(duì)連接翼的靜氣動(dòng)彈性響應(yīng)及靜發(fā)散進(jìn)行了研究。類似于非線性屈曲不穩(wěn)定的位移-突跳概念，提出了非線性結(jié)構(gòu)靜氣動(dòng)彈性發(fā)散的突跳-發(fā)散概念。國(guó)內(nèi)對(duì)連接翼的研究起步較晚且多集中于氣動(dòng)、氣動(dòng)彈性分析，結(jié)構(gòu)優(yōu)化，以及可靠性等。文獻(xiàn)[11]對(duì)連接翼布局實(shí)現(xiàn)直接力控制的可能性進(jìn)行了初步探索。文獻(xiàn)[12]通過實(shí)驗(yàn)研究了五種不同布局形式連接翼的氣動(dòng)特性。文獻(xiàn)[13]研究了連接翼布局的縱向操縱特性，為飛行控制律的設(shè)計(jì)提供了基礎(chǔ)。文獻(xiàn)[14]基于NASTRAN的二次開發(fā)研究了連接翼的非線性顫振問題。文獻(xiàn)[15]采用CFD/CSD耦合方法研究了不同參數(shù)對(duì)連接翼靜氣動(dòng)彈性效應(yīng)的影響規(guī)律。文獻(xiàn)[16]對(duì)比了傳統(tǒng)布局機(jī)翼和連接翼的飛行載荷及顫振。文獻(xiàn)[17]采用同時(shí)考慮靜強(qiáng)度和顫振速度的混合約束剪裁方法設(shè)計(jì)了復(fù)合材料連接翼。文獻(xiàn)[18-21]對(duì)“鉆石背”連接翼的氣動(dòng)特性、靜氣動(dòng)彈性和顫振進(jìn)行了研究。文獻(xiàn)[22]以“鉆石背”連接翼翼尖變形最小為目標(biāo)函數(shù)，采用復(fù)形調(diào)優(yōu)算法和遺傳算法優(yōu)化了連接翼的結(jié)構(gòu)設(shè)計(jì)。文獻(xiàn)[23]采用首次超越破壞準(zhǔn)則建立“鉆石背”連接翼的動(dòng)力可靠性模型，并對(duì)模型的參數(shù)靈敏度進(jìn)行了研究。但目前還未見有連接翼非線性屈曲分析的相關(guān)報(bào)道。研究連接翼的非線性屈曲特性不僅可為結(jié)構(gòu)設(shè)計(jì)和結(jié)構(gòu)優(yōu)化等提供更真實(shí)有效的約束條件，而且也是連接翼飛行器穩(wěn)定性研究的一個(gè)重要方向。本文以一連接翼為研究對(duì)象，從上下翼間夾角、連接翼連接部位的高度、連接翼上下翼間的連接位置，以及連接翼間的載荷分配4方面研究了其非線性屈曲規(guī)律，為更好地理解連接翼結(jié)構(gòu)非線性行為提供參考。

1 非線性屈曲分析的弧長(zhǎng)法

用Co-Rotation有限元格式和三角形平板殼單元可得全局坐標(biāo)系中有限元靜力學(xué)方程為

(1)

(2)

(3)

2 非線性屈曲分析

2.1程序驗(yàn)證

先用標(biāo)準(zhǔn)算例驗(yàn)證自編非線性有限元程序的正確性。一圓柱薄殼如圖1所示，兩邊直線段鉸支，兩邊圓弧段自由。中點(diǎn)A作用集中載荷P=3 000 N。圓柱段半徑R=2 540 mm，弧長(zhǎng)對(duì)應(yīng)的圓心角θ=0.2 rad，直線段長(zhǎng)度L=508 mm。材料彈性模量E=310 275 N/mm2，厚度h=6.35 mm，泊松比ν=0.3。用基于弧長(zhǎng)法的共旋有限元理論求解，所得點(diǎn)A處的z向載荷因子-位移曲線如圖2所示。由圖2可知：計(jì)算結(jié)果與文獻(xiàn)[3]吻合很好，說明計(jì)算程序正確。

2.2翼夾角的影響

為考察連接翼張開夾角對(duì)結(jié)構(gòu)屈曲特性的影響，設(shè)計(jì)了如圖3所示的連接翼結(jié)構(gòu)。圖3中：hjoint為連接部分的高度；θ為上下翼間的夾角，計(jì)算中取θ=20°，30°，40°；連接翼的展向長(zhǎng)度b和根弦長(zhǎng)度c隨θ而變。設(shè)翼面需承載總載荷25 N，以均勻分布的定常壓強(qiáng)載荷作用于翼面，考察在此載荷作用下點(diǎn)P1、P2的位移，所得不同θ下點(diǎn)P1、P2的載荷因子-位移曲線分別如圖4、5所示，θ=20°時(shí)連接翼變形形狀如圖6所示。

由圖4可知：θ=20°時(shí)位移-載荷曲線存在明顯的非線性屈曲失穩(wěn)現(xiàn)象，結(jié)構(gòu)響應(yīng)在臨界載荷之前呈現(xiàn)剛度軟化現(xiàn)象，當(dāng)載荷因子到達(dá)臨界載荷后由于發(fā)生屈曲，結(jié)構(gòu)不能以原來的外形承受更多載荷繼而發(fā)生卸載過程，位移發(fā)生突跳，之后進(jìn)入結(jié)構(gòu)剛化階段；與θ=20°相比，θ=30°時(shí)結(jié)構(gòu)的位移-載荷因子曲線的載荷卸載過程很弱，不平衡路徑不明顯；θ=40°時(shí)在整個(gè)響應(yīng)過程中只存在結(jié)構(gòu)軟化和剛化兩個(gè)階段，無(wú)卸載過程，未發(fā)生位移突跳型屈曲。θ=20°，30°時(shí)的臨界屈曲載荷因子λ分別為0.257，0.205。

由圖5可知：在臨界載荷時(shí)點(diǎn)P2發(fā)生垮塌，結(jié)構(gòu)以新的形式繼續(xù)承受載荷。

由圖6可知：θ=20°時(shí)點(diǎn)P2由于屈曲結(jié)構(gòu)無(wú)法承受載荷而發(fā)生垮塌。

2.3連接翼高度的影響

不同類型飛行器的上下翼在空間上的連接高度hjoint往往相差較大。對(duì)大型高空長(zhǎng)航時(shí)飛行器來說，連接部分的高度相對(duì)較高，誘導(dǎo)阻力較小，但對(duì)“鉆石背”類型的連接翼結(jié)構(gòu)，因需考慮結(jié)構(gòu)的可折疊性，上下翼間的連接高度一般很小，甚至上下翼就在一個(gè)平面內(nèi)，故有必要分析hjoint對(duì)連接翼屈曲特性的影響。hjoint分別為10，20，30 mm時(shí)點(diǎn)P1載荷因子-位移曲線如圖7所示。由圖7可知：結(jié)構(gòu)的屈曲載荷因子隨連接高度的增加而逐漸增大，位移-載荷曲線模式呈現(xiàn)出明顯的軟化-不平衡路徑-剛化的位移突跳屈曲模式。

2.4連接位置的影響

考慮上下翼連接位置對(duì)連接翼屈曲特性的影響，連接翼結(jié)構(gòu)如圖8所示。設(shè)整個(gè)連接翼承受的總載荷25 N，上下翼分別承受載荷12.5 N，計(jì)算此時(shí)結(jié)構(gòu)在不同載荷-位移屈曲路徑的變化規(guī)律，結(jié)果如圖9所示。由圖9可知：隨著連接位置逐漸向根弦處靠近，連接翼的整體剛度逐漸加強(qiáng)，其位移-載荷因子曲線中沒有非平衡路徑。點(diǎn)P2的載荷因子-位移曲線如圖10所示。L=375 mm時(shí)不同載荷因子下的連接翼形狀如圖11所示。由計(jì)算結(jié)果可知：三種連接位置均未發(fā)生位移突跳型屈曲，但圖9中的載荷因子-位移曲線曲率發(fā)生了變化，由凸函數(shù)變?yōu)榘己瘮?shù)，表明結(jié)構(gòu)剛度經(jīng)歷了軟化至剛化的過程。由圖11可知：結(jié)構(gòu)從軟化段進(jìn)入剛化段后上翼點(diǎn)P2附加發(fā)生垮塌，形狀發(fā)生改變后繼續(xù)承受載荷。

2.5上下翼載荷分配的影響

連接翼結(jié)構(gòu)的氣動(dòng)優(yōu)點(diǎn)是可在不增加誘導(dǎo)阻力情況下改變機(jī)翼的載荷分布，因此有必要研究不同載荷分布條件下連接翼的屈曲行為。在機(jī)翼總載荷25 N條件下，考慮5種分布情況進(jìn)行分析，分別為：Fsy=25 N，F(xiàn)xy=0 N；Fsy=17.5 N，F(xiàn)xy=7.5 N；Fsy=12.5 N，F(xiàn)xy=12.5 N；Fsy=7.5 N，F(xiàn)xy=17.5 N；Fsy=0 N，F(xiàn)xy=25 N。此處：Fsy，F(xiàn)xy分別為作用于上翼和下翼的總載荷。所得點(diǎn)P1的位移-載荷因子曲線如圖12所示。由圖12可知：結(jié)構(gòu)的非線性屈曲臨界因子隨上翼載荷增大而逐漸增加，連接翼上翼載荷越大，整個(gè)結(jié)構(gòu)的剛度也越大。

3 結(jié)束語(yǔ)

本文以一連接翼為研究對(duì)象，采用三角形平板殼單元對(duì)其進(jìn)行空間離散，用共旋有限元方法建立非線性靜力學(xué)方程，通過弧長(zhǎng)法求解其在不同載荷作用下的非線性屈曲行為。在不同幾何參數(shù)條件下獲得的連接翼非線性屈曲的一般規(guī)律如下。第一，連接翼夾角的影響是：增加連接翼的夾角可降低位移突跳型屈曲的臨界載荷因子，但當(dāng)角度超過一定程度后連接翼反而不會(huì)發(fā)生位移突跳現(xiàn)象。說明發(fā)生位移突跳型屈曲存在一定范圍，只要幾何參數(shù)處于該范圍之外就不會(huì)發(fā)生此現(xiàn)象。在結(jié)構(gòu)軟化階段，連接翼角度愈大剛度越??；在剛化階段其規(guī)律則相反，角度越大剛度越大。第二，連接翼高度的影響是：增加連接翼高度可提高位移突跳型屈曲的臨界載荷因子；降低連接翼高度會(huì)增大連接翼的誘導(dǎo)阻力，同時(shí)也更易發(fā)生非線性屈曲失穩(wěn)。第三，連接翼連接位置的影響是：減小下翼和上翼連接處到根弦的距離可提高連接翼的剛度，減小翼稍的位移變形。這是因?yàn)橄乱淼拈L(zhǎng)度減短后其抗彎剛度得到提高，從而減小了對(duì)上翼的展向壓縮，提高屈曲失穩(wěn)的載荷因子。第四，上下翼載荷分配的影響是：由計(jì)算結(jié)果可知，對(duì)上翼分配較大載荷可提高屈曲臨界載荷因子，這是因?yàn)橄乱沓惺芨噍d荷時(shí)對(duì)上翼沿展向的壓縮就越大，從而促使上翼發(fā)生屈曲。上翼受載荷越大，其展向壓縮越小因此發(fā)生屈曲的臨界載荷也越大。隨著載荷因子的提高，連接翼結(jié)構(gòu)一般都會(huì)經(jīng)歷軟化-不平衡路徑-剛化三個(gè)階段，若無(wú)不平衡路徑階段則認(rèn)為結(jié)構(gòu)未發(fā)生位移突跳型屈曲。隨著載荷增加，上翼受壓后一般都會(huì)在軟化階段維持一個(gè)形狀，繼續(xù)施加載荷，上翼翼根處受壓垮塌后變?yōu)榱硪粋€(gè)形狀則進(jìn)入剛化階段。

[1] OLIVIERO F, FREDIANI A, RIZZO E. Design of an airfreight system based on innovative PrandtlPlane aircraft[C]// AIAA SciTech, 56thAIAA/ASCE/AHS/ASC Structures, Structural Dynamics and Materials Conference. Kissimmee: AIAA, 2015: 1-23.

[2] KROO I, SMITH S. Aerodynamic and structural studies for joined-wing aircraft[J]. Journal of Aircraft, 1991, 28(1): 74-81.

[3] HAJELA P, CHEN J L. Optimum structural sizing of conventional cantilever and joined wing configurations using equivalent beam models[R]. AIAA, 86-2653, 1986.

[4] SELBERG B P, CRONIN D L. Aerodynamic-structural optimization of positive/negative stagger joined wing configurations[R]. AIAA, 86-2626, 1986.

[5] SMITH S C, CLIFF S E. The design of a joined-wing flight demonstrator aircraft[R]. AIAA, 87-2930, 1987.

[6] STEVEN H C, CHRISTOPHER M B. Planform effects on high speed civil transport design[R]. AIAA, 88-4487, 1988.

[7] GALLMAN J W, KROO I M, SMITH S C. Design synthesis and optimization of joined-wing transports[R]. AIAA, 90-3197, 1990.

[8] LIVNE E. Aeroelasticity of joined-wing airplane configurations-past work and future challenges——a survey[R]. AIAA, 2001-1370, 2001.

[9] DEMASI L, CAVALLARO R, RAZóN A M. Postcritical analysis of Prandtl plane joined-wing configurations[J]. AIAA Journal, 2013, 51(1): 161-177.

[10] DEMASI L, CAVALLARO R, BERTUCCELLI F. Post-critical analysis of joined wings: the concept of snap-divergence as a characterization of the instability[C]// 5thAIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics and Materials Conference. Boston: AIAA, 2013: 1-32.

[11] 潘家正, 呂慶風(fēng), 周欲曉. 連接翼布局直接力控制可能性的初步探索[J]. 飛行力學(xué), 1996, 14(1): 48-53.

[12] 李光里, 李國(guó)文, 黎軍, 等. 連接翼布局氣動(dòng)特性研究[J]. 空氣動(dòng)力學(xué)學(xué)報(bào), 2006, 24(4): 513-519.

[13] 劉杰, 張曙光. 連接翼布局縱向控制特性[J]. 航空學(xué)報(bào), 2008, 29(SI): 91-96.

[14] 冉玉國(guó), 李秋彥, 金偉, 等. 大展弦比連接翼飛機(jī)考慮預(yù)載作用下的非線性模態(tài)與顫振特性研究[C]// 第十一屆全國(guó)空氣彈性學(xué)術(shù)會(huì)議, CARS-2009. 昆明: 中國(guó)空氣動(dòng)力學(xué)會(huì)空氣動(dòng)力彈性專業(yè)委員會(huì), 2009: 28-33.

[15] 張書俊, 王運(yùn)濤, 孟德虹. 大展弦比連接翼靜氣動(dòng)彈性研究[J]. 空氣動(dòng)力學(xué)學(xué)報(bào), 2013, 31(2): 170-174.

[16] 張波成, 萬(wàn)志強(qiáng), 楊超. 連翼布局飛行器飛行載荷與顫振分析[J]. 工程力學(xué), 2010, 27(8): 229-240.

[17] 李少華, 楊智春, 谷迎松, 等. 復(fù)合材料連翼的氣動(dòng)彈性剪裁研究[J]. 西北工業(yè)大學(xué)學(xué)報(bào), 2008, 26(3): 292-296.

[18] 吳小勝, 雷娟棉, 胡俊. 鉆石背彈翼的靜氣動(dòng)彈性研究[J]. 北京理工大學(xué)學(xué)報(bào), 2010, 30(12): 1387-1390.

[19] 吳小勝, 吳甲生. “鉆石背”彈翼顫振的數(shù)值計(jì)算[J]. 北京理工大學(xué)學(xué)報(bào), 2007, 27(5): 385-389.

[20] 雷娟棉, 吳甲生. 鉆石背彈翼外形參數(shù)對(duì)氣動(dòng)特性的影響[J]. 北京理工大學(xué)學(xué)報(bào), 2006, 26(11): 945-949.

[21] 雷娟棉, 吳甲生. 鉆石背彈翼氣動(dòng)特性風(fēng)洞實(shí)驗(yàn)研究[J]. 兵工學(xué)報(bào), 2007, 28(7): 893-896.

[22] 許兆慶, 吳軍基, 薛曉中, 等. 巡飛彈鉆石背折疊翼結(jié)構(gòu)優(yōu)化設(shè)計(jì)[J]. 機(jī)械設(shè)計(jì)與研究, 2011, 27(1): 87-90.

[23] 黃益民, 葛森, 吳煒, 等. 鉆石背彈翼抗風(fēng)動(dòng)力可靠性及靈敏度分析[J]. 機(jī)械科學(xué)與技術(shù), 2013, 32(3): 416-420.

[24] 余學(xué)芳. 混凝土結(jié)構(gòu)裂縫的有限元模擬[D]. 杭州: 浙江大學(xué), 2012.

[25] CRISFIELD M. Non linear finite element analysis of solid and structures[M]. New York: Wiely, 1991: 1-345.

[26] 楊勁松, 夏品奇. 薄殼大轉(zhuǎn)動(dòng)、小應(yīng)變幾何非線性分析共旋有限元法[J]. 中國(guó)科學(xué): 科學(xué)技術(shù), 2012, 42(11): 1295-1304.

[27] FELIPPA C A. A study of optimal membrane triangles with drilling freedoms[J]. Comput Methods Appl Mech Engrg, 2003, 192: 2125-2168.

[28] BATOZ J L, BATHE K J, HO L W. A study of three-node triangular plate bending elements[J]. Int Numer Methods Eng, 1980, 15: 1771-1812.

[29] SZE K Y, LIU X H, LO S H. Popular benchmark problems for geometric nonlinear analysis of shells[J]. Finite Elements Anal Design, 2004, 40: 1551-1569.

AnalysisonNonlinearBucklingofJoined-WingConfigurations

DOU Yi-bin, ZHANG Xiao-hong

(Shanghai Electromechanical Engineering Institute, Shanghai 201109, China)

The nonlinear buckling behavior of a joined-wing configuration was analyzed by the simulation with the nonlinear finite element method, and the effect of varied geometry parameters as well as wing’s load repartition on nonlinear buckling was studied. The joined-wing was modeled with triangular planer shell element based on Co-Rotational framework. The snap-through phenomenon and post-critical pattern was investigated with the adoption of the arc-length technique. This study showed that the typical load-displacement curve of joined-wing was composed by three branches which were stiffness softening, snap-buckling unstable and stiffness stiffening. The buckling load reduced gradually when the angle between the upper and lower wings increased. And there was a critical angle existed that no snap-through phenomenon occurred for joined-wing with bigger angles than this critical angle. Increasing the height of joint was beneficial for increasing the snap-buckling load of the joined-wing. The lift repartition and its effects on buckling showed that an improvement of the structural stability of the joined-wing airplane would be achieved by increasing the amount of load on the upper wing and reducing the load on the lower wing. While keeping the angle between the upper and lower wings unchanging, the decrease of the span of lower wing improved the structural response by increasing the snap-buckling load. And the snap-through phenomenon was not occurred.

joined-wing; nonlinear buckling; Co-Rotational framework; snap-through; arc-length technique; stiffness softening; stiffness stiffening; critical included angle

1006-1630(2017)05-0040-06

2017-03-25；

2017-06-08

竇怡彬(1984—)，男，博士，主要研究方向?yàn)闅鈩?dòng)彈性力學(xué)。

V211.3

A

10.19328/j.cnki.1006-1630.2017.05.006