Simscape Multibody三维物理仿真：从滑块单摆掌握关节与坐标系设计

做Simscape仿真这几年，我最大的感触是：凡是跟“三维”沾边的物理建模，最难的往往不是软件操作，而是你脑子里有没有一套清晰的物理模型。前两篇我们聊了基础刚体和坐标系的搭建逻辑，这一篇我挑了个特别典型的系统——滑块单摆，来讲透三维物理仿真中关节、坐标系、自由度和求解器之间怎么配合。滑块单摆看着简单，实际上是一个被低估的教学模型。它同时包含了平动（滑块）和转动（摆），而且在三维环境下，摆还多了空间摆动的自由度，很多基础的仿真概念都能在这个模型上一一验证。

这篇东西适合两类人：一是刚接触Simscape Multibody、想从二维思路转向三维思路的初学者；二是已经能搭出简单模型，但对关节选型、坐标系设计、求解稳定性这些“看不见的细节”还没有系统认识的用户。我尽量把每一步操作背后的物理依据和逻辑讲清楚，你不仅能照着搭，还能知道为什么要这么搭。

1. 为什么拿滑块单摆当三维仿真的“练功桩”

1.1 平面里的滑块摆，其实没你想的那么无聊

先回顾一下经典力学里的滑块单摆：一个质量为M的滑块放在水平轨道上，一根长度为L的轻杆从滑块底部或侧面的一个铰点垂下，杆末端固定一个质量为m的摆锤。摆一旦偏离竖直位置开始晃动，滑块就会在轨道上来回滑动。这个系统在经典力学里非常有分量，因为它是平动与转动强耦合的典型例子。

如果在二维平面里研究它，我们只需要两个坐标：滑块的水平位移x，以及摆杆与竖直方向的夹角θ。两个自由度，两条运动方程，拉格朗日方程一列，马上就能求解。但问题在于，二维模型里我们默认了一件事：摆只能在某个平面内运动，轨道也只在这个平面内存在。现实中所有的物理装配都有厚度、有宽度、有安装高度差。滑块在导轨上滑行时，导轨的平面度和垂向刚度会影响运动特性；摆杆如果稍有初始侧向偏移，就不再是平面摆动，而是空间摆动。这些效应，二维模型全部无法表达。

1.2 三维环境带来的是自由度质变，而不是视觉升级

很多人以为“三维建模”就是把二维模型加个厚度、上个色，让视觉得到提升。这是最大的误解。三维环境真正的改变在于自由度的数量从“被约束”变成“被选择”。

同样是单摆，二维里只能用Revolute Joint（旋转关节），转轴被限制成一条固定的轴。而在三维空间中，你可以选择三种关节来连接摆杆：

同一套物理结构，仅仅因为关节类型不同，系统呈现出来的动力学行为就完全不同。三维滑块单摆的精髓，就在于你可以自由选择要保留多少自由度，然后用一套通用的建模方法把这些自由度串起来。这是Simscape Multibody相比二维Simscape（比如Simscape Foundation库里的Translational/Rotational元件）最核心的差异。

1.3 本系列第三篇的定位：从“能搭”进化到“会搭”

前两篇我们把Simscape Multibody的界面、基本刚体创建、Rigid Transform坐标变换、以及单个刚体的仿真流程走了一遍。这一篇开始，我们要面对真正的系统级装配——多个刚体之间通过关节约束互相连接，形成一个完整的机构。

选择滑块单摆自作为第三个项目，还有一个现实原因：它足够简单，但结构上完整包含了三维机械仿真中几乎所有的核心要素：

• 地面基座的固定与约束
• 滑块沿单一方向的平动自由度约束
• 悬挂点在滑块上的坐标偏移
• 摆杆的空间转动自由度
• 初始条件的注入方式
• 多刚体之间的动量和能量传递

你可以把滑块单摆看成三维机械仿真的“最小完备系统”。后面无论你要做机械臂、四足机器人还是悬架系统，本质上都是在这个基础上增加更多刚体、更多关节、更多约束。现在把底层逻辑啃透，后面搭什么都能快很多。

2. 关节拓扑与坐标系设计：模型搭不起来，八成是这一步错了

2.1 先画拓扑图，再打开Simulink

我在前面反复强调过一句话：在Simscape Multibody里，模型不是画出来的，是先想出来的。很多新手一上来就拖模块，拖完发现装配报错，再回去检查，折腾半天也不知道哪个环节出了问题。正确的做法是先拿一张纸把拓扑结构画出来。

