4. 传动系统动态响应建模方法
传动系统由伺服电机、联轴器、减速器、丝杠、齿轮副、负载等多个惯性与弹性环节串联而成，其动态响应既包含各轴系的扭转振动，也包含传动链的弹性变形与间隙非线性。要准确预测系统在指令输入或扰动下的响应，必须建立能够反映这些物理特征的多自由度动力学模型。本节围绕多自由度集中参数建模、拉格朗日方程与牛顿-欧拉方法、状态空间表达、传递函数与频响函数推导，以及刚体模态与弹性模态的区分展开。
一、多自由度集中参数建模
集中参数建模的核心思想是把连续分布的轴、丝杠等弹性体离散为若干具有集中惯量和集中刚度的节点。以典型的“电机—联轴器—减速器—丝杠—工作台”传动链为例，可将电机转子、减速器输入轴、减速器输出轴、丝杠、工作台分别视为独立惯量，相邻惯量之间用扭转弹簧和阻尼器连接。设第 i 个惯量为 Ji，第 i 与第 i+1 个惯量之间的扭转刚度为 ki、阻尼为 ci，则系统具有 n 个自由度，广义坐标取各惯量的角位移 θi。
建模时需注意三点：一是惯量与刚度的折算，减速器两侧的惯量、刚度应按传动比平方折算到同一侧；二是丝杠的轴向刚度与扭转刚度需分别考虑，工作台质量折算为丝杠的等效转动惯量；三是阻尼通常难以精确获得，工程上常按等效粘性阻尼或实验辨识确定。集中参数模型的自由度数量应根据关心的频段选择，一般要求最高分析频率不超过模型最低弹性模态频率的 2～3 倍。
二、拉格朗日方程与牛顿-欧拉方法
拉格朗日方程从能量角度建立动力学方程，适合多自由度、约束复杂的系统。定义拉格朗日函数 L = T − V，其中 T 为系统动能，V 为势能，则第 i 个广义坐标的方程为：
d/dt(∂L/∂θ̇_i) − ∂L/∂θ_i = Q_i其中 Qi 为非保守力对应的广义力，包括电机驱动力矩、负载阻力矩和阻尼耗散力。对于前述传动链，动能 T = ½ΣJiθ̇i²，势能 V = ½Σki(θi+1−θi)²，代入后即可得到 n 个耦合的二阶微分方程，写成矩阵形式为：
J θ̈ + C θ̇ + K θ = T牛顿-欧拉方法则对每个惯量单独列写力矩平衡方程，再通过相邻惯量间的弹性力矩建立联系。两种方法结果一致，但拉格朗日方程更适合系统化推导，牛顿-欧拉方法更直观、便于逐段校核。工程实践中常先用牛顿-欧拉方法建立各段方程，再用拉格朗日方程统一整理为矩阵形式。
三、状态空间表达
将二阶矩阵方程 Jθ̈ + Cθ̇ + Kθ = T 转换为一阶状态空间形式，取状态向量 x = [θT, θ̇T]T，则：
ẋ = [ 0 I ] x + [ 0 ] u
[ −J⁻¹K −J⁻¹C ] [ J⁻¹B ]输出方程根据关注量选取，例如测量电机端角速度时 y = [0, 1, 0, 0, …]x。状态空间模型的优势在于可直接用于现代控制算法设计，如 LQR、状态观测器和模型预测控制，也便于在 Simulink 等环境中搭建仿真。对于 4 自由度系统，状态矩阵为 8×8，特征值即系统的复模态频率与阻尼比。
四、传递函数与频响函数推导
对状态空间模型取拉普拉斯变换，令初始条件为零，可得传递函数矩阵 G(s) = C(sI − A)−1B + D。以电机力矩到工作台位移的传递函数为例，其分母为系统的特征多项式 det(s2J + sC + K)，分子则反映输入输出之间的耦合路径。将 s = jω 代入即得频响函数 H(jω)，其幅频特性曲线上的峰值对应系统各阶固有频率，相位在共振区发生 −180° 左右的快速变化。
工程上常用扫频激励或锤击法获取实测频响，与模型频响对比以修正刚度、阻尼参数。若模型预测的第一阶弹性模态频率与实测偏差超过 10%，应优先检查联轴器刚度、丝杠轴向刚度和工作台折算惯量的取值。
五、刚体模态与弹性模态的区分
当系统不受外部约束、仅存在内部弹性连接时，整体旋转或平移对应的模态频率为零，称为刚体模态。刚体模态下各惯量之间没有相对变形，传动链像一个整体同步运动。弹性模态则对应各惯量之间的相对扭转或轴向变形，频率大于零，是振动控制的主要对象。
区分方法有两种：一是观察特征值，零特征值对应刚体模态；二是观察振型，刚体模态中各节点位移同向且比例一致，弹性模态中则出现节点或反相。实际伺服系统中，电机端通常有位置闭环，刚体模态被控制器约束，但弹性模态仍会激发机械谐振。因此，在建模时若只关心谐振抑制，可先约束刚体自由度，仅保留弹性模态；若需分析整定时间与跟随误差，则必须保留刚体模态。
综上，多自由度集中参数建模为传动系统动态响应分析提供了统一框架。拉格朗日方程与牛顿-欧拉方法分别从能量和力平衡角度导出动力学方程，状态空间表达便于控制设计，传递函数与频响函数则连接了模型与实验。正确区分刚体模态与弹性模态，是后续振动抑制与控制器参数整定的前提。