第07课 - 牛顿-欧拉法

运动学告诉"怎么动"，动力学告诉"需要多大的力"。牛顿-欧拉法通过正向递推速度/加速度、反向递推力/力矩，高效计算关节驱动力矩。

📐 核心公式

💻 2R臂牛顿-欧拉法实现

Python

import numpy as np
class NE2R:
    def __init__(self,L1=1.0,L2=0.8,m1=5.0,m2=3.0,lc1=0.5,lc2=0.4,I1=0.5,I2=0.2,g=9.81):
        self.p={'L1':L1,'L2':L2,'m1':m1,'m2':m2,'lc1':lc1,'lc2':lc2,'I1':I1,'I2':I2,'g':g}
    def inv_dyn(self,q,qd,qdd):
        p=self.p; t1,t2=q; td1,td2=qd; tdd1,tdd2=qdd
        c2,s2=np.cos(t2),np.sin(t2)
        # 惯性+科里奥利+重力
        h=p['m2']*p['L1']*p['lc2']*s2
        M11=p['m1']*p['lc1']**2+p['I1']+p['m2']*(p['L1']**2+p['lc2']**2+2*p['L1']*p['lc2']*c2)+p['I2']
        M12=p['m2']*(p['lc2']**2+p['L1']*p['lc2']*c2)+p['I2']
        M22=p['m2']*p['lc2']**2+p['I2']
        M=np.array([[M11,M12],[M12,M22]])
        C=np.array([[-h*td2,-h*(td1+td2)],[h*td1,0]])
        G=np.array([(p['m1']*p['lc1']+p['m2']*p['L1'])*p['g']*np.cos(t1)+p['m2']*p['lc2']*p['g']*np.cos(t1+t2),
                     p['m2']*p['lc2']*p['g']*np.cos(t1+t2)])
        return M@qdd+C@qd+G
arm=NE2R()
q=np.array([np.pi/4,np.pi/6]); qd=np.array([1.0,0.5]); qdd=np.array([0.5,0.3])
print("牛顿-欧拉法逆动力学"); print("="*55)
tau_grav=arm.inv_dyn(q,np.zeros(2),np.zeros(2))
tau_full=arm.inv_dyn(q,qd,qdd)
print(f"重力力矩: τ=({tau_grav[0]:.2f},{tau_grav[1]:.2f}) N·m")
print(f"全力矩: τ=({tau_full[0]:.2f},{tau_full[1]:.2f}) N·m")
print("\n重力力矩扫描(θ₂=0):")
for d in [0,30,60,90]:
    q_t=np.array([np.radians(d),0]); tau=arm.inv_dyn(q_t,np.zeros(2),np.zeros(2))
    print(f"  θ₁={d:3d}°→τ=({tau[0]:+8.2f},{tau[1]:+8.2f})")

运行结果

牛顿-欧拉法逆动力学 ======================================================= 重力力矩: τ=(41.20,3.05) N·m 全力矩: τ=(44.72,4.71) N·m 重力力矩扫描(θ₂=0): θ₁= 0°→τ=( +65.73, +11.77) θ₁= 30°→τ=( +56.92, +10.19) θ₁= 60°→τ=( +32.86, +5.89) θ₁= 90°→τ=( +0.00, +0.00)

🔬 工程实践要点

动力学实施注意事项

参数辨识：实际机器人的质量、质心、惯量参数难以精确获取，需要通过运动激励和最小二乘辨识
实时计算：逆动力学需要1kHz计算，Newton-Euler递推比Lagrange更高效（O(n) vs O(n³)）
摩擦补偿：关节摩擦（Coulomb + 粘性）是最大的未建模项，需要专门的摩擦模型
柔性关节：谐波减速器的柔性会引入谐振，需要弹性关节模型
重力补偿精度：5%的重力补偿误差可能导致数毫米的位置偏差

📖 延伸阅读

本课在知识体系中的位置

本课程按"运动学→动力学→力控制→感知规划→实战项目"组织，每课都建立在前面知识基础上。理解本课后，你将：

深入理解机械臂系统的核心原理
独立实现和验证关键算法
为后续课程打下坚实基础
将理论应用到实际工程项目

🎯 学习目标回顾

📋 本课小结

关键要点

本课的核心概念和公式已在仿真中验证
Python实现提供了从理论到代码的完整路径
课后练习将帮助你深化理解和应用能力

