日韩精品一区蜜桃av在线-久久久久久久精品成人热入口-日本中文字幕二区1080p-久久久久女人精品毛片九一

聲明：本文為《現代防御技術(shù)》雜志社供《中國軍工網(wǎng)》獨家稿件。未經(jīng)許可，請勿轉載。

作者簡(jiǎn)介：高慶豐（1979-），男，內蒙古呼和浩特人，助工，碩士，主要從事導彈總體技術(shù)研究。
通信地址：100854北京142信箱30分箱
高慶豐1，劉莉2，陳羅婧2
（1．中國航天科工集團公司 二院二部，北京100854；2.北京理工大學(xué) 飛行器工程系，北京100081）

摘要：在彈體坐標系和準彈體坐標系中建立了旋轉飛行器角運動(dòng)數學(xué)模型。應用李亞普諾夫第一近似理論和勞斯-霍爾維茨方法導出了旋轉飛行器的非線(xiàn)性運動(dòng)穩定性判據，這個(gè)判據可應用于有控旋轉導彈的運動(dòng)穩定性分析，也可應用到炮彈和火箭彈上。
關(guān)鍵詞：旋轉飛行器；非線(xiàn)性；穩定性
中圖分類(lèi)號：V412；TJ415；O242.2 TJ7611＋3;文獻標識碼：A文章編號：1009086X（2006）01001905

Criteria for the nonlinear dynamic stability of rotative vehicles
GAO Qingfeng1，LIU Li2，CHEN Luojing2
（1．The Second System Design Department of the Second Research Academy of CASIC，Beijing 100854，China;
2．Beijing Institute of Technology，Department of Flight Vehicle Engineering，Beijing 100081，China）

Abstract:The angular motion mathematical model of rotative vehicles is established in the body coordinate system and the quasibody coordinate system.Criteria for the nonlinear dynamic stability of rotative vehicles are derived by the Liapunov′s first method and the RouthHurwitz stability criterion,this criteria can be applied to the dynamic stability analysis of controlled rotative missiles,projectiles and rockets.
Key words:Rotative vehicle；Nonlinear；Stability

