第三章二自由度系统

二自由度系统第三章113.1

引言2kcm建模方法1：将车、人等全部作为一个质量考虑，并考虑弹性和阻尼。要求：对轿车的上下振动进行动力学建模。例子：轿车行驶在路面上会产生上下振动。缺点：模型粗糙，没有考虑人与车、车与车轮之间的相互影响。优点：模型简单；分析：人与车、车与车轮、车轮与地面之间的运动存在耦合。二自由度系统振动3k2c2m车m人k1c1建模方法2：车、人的质量分别考虑，并考虑各自的弹性和阻尼。优点：模型较为精确，考虑了人与车之间的耦合；缺点：没有考虑车与车轮之间的相互影响。二自由度系统振动4m人k1c1k2c2mk3c3k2c2k3c3m车m轮m轮建模方法3：车、人、车轮的质量分别考虑，并考虑各自的弹性和阻尼。优点：分别考虑了人与车、车与车轮之间的相互耦合，模型较为精确.问题：如何描述各个质量之间的相互耦合效应？二自由度系统振动53.2运动微分方程二自由度系统振动6运动微分方程例1：双质量弹簧系统，两质量分别受到激振力。不计摩擦和其他形式的阻尼。试建立系统的运动微分方程。二自由度系统振动/运动微分方程m1m2k3k1k2F1(t)F2(t)c1c2c37解：的原点分别取在的静平衡位置。

建立坐标：设某一瞬时：上分别有位移加速度受力分析：x1二自由度系统振动/运动微分方程速度F1(t)k1x1k2(x1-x2)m1k2(x2-x1)m2k3x2F2(t)x2m1m2k3k1k2F1(t)F2(t)c1c2c38建立方程：矩阵形式：F1(t)F2(t)二自由度系统振动/运动微分方程k1x1k2(x1-x2)m1k2(x2-x1)m2k3x29如同在单自由度系统中所定义的，在多自由度系统中也将质量、刚度、位移、加速度及力都理解为广义的。二自由度系统振动/运动微分方程m1m2k3k1k2F1(t)F2(t)c1c2c310可统一表示为：作用力方程位移向量加速度向量质量矩阵刚度矩阵激励力向量若系统有n个自由度，则各项皆为

n

维矩阵或列向量二自由度系统振动/运动微分方程位移向量刚度矩阵11式中：二自由度系统振动/运动微分方程121.质量矩阵的形成,系统的动能可以表示为

[M]即为所求的质量矩阵，显然为对称矩阵。二自由度系统振动/能量法用能量法确定振动系统的[M]、[K]、[C]13二自由度系统振动/能量法142.刚度矩阵的形成系统的势能可写为

[K]即为所求的刚度矩阵，也是对称矩阵。二自由度系统振动/能量法15二自由度系统振动/能量法163.阻尼矩阵的形成线性阻尼（黏滞阻尼）的耗能函数可写为[C]即为所求的阻尼矩阵，也是对称矩阵。二自由度系统振动/能量法17二自由度系统振动/能量法18例2：两自由度系统摆长

l，无质量，无阻尼做微摆动求：质量矩阵、刚度矩阵和运动微分方程xm1k1k2多自由度系统振动/多自由度系统的动力学方程19质量矩阵和刚度矩阵的正定性质n阶方阵A

正定并且等号仅在时才成立。

是指对于任意的

n维列向量y，总有成立。如果时，等号也成立，那么称矩阵A

是半正定的。

根据分析力学的结论，对于定常约束系统：动能：势能：二自由度系统振动/质量矩阵和刚度矩阵的正定性质20质量矩阵和刚度矩阵的正定性质n阶方阵A

正定并且等号仅在时才成立

是指对于任意的

n维列向量y，总有成立如果时，等号也成立，那么称矩阵A

是半正定的。

动能：除非所以，正定二自由度系统振动/质量矩阵和刚度矩阵的正定性质21质量矩阵和刚度矩阵的正定性质n阶方阵A

正定并且等号仅在时才成立

是指对于任意的

n维列向量y，总有成立如果时，等号也成立，那么称矩阵A

是半正定的。

势能：对于仅具有稳定平衡位置的系统，势能在平衡位置上取极小值。V>0当各个位移不全为零时，

K正定对于具有随遇平衡位置的系统，存在刚体位移。对于不全为零的位移存在V

