主权项 |
一种转动部件对挠性动力学影响的仿真方法,包括:步骤1)、对大型转动部件的特性进行分析,建立所述大型转动部件的动力学模型;步骤2)、对转动部件所带的挠性附件进行结构动力学分析,经过模态截断后生成线性化的模态坐标下描述的振动模型,结合刚体动力学模型建立带有挠性附件的多体动力学模型;步骤3)、在步骤1)和步骤2)所建立的模型的基础上建立挠性卫星姿态动力学与控制系统模型,利用该模型进行仿真;在所述的步骤1)中,所述大型转动部件的动力学模型为所述转动部件的动不平衡量、静不平衡量对整星质心的干扰力F和力矩T; <mrow> <mi>F</mi> <mo>=</mo> <msubsup> <mi>ω</mi> <mi>r</mi> <mn>2</mn> </msubsup> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>B</mi> <mi>j</mi> </msub> <mi>cos</mi> <msub> <mi>β</mi> <mi>j</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>B</mi> <mi>j</mi> </msub> <mi>sin</mi> <msub> <mi>β</mi> <mi>j</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> <mi>T</mi> <mo>=</mo> <msubsup> <mi>ω</mi> <mi>r</mi> <mn>2</mn> </msubsup> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>B</mi> <mi>d</mi> </msub> <mi>sin</mi> <msub> <mi>β</mi> <mi>d</mi> </msub> <mo>-</mo> <msub> <mi>z</mi> <mi>a</mi> </msub> <msub> <mi>B</mi> <mi>j</mi> </msub> <mi>cos</mi> <msub> <mi>β</mi> <mi>j</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>B</mi> <mi>d</mi> </msub> <mi>cos</mi> <msub> <mi>β</mi> <mi>d</mi> </msub> <mo>+</mo> <msub> <mi>z</mi> <mi>a</mi> </msub> <msub> <mi>B</mi> <mi>j</mi> </msub> <mi>sin</mi> <msub> <mi>β</mi> <mi>j</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>x</mi> <mi>a</mi> </msub> <msub> <mi>B</mi> <mi>j</mi> </msub> <mi>cos</mi> <msub> <mi>β</mi> <mi>j</mi> </msub> <mo>-</mo> <msub> <mi>y</mi> <mi>a</mi> </msub> <msub> <mi>B</mi> <mi>j</mi> </msub> <mi>sin</mi> <msub> <mi>β</mi> <mi>j</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> </mrow>其中,Bj为静不平衡量大小,βj为静不平衡量相位,Bd为动不平衡量大小,βd为动不平衡量相位,[xa,ya,za]为参考点在本体系坐标,ωr为转动部件转速;所述的多体动力学模型包括: <mrow> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mi>m</mi> <msub> <mover> <mi>v</mi> <mo>·</mo> </mover> <mi>c</mi> </msub> <mo>+</mo> <mi>BT</mi> <mover> <mi>η</mi> <mrow> <mo>·</mo> <mo>·</mo> </mrow> </mover> <mo>=</mo> <mi>F</mi> </mtd> </mtr> <mtr> <mtd> <msub> <mi>J</mi> <mi>o</mi> </msub> <mo>·</mo> <mover> <mi>ω</mi> <mo>·</mo> </mover> <mo>+</mo> <mi>ω</mi> <mo>×</mo> <mrow> <mo>(</mo> <msub> <mi>J</mi> <mi>o</mi> </msub> <mo>·</mo> <mi>ω</mi> <mo>)</mo> </mrow> <mo>+</mo> <mi>BR</mi> <mover> <mi>η</mi> <mrow> <mo>·</mo> <mo>·</mo> </mrow> </mover> <mo>=</mo> <mi>T</mi> </mtd> </mtr> <mtr> <mtd> <mover> <mi>η</mi> <mrow> <mo>·</mo> <mo>·</mo> </mrow> </mover> <mo>+</mo> <mn>2</mn> <mi>ξΩ</mi> <mover> <mi>η</mi> <mo>·</mo> </mover> <mo>+</mo> <msup> <mi>Ω</mi> <mn>2</mn> </msup> <mi>η</mi> <mo>+</mo> <msup> <mi>BR</mi> <mi>T</mi> </msup> <mover> <mi>ω</mi> <mo>·</mo> </mover> <mo>+</mo> <msup> <mi>BT</mi> <mi>T</mi> </msup> <msub> <mover> <mi>v</mi> <mo>·</mo> </mover> <mi>c</mi> </msub> <mo>=</mo> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> </mrow>其中,ω为角速度,vc为质心线速度,η为模态变量,Ω为挠性附件振动基频,ξ为挠性附件振动阻尼系数,BR、BT分别为挠性附件相对质心的转动、平动耦合系数,F、T分别为转动部件动、静不平衡引起干扰力和力矩。 |