雷萍 劉冬平

摘要：為了研究貨艙門收放機(jī)構(gòu)的載荷及運(yùn)動特性，根據(jù)力矩平衡原理、伯努利方程及轉(zhuǎn)動定理推導(dǎo)出該收放機(jī)構(gòu)的力學(xué)和運(yùn)動學(xué)模型;利用Simulink建立了貨艙門收放機(jī)構(gòu)的計(jì)算模型，通過仿真得到了貨艙門收放機(jī)構(gòu)載荷、速度、加速度變化等重要參數(shù)。仿真結(jié)果表明，該模型有效地模擬了收放機(jī)構(gòu)載荷及運(yùn)動特性，為貨艙門收放機(jī)構(gòu)的設(shè)計(jì)提供了參考。

關(guān)鍵詞：Simulink;貨艙門;收放機(jī)構(gòu);建模與仿真

中圖分類號：V19 文獻(xiàn)標(biāo)識碼：A

運(yùn)輸類飛機(jī)貨艙門是機(jī)身中重要且特殊的運(yùn)動部件，又是設(shè)計(jì)復(fù)雜、功能要求嚴(yán)格、有一定代表性的部件[1]，收放機(jī)構(gòu)是貨艙門上典型的動力支持機(jī)構(gòu)，合理地確定艙門收放機(jī)構(gòu)的載荷及運(yùn)動特性，對艙門收放機(jī)構(gòu)的設(shè)計(jì)及強(qiáng)度校核具有重要意義，并最終影響設(shè)計(jì)的合理性和結(jié)構(gòu)的安全性。

目前，對貨艙門收放機(jī)構(gòu)的設(shè)計(jì)傳統(tǒng)方法是結(jié)構(gòu)力學(xué)、運(yùn)動學(xué)等理論進(jìn)行計(jì)算，通過試驗(yàn)手段進(jìn)行測試，然后進(jìn)行性能改進(jìn)，其過程繁瑣、設(shè)計(jì)周期長，而利用Simulink工具搭建貨艙門收放機(jī)構(gòu)的仿真模型[2]，對貨艙門收放機(jī)構(gòu)的載荷及運(yùn)動特性進(jìn)行仿真計(jì)算，可得到貨艙門收放機(jī)構(gòu)載荷、速度、加速度等重要參數(shù)的變化規(guī)律，為收放機(jī)構(gòu)的載荷及運(yùn)動分析提供一種快速有效的方法。該模型構(gòu)建的貨艙門收放機(jī)構(gòu)參數(shù)化設(shè)計(jì)工具[3]，便于不同工況下收放機(jī)構(gòu)的設(shè)計(jì)，具有一定的工程意義。

1 貨艙門的運(yùn)動分析

1.1 貨艙門運(yùn)動機(jī)構(gòu)的受力分析

貨艙門在開啟過程中靠電機(jī)驅(qū)動貨艙門收放機(jī)構(gòu)的運(yùn)動控制貨艙門的運(yùn)動軌跡，貨艙門受力分析如圖1所示。

圖1中O點(diǎn)為貨艙門相對于門框的轉(zhuǎn)軸，A點(diǎn)為貨艙門收放機(jī)構(gòu)與機(jī)身鉸接點(diǎn)，C點(diǎn)為貨艙門收放機(jī)構(gòu)與貨艙門的連接點(diǎn)，也是貨艙門收放機(jī)構(gòu)的驅(qū)動端，C點(diǎn)沿X方向?yàn)樨浥撻T收放機(jī)構(gòu)的傳動線系洲BC為貨艙門收放機(jī)構(gòu)，G為貨艙門的重心，F(xiàn)點(diǎn)為風(fēng)載作用于貨艙門上點(diǎn)，ABCO構(gòu)成四桿機(jī)構(gòu)，貨艙門開啟過程中O、A兩點(diǎn)固定不變，B、C兩點(diǎn)隨收放機(jī)構(gòu)的運(yùn)動而轉(zhuǎn)動，C點(diǎn)的旋轉(zhuǎn)運(yùn)動實(shí)現(xiàn)貨艙門的收放。

1.2 貨艙門收放機(jī)構(gòu)運(yùn)動參數(shù)確定[4～6]