＝0K半正定二自由度系统振动/质量矩阵和刚度矩阵的正定性质22振动问题中主要讨论K阵正定的系统及K阵半正定的系统，前者称为正定振动系统，后者称为半正定振动系统。二自由度系统振动/质量矩阵和刚度矩阵的正定性质23耦合与坐标变换矩阵中非零的非对角元元素称为耦合项。质量矩阵中出现耦合项称为惯性耦合。刚度矩阵或柔度矩阵中出现耦合项称为弹性耦合。以两自由度系统为例:不存在惯性耦合存在惯性耦合二自由度系统振动/耦合24如果系统仅在第一个坐标上产生加速度不出现惯性耦合时，一个坐标上产生的加速度只在该坐标上引起惯性力.同理，不出现弹性耦合时，一个坐标上产生的位移只在该坐标上引起弹性恢复力；而出现弹性耦合时，一个坐标上产生的位移还会在别的坐标上引起弹性恢复力.耦合的表现形式取决于坐标的选择耦合非耦合出现惯性耦合时，一个坐标上产生的加速度还会在别的坐标上引起惯性力.二自由度系统振动/质量矩阵和刚度矩阵的正定性质253.3不同坐标系下的运动微分方程26自由度与广义坐标

在任意坐标系中，要确定一个物体的位置所确定独立坐标的数目，称为这个物体的运动自由度。比如：在空间作任意运动的质点具有三个自由度；确定一个刚体在空间的位置，则需要六个参数，因而刚体作一般运动时具有六个运动自由度。

为了完全确定物体的位置而选定的任意一组彼此独立的坐标参数，称为这个物体的广义坐标。在选定坐标时，除去直角坐标Ｘ、Ｙ、Ｚ之外，我们也可以用角度φ、θ及从物体中的一点到某些固定点的距离等参数来确定物体在空间的位置。内容回顾27以汽车的二自由度振动模型为例汽车板簧以上部分被简化成为一根刚性杆，具有质量m和绕质心的转动惯量Ic。质心位于C点。分别在A点和B点与杆相联的弹性元件k1、k2为汽车的前，后板簧。Icm二自由度系统振动/不同坐标系的运动微分方程28采用能量表达式来建立其运动微分方程只考虑杆的竖向运动和绕质心的转动。系统的动能和势能为:只需要取定其中两个坐标，而将其他两个消去

二自由度系统振动/不同坐标系的运动微分方程291．取广义坐标为yA,q

yC和yB可用yA和q表示为:

在yA和q表示下系统的动能和势能为表示为:

二自由度系统振动/不同坐标系的运动微分方程30因此在

yA

和q坐标下的运动微分方程为:二自由度系统振动/不同坐标系的运动微分方程312．取广义坐标为yc和q

yA，yB用yc和q表示为为由yc，q到yA，yB的变换矩阵。

二自由度系统振动/不同坐标系的运动微分方程32在yc

，q下系统的动能和势能为[M][K]二自由度系统振动/不同坐标系的运动微分方程33运动微分方程为当时方程将存在弹性耦合。则刚度矩阵为对角矩阵方程已经解耦。这时系统垂直方向的运动与绕质心的转动独立。当[M][K]二自由度系统振动/不同坐标系的运动微分方程343.取广义坐标为yA

，

yB

yc和q可用yA和yB表示为：变换矩阵为:二自由度系统振动/不同坐标系的运动微分方程35在yA和yB下的质量矩阵为当时，方程存在惯性耦合。[M]在yA和yB下的的运动微分方程为二自由度系统振动/不同坐标系的运动微分方程36当时，质量矩阵为对角矩阵，刚度矩阵也为对角矩阵，运动微分方程已经解耦。此时有这时，A点和B点的振动相互独立，对于汽车来说，就是前悬架和后悬架的振动相互独立，系统就如同两个相互独立，没有联系的单自由度系统，在汽车理论中，称为悬挂质量分配系数。此时，系统的运动成为2个相互独立的单自由度系统。二自由度系统振动/不同坐标系的运动微分方程37如果广义坐标{x}和{y}之间有变换关系：在{x}，{y}下的刚度矩阵分别为[K]和[K1]，则由于系统势能大小与广义坐标的选取无关，有:从而得到坐标变换二自由度系统振动/不同坐标系的运动微分方程38同样，系统动能和能量耗散函数的大小也与广义坐标的选取无关，可以得到两个坐标系{x}和{y}的质量矩阵[M]、[M1]和阻尼矩阵[C]、[C1]之间的关系：{x}{y}二自由度系统振动/不同坐标系的运动微分方程39总结系统的质量矩阵、刚度矩阵及阻尼矩阵的具体形式/运动微分方程与所选取的描述系统振动的广义坐标有关，合适的广义坐标能够解除方程的耦合。由于不同广义坐标之间存在着线性变换关系，所以，方程解耦的问题就归结为寻找一个合适的线性变换矩阵[u]，使变换后系统的质量矩阵，阻尼矩阵和刚度矩阵成为对角矩阵，从而求得微分方程的解x。二自由度系统振动/不同坐标系的运动微分方程40研究汽车上下振动和俯仰振动的力学模型。表示车体的刚性杆AB的质量为m，杆绕质心C的转动惯量为Ic。悬挂弹簧和前后轮胎的弹性用刚度为k1和k2的两个弹簧来表示。写出车体微振动的微分方程。选取D点的垂直位移

