楊易鑫,印寅,聶宏,魏小輝
南京航空航天大學 機械結構力學及控制國家重點實驗室,南京 210016
作為起落架收放系統的重要組成部分,撐桿式鎖機構是保證飛機安全起降必不可少的機構之一[1],可將起落架牢固地鎖定于收放位置,承擔飛行時的前起落架重力并傳遞滑跑時的地面載荷[2]。接近上鎖時,由開鎖的不穩定狀態瞬間跳躍至鎖死的穩定狀態(解鎖過程相反),是鎖機構的顯著特征之一,如圖1所示。根據分岔理論[3-4],鎖機構穩定性改變的瞬間即出現了分岔。對鎖機構性能的研究自然集中于分岔出現的臨界點(分岔點)的研究上。工業上一般通過動力學仿真研究起落架收放機構性能[5-6]。在鎖機構設計方面,分析機構穩定性隨參數緩慢改變的變化趨勢時多采用多體動力學軟件,如LMS、ADAMS等,執行多時域仿真[7]。但該方法通常模型復雜,并且選定的參數每改變一次,動力學模型便需重新仿真一次,如此會耗費大量計算時間,尤其考慮到影響參數較多的情況,計算量會呈指數倍增長,設計成本高昂;同時僅僅依靠仿真結果難以解釋鎖機構在解鎖與上鎖瞬間的跳躍現象,也難以捕捉分岔點并分析不同參數對分岔點的影響情況。這在一定程度上限制了鎖機構設計水平的發展[8]。
分岔分析有助于解決上述問題。它基于分岔理論,采用數值延拓法[9],從已知的非線性常微分方程組平衡解出發,在一定誤差范圍內,延拓出所有平衡解隨選定參數變化的軌跡曲線從而找出分岔點。它的一個顯著優勢是,通過COCO(Computational Continuation Core)[10]等數值延拓程序包,能夠以模型任意某狀態變量為延拓參數,無需重構方程組,直接延拓出系統關鍵特性(如分岔點)隨不同參數的變化趨勢。這樣的靈活性,不僅能夠快速分析各種參數對鎖機構鎖定性能的影響并指導機構設計,而且有助于深入了解機構動力學特性,闡述非線性行為。在飛機起落架研發初期可成為強有力的分析手段。

圖1 鎖機構跳躍現象示意圖Fig.1 Schematic of jumping phenomenon of lock mechanism
分岔分析在航空領域已有廣泛運用[11],例如起落架擺振問題,Terkovics[12]、Thota[13]以及陳大偉[14]等便基于飛機前進速度和作用于起落架上的垂向力2個參數,采用分岔分析方法,考察了輪胎充氣壓力、橫向轉彎等因素對前起落架擺振穩定區域及其邊界范圍的影響。在飛機地面操縱問題上,Coetzee[15]和Krauskopf[16]等順利聯合AUTO[17]與SimMechanics模型,使用該方法分析了飛機滑行速度和轉彎角度對A320地面機動穩定性的影響。此外,在飛行動力學[18]尤其是軍機控制[19]領域,分岔分析也發揮了重要作用。
在鎖機構分岔分析方面,學者Knowles等[20-21]開展了一系列研究,以前起落架撐桿鎖機構為對象,先后考察了鎖彈簧各參數對下位上鎖性能的影響[22];分析了彈簧剛度對上位上鎖分岔點的的影響并給出了成功上鎖所需的臨界彈簧剛度值[23],但未考慮彈簧各參數對解鎖性能的影響。Yin等[24]針對三維單側撐桿主起落架,研究了下位解鎖力對解鎖分岔點的影響并優化了臨界解鎖力,但忽略了上位解鎖情況,并且改變彈簧參數可能會對上、下位解鎖產生不同,甚至相反的影響。
為彌補相關研究的不足,本文針對某前起落架撐桿式鎖機構,綜合考慮上、下位解鎖過程,運用分岔分析方法研究前起落架收放過程中,解鎖作動筒作用力對鎖機構性能(鎖撐桿運動軌跡分岔點)的影響,定義解鎖所需的最小作用力(臨界解鎖力)與臨界解鎖角度,并分析彈簧各參數對臨界解鎖力的影響;最終以下位解鎖為例,優化臨界解鎖力。
撐桿式鎖機構是將鎖桿支撐于阻力桿折疊處,上鎖后使阻力桿無法折疊以承受拉壓載荷的鎖機構。某前起落架(Noise Landing Gear,NLG)采用上下位鎖合一的支承式撐桿鎖機構,即上下位一體式鎖機構。其鎖桿與機體相連,剛度較好,上鎖后不會因受載變形而開鎖[24]。結構如圖2所示,該機構由上、下鎖撐桿,鎖彈簧與解鎖作動筒組成,鎖撐桿間裝有止動裝置(擋板和擋塊)。其工作原理為:上鎖時,在彈簧力與鎖桿自身重力共同作用下,上、下鎖撐桿向下折疊越過中心位置(兩鎖桿夾角達180°)后,擋板與擋塊相互接觸阻礙桿件繼續運動,完成上鎖,將前起落架固定于放下或收上位置;解鎖時,解鎖作動筒施加載荷克服彈簧力,拉動(下位解鎖)或推動(上位解鎖)上、下鎖撐桿向上折疊越過中心位置完成開鎖。