滑块单摆的标准拓扑链路是：

World Frame → Rigid Transform → Ground（地面基座） ↓ Prismatic Joint（棱柱关节） ↓ Slider（滑块） ↓ Rigid Transform（悬挂点偏移） ↓ Spherical Joint（球关节） ↓ Pendulum Rod（摆杆组件）

这条链路里每个模块都有不可替代的作用。World Frame是全局坐标系的原点。Rigid Transform是“空间变换器”，用来定义两个连接坐标系之间的固定相对位置。Ground刚体提供物理质量属性和可视化外形。关节则负责定义相邻两个刚体之间的相对运动能力。

2.2 关节选型：为什么是棱柱关节+球关节

滑块单摆的两个关节选型，分别对应两种不同的运动约束。

滑块沿直线导轨运动，在机械上对应的是棱柱约束。Simscape Multibody里的Prismatic Joint正好提供这种运动能力——允许两个相邻坐标系沿某一轴方向相对平移，同时约束其他五个自由度（另外两个平移、三个转动）。通过设置关节的参数，你可以指定平移轴的方向在模型坐标系中的朝向。

摆杆连接滑块的方式，如果完全按三维来做，选Spherical Joint最合适。球关节允许三个转动自由度全部保留，摆杆可以朝任何方向晃动。这个选择的意义在于：它让系统的动态行为更丰富。如果你设定初始摆角只在XZ平面内，那么系统退化成一个平面摆；如果你给一个侧向初始速度，摆杆就会走出三维空间中的菊花状轨迹，滑块也会在导轨上来回滑动，整个系统呈现出一个二维模型永远给不了的耦合运动。

注意：球关节在一个模型里是被当作3个独立的转动自由度处理的，这意味着它的方程在数值求解上比旋转关节要“硬”一些。如果你的摆杆是刚性的实心杆，倒是没什么问题；但如果后续要做柔性体，球关节处的约束方程就会变得很复杂。

2.3 坐标系与刚体变换：悬挂点的偏移必须想清楚

这是我最想强调的一点。在Simscape Multibody里，每个刚体、每个关节都有自己独立的坐标系。刚体的本体系（Body Frame）默认位于质心，关节的坐标系则定义了约束轴的原点和方向。

滑块单摆模型中最容易犯错的地方在滑块与摆杆之间的悬挂点位置。滑块实体本身有几何尺寸（比如长100mm、宽60mm、高40mm），它的Body Frame在质心。但摆杆并不应该从质心垂下，而是应该从滑块底面上某个铰接点垂下。这个物理事实体现在模型上，就是需要在滑块与球关节之间插入一个Rigid Transform，用来描述“从滑块质心到球关节铰点”的平移变换。

比如滑块沿X轴放置，质心位于坐标系原点，铰点在滑块底部前端，位置为(40, 0, -20)mm——也就是沿X方向偏40mm、沿Z方向向下20mm。这个偏移如果不加，球关节会直接位于滑块质心。模型也许也能跑起来，但物理意义就完全错了：摆杆变成了从滑块中间穿出来的一根杆，动态响应自然对不上。

在设计Rigid Transform时，问自己三个问题：

1. 这个变换描述的是哪两个坐标系之间的关系？
2. 偏移量的大小和方向是否符合实际装配关系？
3. 变换是否引入了旋转？如果有，旋转轴是否符合物理约束？

这三个问题在滑块单摆模型里只有第一个变换（World到ground base的定位）和第二个变换（滑块质心到悬挂点）需要回答。等你开始搭六轴机械臂时，追问这三个问题的价值就会完全体现出来。

2.4 初始条件设定：直接决定你看到的是“平面摆”还是“空间摆”

在Simulink里，关节模块的配置界面中有Initial Conditions选项卡。在这里可以给滑块设置初始位置和速度，给摆杆设置初始枢轴角度和角速度。

对于滑块单摆，最经典的初始条件组合是：

• 滑块初始位置：x=0，初始速度：0
• 摆杆初始角度：绕Y轴偏转30度（即在XZ平面内偏离竖直方向），初始角速度：0

这种情况下，系统退化为二维平面摆，滑块会沿X轴往复运动，摆杆只在XZ平面内摆动。这个场景适合验证模型的基本动力学是否正确。