1引言
旋轉飛行器是指在飛行過(guò)程中，繞其縱軸自旋的一類(lèi)飛行器，通常包括小型防空導彈、反坦克導彈、炮彈和火箭彈等。
反坦克導彈有無(wú)控的起始飛行段，對這類(lèi)導彈進(jìn)行設計時(shí)，要對其彈體的動(dòng)態(tài)特性提出穩定性要求。如果用經(jīng)典的方法設計制導系統，也往往首先要研究彈體穩定性問(wèn)題。而對于炮彈和火箭彈，運動(dòng)穩定性更是首要的［1］。
2符號說(shuō)明
a1，b1為與升力和側向力有關(guān)的動(dòng)力系數；a2，b2為與馬格努斯力有關(guān)的動(dòng)力系數；a3，b3為與俯仰力矩和偏航力矩有關(guān)的動(dòng)力系數；a4，b4為與馬格努斯力矩有關(guān)的動(dòng)力系數；a5，b5為與阻尼力矩有關(guān)的動(dòng)力系數；a6，b6為與轉速有關(guān)的動(dòng)力系數；α，β為彈體坐標系中的攻角和側滑角；αf，βf準彈體坐標系中的攻角和側滑角；ωz1，ωy1為彈體坐標系中的俯仰和偏航角速度；ωzf，ωyf為準彈體坐標系中的俯仰和偏航角速度；P0為發(fā)動(dòng)機推力；q為動(dòng)壓；S為參考面積；L為參考長(cháng)度；m為質(zhì)量；v為速度；ωx為彈體繞縱軸的旋轉角速度；ξ~為復攻角；cαy，cα3y為線(xiàn)性和立方升力系數導數；cαz，cα3z為線(xiàn)性和立方側向力系數導數；cβy，cβz為氣動(dòng)交叉力系數導數；mαz，mα3z為線(xiàn)性和立方俯仰力矩系數導數；mβy，mβ3y為線(xiàn)性和立方偏航力矩系數導數；mαy，mβz為馬格努斯力矩系數導數；mωzz，mωyy為俯仰阻尼和偏航阻尼力矩系數；Jx，Jy，Jz為相對彈體坐標系各軸的轉動(dòng)慣量；t為時(shí)間；s為彈道弧長(cháng)；σ為穩定性系數。
3旋轉飛行器角運動(dòng)數學(xué)模型
有一類(lèi)小型防空導彈，彈體在飛行中以一定的角速度繞自身縱軸旋轉，采用單通道控制，由一對舵面同時(shí)控制導彈的俯仰運動(dòng)和偏航運動(dòng)，因此其氣動(dòng)外形是面對稱(chēng)的。
為了便于討論，建立與彈體固聯(lián)的彈體坐標系Ox1y1z1和不隨彈體旋轉的準彈體坐標系Oxfyf zf，它們都以彈體質(zhì)心為坐標原點(diǎn)，彈體坐標系Ox1軸與彈體縱軸重合，向前為正，Oy1軸垂直于Ox1軸及舵軸。Oz1軸與Ox1軸和Oy1軸形成右手系。準彈體坐標系Oxf軸與Ox1軸重合，Oyf 軸垂直于Oxf 軸指向上，Ozf 軸與Oxf 軸和Oyf軸形成右手系。
現代防御技術(shù)·導彈技術(shù)高慶豐，劉莉，陳羅婧：旋轉飛行器非線(xiàn)性運動(dòng)穩定性判據現代防御技術(shù)2006年第34卷第1期忽略重力的影響，以彈體坐標系表征的自由運動(dòng)中力和力矩的平衡方程分別為［2］α·
β·＋a1〖〗a2＋ωx
－（b2＋ωx）〖〗b1α
β－ωz1
ωy1＝0，（1）
ω·z1
ω·y1＝a3〖〗a4
－b4〖〗b3α
β＋
ωx－a6〖〗a5
b5〖〗－（b6－ωx）ωy1
ωz1，（2）式（1）和式（2）中：a1＝［P0＋qS（cαy＋cα3yα2）］/mv， b1＝［P0－qS（cαz＋cα3zα2）］/mv，
a2＝qScβz/mv， b2＝qScβy/mv，
a3＝qSL(mαz＋mα3z)/Jz, b3＝qSL（mβy＋mβ3y）/Jy，
a4＝qSLmαy/Jz， b4＝qSLmβz/Jy，
a5＝qSL2mωzz/Jzv， b5＝qSL2mωyy/Jyv，
a6＝（Jx/Jz）ωx， b6＝（Jx/Jy）ωx 通過(guò)坐標變換，以準彈體坐標系表征的自由運動(dòng)中力和力矩的平衡方程分別為［2］ω·yf
ω·zf＋－（a5＋b5）〖〗2＋（a6－b6）〖〗2sin（2ωxt）＋（a5－b5）〖〗2cos（2ωxt）〖〗（a6＋b6）〖〗2＋（a5－b5）〖〗2sin（2ωxt）－（a6－b6）〖〗2cos（2ωxt）
－（a6＋b6）〖〗2＋（a5－b5）〖〗2sin（2ωxt）＋（a6－b6）〖〗2cos（2ωxt）〖〗－（a5＋b5）〖〗2＋（b6－a6）〖〗2sin（2ωxt）＋（a5－b5）〖〗2cos（2ωxt）
ωyf
ωzf＋（a4＋b4）〖〗2＋（a3－b3）〖〗2sin（2ωxt）－（a4－b4）〖〗2cos（2ωxt）〖〗－（a3＋b3）〖〗2－（b3－a3）〖〗2cos（2ωxt）
－（a3＋b3）〖〗2－（a3－b3）〖〗2cos（2ωxt）〖〗－（a4＋b4）〖〗2＋（a3－b3）〖〗2sin（2ωxt）－（a4－b4）〖〗2cos（2ωxt）αf
βf＝0，（3）
α·f
β·f＋（a1＋b1）〖〗2－（b1－a1）〖〗2cos（ωxt）〖〗（a2＋b2）〖〗2＋（b1－a1）〖〗2sin（2ωxt）＋（a2－b2）〖〗2cos（2ωxt）
－（a2＋b2）〖〗2＋（b1－a1）〖〗2sin（2ωxt）－（a2－b2）〖〗2cos（2ωxt）〖〗（a1＋b1）〖〗2－（b1－a1）〖〗2cos（ωxt）·
αf
βf＋0-1
-10ωyf
ωzf＝0 （4）對于面對稱(chēng)導彈，由于a1≠b1，a2≠b2，a3≠b3，a4≠b4，a5≠b5，a6≠b6，因而含有sin（2ωxt）和cos（2ωxt）項的系數不為0，角頻率為2ωx的擺動(dòng)將不可避免的存在，考慮到擺動(dòng)的幅值不大，且在彈體旋轉一周所產(chǎn)生的平均效應為0。因此，當彈體旋轉頻率遠大于彈體擾動(dòng)運動(dòng)頻率時(shí)，完全可把含sin（2ωxt）和cos（2ωxt）的項略去不計［2］。所以，式（3）可變?yōu)槭剑?），式（4）可變?yōu)槭剑?）。ω·yf
ω·zf＝a5＋b5〖〗2〖〗－（a6＋b6）〖〗2
a6＋b6〖〗2〖〗a5＋b5〖〗2ωyf
ωzf＋－（a4＋b4）〖〗2〖〗a3＋b3〖〗2
a3＋b3〖〗2〖〗－（a4＋b4）〖〗2αf
βf，（5）
α·f
β·f＋a1＋b1〖〗2〖〗a2＋b2〖〗2
－（a2＋b2）〖〗2〖〗a1+b1〖〗2αf
βf＋0〖〗-1
-1〖〗0ωyf
ωzf＝0 （6）對式（6）求導，可得ω·yf
ω·zf＝β¨f
α¨f＋－（a2＋b2）〖〗2〖〗a1＋b1〖〗2
a1＋b1〖〗2〖〗a2＋b2〖〗2α·f
β·f,（7）令式（5）和式（7）右端相等，同時(shí)代入式（6）可得α¨f＋（a1＋b1）－（a5＋b5）〖〗2α·f－（a3＋b3）〖〗2＋（a1＋b1）（a5＋b5）〖〗4－（a2＋b2）（a6＋b6）〖〗4αf－
－（a2＋b2）＋（a6＋b6）〖〗2β·f－（a4＋b4）〖〗2＋（a1＋b1）（a6＋b6）〖〗4＋（a2＋b2）（a5＋b5）〖〗4βf＝0,（8）
β¨f＋（a1＋b1）－（a5＋b5）〖〗2β·f－（a3＋b3）〖〗2＋（a1＋b1）（a5＋b5）〖〗4－（a2＋b2）（a6＋b6）〖〗4βf+
－（a2＋b2）＋（a6＋b6）〖〗2α·f＋（a4＋b4）〖〗2＋（a1＋b1）（a6＋b6）〖〗4＋（a2＋b2）（a5＋b5）〖〗4αf＝0 （9）將式（8）乘以虛數i再與式（9）相加，可得復數表達式ξ~¨＋（A－iB）ξ~·－（C＋iD）ξ~＝0,（10）式（10）中，定義A＝［（a1＋b1）－（a5＋b5）］/2,
B＝［－（a2＋b2）＋（a6＋b6）］/2,
C＝（a3＋b3）/2＋（a1＋b1）（a5＋b5）/4－（a2＋b2）（a6＋b6）/4,
D＝（a4＋b4）/2＋（a1＋b1）（a6＋b6）/4＋（a2＋b2）（a5＋b5）/4,
ξ~＝βf＋iαf 式（10）為復攻角在t域的微分方程。
由于式（10）的系數與速度有關(guān)，這是一變系數微分方程，為使t域角運動(dòng)微分方程系數的時(shí)變性減弱，對式(10)進(jìn)行數學(xué)變換，有［3］df（x）〖〗dt＝df（x）〖〗dsds〖〗dt＝vdf（x）〖〗ds，
df（x）2〖〗d2t＝v2df（x）2〖〗d2s＋v·df（x）〖〗ds （11）應用式（11），將式（10）變?yōu)棣蝵″＋［（v·/v＋A－iB）/v］ξ~′－
（1/v）2（C＋iD）ξ~＝0，（12）式（12）中，定義H＝（v·/v＋A）/v,
P＝B/v,
M＝C/v2,
PT＝D/v2,式（12）可寫(xiě)為ξ~″＋（H－iP）ξ~′－（M＋iPT）ξ~＝0 （13）式（13）為復攻角在s域的微分方程，將H，P，M，PT展開(kāi)，并忽略氣動(dòng)交叉項的影響，J＝Jy≈Jz，可得H（α2）＝ρS〖〗2m（cαy＋cα3yα2）－（cαz＋cα3zα2）〖〗2－mL2〖〗J（mωzz＋mωyy）〖〗2－cx－2mgsin θ〖〗ρSv2＋4P0〖〗ρSv2，
M（α2）＝ρSL〖〗2J（mαz＋mα3zα2）＋（mβy＋mβ3yα2）〖〗2＋P0L〖〗mv2（mωzz＋mωyy）〖〗2－Jxωx〖〗mvL（cβy＋cβz）〖〗2,
P＝Jx〖〗Jωx〖〗v,
T（α2）＝ρSLv〖〗4Jxωx（mαy＋mβz）＋ρS〖〗4m［（cαy＋cα3yα2）－（cαz＋cα3zα2）］＋P0〖〗mv2 反坦克導彈、炮彈和火箭彈為軸對稱(chēng)旋轉飛行器，軸對稱(chēng)旋轉飛行器是面對稱(chēng)旋轉飛行器的特例。對于軸對稱(chēng)旋轉飛行器，有a1＝b1，a2＝b2，a3＝b3，a4＝b4，〖〗a5＝b5，a6＝b6，復攻角在s域的微分方程也為ξ~″＋（H－iP）ξ~′－（M＋iPT）ξ~＝0,（14）式（14）中:H（α2）＝ρS〖〗2m（cαy＋cα3yα2）－mL2〖〗Jmωzz－cx－2mgsin θ〖〗ρSv2＋4P0〖〗ρSv2,
M（α2）＝ρSL〖〗2J（mαz＋mα3zα2）＋P0L〖〗mv2mωzz－Jxωx〖〗mvLcβz,
P＝Jx〖〗Jωx〖〗v,
T（α2）＝ρSLv〖〗2Jxωxmαy＋ρS〖〗2m（cαy＋cα3yα2）＋P0〖〗mv2 4旋轉飛行器運動(dòng)穩定性判據
復攻角非線(xiàn)性微分方程（13）在實(shí)數域的攻角非線(xiàn)性微分方程組為α″f＋H（α2）α′f－M（α2）αf＋Pβ′f＋PT（α2）βf＝0,
β″f＋H（α2）β′f－M（α2）βf－Pα′f-PT（α2）αf＝0 （15）非線(xiàn)性微分方程組（15）是一個(gè)含有3個(gè)非線(xiàn)性函數，1個(gè)常數的四階微分方程組，判斷其解的穩定性十分困難［4］。應用李亞普諾夫第一近似理論，對于非線(xiàn)性微分方程組，如果其線(xiàn)性化微分方程組之特征方程的所有特征根均有負實(shí)部，則非線(xiàn)性微分方程組的原點(diǎn)漸進(jìn)穩定［5］。所以，可通過(guò)線(xiàn)性化