圖2 前起落架及上下位一體式鎖機構結構示意圖Fig.2 Schematic of nose landing gear and combined uplock/downlock mechanism
考慮到前起落架真實機構零件眾多,結構連接復雜,為便于模擬計算同時不使模型過于理想偏離實際,對前起落架機構作出如下假設:① 將緩沖器外筒,活塞桿,輪軸,機輪以及扭力臂等簡化為一個整體構件,統稱為主支柱。主支柱作為前起落架收放的主要旋轉部件,其重心通過計算模擬進行等效;② 各部件間為理想約束;③ 忽略各運動副連接的間隙與摩擦;④ 收放作動筒、解鎖作動筒和鎖彈簧等機構運動對靜力學模型無影響,視為力元而不作為運動部件處理。
適當簡化前起落架機構后,結構如圖3所示。機構主要由主支柱(L1)、上阻力桿(L2)、下阻力桿(L3)、下鎖撐桿(L4)、上鎖撐桿(L5)組成,桿件間通過平面旋轉副連接(連接點由B、C、D表示)。其中,上阻力桿、上鎖撐桿、主支柱以及彈簧與機體連接點(分別對應A、E、O、O1點)假定為固定于基座上的點。解鎖作動筒作用力(Ful)作用于上鎖撐桿(H點);收放作動筒作用力(Fa)作用于主支柱(K點)。建立靜力學模型時所參照的全局坐標系,以O點為原點;x軸方向水平向右逆于航向;y軸垂直于x軸方向向上。
前起落架第i根桿件主要包含4個元素:
Li=(xi,yi,θi,li)
(1)
式中:各參數如圖4所示。(xi,yi)表示桿件重心Gi的全局坐標;θi表示桿件軸線(實線所示,如AB)與x軸的夾角;li表示桿件長度。需要注意的是,各桿件并非對稱結構,重心G與軸線存在偏差,重心與端點間的連線(虛線所示,如AG2、BG2,其長度分別為l21、l22)與軸線間存在夾角β(如β21、β22);主支柱較為特殊,l11~l13分別為OG1、CG1、OC連線長度,β12為OG1和CG1夾角,β13為OG1和OC夾角。考慮到桿件長度為固定

圖3 前起落架收放狀態簡圖Fig.3 Schematic of NLG deployed and retracted

圖4 前起落架各桿件結構參數簡圖Fig.4 Structural parameter diagrams of NLG bars
值,如此每根桿件的狀態變量為:Xi=[xi,yi,θi]T,機構共包含15個狀態變量,根據結構連接關系可知運動自由度為1(以主支柱收放角度θstrut為驅動自由度),因此需14個約束方程描述機構的幾何約束關系,如式(2)所示。式中(Ax,Ay)、(Ex,Ey)分別為A、E點坐標。
(2)