如果你想看真正的三维效果，给摆杆绕X轴也加一个角速度分量，比如初始角速度为(0.5, 0, 0) rad/s。这时候摆杆会绕悬挂点做锥形摆动，滑块的运动也会变得更有趣——不仅有跟随摆杆的往复运动，还会因摆杆的空间指向变化产生微妙的浮动轨迹。

我建议你两种初始条件都试一遍。先在平面模式验证物理正确性，再切换到三维模式观察空间耦合效应。这也是调模型的标准思路：从简单到复杂，逐步增加自由度，每一步都确认无误后再进入下一步。

3. Simscape Multibody建模实操：从零开始搭一个能跑的滑块摆

下面进入实操环节。我会用Simscape Multibody的标准模块库，一步步把滑块单摆搭出来。假定你已经会打开Simulink，会新建Simscape Multibody模型，前两篇讲过的基础操作我就不重复了。

3.1 模型基本框架

在MATLAB命令行输入sm_new，会自动创建一个新的Simscape Multibody模型模板，包含Simulink-PS Converter、PS-Simulink Converter、Solver Configuration等基础配置模块。这套骨架是每个多体模型都需要的：Solver Configuration告诉求解器模型里有多少状态并用什么方式求解，两个Converter负责Simulink信号域与物理域之间的信号转换。

建议在搭建之前，先把模型里的模块重命名。这一点很多人忽略，但等模型复杂到一定程度，模块名字就是你看图的导航。比如把Solid模块直接命名为Slider_Solid，把Joint命名为Prismatic_Joint_Slider，一眼就能看出每个模块的物理含义。

3.2 地面基座部分

从Simscape > Multibody > Bodies库中拖入一个Solid模块，作为地面基座。双击打开，在Geometry选项卡里选择形状为Box，输入尺寸。我这里用一个宽300mm、深100mm、厚50mm的长方体作为轨道底座。在Visualization选项卡里设置颜色，让轨道和滑块区分开——我习惯用浅灰色代表固定底座，深蓝色代表运动部件。

在Solid模块之前，需要有坐标系连接。我通常的做法是：

• 拖入World Frame模块
• 在World Frame和Ground Solid之间插入Rigid Transform，让地面基座的夹角落在全局坐标系原点
• 如果地面基座底部跟世界坐标系平面对齐，就不需要额外的旋转偏移，只平移(0, 0, -25)mm让顶面正好落在Z=0平面上

这里要说明一下为什么用Rigid Transform而不是直接把Solid连在World Frame上。Solid模块的输出端口是它自身的坐标系原点。直接连World Frame意味着Solid的自身坐标系与全局坐标系完全重合。如果Solid的几何中心正好在原点，这没问题；但更多情况下，你需要先让Solid偏移到合理位置，再让下游关节连接在正确的位置上。Rigid Transform就是干这个的。

3.3 滑块与棱柱关节

滑块实体同样是Solid模块，尺寸设为80mm×60mm×40mm，材料选铝。材料在Visualization选项卡的Physical Properties里会自动生成质量、密度和惯量矩阵。如果你的模型对精度要求高，建议手动输入质量属性，而不要让软件按几何体自动计算——后面我会专门讲这个坑。

Slider Solid之前放一个Prismatic Joint。双击关节模块，在Parameters选项卡里可以设置：

• 平移轴方向：默认是X轴，我保持默认，让滑块沿X方向滑动
• 平移范围[lower, upper]：比如限制在[-100, 100]mm，模拟轨道两端限位
• 阻尼系数：如果不想让滑块无限滑动，可以设一个较小的值比如0.1 N·s/m
• 摩擦模型：先不启用，保持理想光滑状态，方便后续验证动量守恒

关节的连接顺序是：一个端口连接Ground侧的坐标系，另一个端口连接滑块Solid。注意Simscape Multibody关节的两个端口是有物理方向性的，F端口接“基座侧”，M端口接“运动侧”。以Prismatic Joint为例，F侧的坐标系作为参考，M侧坐标系沿轴方向运动。这个方向不能搞反，否则滑块的运动方向会跟你的预期相反。

3.4 悬挂点与球关节

滑块Solid之后面的Rigid Transform，这一步就是前面反复强调的悬挂点偏移。

变换参数设置如下：

• 平移：X方向40mm（让悬挂点位于滑块的前端面附近）、Y方向0mm、Z方向-20mm（从滑块质心降落到底面）
• 旋转：保持单位矩阵，方向不旋转