微分方程組的穩定性來(lái)判斷非線(xiàn)性微分方程組的穩定性。非線(xiàn)性微分方程組（15）的線(xiàn)性化微分方程組為α″f＋Hα′f－Mαf＋Pβ′f＋PTβf＝0,
β″f＋Hβ′f－Mβf－Pα′f-PTαf＝0,（16）式（16）中:H＝ρS〖〗2m（cαy－cαz）〖〗2－mL2〖〗J（mωzz＋mωyy）〖〗2－cx－2mgsin θ〖〗ρSv2＋4P0〖〗ρSv2,
M＝ρSL〖〗2J（mαz＋mβy）〖〗2＋P0L〖〗mv2（mωzz＋mωyy）〖〗2－Jxωx〖〗mvL（cβy＋cβz）〖〗2,
P＝Jx〖〗Jωx〖〗v,
T＝ρSLv〖〗4Jxωx（mαy＋mβz）＋ρS〖〗4m（cαy－cαz）＋P0〖〗mv2 
首先，通過(guò)數學(xué)變換（11），微分方程組（16）的系數時(shí)變性減弱了。其次，采用系數凍結法，在不長(cháng)的彈道區間內，微分方程組（16）的系數可視為常數。所以，線(xiàn)性化微分方程組（16）可看作線(xiàn)性定常系統，其特征方程為λ4＋h1λ3＋h2λ2＋h3λ＋h4＝0,（17）式（17）中:h1＝2H；h2＝P2＋H2－2M；h3＝2（P2T－MH）；h4＝M2＋（PT）2。
由勞斯-霍爾維茨方法可知，特征方程（17）全部根的實(shí)部都為負值的充要條件是下列條件成立［6］：h1＞0，h2＞0，h3＞0，h4＞0,
h1h2-h3＞0,
h1h2-h3＞h21h4/h3 （18）根據條件（18)，代入hi關(guān)系式后，整理可得H＞0,
P2T－MH＞0,
（P2＋H2/P2）［H（P2T－MH）－（PT）2］＞0 （19）式（19）中，第二式由第一、第三式成立而自然滿(mǎn)足，與式（19）等價(jià)的條件變?yōu)镠＞0,
H（P2T－MH）－（PT）2＞0 （20）式（20）中的第二式可化為0＜σ＝（2PT－PH）2〖〗H2（P2－4M）＜1 （21） 所以，面對稱(chēng)旋轉飛行器運動(dòng)穩定性判據為0＜σ＜1，由σ不僅可判斷是否滿(mǎn)足運動(dòng)穩定性判據，而且能夠反映出穩定性的好壞，σ越小，穩定性越好［7］。同理，可得到軸對稱(chēng)旋轉飛行器的運動(dòng)穩定性判據也為0＜σ＝（2PT－PH）2〖〗H2（P2－4M）＜1,（22）式（22）中：