圖5 前起落架各桿件受力簡圖Fig.5 Force diagram of NLG


解鎖作動筒作用力及其力矩(以E點為計算中心):
(3)
收放作動筒作用力及其力矩(以O點為計算中心):
(4)
鎖彈簧力及其力矩(以E點為計算中心):
(5)
各參數可參考圖6。式(3)中,(xu1,yu1)、(xu2,yu2)分別為解鎖作動筒兩端點O2、H坐標,β54為HE、DE間夾角;式(4)中,(xa1,ya1)、(xa2,ya2)分別為收放作動筒兩端點O3、K坐標,l14為OK長度,β14為OK、OG1間夾角;式(5)中,(xs1,ys1)、(xs2,ys2)分別為彈簧兩端點O1、F坐標,β53為FE、DE間夾角;k為鎖彈簧剛度;lu為鎖彈簧原長。

圖6 彈簧力、解鎖力、收放力計算參數示意圖Fig.6 Parameter diagram for calculating spring force, unlock actuator force and retraction actuator force
上阻力桿、上鎖撐桿、主支柱的力矩平衡方程分別以A、E、O點為計算中心;下阻力桿和下鎖撐桿的力矩平衡方程以B點為計算中心。最終力矩平衡方程為
(6)
力平衡方程包括下阻力桿和下鎖撐桿在x、y方向的力平衡方程:
(7)
整合所有力和力矩平衡方程,可寫成矩陣形式:
AF-B=0
(8)
式中:A為力系數矩陣;F為狀態變量矢量;B為剩余項矢量。具體見附錄A。
為使前起落架順利、平穩地完成收、放動作,需要解鎖作動筒和收放作動筒按照一定順序驅動起落架運動。以收起過程為例,放下過程與之相似。首先由解鎖作動筒拉動鎖撐桿實現下位解鎖,隨后在收放作動筒作用下收起前起落架,至接近收上位置時依靠彈簧拉力與鎖桿重力實現上位上鎖。通常解鎖作動筒會在起落架成功解鎖后泄壓。但需要說明的是,文獻[23]指出,解鎖作動筒必須在收放作動筒載荷足夠支撐起落架重量后泄壓。因此,在解鎖作動筒作用力不干擾上位上鎖的前提下,為保證不影響起落架正常收起,起落架收上后解鎖作動筒再泄壓。整個收放過程中作動筒工作流程如圖7所示,前起落架在不同作用力下的位置狀態由點C、RU、RD、EU、ED表示。
不難發現,解鎖作動筒作用力Ful和收放作動筒作用力Fa是起落架收放過程中的關鍵因素。本節將基于數值延拓算法對前起落架鎖機構進行分岔分析,于MATLAB中求解靜力學模型(式(2)、式(8))平衡解并進行數值延拓,得到Fa和Ful對鎖機構運動軌跡及其穩定性的影響。同時采用LMS動力學軟件對鎖機構運動軌跡進行傳統動力學仿真,并將仿真軌跡與延拓軌跡進行比較以驗證結果準確性。各結構參數值見附錄B。為便于區分鎖撐桿向上與向下折疊運動,采用過中心角度θov描述鎖機構狀態。該參數表示上、下鎖撐桿角度之差,即θov=θ4-θ5,定義上下鎖撐桿伸展成一條直線狀態時θov=0°。θov<0°時鎖撐桿向下折疊,處于下過中心狀態,θov>0°時鎖撐桿向上折疊,處于上過中心狀態。對于單自由度的前起落架機構,在鎖撐桿折疊方向確定的情況下,主支柱收放角度θstrut和鎖撐桿過中心角度θov一一對應。