下一步

完成练习后，继续下一课的学习。本课内容将在后续课程中被反复使用和扩展。

📚 动力学理论基础

牛顿-欧拉递推的物理意义

牛顿-欧拉法的第一步递推（正向）从基座向末端传播运动量：

每个连杆的速度和加速度由前一个连杆的运动加上当前关节的运动确定
科里奥利加速度和离心加速度是运动耦合的体现
重力加速度作为基座加速度加入递推

第二步递推（反向）从末端向基座传播力：

每个连杆的惯性力等于质量乘加速度
关节力矩由力矩平衡方程确定
相邻连杆间的力传递遵循牛顿第三定律

拉格朗日法的优势

虽然拉格朗日法推导复杂，但它有独特优势：

自动获得惯性矩阵M(q)的解析表达式
斜对称性Ṁ-2C可以直接验证
便于设计基于能量的控制器
适合符号推导和参数辨识

关节摩擦模型

实际关节的摩擦力远比简单阻尼模型复杂，常用模型包括：

Coulomb摩擦：F = μ·N·sign(ω)，与速度方向有关
粘性摩擦：F = b·ω，与速度成正比
Stribeck效应：低速时摩擦力下降
LuGre模型：统一的动态摩擦模型

摩擦补偿是高精度运动控制的必要环节。

🏭 工业应用实例

本课知识在实际工程中的应用

本课涉及的技术在以下场景中有重要应用：

计算力矩控制器设计、前馈补偿精度、实时动力学计算、重力补偿标定——核心算法的直接应用
不同应用场景对算法精度、速度和鲁棒性有不同要求
工业实践中通常需要在理想模型和实际约束之间做权衡
从实验室仿真到工业部署需要大量的工程化工作

理解理论与实践的差距，是成为优秀机器人工程师的关键。

❓ 常见问题解答

Q1: 本课的算法在实际机器人上能直接使用吗？

仿真验证的算法通常需要以下适配才能部署到实际机器人：参数标定、实时性优化、异常处理、安全保护。仿真与现实的差距（sim-to-real gap）是机器人领域的重要挑战。

Q2: Python实现的效率够用吗？

Python适合算法验证和原型开发，但生产级控制器通常用C++实现。Python的NumPy/SciPy计算效率约为C++的1/10到1/100，对于1kHz控制频率可能不够。可以使用Cython、Numba或直接C++重写来加速。

Q3: 如何从仿真过渡到实际机器人？

推荐步骤：(1) Python仿真验证算法正确性；(2) 添加传感器噪声和延迟模型测试鲁棒性；(3) 用C++重写核心计算模块；(4) 在低速度下实际测试；(5) 逐步提高速度和负载。安全永远是第一位的。

🐛 调试技巧

本课代码的常见问题和调试方法

数值不稳定：检查除零、sqrt负数、矩阵奇异等情况，添加小的正则化项
IK不收敛：尝试不同初始值，增大阻尼系数，检查目标可达性
仿真发散：减小时间步长，检查控制器增益，添加速度限幅
结果异常：用简单情况验证（如θ=0、θ=π/2等），对比解析解
性能问题：避免在循环中创建数组，预分配内存，向量化计算

📖 术语表

本课关键术语

术语	英文	定义
正运动学	Forward Kinematics	已知关节角求末端位姿
逆运动学	Inverse Kinematics	已知末端位姿求关节角
雅可比矩阵	Jacobian Matrix	关节速度到末端速度的映射
奇异位形	Singularity	末端失去某些方向运动能力的位形
工作空间	Workspace	末端可达的空间区域
自由度	Degrees of Freedom	独立运动变量的数量
齐次变换	Homogeneous Transform	4x4矩阵表示位姿
阻抗控制	Impedance Control	控制力与位移的动态关系
导纳控制	Admittance Control	输入力输出位移修正
轨迹规划	Trajectory Planning	生成平滑的运动时间函数
力封闭	Force Closure	接触力可抵抗任意外力
碰撞检测	Collision Detection	判断几何体是否相交
路径规划	Path Planning	在障碍物间找到安全路径
视觉伺服	Visual Servoing	基于视觉反馈的运动控制
手眼标定	Hand-Eye Calibration	确定相机与机器人坐标系的变换

牛顿-欧拉法

⚙️ 动力学的意义

两步递推