这个Rigid Transform的物理含义是：从滑块质心坐标系出发，沿滑块本体的X方向走40mm，再沿-Z方向走20mm，到达铰点位置。

接着放Spherical Joint。球关节默认不限制转角范围，三个转动自由度全部自由。在这个关节模块里，Initial Conditions选项卡可以设置初始角位置和角速度。角位置的单位是deg还是rad取决于你的模型单位设置，一般默认是deg，但我建议在模型配置里把单位改成reflected的SI制，统一用rad，减少换算错误。

3.5 摆杆与摆锤

摆杆我用Solid模块做一个圆柱体，直径8mm，长度L=300mm，材料用钢。摆杆的一端连接球关节的M端，另一端连接摆锤。注意，Solid模块自身的坐标系位于几何体的质心，所以从球关节坐标系到摆杆质心之间要有一个Rigid Transform，平移距离是杆长的一半。

完整链路是：

Spherical Joint的M端 → Rigid Transform（沿杆方向偏移L/2） → Pendulum_Link_Solid（质心） → Rigid Transform（沿杆方向偏移L/2） → Pendulum_Bob_Solid（摆锤）

摆锤用一个直径40mm的球体，材质用黄铜或者自定义密度，让它比摆杆重。之所以在杆末端单独加一个摆锤Solid，而不是把杆做成一个带有大端头的整体，是为了方便后面调参——你想改变摆锤质量，直接改Solid参数就行，不需要动整个杆的几何尺寸。

3.6 初版参数汇总

以下是我推荐的一组初值，照着设置就能跑出标准结果：

这组参数下，系统算出来是无阻尼的，能量应该守恒。你可以先用这组数据跑通，再改参数研究不同摆锤质量对滑块运动的影响。

4. 仿真求解设置：让求解器别成为你的天敌

4.1 求解器类型怎么选

Simscape Multibody的模型是微分代数方程（DAE）系统。它同时包含微分方程（刚体的运动学）和代数约束方程（关节的约束）。所以求解器的选择要针对DAE系统的特性。

对于滑块单摆这种无摩擦、刚体连接、约束方程不退缩的模型，MATLAB默认的ode45可以跑通，但容易出现“求解成功后其实精度不够”的隐患。我个人的经验是：优先考虑ode23t或ode15s，这两种求解器对DAE系统的容错性更好，特别是模型存在闭环约束或刚体接触时，它们能稳定处理，不会因为一个误差尖峰直接报错中断。

在实际仿真里，我建议这样配置：

• 求解器类型：Variable-step（变步长）
• 求解器：ode23t
• Max step size：1e-3 s（默认是auto，但建议手动限制，防止运动过快时跳步过多）
• Relative tolerance：1e-4
• Absolute tolerance：1e-6

原理上，变步长求解器会根据当前动态剧烈程度自动调整步长。滑块单摆系统在刚释放的那一瞬，摆杆角速度从0快速上升，变化率很大，如果步长太大，误差会迅速扩散。限制Max step size可以保证在最剧烈的运动段，求解器也不会跳超过一个足够小的时间步。

4.2 刚体惯量设置：一个容易被忽略但影响巨大的细节

前面提到过，Solid模块默认按几何体形状自动计算惯量。这在纯仿真验证里够用，但有两个潜在问题：

第一，如果你的几何体形状不是标准的长方体、圆柱体或球体，而是从CAD里导入的复杂网格，自动计算出的质量属性可能不准确。

第二，自动计算的质量属性基于你选择的材料密度。如果你改了材料或者忘记了材料设置，惯量计算就会跑偏，而关节约束方程的动力学响应极度依赖惯量矩阵的正确性。

所以我给一个通用建议：对于起关键作用的运动部件，打开Solid模块的Inertia选项卡，手动输入质量值、惯量矩阵。惯量矩阵怎么算？对于长方体滑块，连绕X轴的转动惯量是：

Jx = (1/12) * M * (H^2 + D^2)

其中H为高度，D为深度，M为质量。其他轴的公式教科书上都有，不再赘述。关键是你要知道Simscape Multibody里的惯量矩阵是在Solid本体系下定义的，坐标系原点在质心，所以使用公式计算时也要以质心为参考点。

4.3 可视化与数据记录配置