圖7 前起落架作動筒工作流程圖Fig.7 Actuator force schedules of NLG
收放過程鎖機構運動軌跡如圖8所示。其中粗灰色曲線為動力學仿真結果,細線為數值延拓結果:紅色實線為穩定解分支,表示前起落架在收放作動筒作用下的響應(包含了起落架收、放過程),藍色虛線表示不穩定解分支;黑色實線表示前起落架在解鎖作動筒作用下的響應(即解鎖或解鎖作動筒泄壓過程)。分岔點由*表示,標注為SN1~SN4。為清晰地顯示主支柱收放運動與鎖機構折疊運動之間的對應關系,得到θov-θstrut平面內的前起落架運動軌跡如圖9所示,θstrut=0°時前起落架處于收上位置;θstrut=101.2° 時前起落架處于放下位置。
根據圖8、圖9中箭頭所示順序可闡述前起落架收上過程。作動筒按圖7所示流程工作,解鎖階段,Fa=0 kN。起落架僅在解鎖作動筒拉力作用下開鎖,Ful增加至0.43 kN時,鎖撐桿由初始點C1(θov= -20°)運動至解鎖位置RU1(θov=20°),圖9表明在這一過程中主支柱角度幾乎沒有變化;解鎖后Ful保持不變,逐漸增加Fa,θstrut不斷減小,起落架持續收起,上、下鎖撐桿由位置RU1經RU2(θov=80°)、RU3(θov=150°)、RU4(θov=70°)緩慢運動至分岔點SN2處。此時臨近收上位置。可以看出,上述運動過程動力學仿真結果與數值延拓結果吻合較好,而分岔點SN2后2種結果出現偏離,分析數值延拓曲線后不難發現,由于分岔點位于實線和虛線交接處,說明該點兩側平衡解將由穩定狀態變為不穩定狀態。而對于不穩定狀態,任何微小擾動都會使其偏離發散到穩定狀態,因此在準靜態動力學仿真過程中,機構實際不會遵循不穩定解分支運動。SN2點后,鎖撐桿動力學仿真運動軌跡將由上過中心穩定曲線跳躍至下過中心穩定曲線,過中心角度發生突變,到達位置RD1(θov=-40°)。

圖8 前起落架收放循環分岔圖Fig.8 Bifurcation diagram of NLG extension/retraction cycle

圖9 前起落架收上過程鎖機構過中心角度與主支柱角度對應關系分岔圖Fig.9 Bifurcation diagram of NLG retraction cycle with over-centre angle as function of strut angle
需要特別指出的是,為全面考察鎖機構在各過中心角度下的運動軌跡,延拓曲線中鎖桿可越過擋板位置繼續折疊。實際情況中,擋板和擋塊在θov= -2°(該角度也稱為上鎖角度)時互相接觸,阻止鎖撐桿繼續運動,如圖8、圖9中黑色虛線所示。因此,上文所述跳躍至位置RD1的情況在真實收上過程中不會出現,起落架會在位置RD1*處完成上位上鎖,分岔點SN2因此可稱為上位上鎖分岔點。前起落架鎖于收上位置(圖3(b)),收上過程完成。此后解鎖作動筒泄壓(無擋板作用時前起落架運動至RD2位置后泄壓)。收上過程中鎖撐桿不同狀態,都能在圖10中找到相應的機構簡圖。
注意到當允許鎖桿向下折疊至位置RD1后,若逐漸降低Fa,前起落架沿下過中心曲線由RD3(θov= -150°)至RD4(θov= -80°)緩慢放下時,由于存在Ful,出現了分岔點SN3。到達該點后鎖撐桿同樣發生了跳躍現象,類似于C1-RU1的解鎖過程。SN3可稱為下位解鎖分岔點。
前起落架放下過程與收上過程相似,如圖11箭頭所示:上位解鎖時,需要收放作動筒施加載荷克服起落架重力后,由解鎖作動筒完成解鎖。因此在Fa=15 kN時,解鎖作動筒反向施加推力(Ful= -0.3 kN),鎖撐桿由C2(θov= -15°)運動至解鎖位置EU1(θov=8°)完成上位解鎖。此后隨Fa減小,前起落架逐漸放下至EU4(θov=80°),鎖撐桿折疊運動如圖12所示。
動力學仿真軌跡與延拓軌跡同樣有較高的重合度,而當到達分岔點SN1后,動力學軌跡將再次發生上文所述跳躍,鎖撐桿由向上折疊狀態變為向下折疊狀態。同樣地,由于擋板的存在,鎖撐桿實際不會折疊至ED1(θov= -70°),而于ED1*處完成下位上鎖。類似地,稱分岔點SN1為下位上鎖分岔點。放下過程結束后解鎖作動筒泄壓(允許鎖桿自由折疊時,前起落架運動至ED2位置后泄壓)。

