发明名称 一种转动部件对挠性动力学影响的仿真方法
摘要 本发明公开了一种转动部件对挠性动力学影响的仿真方法,包括:对大型转动部件的特性进行分析,建立所述大型转动部件的动力学模型;对转动部件所带的挠性附件进行结构动力学分析,经过模态截断后生成线性化的模态坐标下描述的振动模型,结合刚体动力学模型建立带有挠性附件的多体动力学模型;在之前建立的模型的基础上建立挠性卫星姿态动力学与控制系统模型,利用该模型进行仿真。该方法可用于姿轨控与动力学设计过程中对设计结果进行验证,也可用于在轨运行期间的故障判定和反演。本发明解决了转动部件对动力学影响分析的实际工程问题,取得了提高平台稳定性和可靠性的有益效果。
申请公布号 CN102609561B 申请公布日期 2014.03.05
申请号 CN201110425575.4 申请日期 2011.12.15
申请人 上海卫星工程研究所 发明人 吕旺;薛孝补;杨立峰;朱海江
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 中国和平利用军工技术协会专利中心 11215 代理人 容敦璋
主权项 一种转动部件对挠性动力学影响的仿真方法,包括:步骤1)、对大型转动部件的特性进行分析,建立所述大型转动部件的动力学模型;步骤2)、对转动部件所带的挠性附件进行结构动力学分析,经过模态截断后生成线性化的模态坐标下描述的振动模型,结合刚体动力学模型建立带有挠性附件的多体动力学模型;步骤3)、在步骤1)和步骤2)所建立的模型的基础上建立挠性卫星姿态动力学与控制系统模型,利用该模型进行仿真;在所述的步骤1)中,所述大型转动部件的动力学模型为所述转动部件的动不平衡量、静不平衡量对整星质心的干扰力F和力矩T; <mrow> <mi>F</mi> <mo>=</mo> <msubsup> <mi>&omega;</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>&beta;</mi> <mi>j</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>B</mi> <mi>j</mi> </msub> <mi>sin</mi> <msub> <mi>&beta;</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>&omega;</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>&beta;</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>&beta;</mi> <mi>j</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>B</mi> <mi>d</mi> </msub> <mi>cos</mi> <msub> <mi>&beta;</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>&beta;</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>&beta;</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>&beta;</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>&CenterDot;</mo> </mover> <mi>c</mi> </msub> <mo>+</mo> <mi>BT</mi> <mover> <mi>&eta;</mi> <mrow> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> </mrow> </mover> <mo>=</mo> <mi>F</mi> </mtd> </mtr> <mtr> <mtd> <msub> <mi>J</mi> <mi>o</mi> </msub> <mo>&CenterDot;</mo> <mover> <mi>&omega;</mi> <mo>&CenterDot;</mo> </mover> <mo>+</mo> <mi>&omega;</mi> <mo>&times;</mo> <mrow> <mo>(</mo> <msub> <mi>J</mi> <mi>o</mi> </msub> <mo>&CenterDot;</mo> <mi>&omega;</mi> <mo>)</mo> </mrow> <mo>+</mo> <mi>BR</mi> <mover> <mi>&eta;</mi> <mrow> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> </mrow> </mover> <mo>=</mo> <mi>T</mi> </mtd> </mtr> <mtr> <mtd> <mover> <mi>&eta;</mi> <mrow> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> </mrow> </mover> <mo>+</mo> <mn>2</mn> <mi>&xi;&Omega;</mi> <mover> <mi>&eta;</mi> <mo>&CenterDot;</mo> </mover> <mo>+</mo> <msup> <mi>&Omega;</mi> <mn>2</mn> </msup> <mi>&eta;</mi> <mo>+</mo> <msup> <mi>BR</mi> <mi>T</mi> </msup> <mover> <mi>&omega;</mi> <mo>&CenterDot;</mo> </mover> <mo>+</mo> <msup> <mi>BT</mi> <mi>T</mi> </msup> <msub> <mover> <mi>v</mi> <mo>&CenterDot;</mo> </mover> <mi>c</mi> </msub> <mo>=</mo> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> </mrow>其中,ω为角速度,vc为质心线速度,η为模态变量,Ω为挠性附件振动基频,ξ为挠性附件振动阻尼系数,BR、BT分别为挠性附件相对质心的转动、平动耦合系数,F、T分别为转动部件动、静不平衡引起干扰力和力矩。
地址 200240 上海市闵行区华宁路251号