4、连续体质心计算：体积分法、密度均匀与不均匀情况、对称性简化技巧
连续体的质心计算，说白了就是积分。你想想看，离散点系我们能用加权平均搞定，但到了连续体——比如一块机翼蒙皮、一个发动机叶片——质量是连续分布的，那就得用积分来“求和”。
我个人习惯把这个问题拆成三步走：建立坐标系 → 写出密度函数 → 套体积分公式。嗯，就这么简单。但实际项目中，真正让人头疼的不是公式本身，而是积分怎么算、密度怎么处理、以及能不能偷懒用对称性。
4.1 体积分法：质心计算的基本公式
先看最一般的情况。一个连续体占据空间区域 \( V \)，密度函数为 \( \rho(x,y,z) \)，那么它的总质量是：
M = ∫∫∫_V ρ(x,y,z) dV质心坐标 \( (x_c, y_c, z_c) \) 由下式给出：
x_c = (1/M) ∫∫∫_V x · ρ(x,y,z) dV
y_c = (1/M) ∫∫∫_V y · ρ(x,y,z) dV
z_c = (1/M) ∫∫∫_V z · ρ(x,y,z) dV我在项目中遇到过不少新手，一上来就套公式，结果积分区域搞错了。记住：积分区域 V 必须是整个物体占据的空间，别漏了空心部分。
核心要点：体积分法适用于任意形状、任意密度分布的连续体。但计算量往往较大，需要选择合适的坐标系（直角坐标、柱坐标、球坐标）来简化积分。
4.2 密度均匀的情况：简化到几何中心
如果密度是常数 \( \rho_0 \)，那公式就漂亮多了。总质量 \( M = \rho_0 V \)，质心公式变成：
x_c = (1/V) ∫∫∫_V x dV
y_c = (1/V) ∫∫∫_V y dV
z_c = (1/V) ∫∫∫_V z dV说白了，均匀密度物体的质心就是它的几何中心（形心）。这个结论太有用了。我记得有一次做无人机机翼的质心估算，翼型截面是NACA 0012对称翼型，我直接取弦长的25%位置作为初始质心，后续有限元验证误差不到2%。
常见的均匀形状质心位置，我建议你记在脑子里：
| 形状 | 质心位置 |
|---|---|
| 长方体 | 几何中心 |
| 球体 | 球心 |
| 圆柱体 | 轴线中点 |
| 三角形薄板 | 三条中线的交点（重心） |
| 半圆薄板 | 距直径边 4R/(3π) 处 |
我的小技巧：遇到复杂形状但密度均匀时，先拆成几个标准形状的组合，分别算形心，再用加权平均。这比直接积分快得多。
4.3 密度不均匀的情况：实战中的常态
真实飞行器很少有密度均匀的。复合材料蒙皮、泡沫芯材、局部加强片……密度处处不同。这时候必须老老实实处理 \( \rho(x,y,z) \)。
我建议按以下步骤来：
- 建立密度模型：把密度函数写出来。可能是分段常数、线性变化，或者从材料数据表插值。
- 选择积分顺序：先对密度变化最剧烈的方向积分，往往能简化计算。
- 数值积分兜底：如果解析积分搞不定，直接用高斯求积法或蒙特卡洛法。
举个例子，一个变密度杆，长度 L，密度沿轴向线性变化：\( \rho(x) = \rho_0 (1 + \alpha x) \)。那么：
M = ∫_0^L ρ_0 (1 + αx) · A dx = ρ_0 A (L + αL²/2)
x_c = (1/M) ∫_0^L x · ρ_0 (1 + αx) · A dx
= [ρ_0 A (L²/2 + αL³/3)] / [ρ_0 A (L + αL²/2)]
= (L²/2 + αL³/3) / (L + αL²/2)你看，解析解也不复杂。但如果是二维或三维的不均匀密度，我通常会写一小段Python脚本做数值积分。
避坑指南：我曾经在计算火箭燃料储箱的质心时，忽略了燃料晃动导致的密度分布变化，结果质心偏移了5%。后来我养成了习惯——对于液体或颗粒状内容物，一定要考虑工作状态下的密度分布，而不是静止状态。
4.4 对称性简化技巧：能偷懒就偷懒
这是我最喜欢的一节。飞行器设计里，对称性无处不在。利用好对称性，计算量能减少一半以上。
基本规则：
- 面对称：如果物体关于某个平面对称，且密度分布也对称，那么质心一定在这个平面上。
- 轴对称：如果物体关于某条轴线对称，质心一定在这条轴线上。
- 中心对称：如果物体关于某点对称，质心就是这个点。
你想想看，一个对称的机翼，左右半翼的质心一定在对称面上。那我只需要算一半的积分，然后直接镜像过去就行。
具体操作时，我习惯这样做：
- 先找对称面/轴：把坐标系原点放在对称中心上。
- 判断密度是否对称：如果材料分布也对称，直接断言质心在对称面上。
- 只算非对称方向：比如面对称时，只需要计算垂直于对称面方向的坐标。
实战案例：某型无人机机身，关于XZ平面对称，密度分布也对称。我直接设 y_c = 0，只算 x_c 和 z_c。原本需要三重积分，降成了二重积分，计算时间从半小时缩到两分钟。
但要注意——对称性只适用于密度分布也对称的情况。如果一边是铝合金、另一边是钛合金，那对称面就不管用了。嗯，这个坑我踩过。
4.5 知识体系总览
下面这张图是我自己整理的连续体质心计算知识框架，帮你理清思路：
4.6 一个完整的Python示例
最后，给你一个我常用的Python脚本。它用数值积分计算一个密度不均匀的楔形块的质心。代码不长，但很实用：
import numpy as np
from scipy import integrate
def wedge_mass_center(L, W, H, rho_func):
"""
计算楔形块（底边长L，宽W，高H）的质心
rho_func(x,y,z) 是密度函数
"""
def integrand_mass(z, y, x):
return rho_func(x, y, z)
def integrand_x(z, y, x):
return x * rho_func(x, y, z)
def integrand_y(z, y, x):
return y * rho_func(x, y, z)
def integrand_z(z, y, x):
return z * rho_func(x, y, z)
# 楔形区域：x从0到L，y从0到W，z从0到H*(1 - x/L)
M, err = integrate.nquad(integrand_mass,
[[0, L], [0, W], lambda x: [0, H*(1 - x/L)]])
Mx, err = integrate.nquad(integrand_x,
[[0, L], [0, W], lambda x: [0, H*(1 - x/L)]])
My, err = integrate.nquad(integrand_y,
[[0, L], [0, W], lambda x: [0, H*(1 - x/L)]])
Mz, err = integrate.nquad(integrand_z,
[[0, L], [0, W], lambda x: [0, H*(1 - x/L)]])
return (Mx/M, My/M, Mz/M)
# 示例：密度随高度线性增加
rho = lambda x,y,z: 1.0 + 0.5*z
xc, yc, zc = wedge_mass_center(2.0, 1.0, 0.5, rho)
print(f"质心坐标: ({xc:.3f}, {yc:.3f}, {zc:.3f})")这个脚本我用了好几年，稍微改改就能适配各种形状。你拿过去直接用就行。
公众号:蓝海资料掘金营，微信deep3321