圖10 前起落架收上過程簡圖Fig.10 Sketch of NLG retraction cycle

圖11 前起落架放下過程鎖機構過中心角度與主支柱收放角度對應關系分岔圖Fig.11 Bifurcation diagram of NLG extension cycle with over-centre angle as function of strut angle

圖12 前起落架放下過程簡圖Fig.12 Sketch of NLG extension cycle
不難發現當鎖桿能夠越過擋板繼續向下折疊至位置ED1,并逐漸增加Fa收上前起落架時,在Ful作用下出現了分岔點SN4。鎖撐桿由ED3、ED4位置到達該點后,發生了相似于C2-EU1解鎖過程的跳躍。SN4可稱為上位解鎖分岔點。
通過以上分析可以看出,動力學仿真結果與數值延拓結果在穩定解分支上基本吻合,而由于分岔點后出現不穩定解分支,導致動力學軌跡發生跳躍進而實現起落架上鎖。
由3.1節分析可知,平衡解,尤其是分岔點(SN1~SN4)的出現決定了機構的動力學響應。前起落架鎖機構某些關鍵行為,例如上位上鎖、下位上鎖等,常常與分岔點的出現相關。

圖13 不同解鎖力下鎖撐桿運動軌跡分岔圖Fig.13 Bifurcation diagrams of lock links with different unlock actuator forces
對比圖13不同解鎖作動筒作用力下的鎖撐桿運動軌跡可以發現,Ful變化會導致分岔點的位置發生改變,進而對前起落架機構的收放運動產生影響。本小節將詳細分析不同Ful對分岔點的影響,找到前起落架能夠順利上、下位解鎖的臨界解鎖力。需要說明的是,由于動力學仿真難以準確捕捉分岔點更無法分析其隨不同參數的變化趨勢,而數值延拓方法的準確性已由3.1節得到驗證。因此以下內容均采用數值延拓進行分析。
(1) 上位解鎖
分別選取Ful= -0.2 kN、Ful=0 kN 2種工況進行分析,延拓結果如圖14所示,各線型與標示含義與圖8相同。由放下過程分析可知,為便于上位解鎖,收放作動筒需施加作用力以克服重力。Ful=0 kN的工況下,Fa=15 kN時,只能在圖14(b)的下過中心曲線上發現一個穩定平衡解。這意味著若此時降低Fa放下前起落架,鎖撐桿會向下折疊沿下過中心曲線運動。但如前文所述,擋板阻擋了鎖桿的運動趨勢,該工況實際未完成上位解鎖。為實現解鎖順利放下起落架,需提高Ful使平衡解移動至上過中心穩定曲線。施加推力至0.2 kN(圖14(a))后,Fa=15 kN時對應的平衡解處于上過中心曲線,說明鎖撐桿完成上位解鎖,可以向上折疊使前起落架放下。

圖14 上位解鎖時不同解鎖力下鎖撐桿運動軌跡分岔圖Fig.14 Bifurcation diagrams of lock links with different unlock actuator forces in case of up unlocking

圖15 分岔點SN2、SN4隨解鎖力變化的軌跡線Fig.15 Locus of bifurcation points SN2SN4 varying with unlock actuator force
因此,上位解鎖時可根據不同Ful值將鎖撐桿解鎖狀態分為2類:上位未解鎖狀態和上位解鎖狀態。對比圖14(a)和(b)不難發現,不同解鎖狀態與分岔點SN2、SN4的出現與否緊密相關。對點SN2、SN4進行數值延拓可更直觀地了解分岔點隨Ful的變化趨勢,如圖15中黑色曲線所示。將軌跡線分別投影至Fa-Ful、θov-Ful平面(如圖16、圖17所示)以便清晰地觀察曲線走勢,其中藍色曲線表示上位解鎖曲線。


