要得到整个天体的总质量M，我们只需要把所有薄球壳的质量加起来——也就是从天体中心（r=0）到最外层（r=R，R是天体的最大半径）进行积分，公式就是：
M=∫（从0到R）ρ(r)×4πr²dr
这个积分公式看起来复杂，但只要知道了密度随半径变化的函数ρ(r)，就能通过计算得到总质量。
3.2.2 质能方程与旋转动能的变密度适配
有了变密度下的总质量公式，我们就能很容易地得到质能关系：E=Mc²=∫（从0到R）ρ(r)×4πr²dr ×c²。这个式子告诉我们，变密度天体的总能量，是每个薄球壳能量的累加，更符合天体的实际情况。
再看旋转动能。旋转动能的核心是转动惯量，对于薄球壳微元，转动惯量dI=mr²=ρ(r)×4πr²dr ×r²=ρ(r)×4πr⁴dr。
整个天体的总转动惯量I=∫（从0到R）ρ(r)×4πr⁴dr。
然后代入旋转动能公式Eₖ=1/2Iω²，就能得到变密度下的旋转动能：
Eₖ=1/2×[∫（从0到R）ρ(r)×4πr⁴dr]×ω²
这样计算出来的旋转动能，比理想模型更精准，因为它考虑了不同半径处密度的差异——核心密度大的区域，对转动惯量的贡献也更大。
3.3 变速度模型的积分推导：捕捉天体的动态演化
除了密度，天体的速度也在变化，比如地球的自转速度其实在慢慢变慢（每年大约慢1.8毫秒），太阳的自转速度也随纬度和时间变化。所以我们需要引入变速度模型，让速度成为时间的函数v(t)。
3.3.1 动态力与平动动能
根据牛顿第二定律，动态力F=ma，而加速度a是速度对时间的导数（a=dv/dt），所以动态力可以表示为：
F=m×dv/dt=∫（从0到R）ρ(r)×4πr²dr ×dv/dt
这个式子能反映出，天体受到的动态力，和它的密度分布、速度变化率都相关。
再看平动动能。平动动能是随时间变化的，因为速度v(t)在变。所以变速度下的平动动能的变化量，需要通过积分计算：
ΔEₖ=∫（从t₀到t₁）F×v(t)dt=∫（从t₀到t₁）[∫（从0到R）ρ(r)×4πr²dr ×dv/dt]×v(t)dt
简化后可以得到：
ΔEₖ=1/2×∫（从0到R）ρ(r)×4πr²dr ×[v²(t₁)-v²(t₀)]
其中v(t₀)是初始速度，v(t₁)是t₁时刻的速度。这个公式能精准捕捉天体在不同时间段内平动动能的变化，反映出速度演化对能量的影响。
3.3.2 旋转速度变化的拓展
旋转角速度ω也是时间的函数ω(t)，所以旋转动能随时间的变化可以表示为：
Eₖ(t)=1/2×[∫（从0到R）ρ(r)×4πr⁴dr]×ω²(t)
如果我们知道了角速度随时间的变化规律（比如ω(t)=ω₀e^(-kt)，其中ω₀是初始角速度，k是衰减系数），就能计算出任意时刻天体的旋转动能，完美贴合天体的动态演化过程。
3.4 变密度-变速度耦合模型：核心公式体系汇总
把变密度和变速度结合起来，我们就得到了完整的耦合模型，核心公式总结如下，大家可以参考：
- 总质量：M=∫（从0到R）ρ(r)×4πr²dr
- 质能关系：E=Mc²=∫（从0到R）ρ(r)×4πr²dr ×c²
- 动态力：F=∫（从0到R）ρ(r)×4πr²dr ×dv/dt
- 平动动能变化量：ΔEₖ=1/2×∫（从0到R）ρ(r)×4πr²dr ×[v²(t₁)-v²(t₀)]
- 旋转动能：Eₖ=1/2×[∫（从0到R）ρ(r)×4πr⁴dr]×ω²(t)
这些公式看起来有很多积分，但核心逻辑很简单：把天体拆成无数微元，分别计算每个微元的物理量，再累加起来，就能得到整个天体的精准结果。这就像我们计算一个不规则物体的体积时，把它切成无数个小正方体，再把小正方体的体积加起来，结果会非常精准。
四、天体实证示例：地球与太阳的模型应用
模型好不好用，不能只看理论推导，还得用实际天体数据来验证。我们选择了地球（类地行星）和太阳（恒星）这两个最具代表性的天体，用模型计算它们的质量、旋转动能、万有引力等参数，再和实测数据对比，看看偏差有多大。
4.1 地球场景的实证与量化
4.1.1 地球的核心参数与密度分布函数
首先，我们收集了地球的基础参数（来自公开的天体物理实测数据）：
- 最大径向半径R=6371千米（地球平均半径）
- 核心半径R₁=3480千米（地核的半径，包括内核和外核）
- 核心密度ρ₁=13.1克/立方厘米（地核主要由铁镍合金组成，密度极高）
- 实际总质量M实测=5.97×10²⁴千克
- 赤道自转角速度ω=7.29×10⁻⁵弧度/秒（地球自转一周约24小时）
根据地球的圈层结构，我们把密度分布分成两段，用指数递减函数拟合：
- 地核区（0≤r≤R₁）：ρ₁(r)=13.1×e^(-0.0005r)（r的单位是千米，这个函数表示地核区密度随半径缓慢衰减）
- 地幔-地壳区（R₁<r≤R）：ρ₂(r)=5.5×e^(-0.001r)（地幔和地壳的密度比地核低很多，衰减速度更快）