和绕D点的角位移为坐标。ABCDa1a2el1l2lk1k2多自由度系统振动/多自由度系统的动力学方程41ABCDa1a2el1l2lk1k2ABCDa1a2el1l2lk1k2简化形式多自由度系统振动/多自由度系统的动力学方程42首先求刚度矩阵令：对D点取矩：力平衡：ABCDa1a2el1l2lk1k2CD多自由度系统振动/多自由度系统的动力学方程车体所受外力向D点简化为合力PD

和合力矩MD。微振动，杆质心的垂直位移、杆绕质心的角位移：43令：对D点取矩：力平衡：ABCDa1a2el1l2lk1k2CD刚度矩阵：多自由度系统振动/多自由度系统的动力学方程44求质量矩阵令：ABCDa1a2el1l2lk1k2CD惯性力质心C所受的惯性力：力平衡：力矩平衡：多自由度系统振动/多自由度系统的动力学方程45令：ABCDa1a2el1l2lk1k2质心C所受的惯性力矩：力平衡：对D点取矩：CD惯性力矩惯性力质心C所受的惯性力：质量矩阵：多自由度系统振动/多自由度系统的动力学方程46质量矩阵刚度矩阵运动微分方程多自由度系统振动/多自由度系统的动力学方程：作用在D点的外力合力和合力矩47如果D点选在这样一个特殊位置，使得：ABCDa1a2el1l2lk1k2只存在惯性耦合，而不出现弹性耦合多自由度系统振动/多自由度系统的动力学方程48如果D点选在质心C：ABCDa1a2el1l2lk1k2只存在弹性耦合，而不出现惯性耦合。多自由度系统振动/多自由度系统的动力学方程49问：能否找到这样一种坐标使得系统的运动微分方程既不出现惯性耦合，也不出现弹性耦合？即：若能够，则有：方程解耦，变成了两个单自由度问题。使系统运动微分方程的全部耦合项全部解耦的坐标称为主坐标。多自由度系统振动/多自由度系统的动力学方程50§3.4二自由度

无阻尼自由振动1二自由度系统振动/无阻尼自由振动51

不同广义坐标之间存在着线性变换关系，所以方程解耦的问题就归结为寻找一个合适的线性变换矩阵[u]，使变换后系统的质量矩阵，阻尼矩阵和刚度矩阵成为对角矩阵，从而求得微分方程的解。多自由振动运动方程：当即为无阻尼自由振动，则二自由度系统振动/无阻尼自由振动52

如果存在变换矩阵[u]使方程解耦。即当{x}=[u]{y}时，在{y}下的运动微分方程为：二自由度无阻尼自由振动的运动方程可写为，二自由度系统振动/无阻尼自由振动53其解为：如果初始条件为：则方程的解为：二自由度系统振动/无阻尼自由振动54由此得到方程的解：

可得系统的自由振动是简谐振动，即：二自由度系统振动/无阻尼自由振动55在特殊的初始条件下，方程能够解耦，使变换后系统的质量矩阵，阻尼矩阵和刚度矩阵成为对角矩阵，从而求得微分方程的解。此时，系统的两个自由度以相同的频率做简谐振动。同时达到极值，同时为零。它们之间的相位差为零或。

二自由度系统振动/无阻尼自由振动56例：无阻尼系统

设m1＝m2=m。这是个对称系统，对称点为k1的中点。

m1m2kk1k二自由度系统振动/无阻尼自由振动571.把m1，m2向右移动相同的距离x0，然后同时无初速度地放开，这时初始条件为：在整个振动过程中，弹簧k1不变形，m1和m2受到的力大小、方向均相同，二者的质量又相同，因此它们的速度和位移也相同。这样m1和m2之间的距离始终保持不变，二者就如同一个刚体。