圖16 分岔點軌跡線與上位解鎖曲線(θov-Ful)Fig.16 Locus of bifurcation points and up unlocking curve (θov-Ful)

圖17 分岔點軌跡線與上位解鎖曲線(Fa-Ful)Fig.17 Locus of bifurcation points and up unlocking curve (Fa-Ful)
(2) 下位解鎖
選取Ful=0.1 kN、0.2 kN、0.3 kN、0.4 kN這4種工況分別進行延拓計算,結果如圖18所示。
下位解鎖時,初始狀態沒有收放作動筒作用。當Ful=0 kN,Fa=0 kN時,如圖14(b),此時與上位未解鎖狀態相似,唯一對應的穩定平衡解位于下過中心曲線,實際前起落架未解鎖,無法收起。使初始狀態平衡解移動至上過中心穩定曲線依然是成功解鎖的關鍵。Ful增加至0.1 kN、0.2 kN時,由圖18(a)和圖18(b)可知,均未成功解鎖;進一步增長至0.3 kN后,如圖18(c)所示,上過中心穩定曲線與縱軸相交,意味著起落架能夠解鎖,鎖撐桿隨Fa增加可向上折疊,但折疊幅度(即θov增加幅度)遠不及Fa的增長程度,前起落架實際沒有明顯收上。并且隨后出現了2個分岔點SN1、SN1*將軌跡線分為上、下穩定分支,形如遲滯回線(圖18(a)中下過中心曲線同樣形成了該回線)。Fa需超過SN1*點對應的力值(約5 kN),平衡解從下穩定分支躍至上穩定分支后,鎖撐桿才能大幅折疊,實現前起落架穩定收起。因此,本工況中雖成功解鎖,但遲滯回線的存在證明解鎖仍然不夠徹底。圖18(d)中遲滯回線消失,表明Ful=0.4 kN時前起落架能徹底解鎖并順利收上。
綜上所述,類似地,可將不同Ful下的鎖撐桿下位解鎖狀態分為3類:下位完全解鎖狀態(圖18(d))、下位不完全解鎖狀態(圖18(c),鎖撐桿向上折疊,Fa需增加至一定數值越過分岔點后起落架方能收上)、下位未解鎖狀態(圖18(a)和圖(b))。

與上位解鎖類似,不同解鎖狀態與SN1、SN3點息息相關。同樣對分岔點進行數值延拓,得到其隨不同Ful的變化軌跡,如圖19所示。將曲線分別投影至Fa-Ful、θov-Ful平面(圖20、圖21)可更為清晰地觀察其走勢。
結合圖19、圖20可知,Ful<0.22 kN時,鎖撐桿處于未解鎖狀態,SN1點始終存在,SN3、SN3*點于尖點C3后出現,對應于下過中心遲滯回線的形成,Ful>0.15 kN時SN3*點消失;隨Ful增至0.22 kN后,SN1*點逐漸出現,即形成上過中心遲滯回線,形成過程對應未解鎖狀態向不完全解鎖狀態轉化的過程;隨后2個分岔點再慢慢靠近直至交匯于尖點C2處后最終消失。SN1、SN1*的消失表明鎖撐桿由不完全解鎖狀態轉化為完全解鎖狀態,C2點為轉換的臨界點。注意到圖21中2條軌跡線在θov-Ful平面內的投影出現了一個缺口,圖18(b)代表了該缺口所對應的前起落架延拓曲線,此時鎖撐桿隨Fa增加不發生折疊運動。出現缺口的原因是當Ful逐漸減小至接近未解鎖區域時,Fa會漸漸趨近于無窮大(超出了圖20所示的力值上限)。換言之,在未解鎖區域和不完全解鎖區域的交界處,從上過中心遲滯回線下穩定解分支跳躍至上穩定解分支所需的收放作動筒作用力接近無窮大,遠遠超出了機載液壓系統的能力范圍,起落架因此由不完全解鎖區域過渡至未解鎖區域,無法正常收上。