已知貨艙門收放機(jī)構(gòu)的質(zhì)量m=190kg，貨艙門尺寸為2739mm×2385mm，貨艙門的最大開啟角度為130°，貨艙門在完全關(guān)閉狀態(tài)時，貨艙門初始角度r₀=25.07°，門重心和豎直方向的夾角b₀=21.29°，風(fēng)載作用點(diǎn)和豎直方向的夾角c₀=18.12°，固定連桿長度OA=310mm，活動連桿長度AB=300mm、BC=292mm、OC=280mm，貨艙門重心到門轉(zhuǎn)軸的長度OG=1041.05mm，風(fēng)載作用點(diǎn)到門轉(zhuǎn)軸的長度為OF=1157.58mm。

根據(jù)上述參數(shù)，貨艙門關(guān)閉狀態(tài)到完全打開狀態(tài)AC長度不斷變化，在△OAC中根據(jù)余弦定理可得：式中：r為貨艙門打開過程中轉(zhuǎn)過的角度。

由此推導(dǎo)出四連桿中各角度計(jì)算公式為：

將式（1）、式（4）聯(lián)立可得：

貨艙門在重力作用下圍繞O產(chǎn)生力矩M_g為：

貨艙門在風(fēng)載作用下圍繞O產(chǎn)生的力矩為：

根據(jù)伯努力方程，貨艙門上承受的風(fēng)載為[7]：式中：地面空氣密度ρ=1.225kg/m³，風(fēng)速v=20.6m/s，貨艙門的最大迎風(fēng)面積S_m=6.5325m²。

2 貨艙門收放機(jī)構(gòu)的建模與仿真

2.1 電機(jī)輸出扭矩

電機(jī)輸出扭矩隨著艙門開啟角度及轉(zhuǎn)速在變化，電機(jī)功率扭矩計(jì)算公式為：

電機(jī)轉(zhuǎn)速與貨艙門轉(zhuǎn)速之間的關(guān)系為：

電機(jī)轉(zhuǎn)速與輸出端C點(diǎn)轉(zhuǎn)速關(guān)系為：式中：M_d為電機(jī)輸出扭矩（N·m）;P為電機(jī)額定功率，P=0.2kW;n_d為電機(jī)轉(zhuǎn)速（r/min）;ω_r為貨艙門轉(zhuǎn)速;u為負(fù)載轉(zhuǎn)速與貨艙門轉(zhuǎn)速比;i為收放機(jī)構(gòu)傳動比，i=3612.5。

2.2 貨艙門收放機(jī)構(gòu)的運(yùn)動微分方程

根據(jù)貨艙門的使用要求，貨艙門開啟的時間為30s。不考慮摩擦的情況下，電機(jī)啟動后，貨艙門開啟過程中圍繞口的運(yùn)動微分方程如下：

根據(jù)圖1貨艙門受力簡圖可知：式中：J為貨艙門繞點(diǎn)。的轉(zhuǎn)動慣量J=m×1.04105²kg·m²;M'為電機(jī)對貨艙門轉(zhuǎn)軸O點(diǎn)的力矩;M_g為貨艙門重力繞口點(diǎn)的力矩;M_f為貨艙門風(fēng)載繞O點(diǎn)的力矩;β為貨艙門圍繞O點(diǎn)轉(zhuǎn)動角加速度，ω為角速度（rad/s）;r為角位移（rad）;F₁為連桿AB軸力（N）;T₁為輸出端負(fù)載扭矩（N·m）。

2.3 收放機(jī)構(gòu)Simulink仿真模型建立[8～10]

Simulink是一個非常適用于工程的仿真平臺，能夠快速、準(zhǔn)確地建立仿真模型。貨艙門開啟過程的仿真模型如圖2所示，采用電機(jī)和收放機(jī)構(gòu)運(yùn)動的圖形化模型，將計(jì)算公式用Matlab Function模塊建立。電機(jī)仿真模型輸入?yún)?shù)ω_r，u;輸出參數(shù)T_d，n_d。收放機(jī)構(gòu)運(yùn)動仿真模型輸入?yún)?shù)r，M_d;輸出參數(shù)y，M_g，F(xiàn)₁，T₁，F(xiàn)_fM_f，M，a，v。需要注意的是，這里出現(xiàn)的T_d，y，v分別對應(yīng)于M_d，β，u。

對于存在約束狀態(tài)的變量，要附加適當(dāng)?shù)哪K或者對模塊進(jìn)行參數(shù)設(shè)置。例如貨艙門開啟角度為0°≤r≤130°，則將Saturation模塊參數(shù)的上限設(shè)為130×π/180，下限設(shè)為0;而對于電機(jī)輸出扭矩T_d的約束條件，則在輸出T_d的后面附加Switch模塊，Switch模塊中判斷條件選擇“>門限值”，門限值設(shè)為0.02×π/180，若不滿足判斷條件，選擇輸出Constant模塊中的常數(shù)0.1。