在Simscape Multibody模型的Simulation选项卡里，Mechanics Explorer会自动打开并给出三维动画。默认情况下，Mechanics Explorer会显示所有Solid外形和坐标系图标。建议把坐标系图标关掉，只保留Solid外形，看起来更干净。

为了记录仿真数据，我习惯用Simulink的Signal Logging或者Simscape的Local Solver配置。具体操作是：

在Solver Configuration模块的属性里，勾选“Log all simscape simulation data”，然后在仿真结束后，通过MATLAB命令行查看simlog。比如我想查看滑块位移，可以运行：

logsout = sim('slider_pendulum'); slider_x = logsout.get('Slider');

如果你想要在Simulink Scope里直接看曲线，也可以在Prismatic Joint的输出端口接上Joint Sensor模块，然后经过PS-Simulink Converter转成普通信号，再接到Scope上。Joint Sensor可以同时测量位置、速度、加速度，非常方便。我一般会给滑块和球关节各挂一个Joint Sensor。

5. 仿真结果解读：从能量和动量看系统行为

5.1 动量守恒：滑块为什么会“反向退”

把滑块单摆跑一遍，你会看到这样的现象：摆杆向右摆（沿X正方向）的同时，滑块向左滑（沿X负方向）。这是动量守恒的直接体现。

这个系统的水平方向没有外力（假设轨道无摩擦、空气阻力忽略），因此系统总水平动量恒定。初始静止时，总水平动量为0。那么任意时刻，滑块和摆的水平动量之和必须为0。摆杆向右摆动时，其质心有沿X正方向的速度分量。为了抵消这个动量，滑块必须获得一个沿X负方向的速度，乘积大小正好弥补摆的水平动量。

在Simscape模型里，滑块的速度曲线和摆杆铰点处的水平速度之间应该满足：

M * v_slider + m * v_pendulum_horizontal = 0

其中v_pendulum_horizontal是摆锤质心的水平速度。

你可以把这两个信号都接出来，在MATLAB里做一个线性组合验证。如果结果不恒等于0，说明模型存在额外的动量注入或者数据记录精度不够。这一步是验证模型物理一致性的金标准。

5.2 能量曲线：判断模型和求解器是否“诚实”

一个无阻尼的保守系统，总能量应该恒定。机械能包括：

• 滑块动能：0.5 * M * v_slider²
• 摆杆动能：0.5 * m * v_pendulum² + 0.5 * J_rod * ω²
• 摆锤动能：0.5 * m_bob * v_bob² + 0.5 * J_bob * ω_bob²
• 重力势能：m * g * h_change

在Simscape里，你不需要自己手动计算这些能量。右键点击模型里的Solid模块，选择“Measure”或者使用Simscape的结果查看器，可以直接查看质量属性和能量数据。但如果你想验证总能量是否守恒，建议还是自己在MATLAB后处理里算一遍，结合前面记录的关节位置、速度数据。

如果能量曲线在某个时间段出现锯齿状波动或者单调衰减，排查顺序是：

1. 是否启用了摩擦或阻尼（包括关节默认阻尼）
2. 求解器误差容限是否过松
3. 是否产生了碰撞或约束冲突
4. 初始条件是否有微小速度导致初始动能为0但能量曲线不为常数

5.3 三维轨迹观察：Mechanics Explorer的进阶使用

前面在球关节初始条件里加了一个绕X轴的角速度分量之后，你会在Mechanics Explorer里看到摆锤在空间中走出复杂的三维轨迹。

建议把Mechanics Explorer的Trace选项打开。在可视化配置中，勾选Solid的Trace，可以记录每一个Solid质心的运动轨迹。这样导出后，摆锤部分会显示出一条空间曲线，直观地告诉你系统的三维动力学行为。

三维轨迹的某种形态可以作为系统是否存在混沌的判断依据。滑块单摆在合适的参数下，如果在三维空间释放，运动通常介于规则与不规则之间。改变摆锤质量、摆杆长度或者初始角速度的组合，你可以观察到轨迹形态的显著变化。这种参数敏感性本身就是动力学系统理论里一个值得记录的现象。

6. 踩坑与调试：我实际遇到过的几个问题

6.1 装配报错：刚体变换的坐标系方向误区

我最早搭滑块单摆时，遇到一个非常典型的报错：

Inconsistency in the connectivity of block 'xxx/Rigid Transform'