圖19 分岔點SN1、SN3隨解鎖力變化的軌跡線Fig.19 Loci of bifurcation points SN1SN3 varying with unlock actuator force

圖20 分岔點軌跡線與下位解鎖曲線(Fa-Ful)Fig.20 Loci of bifurcation points and down unlocking curve (Fa-Ful)


圖21 分岔點軌跡線與下位解鎖曲線(θov-Ful)Fig.21 Loci of bifurcation points and down unlocking curve (θov-Ful)

選取不同彈簧剛度、原長,延拓出不同工況下的SN2、SN4點以及SN1、SN1*點軌跡線,并分別提取曲線臨界點C1、C2得到其隨彈簧剛度、原長的變化趨勢,如圖22、圖23所示。可以發現,隨著彈簧剛度增加、原長縮短,上位和下位臨界解鎖力均呈線性增長,意味著解鎖難度增加了,臨界解鎖角度有所上浮進一步證實了這一趨勢。分析原因,增強彈簧剛度或降低彈簧原長都增加了鎖彈簧作用于上鎖撐桿的彈簧力,不利于解鎖的力矩增大,致使解鎖難度加大。

圖22 不同彈簧剛度對臨界解鎖力和臨界解鎖角度 的影響Fig.22 Effects of different spring stiffness on critical unlocking force/angle

圖23 不同彈簧原長對臨界解鎖力和臨界解鎖角度 的影響Fig.23 Effects of different spring unstrained lengths on critical unlocking force/angle
鎖彈簧上安裝點與機身相連,默認為不可調節的固定安裝點,因此僅可改變上鎖撐桿處彈簧的安裝位置。可通過2種方式調節:改變安裝距離sm(EF長度)或安裝角度β(EF與DE夾角),如圖24所示。選取不同sm、β,分別得到臨界點變化趨勢,如圖25、圖26所示。由圖可知,在選定的變化范圍內,增大安裝角度時,減輕了上位解鎖的難度而下位解鎖與之相反;彈簧安裝距離增加時,上、下位解鎖均越來越困難。但相對而言,延長安裝位置更顯著地提高了解鎖難度。


圖24 彈簧不同安裝角度和安裝距離示意圖Fig.24 Schematic of different spring installation angles/distances

圖25 不同彈簧安裝角度對臨界解鎖力、解鎖角度的影響Fig.25 Effects of different spring installation angles on critical unlocking force/angle





圖27 尋優算法結合數值延拓算法優化流程圖Fig.27 Optimization flowchart based on combination of MIGA and COCO


表1 下位臨界解鎖力優化結果

圖28 下位臨界解鎖力優化結果分岔圖Fig.28 Bifurcation diagram of optimization results for critical down unlocking force
本文提出了基于分岔理論,通過數值延拓分析撐桿式鎖機構穩定性并指導參數設計的方法。對前起落架上下位一體式鎖機構進行了分岔分析和優化設計,得到了以下重要結論:
1) 采用分岔分析方法研究鎖機構具有可行性,該方法能快速把握模型穩定性隨選定參數的變化趨勢,準確預測系統動力學行為,提供高質量的設計參數。
2) 根據不同解鎖作動筒作用力,發現了鎖撐桿2種上位解鎖狀態和3種下位解鎖狀態,確定了解鎖所需的最小解鎖力,即臨界解鎖力及解鎖曲線上相應的臨界解鎖角度。
3) 彈簧各結構參數對臨界解鎖角度影響較小,但會顯著影響臨界解鎖力的大小,具體表現為:當增加彈簧剛度、安裝距離,減小彈簧原長時,上、下位解鎖難度均會提高;增加彈簧安裝角度雖然可緩解上位解鎖難度,但同時會使下位解鎖更困難。
4) 以下位解鎖為例,采用多島遺傳全局尋優算法,在下位上鎖能力一定的條件下,優化彈簧參數后,臨界解鎖力下降了34.7%。