各模塊建立起來以后，將它們連接起來，設(shè)置SimulinkParameters，設(shè)置開始時間為0，終止時間為50s，采用算法類型為變步長，算法為ode45對系統(tǒng)進(jìn)行仿真計(jì)算，將需要輸出的信號連入Scope模塊或XY Graph，輸出其隨時間的變化曲線或其隨貨艙門開啟角度的變化曲線。

2.4 仿真結(jié)果

貨艙門收放機(jī)構(gòu)驅(qū)動電機(jī)扭矩、轉(zhuǎn)速、變化如圖3所示，負(fù)載扭矩變化如圖4所示，連桿AB軸力變化如圖5所示，輸出端速度、加速度變化如圖6所示，收放機(jī)構(gòu)開啟艙門所需時間曲線如圖7所示。

從仿真曲線可知：電機(jī)扭矩在0.7～1N·m之間，轉(zhuǎn)速在1380～1550r/min之間工作平穩(wěn);艙門開啟至130°時，連桿AB軸力最大為13000N;艙門開啟至111°時，負(fù)載扭矩最大為1889N·m;貨艙門收放機(jī)構(gòu)將貨艙門開啟過程呈現(xiàn)接近勻速運(yùn)動，艙門開啟至最大位置130°時，需要用時30s，符合實(shí)際情況。

3 結(jié)論

運(yùn)用Simulink系統(tǒng)建模與仿真技術(shù)，對貨艙門收放機(jī)構(gòu)載荷及運(yùn)動特性進(jìn)行仿真分析，具有方便、準(zhǔn)確和直觀的特點(diǎn)，解決了傳統(tǒng)方法計(jì)算繁瑣、設(shè)計(jì)周期長的問題，通過仿真分析，了解貨艙門收放機(jī)構(gòu)載荷及運(yùn)動變化特性，為貨艙門收放機(jī)構(gòu)的設(shè)計(jì)提供了參考、具有一定的工程意義，后續(xù)可就不同參數(shù)變化對收放機(jī)構(gòu)載荷及運(yùn)動特性影響進(jìn)行具體分析，通過調(diào)整參數(shù)、觀察曲線變化，可了解諸參數(shù)對機(jī)構(gòu)的影響，有利于選擇優(yōu)化參數(shù)，設(shè)計(jì)出合理的機(jī)構(gòu)。

參考文獻(xiàn)

[1]Thuis H G S J.Development of a composite cargo door for anaircraft[J].Composite Structures，1999，47：813-825.

[2]薛定宇.基于Matlab/Simulink的系統(tǒng)仿真技術(shù)與應(yīng)用[M].北京：清華大學(xué)出版社，2006.

[3]武娟，任家駿，吳鳳林，等.基于Matlab軟件中Simulink的機(jī)構(gòu)運(yùn)動仿真[J].機(jī)械工程與自動化，2004（3）4-6.

[4]中國民用航空局.CCAR-25-R4運(yùn)輸類飛機(jī)試航標(biāo)準(zhǔn)[S].北京：中國民用航空局，2011.

[5]張德豐，許華興，王旭.Matlab/Simulink建模與仿真實(shí)例精講[M].北京：機(jī)械工業(yè)出版社，2010.

[6]洪嘉振.理論力學(xué)[M].北京：高等教育出版社，2015.

[7]詹光，孫穎.某高空長航時飛機(jī)垂直突風(fēng)過載計(jì)算分析[J].飛機(jī)設(shè)計(jì)，2007，27（6）：7-8.

[8]余知樸，岳珠峰，劉永濤，等.飛行環(huán)境下飛機(jī)貨艙門電液系統(tǒng)的開啟功能分析[J].機(jī)械工程學(xué)報(bào)，2014（3）：51-52.

[9]邱曉林.基于Matlab-Simulink的動態(tài)模型與系統(tǒng)仿真工具[M].西安：西安交通大學(xué)出版社，2003.

[10]許成，王作祥，劉宗杰.基于Simulink的潛空導(dǎo)彈運(yùn)載器水下彈道仿真[J].海軍航空工程學(xué)院學(xué)報(bào)，2012，27（4）：5.