H＝ρS〖〗2mcαy－mL2〖〗Jmωzz－cx－2mgsin θ〖〗ρSv2＋4P0〖〗ρSv2，
M＝ρSL〖〗2Jmαz＋P0L〖〗mv2mωzz－Jxωx〖〗mvLcβz,
P＝Jx〖〗Jωx〖〗v,
T＝ρSLv〖〗2Jxωxmαy＋ρS〖〗2mcαy+P0〖〗mv2 
5結束語(yǔ)
本文得到的旋轉飛行器非線(xiàn)性運動(dòng)穩定性判據具有很強的通用性，對面對稱(chēng)和軸對稱(chēng)的旋轉飛行器都是適用的，不但適合于分析旋轉導彈和主動(dòng)段火箭彈，也適合于分析炮彈和被動(dòng)段火箭彈。
參考文獻：
［1］康景利.變參量彈體運動(dòng)穩定性判據［J］.北京理工大學(xué)學(xué)報，1989，9（2）：35-40.
［2］葉堯卿.便攜式紅外尋的防空導彈設計［M］.北京：宇航出版社，1996.
［3］閆章更，祁載康.射表技術(shù)［M］.北京：國防工業(yè)出版社，2000.
［4］張迅.炮彈非線(xiàn)性運動(dòng)穩定性的勞斯－霍維茨方法［J］.彈箭與制導學(xué)報，1997，（2）：39-43
［5］高為柄.運動(dòng)穩定性基礎［M］ .北京：高等教育出版社，1987.
［6］胡壽松.自動(dòng)控制原理［M］ .北京：國防工業(yè)出版社，1994
［7］周效亮.旋轉彈丸的飛行穩定性設計［J］.兵工學(xué)報彈箭分冊，1983，（4）：8-24.
信度的射擊精度評估方法。事實(shí)上，如何準確地建立武器系統動(dòng)力學(xué)模型，以及模型的VV&A，從而獲得驗前信息，也是一個(gè)關(guān)鍵性的問(wèn)題。
參考文獻：
［1］王國平,芮筱亭,陳衛東.多管武器系統密集度仿真技術(shù)［J］.系統仿真學(xué)報,2004,(5):963-966
［2］芮筱亭,陸毓琪,王國平，等.多管火箭發(fā)射動(dòng)力學(xué)仿真與試驗測試方法［M］.北京:國防工業(yè)出版社,2003
［3］李鵬波,謝紅衛,張金槐.考慮驗前信息可信度時(shí)的Bayes估計［J］.國防科技大學(xué)學(xué)報, 2003，25(4):107-110
［4］張士峰.多元驗前信息的融合方法［J］.飛行器測控學(xué)報,2000,19(1):26-302006年2月〖〗第34卷第1期現代防御技術(shù)〖〗MODERN DEFENCE TECHNOLOGYFeb.