飞行器转动惯量叠加

简介

在飞行器总体设计过程中会涉及到转动惯量张量叠加。对于一个固体推进剂飞行器,已知推进剂装药燃尽后的转动惯量张量(纯结构),需要计算燃烧过程中飞行器整体的转动惯量张量。

计算

以飞行器头部为转动惯量张量的参考点,x轴指向尾喷管

燃烧过程中的转动惯量张量为纯结构与推进剂的叠加
I = I 0 + I p {\bf{I}} = {{\bf{I}}0} + {{\bf{I}}{\rm{p}}} I=I0+Ip

而推进剂的转动惯量张量需要运用平行轴定理计算
I p = [ I x x I x y I x z I y x I y y + h 2 m p I y z I z x I z y I z z + h 2 m p ] \mathbf{I}{\mathrm{p}}=\left[\begin{array}{ccc} I{x x} & I_{x y} & I_{x z} \\ I_{y x} & I_{y y}+h^{2} m_{\mathrm{p}} & I_{y z} \\ I_{z x} & I_{z y} & I_{z z}+h^{2} m_{\mathrm{p}} \end{array}\right] Ip=⎣⎡IxxIyxIzxIxyIyy+h2mpIzyIxzIyzIzz+h2mp⎦⎤

式中 m p m_p mp为当前装药质量, h h h为装药局部坐标系的原点与飞行器头部之间的距离, I x x I_xx Ixx, I x y I_xy Ixy...等为装药局部坐标系下所得的转动惯量,定义式为
I x x = ∫ ( y 2 + z 2 ) d m , I x y = − ∫ x y d m {I_{xx}} = \int {\left( {{y^2} + {z^2}} \right){\rm{d}}m,} {\rm{ }}{I_{xy}} = - \int {xy{\rm{d}}m} Ixx=∫(y2+z2)dm,Ixy=−∫xydm

一般地,由于装药的旋转对称性,非对角线元素均为0

质心位置的叠加如下
c = m 0 c 0 + m p c p m 0 + m p = m 0 c 0 + m p ( h + c p , l o c a l ) m 0 + m p c = {{{m_0}{c_0} + {m_{\rm{p}}}{c_{\rm{p}}}} \over {{m_0} + {m_{\rm{p}}}}} = {{{m_0}{c_0} + {m_{\rm{p}}}\left( {h + {c_{{\rm{p,local}}}}} \right)} \over {{m_0} + {m_{\rm{p}}}}} c=m0+mpm0c0+mpcp=m0+mpm0c0+mp(h+cp,local)

相关推荐
斯内科5 分钟前
四胞胎素数:找出‌个位数分别是 1、3、7、9‌,且‌十位及更高位数字完全相同‌的质数,例如 11、13、17、19
算法·质数·素数·四胞胎素数
Hello.Reader36 分钟前
算法基础(十二)——主方法:快速求解常见递归式
算法
小O的算法实验室1 小时前
2026年IEEE TITS,面向按需外卖配送调度的特定问题知识与基于学习元启发式算法,深度解析+性能实测
算法·论文复现·智能算法·智能算法改进
加勒比海带661 小时前
目标检测算法——农林行业数据集汇总附下载链接【Plant】
大数据·图像处理·人工智能·算法·目标检测
洛水水1 小时前
【力扣100题】23. 螺旋矩阵
算法·leetcode·矩阵
影sir1 小时前
不同测试数据下,该如何选择算法
算法·深度优先
潇湘散客2 小时前
CAX软件插件化设计实现牛刀小试
c++·算法·图形学·opengl
速易达网络2 小时前
2026,视觉算法正在经历一场静默革命
算法
WBluuue2 小时前
Codeforces 1094 Div1+2(ABCDE)
c++·算法
TENSORTEC腾视科技2 小时前
腾视科技大模型一体机解决方案:低成本私有化落地,重塑行业智能应用新格局
大数据·人工智能·科技·算法·ai·零售·大模型一体机