讨论几种特殊的初始条件下的振动m1m2kk1k二自由度系统振动/无阻尼自由振动58此时等效为一个单自由度系统。

固有频率为：

响应为：并有：即此时系统两个自由度以为频率做简谐振动。同时达到极值，同时为零。它们之间的相位为零。m1m2kk1kkk二自由度系统振动/无阻尼自由振动592．m1向左，m2向右，均移动x0，然后同时无初速度地放开，这时初始条件为:由于系统的对称性，在振动过程中，系统的中点即k1的中点没有运动，就象一个固定点。k1被分成相等的两半，每一半的弹簧的刚度为2k1

m1m2kk1k二自由度系统振动/无阻尼自由振动60即此时系统两个自由度以为频率做简谐振动。同时达到极值，同时为零。它们之间的相位差为。响应为:固有频率:并有：m1m2kk1k2k1m1m2kk2k1二自由度系统振动/无阻尼自由振动613．如果初始条件为：即m1和m2的初始位移为零，而初始速度不为零。可得到：此时，与第一种初始条件下一样，系统两个自由度以为频率做简谐振动。同时达到极值，同时为零。它们之间的相位差为零。

二自由度系统振动/无阻尼自由振动624．如果初始条件为：

, 即m1和m2的初始位移为零，而初始速度不为零，方向相反大小相等，均为可得到：此时，与第二种初始条件下情况相同，系统两个自由度以为频率做简谐振动。同时达到极值，同时为零。它们之间的相位差为

m1m2kk2k1二自由度系统振动/无阻尼自由振动63对于任意的初始条件可分解为如下的四种初始条件之和：①②③④

①②③④二自由度系统振动/无阻尼自由振动64根据叠加原理，系统在任意初始条件下的自由振动响应为：或者写成：m1m2kk1k二自由度系统振动/无阻尼自由振动65

(1)

二自由度无阻尼系统在某些特定的初始条件下，其运动方程可解耦，此时自由振动是简谐振动。

(2)

系统的两个自由度以相同的频率振动，同时达到极值，同时为零，它们之间的相位差为零或，它们的坐标之比是与系统的物理参数有关而与时间无关的常数。

(3)

这种振动为系统的固有振动，固有振动时的频率称为系统的固有频率，坐标之比称为固有振型，简称振型，振型与固有频率是一一对应的。二自由度系统存在两种频率的固有振动，因此有两个固有频率，两个固有振型。二自由度系统在任意初始条件下的无阻尼自由振动是这两个固有振动的线性组合。二自由度系统振动/固有频率和固有阵型66设系统的固有振动时的解为：将上式代入方程式可以得到：

对这一线性齐次代数方程组。要想得到且A1，A2不全为零的解，方程系数矩阵行列式的值只能为零。

二自由度系统振动/固有频率和固有阵型67整理得：解之得到两个根：取正平方根即得到系统的两个固有频率：二自由度系统振动/固有频率和固有阵型68将w1代入方程式，得到：得到u11，u21的比值：这与我们前面利用系统对称性分析得到的结果也是一致的。

二自由度系统振动/固有频率和固有阵型69再将w2代入方程，得：这与我们前面利用系统对称性分析得到的结果同样也是一致的。二自由度系统振动/固有频率和固有阵型70设系统的运动微分方程为：系统的响应有如下形式：，这里：

二自由度系统振动/固有频率和固有阵型71代人方程：两边左乘振型的转置{u}T，并设：得：

上式为一个单自由度系统无阻尼自由振动方程，它的解为：

为简谐振动。

二自由度系统振动/固有频率和固有阵型72代入方程，对于线性齐次代数方程组。振型{u}有非零解的充分必要条件是系数矩阵的行列式为零，即有：称为特征方程或频率方程。

可以求出它们各自对应的振型{u1}＝{u1l，u2l}T，{u2}＝{u12，u22}T。可以看出，固有频率和它所对应的振型完全由质量矩阵和刚度矩阵决定，与外部激励无关，是系统固有的性质。二自由度系统振动/固有频率和固有阵型73取{ul}为振型矩阵[u]的第一列，{u2}为第二列，即得到，{u1}和，{u2}后，可以得到解：二自由度系统振动/固有频率和固有阵型74例：两自由度弹簧－质量系统m2m2kkkx1x2求：固有频率和主振型。多自由度系统振动/多自由度系统的自由振动/模态75解：动力学方程：令主振动：或直接用得：m2