【信息科学与工程学】【制造工程】第十六篇 装备制造系统模型 第一章 力学模型 01
装备制造力学系统模型
|
编号 |
类别 |
领域 |
模型配方 |
定理/算法/模型/方法名称 |
定理/算法/模型/方法的逐步思考推理过程及每一个步骤的数学方程式和参数选择/参数优化 |
精度/密度/误差/强度 |
底层规律/理论定理 |
典型应用场景和各类特征 |
变量/常量/参数列表及说明 |
数学特征 |
时序和交互流程的所有细节/分步骤时序情况及数学方程式 |
流动模型和流向方法的数学描述 |
|---|---|---|---|---|---|---|---|---|---|---|---|---|
|
JX-D1-0001 |
基础定律 |
刚体动力学 |
∑F=ma |
牛顿第二定律 |
1. 概念定义:物体动量变化率等于所受合外力。 |
在宏观低速条件下精度极高,是经典力学的基石。误差主要来源于对“质点”和“刚体”的理想化,以及实际中摩擦、空气阻力等未完全建模的力。 |
牛顿运动定律 |
任何宏观物体的机械运动分析基础。如机床运动部件加速度计算、车辆制动距离估算、机器人关节动力学。特征:向量性、瞬时性、独立性。 |
∑F:合外力向量 (N)。 |
代数:向量方程。 |
时序流程: |
常用于描述动量流。系统动量的变化率等于外力的“流动”。方程 dtdp=F描述了动量 p在时间上的流动,外力 F是流动的驱动力。 |
|
JX-D1-0002 |
本构关系 |
材料力学/弹性力学 |
σ=Eϵ |
胡克定律 (一维) |
1. 基本假设:材料处于线弹性、小变形范围。 |
在小应变(通常 ϵ<0.002)范围内精度高。误差来源于材料的非线性、塑性变形、蠕变等。强度由屈服强度 σs和抗拉强度 σb界定。 |
广义胡克定律 |
杆、梁、轴的拉伸/压缩变形分析;弹簧设计;结构静力分析的基础。特征:线性、可逆、各向同性材料在单轴应力下的响应。 |
σ:正应力 (Pa)。 |
代数:线性方程。 |
静态加载流程: |
描述了“应力-应变”能流的本构关系。在弹性范围内,外界做功以应变能形式存储,dVdU=∫σdϵ=21Eϵ2, 能量流是可逆的。 |
|
JX-D1-0003 |
原理/方法 |
多体系统动力学 |
∑i(Fi−miai)⋅δri=0 |
达朗贝尔原理 (动力学普遍方程) |
1. 静力学原理:虚功原理,∑Fi⋅δri=0(对理想约束)。 |
原理精确。应用中的精度取决于对系统自由度、力和约束建模的准确度。是拉格朗日方程和哈密顿原理的基础。 |
牛顿运动定律 + 虚功原理 |
复杂约束系统(如机械臂、车辆悬挂、机构)的动力学建模。特征:将动力学问题转化为瞬时“平衡”问题,消除了理想约束力。 |
Fi:作用在第 i个质点上的主动力。 |
变分法:基于虚位移(等时变分)。 |
建模与求解流程: |
从“力流”和“动量流”的角度统一了主动力、惯性力和约束。原理本身描述了在虚位移空间上,所有“力流”的功率之和为零,是瞬时能量平衡的表述。 |
|
JX-D1-0004 |
运动方程 |
分析力学/多体动力学 |
dtd(∂q˙j∂L)−∂qj∂L=Qj |
拉格朗日方程 (第二类) |
1. 定义拉格朗日函数:L=T−V, 其中 T为系统动能,V为系统势能。 |
在完整约束系统下精确。数值求解精度取决于积分算法(如龙格-库塔法)。建模误差源于对 T, V, Qj的描述是否准确。 |
哈密顿原理、最小作用量原理 |
复杂多自由度系统(如机器人、航天器、柔性机械臂)的动力学建模与控制。特征:标量形式,自动消除理想约束力,易于计算机符号推导。 |
L:拉格朗日函数 (J)。 |
变分法/泛函极值:核心。 |
系统化建模与仿真流程: |
描述了系统在位形空间(由广义坐标张成)中的“轨迹流动”。哈密顿原理表明,真实轨迹是使“作用量”这个泛函取极值的流线。拉格朗日方程是该流线上每一点满足的局部微分条件。 |
|
JX-D1-0005 |
原理 |
分析力学/最优控制 |
δS=δ∫t1t2L,dt=0 |
哈密顿原理 |
1. 定义作用量:S=∫t1t2L(q,q˙,t)dt。 |
是经典力学的一个基本公理,在宏观低速范围内被认为是精确的。为拉格朗日方程、哈密顿方程、诺特定理等提供基础框架。 |
最小作用量原理 |
作为更基本的原理,用于统一推导连续介质力学、电磁场、量子力学的运动定律。在最优控制(如轨迹规划)中,是庞特里亚金极大值原理的类比。特征:整体性、变分性、协变性。 |
S:作用量 (J·s)。 |
泛函分析/变分法:核心是求泛函 S[q(t)]的极值。 |
原理性流程: |
描述了系统在时空中遵循“最小阻力”或“最经济”路径的流动思想。作用量 S是衡量路径“代价”的全局量,真实运动是其极值路径,类似于“最速降线”或“光程最短”(费马原理)。 |
|
JX-D1-0006 |
场方程/守恒律 |
连续介质力学 (流体/固体) |
∂t∂ρ+∇⋅(ρv)=0 |
连续介质力学控制方程 |
1. 质量守恒:对控制体,质量增加率=流入净通量。微分形式即连续性方程:∂t∂ρ+∇⋅(ρv)=0。 |
是物理学基本守恒律,形式精确。应用精度取决于:1)连续介质假设的尺度有效性;2)本构关系(如材料模型)的准确性;3)边界条件;4)数值求解的离散误差。 |
经典力学守恒定律、连续介质假设 |
流体:空气动力学、流体机械、化工流程。 |
ρ:密度 (kg/m³)。 |
向量分析/场论:大量使用梯度 ∇、散度 ∇⋅、物质导数。 |
CFD/FEM求解通用流程: |
这些方程是质量、动量、能量三种物理量在连续介质中“流动”的严格数学描述。 |
|
编号 |
类别 |
领域 |
模型配方 |
定理/算法/模型/方法名称 |
定理/算法/模型/方法的逐步思考推理过程及每一个步骤的数学方程式和参数选择/参数优化 |
精度/误差/强度 |
底层规律/理论定理 |
典型应用场景和各类特征 |
变量/常量/参数列表及说明 |
数学特征 |
时序和交互流程的所有细节/分步骤时序情况及数学方程式 |
流动模型和流向方法的数学描述 |
|---|---|---|---|---|---|---|---|---|---|---|---|---|
|
JX-D1-0007 |
运动学/微分几何 |
机器人学/机构学 |
v=ω×r+vO |
刚体一般运动的瞬时速度/加速度合成(基点法) |
1. 基点选择:在刚体上任选一点 O(基点),其运动已知(vO,aO)。 |
在刚体假设下精确。实际应用误差来源于关节间隙、构件弹性变形。 |
刚体运动学、向量微积分 |
场景:工业机器人末端执行器速度/加速度分析、机械臂奇异性分析、车辆车轮运动学。 |
vP,aP:点P的绝对速度/加速度。 |
向量代数与几何:核心是叉乘运算。 |
计算流程(以机器人末端为例): |
描述了刚体内运动“流”的传递。速度场在刚体内呈线性分布(螺旋场)。从基点到任意点的速度“流动”是通过旋转“算子” ω×实现的。加速度场则是速度场随时间变化与空间变化的对流。 |
|
JX-D1-0008 |
变换矩阵 |
机器人学/空间运动 |
ATB=[ARB01×3ApBorg1] |
齐次变换矩阵 |
1. 位姿描述:坐标系{B}相对于{A}的位姿包括旋转 ARB和平移 ApBorg。 |
数学表示精确。数值计算中的误差来自浮点数运算和旋转矩阵的正交性保持(需定期正交化)。 |
三维欧几里得几何、线性代数(齐次坐标) |
场景:机器人运动学正解与逆解、数控机床坐标变换、计算机图形学。 |
ATB:从{B}到{A}的 4×4齐次变换矩阵。 |
线性代数/矩阵论:核心是矩阵乘法与求逆。 |
机器人正向运动学流程: |
描述了空间点和坐标系“位姿流”的传递。沿着运动链,点的坐标通过连续的矩阵乘法进行变换,如同在变换的“流”中传递。复合变换是“流”的串联,逆变换是“流”的回溯。 |
|
JX-D1-0009 |
平衡方程 |
结构力学/弹性力学 |
∇⋅σ+b=0 |
弹性静力学平衡方程(柯西方程) |
1. 微元体受力:考虑一个无限小的体积微元 dV。 |
在连续介质和静态假设下精确。求解精度取决于边界条件、几何和本构模型的准确性。 |
牛顿第二定律(静力形式)、散度定理 |
场景:任何静态结构(桥梁、建筑、机床床身、压力容器)的应力分析基础。 |
σ:柯西应力张量(二阶,对称)。 |
向量分析/场论:涉及张量的散度。 |
有限元法求解静力问题流程: |
描述了内力(应力)的流动必须与外力(体力)平衡。方程 ∇⋅σ=−b表明,应力场的“源”是负的体力。应力张量的散度衡量了内力的“净流出率”,它必须抵消体内的“源”以达到平衡。 |
|
JX-D1-0010 |
本构模型 |
塑性力学 |
f(σij)≤0, dϵijp=dλ∂σij∂g |
塑性理论(屈服准则与流动法则) |
1. 弹性界限:定义屈服函数 f(σij), 当 f<0为弹性,f=0为屈服。 |
模型是对复杂物理现象的近似。精度取决于材料参数(屈服强度、硬化模量)的校准和硬化模型的选取。 |
不可逆热力学、最大塑性功原理 |
场景:金属成形(冲压、锻造)、结构极限载荷分析、岩土力学。 |
f:屈服函数。 |
优化/变分不等式:塑性流动服从最大塑性耗散原理。 |
增量塑性计算流程(返回映射算法): |
描述了塑性变形“流”的方向和大小。屈服面定义了弹性域的边界。当应力状态试图“流出”弹性域,塑性应变便开始“流动”,其方向由塑性势梯度决定(通常指向屈服面外法向),大小由一致性条件约束,使应力状态被“拉回”或“推移”至屈服面上。这是一种带约束的“流动”。 |
|
JX-D1-0011 |
运动方程 |
转子动力学 |
Mq¨+(C+G)q˙+Kq=F(t) |
对称刚性转子动力学方程(Jeffcott转子简化模型) |
1. 模型简化:将复杂转子简化为一个在无质量弹性轴中央的刚性圆盘(质量为 m),轴两端简支,考虑轴的横向弯曲刚度 k。 |
模型抓住了转子动力学的核心特征(陀螺效应、临界转速)。精度受限于模型的简化(如忽略轴质量、非线性支撑等)。 |
牛顿-欧拉方程、梁的弯曲理论 |
场景:涡轮机械(汽轮机、压缩机)、发电机、机床主轴等的振动分析与稳定性校核。 |
M:质量矩阵(对角或耦合)。 |
常微分方程组:二阶线性ODE,系数可能时变(若Ω变化)。 |
转子系统振动分析流程: |
\lambda_i |
|
JX-D1-0012 |
控制算法 |
运动控制/自动化 |
u(t)=Kpe(t)+Ki∫0te(τ)dτ+Kddtde(t) |
比例-积分-微分控制器(PID) |
1. 误差定义:e(t)=r(t)−y(t), 其中 r(t)为设定值,y(t)为过程变量。 |
依赖于参数整定。性能指标:超调量、调节时间、稳态误差、抗干扰能力。对线性、单输入单输出系统效果良好,对非线性、强耦合系统需改进。 |
经典控制理论(频域法、根轨迹)、误差驱动调节 |
场景:几乎所有工业过程控制(温度、压力、流量)、运动控制(位置、速度)、机器人关节控制。 |
u(t):控制器输出(控制量)。 |
微积分:核心是比例、积分、微分运算。 |
数字PID控制循环时序: |
描述了基于误差信息的“控制作用流”。比例项是即时响应流,大小与误差成正比。积分项是累积误差流,像水池一样积累过去误差,用于消除静差。微分项是误差变化趋势的预测流,像一个阻尼器,抑制过快的变化。这三股“流”汇合后,共同驱动被控对象向设定点“流动”。 |
|
JX-D1-0013 |
滤波算法 |
信号处理/状态估计 |
x^k−Pk−Kkx^k+Pk+=Fkx^k−1+=FkPk−1+FkT+Qk=Pk−HkT(HkPk−HkT+Rk)−1=x^k−+Kk(zk−Hkx^k−)=(I−KkHk)Pk− |
离散卡尔曼滤波器 |
1. 模型假设:线性动态系统:状态方程 xk=Fkxk−1+wk, 观测方程 zk=Hkxk+vk。wk∼N(0,Qk), vk∼N(0,Rk)。 |
在模型准确且噪声为高斯白噪声的假设下,是最优线性无偏估计。误差由 Pk+量化。模型失配或非高斯噪声会降低性能。 |
概率论(贝叶斯定理)、最优估计理论(最小均方误差) |
场景:导航(GPS/INS融合)、目标跟踪、传感器数据融合、故障诊断。 |
x^:状态估计向量(先验 |
概率与统计:核心是高斯分布的条件期望与协方差更新。 |
卡尔曼滤波递归循环: |
描述了信息(状态估计与不确定性)的融合与更新流。预测步将状态估计和不确定性沿着时间流向前推进,同时添加过程噪声(不确定性增加)。更新步将来自传感器的测量信息流与来自模型的预测信息流进行最优融合。卡尔曼增益 Kk是信息融合的阀门,控制着在多大程度上信任测量。协方差矩阵 P描述了估计不确定性的流动与衰减。 |
|
JX-D1-0014 |
强度理论 |
材料力学/失效分析 |
1. σ1≤σb |
经典强度理论(四大强度理论) |
1. 问题:材料在复杂应力状态(σ1,σ2,σ3)下何时失效?需建立与简单拉伸试验(失效应力 σb或 σs)的联系。 |
理论是近似。第三理论保守;第四理论与多数金属试验数据吻合较好,误差一般在±10%以内。选择取决于材料特性(脆性/塑性)和应力状态。 |
材料破坏机理假设、能量原理 |
场景:机械零件(轴、齿轮、压力容器)的强度设计与校核、安全评估。 |
σ1,σ2,σ3:主应力(σ1≥σ2≥σ3)。 |
代数/不等式:核心是比较表达式与许用值。 |
零件强度校核流程: |
描述了导致材料失效的“力学驱动力”的流动与集中。不同的理论假设了不同的“失效驱动力”流动路径: |
|
JX-D1-0015 |
疲劳模型 |
结构耐久性/可靠性 |
σa=σf′(2Nf)b |
应变-寿命方程(Coffin-Manson方程) |
1. 疲劳过程:材料在循环载荷下,微观裂纹萌生并扩展至断裂。 |
适用于低周疲劳(高应变,寿命 <105次)到高周疲劳的过渡区。在高周区(ϵpa很小),退化为应力-寿命模型。试验数据分散性大,需统计分析。 |
能量耗散原理、循环塑性变形累积损伤 |
场景:承受交变载荷的机械零件(发动机曲轴、飞机起落架、涡轮叶片)的寿命预测。 |
ϵa:总应变幅(控制变量)。 |
幂律/对数坐标线性化:方程两边取对数后,log(ϵa)−log(2Nf)呈分段线性关系。 |
基于局部应变法的疲劳寿命分析流程: |
描述了循环塑性变形导致损伤的“能量/微结构演化流”。塑性应变幅 ϵpa是每次循环中不可逆的“损伤流”的度量,其累积导致裂纹萌生。弹性应变幅 ϵea也贡献损伤,尤其在长寿命区。方程本质是损伤累积速率与应变幅之间的本构关系。总损伤流(导致失效)是弹性与塑性损伤流的叠加。 |
|
JX-D1-0016 |
几何关系 |
机构学/机器人学 |
v=J(q)q˙ |
速度雅可比矩阵及其广义逆 |
1. 微分运动关系:机器人末端操作速度 v=[vpT,ωT]T与关节速度 q˙之间呈线性关系,系数矩阵即为雅可比矩阵 J(q)。 |
基于微分运动学,在瞬时构型下精确。误差来源于关节间隙、连杆柔性。奇异性是固有特性。 |
微分几何、线性代数(矩阵的伪逆) |
场景:机器人轨迹规划中的速度控制、力控制、奇异规避、冗余度解析。 |
v:末端操作速度(6×1)。 |
线性代数:核心是矩阵乘法、求逆/伪逆。 |
冗余机器人逆运动学求解时序(梯度投影法): |
描述了从关节空间速度到任务空间速度的“运动流”映射及其逆映射。雅可比矩阵 J是这个映射的“传导率矩阵”。正运动学是速度从关节向末端的“正向流动”。逆运动学是速度从末端向关节的“反向分配”。伪逆 J+实现了在最小能量(最小关节速度范数)意义下的“最优分配流”。零空间 (I−J+J)代表了不影响末端运动的“内部循环流”,可用于系统内部的优化调节。 |
|
JX-D1-0017 |
振动方程 |
结构动力学/模态分析 |
Mx¨(t)+Cx˙(t)+Kx(t)=F(t) |
多自由度系统振动微分方程 |
1. 系统离散化:将连续结构离散为有限个自由度(如集中质量、有限元节点),位移向量 x(t)=[x1(t),x2(t),...,xn(t)]T。 |
模型精度取决于离散化程度和 M,C,K矩阵的准确性。阻尼矩阵 C通常难以精确获得,常假设为比例阻尼(Rayleigh阻尼):C=αM+βK。 |
拉格朗日力学、离散系统理论 |
场景:建筑抗震分析、车辆NVH分析、机床颤振预测、航空航天结构动态响应。 |
M:n×n质量矩阵(通常对称、正定)。 |
常微分方程组:核心是耦合的线性二阶ODE。 |
模态分析与时程分析流程: |
描述了振动能量在多个自由度间的流动与转化。质量矩阵 M存储动能,刚度矩阵 K存储势能,阻尼矩阵 C耗散能量。外力 F(t)是外部能量输入。模态分析揭示了系统固有的能量流动模式(振型),每个振型对应一种能量在自由度间特定比例的分配方式。系统的总响应是这些模态“能量流”的线性叠加。 |
|
JX-D1-0018 |
稳定性判据 |
控制系统/动力学系统 |
Δ(s)=ansn+an−1sn−1+...+a1s+a0=0 |
劳斯-赫尔维茨稳定性判据 |
1. 系统特征方程:线性时不变系统稳定性由其传递函数分母或状态矩阵特征方程 Δ(s)=0的根(极点)决定。系统稳定的充要条件是所有根具有负实部(左半平面)。 |
判据精确(针对线性系统)。它避免了直接求解高次方程的根,是代数判据。对含参数的系统,可用于分析稳定区间。 |
多项式根的性质、线性系统理论 |
场景:控制系统参数稳定性分析、判断特征方程根的位置、预测系统是否发散。 |
Δ(s):系统特征多项式。 |
代数:基于多项式系数进行表格运算。 |
劳斯稳定性分析流程: |
描述了系统极点(模态)在复平面上的“分布流”。稳定的系统,其所有极点(能量模态)都位于左半平面,意味着任何扰动引起的自由响应都会指数衰减(能量耗散)。劳斯判据通过检验特征多项式系数的组合(劳斯表第一列),间接监控了这些极点是否企图“流入”右半平面(实部为正),右半平面的极点对应指数增长的模式(能量输入超过耗散),导致系统失稳。 |
|
JX-D1-0019 |
流体动力学方程 |
计算流体力学 |
∂t∂ρ+∇⋅(ρu)=0 |
可压缩纳维-斯托克斯方程组 |
1. 质量守恒:同上JX-D1-0006。 |
是流体运动最基本的物理定律,形式精确。数值求解的精度和稳定性受离散格式、湍流模型、网格质量等极大影响。对高雷诺数流动,直接求解(DNS)计算量巨大,常用RANS或LES模型。 |
质量、动量、能量守恒定律、牛顿粘性定律、傅里叶热传导定律 |
场景:航空航天(翼型绕流、发动机燃烧)、汽车工业(外流场、发动机舱冷却)、能源动力(涡轮机械、燃烧室)。 |
ρ:密度。 |
非线性偏微分方程组:对流项 (u⋅∇)u是非线性的来源。 |
有限体积法求解流程(以稳态问题为例): |
完整描述了质量流、动量流和能量流在流体域中的对流、扩散、产生与耗散。 |
|
JX-D1-0020 |
优化算法 |
机械设计/运筹学 |
minx∈Rnf(x) |
梯度下降法(最速下降法) |
1. 问题建模:将工程设计问题(如轻量化、性能最优)转化为含目标函数 f(x)和约束 gi(x),hj(x)的数学优化问题。 |
f(\mathbf{x}{k+1}) - f(\mathbf{x}k) |
< \delta$ 时停止迭代。 |
收敛速度慢(线性收敛),特别是在谷底呈“之”字形振荡。对初始值和步长敏感。是许多高级优化算法(如共轭梯度法、Adam)的基础。 |
微积分(梯度方向是函数最速上升方向)、泰勒展开 |
场景:机器学习模型训练、结构参数优化、控制系统参数整定等无约束或简单约束的连续优化问题。 |
x:n维设计变量向量。 |
微积分/向量分析:核心是梯度。 |
|
JX-D1-0021 |
接触力学模型 |
机械连接/摩擦学 |
FN=Kδ+Dδ˙ |
罚函数法接触模型(含库仑摩擦) |
1. 接触检测:判断从点与主面是否发生穿透。定义法向间隙 gN, 若 gN<0则发生接触。 |
是物理接触的简化模型。精度取决于罚刚度 K(理论上 K→∞精确,但数值上会导致病态)。摩擦模型是对复杂微观滑移的宏观近似。计算效率高,广泛应用于显式动力学。 |
赫兹接触理论(小变形)、库仑摩擦定律、惩罚原理 |
场景:多体动力学中的碰撞与接触(齿轮啮合、轴承滚子-滚道、车门关闭)、机械装配干涉分析。 |
gN:法向间隙(负值表示穿透)。 |
非光滑动力学:摩擦力在滑动开始点不连续。 |
显式动力学中接触力计算流程(每个时间步): |
描述了接触界面间力与运动的“约束流”。法向接触像一根非线性的弹簧-阻尼器,仅在受压时传递“力流”。切向摩擦则像一个带屈服极限的剪切弹簧-阻尼器:在屈服前(粘着),传递弹性“剪力流”;达到屈服极限后(滑动),“剪力流”被限制在最大值 μFN, 多余的相对运动能量被耗散。状态切换点就是“力流”通道从弹性传递切换到限幅耗散的临界点。 |
|
JX-D1-0022 |
波动方程 |
声学/弹性动力学 |
∂t2∂2u=c2∇2u |
经典波动方程 |
1. 物理背景:描述微小扰动在弹性介质(弦、杆、膜、声场)中的传播。 |
|
编号 |
类别 |
领域 |
模型配方 |
定理/算法/模型/方法名称 |
定理/算法/模型/方法的逐步思考推理过程及每一个步骤的数学方程式和参数选择/参数优化 |
精度/误差/强度 |
底层规律/理论定理 |
典型应用场景和各类特征 |
变量/常量/参数列表及说明 |
数学特征 |
时序和交互流程的所有细节/分步骤时序情况及数学方程式 |
流动模型和流向方法的数学描述 |
|---|---|---|---|---|---|---|---|---|---|---|---|---|
|
JX-D1-0022 |
控制原理 |
线性系统理论 |
G(s)=U(s)Y(s)=ansn+...+a1s+a0bmsm+...+b1s+b0 |
传递函数 |
1. 系统线性化:在平衡点附近,对非线性系统进行线性化,得到线性常系数微分方程:any(n)+...+a1y˙+a0y=bmu(m)+...+b1u˙+b0u。 |
仅在系统为线性时不变(LTI)且在零初始条件下精确成立。对于非线性系统,是其在小信号扰动下的近似模型。 |
拉普拉斯变换理论、线性系统理论 |
场景:控制系统分析与综合、频域响应(伯德图、奈奎斯特图)绘制、稳定性分析、控制器设计。 |
G(s):传递函数。 |
复变函数:G(s)是复变量s的有理函数。 |
从物理系统到传递函数的建模流程: |
G(j\omega) |
|
JX-D1-0023 |
稳定性判据 |
频域控制理论 |
Z=N+P |
奈奎斯特稳定性判据 |
1. 原理基础:幅角原理。映射定理:当s沿奈奎斯特路径(包围整个右半平面)顺时针绕行一周时,开环传递函数 G(s)H(s)在 G(s)H(s)平面上环绕点 (−1,j0)的净圈数 N与闭环系统右半平面极点 Z和开环右半平面极点 P满足 Z=N+P。 |
判据精确,适用于线性时不变系统。可以处理延迟环节 e−sτ, 这是劳斯判据难以处理的。图形化方法还能提供稳定裕度(幅值裕度、相位裕度)信息。 |
复变函数(幅角原理)、反馈系统理论 |
场景:判断闭环控制系统稳定性,尤其是含延迟或难以获得特征多项式的系统;评估系统相对稳定性(稳定裕度)。 |
G(s)H(s):开环传递函数。 |
复分析/映射:核心是幅角原理。 |
奈奎斯特稳定性分析流程: |
GH |
|
JX-D1-0024 |
数值方法 |
计算固体力学 |
Ka=F |
有限元法基本方程(线弹性静力) |
1. 强形式到弱形式:将平衡方程 ∇⋅σ+b=0及其边界条件,通过加权残值法(如 Galerkin 法)转化为积分弱形式:∫ΩδϵTσdΩ=∫ΩδuTbdΩ+∫ΓtδuTtdΓ。 |
精度取决于:网格密度(h-收敛)、单元阶次(p-收敛)、几何与物理模型准确性。误差有离散误差、数值积分误差、舍入误差。可通过后验误差估计进行网格自适应。 |
变分原理(最小势能原理或虚功原理)、加权残值法 |
场景:复杂几何结构的静力、动力、热、流固耦合等物理场分析。是CAE软件(如ANSYS, ABAQUS)的核心。 |
K:整体刚度矩阵(对称、正定、稀疏)。 |
变分法/泛函:基于能量最小化或加权残值。 |
有限元分析通用流程: |
描述了外力在离散化结构中的“传递与平衡流”。刚度矩阵 K是结构内部抗力通道的传导率张量,它将节点位移“流” a映射为节点力“流” Ka。方程 Ka=F意味着外部载荷 F流入结构后,通过内部的所有“传导路径”(单元),最终被结构内力 Ka完全平衡。每个单元的贡献是局部“传导路径”的描述,组装过程就是将所有局部路径连接成全局网络。 |
|
JX-D1-0025 |
数值方法 |
计算流体力学 |
∂t∂U+∂x∂F+∂y∂G+∂z∂H=S |
有限体积法基本方程(积分守恒形式) |
1. 守恒定律的积分形式:任何流动的守恒律(质量、动量、能量)均可写成上述散度形式。对固定控制体 V积分,利用散度定理,将体积分转化为通过控制体表面 A的通量积分。 |
精度取决于:网格质量、通量计算格式的阶数、时间积分格式。守恒性好,能精确捕捉激波等间断。湍流模型引入主要建模误差。 |
物理守恒定律(积分形式)、散度定理 |
场景:复杂外形的流体动力学计算,特别是涉及激波、分离流、多相流的场合。是主流CFD软件(如Fluent, STAR-CCM+)的核心方法。 |
U:守恒变量向量 [ρ,ρu,ρv,ρw,ρE]T。 |
守恒型偏微分方程:方程写成散度形式。 |
有限体积法求解流程(稳态问题,隐式求解): |
\mathbf{R}_i |
|
JX-D1-0026 |
磨损模型 |
摩擦学/可靠性 |
V=KHWs |
阿查德磨损定律(Archard) |
1. 假设:磨损是由于表面微凸体接触产生塑性变形和材料迁移所致。真实接触面积与载荷成正比,与材料硬度成反比。 |
是一个高度简化的经验-半经验模型。在中等载荷、稳定滑动条件下,对多种材料组合的磨损体积趋势预测有效。精度严重依赖磨损系数 K的准确获取,而 K本身受多因素影响。 |
接触力学、塑性变形、材料剥离的物理模型 |
场景:机械密封、轴承、齿轮、制动盘、滑轨等存在滑动接触的零件的寿命预估与设计。 |
V:磨损体积(m³)。 |
比例关系:磨损体积与载荷、滑动距离成正比,与硬度成反比。 |
基于磨损模型的寿命预测流程: |
描述了材料以磨屑形式从接触表面“流失”的速率。磨损系数 K可以理解为每次滑动中,真实接触面积内“转化为磨屑的材料比例”。方程 V˙=K(Wv/H)表明,材料流失率 正比于能量输入率(Wv近似为摩擦功率)除以材料抵抗变形的能力(H)。这是一个典型的“输入-耗散”流模型,输入是机械功,耗散是材料损失。 |
|
JX-D1-0027 |
运动规划 |
机器人学/自动化 |
minq(t)∫0tf(∑i=1nwiq˙i2)dt |
关节空间轨迹规划(三次/五次多项式插值) |
1. 问题定义:在关节空间生成一条从起点 qs到终点 qg的光滑轨迹 q(t), 满足运动约束(速度、加速度极限),并优化某种性能指标(如时间最短、能量最小、平滑)。 |
多项式轨迹在路径点上精确通过,且速度/加速度连续。这是对理想路径的近似,未考虑动力学(如惯性力、哥氏力)约束,实际跟踪时可能需要模型预测控制。 |
多项式插值理论、最优控制理论 |
场景:工业机器人点对点运动、机械臂拾放操作、数控机床的轮廓加工。 |
q(t):关节角度向量(n维函数)。 |
多项式代数:轨迹是时间多项式。 |
三次多项式轨迹生成流程: |
描述了关节角度从起始值到目标值的平滑“流动”过程。多项式函数定义了角度随时间变化的“流动曲线”。三次多项式确保了速度的连续性(没有突变),就像定义了“流动”的斜率是平滑变化的。五次多项式进一步确保了加速度(流动曲线的曲率)的连续性,使得“流动”更加平滑,减少了冲击。轨迹规划就是设计这条“流动曲线”的形状,使其满足边界和约束。 |
|
JX-D1-0028 |
流体模型 |
润滑/轴承技术 |
∂x∂(μh3∂x∂p)+∂z∂(μh3∂z∂p)=6U∂x∂h+12∂t∂h |
雷诺方程(流体动压润滑) |
1. 简化假设:基于 Navier-Stokes 方程,针对薄层润滑膜(h≪L)作如下简化:忽略体积力、惯性力;压力沿膜厚方向不变;流速沿膜厚呈抛物线分布(Couette-Poiseuille流);无滑移边界;流体为牛顿流体。 |
_{y=0}=0,v |
_{y=h}=V),得到上述雷诺方程。<br>∗∗3.各项物理意义∗∗:<br>−左边:压力引起的流动(Poiseuille流)的扩散项。<br>−右边第一项:表面切向运动产生的楔形效应(U \partial h/\partial x)。<br>−右边第二项:挤压膜效应(\partial h/\partial t$)。 |
在假设成立的范围内(等粘、不可压、层流、薄膜)是精确的。实际应用中,需考虑端泄、空化、变粘度、表面粗糙度、湍流等修正。 |
纳维-斯托克斯方程、连续方程、薄层流体近似 |
场景:滑动轴承、径向轴承、推力轴承、导轨、活塞环-缸套、齿轮啮合等流体动压润滑分析。 |
p:油膜压力 (Pa)。 |
偏微分方程:椭圆型PDE。 |
|
JX-D1-0029 |
热传导方程 |
传热学/热力学 |
ρcp∂t∂T=∇⋅(k∇T)+q˙v |
傅里叶热传导定律与能量方程 |
1. 傅里叶定律:热流密度向量 q与温度梯度成正比,方向相反:q=−k∇T, k为热导率。 |
在连续介质和各向同性材料中精确。实际误差来源于:k,cp,ρ随温度变化;接触热阻;复杂边界条件(如辐射)的近似。 |
能量守恒定律、傅里叶热传导经验定律 |
场景:发动机缸体热分析、电子设备散热、热处理工艺模拟、建筑保温计算。 |
T:温度 (K)。 |
偏微分方程:扩散方程。 |
瞬态热分析有限元流程: |
描述了热量在介质中的“扩散流”。傅里叶定律 q=−k∇T表明,热流从高温区流向低温区,流速(热流密度)正比于温度梯度(驱动力),比例系数是热导率 k。热传导方程 ρcp∂t∂T=∇⋅(k∇T)+q˙v是热量守恒方程:左边是内能(热量)的储存率,右边第一项是净流入的导热热量(扩散),第二项是内部产热(源)。 |
|
JX-D1-0030 |
动力学算法 |
多体系统动力学 |
[MΦqΦqT0][q¨λ]=[Qγ] |
带拉格朗日乘子的多体系统动力学方程(指标-1 DAE) |
1. 系统描述:用广义坐标 q描述位形,系统受完整约束 Φ(q,t)=0。 |
方程在加速度级精确。数值积分会产生约束漂移(即 Φ=0和 Φ˙=0被破坏),需采用稳定化方法(如Baumgarte修正、投影法)。 |
达朗贝尔原理、拉格朗日乘子法、微分-代数方程理论 |
场景:包含复杂运动副(铰链)的多体系统动力学仿真,如汽车悬架、机器人闭链机构、航天器太阳能帆板。 |
q,q˙,q¨:广义坐标、速度、加速度向量。 |
线性代数:求解一个稀疏的对称/非对称线性系统。 |
多体系统动力学仿真流程(常采用): |
\Phi |
|
JX-D1-0031 |
故障诊断模型 |
状态监测/可靠性工程 |
x(t)=∑k=1KAk(t)ej(2πfkt+ϕk)+n(t) |
信号分解与频谱分析模型(傅里叶级数/变换) |
1. 周期性故障特征:旋转机械的局部故障(如轴承点蚀、齿轮断齿)会产生周期性冲击,激励结构固有频率,在振动信号中表现为以故障特征频率 fc及其谐波为载波,以固有频率 fn为调制的调幅-调频信号。 |
X(f) |
^2$。 |
傅里叶分析对平稳周期信号理论精确。实际中频谱泄露、栅栏效应、噪声会影响精度。对非平稳信号,传统频谱分析会模糊时频特征。 |
傅里叶分析理论、信号处理、调制原理 |
场景:旋转机械(轴承、齿轮箱、风机、泵)的故障检测与诊断。 |
x(t):时域振动/声发射信号。 |
谐波分析:将信号分解为正弦波的叠加。 |
|
JX-D1-0032 |
相似准则 |
实验力学/流体力学 |
ρU2L2F=f(μρUL,gLU,aU,...) |
量纲分析与π定理(白金汉π定理) |
1. 确定相关物理量:分析物理现象,列出所有影响该现象的物理量,共 n个。例如,物体在流体中受到的阻力 F与密度 ρ、流速 U、特征长度 L、粘度 μ、重力加速度 g、声速 a等有关。 |
π定理是精确的数学定理。相似准则的实用性取决于是否正确地选择了所有重要的相关物理量。完全相似(所有π数相等)通常难以实现,需根据主导作用保证部分相似。 |
量纲齐次原理、群论(李群) |
场景:设计模型实验(风洞、水洞),将模型实验结果换算到原型;指导计算流体力学中的无量纲化;整理实验数据得到经验公式。 |
物理量:F,ρ,U,L,μ,g,a,...。 |
线性代数:量纲矩阵的秩与零空间。 |
应用π定理推导阻力公式的步骤: |
提供了在不同尺度或条件下,物理现象之间“相似性”的判据。每个无量纲数(如 Re,Fr,Ma)代表了两种竞争效应的比值。例如,雷诺数 Re=μρUL是惯性力流与粘性力流的比。π定理保证了,如果两个系统的所有无量纲数相同,那么尽管它们的绝对尺寸、速度等不同,但其内部“各种力量(或能量)的流动格局”是几何相似的,因此物理现象相似。模型实验就是基于此,创造与原型具有相同“力流比”的环境。 |
|
JX-D1-0033 |
断裂判据 |
断裂力学/损伤容限 |
KI=Yσπa |
线弹性断裂力学 - 应力强度因子与断裂韧性 |
1. 裂纹尖端场:Irwin等人分析表明,I型(张开型)裂纹尖端附近的应力场具有奇异性,形式为 σij=2πrKIfij(θ)+..., 其中 r,θ为以裂纹尖端为原点的极坐标。 |
在小范围屈服条件下是精确的工程判据。当塑性区较大时,需采用弹塑性断裂力学(如J积分)。KIC的测试有严格标准。 |
线弹性理论(Williams 渐近展开)、能量释放率理论(Griffith理论) |
场景:含缺陷结构(压力容器、管道、飞机结构、焊接接头)的安全评定、疲劳裂纹扩展寿命预测、材料筛选。 |
KI:I型应力强度因子 (MPa·√m)。 |
奇异积分:裂纹尖端应力场具有 1/r奇异性。 |
含裂纹结构安全评估流程: |
描述了导致裂纹扩展的“驱动力”的流动与集中。应力强度因子 KI是裂纹尖端应力奇异性强度的度量,可以理解为“驱动裂纹扩展的能量流”的强度。它汇聚了远场载荷 σ和几何(a,Y)的贡献。断裂韧性 KIC则是材料抵抗这种扩展的能量耗散能力的阈值。断裂判据 KI≥KIC意味着,当“驱动力流”的强度超过材料的“耗散能力”时,富余的能量流将用于创造新表面,驱动裂纹前进。 |
|
JX-D1-0034 |
优化算法 |
机器学习/控制 |
J(θ)=m1∑i=1mL(f(x(i);θ),y(i)) |
Adam 优化算法 |
1. 问题:最小化目标函数 J(θ), 如神经网络损失函数,其中 θ为参数。 |
在实践中对许多非凸优化问题(如训练DNN)非常有效,收敛速度快。其性能依赖于超参数 α,β1,β2的选择。理论收敛性在凸优化框架下得到证明。 |
随机梯度下降、动量法、RMSProp |
场景:深度神经网络训练、强化学习策略优化、大规模机器学习模型参数训练。 |
θ:待优化参数向量。 |
随机优化:基于小批量数据的随机梯度。 |
Adam 算法迭代流程: |
描述了在参数空间中寻找最优点的、智能化的“下降流”。与简单梯度下降的“盲流”不同,Adam引入了两个“记忆流”: |
|
JX-D1-0035 |
可靠性模型 |
系统可靠性工程 |
Rs(t)=∏i=1nRi(t)(串联) |
串联与并联系统可靠性模型 |
1. 可靠性定义:部件或系统在规定的条件下和规定的时间内,完成规定功能的概率,记为 R(t)。 |
在部件失效独立的假设下精确。实际系统中,部件失效可能存在相关性(共因失效),模型会高估可靠度。指数分布假设(恒定失效率)适用于电子产品的随机失效期,不适用于机械磨损期。 |
概率论(独立事件的概率)、布尔代数 |
场景:复杂装备(如飞机、发电厂、通信网络)的可靠性设计、预测与评估;冗余设计(如备份系统)的效果量化。 |
Ri(t):第 i个部件在时刻 t的可靠度。 |
概率:可靠度是概率。 |
系统可靠性评估流程: |
描述了系统功能成功的“概率流”如何通过各部件传递。 |
|
编号 |
类别 |
领域 |
模型配方 |
定理/算法/模型/方法名称 |
定理/算法/模型/方法的逐步思考推理过程及每一个步骤的数学方程式和参数选择/参数优化 |
精度/误差/强度 |
底层规律/理论定理 |
典型应用场景和各类特征 |
变量/常量/参数列表及说明 |
数学特征 |
时序和交互流程的所有细节/分步骤时序情况及数学方程式 |
流动模型和流向方法的数学描述 |
|---|---|---|---|---|---|---|---|---|---|---|---|---|
|
JX-D1-0036 |
控制理论 |
状态空间分析 |
x˙=Ax+Bu |
线性时不变系统的状态空间方程 |
1. 状态变量选择:选择一组能够完全确定系统未来行为的最小变量集 x=[x1,x2,...,xn]T。 |
模型是线性时不变系统在小扰动下的精确或近似描述。精度取决于状态变量选择的完备性和线性化的有效性。 |
线性系统理论、常微分方程理论 |
场景:多输入多输出系统的分析与控制设计、最优控制(LQR)、状态估计(卡尔曼滤波)、现代控制系统设计。 |
x:状态向量 (n×1)。 |
线性代数:核心是矩阵运算。 |
从状态空间到响应的求解流程: |
描述了系统内部状态 x的“演化流”。系统矩阵 A定义了状态变量之间如何相互影响和自演化。输入矩阵 B定义了外部控制信号 u如何“注入”并驱动状态流。输出方程定义了内部状态流中哪些部分以及如何被“观测”为输出 y。状态空间表示清晰地勾勒了信息(状态)在系统内部的流动路径。 |
|
JX-D1-0037 |
动力学算法 |
机器人学/多体动力学 |
正向迭代:iωi=i−1RiT(i−1ωi−1+θ˙ii−1z^i−1)iω˙i=i−1RiT(i−1ω˙i−1+i−1θ¨ii−1z^i−1+i−1ωi−1×(θ˙ii−1z^i−1))iv˙i=i−1RiT[i−1v˙i−1+i−1ω˙i−1×i−1pi+i−1ωi−1×(i−1ωi−1×i−1pi)]+d¨ii−1z^i−1+...iv˙Ci=iω˙i×ipCi+iωi×(iωi×ipCi)+iv˙i逆向迭代:ifi=iFi+i+1Rii+1fi+1ini=iNi+i+1Rii+1ni+1+ipCi×iFi+ipi+1×(i+1Rii+1fi+1)τi=iniT(i−1RiTi−1z^i−1) (旋转关节) |
牛顿-欧拉迭代动力学算法 |
1. 正向迭代(外推):从基座(连杆0)到末端,递推计算每个连杆的角速度、角加速度、线加速度和质心加速度。 |
在刚体假设和运动副理想约束下精确。计算是数值稳定的。实际误差来源于质量、惯性参数标定误差,以及忽略关节摩擦、柔性。 |
牛顿-欧拉方程(刚体动力学)、运动学递推 |
场景:工业机器人、仿人机器人的实时模型计算、动力学仿真、基于模型的控制(如计算力矩控制)。 |
i−1Ri:从连杆 i坐标系到连杆 i−1坐标系的旋转矩阵。 |
向量/张量运算:涉及叉乘、旋转矩阵变换。 |
牛顿-欧拉算法完整时序流程: |
描述了运动、力和力矩在机器人运动链中的双向流动。 |
|
JX-D1-0038 |
机构学公式 |
平面机构学 |
F=3(n−1)−2j1−j2 |
平面机构自由度计算公式(格鲁勃-契贝谢夫公式) |
1. 自由度定义:机构所具有的独立运动的数目,或确定机构位置所需独立广义坐标的数目。 |
公式在理想情况下精确。实际机构中,间隙、弹性变形可能导致“过约束”或“欠约束”,影响实际运动确定性。 |
运动学、约束理论、螺旋理论 |
场景:新机构的概念设计、运动可行性分析、判断机构是否具有确定运动。 |
F:机构的自由度(Degree of Freedom)。 |
计数与代数:基于构件和运动副的计数进行简单代数运算。 |
机构自由度分析流程: |
描述了机构中“可动性”或“运动潜力”的资源分配。每个活动构件贡献3个“运动资源”(自由度)。每个运动副消耗一定数量的“运动资源”以建立构件间的连接(约束),低副消耗2个,高副消耗1个。机架固定消耗3个。机构的总体“可动性” F就是初始总资源减去所有消耗后剩余的独立“运动资源”数量,也就是需要外部提供多少“驱动能量流”才能完全确定整个系统的运动状态。 |
|
JX-D1-0039 |
电磁学方程 |
电机与驱动 |
vaeaTeJdtdωm+Bωm=Raia+Ladtdia+ea=Keωm=Ktia=Te−TL |
直流有刷电机数学模型 |
1. 电枢回路方程:电枢电压 va平衡电阻压降 Raia、电感压降 Ladia/dt和反电动势 ea。反电动势与转速成正比:ea=Keωm, Ke为反电动势常数。 |
模型是线性的,忽略了磁饱和、齿槽效应、电刷压降等非线性因素。在小信号分析和控制设计中有良好精度。参数 Ra,Kt,Ke,J可通过测试获得。 |
电路理论(基尔霍夫电压定律)、电磁感应定律(法拉第定律)、安培力定律、刚体转动定律 |
场景:直流伺服电机控制系统设计与分析、机器人关节驱动建模、数控机床进给系统建模。 |
va,ia:电枢电压 (V) 和电流 (A)。 |
常微分方程组:电和机械的耦合ODE。 |
直流电机速度控制系统仿真流程: |
清晰描述了能量在电、磁、机械域之间的流动与转换。 |
|
JX-D1-0040 |
流体方程 |
空气动力学/流体机械 |
p+21ρv2+ρgh=constant |
伯努利方程(沿流线,无粘、不可压、稳态) |
1. 假设条件:理想流体(无粘)、不可压缩、稳态流动、质量力有势(如重力)、沿同一流线。 |
在假设条件下精确。实际应用误差来源于粘性效应、可压缩性、非稳态性以及流线不同。常用于定性分析和近似估算。 |
牛顿第二定律(欧拉方程)、机械能守恒 |
场景:管道流动测量(文丘里管、孔板)、机翼升力原理的定性解释、水泵和风机进出口分析、空速管原理。 |
p:静压 (Pa)。 |
代数方程:沿流线的能量代数关系。 |
应用伯努利方程分析文丘里管的流程: |
完美诠释了流体机械能沿流线的守恒流动与转换。方程 p+21ρv2+ρgh=C表明,在无粘、不可压、稳态流动中,流线上任意一点的压力能、动能和势能可以互相转换,但总和不变。这就像一条“能量流”管道,截面积变化(A)会导致流速 v变化,从而引起动能 21ρv2和压力能 p之间的此消彼长(文丘里效应)。速度大的地方压力小,速度小的地方压力大,这是伯努利原理的核心,描述了机械能在不同形式间的“流动”与“再分配”。 |
|
JX-D1-0041 |
弹性理论 |
弹性力学 |
{∂x∂σxx+∂y∂τxy+bx=0∂x∂τxy+∂y∂σyy+by=0 |
平面问题基本方程(平衡、几何、物理) |
1. 平面应力与平面应变:针对特定几何和载荷的简化。 |
在假设成立的条件下精确。平面应力适用于厚度远小于其他尺寸的薄板;平面应变适用于长度远大于横截面尺寸的柱体。是许多弹性力学解析解和有限元分析的基础。 |
连续介质力学、线弹性理论 |
场景:薄板弯曲(近似)、长坝、滚子接触、齿轮齿根应力等二维结构分析。 |
u,v:x,y方向的位移分量。 |
偏微分方程组:平衡、几何、本构方程组成PDE组。 |
基于应力函数的求解流程: |
描述了应力、应变、位移在二维弹性体中的分布与传递流。 |
|
JX-D1-0042 |
复合材料力学 |
层合板理论 |
[NM]=[ABBD][ϵ0κ] |
经典层合板理论(CLT)本构关系 |
1. 基本假设:直法线假设(Kirchhoff-Love),层间完美粘结,平面应力状态。 |
对于薄板(跨度/厚度 > 10)和小变形情况精度良好。忽略了横向剪切变形,对厚板或某些铺层(如 ±45∘)误差较大,需用一阶剪切变形理论(FSDT)修正。 |
弹性力学、平板理论、等效单层理论 |
场景:复合材料机翼、车身、船体的结构设计与分析。 |
N,M:单位宽度上的面内合力和合力矩向量 (Nx,Ny,Nxy,Mx,My,Mxy)。 |
线性代数:核心是矩阵方程。 |
层合板分析与设计流程: |
描述了外力在复合材料层合板中的“传递与分配流”。刚度矩阵 [AB;BD]是整个层合板的“传导率张量”,它将宏观的“中面变形流” (ϵ0,κ) 映射为宏观的“力流” (N,M)。 |
|
JX-D1-0043 |
振动理论 |
单自由度系统 |
mx¨+cx˙+kx=F0cos(ωt) |
单自由度系统受迫振动方程 |
1. 系统建模:质量 m、弹簧 k、阻尼 c组成的系统,受简谐激振力 F(t)=F0cos(ωt)。 |
模型精确。是理解振动现象的基础。实际系统可能具有非线性(刚度、阻尼),但线性模型在振幅不大时是良好的近似。 |
牛顿第二定律、线性微分方程理论 |
场景:隔振系统设计、振动传感器原理、结构共振评估、动平衡分析。 |
m,c,k:质量、阻尼系数、刚度系数。 |
二阶线性常微分方程:具有解析解。 |
受迫振动响应分析流程: |
X |
|
JX-D1-0044 |
非线性动力学 |
自激振动 |
x¨−μ(1−x2)x˙+x=0 |
范德波尔方程 |
1. 方程形式:典型的非线性阻尼(或称“负阻尼”)方程。阻尼项为 −μ(1−x2)x˙。 |
x |
< 1时,(1-x^2) > 0,阻尼项为负(-\mu(...)\dot{x}),系统获得能量,振幅增大。<br>−当 |
x |
> 1时,(1-x^2) < 0,阻尼项为正,系统耗散能量,振幅减小。<br>∗∗3.极限环∗∗:上述特性导致系统存在一个稳定的周期解,即极限环。无论初始条件如何,系统最终都会稳定到该周期运动上。<br>∗∗4.参数影响∗∗:\mu是小参数时,极限环近似为正弦波;\mu$ 较大时,极限环是 relaxation oscillation(张弛振荡),波形剧烈变化。 |
方程是自激振动的一个经典数学模型。其解(极限环)的周期和形状可以通过近似方法(如平均法、多尺度法)或数值积分得到。 |
非线性振动理论、动力系统理论(极限环、稳定性) |
场景:电子振荡器、机械系统的摩擦颤振(如刹车尖叫)、心脏的电生理模型、某些化学反应振荡。 |
|
JX-D1-0045 |
热力学定律 |
工程热力学 |
∮TδQ≤0 |
克劳修斯不等式(热力学第二定律的数学表述) |
1. 热力学第二定律:指出热过程的方向性,如热不能自发地从低温传到高温。 |
是热力学的基本定律,普遍成立。在实际计算中,需要准确知道热交换和温度历程。 |
热力学第二定律 |
场景:热机(内燃机、燃气轮机)最大效率分析、制冷循环性能系数计算、化工过程可行性判断、能量系统优化。 |
S:熵 (J/K)。 |
循环积分:对闭合路径的积分。 |
应用克劳修斯不等式分析热机循环: |
描述了热量传递与转换过程中“品质”的退化或“无序度”(熵)的流动。 |
|
JX-D1-0046 |
流体方程 |
边界层理论 |
u∂x∂u+v∂y∂u=−ρ1dxdp+ν∂y2∂2u |
普朗特边界层方程(不可压,稳态,二维) |
1. 边界层概念:高雷诺数下,粘性影响仅局限在物面附近的薄层内,层外为无粘势流。 |
在高雷诺数 (Re>104) 的平板、楔形流等情况下精度很高。对于分离流、强压力梯度情况,可能需要更复杂的模型。 |
纳维-斯托克斯方程、量纲分析与摄动理论 |
场景:飞机机翼、涡轮叶片、船舶外壳的摩擦阻力计算、传热与传质分析、流动分离预测。 |
u,v:边界层内 x(流向)和 y(法向)速度分量。 |
偏微分方程:简化的PDE组。 |
应用卡门动量积分方程近似求解流程: |
_{y=0} = \frac{\mu U}{\delta} f'(0)。<br>6.∗∗代入动量积分方程∗∗:得到一个关于\delta(x)的常微分方程。<br>7.∗∗求解∗∗\delta(x):给定初始条件(如前缘\delta(0)=0),数值积分求解\delta(x),进而得到\theta(x), \tau_w(x)$ 等。 |
|
JX-D1-0047 |
声学方程 |
声学/振动 |
∇2p−c021∂t2∂2p=0 |
声波方程(线性化、无源、均匀介质) |
1. 基本假设:介质为理想气体,声波引起的小扰动(压力 p′, 密度 ρ′, 速度 u)满足:p′≪p0, ρ′≪ρ0, u≪c0。 |
在小振幅声波的假设下是精确的。当声压级很高(如 >140 dB)或介质不均匀、有流动时,需用非线性声学或对流波动方程。 |
流体力学方程的线性化、气体状态方程 |
场景:室内声学设计、噪声控制、声纳、超声检测、麦克风和扬声器建模。 |
p:声压扰动 (Pa)。 |
偏微分方程:双曲型波动方程。 |
平面波、球面波、柱面波基本解推导: |
描述了声压扰动(即机械振动)在弹性介质中的“传播波”。波动方程 ∇2p=c021∂t2∂2p本质是动量守恒与质量守恒耦合的结果。左边 ∇2p代表压力的空间变化(“压力聚集度”),右边代表压力的时间二阶变化(“加速”)。方程表明,空间上某点的压力聚集会产生该点压力的时间变化(反之亦然),这种关系以波速 c0为特征,导致扰动以“波”的形式传播出去。这是一种能量和信息的波动式流动。 |
|
JX-D1-0048 |
电磁学方程 |
宏观电磁学 |
∇⋅D∇×E∇⋅B∇×H=ρf=−∂t∂B=0=Jf+∂t∂D |
麦克斯韦方程组(微分形式) |
1. 高斯电场定律:电位移 D的散度等于自由电荷密度 ρf。描述了电场的有源性。 |
|
编号 |
类别 |
领域 |
定理/定律/原理名称 |
核心表述(公式/叙述) |
定理的逐步思考推理过程、证明思路及关键数学步骤 |
成立条件与适用范围 |
底层规律/理论归属 |
典型应用场景和工程意义 |
变量/参数列表及说明 |
数学与逻辑特征 |
定理的演绎/应用流程步骤 |
物理本质与“流”描述 |
|---|---|---|---|---|---|---|---|---|---|---|---|---|
|
TH-D1-0001 |
静力学 |
刚体力学 |
力的平行四边形法则 |
作用在同一点上的两个力 F1和 F2, 其合力 R可以用以这两个力为邻边构成的平行四边形的对角线来表示。即:R=F1+F2(向量加法)。 |
1. 实验基础:通过弹簧秤实验验证,两个力共同作用的效果与一个力单独作用的效果相同,该力即合力。 |
条件:力作用于同一刚体上的同一点。 |
欧几里得几何、向量代数 |
场景:任何力的合成与分解分析。例如,斜拉桥索力的合成、机械构件上多载荷的等效、运动学中速度与加速度的合成。 |
F1,F2:两个分力矢量。 |
几何性:基于几何作图。 |
1. 作图法:按比例尺画出两个力,做平行四边形,量取对角线长度和方向。 |
描述了力作为一种矢量“流”的叠加原理。两个力矢量如同两条有方向的“作用流”,在作用点汇合。平行四边形法则定义了这两条“流”如何合并为一条等效的“总作用流”。合力矢量代表了两个分力对物体平动效应(动量变化率)的矢量和,是力流在向量空间中的线性叠加。 |
|
TH-D1-0002 |
静力学 |
刚体力学 |
二力平衡原理 |
作用于同一刚体上的两个力,使刚体保持平衡的充分必要条件是:这两个力大小相等、方向相反、且作用在同一条直线上。即:F1=−F2, 且两力共线。 |
1. 必要性证明(反证法):假设两个力平衡但不共线。根据力的可传性,将它们的作用点移到其作用线的交点。此时,两个力构成一个力偶,力偶矩不为零,将使刚体产生转动,与平衡状态(既无平动也无转动)矛盾。因此,它们必须共线。 |
条件:两个力作用于同一刚体。刚体处于平衡状态(静止或匀速直线运动)。 |
牛顿第一定律、刚体静力学公理 |
场景:桁架中杆件受力分析(二力杆假设)、绳索和链的受力、简单的支撑反力分析。 |
F1,F2:作用在刚体上的两个力。 |
充要条件:既是平衡的必要条件,也是充分条件。 |
1. 识别二力构件:在结构(如桁架)中,找出仅在两处受力且不计自重的构件。 |
描述了力流在刚体内部传递达到静力平衡时的一种最简单、最直接的路径。当刚体只受两个力作用且平衡时,这两个力必须构成一条封闭的、无矩的“力流环”。任何偏离共线的情况都会产生力偶矩,打破平衡,引发旋转的“角动量流”。二力平衡原理定义了这种单向、共线的静力传递通道。 |
|
TH-D1-0003 |
静力学 |
刚体力学 |
加减平衡力系原理 |
在作用于刚体的已知力系上,加上或减去任意的平衡力系,并不改变原力系对刚体的作用效应。 |
1. 公理性陈述:此原理是静力学的基本公理之一,无法从更基础的原理推导,其正确性由大量实践验证。 |
条件:仅适用于刚体。对于变形体,力沿作用线移动会改变其内部应力分布,从而改变效应。 |
静力学公理 |
场景:力系的简化、力向一点平移(需附加力偶)、约束反力的分析、复杂载荷的等效替换。 |
无特定公式,是一个操作原理。 |
公理性:是逻辑起点之一。 |
1. 目标:简化力系或改变力的作用点以方便分析。 |
描述了在刚体模型中“力流”的重新分配与组合不影响整体宏观效应的原理。刚体假设下,内部“力流”路径的细节被忽略,只关心净效应。加减平衡力系,相当于在系统的“力流网络”中,增加或移除一个内部自平衡、不对外产生净“流量”(主矢)和“涡流”(主矩)的局部回路,因此不改变整个网络对外的输运特性。 |
|
TH-D1-0004 |
静力学 |
力偶理论 |
力偶等效定理 |
作用在同一平面内的两个力偶,如果它们的力偶矩(大小和转向)相同,则它们彼此等效。即:力偶对刚体的转动效应完全由力偶矩矢量决定。 |
1. 力偶矩定义:力偶矩 M等于力偶中一个力的大小 F与力偶臂 d(两力作用线之间的垂直距离)的乘积,即 M=±Fd, 正负号表示转向(通常逆时针为正)。 |
条件:力偶作用在刚体上。对于平面力系,力偶矩是代数量;对于空间力系,力偶矩是自由矢量。 |
静力学公理、矢量代数 |
场景:机械中转动副的受力分析、机构中力矩的等效替换、复杂力系简化时力偶的处理。 |
M:力偶矩矢量,其大小 M=Fd, 方向垂直于力偶所在平面,按右手螺旋法则确定。 |
自由矢量性:在空间问题中,力偶矩矢量可以滑移,是自由矢量。 |
1. 计算力偶矩:对于给定的力偶,计算其力偶矩 M=F×d(平面内)或确定其力偶矩矢量 M=r×F(空间,其中 r为从一个力作用点到另一个力作用点的矢量)。 |
描述了纯转动作用的“矩流”的等效性。力偶是一种特殊的力系,其“力流”自成闭环(主矢为零),不产生平动效应,只产生“转动效应流”或“角动量流”。力偶矩矢量 M就是这个“转动效应流”的强度与方向度量。该定理指出,无论产生这个“矩流”的具体力对(F,d)如何,只要其“流量”(大小和方向)相同,对刚体的外部转动效应就完全相同。这类似于电路中的电流源,只要电流值相同,其对外电路的效果就相同。 |
|
TH-D1-0005 |
静力学 |
约束与受力分析 |
三力平衡汇交定理 |
刚体在三个互不平行的力作用下处于平衡时,这三个力的作用线必汇交于同一点。 |
1. 证明思路:设三个力 FA,FB,FC分别作用于刚体上 A,B,C三点,且刚体平衡。 |
条件:刚体;仅受三个力作用;三个力不平行;刚体处于平衡状态。 |
二力平衡原理、力的可传性原理 |
场景:三角支架、三铰拱、某些机构在特定位置(受三个力)的受力分析。 |
F1,F2,F3:作用在刚体上的三个力。 |
几何约束:提供了一个力的方向必须满足的几何条件(共点)。 |
1. 识别:判断研究对象是否仅受三个不平行力作用且平衡。 |
描述了三力平衡时,“力流”必须汇聚于一个“节点”。在三个力作用的平衡系统中,任何两个力的合力必须与第三个力共线且反向以维持平衡。这就要求这三个力的作用线必须相交于一点,使得它们的“力流”线能够在一个共同的“汇交点”上实现矢量闭合。这个“汇交点”是系统内“力流”的静力学奇点,所有力的作用线在此交汇,确保了净力矩为零(因为每个力对该点的力矩均为零)。 |
|
TH-D1-0006 |
静力学 |
力系简化 |
力的平移定理 |
作用在刚体上某点的力,可以平行移动到刚体内任意一点,但必须附加一个力偶,该附加力偶的矩等于原力对新作用点的矩。 |
1. 等效变换过程:设力 F作用于点 A。欲将其平移到点 B。 |
条件:适用于刚体。 |
加减平衡力系原理 |
场景:固定端约束反力的简化(得到一个力和一个力偶)、复杂分布载荷的等效、机构动力学中惯性力系的简化。 |
F:原力矢量。 |
等效性:变换前后力系对刚体的外效应(主矢和主矩)不变。 |
1. 选择简化中心 O。 |
揭示了力对刚体作用的双重性(平动与转动)在空间变换下的分解规律。一个力 F从 A点平移到 B点,其“平动效应流” (F) 被完整地带到 B点。但是,其关于 B点的“转动效应流” (M=r×F) 却丢失了。为了完全等效,必须额外补充一个“转动效应流”(即附加力偶矩)来弥补这个差额。该定理本质是力系效应在改变参考点时的变换规则。 |
|
TH-D1-0007 |
静力学 |
力系平衡 |
平面任意力系的平衡方程 |
刚体在平面任意力系作用下平衡的充分必要条件是:力系的主矢和对任意一点的主矩均为零。即: |
1. 必要性:如果刚体平衡(静止或匀速直线平动、无转动),则根据牛顿第一定律和转动定律,其合力(主矢)和合力矩(对任一点的主矩)必须为零。 |
条件:力系作用在同一个刚体上;刚体处于平面一般运动状态(或静止);力系是平面任意力系。 |
牛顿运动定律、力的平移定理 |
场景:求解梁的支座反力、桁架中指定杆件内力、刚架的内力分析、机械结构的静力学设计校核。 |
∑Fx:所有力在 x轴上投影的代数和。 |
代数方程组:三个独立的线性方程(对于线性系统)。 |
1. 确定研究对象:选取需要计算未知力的物体或系统。 |
描述了平面力系达到静力平衡时,“力流”和“矩流”的全局守恒条件。 |
|
TH-D1-0008 |
材料力学 |
应力分析 |
圣维南原理 |
作用于物体表面某一小区域上的力系,如果用一个静力等效的力系(即主矢和主矩相同)来代替,那么这种替换仅在力系作用区域附近产生显著的应力差异,在距离该区域较远处,应力分布的影响可以忽略不计。 |
1. 定性陈述:局部载荷的等效替换,只影响载荷作用区域附近的应力分布,对“远处”的应力影响很小。 |
条件:替换的力系必须是静力等效的(主矢、主矩相同)。物体是线弹性的(通常)。 |
弹性力学、实验归纳 |
场景: |
无具体公式,是一个定性原理。 |
定性原理:非定量,给出趋势性结论。 |
1. 实际问题:一个构件承受复杂的表面力分布。 |
描述了弹性体内应力“流”的局部扰动特性。复杂的表面载荷会在其作用点附近产生复杂的、高梯度的“应力流场”。圣维南原理指出,这种扰动是一种局部衰减的“流”。在远离扰动源的地方,应力“流”的分布模式主要由结构的整体形状和边界条件决定,而不再“记得”载荷的具体分布细节,只“记得”其静力等效的总效应(主矢和主矩)。这类似于流体中局部扰动产生的涡流会随着距离衰减。 |
|
TH-D1-0009 |
材料力学 |
应力应变关系 |
胡克定律(一维) |
在材料的线弹性范围内,应力 σ与应变 ϵ成正比。即:σ=Eϵ, 其中 E为杨氏模量(弹性模量)。 |
1. 实验基础:罗伯特·胡克通过弹簧实验发现“伸长量与力成正比”(Ut tensio, sic vis)。推广到连续介质,即应力与应变成正比。 |
条件:材料处于线弹性阶段;小变形;通常适用于各向同性材料。 |
固体物理学、连续介质力学 |
场景:几乎所有弹性结构的设计与校核计算:杆的拉压、梁的弯曲、轴的扭转、薄壳的变形等。 |
σ:正应力 (Pa)。 |
线性关系:应力与应变成正比。 |
1. 确定载荷和几何:已知外力 F和杆件横截面积 A、原长 L0。 |
描述了材料在弹性范围内,内部“应力流”与“应变流”之间的线性传导关系。杨氏模量 E是这个传导关系的比例系数或阻抗。应力 σ是驱动变形的“广义力”,应变 ϵ是产生的“广义位移”。胡克定律 σ=Eϵ类似于电路中的欧姆定律 V=IR, 其中 E类似于电阻 R, 表征材料抵抗弹性变形的“流阻”。这是一种最简单的、线性的“力流”-“变形流”本构关系。 |
|
TH-D1-0010 |
材料力学 |
梁的弯曲 |
平面弯曲的正应力公式(弯曲公式) |
梁在纯弯曲时,横截面上任一点的正应力 σ与该点到中性轴的距离 y成正比,与截面惯性矩 I成反比。公式为:σ=−IMy。 |
1. 基本假设:平面假设(变形后横截面保持平面)、单向受力假设(纵向纤维间无挤压)、材料线弹性(服从胡克定律)。 |
条件:纯弯曲(横截面上只有弯矩,无剪力);材料线弹性;小变形;平面弯曲(弯矩作用在纵向对称面内)。对于横力弯曲(有剪力),该公式近似成立,误差在工程允许范围内。 |
几何学、平衡条件、胡克定律 |
场景:各种梁式构件(如桥梁、楼板、机床主轴)的强度校核与设计。计算梁在弯矩作用下,截面上最大拉应力和最大压应力。 |
σ:横截面上某点的正应力 (Pa), 拉为正,压为负。 |
线性分布:应力沿截面高度线性分布,中性轴上为零。 |
1. 求支座反力, 画出梁的弯矩图,确定危险截面及其弯矩 M。 |
描述了弯曲弯矩 M如何通过截面几何属性 I转化为横截面上线性分布的“正应力流”。弯矩 M是截面上的“内力矩流”,它需要由截面上的微内力 σdA对中性轴的矩来平衡。公式 σ=−My/I表明,这个“内力矩流”被“分配”到截面的各个微面积上,分配的比例因子就是 y/I。惯性矩 I是截面抵抗弯曲变形能力的度量,I越大,相同的弯矩 M产生的最大应力 σmax越小。这体现了通过合理设计截面形状(增大 I)来“疏导”和“降低”弯曲应力“流”强度的思想。 |
|
TH-D1-0011 |
材料力学 |
梁的变形 |
梁的挠曲线近似微分方程 |
在小变形和线弹性条件下,梁的挠度 w(x)与弯矩 M(x)满足关系:dx2d2w=EIM(x)或 EIw′′=M(x)。 |
1. 曲率公式:由纯弯曲公式,曲率 κ=ρ1=EIM。 |
w'' |
}{(1+(w')^2)^{3/2}},其中w' = dw/dx为转角。<br>∗∗3.小变形近似∗∗:当梁的挠度很小,即w' << 1(转角远小于1弧度)时,(1+(w')^2)^{3/2} \approx 1。因此,曲率可近似为\kappa \approx |
w'' |
。<br>∗∗4.符号约定∗∗:根据弯矩正负和挠曲线凹凸方向(w''的正负)的协调,通常规定:当弯矩M为正(使梁下侧受拉)时,挠曲线向下凸,w''为负。因此,引入符号得到:w'' = -\frac |
|
编号 |
类别 |
领域 |
定理/定律/原理名称 |
核心表述(公式/叙述) |
定理的逐步思考推理过程、证明思路及关键数学步骤 |
成立条件与适用范围 |
底层规律/理论归属 |
典型应用场景和工程意义 |
变量/参数列表及说明 |
数学与逻辑特征 |
定理的演绎/应用流程步骤 |
物理本质与“流”描述 |
|---|---|---|---|---|---|---|---|---|---|---|---|---|
|
TH-D1-0021 |
系统动力学 |
分析力学 |
达朗贝尔原理 |
在质点运动的任一瞬时,作用于质点的主动力 F、约束反力 N和质点的惯性力 FI=−ma在形式上组成平衡力系。即:F+N+FI=0。 |
1. 牛顿第二定律:F+N=ma。 |
条件:适用于惯性参考系。惯性力是虚拟力,没有施力物体,不服从牛顿第三定律。 |
牛顿第二定律、静力平衡原理 |
场景:求解非自由质点系的加速度和约束反力(如滑轮系统、旋转机械的轴承动反力)、转子动平衡分析。 |
F:主动力(真实存在的力)。 |
动静法:将动力学问题转化为静力学问题求解。 |
1. 受力分析:画出研究对象上的所有主动力 F和约束反力 N。 |
描述了在惯性参考系中,真实力与虚拟的惯性力共同构成一种“动态平衡”。惯性力 −ma代表了质点由于惯性而抵抗运动状态改变的“趋势”或“效应流”。达朗贝尔原理本质是将动量变化率 ma 从方程的一侧“搬运”到另一侧,并赋予其“力”的形式,从而将动力学的“因果流”(力产生加速度)重新表述为静力学的“平衡流”(所有力的矢量和为零)。这为理解非惯性系中的“表观力”提供了桥梁。 |
|
TH-D1-0022 |
系统动力学 |
分析力学 |
虚位移原理(虚功原理) |
对于具有理想约束的质点系,其平衡的充分必要条件是:所有主动力在任何虚位移上所作的元功之和等于零。即:∑Fi⋅δri=0。 |
1. 基本概念: |
条件:系统受理想约束;考虑虚位移,而非实位移;适用于静平衡问题。 |
能量守恒、约束力学 |
场景:复杂机构(如连杆机构、差动滑轮)的平衡位置或平衡力求解、结构力学中利用虚功原理求位移(单位载荷法)。 |
Fi:第 i个质点所受的主动力。 |
标量方程:将矢量平衡问题转化为一个标量方程。 |
1. 判断约束:确认系统所受约束为理想约束。 |
描述了系统平衡时,所有主动力在任意可能(虚拟)的位移模式上所做的“功流”之和必须为零。虚位移 δr代表了系统在约束允许下的所有可能运动趋势。虚功原理要求,在这些所有可能的“运动趋势流”中,主动力不能在其中任何一个趋势上做净正功,否则系统就会沿该趋势运动而打破平衡。这本质上是能量原理在约束系统平衡问题上的体现,即系统在平衡位置时,其势能(对于保守系统)取驻值。 |
|
TH-D1-0023 |
系统动力学 |
分析力学 |
拉格朗日方程 |
对于具有完整、理想约束的系统,其运动微分方程可由下式给出: |
1. 从达朗贝尔-拉格朗日原理出发:结合达朗贝尔原理和虚位移原理,得到动力学普遍方程:∑i(Fi−miai)⋅δri=0。 |
条件:系统受完整约束(约束方程不含速度项,或可积分);理想约束;广义坐标独立。 |
达朗贝尔原理、虚功原理、变分法 |
场景:机器人动力学建模、多体系统仿真、航天器姿态动力学、复杂机械系统(如汽车悬架)的振动分析。 |
qj:第 j个广义坐标。 |
二阶常微分方程组:对每个广义坐标得到一个二阶ODE。 |
1. 确定系统:分析自由度,选取合适的广义坐标 q1,q2,...,qn。 |
是能量流在广义坐标空间中的动力学表述。方程 dtd(∂q˙j∂L)−∂qj∂L=Qj∗可以理解为: |
|
TH-D1-0024 |
系统动力学 |
分析力学 |
哈密顿原理 |
在相同的时间间隔 [t1,t2]内,具有相同起止位形的所有可能运动中,真实运动是使哈密顿作用量 S取驻值(通常为极小值)的那一条。即: |
1. 作用量定义:S=∫t1t2Ldt, 它是一个泛函,依赖于路径 q(t)。 |
{t_1}^{t_2} - \int{t_1}^{t_2} \frac{d}{dt}\left(\frac{\partial L}{\partial \dot{q}}\right) \delta q dt。由于端点变分为零,第一项消失。<br>∗∗4.得到欧拉−拉格朗日方程∗∗:代入得\delta S = \int_{t_1}^{t_2} \left[ \frac{\partial L}{\partial q} - \frac{d}{dt}\left(\frac{\partial L}{\partial \dot{q}}\right) \right] \delta q dt = 0。由于\delta q的任意性,被积函数必须为零,即得到拉格朗日方程。<br>∗∗5.原理陈述∗∗:真实运动使S$ 取驻值。 |
条件:系统受完整、理想约束;作用量积分存在;起止位形固定。 |
变分法、最小作用量原理 |
场景:理论推导连续介质力学方程(如梁、板的振动方程)、最优控制理论(庞特里亚金最小值原理的源头)、量子力学的路径积分表述。 |
S:哈密顿作用量,是一个泛函。 |
变分原理:用泛函的驻值条件代替微分方程。 |
1. 确定系统:明确系统的动能 T和势能 V, 构造拉格朗日函数 L=T−V。 |
|
TH-D1-0025 |
系统动力学 |
振动理论 |
单自由度系统无阻尼自由振动方程 |
质量-弹簧系统的运动由以下微分方程描述: |
1. 建立模型:考虑一个质量为 m的物块,连接一个刚度为 k的弹簧,置于光滑平面。以静平衡位置为坐标原点。 |
条件:系统是线性的(弹簧力与位移成正比);无阻尼;无外部激励;小位移(保证线性)。 |
牛顿第二定律、胡克定律 |
场景:任何近似为质量-弹簧系统的振动分析:机床基础的隔振、车辆悬架的垂向振动、精密仪器平台的防振设计。 |
x(t):质量块相对于静平衡位置的位移。 |
二阶线性齐次ODE:标准形式 x¨+ωn2x=0。 |
1. 建模:将实际系统简化为质量-弹簧模型,确定等效 m和 k。 |
描述了动能与势能之间无耗散的周期性流动与转换。质量 m代表惯性,是“动能储存器”;弹簧 k代表弹性,是“势能储存器”。运动方程 mx¨+kx=0表明,加速度(动量的变化率)由位移(势能的度量)驱动。在振动过程中,能量在动能 21mx˙2和势能 21kx2之间来回流动,总量守恒。固有频率 ωn=k/m是这种能量流动节奏的固有速率,由系统的“惯性流阻”和“弹性容抗”共同决定。 |
|
TH-D1-0026 |
系统动力学 |
振动理论 |
单自由度系统有阻尼自由振动方程 |
考虑粘性阻尼的质量-弹簧-阻尼器系统的运动方程为: |
1. 建立模型:在质量-弹簧系统上增加一个粘性阻尼器,阻尼力 Fd=−cx˙, 与速度成正比,方向相反。 |
条件:线性系统;阻尼为粘性阻尼(与速度成正比)。 |
牛顿第二定律、粘性阻尼模型 |
场景:车辆减震器设计、建筑结构在地震后的余振分析、精密仪器需要快速停止晃动的场合(如天平、光学平台)。 |
c:粘性阻尼系数 (N·s/m)。 |
二阶线性齐次ODE:标准形式 x¨+2ζωnx˙+ωn2x=0。 |
1. 建模:确定 m,c,k。 |
描述了能量在动能、势能和热耗散之间的流动。阻尼器 c代表一种“能量耗散器”或“阻尼流阻”,它将机械能(动能)不可逆地转化为热能。运动方程 mx¨+cx˙+kx=0中,cx˙项代表了与速度成正比的耗散力。在振动过程中,系统的总机械能不断通过阻尼器“泄漏”出去,导致振幅指数衰减 e−ζωnt。阻尼比 ζ量化了这种“泄漏”速率相对于系统固有振荡速率 ωn的大小。临界阻尼 (ζ=1) 是“能量流”能够最快耗散而不发生往复振荡的临界状态。 |
|
TH-D1-0027 |
系统动力学 |
振动理论 |
单自由度系统受迫振动(简谐激励)方程 |
系统在简谐激励力 F(t)=F0cos(ωt)作用下的运动方程为: |
1. 建立方程:在阻尼自由振动方程右端加上激励项。 |
条件:线性系统;激励为简谐形式;系统稳定(ζ>0)。 |
线性系统理论、常微分方程 |
场景:旋转机械(如电机、涡轮机)由于不平衡质量引起的振动、车辆在周期性不平路面上的振动、建筑物在风载荷或地震波作用下的响应分析。 |
F0:激励力的幅值 (N)。 |
非齐次线性ODE:标准形式。 |
1. 写出方程:mx¨+cx˙+kx=F0cos(ωt)。 |
描述了外部能量通过激励力持续输入系统,与系统内部的惯性、阻尼和弹性特性相互作用,形成稳态振动的过程。激励力 F0cos(ωt)是一个外部“能量流”源。系统的响应振幅 X和相位 ϕ反映了系统对这个输入“能量流”的“阻抗”特性。频率响应函数 X/(F0/k)类似于电路的增益。共振发生时 (ω≈ωn), 系统的动态阻抗最小,外部能量输入与系统内部能量存储(动能与势能)的节奏同步,导致能量在系统内高效积累(振幅增大)。阻尼则扮演了“能量泄放通道”的角色,限制了共振时的能量积累。相位差 ϕ反映了能量流动的时序关系。 |
|
TH-D1-0028 |
材料力学 |
强度理论 |
最大拉应力理论(第一强度理论) |
引起材料脆性断裂的主要因素是最大拉应力。无论材料处于何种应力状态,只要其最大主应力 σ1达到材料单向拉伸时的强度极限 σb, 材料即发生断裂。强度条件为:σ1≤[σ], 其中 [σ]=σb/nb, nb为安全系数。 |
1. 基本假设:材料的破坏(断裂)是由最大拉应力引起的。适用于脆性材料(如铸铁、陶瓷、岩石)。 |
条件:适用于脆性材料,且最大主应力为拉应力的情况。材料通常拉压强度不同。 |
材料破坏机理、主应力分析 |
场景:铸铁构件(如机床床身、管道)的拉伸强度设计、岩石和混凝土结构的抗拉强度评估。 |
σ1,σ2,σ3:一点处的三个主应力,且 σ1≥σ2≥σ3。 |
单参数准则:只依赖于最大主应力 σ1。 |
1. 应力分析:对危险点进行应力分析,求出其三个主应力 σ1,σ2,σ3。 |
基于一种简单的“力流”中断模型:认为材料内部的结合力(分子、原子间作用力)主要抵抗拉应力。当某点的最大拉应力 σ1超过材料固有的结合强度 σb时,该处的“结合键流”被拉断,裂纹萌生并扩展,导致宏观断裂。该理论将复杂的多向应力状态“映射”到最简单的单向拉伸试验上,认为破坏是由单一方向的“拉应力流”过载引起的。 |
|
TH-D1-0029 |
材料力学 |
强度理论 |
最大拉应变理论(第二强度理论) |
引起材料脆性断裂的主要因素是最大拉应变。无论材料处于何种应力状态,只要其最大拉应变 ϵ1达到材料单向拉伸时的极限拉应变 ϵu(ϵu=σb/E), 材料即发生断裂。强度条件为:σ1−ν(σ2+σ3)≤[σ]。 |
1. 基本假设:材料的破坏是由最大拉应变引起的。也主要针对脆性材料。 |
条件:适用于脆性材料;材料服从广义胡克定律直到破坏;通常拉压强度不同。 |
广义胡克定律、应变能 |
场景:历史上用于砖石、混凝土结构的强度分析。在某些特定组合应力状态下,对脆性材料破坏的预测可能比第一理论更准确。 |
ϵ1:最大主应变(通常对应最大主应力方向)。 |
多参数准则:依赖于三个主应力,通过泊松比耦合。 |
1. 应力分析:求出危险点的三个主应力。 |
从“变形流”的角度解释破坏:认为材料的破坏是由于某个方向的伸长变形(拉应变)超过了材料所能承受的极限。广义胡克定律 $\epsilon_1 = [\sigma |
(定理编号:TH-D1-0040 至 TH-D1-0060)。
|
编号 |
类别 |
领域 |
定理/定律/原理名称 |
核心表述(公式/叙述) |
定理的逐步思考推理过程、证明思路及关键数学步骤 |
成立条件与适用范围 |
底层规律/理论归属 |
典型应用场景和工程意义 |
变量/参数列表及说明 |
数学与逻辑特征 |
定理的演绎/应用流程步骤 |
物理本质与“流”描述 |
|---|---|---|---|---|---|---|---|---|---|---|---|---|
|
TH-D1-0040 |
系统动力学 |
振动理论 |
多自由度系统无阻尼自由振动方程 |
一个具有 n个自由度的线性系统,其无阻尼自由振动由以下矩阵方程描述: |
1. 系统建模:用 n个广义坐标 x1,x2,...,xn描述系统位形。 |
条件:线性系统;无阻尼;无外激励;小变形。 |
拉格朗日方程、线性系统理论 |
场景:建筑结构的地震模态分析、汽车车身的振动特性、飞机机翼的颤振初步分析、多转子轴承系统的临界转速计算。 |
{x(t)}:n维位移列向量。 |
矩阵微分方程:用矩阵形式表示的耦合ODE组。 |
1. 建立模型:确定自由度,建立 [M]和 [K]矩阵(可用拉格朗日法、影响系数法等)。 |
描述了多自由度系统中,振动能量在各模态间可能的分布与流动。质量矩阵 [M]和刚度矩阵 [K]分别描述了系统的惯性分布和弹性耦合。运动方程 [M]{x¨}+[K]{x}={0}表明,加速度 {x¨}与位移 {x}通过 [M]和 [K]相互耦合。求解特征值问题,就是寻找一组特殊的同步运动模式(振型 {ϕ}r),在这些模式下,所有坐标以同一频率 ωr做简谐振动,且惯性力 −ωr2[M]{ϕ}r与弹性恢复力 [K]{ϕ}r处处平衡。每个振型对应一种系统内能量储存与转换的固有空间分布模式。 |
|
TH-D1-0041 |
系统动力学 |
振动理论 |
模态叠加法 |
对于线性阻尼系统,其动力响应可以通过系统的各阶模态振型的线性叠加来表示。即: |
1. 坐标变换:将物理坐标 {x}变换到模态坐标 {q}, 即 {x}=[Φ]{q}。 这里 [Φ]由无阻尼系统的模态振型构成。 |
条件:线性系统;阻尼矩阵可对角化(通常满足比例阻尼假设);已求得无阻尼系统的模态。 |
线性代数(坐标变换、矩阵对角化)、单自由度系统理论 |
场景:结构在地震载荷、风载荷、冲击载荷下的时程响应分析;机械系统在复杂激励下的振动计算;响应谱分析。 |
{q(t)}:模态坐标向量。 |
坐标变换:通过线性变换实现解耦。 |
1. 模态分析:求解无阻尼系统特征值问题,得到 ωr和 {ϕ}r, 并质量归一化。 |
提供了一种在“模态空间”中分析系统动力响应的“投影-叠加”流。物理坐标下的复杂运动 {x(t)}, 可以看作是在各个“固有运动模式通道”(模态振型 {ϕ}r)上同时发生的、独立运动的叠加。每个模态通道的“强度”由模态坐标 qr(t)描述。模态叠加法的精髓在于,通过模态变换,将物理空间中相互耦合的“能量流”通道,解耦为一组独立的、单方向的“模态能量流”通道。每个通道有自己的频率 ωr和阻尼 ζr, 外部激励 {F(t)}也被投影到各通道上成为模态力 fr(t)。系统的总响应就是所有这些独立通道响应的“矢量合成流”。 |
|
TH-D1-0042 |
系统动力学 |
振动理论 |
动力吸振器原理 |
在原主系统(质量 m1, 刚度 k1)上附加一个较小的子系统(质量 m2, 刚度 k2, 阻尼 c2), 构成一个两自由度系统。通过合理设计 m2, k2, c2, 可以使主系统在特定激励频率下的振动幅值显著减小,甚至为零。 |
1. 系统模型:主系统受简谐力 F0eiωt激励。附加系统与主系统连接。 |
条件:主系统和吸振器均为线性;激励为简谐或窄带;主系统可近似为单自由度系统。 |
两自由度系统强迫振动理论 |
场景:汽车发动机悬置系统、直升机旋翼的振动控制、精密光学平台隔振、输电线防舞动。 |
m1,k1:主系统的质量和刚度。 |
两自由度系统:系统模型简单,但现象丰富。 |
1. 识别干扰频率 ω:确定需要抑制的主要振动频率。 |
体现了通过引入一个辅助的“能量交换器”,将主系统的振动能量“转移”和“耗散”掉的思想。主系统在外激励下获得能量。附加的吸振器被调谐到激励频率,与主系统发生强烈的相互作用。在理想情况下(无阻尼),吸振器像一个“动力反作用力发生器”,其运动产生的惯性力与弹性力,恰好与外激励力大小相等、方向相反,从而在主系统上完全抵消了激励的作用,使主系统保持静止。外激励输入的能量全部流入吸振器,并在其内部循环。加入阻尼后,这部分能量被吸振器的阻尼器耗散掉,从而更鲁棒地抑制了振动。这是一个典型的“能量转移与耗散”流控制策略。 |
|
TH-D1-0043 |
材料力学 |
应力分析 |
剪切胡克定律 |
在材料的线弹性范围内,剪应力 τ与剪应变 γ成正比。即:τ=Gγ, 其中 G为剪切模量(或切变模量)。 |
1. 实验基础:纯剪切试验(如薄壁圆管扭转)表明,在小变形下,剪应力与剪应变成正比。 |
条件:材料处于线弹性范围;小变形;通常适用于各向同性材料。 |
线弹性理论、连续介质力学 |
场景:圆轴扭转时的应力应变计算、螺栓和铆钉的剪切强度校核、弹性联轴器的扭转变形分析。 |
τ:剪应力 (Pa)。 |
线性本构关系:与正应力的胡克定律形式对称。 |
1. 确定受力状态:判断构件是否处于纯剪切或主要受剪状态。 |
描述了材料抵抗形状畸变(剪切变形)的“刚度”特性。剪切模量 G是材料抵抗“角度变形流”的“流阻”。剪应力 τ是驱动这种角度变形的“广义力”,剪应变 γ是产生的“广义位移”。定律 τ=Gγ表明,要产生一定的角度畸变,所需的剪应力与 G成正比。它是材料内部抵抗晶格滑移或分子链相对转动能力的宏观度量。 |
|
TH-D1-0044 |
材料力学 |
梁的弯曲 |
梁的剪力-弯矩微分关系 |
在分布载荷 q(x)作用的直梁上,剪力 Fs(x)、弯矩 M(x)和载荷集度之间存在以下微分关系: |
1. 微段平衡:从梁中取一微段 dx, 其左右截面内力分别为 Fs,M和 Fs+dFs,M+dM, 其上作用有分布载荷 q(x)(向上为正)。 |
条件:直梁;小变形;平面弯曲;载荷垂直于轴线。 |
微元体平衡、静力学 |
场景:快速绘制梁的内力图(剪力图、弯矩图);根据已知载荷推断内力分布;检验内力图正确性;推导梁的挠曲线微分方程。 |
Fs(x):横截面上的剪力,符号规定:使微段顺时针转为正。 |
微分关系:内力是载荷的积分,载荷是内力的微分。 |
1. 求支座反力。 |
揭示了梁的内力(剪力和弯矩)沿轴线方向的“流动”或“传递”规律。方程 dxdFs=−q(x)是剪力的连续性方程:剪力的变化率等于负的载荷集度。这意味着,分布载荷是“吸收”或“供给”剪力流的“源”或“汇”。方程 dxdM=Fs是弯矩的平衡方程:弯矩的变化率等于剪力。这表明,剪力是产生弯矩变化的“驱动力”,类似于剪力是弯矩的“流”。这两个微分关系共同描述了外部载荷 q(x)如何通过截面剪力的“传导”,最终体现为弯矩的“分布”。 |
|
TH-D1-0045 |
材料力学 |
组合变形 |
斜弯曲应力公式 |
当横向载荷不作用在梁的纵向对称面内时,梁发生斜弯曲。横截面上任一点的正应力为: |
1. 载荷分解:将横向力 F沿截面的两个形心主轴方向 y,z分解,产生弯矩 Mz=Fyx和 My=−Fzx。 |
条件:线弹性材料;小变形;平面弯曲的叠加原理适用;载荷垂直于轴线,但不通过弯曲中心(对开口薄壁截面需注意约束扭转)。 |
叠加原理、平面弯曲理论 |
场景:屋面桁条、桥梁的桥面板、机床导轨等在非对称载荷下的强度计算。 |
My,Mz:弯矩在形心主轴 y和 z方向的分量 (N·m)。 注意符号:My使 z坐标为正值处产生拉应力时为正;Mz使 y坐标为正值处产生拉应力时为正。 |
线性叠加:两个平面弯曲应力的代数和。 |
1. 建立形心主轴坐标系 Oyz。 |
描述了当弯曲载荷不在主轴平面时,弯矩“流”被分解到两个正交的主方向,分别产生应力“流”,然后叠加的机制。合弯矩 M是一个矢量。在主轴坐标系下,它可以分解为 My和 Mz两个分量。每个分量作为一个独立的“弯矩流”,按照各自方向的平面弯曲规律(σ∝M/I和到中性轴的距离)产生正应力分布。总的正应力场是这两个独立应力场的线性叠加。这体现了矢量分解与线性叠加在力学分析中的强大应用。 |
|
TH-D1-0046 |
材料力学 |
组合变形 |
偏心拉伸/压缩应力公式 |
当轴向力 F不通过杆件横截面的形心时,产生偏心拉压。横截面上任一点的正应力为: |
1. 力系简化:将偏心力 F向截面形心简化,得到一个轴向力 F和两个弯矩 Mz=F⋅ey, My=F⋅ez。 |
条件:线弹性材料;小变形;截面至少有一个对称轴(否则可能引起弯曲)。 |
叠加原理、轴向拉压与弯曲理论 |
场景:厂房立柱、桥梁墩台、机械设备地脚螺栓的应力计算、土建结构中避免截面出现拉应力的设计(如水库重力坝)。 |
F:轴向力 (N), 拉力为正。 |
线性叠加:轴向应力与两个方向弯曲应力的叠加。 |
1. 确定内力:将偏心力向形心简化,得到轴力 F和弯矩 My, Mz。 |
描述了轴向力“流”和弯曲力矩“流”共同作用下的正应力合成。偏心力 F的作用可以分解为经过形心的均匀拉伸/压缩流 (F/A) 和绕主轴的弯曲流 (My,Mz)。总的正应力场是均匀场和两个线性分布场的叠加。中性轴 σ=0是总应力场中正应力“汇”与“源”平衡的界线。形心核心是这样一个区域:当压力作用在其中时,由弯曲产生的拉应力“流”处处小于由轴向压力产生的压应力“流”,从而保证整个截面上只有压应力“流”,这对于“抗拉能力弱”的材料意味着“力流”通道不会被拉应力中断。 |
|
TH-D1-0047 |
材料力学 |
能量法 |
卡氏第一定理 |
对于线弹性结构,若其应变能 U表示为载荷 Fi的函数,则应变能对任一载荷 Fi的偏导数,等于该载荷作用点沿其作用方向的位移 δi。即: |
1. 应变能:线弹性结构的应变能是载荷的二次齐次函数。例如,对拉伸杆 U=F2L/(2EA)。 |
条件:线弹性结构;小变形;应变能必须表示为载荷 Fi的显函数;位移 δi是 Fi作用点沿其方向的位移。 |
功能原理、线弹性理论 |
场景:计算静定或超静定结构在特定点、特定方向的位移。例如,求梁上某点的挠度或转角,求桁架某节点的位移。 |
U:结构的应变能 (J), 是外力 F1,F2,...,Fn的函数。 |
导数关系:位移是应变能对力的偏导数。 |
1. 写出应变能表达式:用所有外力 Fi表示结构的总应变能 U。例如,对由 m根杆组成的桁架,U=∑j=1m2EjAjNj2Lj, 其中 Nj是第 j根杆的轴力,需表示为 Fi的函数。 |
反映了应变能对广义力的变化率等于相应的广义位移。可以将应变能 U视为储存在结构中的“弹性势能场”。广义力 Fi是这个场的“广义坐标”。卡氏第一定理 δi=∂U/∂Fi表明,广义位移 δi是“势能场”在力 Fi方向上的“梯度”。换句话说,如果你轻微地推动某个载荷 Fi(改变其大小),系统储存的弹性势能的变化率,正好等于该载荷作用点已经发生的位移。这揭示了力与位移在能量意义上的共轭关系。 |
|
TH-D1-0048 |
材料力学 |
能量法 |
卡氏第二定理 |
对于线弹性结构,若其应变能 U表示为位移 δi的函数,则应变能对任一位移 δi的偏导数,等于该位移对应的力 Fi。即: |
1. 与第一定理的比较:卡氏第一定理中 U=U(Fi), 第二定理中 U=U(δi)。 由于线弹性结构力与位移成正比,二者存在一一对应的 Legendre 变换关系。 |
条件:线弹性结构;小变形;应变能必须表示为位移 δi的显函数。 |
余能原理、Legendre变换 |
场景:作为理论工具,用于推导结构力学中的其他方法(如单位载荷法、力法正则方程)。在分析以位移为基本未知量的有限元法中也有体现。 |
U:结构的应变能,是广义位移 δ1,δ2,...,δn的函数。 |
导数关系:力是应变能对位移的偏导数。 |
1. (理论上)将 U表示为位移的函数:这通常需要知道力-位移关系和结构的本构方程,往往很复杂。 |
与第一定理对偶,描述了应变能对广义位移的变化率等于相应的广义力。如果将结构视为一个“弹性系统”,其广义位移 δi是描述系统位形的“坐标”,那么应变能 U(δi)就是该系统的“势能函数”。卡氏第二定理 Fi=∂U/∂δi正是势能梯度等于力这一基本原理的表述。广义力 Fi是试图改变广义位移 δi的“广义推力”,其大小等于势能场在该方向的“坡度”。这为基于位移的有限元法提供了理论基础:单元刚度矩阵本质上就是应变能对节点位移求二阶导的结果。 |
|
TH-D1-0049 |
材料力学 |
能量法 |
莫尔定理(单位载荷法) |
计算线弹性结构在某点、沿某方向的位移 Δ的公式为: |
1. 虚功原理:在结构实际位移状态上,施加一组虚设的单位载荷,则外力虚功等于内力虚功。 |
条件:线弹性结构;小变形;单位载荷引起的结构响应与实际载荷引起的响应无关(适用叠加原理)。 |
虚功原理、线弹性本构关系 |
场景:计算梁或刚架指定点的挠度和转角;求桁架节点的位移;校核结构的刚度;解超静定结构时的变形协调条件计算。 |
Δ:所求的广义位移(线位移或角位移)。 |
积分计算:位移表达为定积分。 |
1. 求实际内力:在实际载荷作用下,求出结构各段的内力方程 M(x),N(x)等。 |
是虚功原理在线弹性位移计算中的具体实现。其核心思想是:在真实的变形状态(由实际载荷引起)上,施加一个虚拟的“单位力流”,并计算这个“虚拟力流”在整个真实变形场上所做的“虚功”。根据虚功原理,这个“外力虚功”等于“内力虚功”。内力虚功是虚拟内力 N,M在真实微元变形 d(ΔL),dθ上做功的积分。由于真实变形由实际内力通过本构关系决定,最终公式表现为两个内力场乘积的积分。这本质上是一种“标量积”或“投影”:将实际的内力(变形)场“投影”到单位虚载荷引起的内力场上,其“投影长度”就是所求的位移 Δ。图乘法则是这个投影在特定条件下的几何直观体现。 |
|
TH-D1-0050 |
材料力学 |
超静定问题 |
力法(柔度法)正则方程 |
用于求解 n次超静定结构。基本思想是解除多余约束,代之以多余未知力 X1,X2,...,Xn, 根据原结构在解除约束处的位移协调条件建立方程: |
1. 确定超静定次数, 选择基本静定系(静定基),并解除相应的多余约束,用多余未知力 Xi代替。 |
条件:线弹性结构;小变形。 |
叠加原理、位移协调条件、线弹性理论 |
场景:求解多跨连续梁、超静定刚架、桁架的内力分布;计算复杂结构的支座反力;为后续的强度、刚度和稳定性分析提供基础。 |
Xi:第 i个多余未知力(广义力)。 |
线性方程组:关于多余未知力的线性方程组。 |
1. 判断超静定次数 n。 |
是位移协调条件在力空间中的表达。超静定结构的内力分布不仅满足平衡条件,还必须满足变形协调条件。力法将后者表达为关于多余未知力 Xi的方程。柔度系数 δij反映了结构在单位力 Xj作用下,在 i方向产生的“位移柔量”,它描述了结构“力流”通道的“柔顺性”。正则方程 ∑δijXj=−ΔiP要求,所有多余力 Xj在其自身作用点引起的位移总和,必须恰好抵消外载荷在该点引起的位移,以满足原结构的约束条件。这相当于在结构的“变形流场”中,用一组内部的“自平衡力流” Xi来“修正”由外载荷引起的变形,使其与外部约束相容。 |
|
TH-D1-0051 |
材料力学 |
压杆稳定 |
欧拉公式(理想压杆临界力) |
对于一端固定、另一端自由的细长理想压杆,其临界压力为: |
1. 理想模型:等截面直杆;材料线弹性;压力严格沿轴线;初始无弯曲;小变形。 |
条件:理想压杆(无初始缺陷);线弹性材料;小变形;细长杆(满足欧拉公式适用条件 λ≥λP)。 |
梁的弯曲微分方程、特征值问题 |
场景:钢结构立柱、桥梁桁架中的压杆、机械工程中的顶杆、活塞杆的稳定性设计。 |
Fcr:临界压力 (N), 超过此值,杆将失稳。 |
特征值问题:临界力是微分方程边值问题的特征值。 |
1. 判断杆端约束, 确定长度因数 μ。 |
描述了压杆在轴向压力下,从直线平衡形态向弯曲平衡形态“分岔”的临界条件。压力 F产生弯矩 −Fw, 而弯曲变形产生恢复弯矩 EIw′′。平衡方程 EIw′′+Fw=0表明,压力 F扮演了“负刚度”的角色。当 F较小时,正的弯曲刚度占优,直线形态稳定。当 F达到临界值 Fcr=π2EI/(μL)2时,压力产生的“负刚度”效应恰好抵消了弯曲刚度,系统处于随遇平衡的临界状态。超过 Fcr, “负刚度”占优,微小的扰动就会使弯曲变形迅速增大,导致失稳。欧拉公式给出了这个“刚度竞争”的平衡点。 |
|
TH-D1-0052 |
材料力学 |
动载荷 |
惯性力法(动静法)求动应力 |
对于做等加速直线运动或等角速转动的构件,可以通过施加惯性力,将动力学问题在形式上转化为静力学问题,然后按静载荷方法计算应力(动应力)。 |
1. 达朗贝尔原理应用:对构件中的每个质点,施加惯性力 −aidm。 |
条件:构件作等加速直线运动或等角速度转动;材料仍处于线弹性范围;不考虑由加速度变化引起的弹性振动。 |
达朗贝尔原理、静力学方法 |
场景:起重机加速吊装时钢丝绳的应力计算、旋转飞轮或涡轮机转子的应力分析、火箭发射时部件的过载分析。 |
a:线加速度 (m/s²)。 |
动静法转化:通过加惯性力化为静力问题。 |
1. 运动分析:确定构件各点的加速度。 |
是达朗贝尔原理在动应力计算中的直接应用。将惯性力 −ρa视为作用在构件上的“虚拟体积力场”,这个力场与真实外力场和内部应力场共同构成“瞬时静力平衡”。这样,就将由加速度引起的、分布在整个构件上的“动量变化效应流”,转化为一个虚构的、静态的“体积力流”来处理。动应力就是真实外力与这个虚构惯性力共同作用产生的应力。动荷因数 Kd衡量了惯性力“流”相对于静重力“流”的强度。 |
|
TH-D1-0053 |
材料力学 |
动载荷 |
冲击动荷因数 |
当重物以速度 v垂直冲击被冲击构件时,冲击动荷因数 Kd近似为: |
1. 基本假设:冲击物为刚体;被冲击物为线弹性体(质量忽略或等效);冲击过程无能量损失(能量守恒);冲击后两者一起运动。 |
条件:冲击物为刚性;被冲击物为线弹性;忽略冲击过程中的其他能量损失(如热、声);被冲击物质量可忽略,或冲击持续时间远大于被冲击物的自振周期。 |
机械能守恒定律、线弹性假设 |
场景:打桩机效率分析、冲压加工中模具的强度校核、汽车碰撞时保险杠的吸能设计、重物坠落对平台或地板的冲击评估。 |
Q:冲击物的重量 (N)。 |
能量法:基于能量守恒推导。 |
1. 计算静位移 Δst:在冲击点,沿冲击方向,施加一个大小等于冲击物重量 Q的静力,计算该点的位移。 |
描述了冲击过程中,冲击物的动能或势能突然转化为被冲击物的应变能,导致瞬时载荷和应力急剧放大的现象。Kd是这个能量转换过程中“冲击效应放大系数”。静位移 Δst反映了被冲击系统的“柔度”或“吸能潜力”。Δst越大,系统在冲击时能通过更大的变形来吸收能量,从而降低冲击力峰值,即“以柔克刚”。能量守恒方程 Q(h+Δd)=21kΔd2是“能量流入”等于“能量储存”的平衡。冲击问题的核心是能量的快速流动与再分配。 |
|
TH-D1-0054 |
材料力学 |
交变应力 |
疲劳极限(持久极限) |
在交变应力作用下,材料经过无限次应力循环而不发生疲劳破坏的最大应力值,称为疲劳极限或持久极限,记为 σr。下标 r表示循环特征,如对称循环 (r=−1) 的疲劳 |
(定理编号:TH-D1-0054 至 TH-D1-0074)。
|
编号 |
类别 |
领域 |
定理/定律/原理名称 |
核心表述(公式/叙述) |
定理的逐步思考推理过程、证明思路及关键数学步骤 |
成立条件与适用范围 |
底层规律/理论归属 |
典型应用场景和工程意义 |
变量/参数列表及说明 |
数学与逻辑特征 |
定理的演绎/应用流程步骤 |
物理本质与“流”描述 |
|---|---|---|---|---|---|---|---|---|---|---|---|---|
|
TH-D1-0054 |
材料力学 |
交变应力 |
疲劳极限(持久极限) |
在交变应力作用下,材料经过无限次应力循环而不发生疲劳破坏的最大应力值,称为疲劳极限或持久极限,记为 σr。下标 r表示循环特征,如对称循环 (r=−1) 的疲劳极限记为 σ−1。对于钢铁材料,通常以 107次循环未破坏对应的应力作为疲劳极限。 |
1. 疲劳试验:通过旋转弯曲疲劳试验机或拉压疲劳试验机,对一组标准光滑试件施加不同应力幅 σa的对称循环应力,记录导致破坏的循环次数 N。 |
条件:标准光滑小试件;恒幅对称循环应力;室温、空气环境。 |
材料科学、损伤累积理论 |
场景:轴、齿轮、弹簧、连杆、桥梁、飞机结构等承受交变载荷的零部件设计。 |
σr:循环特征为 r时的疲劳极限 (Pa)。 |
经验性/实验性:由试验测定,无严格数学推导。 |
1. 制备标准试件。 |
代表了材料抵抗在交变应力下裂纹萌生和早期扩展的“内在抗力阈值”。当外部交变应力“流”的幅值低于这个阈值时,其能量不足以驱动材料内部缺陷(如位错、夹杂物)处发生不可逆的滑移和微裂纹形成,或者不足以使已有微裂纹克服扩展阻力。因此,能量“流”在材料内部循环往复,但不会造成损伤的累积和宏观破坏。疲劳极限是材料微观结构(晶粒、相界、缺陷)对循环滑移阻碍能力的宏观体现。对于钢铁材料,其明显的疲劳极限可能与循环硬化及裂纹尖端的闭合效应有关。 |
|
TH-D1-0055 |
材料力学 |
交变应力 |
疲劳寿命的 Miner 线性累积损伤理论 |
在变幅交变应力作用下,疲劳损伤可以线性累加。若在应力水平 Si下循环 ni次,该应力水平下的疲劳寿命为 Ni(由 S-N 曲线确定),则损伤率 Di=ni/Ni。当各应力水平造成的损伤率之和等于 1 时,发生疲劳破坏。即: |
1. 基本假设: |
条件:变幅载荷;各应力水平下的损伤可线性叠加;忽略载荷间的相互作用(如过载引起的残余应力、裂纹闭合等)。 |
损伤力学、经验假设 |
场景:汽车底盘、飞机机翼、风力发电机叶片等在随机载荷下的寿命预估;机械零件的定期检修周期制定。 |
Si:第 i个应力水平(通常指应力幅)。 |
线性假设:损伤与循环次数成正比,且可线性叠加。 |
1. 获取载荷谱:通过实测或标准,得到构件危险点的应力时间历程。 |
将疲劳破坏视为一种“损伤量”的线性累积过程。每个应力循环都像是一滴“损伤水滴”,滴入一个“损伤桶”中。不同应力水平下的循环,其“水滴”大小不同(1/Ni)。Miner 理论假设这些“水滴”可以简单地累加,且“桶”的容量是固定的(D=1)。当累积的“损伤流体”达到桶的容量时,破坏发生。这是一种高度简化的模型,它忽略了损伤累积可能存在的非线性(例如,大应力“水滴”可能会改变“桶”的容量或后续“水滴”的累积效率),但为处理复杂的载荷“流”提供了一条实用的量化路径。 |
|
TH-D1-0056 |
材料力学 |
强度理论 |
最大拉应力理论(第一强度理论) |
认为材料的破坏是由最大拉应力引起的。当构件危险点处的最大拉应力 σ1达到材料单向拉伸时的极限应力 σb时,发生脆性断裂。强度条件为: |
1. 破坏假设:无论应力状态如何,只要最大主应力 σ1达到材料极限,即发生断裂。 |
条件:适用于脆性材料(如铸铁、混凝土、岩石);且危险点处于拉应力状态(通常 σ1>0)。 |
经验性破坏准则 |
场景:铸铁构件受拉或受弯时的强度校核;混凝土梁受拉区开裂判断;地质岩层拉伸断裂分析。 |
σ1,σ2,σ3:主应力,且 σ1≥σ2≥σ3。 |
单参数准则:只依赖于最大主应力 σ1。 |
1. 应力分析:确定构件危险点的应力状态,计算其主应力 σ1,σ2,σ3。 |
认为材料的破坏是由“拉应力流”通道被拉断所致。该理论假设材料内部存在微观缺陷或薄弱面,当最大主应力方向(σ1方向)的拉应力“流”强度超过材料内部结合力所能承受的极限时,缺陷会扩展并连接,导致宏观上的脆性断裂。它忽略了压应力“流”可能对缺陷闭合产生的有利影响,也忽略了剪切“流”的作用,是一种最直接的“拉断”模型。 |
|
TH-D1-0057 |
材料力学 |
强度理论 |
最大拉应变理论(第二强度理论) |
认为材料的破坏是由最大拉应变引起的。当构件危险点处的最大拉应变 ϵ1达到材料单向拉伸时的极限应变 ϵu时,发生破坏。利用广义胡克定律,强度条件可写为: |
1. 破坏假设:无论应力状态如何,最大拉应变是破坏主因。 |
条件:线弹性材料直到破坏;小变形;适用于部分脆性材料,特别是 σ2, σ3为压应力的情况。 |
经验性破坏准则、广义胡克定律 |
场景:早期用于铸铁构件在复杂应力状态下的强度计算,如受内压的脆性材料管道。 |
ϵ1:最大主应变。 |
三参数准则:依赖于三个主应力,但以组合形式出现。 |
1. 应力分析:计算危险点的主应力 σ1,σ2,σ3。 |
认为破坏是由“拉伸变形流”超过容限引起的。它不仅关注最大主应力方向的直接拉伸,还考虑了由于泊松效应,其他方向的应力(σ2, σ3)会间接影响该方向的应变。例如,压应力 σ2, σ3会产生横向膨胀,从而增加 σ1方向的拉应变。因此,该理论认为,即使 σ1相同,处于双向受压状态的点比单向受拉的点更不容易因拉应变而破坏,因为压应力抑制了横向膨胀。这是一种基于“变形协调”的破坏观点。 |
|
TH-D1-0058 |
材料力学 |
强度理论 |
最大切应力理论(第三强度理论,Tresca 准则) |
认为材料的破坏(屈服)是由最大切应力引起的。当构件危险点处的最大切应力 τmax达到材料单向拉伸屈服时的最大切应力 τs时,发生塑性屈服。强度条件为: |
1. 破坏假设:塑性屈服由最大切应力控制。 |
条件:适用于塑性材料(如大多数金属);材料屈服后即认为失效;通常用于静载荷下的强度设计。 |
塑性力学、剪切滑移机制 |
场景:轴类零件的强度设计(弯扭组合);压力容器筒体的强度校核;塑性材料构件在复杂应力状态下的屈服判断。 |
τmax:最大切应力。 |
双参数准则:只依赖于最大和最小主应力 σ1和 σ3, 忽略中间主应力 σ2的影响。 |
1. 应力分析:计算危险点的主应力 σ1,σ2,σ3。 |
认为材料的塑性屈服是“剪切滑移流”达到临界值的结果。该理论基于金属晶体学的位错滑移机制:当作用在某个滑移面上的切应力达到临界值时,位错开始大规模运动,导致宏观屈服。最大切应力 τmax代表了最有可能发生滑移的平面上的切应力驱动“流”。准则 σ1−σ3=σs表明,屈服只取决于最大主应力与最小主应力之差,即最大剪应力幅值,而与静水压力(σ1+σ2+σ3)/3)无关。这符合金属材料屈服对静水压力不敏感的观察。 |
|
TH-D1-0059 |
材料力学 |
强度理论 |
形状改变比能理论(第四强度理论,von Mises 准则) |
认为材料的破坏(屈服)是由形状改变比能(畸变能密度)引起的。当构件危险点处的形状改变比能 ud达到材料单向拉伸屈服时的形状改变比能 uds时,发生塑性屈服。强度条件为: |
1. 能量分解:应变比能 u可分解为体积改变比能 uv和形状改变比能 ud:u=uv+ud。 |
条件:适用于塑性材料;材料屈服后即认为失效;通常用于静载荷下的强度设计。 |
塑性力学、能量原理 |
场景:复杂应力状态下塑性材料构件的精确强度分析;有限元分析软件中默认的金属材料屈服准则;航空航天、汽车等高科技领域的关键部件设计。 |
ud:形状改变比能(畸变能密度)(J/m³)。 |
三参数准则:依赖于所有三个主应力之差,考虑了中间主应力的影响。 |
1. 应力分析:计算危险点的主应力 σ1,σ2,σ3。 |
认为材料的塑性屈服是“畸变能流”密度达到临界值的结果。该理论从能量角度出发,认为引起材料形状改变(畸变)的那部分弹性应变能,是驱动微观位错运动和宏观塑性变形的“能量流”。体积改变部分(静水压力)只改变材料密度,不引起形状变化,因此与屈服无关。准则的数学形式 J2=k(其中 J2是应力偏张量的第二不变量)意味着,屈服发生在应力偏张量的“模长”达到某一常数时,这与金属晶体多滑移系的宏观统计行为相符,能更好地拟合实验数据。 |
|
TH-D1-0060 |
材料力学 |
强度理论 |
莫尔强度理论 |
认为材料的破坏(断裂或屈服)主要是由于某个截面上的切应力达到极限值,但该极限值与该截面上的正应力有关。其强度条件可表示为: |
1. 基本思想:材料在不同应力状态下的破坏,是由于某一斜截面上的切应力 τ和正应力 σ共同作用的结果。破坏条件可写为 τ=f(σ), 函数 f由实验确定。 |
条件:适用于拉压强度不等的材料,如铸铁、混凝土、岩石等脆性材料。 |
经验性破坏准则、图解法 |
场景:铸铁构件在拉-压复合应力下的强度校核;混凝土结构在双向应力下的破坏分析;岩石地基的承载力计算。 |
σb:材料的抗拉强度极限。 |
\sigma_{bc} |
> \sigma_b,故[\sigma_c] > [\sigma_t]$。 |
双参数准则:依赖于 σ1和 σ3, 但系数由材料拉压强度比决定。 |
|
TH-D1-0061 |
系统动力学 |
分析力学 |
哈密顿原理 |
在相同时间内,保守系统从初始位形到最终位形的所有可能运动中,真实运动是使哈密顿作用量 S取驻值(通常是极小值)的那一条。即: |
1. 作用量定义:S=∫t1t2Ldt, 它是路径 q(t)的泛函。 |
条件:完整系统;对于保守系统,L=T−V;对于非保守系统,需已知广义力。变分路径满足端点固定。 |
变分法、最小作用量原理 |
场景:推导复杂离散系统或连续系统的运动方程;弹性力学中平衡方程和自然边界条件的统一推导;量子力学中费曼路径积分的经典对应。 |
S:哈密顿作用量,是泛函。 |
积分变分原理:用作用量积分 S的驻值条件代替微分方程。 |
1. 确定系统动能 T和势能 V, 构造拉格朗日函数 L=T−V。 |
是物理学中“最小作用量原理” 在力学中的表述。它断言:自然界的真实运动总是选择一条使某个称为“作用量”的积分量取极值(通常是极小值)的路径。作用量 S是拉格朗日函数 L(动能与势能之差)对时间的积分。哈密顿原理 δS=0意味着,在从起点到终点的所有“虚拟路径流”中,真实路径是使这个“作用量流”的总和取驻值的那一条。这就像光走最短时间路径一样,力学系统的运动也遵循某种“最经济”或“最平稳”的法则。它从全局和历史的视角,用一条简洁的数学原理概括了牛顿力学的所有运动定律。 |
|
TH-D1-0062 |
系统动力学 |
分析力学 |
哈密顿正则方程 |
对于保守系统,引入广义动量 pj=∂L/∂q˙j和哈密顿函数 H(q,p,t)=∑jpjq˙j−L, 则系统的运动由以下一对一阶方程描述: |
1. Legendre 变换:从以 (q,q˙)为自变量的拉格朗日函数 L, 变换到以 (q,p)为自变量的哈密顿函数 H:H(q,p,t)=pq˙−L(q,q˙,t), 其中 q˙要通过关系 p=∂L/∂q˙用 p表示。 |
条件:完整系统;变换 p=∂L/∂q˙可逆(即 Hessian 矩阵 ∂2L/∂q˙i∂q˙j非奇异);系统保守或广义力有势。 |
勒让德变换、辛几何 |
场景:天体力学中的轨道计算;分子动力学模拟;控制理论中的最优控制(哈密顿-雅可比-贝尔曼方程);量子力学的经典极限。 |
qj:广义坐标。 |
一阶方程组:将二阶的拉格朗日方程化为一阶方程组。 |
1. 写出拉格朗日函数 L(q,q˙,t)。 |
实现了从位形空间 (q)到相空间 (q,p)的转变,并将动力学表述为相空间中的“辛流”。正则方程 q˙=∂H/∂p,p˙=−∂H/∂q定义了一个相空间中的矢量场。这个矢量场由哈密顿函数 H的梯度决定,但其结构是“旋转”的:位置 q的变化率由动量 p方向的梯度给出,而动量的变化率由位置 q方向的负梯度给出。这种结构保证了相空间体积在流动下保持不变(刘维尔定理),并揭示了力学系统内在的辛几何特性。哈密顿函数 H是驱动这个“辛流”的“能量源”。 |
|
TH-D1-0063 |
系统动力学 |
振动理论 |
单自由度有阻尼系统自由振动解 |
系统方程:mx¨+cx˙+kx=0。 其解的形式取决于阻尼比 ζ=c/(2mk)=c/(2mωn)。 |
1. 特征方程:设解 x=est, 代入方程得特征方程 ms2+cs+k=0。 |
(定理编号:TH-D1-0064 至 TH-D1-0084)。
|
编号 |
类别 |
领域 |
定理/定律/原理名称 |
核心表述(公式/叙述) |
定理的逐步思考推理过程、证明思路及关键数学步骤 |
成立条件与适用范围 |
底层规律/理论归属 |
典型应用场景和工程意义 |
变量/参数列表及说明 |
数学与逻辑特征 |
定理的演绎/应用流程步骤 |
物理本质与“流”描述 |
|---|---|---|---|---|---|---|---|---|---|---|---|---|
|
TH-D1-0064 |
系统动力学 |
随机振动 |
随机过程的均方值、自相关函数与功率谱密度关系(Wiener-Khinchin定理) |
对于一个宽平稳随机过程 X(t),其自相关函数 RXX(τ)=E[X(t)X(t+τ)]与功率谱密度函数 SXX(ω)构成一对傅里叶变换对: |
1. 定义自相关函数:衡量随机信号自身在不同时刻的相关性,对于平稳过程,只依赖于时间差 τ。 |
条件:X(t)是宽平稳随机过程;自相关函数绝对可积 $\int_{-\infty}^{\infty} |
R_{XX}(\tau) |
d\tau < \infty$。 |
随机过程理论、傅里叶分析、信号处理 |
场景:分析风载荷、路面不平度、海浪、地震动等随机激励的频域特性;计算线性系统在随机激励下的输出功率谱密度和均方响应;结构随机疲劳分析。 |
X(t):平稳随机过程。 |
傅里叶变换对:SXX(ω)与 RXX(τ)互为傅里叶变换。 |
|
TH-D1-0065 |
系统动力学 |
随机振动 |
线性系统在平稳随机激励下的响应功率谱密度 |
对于一个线性时不变系统,其频率响应函数为 H(ω)。当输入是一个平稳随机过程 X(t), 其功率谱密度为 SXX(ω)时,输出 Y(t)也是一个平稳随机过程,且其功率谱密度 SYY(ω)和输入-输出互谱密度 SXY(ω)分别为: |
H(\omega) |
^2 S{XX}(\omega)<br>S{XY}(\omega) = H(\omega) S_{XX}(\omega)$ |
1. 系统输入输出关系:在频域,输出傅里叶变换 Y(ω)=H(ω)X(ω)。 |
Y_T(\omega) |
^2]}{2T}。其中Y_T(\omega)是y_T(t)(y(t)的截断)的傅里叶变换。代入Y_T(\omega)=H(\omega)X_T(\omega),得:<br>S{YY}(\omega) = \lim{T\to\infty} \frac{E[ |
H(\omega)X_T(\omega) |
^2]}{2T} = |
H(\omega) |
|
TH-D1-0066 |
系统动力学 |
非线性振动 |
相平面与奇点分类 |
对于二阶自治系统 x¨+f(x,x˙)=0, 令 y=x˙, 可化为一阶方程组: |
1. 相平面构建:以位移 x和速度 y=x˙为坐标轴,构成相平面。系统的一个状态对应相平面中的一个点,运动过程形成一条轨迹(相轨迹)。 |
J - \lambda I |
= 0得到特征根\lambda_{1,2}。根据特征根的类型判断奇点类型:<br>−∗∗结点∗∗:\lambda_1, \lambda_2为同号实数。稳定结点(负实部),不稳定结点(正实部)。<br>−∗∗焦点∗∗:\lambda_1, \lambda_2为共轭复数。稳定焦点(实部负),不稳定焦点(实部正)。<br>−∗∗中心∗∗:\lambda_1, \lambda_2为纯虚数。<br>−∗∗鞍点∗∗:\lambda_1, \lambda_2$ 为异号实数。 |
条件:系统是二阶、自治的(f不显含时间 t);在平衡点附近可线性化。 |
微分方程定性理论、动力系统 |
场景:分析单自由度非线性振子的稳定性(如 Duffing 振子、Van der Pol 振子);研究简单生态模型(捕食者-被捕食者);电路系统(如 LC 振荡器)的相图分析。 |
x:系统位移(或广义坐标)。 |
几何方法:用图形(相图)代替解析解。 |
|
TH-D1-0067 |
系统动力学 |
非线性振动 |
极限环 |
在相平面中,一个孤立的闭合轨迹称为极限环。它对应于非线性自治系统的周期解。极限环可以是稳定的、不稳定的或半稳定的。 |
1. 存在性:极限环是二阶(或更高维)非线性自治系统特有的现象,线性系统不可能有孤立的闭合轨迹(线性系统的闭合轨迹是连续族,如中心点周围的无穷多同心环)。 |
条件:非线性自治系统;通常是二阶系统,在高维系统中也可能存在极限环(周期轨道)。 |
动力系统理论、非线性微分方程 |
场景:机械系统中的自激振动,如机床颤振、输电线舞动、刹车尖叫;电子振荡器(如 Van der Pol 振荡器);生物节律(如心脏跳动、神经元放电);化学反应中的振荡现象。 |
极限环:相平面中的孤立闭合轨迹。 |
孤立闭轨:与线性系统的中心点周围连续闭轨族不同。 |
1. 建立系统模型:得到非线性自治微分方程。 |
代表了非线性系统中一种动态平衡的“闭合能量流”。在稳定极限环上,系统在一个周期内,由非线性机制提供的能量输入与由阻尼机制导致的能量耗散达到精确平衡。例如,在 Van der Pol 振子中,当振幅小时,非线性项提供负阻尼,相当于“泵入”能量,使振幅增大;当振幅大时,非线性项表现为正阻尼,“耗散”能量,使振幅减小。最终,系统状态被“吸引”到一个特定的闭合轨道上,形成稳定的周期振荡。这个闭合轨道就是极限环,它像一个吸引子,将附近的运动状态都吸引过来。极限环上的运动是自维持的,不需要外部周期激励,其振幅和频率由系统内部参数决定。 |
|
TH-D1-0068 |
系统动力学 |
分析力学 |
刘维尔定理 |
在哈密顿力学中,保守系统在相空间中的概率分布函数 ρ(q,p,t)沿着系统的运动轨迹是常数,即其物质导数为零: |
1. 相空间与系综:考虑大量性质完全相同、但初始条件各异的系统构成的系综。系综在相空间中的分布用概率密度函数 ρ(q,p,t)描述。 |
条件:系统是保守的,可由哈密顿量 H描述;相空间概率分布函数 ρ足够光滑。 |
哈密顿力学、统计力学 |
场景:统计力学中推导平衡态分布;分子动力学模拟中验证算法是否保辛(辛算法能更好地保持刘维尔定理);研究动力系统长期演化行为。 |
ρ(q,p,t):相空间中的概率密度函数。 |
守恒律:概率密度沿轨迹守恒。 |
1. 确定系统的哈密顿量 H(q,p)。 |
描述了保守哈密顿系统在相空间中“概率流体”演化的不可压缩性。想象相空间像一个充满“概率粒子”的流体,每个粒子代表系统的一个可能状态。哈密顿正则方程定义了该流体的速度场。刘维尔定理指出,这个速度场的散度为零,意味着流体是不可压缩的。因此,在流动过程中,任何一块相空间体积元虽然形状会扭曲,但其体积保持不变。同时,流体中“概率密度” ρ的“物质导数”为零,意味着跟随一个“概率粒子”运动,其周围的概率密度不会改变。这保证了相空间中的概率分布像刚体一样被“搬运”,而不会发生扩散或集中,是统计力学中各态历经假说和微正则分布成立的动力学基础。 |
|
TH-D1-0069 |
系统动力学 |
分析力学 |
泊松括号 |
对于相空间中的任意两个动力学变量 A(q,p,t)和 B(q,p,t), 其泊松括号定义为: |
1. 定义:泊松括号是一种双线性、反对称的微分运算,满足雅可比恒等式 {A,{B,C}}+{B,{C,A}}+{C,{A,B}}=0。 |
条件:系统可由正则变量 (q,p)描述;函数 A,B关于 (q,p)连续可微。 |
辛几何、李代数 |
场景:简洁表达哈密顿系统的运动方程;判断物理量是否为守恒量;研究动力系统的对称性(诺特定理);在摄动理论中研究近可积系统;量子力学中经典对应的建立。 |
A,B:相空间中的任意动力学函数(物理量)。 |
双线性:{aA+bB,C}=a{A,C}+b{B,C}。 |
1. 给定两个相空间函数 A(q,p,t)和 B(q,p,t)。 |
泊松括号 {A,B}度量了两个物理量 A和 B在相空间“流动”方向上的“不可对易性”或“相互影响”。它可以理解为:沿着由 B生成的哈密顿流移动,物理量 A的变化率。特别地,{A,H}就是 A沿着由系统哈密顿量 H生成的真实运动轨迹的变化率(若 A不显含时间)。基本关系 {qi,pj}=δij反映了坐标和其共轭动量之间的经典对易关系,是相空间基本辛结构的体现。泊松括号为零意味着两个量在运动中是“对易”的,其中一个量的守恒不会因另一个量的变化而改变。这种代数结构深刻地反映了力学系统的内在对称性和约束。 |
|
TH-D1-0070 |
系统动力学 |
连续系统振动 |
一维波动方程 |
描述一维均匀弹性介质(如弦、杆)中波传播的基本方程: |
1. 模型建立:以均匀弦的横向振动为例。取微段 dx, 设弦的线密度为 ρ, 张力为 T。微段两端张力在垂直方向的分量差提供恢复力。 |
_x,右端垂直力T \frac{\partial u}{\partial x} |
_{x+dx}。合力为T \left( \frac{\partial u}{\partial x} |
_{x+dx} - \frac{\partial u}{\partial x} |
_x \right) \approx T \frac{\partial^2 u}{\partial x^2} dx。<br>∗∗3.运动方程∗∗:根据牛顿第二定律,(\rho dx) \frac{\partial^2 u}{\partial t^2} = T \frac{\partial^2 u}{\partial x^2} dx,两边除以\rho dx,得\frac{\partial^2 u}{\partial t^2} = \frac{T}{\rho} \frac{\partial^2 u}{\partial x^2}。令c^2 = T/\rho,即得一维波动方程。<br>∗∗4.杆的纵向振动∗∗:类似推导,用应力应变关系\sigma = E\epsilon和牛顿定律,得到相同形式的方程,其中c^2 = E/\rho$。 |
条件:均匀介质;小变形;忽略阻尼和外部激励;张力 T恒定(对于弦)。 |
牛顿第二定律、弹性力学 |
场景:弦乐器(吉他、钢琴)的振动分析;地震波在地层中的一维传播模型;电缆或绳索的波动;杆状结构的纵向冲击响应。 |
|
TH-D1-0071 |
系统动力学 |
连续系统振动 |
梁的横向自由振动微分方程(欧拉-伯努利梁) |
对于细长、等截面梁的横向小振幅自由振动,忽略剪切变形和转动惯量,其运动方程为: |
1. 受力分析:从梁中取微段 dx。截面上的内力有剪力 Fs和弯矩 M。根据材料力学,弯矩与曲率关系:M=−EI∂x2∂2w(小变形下,曲率近似为 w′′)。 |
(定理编号:TH-D1-0071 至 TH-D1-0091)
|
编号 |
类别 |
领域 |
定理/定律/原理名称 |
核心表述(公式/叙述) |
定理的逐步思考推理过程、证明思路及关键数学步骤 |
成立条件与适用范围 |
底层规律/理论归属 |
典型应用场景和工程意义 |
变量/参数列表及说明 |
数学与逻辑特征 |
定理的演绎/应用流程步骤 |
物理本质与“流”描述 |
|---|---|---|---|---|---|---|---|---|---|---|---|---|
|
TH-D1-0071 |
系统动力学 |
连续系统振动 |
梁的横向自由振动微分方程(欧拉-伯努利梁) |
对于细长、等截面梁的横向小振幅自由振动,忽略剪切变形和转动惯量,其运动方程为: |
1. 受力分析:从梁中取微段 dx。截面上的内力有剪力 Fs和弯矩 M。根据材料力学,弯矩与曲率关系:M=−EI∂x2∂2w(小变形下,曲率近似为 w′′)。 |
条件:细长梁(长度远大于截面尺寸);材料线弹性;小变形(曲率近似为 w′′);平截面假定;忽略剪切变形和截面转动惯量。 |
弹性力学、达朗贝尔原理 |
场景:桥梁、楼板、飞机机翼、转子叶片、精密仪器悬臂等结构的自由振动分析;确定结构的固有频率和振型。 |
w(x,t):梁在位置 x处、时间 t时的横向位移。 |
四阶空间导数:反映了弯曲变形与内力(弯矩)的关系。 |
1. 建立方程:根据模型写出 EIw′′′′+ρAw¨=0。 |
描述了梁在弯曲振动中“惯性流”与“弹性恢复力流”的动态平衡。方程左边第一项 EIw′′′′代表由于梁的弯曲变形而产生的弹性恢复力流密度。它是弯矩的二阶导数(即剪力梯度),最终表现为分布力。第二项 ρAw¨代表单位长度梁的惯性力流密度(质量乘以加速度)。方程表明,在自由振动中,每一点上因弯曲变形产生的净弹性力,完全用于平衡该点质量的惯性力。这是一个波动方程,其解对应着特定频率 ωn和空间形态 Wn(x)的驻波,即模态。不同模态代表了不同频率下,“惯性流”和“弹性恢复力流”在梁空间分布上达到共振平衡的状态。 |
|
TH-D1-0072 |
系统动力学 |
连续系统振动 |
模态叠加法 |
对于线性连续系统(如梁、板、壳),其受迫振动的响应可以表示为各阶固有模态(振型)的线性叠加: |
1. 模态展开假设:基于线性系统振动的完备性,任何满足边界条件的位移场 w(x,t)可展开为固有振型 Wn(x)的级数。 |
条件:系统是线性的;阻尼是比例阻尼或经典阻尼,以保证解耦;固有振型 Wn(x)已知且构成完备集。 |
线性代数、特征值问题、振动理论 |
场景:建筑结构在地震载荷下的响应分析;飞机机翼在湍流中的振动;汽车车身在路面激励下的动态响应;任何复杂线性结构的瞬态或稳态动力分析。 |
Wn(x):第 n阶模态振型函数(特征函数)。 |
坐标变换:从物理坐标 w(x,t)变换到模态坐标 qn(t)。 |
1. 进行模态分析:求解系统的无阻尼自由振动方程,得到固有频率 ωn和振型 Wn(x)。 |
体现了线性系统动力响应可以分解为一系列独立“模态通道”的贡献。每个模态 Wn(x)定义了系统的一种特定的、空间固定的振动形态。广义坐标 qn(t)则描述了该形态随时间参与的“强度”或“振幅”。模态叠加法的核心在于正交性,它保证了不同模态之间的能量不会相互交换,就像一个多维空间中的正交基。外部激励 f(x,t)通过投影到各个模态上(计算广义力 Fn(t))来“激发”这些模态。每个被激发的模态就像一个独立的单自由度振子,按照自己的固有频率 ωn和阻尼比 ζn演化。最终观测到的物理振动,是所有被激发模态的振动“流”在空间上的线性叠加。这种方法将复杂的连续系统振动,转化为多个简单振动的合成。 |
|
TH-D1-0073 |
系统动力学 |
多体系统动力学 |
凯恩方程 |
对于具有 n个自由度的多体系统,凯恩方程表述为: |
1. 引入广义速度:选择广义坐标 qk和广义速度 q˙k。偏速度定义为质心速度对广义速度的偏导数:vi(k)=∂q˙k∂vi。 偏角速度类似定义。 |
条件:系统约束是理想的;需要正确计算偏速度和偏角速度。 |
分析力学、达朗贝尔原理 |
场景:机器人动力学建模与控制;航天器姿态动力学;车辆多体动力学仿真;复杂机构(如并联机床)的运动方程自动生成。 |
qk:第 k个广义坐标。 |
功率形式:方程本质上是功率的平衡(主动力功率与惯性力功率之和为零)。 |
1. 选择广义速度 q˙k。 |
凯恩方程是达朗贝尔原理在广义速度空间(切空间)上的投影。它将系统的动力学平衡条件,表达为在每一个独立的“运动方向”(由偏速度 vi(k)和 ωi(k)张成)上,主动力的“功率流”与惯性力的“功率流”之和必须为零。偏速度可以理解为:当仅让第 k个广义速度发生单位变化,而其他广义速度保持不变时,各个刚体所产生的虚拟运动模式。凯恩方程巧妙地避开了约束力的直接计算,因为对于理想约束,约束反力在这些“虚拟运动模式”上不做功(功率为零)。因此,方程直接建立了驱动系统的“主动功率流”与抵抗运动的“惯性功率流”之间的平衡关系,高效地描述了系统广义动量的变化率。 |
|
TH-D1-0074 |
系统动力学 |
多体系统动力学 |
罗伯逊-维滕堡方法(Roberson-Wittenburg Method) |
一种基于图论和铰接体理论的系统化方法,用于推导复杂多体系统的动力学方程。其核心是利用通路矩阵和变换矩阵来描述系统中刚体之间的相对运动,并最终组装成系统的整体质量矩阵和力向量,方程形式为: |
1. 系统拓扑描述:将多体系统抽象为一个有根树(或带环图)。每个刚体是一个顶点,铰(关节)是连接顶点的边。定义每个体相对于其内接体(向根方向)的铰坐标 qi。 |
条件:系统由刚体通过理想铰连接而成;铰的特性(自由度、运动学关系)已知。 |
图论、旋量理论、分析力学 |
场景:工业机器人逆动力学计算;车辆平顺性与操纵稳定性仿真;卫星编队或空间机械臂的动力学建模;生物力学中人体运动的模拟。 |
q:系统的广义坐标向量(通常由各铰的相对坐标组成)。 |
系统化/程式化:有固定的步骤(运动学递推、动力学组装)。 |
1. 系统拓扑描述:定义刚体、铰、建立有根树结构,编号。 |
该方法将复杂多体系统视为一个由“运动流”和“力流”网络构成的拓扑结构。前向运动学递推沿着树的枝干从根到叶传递“运动流”:每个铰的相对运动(广义速度 q˙i)像水流一样,从父物体流向子物体,叠加并变换,最终决定每个末端执行器的运动。后向动力学递推则相反,它从叶到根传递“力流”:末端物体的惯性效应和外力,通过各个铰关节,反向传递并叠加,最终在基座或广义坐标空间形成总的广义惯性力和广义主动力。通路矩阵和变换矩阵精确地描述了这种“运动流”和“力流”的传递路径与变换关系。最终方程 Mq¨+C=Q是整个系统在广义坐标空间中的“动量平衡”表述,其中 M描述了系统的惯性分布如何影响广义加速度,C反映了由于系统内部运动(科里奥利效应、离心效应)产生的“惯性耦合流”。 |
|
TH-D1-0075 |
系统动力学 |
分析力学 |
哈密顿原理 |
在相同的时间端点 t1和 t2,以及相同的边界条件下,保守系统(或有势力系统)的真实运动路径,使作用量 S取驻值(通常为极小值): |
1. 作用量定义:定义泛函 S[q(t)]=∫t1t2L(q(t),q˙(t),t)dt, 它依赖于系统的路径 q(t)。 |
{t_1}^{t_2} - \int{t_1}^{t_2} \frac{d}{dt}\left( \frac{\partial L}{\partial \dot{q}} \right) \delta q dt。<br>由于端点变分为零,边界项消失。<br>∗∗4.驻值条件∗∗:代入得\delta S = \int_{t_1}^{t_2} \left[ \frac{\partial L}{\partial q} - \frac{d}{dt}\left( \frac{\partial L}{\partial \dot{q}} \right) \right] \delta q dt = 0。<br>由于\delta q是任意的,根据变分法基本引理,被积函数必须为零:<br>\frac{d}{dt}\left( \frac{\partial L}{\partial \dot{q}} \right) - \frac{\partial L}{\partial q} = 0$。 |
条件:系统是完整的(约束可积);主动力有势(或部分有势,非有势部分可单独处理);比较路径具有相同的端点时间和状态。 |
变分法、最小作用量原理 |
场景:从第一性原理推导复杂机械系统、电路系统、电磁场系统的运动方程或场方程;有限元法中的弱形式推导;最优控制理论(如庞特里亚金最小值原理与之相关)。 |
S:作用量,是一个泛函。 |
变分原理:真实路径使某个泛函(作用量)取驻值。 |
1. 确定系统:明确系统的广义坐标 qi。 |
|
TH-D1-0076 |
系统动力学 |
分析力学 |
诺特定理 |
对于一个力学系统,如果其作用量在某种连续变换下保持不变(即具有某种对称性),那么必然存在一个相应的守恒量。 |
1. 对称变换:考虑依赖于参数 ϵ的连续变换 qi(t)→Qi(t,ϵ), 且当 ϵ=0时,Qi(t,0)=qi(t)。其生成元为 $\Phi_i = \frac{\partial Q_i}{\partial \epsilon} |
{\epsilon=0}。<br>∗∗2.作用量不变性∗∗:假设在该变换下,拉格朗日量可能变化一个全导数项,即L(Q, \dot{Q}, t) = L(q, \dot{q}, t) + \epsilon \frac{dG}{dt} + O(\epsilon^2)。则作用量S的变分为\delta S = \epsilon (G(t_2)-G(t_1)),仅依赖于端点。对于端点固定的变分,\delta S=0要求G(t_2)=G(t_1),但更一般地,作用量的不变性意味着运动方程的不变性。<br>∗∗3.推导守恒量∗∗:计算拉格朗日量在变换下的变分\delta L,并利用运动方程(欧拉−拉格朗日方程)进行化简。经过推导可得:<br>\frac{d}{dt} \left( \sum_i \frac{\partial L}{\partial \dot{q}i} \Phi_i - G \right) = 0。<br>因此,量I = \sum_i p_i \Phi_i - G是守恒的,其中p_i = \frac{\partial L}{\partial \dot{q}i}为广义动量。<br>∗∗4.经典对应∗∗:<br>−∗∗时间平移不变性∗∗(t \to t+\epsilon):若L不显含时间,则能量H = \sum_i p_i \dot{q}i - L守恒。<br>−∗∗空间平移不变性∗∗(x_i \to x_i + \epsilon):若L在某个方向平移下不变,则该方向的线动量守恒。<br>−∗∗空间旋转不变性∗∗:若L$ 在旋转下不变,则角动量守恒。 |
条件:系统可由作用量原理描述;存在连续对称变换使作用量保持不变(至多相差一个边界项)。 |
群论、变分法、理论物理 |
场景:在粒子物理中,由规范对称性导出电荷守恒;在广义相对论中,由时空对称性导出能量-动量守恒;在经典力学中,快速判断系统的守恒量(能量、动量、角动量)。 |
ϵ:无穷小变换参数。 |
深刻对应:对称性 ↔守恒律。 |
1. 识别对称性:分析系统的拉格朗日量 L, 找出其在何种连续变换下保持不变(或仅改变一个全导数)。 |
|
TH-D1-0077 |
系统动力学 |
稳定性理论 |
李雅普诺夫直接法(第二方法) |
对于自治系统 x˙=f(x), f(0)=0(原点为平衡点)。如果存在一个标量函数 V(x), 满足: |
1. 基本思想:不直接求解微分方程,而是构造一个类似于“能量”的标量函数 V(x)(李雅普诺夫函数),通过分析 V及其沿系统轨迹的导数 V˙的符号来判断稳定性。 |
好的,为您补充系统动力学部分的核心定理、定律和原理,编号从 TH-D1-0078 至 TH-D1-0091。这些内容构成了从经典力学基础到现代动力系统理论的桥梁。
以下是补充内容:
|
编号 |
类别 |
领域 |
定理/定律/原理名称 |
核心表述(公式/叙述) |
定理的逐步思考推理过程、证明思路及关键数学步骤 |
成立条件与适用范围 |
底层规律/理论归属 |
典型应用场景和工程意义 |
变量/参数列表及说明 |
数学与逻辑特征 |
定理的演绎/应用流程步骤 |
物理本质与“流”描述 |
|---|---|---|---|---|---|---|---|---|---|---|---|---|
|
TH-D1-0078 |
系统动力学 |
经典力学 |
牛顿第二定律 |
物体动量随时间的变化率等于作用在该物体上的合外力:F=dtdp=dtd(mv)。对于质量不变的系统,简化为 F=ma。 |
1. 实验基础:基于伽利略斜面实验和牛顿本人的研究,总结出力与加速度的瞬时正比关系。 |
条件:适用于宏观、低速(远低于光速)的物体;在惯性参考系中成立;物体可视为质点或刚体(平动)。 |
实验归纳、公理化体系 |
场景:所有机械系统的运动分析和受力计算;结构动力学中惯性力的计算(达朗贝尔原理);车辆、飞机、航天器的轨迹规划与控制。 |
F:作用在物体上的合外力(矢量)。 |
因果律:力是原因,加速度是瞬时结果。 |
1. 受力分析:隔离研究对象,画出所有外力。 |
牛顿第二定律描述了“动量流”的守恒关系。外力 F是输入系统的动量流率,即单位时间内注入系统的动量。物体的动量 p是系统存储的动量存量。定律 F=dp/dt表明:外力引起的动量流入率,等于系统内动量存量的变化率。当质量不变时,动量变化体现为速度变化(加速度)。这一定律是“流平衡”思想在力学中最原始的体现:推动一个物体,就是在持续地向它输送动量。 |
|
TH-D1-0079 |
系统动力学 |
分析力学 |
达朗贝尔原理 |
在任意瞬时,作用于一个质点系上的所有主动力 Fi和惯性力 −miai在系统的任何虚位移 δri上所做的虚功之和为零: |
1. 牛顿第二定律:对每个质点有 Fi+Ri=miai, 其中 Ri是约束力。 |
条件:约束是理想的(约束力在虚位移上不做功);虚位移满足瞬时约束条件。 |
牛顿力学、虚功原理 |
场景:复杂机械系统(如机构、机器人)的动力学建模;多体系统动力学方程的推导;结构动力学中动载荷的分析。 |
Fi:作用在质点 i上的主动力。 |
动静法:将动力学问题转化为静力平衡问题。 |
1. 分析系统:确定所有主动力 Fi。 |
达朗贝尔原理是“力流”平衡的瞬时表述。它将真实的主动力 Fi和虚构的惯性力 −miai放在同等地位,认为在每一瞬间,这些“力”在系统可能的虚位移模式上所做的“功率流”之和为零。惯性力 −ma代表了质量对运动状态改变的“惯性抵抗流”。原理的本质是:在由真实力和惯性力共同构成的“广义力场”中,系统在任意虚位移方向上都没有净的“功率流”。这保证了动力学演化路径与约束的兼容性。 |
|
TH-D1-0080 |
系统动力学 |
分析力学 |
虚功原理 |
对于一个处于平衡状态的静力学系统,所有主动力在系统的任何一组虚位移上所做的虚功之和为零: |
1. 平衡条件:系统平衡时,每个质点合力为零:Fi+Ri=0, 其中 Ri为约束力。 |
条件:系统处于静力学平衡;约束是理想的、定常的;虚位移是无穷小、满足约束的任意位移。 |
静力学、变分原理 |
场景:机构静力分析(求平衡位置、约束反力);结构力学中利用虚位移法求内力;柔性机构的平衡构型分析。 |
Fi:作用在质点 i上的主动力。 |
标量化:将矢量平衡问题转化为标量虚功计算。 |
1. 确认平衡:明确系统处于静止状态。 |
虚功原理描述了静力学系统中“能量流”的平衡。虚位移 δri是系统在平衡位置附近所有可能发生的微小“构型流”。主动力 Fi在这些虚位移上做的虚功 δW, 代表了如果系统真的沿该方向发生位移,外力将输入系统的能量流。原理 δW=0表明:在平衡位置,沿任何可能的虚位移方向,净能量流均为零。系统没有自发沿任何方向运动的趋势,所有能量输入通道都被“堵塞”或平衡。这是静力学平衡最本质的“流”表述。 |
|
TH-D1-0081 |
系统动力学 |
分析力学 |
哈密顿原理 |
在所有满足固定端点条件 q(t1)=q1, q(t2)=q2的可能运动路径 q(t)中,真实发生的路径是使作用量 S取驻值(通常为极小值)的那一条: |
1. 作用量定义:定义泛函 S[q(t)]=∫t1t2L(q,q˙,t)dt。 |
{t_1}^{t_2} - \int{t_1}^{t_2} \frac{d}{dt}\left( \frac{\partial L}{\partial \dot{q}} \right) \delta q dt = - \int{t_1}^{t_2} \frac{d}{dt}\left( \frac{\partial L}{\partial \dot{q}} \right) \delta q dt。<br>∗∗4.得到欧拉−拉格朗日方程∗∗:代入得\delta S = \int{t_1}^{t_2} \left[ \frac{\partial L}{\partial q} - \frac{d}{dt}\left( \frac{\partial L}{\partial \dot{q}} \right) \right] \delta q dt = 0。由于\delta q任意,故被积函数必须为零:\frac{\partial L}{\partial q} - \frac{d}{dt}\left( \frac{\partial L}{\partial \dot{q}} \right) = 0$。这就是系统的运动方程。 |
条件:系统是完整、保守的;作用量积分存在;端点固定;拉格朗日量 L足够光滑。 |
变分法、最小作用量原理 |
场景:推导复杂机械系统的运动方程;经典力学到量子力学的过渡(路径积分);最优控制理论(变分法);有限元法中的能量原理。 |
S:作用量,是一个泛函。 |
变分原理:真实路径由泛函驻值条件决定。<br全局性:比较整条路径,而非瞬时状态。 |
1. 建立拉格朗日量:根据系统动能 T和势能 V, 写出 L(q,q˙,t)=T−V。 |
|
TH-D1-0082 |
系统动力学 |
分析力学 |
拉格朗日方程 |
对于具有 n个广义坐标 qi的完整系统,运动方程为: |
1. 从哈密顿原理导出:这是最优雅的推导。由哈密顿原理 δ∫Ldt=0, 经过变分计算(如 TH-D1-0081 所述),直接得到上述方程。 |
条件:系统是完整的(约束方程不含速度,或可积分);理想约束;拉格朗日量 L足够光滑。 |
变分原理、分析力学 |
场景:机器人动力学建模与控制;多体系统仿真;航天器姿态动力学;车辆悬架系统分析;任何复杂机械系统的运动方程推导。 |
qi:第 i个广义坐标。 |
二阶常微分方程组:关于 qi(t)的方程。 |
1. 选择广义坐标:确定能唯一确定系统位形且独立的坐标 q1,...,qn。 |
拉格朗日方程描述了广义坐标空间中“广义动量流”与“广义力势”的平衡。定义广义动量 pi=∂L/∂q˙i, 它是系统在广义方向 qi上的“运动趋势存量”。方程 dtdpi=∂qi∂L+Qi表明:广义动量 pi的变化率(广义动量流),由两部分驱动:一是拉格朗日量对坐标的偏导 ∂L/∂qi(来自保守势的“广义力势流”),二是非保守广义力 Qi(外部输入的“广义力流”)。动能 T贡献了广义动量,势能 V贡献了广义力。这个方程统一了各种力学系统的演化规律,其核心是广义能量(拉格朗日量)在广义坐标和速度空间中的“梯度流”平衡。 |
|
TH-D1-0083 |
系统动力学 |
分析力学 |
哈密顿正则方程 |
通过勒让德变换,将拉格朗日方程化为关于广义坐标 qi和广义动量 pi的一阶对称方程组: |
1. 定义广义动量:pi=∂q˙i∂L。 |
条件:变换 pi=∂L/∂q˙i可逆(即 Hessian 矩阵 ∂2L/∂q˙i∂q˙j非奇异);系统完整。 |
勒让德变换、辛几何 |
场景:天体力学和航天轨道计算;分子动力学模拟;量子力学的经典对应;最优控制(哈密顿体系);动力系统理论的研究。 |
qi:广义坐标。 |
一阶对称方程组:关于 q和 p的方程形式对称优美。 |
1. 建立拉格朗日量:L=T−V。 |
哈密顿正则方程描述了相空间中“状态流”的哈密顿梯度场。相空间点 z=(q,p)代表系统的完整状态。哈密顿量 H是这个空间上的一个标量场(能量场)。方程 q˙=∂H/∂p, p˙=−∂H/∂q表明:状态点 z的运动速度,由哈密顿量 H的梯度场 ∇H经过一个固定的旋转(辛矩阵 J)所决定。这个“流” z˙=J∇H就是哈密顿流。它具有关键性质:1. 沿流 H守恒(能量守恒);2. 流保持相空间的辛结构(面积/体积守恒)。因此,系统的演化可以看作相空间中一个不可压缩的、沿等能面流动的“流体”。这是对动力学最几何化、最深刻的描述之一。 |
|
TH-D1-0084 |
系统动力学 |
稳定性理论 |
李雅普诺夫直接法(第二方法) |
设 x=0是系统 x˙=f(x)的一个平衡点。如果存在一个标量函数 V(x), 满足: |
1. 能量类比:李雅普诺夫函数 V(x)可以看作系统的一种“广义能量”。稳定性要求能量为正,且沿轨迹不断衰减(V˙<0)。 |
条件:系统是自治的(或可化为自治);平衡点位于原点;能构造出满足条件的李雅普诺夫函数 V(x)。 |
稳定性理论、能量法 |
场景:机器人控制器的稳定性设计(如计算力矩控制);电力系统暂态稳定性分析;航空航天器姿态稳定控制;化学反应器的稳定性分析。 |
x:系统状态向量。 |
充分条件:找到 V函数则稳定,找不到不能断定不稳定。 |
1. 确定平衡点:通常平移至原点 x=0。 |
李雅普诺夫直接法描述了系统状态空间中“能量流”的耗散过程。李雅普诺夫函数 V(x)定义了状态空间中的一个“能量盆地”。V正定意味着原点处于盆地底部。V˙(x)=∇V⋅f(x)是 V沿系统轨迹方向的方向导数,代表了能量沿轨迹的变化率,即“能量流”。负定 V˙<0意味着无论系统从盆地中的哪一点出发,其“能量”都严格单调递减,状态点必然沿着能量下降的方向“流”向盆地底部(原点),从而实现渐近稳定。负半定 V˙≤0意味着能量不增加,状态点不会跑出某个等能面,但不保证一定流向原点,可能停留在某个不变集上(稳定但非渐近)。该方法的核心思想是:稳定性取决于能否找到一个“能量函数”,其沿系统动力学的“流”是耗散的。 |
|
TH-D1-0085 |
系统动力学 |
稳定性理论 |
李雅普诺夫间接法(第一方法) |
对于非线性系统 x˙=f(x), 设 x=0是平衡点,f连续可微。计算雅可比矩阵 A=Df(0)。则: |
1. 线性化:在平衡点附近将系统展开:x˙=Ax+g(x), 其中 g(x)包含高阶项,且 ∥g(x)∥/∥x∥→0当 x→0。 |
条件:f在平衡点处连续可微;只给出平衡点局部稳定性的信息。 |
线性系统理论、特征值分析 |
场景:快速评估控制系统在工作点附近的局部稳定性;判断机器人平衡姿态的稳定性;分析电力系统运行点的静态稳定性;化学反应平衡点的稳定性初步判断。 |
x:系统状态向量。 |
局部性:结论仅适用于平衡点足够小的邻域内。 |
1. 计算雅可比矩阵:$A = \left. \frac{\partial \mathbf{f}}{\partial \mathbf{x}} \right |
_{\mathbf{x}=0}。<br>∗∗2.求特征值∗∗:计算A的特征值\lambda_1, ..., \lambda_n。<br>∗∗3.检查实部符号∗∗:<br>−若所有\text{Re}(\lambda_i) < 0,则原点∗∗局部渐近稳定∗∗。<br>−若存在某个\text{Re}(\lambda_i) > 0,则原点∗∗不稳定∗∗。<br>−若所有\text{Re}(\lambda_i) \le 0,且存在\text{Re}(\lambda_i)=0$, 则无法判断,需用中心流形定理或直接法。 |
|
TH-D1-0086 |
系统动力学 |
动力系统 |
庞加莱-本迪克松定理 |
在平面上,如果一个动力系统的轨迹被限制在一个有界闭区域 D内,且该区域不包含任何平衡点,那么这条轨迹要么本身是一个周期轨道(极限环),要么随时间趋近于一个周期轨道(即 ω-极限集是一个周期轨道)。 |
1. 平面限制:定理的关键前提是系统状态空间是二维的(平面)。在二维中,轨迹的拓扑可能性有限。 |
条件:系统是连续时间自治系统;状态空间是二维平面(或二维流形);轨迹进入并停留在某个有界闭区域 D内;D内不含系统的平衡点。 |
平面拓扑、微分方程 |
场景:分析非线性电路(如 Van der Pol 振子)中极限环的存在性;研究生态学中捕食者-猎物模型的周期振荡;证明某些机械振动系统(如自激振动)必然产生周期运动。 |
D:平面上的一个有界闭区域(紧致集)。 |
二维特性:强烈依赖于平面拓扑,在高维不成立。 |
1. 建立模型:获得系统的微分方程,确认是二维自治系统。 |
好的,继续为您系统性地梳理和呈现系统动力学部分的核心定理、定律和原理,编号从 TH-D1-0087 至 TH-D1-0091。这些内容涵盖了非线性系统稳定性分析、几何理论及互联系统稳定性的关键原理。
以下是补充内容:
|
编号 |
类别 |
领域 |
定理/定律/原理名称 |
核心表述(公式/叙述) |
定理的逐步思考推理过程、证明思路及关键数学步骤 |
成立条件与适用范围 |
底层规律/理论归属 |
典型应用场景和工程意义 |
变量/参数列表及说明 |
数学与逻辑特征 |
定理的演绎/应用流程步骤 |
物理本质与“流”描述 |
|---|---|---|---|---|---|---|---|---|---|---|---|---|
|
TH-D1-0087 |
系统动力学 |
非线性动力学 |
庞加莱-本迪克森定理 |
考虑一个在平面(二维)上的连续动力系统 x˙=f(x), x∈R2, 其中 f连续可微。如果系统的正半轨 γ+(从某点出发的未来轨迹)有界且不趋于任何平衡点,那么该正半轨的极限集 ω(γ+)要么是一个周期轨道(极限环),要么包含一个同宿轨道或异宿环。特别地,如果一个有界闭区域 D是正向不变的(从 D内出发的轨迹始终留在 D内),并且 D内不含平衡点,则 D内至少存在一个周期轨道(极限环)。 |
1. 平面拓扑限制:在二维平面上,若一条有界轨迹不自交且不趋于平衡点,由于平面 Jordan 曲线定理的限制,它只能以越来越接近某个闭曲线的方式运动。 |
条件:系统是二维的、自治的、连续可微的;轨迹有界且不趋于平衡点;或者存在一个不含平衡点的有界正向不变集。 |
微分方程定性理论、平面拓扑学 |
场景:证明生物种群竞争模型(如 Lotka-Volterra 模型)的修正版本中存在周期振荡;分析化学反应器中是否存在稳定的周期性状态;研究生态系统中捕食者-猎物数量的周期性波动;证明某些机械或电子振荡器(如 Van der Pol 振子)极限环的存在性。 |
x:二维状态向量 (x,y)。 |
二维性:结论强烈依赖于二维平面的拓扑性质。 |
1. 构造正向不变集:找到一个有界闭区域 D, 并证明从 D内出发的轨迹永远不会离开 D(常用方法是证明在 D的边界上,向量场都指向内部)。 |
庞加莱-本迪克森定理描述了二维有界“状态流”的必然归宿。在二维平面上,状态空间像一张无限延伸的橡胶膜。有界正向不变区域 D就像是膜上一个封闭的“池塘”,水流(向量场)在池塘边界处都指向内部,因此一旦进入池塘的水滴(状态点)就再也出不去。由于池塘内没有静止的“漩涡中心”(平衡点),水滴不能永远停在一个点。在二维的约束下,水滴的运动选项非常有限:它不能像在三维中那样无限复杂地缠绕(混沌)。最终,水滴的轨迹要么自己首尾相接形成一个闭合的环形水道(极限环),要么无限逼近于这样一个环形水道。定理的本质是:在二维有界且无平衡点的流场中,任何轨迹必然渐近地趋于一个周期性的“循环流”。这个循环流就是极限环,它代表了系统的一种稳定的周期运动模式。 |
|
TH-D1-0088 |
系统动力学 |
非线性动力学 |
哈特曼-格罗布曼定理 |
考虑非线性自治系统 x˙=f(x), 设 xe是一个双曲平衡点,即其雅可比矩阵 A=Df(xe)的所有特征值的实部均不为零。则在平衡点 xe的某个邻域 U内,存在一个同胚 h:U→V(V是 Rn中包含原点的开集),使得非线性流 ϕt与线性化系统 y˙=Ay的流 ψt在此邻域内拓扑共轭。即,对于所有 x∈U和足够小的 t, 有 h(ϕt(x))=ψt(h(x))。换言之,非线性系统在双曲平衡点附近的相图与其线性化系统的相图拓扑等价。 |
1. 线性化与误差:在平衡点附近,系统可写为 x˙=Ax+g(x), 其中 g(x)=o(∥x∥)为高阶非线性项。 |
条件:平衡点是双曲的(线性化矩阵无非零实部的特征值);向量场 f至少是 C1的(连续可微)。 |
微分方程几何理论、动力系统 |
场景:快速判断非线性系统在平衡点附近的局部稳定性。如果线性化系统是渐近稳定的(所有特征值实部为负),则原非线性系统在该点也是局部渐近稳定的;如果线性化系统有正实部特征值,则原系统不稳定。在机器人平衡点稳定性分析、电力系统暂态稳定性的初步判断、化学反应平衡点分析中广泛应用。 |
xe:非线性系统的平衡点。 |
局部性:结论仅在平衡点的小邻域内成立。 |
1. 求平衡点:解 f(xe)=0。 |
哈特曼-格罗布曼定理揭示了在双曲平衡点附近,非线性“流”的局部拓扑结构完全由其线性近似主导。将非线性系统在平衡点附近展开,x˙=Ax+g(x), 其中 Ax是线性主导流,g(x)是高阶扰动流。双曲性条件(特征值实部非零)意味着线性主导流在稳定和不稳定方向上都有指数级的收缩或扩张速率,这个速率远大于高阶扰动流的影响(因为 ∥g(x)∥/∥x∥→0当 x→0)。因此,在足够小的邻域内,线性流的指数行为是“强硬”的,它强行规定了轨迹进出平衡点的大致方向(稳定/不稳定流形)。高阶扰动 g(x)只能对这些轨迹进行“弯曲”或“扭曲”,但无法改变它们连接和分离的整体拓扑方式。定理中的同胚 h正是这样一个“橡皮泥映射”,它可以将被非线性扰动弯曲了的轨迹,连续地拉伸和压缩回线性流的直线轨迹,从而在拓扑意义上将两者等同起来。简言之,在双曲点附近,非线性流的“骨架”是由线性项铸就的,非线性项只是附着在上面的“肌肉”。 |
|
TH-D1-0089 |
系统动力学 |
稳定性理论 |
拉萨尔不变性原理 |
考虑自治系统 x˙=f(x), 设 Ω是一个紧致(有界闭)正不变集。设 V:Ω→R是一个连续可微的函数,且在 Ω上满足 V˙(x)=∇V(x)⋅f(x)≤0。令 E={x∈Ω∣V˙(x)=0}, 并令 M为 E中的最大不变集。则从 Ω内任一点出发的解,当 t→∞时,都趋近于 M。即,ω-极限集包含于 M。 |
1. 李雅普诺夫函数性质:由于 V˙≤0, V沿轨迹非增。又因为 V在紧集 Ω上连续,故有下界。因此,V(x(t))单调递减且有下界,必收敛到某个常数 c。 |
条件:系统是自治的;存在一个紧致正不变集 Ω;存在一个 C1函数 V, 在 Ω上 V˙≤0。 |
李雅普诺夫稳定性理论、动力系统 |
场景:分析具有耗散但 V˙仅负半定的物理系统,如带摩擦的机械系统(V˙仅在速度为零时为零);自适应控制系统;梯度下降算法;化学反应网络的质量守恒系统。 |
Ω:状态空间中的一个紧致(有界闭)集,且是正不变的(从 Ω内出发的轨迹始终留在 Ω内)。 |
弱条件:只要求 V˙≤0, 而非负定。 |
1. 构造正不变集:通常利用 V的等值线或水平集来构造一个紧致正不变集 Ω(例如,Ω={x∣V(x)≤c})。 |
拉萨尔不变性原理描述了耗散系统中“能量流”最终停滞的区域。V函数可以看作系统的“能量”或“耗散函数”。条件 V˙≤0意味着系统是耗散的,能量不会增加。集合 E={x∣V˙=0}是能量流为零的“停滞带”,在这里系统没有耗散。然而,轨迹可能只是瞬时穿过 E。拉萨尔原理指出,由于能量 V不断下降(或不变)且有下界,系统最终必须进入一个状态,使得 V不再变化,即轨迹必须永远停留在 V˙=0的区域内。但轨迹本身是动态的,所以它最终必须进入 E中一个动态不变的子集 M——即一个一旦进入就永远不会离开 E的集合。这个集合 M就是系统所有可能的终态的集合。例如,在带摩擦的摆中,V是机械能,V˙≤0(摩擦耗散),V˙=0当且仅当角速度为零(摆瞬时静止)。但摆不会永远停在任意角度,因为重力会使其运动。最终,摆会停在最低点(平衡点)。这里,E是所有速度为零的状态(一条线),而最大不变集 M只是这条线上的一个点(最低点)。原理的本质是:在耗散系统(V˙≤0)中,状态流最终会被吸引到能量不再耗散(V˙=0)且能自我维持(不变)的“动态平衡带” M上。 |
|
TH-D1-0090 |
系统动力学 |
控制理论 |
无源性定理 |
考虑一个系统 H, 其输入为 u(t), 输出为 y(t)。设系统初始状态存储的能量为 V(x(0))≥0。如果存在一个连续可微的半正定函数 V(x)(存储函数),使得对于所有 t≥0和所有输入 u, 满足耗散不等式: |
1. 无源性定义:耗散不等式 V(t)−V(0)≤∫uTydτ物理上表示:系统内部存储能量的增加不超过外部输入供给的能量(uTy是瞬时输入功率)。 |
条件:系统是因果的;存储函数 V(x)连续可微、半正定;耗散不等式成立。对于稳定性结论,通常还需要系统是零状态可观测等条件。 |
耗散系统理论、非线性控制 |
场景:机器人力/位置控制中的阻抗控制设计;电力电子变换器的稳定性分析;网络化控制系统的协议设计;自适应控制中保证全局稳定性的工具;机械系统的能量整形控制(如基于无源性的控制)。 |
u(t):系统的输入向量。 |
能量不等式:核心是一个关于能量的积分不等式。 |
1. 验证子系统无源性:为待互联的子系统 H1和 H2分别找到存储函数 V1, V2, 并验证其满足耗散不等式。 |
无源性定理描述了基于“能量流”互联的系统的稳定性本质。每个无源子系统都可以看作一个“能量处理器”:它从外部输入 u接收功率流 (uTy),一部分能量被存储起来(V增加),一部分可能被耗散。耗散不等式 V(t)−V(0)≤∫uTydτ意味着系统是能量消耗型的,它产生的能量不会多于输入的能量。当两个这样的系统通过负反馈互联时,一个系统的输出作为另一个系统的输入。无源性定理表明,这种互联不会产生“能量创造”的恶性循环。从能量角度看,反馈互联相当于将两个系统的能量端口连接起来。总存储能量 V=V1+V2的变化率满足 V˙≤u1Ty1+u2Ty2。代入反馈连接 u1=r−y2, u2=y1, 得到 V˙≤rTy1。外部输入 r提供的功率流 rTy1是唯一能增加系统总能量的来源。系统内部两个子系统之间的能量交换 (−y2Ty1+y1Ty2=0) 是守恒的,互相抵消了。因此,整个闭环系统对外部而言,仍然表现为一个无源的、能量有界的系统。这就像两个遵守能量守恒的物体耦合在一起,整体也必然遵守能量守恒。如果其中一个子系统是严格无源的(内部有耗散),那么互联系统整体就具有能量耗散机制,从而能镇定外部输入的影响,实现稳定性。 |
|
TH-D1-0091 |
系统动力学 |
控制理论 |
小增益定理 |
考虑两个系统 H1和 H2的反馈互联。设 H1:L2e→L2e具有 L2增益 γ1≥0, H2:L2e→L2e具有 L2增益 γ2≥0。即,存在常数 β1,β2使得对于所有 T≥0, |
1. 反馈互联方程:互联系统方程为 y1=H1(u1), y2=H2(u2), u1=r1+y2, u2=r2−y1(或其他符号约定)。 |
条件:系统 H1, H2是因果的、有限增益 L2稳定的;小增益条件 γ1γ2<1成立。对于非线性时变系统,定理也成立,是输入输出稳定性的核心结果之一。 |
输入输出稳定性理论、算子理论 |
场景:鲁棒控制中基于 H∞范数的控制器设计(小增益定理是 H∞控制的理论基础);分析具有不确定性的反馈系统的鲁棒稳定性;网络控制系统中考虑时延和数据包丢失的稳定性分析;多机器人系统协同控制的稳定性证明。 |
H1,H2:待互联的系统,视为输入输出算子。 |
x(t) |
^2 dt}。<br>|\cdot |T:截断信号的\mathcal{L}2范数,|x_T|_2 = \sqrt{\int_0^T |
x(t) |
(定理编号:TH-D1-0092 至 TH-D1-0102)
|
编号 |
类别 |
领域 |
定理/定律/原理名称 |
核心表述(公式/叙述) |
定理的逐步思考推理过程、证明思路及关键数学步骤 |
成立条件与适用范围 |
底层规律/理论归属 |
典型应用场景和工程意义 |
变量/参数列表及说明 |
数学与逻辑特征 |
定理的演绎/应用流程步骤 |
物理本质与“流”描述 |
|---|---|---|---|---|---|---|---|---|---|---|---|---|
|
TH-D1-0092 |
系统动力学 |
分析力学 |
哈密顿-雅可比方程 |
对于哈密顿系统,存在一个母函数 S(q,α,t),称为哈密顿主函数,它是如下一阶非线性偏微分方程的解: |
1. 动机:寻找一个正则变换,使得新哈密顿量 K恒为零,则新坐标 α和新动量 β均为常数,运动方程平凡可解。 |
条件:系统可由哈密顿量描述;哈密顿-雅可比方程可解(通常通过分离变量法)。 |
偏微分方程理论、正则变换 |
场景:天体力学中二体问题、开普勒轨道的精确求解;量子力学中 WKB 近似(经典作用量);最优控制理论(Hamilton-Jacobi-Bellman 方程)。 |
S(q,α,t):哈密顿主函数,也是作用量。 |
一阶非线性 PDE:方程关于 S及其偏导是非线性的。<br特征线法:该 PDE 的特征线方程就是原哈密顿正则方程。 |
1. 写出哈密顿量 H(q,p,t)。 |
哈密顿-雅可比方程描述了作用量 S在状态空间和时间中的“传播”规律。主函数 S可以看作是系统从某个初始状态演化到当前状态 (q,t)所累积的作用量。H-J 方程 ∂S/∂t+H(q,∇qS,t)=0类似于一个动力学波动方程,决定了作用量“波前”的传播。其中,哈密顿量 H扮演了“色散关系”的角色。方程的解 S(q,α,t)定义了一个在状态空间 (q,t)中移动的“作用量场”。系统的真实运动轨迹,就是这个场的“梯度线”(由 p=∇qS定义),垂直于等 S面。常数 α,β则标记了不同的“作用量射线”。这深刻揭示了力学与几何光学的类比:粒子运动就像作用量波的射线。 |
|
TH-D1-0093 |
系统动力学 |
非线性动力学 |
庞加莱映射 |
对于一个连续时间动力系统,在其周期轨道附近选取一个超曲面 Σ(庞加莱截面),该截面横截于(即不相切于)流。系统轨迹每次从同一方向穿过 Σ的点构成一个序列 {x0,x1,x2,...}。定义映射 P:Σ→Σ, 使得 xk+1=P(xk)。这个离散映射 P称为庞加莱映射。 |
1. 截面选取:对于 n维连续系统,选取一个 (n−1)维的超曲面 Σ, 通常要求它横截于流,即轨迹在交点处不与之相切。对于周期轨道,常选取通过轨道上某点的截面。 |
条件:系统是连续时间自治或周期非自治系统;能定义合适的横截截面。 |
动力系统理论、微分方程 |
场景:分析非线性振子(如 Duffing, Van der Pol)的周期解稳定性;研究转子动力学中的油膜涡动、喘振;分析电力系统的极限环振荡;混沌系统的数值识别(通过观察庞加莱截面的点集结构)。 |
Σ:庞加莱截面,一个 (n−1)维超曲面。 |
降维:将 n维连续流转化为 (n−1)维离散映射。 |
1. 选择截面:根据系统特性,选择一个合适的庞加莱截面 Σ(如 θ=0平面对于周期激励系统)。 |
庞加莱映射是在状态空间的某个“观察窗”Σ上,对连续动力系统“流”进行“频闪采样”的结果。它抛弃了轨迹在截面之间流动的细节,只记录每次穿过这个“窗口”的瞬间状态。这样,连续的、复杂的轨迹演化,被简化为一个离散的、低维的点列演化 xk+1=P(xk)。不动点 x∗意味着轨迹每次都以相同的状态穿过截面,这对应着一个周期运动。映射 P的雅可比矩阵 DP描述了截面附近状态点的线性化演化,其特征值决定了这个不动点的稳定性:如果所有特征值都“收缩”(模<1),意味着邻近的轨迹在每次穿越时都更靠近不动点,周期轨道稳定;反之则不稳定。庞加莱映射将“流”的稳定性问题,转化为“映射”的稳定性问题,是动力系统理论的核心简化工具。 |
|
TH-D1-0094 |
系统动力学 |
非线性动力学 |
梅尔尼科夫方法 |
用于判断平面哈密顿系统在加入小周期阻尼和激励后,其同宿轨(或异宿轨)是否会发生横截相交,从而预测 Smale 马蹄意义下的混沌出现。 |
1. 未扰系统:考虑一个平面哈密顿系统 q˙=f(q), 它有一个双曲鞍点 S, 以及从 S出发又回到 S的同宿轨 qh(t)。 |
\mathbf{f}(\mathbf{q}^h(0)) |
}。因此,M(t_0)的零点对应W^s和W^u的交点。<br>∗∗5.横截相交∗∗:若M(t_0)$ 有简单零点,则交点是非退化的、横截的。根据 Smale-Birkhoff 同宿定理,这意味着映射存在马蹄,从而系统具有复杂的动力学,包括无穷多个周期轨道和混沌不变集。 |
条件:未扰系统是平面哈密顿系统,存在双曲鞍点和同宿(或异宿)轨;扰动是小量的、光滑的、周期的;方法给出的是必要条件(横截相交导致混沌,但混沌不一定需要横截相交)。 |
微分方程摄动理论、动力系统 |
场景:分析阻尼和周期外力作用下的 Duffing 振子(双势阱)何时会出现混沌;研究参数激励摆的混沌阈值;分析具有非线性刚度的机械系统在周期激励下的复杂动力学。 |
qh(t):未扰系统的同宿轨解。 |
摄动法:基于小参数 ϵ展开。 |
|
TH-D1-0095 |
系统动力学 |
非线性动力学 |
中心流形定理 |
对于一个动力系统 x˙=f(x), 设 x=0是一个平衡点,其雅可比矩阵 Df(0)的特征值实部可分为三类:负实部(稳定)、零实部(中心)、正实部(不稳定)。则存在局部不变的流形 Wc(中心流形)、Ws(稳定流形)、Wu(不稳定流形),分别与中心特征空间、稳定特征空间、不稳定特征空间相切。中心流形 Wc的维数等于中心特征值的个数(按重数计),且系统在 Wc上的动力学决定了平衡点附近的长期行为。 |
1. 线性化与特征空间:在平衡点线性化,x˙=Ax, 其中 A=Df(0)。将状态空间分解为稳定子空间 Es、中心子空间 Ec、不稳定子空间 Eu, 分别由对应特征值的特征向量张成。 |
条件:f是 Cr光滑的 (r≥1);平衡点是非双曲的(即存在实部为零的特征值)。 |
微分方程几何理论、不变流形 |
场景:分析结构失稳后的后屈曲行为;研究非线性振动系统在临界阻尼参数附近的动力学;电力系统电压稳定性分析;化学反应动力学中临界点的分析。 |
x:系统状态向量。 |
降维:将系统动力学约化到低维流形上。 |
1. 线性化:计算平衡点 x=0处的雅可比矩阵 A, 求出其特征值,识别中心模态。 |
中心流形定理揭示了在非双曲平衡点附近,高维动力系统的“慢流形”支配着长期演化。线性化矩阵 A的特征值实部决定了模态的“时间尺度”:负实部对应快衰减的稳定模态,正实部对应快增长的不稳定模态,而零实部对应中性或慢变模态。中心流形 Wc就是由这些慢变模态张成的局部不变曲面。定理表明,所有快变模态(稳定和不稳定)的演化,在局部上都“锁定”或“从属于”慢变模态,即可以表示为慢变量的函数 v=h(u)。因此,系统的长期行为(如是否稳定、是否分岔出周期解)完全由限制在 Wc上的慢变子系统决定。这就像在湍急的河流(快变模态)中,有一个缓慢移动的漩涡中心(中心流形),最终决定漂浮物长期轨迹的,正是这个漩涡的结构,而不是瞬间的湍流。中心流形是连接平衡点与更复杂动力学(如极限环)的“桥梁”。 |
|
TH-D1-0096 |
系统动力学 |
非线性动力学 |
规范形理论 |
对于一个非线性动力系统 x˙=f(x), 在平衡点 x=0附近,可以通过一系列近恒等非线性坐标变换 x=y+h(y), 将系统化简为尽可能简单的形式 y˙=Λy+g(y), 其中 Λ是雅可比矩阵的 Jordan 形,g(y)只包含共振项(即无法通过非线性变换消除的项)。这种简化后的系统称为规范形。 |
1. 泰勒展开:在平衡点附近将 f(x)展开为泰勒级数:x˙=Ax+f(2)(x)+f(3)(x)+..., 其中 f(k)是 k次齐次多项式向量。 |
条件:f是光滑的;平衡点处的线性化矩阵 Λ的特征值已知;通常用于局部分析。 |
李代数、幂级数、动力系统 |
场景:分析非线性振子在平衡点附近的动力学;计算霍普夫分岔的正则形系数以判断分岔类型(超临界或亚临界);研究具有对称性的系统(如哈密顿系统)的规范形;简化控制器设计模型。 |
x,y:变换前后的状态向量。 |
近恒等变换:变换在原点附近接近恒等映射。 |
1. 线性化与 Jordan 化:计算 A=Df(0), 求其特征值和特征向量,将系统通过线性变换化为 z˙=Λz+...。 |
规范形理论旨在通过“扭曲”状态空间的坐标,来“拉直”非线性系统的流线,从而消除所有不必要的非线性耦合,只留下反映系统内在共振结构的本质非线性项。线性变换已经将对角化(或 Jordan 化),消除了线性耦合。非线性变换 x=y+h(y)则试图进一步“校正”坐标,使得在新坐标 y下,动力学的表达式中,那些可以通过非线性坐标变换补偿掉的项都被消除。共振项是那些无法被消除的,因为它们对应的非线性相互作用与线性部分的特征频率匹配,从而产生了持久的、累积的效应。例如,在霍普夫分岔中,特征值为纯虚数 ±iω, 则 (λ1+λ2)=0导致三次项 y12y2或 y1y22可能是共振的(实际上需要更精确的条件)。规范形就像给非线性系统做了一个“外科手术”,切除了所有“良性”的非线性组织,只留下关键的“病灶”(共振项),从而让我们能清晰地诊断系统的动力学疾病。 |
|
TH-D1-0097 |
系统动力学 |
控制理论 |
绝对稳定性准则(圆判据、波波夫准则) |
针对 Lur'e 型非线性系统:x˙=Ax+Bu, y=Cx, u=−ϕ(y), 其中 ϕ(⋅)是满足扇形条件 [k1,k2]的静态非线性函数(即 k1y2≤yϕ(y)≤k2y2)。 |
1. 系统描述:研究一个线性部分 G(s)与一个满足扇形条件的非线性反馈 ϕ(⋅)构成的闭环系统。绝对稳定性指对扇形内所有可能的 ϕ, 平衡点都是全局渐近稳定的。 |
条件:线性部分 (A,B,C)可控可观或最小实现;非线性函数 ϕ满足扇形条件 [k1,k2], 且 ϕ(0)=0;通常假设 A是 Hurwitz 的(开环稳定),但圆判据也可处理部分不稳定情况。 |
频域稳定性理论、正实系统 |
场景:分析具有饱和、死区、继电器等非线性特性的控制系统稳定性;设计鲁棒控制器,保证在非线性不确定性下的全局稳定;航空航天器中含非线性执行器的控制系统分析。 |
G(s):线性部分的传递函数。 |
频域判据:稳定性条件通过传递函数的频率响应图形表达。 |
1. 确定扇形条件:根据非线性特性 ϕ(y), 确定其所属的扇形 [k1,k2]。 |
绝对稳定性准则建立了一个“安全边界”,只要线性部分的频率响应 G(jω)不穿越这个边界,那么无论非线性反馈 ϕ在扇形内如何“扭曲”信号,系统的总“能量流”总是耗散的,从而保证全局稳定。圆判据中的“禁止圆”可以理解为:非线性反馈 ϕ在扇形 [k1,k2]内变化,相当于在环路中引入了一个时变的、但有界的“增益”。为了保证闭环稳定,线性部分 G(s)的 Nyquist 曲线必须与这个增益变化范围所对应的临界区域(圆)保持安全距离。如果 G(s)的曲线始终位于圆外,意味着环路增益的相位和幅值组合永远不会达到-180°且增益大于1的条件,从而不会发生振荡失稳。波波夫准则则通过引入一个频率加权 1+qs, 相当于在时域的李雅普诺夫函数中增加了一个“交叉项”,从而可能扩大可判定的稳定扇形。这些准则将非线性系统的稳定性,巧妙地转化为对线性部分频率响应的几何约束。 |
|
TH-D1-0098 |
系统动力学 |
控制理论 |
无源性定理 |
一个动态系统如果对于所有 t≥0和所有输入输出对 (u(t),y(t)), 存在一个存储函数 S(x)≥0使得以下耗散不等式成立: |
1. 无源性定义:源于电路网络理论。一个无源系统不能产生能量,只能存储或耗散能量。不等式 S˙≤uTy表明存储函数 S的变化率不超过输入功率。 |
条件:系统是耗散的,满足上述不等式;存储函数 S(x)正定且径向无界(对于全局稳定性);可检测性条件用于保证渐近稳定性而不仅仅是稳定性。 |
耗散系统理论、李雅普诺夫稳定性 |
场景:机器人阻抗控制和力控制(机器人动力学通常是无源的);电力系统稳定性分析(发电机、负荷模型的无源性);网络拥塞控制(将数据流建模为无源系统);设计基于能量的控制器(如无源化控制)。 |
u(t):系统输入向量。 |
能量不等式:核心是一个积分不等式或微分不等式。 |
1. 确定输入输出对:为系统定义有物理意义的输入 u和输出 y, 通常使得 uTy具有功率量纲。 |
好的,继续为您系统性地梳理和呈现装备制造与力学系统中的核心定理、定律和原理。本批次将继续严格遵循您的要求,优先考虑系统动力学内容,并继续使用详细格式。
以下是第九部分(定理编号:TH-D1-0099 至 TH-D1-0102)。
|
编号 |
类别 |
领域 |
定理/定律/原理名称 |
核心表述(公式/叙述) |
定理的逐步思考推理过程、证明思路及关键数学步骤 |
成立条件与适用范围 |
底层规律/理论归属 |
典型应用场景和工程意义 |
变量/参数列表及说明 |
数学与逻辑特征 |
定理的演绎/应用流程步骤 |
物理本质与“流”描述 |
|---|---|---|---|---|---|---|---|---|---|---|---|---|
|
TH-D1-0099 |
系统动力学 |
非线性动力学 |
Hopf 分岔定理 |
考虑一个依赖于参数 μ的动力系统 x˙=f(x,μ), 设 x=0对所有 μ都是平衡点。在 μ=0处,假设线性化矩阵 A(μ)=Dxf(0,μ)有一对共轭复特征值 λ(μ),λˉ(μ), 满足: |
_{\mu=0} = d \ne 0(特征值实部以非零速率穿过虚轴)。<br>3.∗∗非退化条件∗∗:系统在\mu=0处的三阶非线性项满足某个系数l_1(0) \ne 0(第一李雅普诺夫系数)。<br>则系统在\mu=0处经历Hopf分岔。当d > 0且l_1(0) < 0时,为∗∗超临界Hopf分岔∗∗:当\mu从负变为正时,平衡点失稳,同时产生一个稳定的小振幅极限环。当l_1(0) > 0$ 时,为亚临界 Hopf 分岔:产生不稳定的极限环。 |
1. 中心流形约化:由于在 μ=0时有一对纯虚特征值,根据中心流形定理,系统动力学可约化到一个二维中心流形上,其方程可写为复数形式: |
z |
^2 + \text{高阶项}。<br>其中\alpha(\mu)=\text{Re}(\lambda(\mu)),\omega(\mu)=\text{Im}(\lambda(\mu))。<br>∗∗3.极坐标分析∗∗:令z = r e^{i\theta},代入规范形,分离实部虚部,得到关于振幅r的方程:<br>\dot{r} = \alpha(\mu) r + l_1(\mu) r^3 + \cdots,\dot{\theta} = \omega(\mu) + \cdots。<br>∗∗4.分岔分析∗∗:对于小的\mu和r,截断到三阶项。平衡点r=0对应原系统的平衡点。非零平衡点满足\alpha(\mu) + l_1(\mu) r^2 = 0,即r = \sqrt{-\alpha(\mu)/l_1(\mu)}(要求\alpha(\mu)/l_1(\mu) < 0)。结合横截条件\alpha'(0)=d \ne 0$, 即可分析分岔方向与极限环稳定性。 |
条件:系统光滑(至少 C3);平衡点处线性化矩阵有一对单重纯虚特征值;横截条件和非退化条件满足。 |
动力系统分岔理论、中心流形定理 |
场景:机翼颤振(当气流速度超过临界值,平衡姿态失稳,产生周期振动);电子振荡器(如 Van der Pol 振子)起振过程;激光器阈值行为;化学反应中的周期性振荡(如 BZ 反应);生物神经元放电的 Hopf 分岔模型。 |
μ:分岔参数。 |
|
TH-D1-0100 |
系统动力学 |
非线性动力学 |
中心流形定理 |
考虑系统 x˙=Ax+f(x,y), y˙=By+g(x,y), 其中 x∈Rn, y∈Rm, A的特征值实部为零(中心子空间 Ec), B的特征值实部为负(稳定子空间 Es)。假设 f,g光滑且在原点的阶数大于1。则存在一个局部不变流形(中心流形)Wlocc, 可表示为 y=h(x), 其中 h(0)=0, Dh(0)=0。系统在中心流形上的约化动力学为: |
1. 空间分解:将系统的状态空间分解为中心子空间 Ec(对应零实部特征值)和稳定/不稳定子空间 Es/u。 |
条件:系统在平衡点处光滑(至少 Cr, r≥1);线性化矩阵可对角化或分块对角化,特征值实部分为零、负(和正)三部分。 |
微分方程几何理论、不变流形理论 |
场景:分析高维系统的静态分岔(如鞍结分岔、跨临界分岔)和动态分岔(如 Hopf 分岔);简化柔性结构(如板、壳)在临界载荷附近的屈曲后行为模型;降阶处理具有慢变参数的快慢系统;研究化学反应网络中关键物种的动力学。 |
x:中心变量(对应中心特征值)。 |
降维:将系统维数从 n+m降至 n(中心子空间维数)。 |
1. 线性化与分解:在平衡点线性化,计算特征值,将状态空间分解为中心子空间和稳定子空间,并引入相应坐标 (x,y)。 |
中心流形定理描述了高维动力系统中“慢流形”对长期动力学的支配作用。在平衡点附近,状态空间可以分解为“快子空间” Es(特征值实部为负)和“慢子空间” Ec(特征值实部为零)。快子空间上的状态变量 y代表系统的“快模态”,它们像被强阻尼吸引,会指数快速地松弛到一个由慢变量 x决定的“准平衡”状态,即 y=h(x)。这个关系定义了一个中心流形 Wc, 它是整个状态空间中一个低维的、局部不变的“慢流通道”。定理表明,所有轨迹都会迅速被吸引到这个慢流通道上,并在其上缓慢演化。因此,系统的长期行为(稳定性、振荡、分岔)完全由这个通道上的“慢流”动力学决定。这就像河流的主航道决定了整体的流向,而岸边的小漩涡(快模态)虽然存在,但不影响主航道的趋势。中心流形定理通过识别并约化到这条“主航道”,实现了对高维系统本质动力学的降维理解和分析。 |
|
TH-D1-0101 |
系统动力学 |
非线性动力学 |
规范形理论 |
对于一个在平衡点附近的光滑动力系统 x˙=Ax+F(x), 其中 F(x)包含非线性项,规范形理论旨在通过一系列近恒等坐标变换 x=y+P(y)(P为多项式),将系统简化为尽可能简单的形式 y˙=Ay+G(y), 其中 G(y)仅包含与线性部分 A共振的非线性项。这种简化后的系统称为规范形。 |
1. 同调方程:设变换为 x=y+Pk(y), 其中 Pk是 k次齐次多项式。将变换代入原方程,并假设 y˙=Ay+Gk(y)+⋯, 其中 Gk是待定的 k次项。在 k次项上匹配系数,得到所谓的同调方程: |
条件:系统在平衡点处足够光滑(Cr, r足够高);线性化矩阵 A可对角化(或至少可化为 Jordan 型)。 |
微分方程、李代数、共振理论 |
场景:分析非线性振子的内部共振(如两个频率满足 ω1≈2ω2);简化分岔(如 Hopf 分岔、双 Hopf 分岔)附近的系统方程,以便于分析;研究哈密顿系统的非线性振动;作为数值计算规范形的算法基础。 |
x,y:变换前后的状态向量。 |
简化:将系统化为含有最少非线性项的形式。 |
1. 线性化与对角化:在平衡点线性化,求得矩阵 A及其特征值。 |
规范形理论是对非线性系统进行“动力学的代数化简”,目的是找到一组“最自然的”局部坐标,使得系统的方程呈现出最简洁的代数结构。线性部分 A定义了一个“主导流”。非线性项 F(x)可以看作对这个线性流的扰动。通过近恒等变换 x=y+P(y), 我们试图“旋转”或“弯曲”坐标架,以对齐这个主导流,从而尽可能多地抵消掉非线性扰动。同调方程 LA(Pk)=Fk−Gk正是这个对齐过程的数学表述:我们希望通过选择变换 Pk, 将非线性项 Fk的大部分“吸收”掉,只剩下无法被吸收的残余部分 Gk。那些无法被消除的项 Gk, 即共振项,是与线性流的本征模式发生长期相互作用的项。例如,在振荡系统中,如果两个模式的频率满足 mω1+nω2≈0(m,n为整数),那么相应的非线性项会导致两个模式间持续的能量交换,无法通过坐标变换消除。因此,规范形剥离了所有“无关紧要”的非线性细节,只保留了决定系统长期动力学拓扑类型的核心共振相互作用,是洞察非线性系统本质的“数学显微镜”。 |
|
TH-D1-0102 |
系统动力学 |
控制理论 |
内模原理 |
对于一个控制系统,如果要在存在外部信号(参考输入或扰动)的情况下实现渐近跟踪或扰动抑制,那么控制器中必须包含一个能再现该外部信号动态模型的子系统,这个子系统称为内模。更形式化地:考虑被控对象 P(s)和控制器 C(s)组成的单位反馈系统。设外部信号 r(t)(参考)或 d(t)(扰动)由一自治微分方程 w˙=Sw, r=Qw(或 d=Qw)生成。则闭环系统能实现渐近跟踪(e(t)=r(t)−y(t)→0)的充分必要条件是: |
1. 必要性思路(反证):假设控制器不包含内模,但实现了渐近跟踪。考虑误差 e(t)=r(t)−y(t)的 Laplace 变换 E(s)。由终值定理,稳态误差为零要求 E(s)的所有极点都在左半平面。但 R(s)(参考输入的 Laplace 变换)的极点位于虚轴或右半平面(对应持续信号)。通过分析闭环传递函数,可以证明,如果控制器没有与 R(s)极点对消的极点,那么 R(s)的这些极点将成为 E(s)的极点,导致稳态误差不为零或无穷。因此,控制器必须“内嵌”参考信号的模型。 |
条件:被控对象和控制器是线性时不变的;外部信号由自治线性系统生成;闭环系统可镇定且可检测。 |
线性系统理论、伺服控制 |
场景:高精度运动控制(如光盘驱动器、机器人轨迹跟踪),要求对正弦参考信号无静差跟踪;主动噪声与振动控制,抑制已知频率的周期性扰动;电网中逆变器的控制,以产生高质量的正弦波;过程控制中,抑制由往复泵等引起的周期性扰动。 |
P(s):被控对象的传递函数。 |
结构性条件:对控制器的极点配置提出了明确要求。 |
1. 建模外部信号:确定参考或扰动信号的模型,如阶跃 (1/s)、斜坡 (1/s2)、正弦 (ω/(s2+ω2)) 或其组合。得到其最小多项式或状态空间模型 (S,Q)。 |
内模原理揭示了实现精确伺服控制的本质:控制器必须内部“模拟”外部世界的动态,以生成一个完全对消的信号。外部信号(参考或扰动)是一个持续的“激励流”。要使其对输出不产生净影响,系统必须产生一个与之大小相等、方向相反的“补偿流”。内模原理指出,产生这种精确补偿流的唯一可靠方法,是在控制器内部复制一个能生成该外部信号的“虚拟发生器”(内模)。这个内模在闭环中运行,其输出被反馈并与实际外部信号比较。由于内模的动态与外部信号相同,在稳态时,内模的状态会自发调整到恰好产生一个输出,该输出通过前向通道后,在系统输出端完全抵消外部信号的影响。例如,积分器 (1/s) 是阶跃信号 (1/s) 的内模,它能持续产生一个常值输出以抵消常值扰动。因此,内模原理是“以模型对抗模型”的控制哲学:在控制器内部植入外部环境的动态模型,利用模型间的动态匹配,在闭环内实现精确的对消,从而达成完美的跟踪或抑制。 |
(定理编号:TH-D1-0103 至 TH-D1-0107)
|
编号 |
类别 |
领域 |
定理/定律/原理名称 |
核心表述(公式/叙述) |
定理的逐步思考推理过程、证明思路及关键数学步骤 |
成立条件与适用范围 |
底层规律/理论归属 |
典型应用场景和工程意义 |
变量/参数列表及说明 |
数学与逻辑特征 |
定理的演绎/应用流程步骤 |
物理本质与“流”描述 |
|---|---|---|---|---|---|---|---|---|---|---|---|---|
|
TH-D1-0103 |
系统动力学 |
连续系统振动 |
铁木辛柯梁理论 |
考虑剪切变形和转动惯量的梁横向振动方程: |
1. 运动学关系:铁木辛柯梁理论中,截面转角 ψ由两部分组成:弯曲引起的转角 ∂wb/∂x和剪切变形引起的附加转角 γ。总转角 ψ=∂w/∂x−γ, 其中 γ是剪切角。 |
条件:梁的横截面在变形后仍保持平面,但不一定垂直于中性轴(考虑剪切变形);材料线弹性;小变形;考虑截面转动惯量。 |
弹性力学、达朗贝尔原理 |
场景:高速旋转的涡轮叶片、短粗的机床主轴、复合材料层合梁、微梁(MEMS)的动力学分析;地震工程中剪力墙的动力响应。 |
w(x,t):横向位移。 |
耦合方程组:两个二阶偏微分方程耦合在一起。 |
1. 建立方程:根据模型写出铁木辛柯梁的两个耦合运动方程。 |
铁木辛柯梁理论描述了梁在振动中“平动惯性流”、“转动惯性流”、“弯曲弹性恢复力流”和“剪切弹性恢复力流”之间的动态平衡与耦合。第一个方程是平动平衡:左侧 ρAw¨是平动惯性力流密度;右侧第一项 ∂x∂(κGA(∂w/∂x−ψ))是剪切力梯度,即由截面剪切变形产生的净弹性力流。第二个方程是转动平衡:左侧 ρIψ¨是转动惯性力矩流密度;右侧第一项 ∂x∂(EI∂ψ/∂x)是弯矩梯度,即由弯曲变形产生的净弹性力矩流;右侧第二项 κGA(∂w/∂x−ψ)是剪力产生的力矩。关键耦合项是剪切角 γ=∂w/∂x−ψ, 它连接了平动和转动方程。当剪切刚度 κGA→∞时,γ→0, 即 ψ→∂w/∂x, 第二个方程退化为 EI∂4w/∂x4+ρAw¨−ρI∂t2∂2(∂x2∂2w)=0, 即考虑转动惯量的 Rayleigh 梁方程;若进一步忽略转动惯量,则完全退化为欧拉-伯努利梁方程。该理论更完整地刻画了梁内两种能量存储(弯曲应变能、剪切应变能)和两种动能(平移动能、转动动能)的交互与传播。 |
|
TH-D1-0104 |
系统动力学 |
非线性动力学 |
同宿分岔与异宿分岔 |
在动力系统中,当参数变化时,若稳定流形 Ws和不稳定流形 Wu发生拓扑结构上的碰撞和横截相交,则称为全局分岔。具体地: |
1. 背景:在平面系统中,鞍点的稳定流形 Ws和不稳定流形 Wu通常是一维曲线。当参数变化时,这些流形的相对位置可能改变。 |
条件:系统存在双曲鞍点;参数变化导致稳定/不稳定流形发生拓扑连接;通常发生在非线性系统中。 |
动力系统分岔理论、微分方程几何理论 |
场景:Duffing 振子在双势阱中,当阻尼和激励参数变化时,可能出现同宿分岔,导致混沌运动;电力系统中发电机摇摆方程的同宿分岔分析;化学反应动力学中多稳态系统间的异宿连接;生态学模型中种群灭绝的临界现象。 |
μ:系统参数。 |
全局性:分岔涉及相空间中大范围轨道的拓扑变化,而不仅仅是平衡点或周期轨道的局部稳定性变化。 |
1. 寻找鞍点:对于给定的参数,找到系统的所有双曲鞍点平衡点。 |
同宿/异宿分岔描述了系统相空间中“不变流形”的全局连接与断裂,这直接改变了系统“状态流”的全局通道。在分岔前,鞍点的稳定流形 Ws和不稳定流形 Wu是分离的,它们像两条高速公路,引导着状态点的流向。同宿分岔发生时,Ws和 Wu恰好连接,形成一条封闭的“超级回路”:状态可以从鞍点沿 Wu出发,在相空间中漫游一大圈后,又精确地沿着 Ws回到同一个鞍点。这创造了一个无限长的 transient 过程。当参数越过临界值,这个精确的连接被打破,Ws和 Wu要么分离,要么发生横截相交。横截相交是动力学的“风暴眼”:它意味着 Wu会反复折叠并穿过 Ws, 在相交点附近形成一个极其复杂的、类似编织的结构(同宿缠结)。在这个缠结中,存在无穷多的周期轨道和敏感依赖于初值的轨迹,即混沌。因此,同宿/异宿分岔是系统从规则流动(如趋向平衡点或周期轨道)向湍流式、不可预测的混沌流动转变的关键阀门。 |
|
TH-D1-0105 |
系统动力学 |
控制理论 |
小增益定理 |
考虑两个子系统 H1和 H2的反馈互联。设 H1和 H2的输入输出增益(L2增益或其它范数定义的增益)分别为 γ1和 γ2。如果闭环系统是良定的,且满足小增益条件: |
1. 系统设置:考虑标准反馈结构:e1=u1−H2(e2), e2=u2+H1(e1)。 u1,u2是外部输入,e1,e2是子系统输入,y1=H1(e1), y2=H2(e2)是输出。 |
条件:子系统 H1和 H2是有限增益稳定的(即其增益 γ1,γ2有限);反馈互联是良定的(解存在唯一);小增益条件 γ1γ2<1成立。 |
泛函分析、输入输出稳定性理论 |
场景:鲁棒控制器设计(将模型不确定性建模为一个有界增益的扰动环节);网络控制系统稳定性分析(将网络延迟、丢包建模为动态环节);大系统分解与互联稳定性分析;非线性系统(如 Lur'e 系统)的绝对稳定性分析(可与圆判据联系)。 |
H1,H2:子系统(算子)。 |
输入输出特性:基于算子的增益,而非状态空间模型。 |
1. 系统建模:将待分析的系统分解为两个(或多个)子系统 H1和 H2的反馈互联形式。 |
小增益定理描述了一个反馈环路中“信号放大”能力的制约关系。每个子系统 Hi都可以看作一个“放大器”,其增益 γi描述了它最大能将输入信号的“强度”(用范数度量)放大多少倍。当两个这样的放大器连接成一个闭环时,信号会在环路中循环放大。小增益条件 γ1γ2<1意味着环路的总放大倍数小于1。这就像一个麦克风和扬声器组成的反馈系统:如果麦克风的增益 γ1和扬声器的增益 γ2的乘积小于1,那么任何初始的微小声音在环路中每循环一次就被衰减一次,最终会消失,系统稳定(不啸叫)。反之,如果乘积大于1,微小的噪声会被不断放大,导致信号强度趋于无穷(发散),系统不稳定(啸叫)。因此,小增益定理为反馈系统的稳定性设定了一个根本性的限制:环路的总体放大能力必须小于1,才能确保信号流不会在循环中无限增长。这是一种基于“能量”或“信号强度”流观点的稳定性判据。 |
|
TH-D1-0106 |
系统动力学 |
多体系统动力学 |
铰接体算法 |
一种用于计算树状结构多刚体系统正向动力学(已知力和力矩,求加速度)的高效递归算法。其核心思想是将系统视为由基体向外递归组合的“铰接体”,并利用铰接体惯性的概念。算法复杂度为 O(N), N为刚体数目。主要包含两个递归过程: |
1. 铰接体惯性定义:体 i的铰接体惯性 IiA定义为:将体 i及其所有子体视为一个刚性子树时,该子树在体 i铰接点处的等效惯性。它是一个 6×6的空间惯性矩阵。 |
条件:系统为树状拓扑结构(无闭环);铰是理想、无摩擦的;需要已知系统的运动学模型(连接关系、惯性参数)和所受外力。 |
图论、空间算子代数、牛顿-欧拉方程 |
场景:工业机器人、仿人机器人、机械臂的实时动力学计算与控制;虚拟现实和游戏中的角色动画;航天器多体附件(如太阳帆板、机械臂)的动力学仿真;车辆多体动力学(悬架、转向系统)。 |
IiA:体 i的铰接体惯性(6x6 空间矩阵)。 |
递归计算:算法分向内、向外递归遍历系统树。 |
1. 初始化:给定系统模型(拓扑、惯性、铰类型)、状态(q,q˙)和外力。 |
铰接体算法体现了多体系统动力学中“惯性”和“力”的递归传播与累积。它将整个多体系统视为由嵌套的“铰接体”组成。向内递归是一个“信息聚合”过程:从末梢开始,每个刚体将其自身的惯性 Ii和受到的外力,与所有子体的等效惯性(即子体的铰接体惯性 IjA)和力进行组合,形成以该体为根的子树的总等效惯性 IiA和总等效偏力 piA。这就像计算一个复杂机构的“有效质量”和“有效受力”,从末端传递到根部。向外递归则是一个“求解与传播”过程:从已知加速度的基体开始,利用根部已组装好的等效动力学方程,求解出第一个铰的加速度 q¨1;然后,这个加速度作为已知条件,结合下一级的等效动力学方程,求解出下一个铰的加速度,如此向外传播,直到末端。整个过程避免了构建和求解庞大的整体系统方程,而是通过局部递归计算,高效地完成了“力流”和“运动流”在整个树状网络中的正向和反向传播。 |
|
TH-D1-0107 |
系统动力学 |
分析力学 |
诺特定理 |
对于一个力学系统,如果其作用量 S=∫Ldt在某种连续变换下保持不变(即具有对称性),那么必然存在一个相应的守恒量。 |
1. 对称性假设:考虑依赖于小参数 ϵ的连续变换。作用量 S=∫t1t2Ldt在该变换下不变,即 δS=0。 |
条件:系统可由拉格朗日量描述;作用量在某个连续变换下具有不变性(对称性);变换是无穷小、连续的。 |
变分原理、群论 |
场景:分析力学系统的基本守恒律(能量、动量、角动量)的来源;在粒子物理中,由规范不变性导出电荷守恒;在广义相对论中,由微分同胚不变性导出能量-动量张量守恒;在工程中,分析复杂机械系统(如卫星、机器人)是否存在守恒量以简化模型。 |
L:系统的拉格朗日函数。 |
一一对应:连续对称性 ⇔守恒量。 |
1. 识别对称性:分析系统的拉格朗日量 L, 找出其在何种连续变换下,作用量 S保持不变(或仅改变边界项)。常见的对称性包括时间平移、空间平移、旋转等。 |
诺特定理揭示了物理系统的对称性(不变性)与守恒律之间深刻的“流”守恒关系。系统的每一种连续对称性,都对应着相空间中一种特定的“流”,而守恒量 C就是这个流的“通量”或“荷”。例如,时间平移对称性意味着物理规律不随时间原点变化。这产生了一个“时间流”,其对应的守恒量是能量 H。能量守恒意味着随着时间演化,系统在时间方向上的“平移流”是守恒的。空间平移对称性意味着物理规律在空间中是均匀的。这产生了一个“空间平移流”,其对应的守恒量是动量 p。动量守恒意味着系统在空间平移下的“流”是守恒的。空间旋转对称性意味着物理规律是各向同性的。这产生了一个“旋转流”,其对应的守恒量是角动量 L。诺特守恒量 C的表达式 ∑piζi−Hτ+F可以理解为:广义动量 pi沿着对称变换方向 ζi的分量之和,减去能量 H乘以时间变换 τ,再加上由于拉格朗日量变化产生的边界项 F。这个量在时间演化中保持不变,意味着对称变换所生成的“流”在动力学的“相流”下是守恒的。这是连接几何(对称性)与动力学(守恒律)的桥梁。 |
模块一:多体系统动力学与运动学基础 (公式 1-15)
此模块描述系统各部件的大范围刚体运动与约束关系,是控制的基础。
-
刚体质心运动方程:mir¨i=Fiext+∑Ficonstraint, 其中 mi为质量,ri为质心位置,Fiext为外力,Ficonstraint为约束力。
-
刚体绕质心转动方程 (欧拉方程):Iiω˙i+ωi×(Iiωi)=Tiext+∑Ticonstraint, Ii为惯性张量,ωi为角速度。
-
齐次变换矩阵:T=[R0Tp1], 描述坐标系间的旋转 R和平移 p,用于运动链正解。
-
Denavit-Hartenberg (D-H) 参数法:标准化串联运动链建模,关节 i到 i+1的变换:Tii+1=Rotz(θi)Transz(di)Transx(ai)Rotx(αi)。
-
运动学正解:末端位姿 Xend=T10T21⋯Tnn−1X0。
-
运动学逆解 (雅可比矩阵法):X˙=J(q)q˙, 其中 J为几何雅可比矩阵,q为关节坐标。
-
速度级逆解:q˙=J†X˙+(I−J†J)q˙0, J†为伪逆,(I−J†J)投影到零空间,q˙0用于优化。
-
加速度级逆解:q¨=J†(X¨−J˙q˙)。
-
闭链系统约束方程:Φ(q,t)=0。
-
约束方程的雅可比矩阵:Φqq˙+Φt=0, Φq为约束雅可比。
-
带拉格朗日乘子的系统方程:Mq¨+ΦqTλ=Q, λ为约束力乘子。
-
最小坐标形式方程:通过独立坐标 y表示 q, 消除约束力:M∗(y)y¨=Q∗(y,y˙)。
-
牛顿-欧拉递推动力学算法:从基座到末端递推计算速度/加速度,再从末端到基座递推算力/力矩,高效计算逆动力学。
-
计算力矩控制律:τ=M(q)(q¨d+Kve˙+Kpe)+C(q,q˙)q˙+G(q), 实现非线性解耦。
-
运动学标定模型:ΔX=JξΔξ, ξ为 D-H 参数误差向量,用于几何误差补偿。
模块二:柔性多体动力学与结构振动 (公式 16-35)
此模块描述由于结构柔性导致的变形、振动及其与刚体运动的耦合,是纳米精度的关键制约因素。
-
浮动坐标系下的柔性体运动方程:mttsymmtrmrrmtfmrfmffr¨ω˙q¨f+0000−ω~Iω000Kffrωqf=FtTrFf, qf为柔性模态坐标。
-
单元局部变形位移场:ue(x,t)=N(x)qfe(t), N为形函数矩阵。
-
单元质量矩阵:Me=∫VeρNTNdV。
-
单元刚度矩阵:Ke=∫VeBTDBdV, B为应变-位移矩阵,D为材料本构矩阵。
-
单元阻尼矩阵 (瑞利阻尼):Ce=αMe+βKe。
-
整体有限元组装:Mq¨+Cq˙+Kq=F。
-
模态坐标变换:q=Φη, Φ为模态振型矩阵,η为模态坐标。
-
模态质量、刚度、阻尼矩阵:Mm=ΦTMΦ=I, Km=ΦTKΦ=Ω2, Cm=ΦTCΦ=diag(2ζiωi)。
-
解耦的模态方程:η¨i+2ζiωiη˙i+ωi2ηi=ϕiTF。
-
残余柔度:Gr=K−1−∑i=1nmωi2ϕiϕiT, 描述被截断高阶模态的静态影响。
-
动态子结构法 (Craig-Bampton):界面坐标 qj与固定界面模态坐标 ηi:q=[qjqi]=[IΨij0Φi][qjηi]。
-
谐响应分析:(−ω2M+iωC+K)Q(ω)=F(ω)。
-
瞬态响应 (Newmark-β法):q¨t+Δt=a0(qt+Δt−qt)−a2q˙t−a3q¨t, q˙t+Δt=q˙t+a6q¨t+a7q¨t+Δt。
-
振动主动控制力:Fact=−Gvq˙−Gdq, Gv, Gd为速度/位移反馈增益。
-
压电作动器出力:Fpiezo=KaVa, Ka为作动器常数矩阵,Va为驱动电压。
-
应变传感器输出:Vs=Ksϵ, ϵ为应变。
-
基于模态观测器的状态估计:x^˙m=Amx^m+BmF+L(y−Cmx^m), xm=[η;η˙]。
-
LQR 振动控制:最小化 J=∫0∞(xmTQxm+uTRu)dt, 得最优控制 u=−Klqrxm。
-
频域振动抑制 (陷波滤波器):Hnotch(s)=s2+2ζn2ωns+ωn2s2+2ζn1ωns+ωn2, ζn2>ζn1, 用于抑制特定频率 ωn的振动。
-
振动功率流:P=FTq˙, 用于评估振动能量传递路径。
模块三:精密运动控制与伺服驱动 (公式 36-55)
此模块描述执行器、传感器与控制器构成的闭环系统,实现高精度轨迹跟踪。
-
永磁同步直线电机 (PMLSM) 推力方程:F=kfiq, kf为推力常数,iq为 q 轴电流。
-
PMLSM 电压方程 (d-q轴):ud=Rid+Lddtdid−ωeLqiq, uq=Riq+Lqdtdiq+ωe(Ldid+ψf)。
-
机械运动方程:mx¨=F−Ffric−Fdist。
-
三环伺服控制结构 (位置/速度/电流环):位置环输出速度指令,速度环输出电流指令,电流环控制电机电流。
-
PID 控制律:u(t)=Kpe(t)+Ki∫0te(τ)dτ+Kddtde(t)。
-
前馈控制:uff=Mx¨d+Cx˙d+G, 补偿已知动力学。
-
摩擦模型 (LuGre):dtdz=v−g(v)∣v∣z, Ffric=σ0z+σ1dtdz+σ2v, g(v)=Fc+(Fs−Fc)e−(v/vs)2。
-
摩擦前馈补偿:ufric_ff=F^fric(x˙d)。
-
重复控制:CRC(z)=1−Q(z)z−NQ(z)z−N, 用于抑制周期性误差,N为周期采样数。
-
迭代学习控制 (ILC):uk+1(t)=uk(t)+L∗ek(t+λ), k为迭代次数,L为学习增益,λ为超前时间。
-
扰动观测器 (DOB):d^=Q(s)[Fcmd−(Mns2+Cns)x], 估计并补偿扰动 d,Q(s)为低通滤波器。
-
加速度反馈:在速度环内引入加速度反馈,等效增加系统阻尼:Haf(s)=1+(Ka/Kv)s1。
-
位置环开环传递函数:Gol(s)=sKv⋅1+ωvs1⋅s2+2ζωns+ωn2ωn2, Kv为速度增益,ωv为速度环带宽,ωn为机械共振频率。
-
闭环带宽估算:ωc≈Kv(当机械谐振远高于 ωc时)。
-
跟踪误差估算 (恒速):ess=Kvv。
-
跟踪误差估算 (恒加速):ess=Kaa, Ka=Kv⋅ωv。
-
量化误差影响:位置分辨率 δx=2nL, L为量程,n为编码器位数。速度估计噪声 σv≈Tsδx32, Ts为采样周期。
-
延时补偿 (Smith Predictor):Csp(s)=1+C(s)G(s)(1−e−sTd)C(s), Td为延时。
-
自适应摩擦补偿:在线辨识 σ^0, σ^1, σ^2, Fc, Fs等参数。
-
基于干扰观测器的鲁棒控制:u=unom−d^/b0, b0为控制增益标称值。
模块四:热-力-流耦合动力学与误差补偿 (公式 56-70)
此模块描述温度场、应力场、流体场及其相互作用导致的时变误差。
-
瞬态热传导方程 (傅里叶定律):ρcp∂t∂T=∇⋅(k∇T)+q˙v, q˙v为内热源。
-
热对流边界条件:−k∂n∂T=h(T−T∞)。
-
热辐射边界条件:−k∂n∂T=ϵσ(T4−Tsurr4)。
-
热-结构耦合方程:Ku=Fth, Fth=∫VBTDϵthdV, ϵth=αΔT。
-
热致变形 (线性近似):ΔL=αLΔT。
-
主轴热漂移模型:δthermal(t)=∑iAi(1−e−t/τi), τi为热时间常数。
-
基于热电偶/热成像的温度场测量:T(x)=f(Vtc,Iir)。
-
热误差多元线性回归模型:δtherm=aTTkey+b, Tkey为关键点温度向量。
-
热误差神经网络补偿模型:δtherm=NN(Tkey,t,q˙...;W)。
-
计算流体动力学 (CFD) 控制方程 (N-S):∂t∂(ρu)+∇⋅(ρuu)=−∇p+∇⋅τ+ρg。
-
冷却液对流换热系数:Nu=kfhL=CRemPrn。
-
气浮导轨/主轴的压力分布 (雷诺方程):∂x∂(μph3∂x∂p)+∂y∂(μph3∂y∂p)=6U∂x∂(ph)+12∂t∂(ph)。
-
气膜刚度:Kair=−∂h∂F。
-
热-流-固强耦合迭代求解流程:[Tn+1,un+1,pn+1]=Solver(Tn,un,pn,BCs)。
-
实时热误差补偿值:ΔCcomp=−δtherm_predicted, 输入 CNC 系统。
模块五:多自由度协同与耦合控制 (公式 71-85)
此模块描述数千个自由度(如多轴机床、多机器人系统)间的协同、同步与解耦控制。
-
任务空间协同控制律:τ=JT(q)Ftask+N(q)τnull, Ftask为任务空间力,N投影到零空间。
-
主从同步控制 (位置偏差):τs=Kp(qm−qs)+Kd(q˙m−q˙s)。
-
交叉耦合控制 (CCC):对于双轴轮廓跟踪,定义轮廓误差 ε, 控制律修正各轴指令以最小化 ε。
-
轮廓误差估计 (直线/圆弧):对于直线路径,ε=(x˙d)2+(y˙d)2−exsinθ+eycosθ⋅(x˙dey−y˙dex)的近似。
-
自适应交叉耦合控制:uccc=Kc(t)ε, K˙c=ΓεeTJcT。
-
分散控制与耦合项视为扰动:对第 i个子系统,ui=Ci(s)(ri−yi), 耦合 di=∑j=igij(yj)。
-
集中控制 (状态反馈):u=−Kx, x为整个大系统状态,K通过 LQR 等求解。
-
分布式一致性协议:x˙i=ui, ui=∑j∈Niaij(xj−xi), 使所有 xi趋于一致。
-
包含领导者的协同:ui=∑j∈Niaij(xj−xi)+gi(x0−xi), x0为领导者状态。
-
基于图论的协同稳定性判据:系统矩阵 A−L⊗BK的特征值,L为拉普拉斯矩阵。
-
操作空间控制 (OSC):Γ=JTFtask+(I−JTJT†)τposture, 同时控制末端力与姿态。
-
零空间投影:τnull=(I−J†J)y, y为优化目标梯度,用于避障、能量优化等。
-
阻抗控制:Mde¨+Cde˙+Kde=Fext, 调节末端与环境交互的动态特性。
-
导纳控制:Mdx¨c+Cdx˙c+Kdxc=Fext, xc为位置指令修正量。
-
多速率控制:不同环路(位置、振动、热补偿)以不同采样率 fs1, fs2, fs3运行,通过插值/抽取同步。
模块六:鲁棒控制、自适应控制与智能控制 (公式 86-95)
此模块处理系统参数不确定性、非线性及时变特性,确保极端条件下的精度。
-
H∞ 控制标准问题:寻找控制器 K使闭环系统稳定且 ∥Tzw∥∞<γ, Tzw为扰动 w到性能输出 z的传递函数。
-
混合灵敏度 H∞ 设计:最小化 ∥[WSSWTT]∥∞, S=(I+GK)−1为灵敏度函数,T=I−S为补灵敏度函数。
-
μ 综合 (D-K迭代):处理结构不确定性 Δ, 最小化 ∥μΔ(M)∥∞, M为广义对象。
-
滑模控制 (SMC):定义滑模面 s=e˙+λe, 控制律 u=ueq−Ksgn(s), ueq为等效控制。
-
自适应控制 (模型参考 MRAC):对象 x˙p=Apxp+Bpu, 参考模型 x˙m=Amxm+Bmr, 控制律 u=ΘTω, 参数自适应律 Θ˙=−ΓωeTPBp。
-
自抗扰控制 (ADRC):u=b01(u0−f^), 其中 f^由扩张状态观测器 (ESO) 估计:z˙1=z2+β1(y−z1), z˙2=z3+β2(y−z1)+b0u, z˙3=β3(y−z1)。
-
模糊逻辑控制:u=∑i=1Nμi(e)∑i=1Nμi(e)⋅ui, μi为隶属度函数。
-
神经网络逆动力学控制:u=NNinv(qd,q˙d,q¨d;W), 网络训练逼近系统逆模型。
-
强化学习控制策略:策略 π(u∣s), 通过最大化累积奖励 R=∑γtrt来优化,s为状态,u为动作。
-
基于数据的预测控制 (MPC):在时域 N上,每一步求解优化问题:minu∑k=0N−1∥xk+1−xref∥Q2+∥uk∥R2, 施加动力学 xk+1=f(xk,uk)约束。
模块七:测量、传感与状态估计 (公式 96-100)
此模块描述纳米级信息的获取与融合,是闭环控制的前提。
-
激光干涉仪测长:ΔL=N⋅2nλ, N为条纹计数,λ为真空波长,n为空气折射率(需 Edlen 公式修正)。
-
电容传感器测微距:C=ϵ0ϵrdA, 通过测 C反推间隙 d。
-
多传感器数据融合 (卡尔曼滤波):预测:x^k∣k−1=Fkx^k−1∣k−1+Bkuk, Pk∣k−1=FkPk−1∣k−1FkT+Qk。更新:Kk=Pk∣k−1HkT(HkPk∣k−1HkT+Rk)−1, x^k∣k=x^k∣k−1+Kk(zk−Hkx^k∣k−1), Pk∣k=(I−KkHk)Pk∣k−1。
-
基于观测器的扰动估计:设计观测器 x^˙=Ax^+Bu+L(y−y^), d^=Vx^。
-
测量不确定度合成:Utotal=∑i(ciui)2, ci为灵敏度系数,ui为各不确定度分量。
从宏观运动规划、中观柔性振动与热变形控制、到微观纳米级误差补偿的完整技术体系。要驾驭数千自由度的超精密系统,关键在于模块化建模、多速率分层控制、以及基于多物理场模型的主动补偿。实际应用中,需根据具体机床结构(龙门式、并联式等)和工艺(车、铣、磨、光刻)从中选取和组合相应的公式群,构建专用的数字孪生与控制平台。
|
编号 |
类别 |
领域 |
定理/定律/原理名称 |
核心表述(公式/叙述) |
定理的逐步思考推理过程、证明思路及关键数学步骤 |
成立条件与适用范围 |
底层规律/理论归属 |
典型应用场景和工程意义 |
变量/参数列表及说明 |
数学与逻辑特征 |
定理的演绎/应用流程步骤 |
物理本质与“流”描述 |
|---|---|---|---|---|---|---|---|---|---|---|---|---|
|
TH-D1-0108 |
系统动力学 |
多体动力学 |
柔性多体系统动力学-浮动坐标系法 |
柔性体Bi的运动由大范围刚体运动(参考坐标系{Oi})和弹性变形(模态坐标qfi)叠加描述: |
1. 运动学描述:将柔性体上任意点P的位置分解为参考系运动与弹性变形的叠加,弹性变形用模态振型Φ展开。 |
条件:适用于小变形、线性弹性材料;模态截断合理;参考坐标系运动已知或与变形耦合求解。 |
计算多体动力学、结构动力学 |
场景:高速加工中心主轴-刀柄系统的颤振分析;工业机器人高速运动时的末端抖动预测;卫星太阳能帆板展开动力学;精密光学定位平台的微振动传递分析。 |
rOi:柔性体i参考坐标系原点在全局系中的位置。 |
强耦合非线性:方程高度非线性,源于刚体运动的大位移和姿态变化。 |
1. 柔性体建模:对每个柔性部件进行有限元分析,提取其低阶模态振型Φ和频率ω。 |
该模型描述了“刚体运动流”与“弹性变形波”的相互激发与耦合。参考坐标系定义了部件大范围运动的“主干流”。弹性变形是叠加在主干流上的“涟漪”或“波动”。当主干流加速或转向时(科氏加速度、向心加速度),其惯性力会“拍打”弹性体,激发出特定的波动模式(模态),这就是动力刚化效应。反之,这些弹性波动产生的动态力又会反作用于主干流,影响其运动轨迹。在超精密系统中,即使微小的弹性波动(纳米级)也会被敏感地反映在末端精度上。此模型的核心在于将这两种“流”统一在一个动力学框架内,精确计算它们的能量交换与相互调制。 |
|
TH-D1-0109 |
系统动力学 |
多体动力学 |
约束系统的微分-代数方程(DAE)指标-3问题 |
受完整约束的多体系统动力学可表述为指标-3的微分-代数方程组: |
1. 问题来源:约束方程Φ=0是对位置的约束。数值积分中,微小的舍入误差会导致解偏离约束流形,产生“违约”。 |
条件:系统受完整约束;质量矩阵M正定;约束雅可比Φq行满秩。 |
计算多体动力学、数值分析 |
场景:并联机床(如Stewart平台)的实时动力学仿真;复杂装配体(如汽车悬挂)的多体动力学分析;考虑关节间隙的机器人动力学建模。 |
q:系统的广义坐标向量。 |
混合系统:微分方程与代数方程紧密耦合。 |
1. 系统建模:建立系统的运动学模型,列出所有完整约束方程Φ(q,t)=0。 |
该问题揭示了物理世界“硬约束”与数值计算“软误差”之间的根本矛盾。物理约束(如两个部件必须始终通过铰链连接)定义了状态空间中的一个低维流形(约束流形)。系统的真实动力学轨迹严格位于此流形上。然而,数值积分过程如同一个视力模糊的登山者,试图沿着山脊(约束流形)行走。每走一步(积分步长)都会因为计算误差而稍微偏离山脊。对于指标-3的DAE,这种偏离不会被自动纠正,反而会累积放大,最终使登山者坠入山谷(解违约)。指标降阶和违约稳定化本质上是为登山者提供“纠偏反馈”:要么通过预测更准确的落脚点(利用加速度约束),要么在偏离时施加一个将其拉回山脊的力(Baumgarte稳定项)。在精密系统分析中,确保数值解严格满足约束,是准确计算关节反力、接触力的基础,这些力直接影响变形和精度。 |
|
TH-D1-0110 |
系统动力学 |
多体动力学 |
广义坐标分区法 |
对于受完整约束的系统,广义坐标q∈Rn, 独立约束m个,则可将q划分为独立坐标u∈Rn−m和非独立坐标v∈Rm, 即q=[u;v]。约束方程Φ(q)=0可(至少在局部)唯一确定v=g(u)。系统动力学可降阶为仅关于u的常微分方程: |
1. 隐函数定理:由约束方程Φ(q)=0和约束雅可比Φq行满秩,根据隐函数定理,可局部将m个坐标用其余n−m个坐标表示。 |
条件:约束是完整且独立的;划分q=[u;v]需保证对应的Φv非奇异(以保证v=g(u)存在);通常适用于开链或简单闭链系统。 |
分析力学、计算多体动力学 |
场景:串联机器人(如6轴关节型机器人)的实时动力学计算;简单平面机构(如曲柄滑块)的动力学仿真;作为其他高级算法(如计算力矩控制)的基础。 |
u:独立广义坐标向量(系统自由度)。 |
降维:将n维DAE降为(n−m)维ODE。 |
1. 坐标划分:根据系统拓扑结构,选择一组独立的广义坐标u。 |
该方法本质上是在约束流形上建立一套“局部坐标系”。独立坐标u就是这套局部坐标系的坐标。约束流形本身是弯曲的,但用u来描述时,我们是在其切空间里工作,从而将复杂的弯曲流形上的动力学,转化为平坦的切空间中的动力学。矩阵P就是连接流形与切空间的“映射”或“投影算子”。它将独立坐标空间中的运动“拉回”到全坐标空间。通过投影消去约束力,意味着我们只关心在流形切向上的动力学(由主动力和惯性决定),而自动满足了法向的约束条件。对于精密系统,这使我们能专注于自由度的本质动力学,避免处理数值上棘手的约束力,从而更清晰、更高效地分析影响精度的主要惯性力和耦合效应。 |
|
TH-D1-0111 |
系统动力学 |
多体动力学 |
旋量理论与对偶四元数 |
刚体的空间运动(螺旋运动)可用对偶四元数 q^=q+ϵq0紧凑表示,其中q为单位四元数表示旋转,q0=21tq, t为平移向量的四元数形式,ϵ为对偶单位(ϵ2=0)。空间速度(运动旋量)可表示为对偶四元数的导数:V^=2q^∗q^˙, 其中V^=ω+ϵ(v+ω×r)。动力学方程也可用此形式表达。 |
1. 旋量基础:空间速度(运动旋量)V=[ω;v]和空间力(力旋量)F=[f;τ]是描述刚体运动与受力的自然数学对象。 |
条件:适用于描述三维空间中的刚体运动;对偶四元数需满足单位约束。 |
旋量理论、几何代数 |
场景:并联机床(如Stewart平台)的正逆运动学求解与轨迹规划;多机器人协同搬运的协调控制;虚拟现实中的平滑位姿插值;基于视觉的精密位姿测量与伺服。 |
q^:对偶四元数,表示刚体位姿。 |
几何性:紧密贴合刚体运动的几何本质。 |
1. 位姿表示:将刚体的旋转矩阵R和平移向量p转换为对偶四元数q^。 |
该理论将刚体的空间运动视为绕空间某一轴线的螺旋运动(旋转+沿该轴的平移)。对偶四元数是描述这种螺旋运动的“语言”。它不像矩阵那样将旋转和平移分开处理,而是将其视为一个不可分割的几何实体。空间速度(运动旋量)V描述了刚体上所有点的速度场,如同一个“速度流场”。这个流场可以分解为绕某轴的旋转“涡流”ω和沿某方向的平移“均匀流”v。对偶四元数的微分关系q^˙=21q^V^, 正是描述了位姿“流”q^的变化率由瞬时速度“流”V^驱动。在超精密多轴运动控制中,使用这种表示法可以保证在大角度机动时,位姿插值和计算依然光滑、无奇异性,从而避免因数学表示缺陷引入的虚假误差。 |
|
TH-D1-0112 |
系统动力学 |
多体动力学 |
凯恩方程 |
对于具有n个自由度的系统,选择n个广义速率 ur(r=1,...,n), 它们与广义速度线性相关。定义偏速度 vr=∂ur∂v和偏角速度 ωr=∂ur∂ω。凯恩方程表述为: |
1. 达朗贝尔原理:凯恩方程源于达朗贝尔原理的投影形式。达朗贝尔原理:∑(Fi−miai)⋅δri=0。 |
条件:系统受理想约束;偏速度/偏角速度计算正确。 |
分析力学、多体动力学 |
场景:空间机械臂(多个刚体铰接)的动力学建模与实时控制;并联机器人的逆动力学计算;考虑非完整约束(如轮式机器人)的系统动力学。 |
ur:广义速率,是广义速度的线性组合。 |
速度空间投影:在速度空间而非位形空间工作。 |
1. 定义广义速率:根据系统自由度,选择一组线性无关的广义速率ur。 |
凯恩方程可以理解为在“速度流”的层面上直接建立动力学平衡。它将达朗贝尔原理中的“力平衡”投影到了由“偏速度”张成的“功率流通道”上。每个广义速率ur代表一个独立的运动模式(如某个关节的转速),其对应的偏速度vr,i描述了该运动模式在整个系统中引发的速度场分布。广义主动力Fr是外界输入给第r个运动模式的“功率流”。广义惯性力Fr∗是系统惯性抵抗该运动模式变化所“吸收”的负功率流。凯恩方程Fr∗+Fr=0意味着,对于每一个运动模式,输入功率与惯性吸收功率瞬时平衡。这种方法绕过了位形空间的复杂几何,直接在功率流层面进行结算,对于由多个部件通过复杂方式连接而成的系统,这种“分模式结算”的思路非常清晰和高效。 |
|
TH-D1-0113 |
系统动力学 |
多体动力学 |
虚功原理 |
对于一个处于平衡状态的力学系统,所有作用在系统上的力(包括约束力,如果是理想约束)在任何一组与约束相容的虚位移上所做的虚功之和为零: |
1. 平衡条件:对于静力学,平衡意味着合力为零。虚功原理是这一条件的积分(功)形式。 |
条件:适用于受理想约束的系统;虚位移必须与约束相容;原理适用于静力学和动力学。 |
分析力学、连续介质力学 |
场景:静力学结构分析(如桁架、梁);柔性多体系统动力学方程的推导;有限元法中单元刚度矩阵的推导;机器人静力学分析(如计算末端力与关节力矩的关系)。 |
Fi:作用在第i个质点或刚体上的主动力。 |
标量化:将矢量平衡问题转化为标量方程。 |
1. 确定系统与主动力:明确研究对象及所有主动力Fi。 |
虚功原理揭示了力学系统平衡的能量本质。它说:一个系统平衡,当且仅当任何微小的、允许的“试探性位移”都不会引起净的能量输入或输出。虚位移δri可以看作是在约束许可下,系统状态的一种“试探性波动”。主动力Fi在这个波动上做的功Fi⋅δri, 是外界试图通过这个波动输入系统的能量。原理要求所有这些能量输入的总和为零,意味着系统处于一种“能量稳定”状态——没有任何波动方向能让你白做功。对于动力系统,惯性力−miai的虚功代表了系统动能变化的负值。因此,动力学的虚功原理意味着:主动力输入的虚功率,等于系统动能变化的虚功率。这是一种瞬时功率平衡。这种方法的美妙在于,它完全避开了复杂约束力的具体方向和大小的分析,只关心在“运动可能的方向”上的能量交换,抓住了问题的核心。 |
|
TH-D1-0114 |
系统动力学 |
多体动力学 |
拉格朗日方程(第二类) |
对于具有n个自由度的完整系统,其运动由n个广义坐标 q1,...,qn描述。定义拉格朗日函数 L=T−V, 其中T为系统动能,V为系统势能。则系统的运动方程为: |
1. 从虚功原理出发:对于完整系统,结合达朗贝尔原理的虚功方程为:∑i(Fi−miai)⋅δri=0。 |
条件:系统必须是完整系统(约束可积,能完全用广义坐标描述);广义坐标独立;动能T和势能V需能明确表达。 |
分析力学 |
场景:多自由度机械振动系统(如机床床身模态分析)的建模;机器人动力学方程的推导;卫星姿态动力学建模;任何受完整约束的复杂机械系统。 |
qj:第j个广义坐标。 |
能量法:基于标量能量函数,而非矢量力。 |
1. 确定自由度与广义坐标:分析系统,选择一组独立的广义坐标qj。 |
拉格朗日方程是最小作用量原理(哈密顿原理) 的微分形式。哈密顿原理要求真实轨迹使作用量 S=∫t1t2Ldt取极值。对作用量变分,利用分部积分和端点固定条件,即可导出拉格朗日方程。物理上,∂q˙j∂L是广义动量,dtd项是其变化率;∂qj∂L是广义力(保守部分)。方程表明:广义动量的变化率,等于作用在系统上的总广义力(保守+非保守)。这实质上是牛顿第二定律在广义坐标下的推广。 |
|
编号 |
类别 |
领域 |
定理/定律/原理名称 |
核心表述(公式/叙述) |
定理的逐步思考推理过程、证明思路及关键数学步骤 |
成立条件与适用范围 |
底层规律/理论归属 |
典型应用场景和工程意义 |
变量/参数列表及说明 |
数学与逻辑特征 |
定理的演绎/应用流程步骤 |
物理本质与“流”描述 |
|---|---|---|---|---|---|---|---|---|---|---|---|---|
|
TH-D1-0115 |
系统动力学 |
结构动力学 |
Craig-Bampton 模态综合法 |
将复杂柔性结构划分为若干子结构。对于每个子结构,其物理坐标 u分为界面坐标 uj和内部坐标 ui。通过固定界面主模态 Φk和约束模态 Ψc进行坐标变换: |
1. 子结构划分:基于物理或几何边界,将大型结构划分为若干子结构。 |
条件:子结构内部线性;界面连接为刚性或可通过弹簧-阻尼单元模拟;固定界面模态和约束模态能充分表征子结构动态特性。 |
结构动力学, 模型降阶 |
场景:超精密机床整机动态特性分析;航天器太阳能帆板展开动力学仿真;汽车 NVH(噪声、振动与声振粗糙度)分析。 |
ui, uj:子结构内部和界面物理自由度坐标。 |
混合坐标:同时包含模态坐标(ηk)和物理坐标(uj)。 |
1. 子结构建模:对每个子结构进行有限元离散,得到 M, K, 并划分内部和界面自由度。 |
该方法将结构的动态响应视为“全局弹性波”与“局部静态变形”的叠加。固定界面主模态 Φk描述了子结构在“锚定”状态下的固有振动模式,如同被固定边界的鼓面振动模式。约束模态 Ψc则描述了当界面被“推动”时,子结构内部产生的准静态变形场,如同用手按压鼓面边缘时鼓面的整体凹陷。Craig-Bampton 基将这两种效应解耦并组合:模态坐标 ηk捕捉了内部的动态“涟漪”,而界面坐标 uj直接驱动了静态的“形状改变”。在整机动力学中,这种方法允许我们将复杂的结构“流”(振动波传播)分解为各子结构内部的局部动态“流”和通过界面传递的全局运动“流”,从而高效且高保真地模拟系统的整体动力学行为。 |
|
TH-D1-0116 |
系统动力学 |
接触力学 |
赫兹接触理论 |
两个弹性体在点或线接触时,由于局部变形,接触区域为椭圆(点接触)或矩形(线接触)。接触力 F与法向接近量 δ(两物体在接触点处的相对弹性趋近量)呈 3/2 次方关系: |
1. 基本假设:接触区域尺寸远小于物体曲率半径;材料各向同性、线弹性;接触表面光滑、无摩擦;变形为小变形。 |
条件:材料线弹性;接触区域小;表面光滑;无摩擦或滑动很小;变形远小于物体尺寸。 |
弹性力学, 接触力学 |
场景:精密轴承中球与滚道的接触刚度计算;齿轮齿面接触应力分析;三坐标测量机测头与工件的接触变形估计;微机电系统中微接触的力学分析。 |
F:法向接触力。 |
非线性:力-变形关系为 3/2 次幂律,刚度随载荷增大而增大(硬化)。 |
1. 确定接触几何:计算两接触体的主曲率半径,确定是点接触还是线接触,计算等效曲率半径 R。 |
赫兹接触描述了两个弹性体在微观尺度上的“拥抱”。当它们靠近时,并非在数学点上接触,而是由于弹性变形,形成一个微小的接触斑。接触力 F本质上是这个接触斑上分布压力的积分。法向趋近量 δ是两物体心“想要”靠近的距离,而材料的弹性抵抗了这种靠近。3/2 次方关系源于几何与弹性的耦合:随着载荷增加,接触斑面积(∝a2∝δ)增大,使得平均压力增长慢于载荷线性增长,从而导致刚度随载荷增加。在精密机械中,这种非线性刚度意味着系统的动态特性(如固有频率)会随着预紧力或工作载荷的变化而轻微漂移,是影响精度稳定性的一个因素。 |
|
TH-D1-0117 |
系统动力学 |
摩擦学 |
LuGre 摩擦模型 |
该模型用内部鬃毛(bristle)的平均变形 z来描述摩擦的微观预滑动行为。动力学方程为: |
v |
}{g(v)} z<br>g(v) = F_c + (F_s - F_c) e^{-(v/v_s)^\gamma}<br>F{fric} = \sigma_0 z + \sigma_1 \frac{dz}{dt} + \sigma_2 v<br>其中,v为相对滑动速度,F{fric}$ 为总摩擦力。 |
1. 鬃毛假设:将接触面建模为相互啮合的弹性鬃毛。相对运动导致鬃毛弯曲,产生摩擦力。当鬃毛变形超过一定限度时发生宏观滑动。 |
v |
}{g(v)} z项代表鬃毛的“松弛”或“断裂”,其速率与速度成正比。函数g(v)决定了稳态变形。<br>∗∗3.稳态特性∗∗:当\frac{dz}{dt}=0时,稳态变形z{ss} = g(v) \text{sgn}(v)。代入摩擦力公式,得稳态摩擦−速度曲线:F{ss}(v) = [\sigma_0 g(v) + \sigma_2 v] \text{sgn}(v)。通过设计g(v)可模拟Stribeck效应。<br>∗∗4.摩擦力构成∗∗:总摩擦力由三部分构成:鬃毛的弹性力\sigma_0 z(主导预滑动和静摩擦)、鬃毛变形的阻尼力\sigma_1 \frac{dz}{dt}(主导滞后和摩擦记忆)、粘性摩擦\sigma_2 v。<br>∗∗5.参数辨识∗∗:模型参数(F_c, F_s, v_s, \sigma_0, \sigma_1, \sigma_2, \gamma$) 需要通过实验数据(预滑动位移、摩擦滞后环、稳态曲线)拟合得到。 |
条件:适用于描述低速、高精度运动中的复杂摩擦现象,特别是预滑动位移、摩擦滞后、变静摩擦、Stribeck效应等。 |
摩擦学, 控制系统 |
场景:数控机床进给轴的低速爬行(stick-slip)现象分析与抑制;高精度转台/回转轴的静态摩擦补偿;机器人关节的精确力/位置控制;精密运动平台(如光刻机工件台)的微动定位控制。 |
|
TH-D1-0118 |
系统动力学 |
热力学/误差补偿 |
机床热误差建模 - 多元线性回归 |
机床关键点(如主轴端、刀尖点)的热致位移误差 ΔE与多个温度测点 Ti的测量值之间存在近似线性关系,可通过多元线性回归建立模型: |
1. 热源与热传递分析:识别机床主要热源(主轴电机、轴承、导轨丝杠、环境等)和热传递路径。 |
条件:假设热误差与温度场之间在测量范围内呈线性关系;温度测点能有效表征导致误差的主要热模态;测量数据足够多且覆盖典型工况。 |
统计学, 热误差工程 |
场景:加工中心、车床等数控机床的热误差实时补偿;三坐标测量机的温度补偿;精密光学设备的热稳定性提升。 |
ΔE:热致位移误差(例如,在 X, Y, Z 方向上的偏差)。 |
线性假设:模型形式简单,参数物理意义明确(βi可理解为热膨胀系数或结构柔度的某种组合)。 |
1. 实验设计:规划机床的加热实验,包括空载运行、切削负载、环境变化等多种工况。 |
该模型将复杂的热-结构耦合过程视为一个“黑箱”或“灰箱”。温度测点 Ti是黑箱的输入,热误差 ΔE是输出。多元线性回归试图找到一组权重 βi, 使得输入的线性组合能最佳拟合输出。物理上,βi综合反映了从第 i个热源到误差测量点之间的“热传递增益”和“机械转换增益”。前者取决于热传导/对流路径,后者取决于结构刚度与热膨胀系数。这种方法回避了求解复杂的三维瞬态热传导方程和热弹性方程,而是通过实验数据直接“学习”输入与输出之间的统计映射关系。其实时补偿的本质是构建一个与热变形“流”大小相等、方向相反的“补偿流”, 以抵消其对精度的影响。 |
|
TH-D1-0119 |
系统动力学 |
热力学/误差补偿 |
机床热误差建模 - 状态空间模型 |
将机床的热动态和热致变形视为一个动态系统,用状态空间方程描述: |
1. 热系统动态特性:机床热动态是一个分布参数系统,但可用集中参数的线性时不变(LTI)状态空间模型近似,以捕获其主要动态模态(如不同时间常数的升温过程)。 |
条件:假设热动态在辨识工况范围内可近似为线性时不变系统;输入输出数据信噪比足够高;激励信号能充分激发系统动态。 |
控制理论, 系统辨识 |
场景:高精度机床、光刻机工件台等对热动态补偿要求极高的设备;需要预测未来一段时间热误差的趋势,以进行前瞻性补偿。 |
x:状态向量,维度 n为模型阶次。 |
动态模型:能描述热过程的瞬态和稳态响应。 |
1. 实验设计与数据采集:设计能激发系统热动态的输入激励序列(如变转速、间歇运行),同步采集输入 u(t)和输出 y(t)数据。 |
状态空间模型将机床的热动态抽象为一个“热状态流”网络。状态向量 x代表了网络中各个“储能节点”的热能水平。矩阵 A描述了这些节点之间热能交换的规律(类似于热阻网络),决定了系统的自然热模态(快慢不同的升温过程)。矩阵 B描述了外部热源(输入 u)如何向这些节点注入“热流”。矩阵 C描述了如何从这些内部热能状态“读出”我们关心的可测量(温度、变形)。热误差是这种内部热状态流在机械结构上引发的“溢出”效应。通过辨识得到这个网络模型,我们就可以用少数传感器(输出)来反推整个网络的状态(观测器),并预测其未来的“溢出”(误差预测),从而实现精准的动态补偿。 |
|
TH-D1-0120 |
系统动力学 |
控制理论 |
自适应摩擦补偿 (基于 Lyapunov 稳定性) |
针对含摩擦的系统 x˙=v, mv˙=u−Ffric+d, 其中 Ffric用参数化模型(如 LuGre)描述。设计控制律 u和参数自适应律,使得跟踪误差 e=x−xd渐近收敛于零。以简化模型 Ffric=θTϕ(v)为例,其中 θ为未知参数向量,ϕ(v)为已知的回归函数向量。设计: |
1. 系统与误差方程:定义跟踪误差 e=x−xd, 其导数 e˙=v−x˙d。将系统动力学方程写为误差形式:me¨=u−Ffric+d−mx¨d。 |
条件:摩擦模型可参数化为线性参数化形式 Ffric=θTϕ(v);参数 θ为常数或慢变;扰动 d有界;回归函数 ϕ(v)已知。 |
自适应控制, 非线性控制 |
场景:高精度伺服系统中,摩擦参数随温度、润滑条件、磨损而缓慢变化时的在线补偿;机器人关节的摩擦自适应;精密光学元件的定位控制。 |
x,v:系统位置和速度。 |
Lyapunov 综合法:控制律和自适应律基于稳定性理论一体化设计。 |
1. 摩擦模型参数化:将采用的摩擦模型(如简化 LuGre, 库仑+粘性+Stribeck)写成线性参数化形式 Ffric=θTϕ(v)。 |
自适应摩擦补偿是一个“边干边学” 的过程。控制器最初基于一个粗略的摩擦参数估计 θ^(0)进行补偿。由于估计不准,会产生跟踪误差 e和 e˙。自适应律 θ^˙=−Γϕ(v)e˙像一个“误差驱动学习器”:当存在速度跟踪误差 e˙时,它就根据当前的运动状态 ϕ(v)按方向调整参数估计 θ^。调整的目的是使基于新估计的补偿力能更好地抵消真实摩擦力,从而减少未来的跟踪误差。Lyapunov 函数 V是衡量系统总“不理想程度”的能量函数,包含机械跟踪误差能量和参数估计误差能量。自适应律的设计保证了在控制作用下,这个总能量函数 V至少不会增加(V˙≤0),从而引导系统状态和参数估计共同趋向于理想平衡点。这相当于在控制“流”中引入了一个自适应的“摩擦补偿流”,该补偿流能自动调整其大小和形态,以匹配真实变化的摩擦“阻力流”。 |
|
TH-D1-0121 |
系统动力学 |
控制理论 |
迭代学习控制 (ILC) 基本算法 |
对于执行重复轨迹的任务,ILC 利用前一次迭代的跟踪误差信息来修正当前迭代的控制输入,以期在迭代域上逐步提高跟踪精度。基本 P 型 ILC 更新律为: |
1. 问题设定:系统在有限时间区间 t∈[0,T]内重复运行。每次迭代从相同的初始条件开始。系统动态可能不完全已知,但具有重复性。 |
1 - L(e^{j\omega}) G(e^{j\omega}) |
< 1对所有频率\omega成立。<br>∗∗4.超前时间\lambda∗∗:引入\lambda是为了补偿系统的动态延迟。对于因果系统,用e_k(t+\lambda)$ 意味着使用“未来”的误差信息来修正当前输入,这需要离线处理或存储上一批次的数据。 |
条件:系统运行任务必须是严格重复的(相同轨迹、相同初始条件);系统动态和扰动在迭代间具有可重复性;需要存储整条轨迹的误差和控制量数据。 |
学习控制, 迭代优化 |
场景:数控机床对复杂零件轮廓的重复精加工,以消除由伺服滞后、摩擦、背隙等引起的重复性轮廓误差;SCARA 机器人执行高速拾放操作;光盘驱动器光头定位。 |
k:迭代批次索引。 |
AtomGit 是由开放原子开源基金会联合 CSDN 等生态伙伴共同推出的新一代开源与人工智能协作平台。平台坚持“开放、中立、公益”的理念,把代码托管、模型共享、数据集托管、智能体开发体验和算力服务整合在一起,为开发者提供从开发、训练到部署的一站式体验。
更多推荐


所有评论(0)