3、离散质量系统质心计算：质点系质心公式、加权平均法、多组件系统的质心合成

各位同学，今天我们来聊聊离散质量系统的质心计算。说实话，这是飞行器设计中我最常打交道的内容之一。你想想看，一架飞机或者一枚火箭，本质上就是一堆零部件的集合——发动机、机翼、油箱、航电设备……每个部件都有自己的重量和位置。我们要做的，就是把它们当成一个个质点，算出整个系统的质心。

3.1 质点系质心公式——最根本的出发点

先看最基础的情况。假设有 n 个质点，每个质点的质量是 m_i，位置向量是 r_i。那么整个质点系的质心位置 r_c 就是：

rc = (Σ mi · ri) / (Σ mi)

这个公式看着简单，但我在项目中见过不少人用错。说白了，这就是一个加权平均——质量越大的部件，对质心的「话语权」就越大。拆开来看：

• 分子：每个质点的质量乘以它的位置，然后求和。这叫「质量矩」或「静矩」。
• 分母：所有质点的质量之和，也就是总质量。

嗯，这里要注意：位置向量 r_i 必须相对于同一个坐标系。我见过有人把不同部件的坐标搞混了，结果算出来的质心偏了十万八千里。所以，开工之前先统一坐标系，这是铁律。

3.2 加权平均法——工程中的实用技巧

加权平均法，其实就是把质心公式换了个说法。但在工程实践中，它特别好用。为什么呢？因为我们可以把每个部件的数据整理成表格，然后逐项计算。

我个人习惯用下面的表格来组织数据：

然后质心坐标就是：

xc = 738 / 250 = 2.952 m
yc = 3.0 / 250 = 0.012 m
zc = 133 / 250 = 0.532 m

你看，是不是很直观？每一列加起来，最后一除，完事。我在做某型无人机质心校核时，就是用这个表格，配合 Excel 或者 Python，几分钟就能搞定。

3.3 多组件系统的质心合成——从局部到整体

实际工程中，我们很少直接面对一堆散乱的质点。更多的情况是：先有子组件（比如机翼、机身、尾翼），每个子组件已经算好了自己的质心和质量。然后我们要把这些子组件「合成」成整个飞行器的质心。

这其实就是递归应用质心公式。假设有 k 个子组件，第 j 个子组件的质量为 M_j，质心位置为 R_j，那么：

Rc = (Σ Mj · Rj) / (Σ Mj)

和之前一模一样，只是把「质点」换成了「子组件」。但这里有个坑——子组件的质心位置 R_j 必须相对于同一个全局坐标系。我曾经在项目中遇到过这样的情况：机翼组给的是局部坐标（以机翼根部为原点），机身组给的是全局坐标，结果合成时忘了转换，算出来的质心偏了半个机身长度。嗯，从那以后，我每次做合成前都会先检查一遍坐标基准。

3.4 知识体系总览

为了让你更直观地理解本章的知识结构，我画了一张流程图：

import numpy as np

def centroid_discrete(masses, positions):
    """
    计算离散质量系统的质心
    
    参数:
        masses: 质量列表 (numpy array)
        positions: 位置矩阵 (n x 3)
    
    返回:
        centroid: 质心坐标 (x, y, z)
    """
    total_mass = np.sum(masses)
    mass_moments = np.sum(masses[:, np.newaxis] * positions, axis=0)
    centroid = mass_moments / total_mass
    return centroid

# 示例：三个部件
masses = np.array([120, 85, 45])
positions = np.array([
    [2.5, 0.0, 0.8],
    [4.2, 0.3, -0.2],
    [1.8, -0.5, 1.2]
])

cg = centroid_discrete(masses, positions)
print(f"系统质心: x={cg[0]:.3f} m, y={cg[1]:.3f} m, z={cg[2]:.3f} m")
# 输出: 系统质心: x=2.952 m, y=0.012 m, z=0.532 m