刚体的质量几何

系统对轴的惯性矩

  设PνP_{\nu}Pν点到某个轴uuu的距离等于ρν\rho_{\nu}ρν,称
Ju=∑ν=1Nmνρν2(1)J_{u}=\sum_{\nu=1}^{N} m_{\nu} \rho_{\nu}^{2}{\tag1}Ju=ν=1Nmνρν2(1)

是系统相对于uuu的惯性矩。可见,当系统为刚体时,系统对于 相对刚体固定的轴 的惯性矩是固定的,不随坐标系的选取而改变。

(图1)
坐标表示

  为了进一步量化分析,下面引入坐标系。假定OOO为刚体内的确定点,OxyzOxyzOxyz是固联在刚体上的坐标系。轴uuu经过原点,其与OxOxOx,OyOyOy,OzOzOz轴夹角的余弦分别为α\alphaα,β\betaβ,γ\gammaγ,如图(1)所示。于是
Ju=∑ν=1Nmνρν2=∑ν=1Nmν[(xν2+yν2+zν2)−(xνα+yνβ+zνγ)2]J_{u}=\sum_{\nu=1}^{N} m_{\nu} \rho_{\nu}^{2}=\sum_{\nu=1}^{N} m_{\nu}\left[\left(x_{\nu}^{2}+y_{\nu}^{2}+z_{\nu}^{2}\right)-\left(x_{\nu} \alpha+y_{\nu} \beta+z_{\nu} \gamma\right)^{2}\right]Ju=ν=1Nmνρν2=ν=1Nmν[(xν2+yν2+zν2)(xνα+yνβ+zνγ)2]

经过整理,得
Ju=Jxα2+Jyβ2+Jzγ2−2Jxyαβ−2Jxzαγ−2Jyzβγ(2)J_{u}=J_{x} \alpha^{2}+J_{y} \beta^{2}+J_{z} \gamma^{2}-2 J_{x y} \alpha \beta-2 J_{x z} \alpha \gamma-2 J_{y z} \beta \gamma{\tag2}Ju=Jxα2+Jyβ2+Jzγ22Jxyαβ2Jxzαγ2Jyzβγ(2)

其中引入了下面的记号
Jx=∑ν=1Nmν(yν2+zν2),Jy=∑ν=1Nmν(zν2+xν2),Jz=∑ν=1Nmν(xν2+yν2)J_{x}=\sum_{\nu=1}^{N} m_{\nu}\left(y_{\nu}^{2}+z_{\nu}^{2}\right), \quad J_{y}=\sum_{\nu=1}^{N} m_{\nu}\left(z_{\nu}^{2}+x_{\nu}^{2}\right), \quad J_{z}=\sum_{\nu=1}^{N} m_{\nu}\left(x_{\nu}^{2}+y_{\nu}^{2}\right)Jx=ν=1Nmν(yν2+zν2),Jy=ν=1Nmν(zν2+xν2),Jz=ν=1Nmν(xν2+yν2)

(3){\tag3}(3)

显然(2)和(3)不依赖于轴uuu的选择。(2)中各量称为轴惯性矩(也称对轴的转动惯量)。JxJ_{x}Jx是对OxOxOx轴的惯性矩,以此类推。(3)中各量称为离心惯性矩(也称惯量积)。轴惯性矩是系统绕相应轴转动时的惯性度量,离心惯性矩可以理解为系统质量非平衡性的度量,它们描述系统相对坐标平面质量分布的非对称性。
  轴惯性矩和离心惯性矩对于不同的点OOO是不同的。当坐标系绕同一点OOO变化时,它们也会改变。

惯量椭球

  惯量椭球是表征刚体相对某个点质量分布的物理量。
  假定OOO是刚体内的确定点。对任意过点OOO的任意轴uuu,在其轴上都有点NNN,使得ONONON的距离与对uuu的惯性矩的平方根成反比,即
ON=1/Ju(4)O N=1 / \sqrt{J_{u}}{\tag4}ON=1/Ju(4)

改变uuu的方向,则对应的ONONONJu\sqrt{J_{u}}Ju随之改变。由各个方向上NNN点在空间中形成会形成封闭曲面。这一封闭曲面与刚体相固联的,其形状与坐标系的选择无关。

坐标表示

  为了进一步分析该曲面,下面引入坐标系。和上文一样,假定OOO为刚体内的确定点,OxyzOxyzOxyz是固联在刚体上的坐标系。轴uuu经过原点,其与OxOxOx,OyOyOy,OzOzOz轴夹角的余弦分别为α\alphaα,β\betaβ,γ\gammaγ。令NNN的几何坐标(x,y,z)(x,y,z)(x,y,z),该坐标满足的方程即是上述的曲面方程。根据式(4),易得下式成立