排查半天，发现问题是Rigid Transform的旋转参数被我不小心设成一个非零旋转矩阵，导致下游关节的坐标系方向和全局坐标系不平行。在Simscape Multibody里，关节的状态（位置、速度）都是相对于相关坐标系定义的。如果坐标系方向不对，关节装配时会产生一个隐含的不连续，求解器就会拒绝继续。

处理办法：检查每个Rigid Transform的旋转参数，确认它是否与物理装配一致。滑块到悬挂点这个变换，应该是一个纯平移变换（旋转矩阵为单位阵）。如果确实需要旋转，一定画出坐标系的旋转轴和方向，别凭感觉输入。

6.2 模型跑不通：初值问题导致数值奇异

还有一种情况是模型一运行就报：

Singularity in the solution of the algebraic constraints

我记得有一次是因为球关节的初始角度设为90度——摆杆水平放置。这时候摆锤的重力方向正好与摆杆垂直，系统的几何约束处于临界状态。虽然是合法位置，但在数值求解上容易碰到奇异点。

处理办法：初始角度不要给得太极端，先从20-30度开始。如果你一定要从水平位置启动，可以把求解器改成ode23t并把Relative tolerance收紧到1e-6。如果还不行，就检查摆杆是否允许穿过地面，必要时在模型中增加一个Solid代表地面去约束摆杆的行程边界。

6.3 动画卡顿与可视化异常

Mechanics Explorer在模型刚体较多时会卡顿。滑块单摆只有两三个Solid，一般不会卡。但如果你进行了大规模的参数化扫描，在循环里反复打开Mechanics Explorer就会拖慢速度。

处理办法：批量仿真时，在sim命令中设置MechanicsExplorerEnable为false，关掉可视化，只记录数据。需要看结果时再单独运行一次带可视化的仿真。具体代码是：

set_param('model', 'SimscapeLogType', 'all'); set_param('model', 'SimMechanicsOpenAnimationOnUpdate', 'off');

6.4 滑块冲出轨道：限位设置与接触建模

棱柱关节设了平移范围[-100, 100]mm之后，如果初始速度足够大，滑块会撞上限位。默认情况下，限位是硬约束，Simscape会直接阻止滑块越界。但如果限位设置得太靠近初始位置，你会发现模型运行一小段就中断，报出约束冲突。

处理办法：如果只是做动力性验证，直接把限位范围给得足够大，或者在初始条件里把速度设小一点。如果确实要研究限位碰撞，就需要在限位处加弹簧阻尼器或者导入接触力模型，不要依赖关节默认的硬限位。

6.5 利用Simscape的日志信息进行快速调错

最后分享一个调试技巧。遇到模型不收敛或报错时，打开Simulink的“Diagnostics Viewer”窗口查看警告信息。Simscape会报告具体哪个模块、哪个物理变量出现了问题。还有，在MATLAB命令行运行：

sim('slider_pendulum', 'SimscapeLogType', 'all')

仿真结束后，用simlog这个变量检查每个关节的状态变量曲线，跟着物理量走一遍，通常能很快定位到异常节点。

我在实际搭建滑块单摆时，踩过最深的一个坑是忽略了球关节的默认阻尼参数。Simscape Multibody里的关节模块，默认阻尼系数不是0，而是有一个很小的数值。这个默认值在大多数情况下无所谓，但如果你在做能量守恒验证，最后会发现总能量莫名其妙凭空消失了一点点。解决方法是把球关节和棱柱关节的内摩擦参数全部显式设为0，再去验证能量守恒。

另外，给模型加弹簧这个扩展方向也值得试试。在滑块端加一个线性弹簧将滑块拉回中心位置，系统就从“自由漂移”变成了“受迫振荡”。加上弹簧后，滑块单摆的系统固有频率会变得更加明确，你可以在MATLAB里做一次FFT频谱分析，观察系统响应里的频率分量，进一步验证建模的准确性。具体的弹簧参数可以在Joint模块的Actuation选项卡里配置内部弹簧，这种方式比外部单独连一个Spring & Damper物理元件更简洁，也不会引入额外的建模复杂度。

滑块单摆作为三维物理仿真的练手项目，最大的价值在于：它把所有关键概念都浓缩在一个小模型里。关节、坐标变换、自由度、求解器、能量验证、可视化——每一个环节你都能在这个模型上获得第一手体验。把这一篇的模型吃透，后面真正做复杂机械系统仿真时，你会省掉很多因为基础概念不牢而付出的试错成本。