发明名称 一种使用多级隔振平台提高光学载荷成像质量的方法
摘要 本发明涉及一种提高光学载荷成像质量的方法,特别涉及一种使用多级隔振平台提高光学载荷成像质量的方法,属于高频振动控制领域。一种隔振平台包括:上平台、下平台以及连接上平台和下平台的支杆。所述的一种隔振平台可以安装在飞轮或者控制力矩陀螺等卫星执行机构和卫星星体之间,或光学有效载荷和卫星星体之间;控制力矩陀螺、安装架和隔振平台的上平台共同组成了上平台系统;当需要在光学有效载荷和卫星星体加装隔振平台时,光学有效载荷的底座和隔振平台的上平台固连,同样能够组成一套上平台系统。所有隔振平台的下平台都固定连接在卫星星体上,共同组成了下平台系统。解决了星上光学有效载荷成像精度和稳定度低的问题。
申请公布号 CN102759927B 申请公布日期 2014.05.14
申请号 CN201210275120.3 申请日期 2012.08.03
申请人 北京理工大学 发明人 张尧;张景瑞;翟光
分类号 G05D1/08(2006.01)I;G05B13/04(2006.01)I 主分类号 G05D1/08(2006.01)I
代理机构 代理人
主权项 1.一种使用多级隔振平台提高光学载荷成像质量的方法,其特征在于包括如下步骤:步骤一:根据卫星所携带的飞轮(5)或者控制力矩陀螺(6)的个数来确定使用隔振平台(10)的个数,做到每个飞轮(5)或者每个控制力矩陀螺(6)都使用对应大小的隔振平台(10)来隔离它们造成的高频振动;步骤二:选用任意一个隔振平台(10),该隔振平台(10)共有N个支杆,求解第i个支杆在惯性坐标系下表示的支杆向量S<sub>i</sub>,其具体的表达式如下:S<sub>i</sub>=t+A<sub>eu</sub>p<sub>i</sub>-(b+r<sub>db</sub>+A<sub>ed</sub>q<sub>i</sub>)(1)式(1)中t和b表示惯性坐标系中心分别到上平台系统坐标系中心、卫星星体坐标系中心的矢量列阵;A<sub>eu</sub>和A<sub>ed</sub>是分别从上平台系统坐标系和下平台坐标系到卫星星体坐标系的坐标转换矩阵;p<sub>i</sub>和q<sub>i</sub>表示上平台质心分别到上平台第i个连接点处和下平台第i个连接点处的矢量列阵;r<sub>db</sub>是从卫星星体坐标系原点到下平台坐标系原点的矢量列阵;步骤三:对式(1)求模,可以解得第i个支杆(2)的杆长l<sub>i</sub>,表达形式如下l<sub>i</sub>=||S<sub>i</sub>||(2)步骤四:第i个支杆方向的单位向量定义为s<sub>ui</sub>,计算s<sub>ui</sub>,其表达形式如下所示s<sub>ui</sub>=S<sub>i</sub>/l<sub>i</sub>(3)步骤五:求取第i个支杆的伸缩速度l<sub>i</sub>,其表达形式分别如下<maths num="0001"><![CDATA[<math> <mrow> <msub> <mover> <mi>l</mi> <mo>&CenterDot;</mo> </mover> <mi>i</mi> </msub> <mo>=</mo> <msubsup> <mi>s</mi> <mi>ui</mi> <mi>T</mi> </msubsup> <mo>=</mo> <mo>[</mo> <msub> <mi>v</mi> <mi>p</mi> </msub> <mo>+</mo> <msub> <mi>A</mi> <mi>eu</mi> </msub> <msubsup> <mi>&omega;</mi> <mi>b</mi> <mo>&times;</mo> </msubsup> <msub> <mi>p</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>A</mi> <mi>eb</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>v</mi> <mi>b</mi> </msub> <mo>+</mo> <msubsup> <mi>&omega;</mi> <mi>b</mi> <mo>&times;</mo> </msubsup> <msub> <mi>r</mi> <mi>ab</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>A</mi> <mi>ed</mi> </msub> <msubsup> <mi>&omega;</mi> <mi>b</mi> <mo>&times;</mo> </msubsup> <msub> <mi>q</mi> <mi>i</mi> </msub> <mo>]</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow></math>]]></maths>式(4)中A<sub>eb</sub>是从卫星星体坐标系到惯性坐标系的坐标转换矩阵;v<sub>p</sub>和ω<sub>p</sub>分别是上平台系统的速度和角速度;v<sub>b</sub>和ω<sub>b</sub>分别是卫星星体的速度和角速度;上标“×”表示矢量列阵的反对称斜方阵;步骤六:每个支杆中间段(4)被简化成弹簧和阻尼并联的一个系统,令k<sub>i</sub>表示它的刚度参数;c<sub>i</sub>表示它的阻尼参数;则可求取隔振平台第i个支杆作用在上平台连接点处的力F<sub>si</sub>,其具体表达如下:<maths num="0002"><![CDATA[<math> <mrow> <msub> <mi>F</mi> <mi>si</mi> </msub> <mo>=</mo> <msub> <mi>k</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>l</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>l</mi> <mrow> <mn>0</mn> <mi>i</mi> </mrow> </msub> <mo>)</mo> </mrow> <msub> <mi>s</mi> <mi>ui</mi> </msub> <mo>+</mo> <msub> <mi>c</mi> <mi>i</mi> </msub> <msub> <mover> <mi>l</mi> <mo>&CenterDot;</mo> </mover> <mi>i</mi> </msub> <msub> <mi>s</mi> <mi>ui</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow></math>]]></maths>式中,l<sub>i</sub>是隔振平台(10)第i个支杆的杆长;l<sub>0i</sub>是第i个支杆(2)的标称长度;在计算F<sub>si</sub>的时候,需将式(1)、式(2)、式(3)式(4)的表达式带入到式(5)中;隔振平台(10)第i个支杆(2)作用在下平台连接点处的力F<sub>ui</sub>和第i个支杆(2)作用在上平台(1)连接点处的力F<sub>si</sub>互为反作用力,因此第i个支杆(2)作用在下平台(3)连接点处的力F<sub>ui</sub>如下式所示F<sub>ui</sub>=-F<sub>si</sub>(<sub>6</sub>)步骤七:由得出的各个支杆对上平台(1)连接点和对下平台(3)连接点的力,根据力的相互作用原理,可推得出含有一套隔振平台(10)的整星动力学方程;该方程包含了上平台系统的动力学方程和下平台系统的动力学方程,具体如下所示:<maths num="0003"><![CDATA[<math> <mrow> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <msub> <mi>m</mi> <mi>p</mi> </msub> <msub> <mover> <mi>v</mi> <mo>&CenterDot;</mo> </mover> <mi>p</mi> </msub> <mo>-</mo> <msubsup> <mi>S</mi> <mi>p</mi> <mo>&times;</mo> </msubsup> <msub> <mover> <mi>&omega;</mi> <mo>&CenterDot;</mo> </mover> <mi>p</mi> </msub> <mo>=</mo> <mo>-</mo> <msubsup> <mi>A</mi> <mi>eu</mi> <mi>T</mi> </msubsup> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msub> <mi>F</mi> <mi>si</mi> </msub> <mo>+</mo> <msub> <mi>F</mi> <mi>d</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>I</mi> <mi>p</mi> </msub> <mrow> <msub> <mover> <mi>&omega;</mi> <mo>&CenterDot;</mo> </mover> <mi>p</mi> </msub> <mo>+</mo> <msubsup> <mi>&omega;</mi> <mi>p</mi> <mo>&times;</mo> </msubsup> <msub> <mi>I</mi> <mi>p</mi> </msub> <msub> <mi>&omega;</mi> <mi>p</mi> </msub> <mo>+</mo> <msubsup> <mi>S</mi> <mi>p</mi> <mo>&times;</mo> </msubsup> <msub> <mover> <mi>v</mi> <mo>&CenterDot;</mo> </mover> <mi>p</mi> </msub> <mo>=</mo> <msub> <mi>T</mi> <mi>c</mi> </msub> <mo>-</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msubsup> <mi>p</mi> <mi>i</mi> <mo>&times;</mo> </msubsup> <msubsup> <mi>A</mi> <mi>eu</mi> <mi>T</mi> </msubsup> <msub> <mi>F</mi> <mi>si</mi> </msub> <mo>+</mo> <msub> <mi>T</mi> <mi>d</mi> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow></math>]]></maths><maths num="0004"><![CDATA[<math> <mrow> <mfenced open='{' close='-'> <mtable> <mtr> <mtd> <msub> <mi>m</mi> <mi>b</mi> </msub> <msub> <mover> <mi>v</mi> <mo>&CenterDot;</mo> </mover> <mi>b</mi> </msub> <mo>-</mo> <msubsup> <mi>S</mi> <mi>b</mi> <mo>&times;</mo> </msubsup> <msub> <mover> <mi>&omega;</mi> <mo>&CenterDot;</mo> </mover> <mi>b</mi> </msub> <mo>=</mo> <mo>-</mo> <msubsup> <mi>A</mi> <mi>eb</mi> <mi>T</mi> </msubsup> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msub> <mi>F</mi> <mi>ui</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>I</mi> <mi>b</mi> </msub> <mrow> <msub> <mover> <mi>&omega;</mi> <mo>&CenterDot;</mo> </mover> <mi>b</mi> </msub> <mo>+</mo> <msubsup> <mi>&omega;</mi> <mi>b</mi> <mo>&times;</mo> </msubsup> <msub> <mi>I</mi> <mi>b</mi> </msub> <msub> <mi>&omega;</mi> <mi>b</mi> </msub> <mo>+</mo> <msubsup> <mi>S</mi> <mi>b</mi> <mo>&times;</mo> </msubsup> <msub> <mover> <mi>v</mi> <mo>&CenterDot;</mo> </mover> <mi>b</mi> </msub> <mo>=</mo> <mo>-</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msup> <mrow> <mo>(</mo> <msub> <mi>A</mi> <mi>bd</mi> </msub> <msub> <mi>q</mi> <mi>i</mi> </msub> <mo>+</mo> <msub> <mi>r</mi> <mi>db</mi> </msub> <mo>)</mo> </mrow> <mo>&times;</mo> </msup> <msup> <mrow> <mo>(</mo> <msubsup> <mi>A</mi> <mi>eb</mi> <mi>T</mi> </msubsup> <msub> <mi>F</mi> <mi>ui</mi> </msub> <mo>)</mo> </mrow> </msup> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow></math>]]></maths>式(7)是上平台系统的动力学方程;式(8)是下平台系统的动力学方程;式中m<sub>p</sub>和I<sub>p</sub>分别是上平台系统的质量和惯量;m<sub>b</sub>和I<sub>b</sub>分别是卫星星体(9)的质量和惯量;S<sub>p</sub>和S<sub>b</sub>分别是上平台系统和卫星星体(9)的静矩;F<sub>d</sub>和T<sub>d</sub>分别是执行机构自身产生的扰动力和扰动力矩;T<sub>c</sub>是执行机构产生的有效输出力矩;A<sub>bd</sub>是从隔振平台下平台坐标系到卫星星体坐标系的坐标转换矩阵;将得到的含有单个隔振平台的整星动力学模型扩展成含有多级隔振平台的整星动力学模型;对于含有多级隔振平台的整星动力学模型,是由多个上平台系统动力学方程和一个下平台系统动力学方程组成;首先将式(7)推广成为第j个上平台系统动力学方程,具体如下所示:<maths num="0005"><![CDATA[<math> <mrow> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <msub> <mi>m</mi> <mi>pj</mi> </msub> <msub> <mover> <mi>v</mi> <mo>&CenterDot;</mo> </mover> <mi>pj</mi> </msub> <mo>-</mo> <msubsup> <mi>S</mi> <mi>pj</mi> <mo>&times;</mo> </msubsup> <msub> <mover> <mi>&omega;</mi> <mo>&CenterDot;</mo> </mover> <mi>pj</mi> </msub> <mo>=</mo> <mo>-</mo> <msubsup> <mi>A</mi> <mi>cuj</mi> <mi>T</mi> </msubsup> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msub> <mi>F</mi> <mi>sji</mi> </msub> <mo>+</mo> <msub> <mi>F</mi> <mi>dj</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>I</mi> <mi>pj</mi> </msub> <mrow> <msub> <mover> <mi>&omega;</mi> <mo>&CenterDot;</mo> </mover> <mi>pj</mi> </msub> <mo>+</mo> <msubsup> <mi>&omega;</mi> <mi>pj</mi> <mo>&times;</mo> </msubsup> <msub> <mi>I</mi> <mi>pj</mi> </msub> <msub> <mi>&omega;</mi> <mi>pj</mi> </msub> <mo>+</mo> <msubsup> <mi>S</mi> <mi>pj</mi> <mo>&times;</mo> </msubsup> <msub> <mover> <mi>v</mi> <mo>&CenterDot;</mo> </mover> <mi>pj</mi> </msub> <mo>=</mo> <msub> <mi>T</mi> <mi>jc</mi> </msub> <mo>-</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msubsup> <mi>p</mi> <mi>ji</mi> <mo>&times;</mo> </msubsup> <msubsup> <mi>A</mi> <mi>euj</mi> <mi>T</mi> </msubsup> <msub> <mi>F</mi> <mi>sji</mi> </msub> <mo>+</mo> <msub> <mi>T</mi> <mi>dj</mi> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow></math>]]></maths>其中,各变量的物理含义参看方程(7)中的对应解释;在方程(7)的应用过程中,仅有当执行机构安装在隔振平台上平台的时候,才有执行机构自身产生的扰动力和扰动力矩以及执行机构产生的有效输出力矩项;当光学有效载荷(8)安装在隔振平台(10)上平台(1)的时候,不含有这些扰动力、扰动力矩以及有效输出力矩项;由于加装了多个隔振平台(10),使得每个隔振平台(10)下平台(3)连接点处的力均要作用在卫星星体(9)上,因此式(8)扩展成如下形式:<maths num="0006"><![CDATA[<math> <mrow> <mfenced open='{' close='-' separators=''> <mtable> <mtr> <mtd> <msub> <mi>m</mi> <mi>B</mi> </msub> <msub> <mover> <mi>v</mi> <mo>&CenterDot;</mo> </mover> <mi>b</mi> </msub> <mo>-</mo> <msubsup> <mi>S</mi> <mi>B</mi> <mo>&times;</mo> </msubsup> <msub> <mover> <mi>&omega;</mi> <mo>&CenterDot;</mo> </mover> <mi>b</mi> </msub> <mo>=</mo> <mo>-</mo> <msubsup> <mi>A</mi> <mi>eb</mi> <mi>T</mi> </msubsup> <mrow> <mo>(</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>M</mi> </munderover> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msub> <mi>F</mi> <mi>uji</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>F</mi> <mi>ext</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>I</mi> <mi>B</mi> </msub> <msub> <mover> <mi>&omega;</mi> <mo>&CenterDot;</mo> </mover> <mi>b</mi> </msub> <mo>+</mo> <msubsup> <mi>&omega;</mi> <mi>b</mi> <mo>&times;</mo> </msubsup> <msub> <mi>I</mi> <mi>B</mi> </msub> <msub> <mi>&omega;</mi> <mi>b</mi> </msub> <mo>+</mo> <msubsup> <mi>S</mi> <mi>B</mi> <mo>&times;</mo> </msubsup> <msub> <mover> <mi>v</mi> <mo>&CenterDot;</mo> </mover> <mi>b</mi> </msub> <mo>=</mo> <mo>-</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>M</mi> </munderover> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msup> <mrow> <mo>(</mo> <msub> <mi>A</mi> <mi>bdj</mi> </msub> <msub> <mi>q</mi> <mi>ji</mi> </msub> <mo>+</mo> <msub> <mi>r</mi> <mi>djb</mi> </msub> <mo>)</mo> </mrow> <mo>&times;</mo> </msup> <mrow> <mo>(</mo> <msubsup> <mi>A</mi> <mi>eb</mi> <mi>T</mi> </msubsup> <msub> <mi>F</mi> <mi>uji</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>T</mi> <mi>ext</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow></math>]]></maths>同方程(8)相比,由于引入了M个隔振平台(10),使得卫星星体(9)的质量和惯量增加,分别用m<sub>B</sub>和I<sub>B</sub>表示;星本体系统的静矩也会发生变化,用S<sub>B</sub>表示;其他变量的物理意义参看方程(8)的对应解释;求解式(9)和式(10),可以得到卫星星体角速度ω<sub>b</sub>或任意上平台系统的角速度ω<sub>pj</sub>;各支杆(2)作用力的求解,即式(1)~式(6)对每个隔振平台(10)都适用,并且各支杆(2)作用力的表达式也均相同,只是在使用过程中为便于表达,第j个隔振平台支杆作用力表达式中会加入下标j,以进行区分不同的隔振平台(10);步骤八:在卫星姿态稳定控制中姿态角度φ、θ和ψ为小角度,可以得知卫星星体的姿态角速度等同于卫星星体角速度ω<sub>b</sub>,上平台系统的姿态角速度等同于上平台系统角速度ω<sub>pj</sub>;对姿态角速度进行时间积分可以解得姿态角度;由此可以得知上平台系统的姿态角度和角速度以及卫星星体的姿态角度和角速度;步骤九:再重复步骤二到步骤八,直到卫星姿态稳定;即可绘制光学有效载荷姿态角度和姿态角速度效果图,可通过效果图判定光学有效载荷成像精度和稳定度的提高程度;所述的隔振平台(10)包括上平台(1)、下平台(3)以及连接上平台(1)和下平台(3)的支杆(2);所述的隔振平台(10)安装在飞轮(5)或者控制力矩陀螺(6)卫星执行机构和卫星星体(9)之间,或光学有效载荷(8)和卫星星体(9)之间;当在飞轮(5)和卫星星体(9)之间加入隔振平台(10)的时候,令隔振平台(10)的上平台(1)和飞轮(5)外壳固连,共同组成了一套上平台系统;当在控制力矩陀螺(6)和卫星星体(9)之间使用隔振平(10)台时,需加装一套安装架(7),安装架(7)的一端固连在隔振平台(10)的上平台(1)处,安装架(7)的另一端固连在控制力矩陀螺的框架上,以保证框架沿着轴向转动,此时,控制力矩陀螺(6)、安装架(7)和隔振平台(10)的上平台(1)共同组成了上平台系统;当需要在光学有效载荷(8)和卫星星体(9)加装隔振平台(10)时,光学有效载荷(8)的底座和隔振平台(10)的上平台(1)固连,同样能够组成一套上平台系统;所有隔振平台(10)的下平台(3)都固定连接在卫星星体(9)上,共同组成了下平台系统。
地址 100081 北京市海淀区中关村南大街5号