α=Jux,β=Juyγ=Juz(5)\begin{array}{c} \alpha=\sqrt{J_{u}} x, \quad \beta=\sqrt{J_{u}} y \\ \gamma=\sqrt{J_{u}} z \end{array}{\tag5}α=Jux,β=Juyγ=Juz(5)

代入式(2),得到
Jxx2+Jyy2+Jzz2−2Jxyxy−2Jxzxz−2Jyzyz=1(6)J_{x} x^{2}+J_{y} y^{2}+J_{z} z^{2}-2 J_{x y} x y-2 J_{x z} x z-2 J_{y z} y z=1{\tag6}Jxx2+Jyy2+Jzz22Jxyxy2Jxzxz2Jyzyz=1(6)

  这就是曲面方程在OxyzOxyzOxyz坐标系下的数学表达式,根据解析几何知识判断是一个椭球,这就是惯性椭球。根据前文所述,该曲面(椭球)固联于刚体,其形状与坐标系选择无关。根据解析几何知识,椭球自身有三个方向的主轴。因此,当所选择另一个坐标系Ox′y′z′Ox^{\prime}y^{\prime}z^{\prime}Oxyz与三个方向的主轴重合时,椭球方程可写为
Jx′x′2+Jy′y′2+Jz′z′2=1(7)J_{x^{\prime}} x^{\prime2}+J_{y^{\prime}} y^{\prime2}+J_{z^{\prime}} z^{\prime2}=1{\tag7}Jxx2+Jyy2+Jzz2=1(7)

其中Jx′J_{x^{\prime}}Jx,Jy′J_{y^{\prime}}Jy,Jz′J_{z^{\prime}}Jz被称为主轴惯性矩。如果它们各不相同,则主轴是唯一的。如果椭球是旋转椭球,则主轴选为旋转轴以及赤道平面任意两个垂直的轴。如果三者相等,则所有过OOO点的轴都是惯性主轴。
  惯量椭球形象化地表示出刚体对过OOO点的所有惯量矩分布状况。当有了惯性椭球,对任意轴uuu的惯性矩都等于1/ON21/ON^21/ON2ONONON是连接OOO与椭球交点的线段。相对椭球最短轴的惯性矩最大,相对椭球最长轴的惯性矩最小。

在 Abaqus 中为刚体部件定义质量属性是进行动力学分析的重要步骤,尤其是在涉及刚体运动(如平动和转动)的仿真中。刚体质量属性包括质量(Mass)和转动惯量(Moment of Inertia),它们需要在刚体的参考点(Reference Point)上进行定义。 ### 定义刚体质量属性的条件 当刚体的 **三个平动自由度** 没有被完全约束时,必须为其参考点定义质量;同样,若 **三个转动自由度** 未被完全约束,则必须定义转动惯量。这种定义方式确保在动力学分析中,刚体能够按照牛顿运动定律进行响应[^4]。 ### 设置质量属性的方法 在 Abaqus/CAE 中,可以通过以下步骤为刚体设置质量属性: 1. **创建参考点(Reference Point)** 在 `Model` 下选择 `Reference Point`,输入坐标或选择模型中的某一点作为刚体的参考点。 2. **将部件定义为刚体** 在 `Property` 模块中,选择部件并将其材料行为定义为刚体。可以通过 `Special > Rigid Body` 来指定参考点。 3. **添加质量与转动惯量** 在 `Property` 模块中,选择 `Special > Inertia > Create`,然后选择之前定义的参考点。此时可以输入: - **质量(Mass)** - **转动惯量矩阵(Moment of Inertia)**:包括绕三个坐标轴的转动惯量值(Ixx, Iyy, Izz)以及交叉项(Ixy, Iyz, Izx) ```python # 示例:在 inp 文件中定义刚体质量属性 *Rigid Body, ref point=2, mass=5.0, rotary inertia=10.0 ``` 上述关键字表示参考点为节点编号 2,质量为 5.0,转动惯量为 10.0(适用于绕各轴的惯性矩相同的情况)。 ### 注意事项 - 如果刚体的运动自由度被完全约束,则不需要定义质量或转动惯量。 - 转动惯量的数值应根据刚体的实际几何形状和质量分布计算得出,例如对于一个绕质心旋转的立方体,其绕某轴的转动惯量为 $ I = \frac{1}{6} m a^2 $,其中 $ a $ 是边长。 - 在显式动力学分析(如 ABAQUS/Explicit)中,质量属性的准确性对模拟结果的稳定性与物理合理性至关重要[^1]。 ---
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值