振动理论及工程应用9-第十章-非线性振动课件.ppt

1、第第10章章 非线性振动非线性振动 一般来说，振动系统总是非线性的，线性系一般来说，振动系统总是非线性的，线性系统只是一种简单模型。如果线性理论能反映所要统只是一种简单模型。如果线性理论能反映所要考察的物理现象的定性性质和适当的定量结果，考察的物理现象的定性性质和适当的定量结果，那么就把它当作线性系统来处理；否则，就要研那么就把它当作线性系统来处理；否则，就要研究非线性系统。究非线性系统。在线性系统的研究中可以应用叠加原理，即在线性系统的研究中可以应用叠加原理，即系统对不同激励的响应可以线性相加，而对非线系统对不同激励的响应可以线性相加，而对非线性系统叠加原理不成立，因此对非线性系统的研性系统

2、叠加原理不成立，因此对非线性系统的研究比线性系统要复杂得多。究比线性系统要复杂得多。从研究方法上或是振动过程的变化规律上，从研究方法上或是振动过程的变化规律上，非线性振动与线性振动之间有本质区别。非线性振动与线性振动之间有本质区别。研究非线性振动有两种基本方法研究非线性振动有两种基本方法定性方法：定性方法：定性方法关心的是在已知解的邻域内系统的一定性方法关心的是在已知解的邻域内系统的一般稳定性特征，并非寻求与时间相关的解。般稳定性特征，并非寻求与时间相关的解。定量方法：定量方法：定量方法关心的是运动的时间历程，一般应用定量方法关心的是运动的时间历程，一般应用摄动法来求得这类方程的近似解析解。摄

3、动法来求得这类方程的近似解析解。10.1 非线性振动的例子非线性振动的例子 单摆的有限振幅振动是最简单的一个例子单摆的有限振幅振动是最简单的一个例子 0sin lg 063 lg 对于微小振动，对于微小振动，sin如果振幅不是很小如果振幅不是很小 运动微分方程为运动微分方程为 0 lg 线性系统线性系统非线性系统非线性系统rad/s/t0sin lg 单摆运动特性单摆运动特性 质量质量m在拉紧着的钢丝中的振动。设质量在拉紧着的钢丝中的振动。设质量m附着附着在长度为在长度为2l的钢丝中间，钢丝两端的拉力为的钢丝中间，钢丝两端的拉力为S。当质点。当质点从其平衡位置侧向移动距离从其平衡位置侧向移动距

4、离x时，钢丝产生恢复力，时，钢丝产生恢复力，运动微分方程为运动微分方程为 0sin2 llAESxm 其中其中A,E和和 l分别表示钢丝的横截面分别表示钢丝的横截面积，弹性模量和长度增量；积，弹性模量和长度增量；为钢丝为钢丝与竖直线的偏角。与竖直线的偏角。lxlxll2222 lxxlx 22sin 0233 xlAExlSxm 0sin2 llAESxm 运动微分方程为运动微分方程为 其中其中 如果不再假设位移如果不再假设位移x很小，那么弹簧的弹性恢复很小，那么弹簧的弹性恢复力一般地是位移力一般地是位移x的非线性函数的非线性函数 代入整理得代入整理得一般非线性系统的运动微分方程可表示为一般非

5、线性系统的运动微分方程可表示为 0 xFxm 0 xFx 0 xFx如果如果则称弹性恢复力为硬特性恢复力（称为硬弹簧）；则称弹性恢复力为硬特性恢复力（称为硬弹簧）；则称弹性恢复力为软特性恢复力（称为软弹簧）则称弹性恢复力为软特性恢复力（称为软弹簧）如果如果 0,3 xxxF0 0 例如例如当当时表示硬弹簧；时表示硬弹簧；时表示软弹簧。时表示软弹簧。当当硬弹簧曲线示意图硬弹簧曲线示意图 软弹簧曲线示意图软弹簧曲线示意图 tFxFxxm 如果系统还受到阻力强迫力的作用，则系统的运如果系统还受到阻力强迫力的作用，则系统的运动微分方程为动微分方程为 0,txxFxm 0,txxfx txxFmtxxf

6、,1,在一般情况下，单自由度系统的运动微分方程为在一般情况下，单自由度系统的运动微分方程为或或x 它是它是x和和的非线性函数。的非线性函数。其中其中 如果函数如果函数 f 不显含不显含t，则称这个系统为自治系统，则称这个系统为自治系统，否则称为非自治系统。否则称为非自治系统。10.2 相平面相平面 平衡点平衡点 0,xxfx yxfyyx,yxYyyxXx,式中式中x表示质点的位移，表示质点的位移，表示质点的速度。如表示质点的速度。如果把果把(x,y)看作平面上点的坐标（称为相点）看作平面上点的坐标（称为相点），该平，该平面称为相平面。面称为相平面。xy 设自治系统可表示为设自治系统可表示为或

7、或对于更一般的情形，方程可表示为对于更一般的情形，方程可表示为 微分方程式的一个解微分方程式的一个解x=x(t),y=y(t)对应于相平面对应于相平面上的一条曲线，称为相轨迹，简称轨迹。上的一条曲线，称为相轨迹，简称轨迹。0,0 yx 0,0,SSSSyxYyxX若相平面上的点为若相平面上的点为即即则称点则称点(xs,ys)为方程式的平衡点。为方程式的平衡点。设点设点O(xs,ys)是一个平衡点。令是一个平衡点。令 OOOOyYdxYcyXbxXa,yxYdycxyyxXbyaxx,11 dycxybyaxx ,02 qp 不妨设平衡点不妨设平衡点O为原点，则方程式可写成为原点，则方程式可写成

8、对于线性方程组对于线性方程组特征方程为特征方程为 两个特征根为两个特征根为 qppqpp4214212221 平衡点（平衡点（0，0）有如下类型：）有如下类型：（1）特征值均为负实数）特征值均为负实数(p0,p2 4q0)，则平衡，则平衡点是稳定结点点是稳定结点 21 21稳定结点稳定结点稳定非正常结点稳定非正常结点稳定星形结点稳定星形结点21 （2）两特征值均为正实数）两特征值均为正实数(p0,p2 4q0)，则平，则平衡点是不稳定结点。分别称为不稳定结点，不稳定非衡点是不稳定结点。分别称为不稳定结点，不稳定非正常结点和不稳定星形结点。图形分别与上图相似，正常结点和不稳定星形结点。图形分别与

9、上图相似，但箭头方向相反。但箭头方向相反。（3）特征值为相异实数（）特征值为相异实数（q0），则平衡点称为鞍），则平衡点称为鞍点，如图所示。点，如图所示。（4）特征值为复数，实部为负）特征值为复数，实部为负(p0,4qp2)，则，则平衡点称为稳定焦点，如图所示。平衡点称为稳定焦点，如图所示。稳定焦点稳定焦点鞍点鞍点 （5）特征值为复数，实部为正）特征值为复数，实部为正(p0,4qp2)，则平衡点称为不稳定焦点，此时形状与上图相同，则平衡点称为不稳定焦点，此时形状与上图相同，但箭头方向相反。但箭头方向相反。（6）特征值为纯虚数，则平衡点称为中心，此时）特征值为纯虚数，则平衡点称为中心，此时相迹为

10、封闭的圆，如图所示。相迹为封闭的圆，如图所示。中心点中心点02 qp 平平衡衡点点类类型型示示意意图图 总结以上各种情况，平衡点类型可在总结以上各种情况，平衡点类型可在pq平面平面上简单表示，如图所示。上简单表示，如图所示。例例1 设质量为设质量为m,长为长为l的单摆在具有粘性阻尼的介质的单摆在具有粘性阻尼的介质中运动中运动,阻尼系数为阻尼系数为c,其运动微分方程为其运动微分方程为0sin mgclml 试研究单摆运动的相图试研究单摆运动的相图.0202,mclg 解：解：令令0sin2200 则方程式可写成则方程式可写成 yx ,则方程式（则方程式（b）可表示为）可表示为再令再令 xyyxY

11、yyxXyxsin2,200 0,ynx,2,1,0n 0202,1,0 dcba022002 1 ,1202201 则方程式的平衡点为则方程式的平衡点为对于平衡点（对于平衡点（0，0），按式求得），按式求得特征方程为特征方程为特征值为特征值为 xyyxYyyxXyxsin2,200 当当 1 时，点（时，点（0，0）为稳定结点；）为稳定结点；对于平衡点对于平衡点(，0)，两特征值为相异实数，两特征值为相异实数，点点(，0)为鞍点为鞍点 1 ,1202201 当当 =1 时，点（时，点（0，0）为稳定的非正常结点；）为稳定的非正常结点；当当 1 时，点（时，点（0，0）为稳定焦点。）为稳定焦点

12、。其它平衡点可类似讨论。一般来说其它平衡点可类似讨论。一般来说平衡点（平衡点（0，0），），(2，0)，(4，0)为同类型平衡点。为同类型平衡点。大阻尼大阻尼小阻尼小阻尼临界阻尼临界阻尼 由摆的相图可见，摆的最低位置（由摆的相图可见，摆的最低位置（=0）是稳定的）是稳定的，而摆的最高位置（，而摆的最高位置（=）是不稳定的（同为鞍点）是不稳定的（同为鞍点）。）。摆的相平面示意图摆的相平面示意图稳定结点稳定结点非正常结点非正常结点稳定焦点稳定焦点（0，0）点点 1 =1 1大阻尼大阻尼小阻尼小阻尼临界阻尼临界阻尼10.3 保守系统保守系统 0 xfx ExVx 221保守系统的运动微分方程可以写成

13、保守系统的运动微分方程可以写成将上式积分得将上式积分得或或 设设x=0处势能为零，处势能为零，E是积分常数，它表示质是积分常数，它表示质点的机械能。表明系统的机械能守恒。点的机械能。表明系统的机械能守恒。EVT 保守系统中的周期运动对应于相平面上的闭保守系统中的周期运动对应于相平面上的闭轨线。闭轨线不是孤立的，这是保守系统的一个轨线。闭轨线不是孤立的，这是保守系统的一个特点。闭轨线一个包围一个，在相平面上充满某特点。闭轨线一个包围一个，在相平面上充满某个区域。从物理意义上来看，这意味着，如果有个区域。从物理意义上来看，这意味着，如果有一个周期运动，就有无穷多个周期运动。一个周期运动，就有无穷多

14、个周期运动。周期运动的周期为周期运动的周期为 mxxxVExT02d2 mxxVExT02d4 由于由于E依赖于初始条件，可见周期与初始条件有依赖于初始条件，可见周期与初始条件有关，而对于线性系统，周期与初始条件是无关的。关，而对于线性系统，周期与初始条件是无关的。这是非线性系统与线性系统的一个区别。这是非线性系统与线性系统的一个区别。例例10-3 求质量为求质量为m长为长为l的单摆，在作无阻尼自由振的单摆，在作无阻尼自由振动时的周期动时的周期T。解解 单摆的运动微分方程为单摆的运动微分方程为0sin lg E cos22020 lg 20 令令，并对式积分得，并对式积分得 AA 0为摆幅，则

15、为摆幅，则 AAAAT02200202sin2sind2coscos2d4 单摆振动得周期为单摆振动得周期为 sin2sin2sinA 2020sin2sin1d4 AT设设则上式可写成则上式可写成令令 12sin Ak，则上式展开后积分得，则上式展开后积分得 42020442206494112dsin83sin2114kkkkT 348122sinAAAk 42030721116112AAT 以以代入上式得代入上式得 若将式直接进行数值积分求其数值解，其时间历若将式直接进行数值积分求其数值解，其时间历程如图所示，它与解析结果表达式是基本一致的。由程如图所示，它与解析结果表达式是基本一致的。由

16、此可知，非线性振动的周期与振幅的大小有关。此可知，非线性振动的周期与振幅的大小有关。时间历程示意图时间历程示意图 4206494112kkT 即即10.4 非保守系统非保守系统 由于阻尼的存在，其总能量将随时间而变化，由于阻尼的存在，其总能量将随时间而变化，这一类系统称为非保守系统。这一类系统称为非保守系统。若系统的总能量随时间增加而单调减少，运若系统的总能量随时间增加而单调减少，运动最终趋向一个稳定的平衡位置，称为耗散系统。动最终趋向一个稳定的平衡位置，称为耗散系统。耗散系统中不存在周期运动。若在有阻尼的情耗散系统中不存在周期运动。若在有阻尼的情况下实现周期振动，必须不断地补充能量。况下实现

17、周期振动，必须不断地补充能量。若用外加周期力来补充能量，称为强迫振动。若用外加周期力来补充能量，称为强迫振动。若系统内部存在非振荡性能源，由它引起的周若系统内部存在非振荡性能源，由它引起的周期性振动，称为自激振动，简称自振。期性振动，称为自激振动，简称自振。10.4.1 极限环极限环 若在非线性系统中存在着封闭的相轨线，它对应若在非线性系统中存在着封闭的相轨线，它对应于系统的周期运动。于系统的周期运动。在保守系统中，封闭相轨线包围着奇数个奇点，在保守系统中，封闭相轨线包围着奇数个奇点，这些奇点是中心和鞍点。形成层层相套的连续系统。这些奇点是中心和鞍点。形成层层相套的连续系统。在非保守系统中，存

18、在着一条孤立的简单的封闭在非保守系统中，存在着一条孤立的简单的封闭相规线，称为极限环。相规线，称为极限环。极限环的振幅仅取决于系统的参数，而保守系统极限环的振幅仅取决于系统的参数，而保守系统封闭相轨线的振幅决定于初始能量。封闭相轨线的振幅决定于初始能量。在非保守系统中必须是经过一整周后，能量的增在非保守系统中必须是经过一整周后，能量的增量为零。量为零。稳定的极限环稳定的极限环不稳定的极限环图不稳定的极限环图 极限环极限环L有一个环形邻域，当随着时间的增有一个环形邻域，当随着时间的增长，所有轨线都从内部和外部无限趋近极限环长，所有轨线都从内部和外部无限趋近极限环L时，则称此极限环为稳定的，时，则

19、称此极限环为稳定的，反之，都远离极限环反之，都远离极限环L，则称此极限环为，则称此极限环为不稳定的，不稳定的，因此，极限环是否存在，是否唯一，是否稳定，以因此，极限环是否存在，是否唯一，是否稳定，以及其位置如何，都要根据具体问题认真分析讨论。及其位置如何，都要根据具体问题认真分析讨论。若属于若属于L外部的（内部的）趋于外部的（内部的）趋于L，内部的（，内部的（外部的）离开外部的）离开L，这种极限环称为半稳定的，这种极限环称为半稳定的，半稳定的极限环半稳定的极限环10.4.2 一个由摩擦引起的自激振动现象一个由摩擦引起的自激振动现象10 xSvr 若若 为皮带的速度，并为为皮带的速度，并为常数，

20、设常数，设x1为为M离开原点的离开原点的距离，则皮带相对于距离，则皮带相对于M的速的速度为度为0S 10 xSFvFr 1011xSFkxxm 摩擦力摩擦力F是皮带对物体的相对速度是皮带对物体的相对速度vr的非线性函数的非线性函数物块物块M的运动微分方程为的运动微分方程为设物块的平衡位置为设物块的平衡位置为x0，以以 代入上式得代入上式得000 xx 001SFkx kSFxxxx0101 0 kxxGxm xSgSgxG 00则函数则函数 的图形如图所示。的图形如图所示。xG 以以x0为新的坐标原点，其位移为为新的坐标原点，其位移为x，由坐标变换，由坐标变换则运动方程式可表示为则运动方程式可

21、表示为其中令其中令 若将此系统的运动在相平若将此系统的运动在相平面面(x，y)上表示，则由任意初上表示，则由任意初始状态确定的相点始状态确定的相点P0(x0,y0)出出发，绘制成相轨迹。发，绘制成相轨迹。它是具有恒定频率和恒定振幅的周期运动，称它是具有恒定频率和恒定振幅的周期运动，称为自激振动。工程中称为颤振。为自激振动。工程中称为颤振。由直线由直线P3P1和曲线和曲线P1P P3构成的封闭曲线构成的封闭曲线S，成为一个，成为一个稳定的极限环，稳定的极限环，10.4.3 中华文物龙洗的振动机理中华文物龙洗的振动机理 假设流体为无粘、无旋和不可压缩。在摩擦作用假设流体为无粘、无旋和不可压缩。在摩

22、擦作用下，壳体的振动很小，完全在线性范围内。设流体下，壳体的振动很小，完全在线性范围内。设流体的运动也属于线性范围，所以流体速度势表示为的运动也属于线性范围，所以流体速度势表示为 忽略厚度的不均匀忽略厚度的不均匀性，将龙洗简化为一个性，将龙洗简化为一个铜制的近似旋转薄壳，铜制的近似旋转薄壳，靠近边缘有两个供手搓靠近边缘有两个供手搓动摩擦的耳。动摩擦的耳。tTzrtzr,0,2 zr trgftHz,nwn zr,由势流理论，由势流理论，应满足应满足Laplace方程，即方程，即流场的边界条件是流场的边界条件是 00,0,zzrrA 00,zzrrB 把两耳上的搓动视为在壳上两点把两耳上的搓动视

23、为在壳上两点上的搓动。并搓动的速度相同上的搓动。并搓动的速度相同 搓动速度为搓动速度为0S，3111sgnrrrrvcvbvavF 干摩擦力函数为干摩擦力函数为 0,221 tuhwvuL 000000222 ,0,z,z,rrVSFNzzrrVSFNtvhwvuLBBAA tTntwhwvuL 0223,壳体的振动方程为壳体的振动方程为 ,izrWzrVzrUUiiiTi21 ,1,iiitTUtzru 21 2,itPtTtTiiii 将耦合系统的湿模态，记为将耦合系统的湿模态，记为方程式的解表示为模态的线性组合方程式的解表示为模态的线性组合由湿模态的正交性，模态坐标下的运动方程为由湿模态

24、的正交性，模态坐标下的运动方程为 wvuu,T 振动位移表示为振动位移表示为 其中其中Pi(t)为模态广义力。为模态广义力。平面是镜面对称的，模态可分为三类：平面是镜面对称的，模态可分为三类：（）对）对xz平面对称；平面对称；（）对）对xz平面反对称，对平面反对称，对yz平面也反对称；平面也反对称；（）对）对xz平面反对称，对平面反对称，对yz平面对称。平面对称。21,0I,itPi 21 sgn2II3010113IIII01IIII01IIII01II,iVScSbaTVS c TVSbTVSaNtPiiiiiiiAi 21 sgn2IIIi3010113IIIiIII01IIIIII01

25、IIIIII01III,iVScSbaTVS c TVSbTVSaNtPiiiiiAi 对于第对于第I类，模态在类，模态在A、B点处周向位移为零。点处周向位移为零。设设NA=NB，计算非线性模态广义力得计算非线性模态广义力得第一种情形，双耳同向搓动第一种情形，双耳同向搓动第二种情形，双耳反向搓动第二种情形，双耳反向搓动 模态广义力是广义速度的非线性函数，此系统是模态广义力是广义速度的非线性函数，此系统是一类多自由度自激振动系统。一类多自由度自激振动系统。第一种情形只产生第第一种情形只产生第II类模态的自激振动；类模态的自激振动；第二种情形只产生第第二种情形只产生第III类模态的自激振动。类模态

26、的自激振动。方程式是多自由度的非线性振动方程。直接求解析方程式是多自由度的非线性振动方程。直接求解析解较困难，首先以第一种情形为例进行数值解。首解较困难，首先以第一种情形为例进行数值解。首先对时间先对时间t进行坐标变换，引入新变量进行坐标变换，引入新变量，设，设1II1IIdtd t 21II22221IIdddd dddd TtTTtT 2II30101131112IIII21II30101131111IIII1sgn2sgn2VScSbaScSbSaNKTTVScSbaScSbSaNTTrrrArrrA 2II2II1II1II0TKVTVSSr 21II22II K269813.01II

27、 V0448042.02II V15.8531II 代入方程式，取前两阶湿模态，可得代入方程式，取前两阶湿模态，可得其中其中由模态实验结果，湿模态参数为由模态实验结果，湿模态参数为68.31692II 计算结果（相平面图）如图所示。计算结果（相平面图）如图所示。第一阶模态的数值均在第二阶模态相应数值的第一阶模态的数值均在第二阶模态相应数值的3以以下，并且都是稳定解。结论：在此多自由度非线性振下，并且都是稳定解。结论：在此多自由度非线性振动系统中，自激振动主要以同类最低阶湿模态为主，动系统中，自激振动主要以同类最低阶湿模态为主，高阶湿模态可忽略不计。高阶湿模态可忽略不计。0100SS 010SS

28、 01S 振幅的最大值发生在振幅的最大值发生在 点。点。以搓动速度为横坐以搓动速度为横坐标，振幅为纵坐标，标，振幅为纵坐标，结果曲线如图所示。结果曲线如图所示。它表示当粘滑系数一它表示当粘滑系数一定时振幅与搓动速度定时振幅与搓动速度的变化关系的变化关系0100SS 时，振幅随搓动速度的增大而增大时，振幅随搓动速度的增大而增大010SS 时，振幅随搓动速度的增大而减小时，振幅随搓动速度的增大而减小 定常解振幅值与粘滑系数和搓动速度有关。要达到定常解振幅值与粘滑系数和搓动速度有关。要达到较大值，使水珠喷得高，不但调整搓动速度，还要选较大值，使水珠喷得高，不但调整搓动速度，还要选择适当的粘滑系数。（

29、在搓动之前要用肥皂洗手）择适当的粘滑系数。（在搓动之前要用肥皂洗手）第一种情形，龙洗振动表现的是第第一种情形，龙洗振动表现的是第类中第一阶模类中第一阶模态。它有四个振动的波峰处，此处的振动最大，而与之态。它有四个振动的波峰处，此处的振动最大，而与之交错的节线将水面分成四部分，因此就形成了水珠在四交错的节线将水面分成四部分，因此就形成了水珠在四个部位被喷起的景观。个部位被喷起的景观。第二种情形，龙洗振动表现的是第第二种情形，龙洗振动表现的是第类第一阶模态类第一阶模态，振动的六个波峰与节线相互交错将水面分成六个区，振动的六个波峰与节线相互交错将水面分成六个区域，形成了水珠在六个部位被喷起的壮观场景

30、。域，形成了水珠在六个部位被喷起的壮观场景。10.5 摄动法摄动法 众所周知，只有为数很少的非线性微分方程可求得精众所周知，只有为数很少的非线性微分方程可求得精确的解析解。在一般情况下，只能用近似方法求解，将微确的解析解。在一般情况下，只能用近似方法求解，将微分方程的解近似地表为幂级数形式，其近似程度取决于级分方程的解近似地表为幂级数形式，其近似程度取决于级数解中所取得的项数。数解中所取得的项数。如果表示系统的微分方程中的非线性项是个小量，如果表示系统的微分方程中的非线性项是个小量，那末这个小的非线性项就称为摄动。摄动项一般用一个那末这个小的非线性项就称为摄动。摄动项一般用一个小参数小参数 标

31、出。我们通常是求小参数标出。我们通常是求小参数 的幂级数形式的解的幂级数形式的解。这种求解方法统称为摄动法，或称小参数法。这种求解方法统称为摄动法，或称小参数法。非线性自治系统的方程为),(20 xxfxx 派生系统为020 xx)()()(),(2210txtxtxtx设非线性自治系统的方程为 代入方程，由于等号两端的同幂项的系数必须相等，于是整理后有 xxxfxxxxfxxxxxfxxxx),(),(),(000100122020012010200由此，得到一系列关于函数的微分方程。这样，解非线性微分方程就代之以解一系列线性微分方程，而且可以递归求解。可得精确到n的近似解)(12210nn

32、nOxxxxx在应用基本摄动法时，在解中经常出现久期项，若只保留前几项，它将随时间无限增长，但是，如果取无限多项，那末级数必定是收敛的，但我们又不能求出无限多项。10.6 平均法平均法 在一般情况下，非线性振动微分方程很难求得精在一般情况下，非线性振动微分方程很难求得精确的解析解，只能用近似方法求得近似解析解，平均确的解析解，只能用近似方法求得近似解析解，平均法就是的其中一种有效的方法。法就是的其中一种有效的方法。所谓平均法就是将以位移为未知量的振动方程，所谓平均法就是将以位移为未知量的振动方程，化成以振幅、相位为未知量的标准方程组，因为振幅化成以振幅、相位为未知量的标准方程组，因为振幅和相位

33、的导数都是和相位的导数都是O()量级的周期函数，因此，可用量级的周期函数，因此，可用一个周期的平均值代替它，故称其为平均法。一个周期的平均值代替它，故称其为平均法。平均法有多种形式，下面以平均法有多种形式，下面以KBM法为例加以说明。法为例加以说明。10.6.1 自治系统的平均法自治系统的平均法一个自由度系统的自由振动方程为一个自由度系统的自由振动方程为其中其中 为正的小参数；为正的小参数；x为振动位移；为振动位移；txxfdd,为非线性函数。为非线性函数。如如=0，可得其解为，可得其解为 t其中其中a、为起始条件确定的常数为起始条件确定的常数 txxfxtxdd,dd222 （1）sinco

34、saxax（2）当当0时，即有非线性干扰力存在，则时，即有非线性干扰力存在，则a、将为将为时间时间t的函数，现研究是什么函数时的函数，现研究是什么函数时 sinddcosdddd tatatx微分式（微分式（2）的第一式）的第一式令上式与式（令上式与式（2）的第二式相等，则有）的第二式相等，则有 sincosaxax txxfxtxdd,dd222 满足满足（1）（2）0sinddcosdd tata（3）cosddsindddd22 tatatx sin,coscosddsinddaaftata 微分式（微分式（2）的第二式）的第二式代入式（代入式（1），则），则 （4）从从(3)、(4)中

35、解出中解出称为标准方程组称为标准方程组 ,sinsin,cosddaaafta ,sinsin,cosddaaafat （5 5）设函数设函数a、由平稳变化的项由平稳变化的项y、和小振动相叠和小振动相叠加而成，取加而成，取KB变换为变换为 ,221221ytVytVytUytUya yZyZyZtyYyYyYty3322133221dddd 并要求新变量并要求新变量y、的导数为的导数为其中其中Y1、Y2、Y3、Z1、Z2、Z3不显含不显含t，U1、U2、V1、V2为为 的以的以2 为周期的周期函数及为为周期的周期函数及为t的以的以T为周期的为周期的周期函数。周期函数。将上式代入式（将上式代入式

36、（5 5）得）得 1203322123322122223322113322111332211203323122332212222332211332211133221,ytZZZVYYYyVtVZZZVYYYyVtVZZZytZZZUYYYyUtUZZZUYYYyUtUYYY 011011 tVZtUY 12111121211112 tVZVYyVZtUZUYyUY其中其中 为将为将 在在=0点展成泰勒点展成泰勒级数时级数时 的系数。的系数。1010,令上面等式令上面等式两端两端 同次方的系数相等同次方的系数相等，则可确定，则可确定Y、Z、V、U的微分方程式，下面只写出的微分方程式，下面只写出

37、前二阶的系数方程前二阶的系数方程 TTtTZtTY020001020001d21d1d21d1 根据根据Y、Z不显含不显含t的条件，可用下列方法确定的条件，可用下列方法确定和和其中其中T为振动周期。为振动周期。TTtZVYyVTZtZUYyUTY01111120111112d1d1 ttVtZVUtYU010101010111dd ttVtZUYyUZVUtZUYyUYU02011112020201111212dd 而而和和在运算过程中，将不显含在运算过程中，将不显含t t的函数当作常数。的函数当作常数。yZtyYtytyx11ddddcos ,cos11ytVytUyatax第一次近似解第一

38、次近似解第二次近似解第二次近似解 ttVtZytVUtYytU010101010111d,d,yZyZtyYyYty221221dddd 其中其中 和和 由此可知，如将式写成由此可知，如将式写成 的的m阶级数，当有关阶级数，当有关的系数确定后，则可求得其的系数确定后，则可求得其m次近似解。次近似解。vtEtxxfxtxsindd,dd1222 如果在振动系统上作用有外干扰力如果在振动系统上作用有外干扰力Esinvt，则单自由度非线性系统强迫振动的方程式则单自由度非线性系统强迫振动的方程式10.6.2 非共振情况的平均法非共振情况的平均法在线性系统中，定常解的频率与干扰频率相同，在线性系统中，定

39、常解的频率与干扰频率相同，并且定常解和起始条件无关。在非线性系统中，除并且定常解和起始条件无关。在非线性系统中，除含有和干扰力相同的频率成分外，还有高频率成分含有和干扰力相同的频率成分外，还有高频率成分存在，并且与起始条件有关。存在，并且与起始条件有关。如果干扰力频率和派生系统（当如果干扰力频率和派生系统（当=0时的系统时的系统）的固有频率差值较大，则为非共振情况。）的固有频率差值较大，则为非共振情况。当当=0时，设其解为时，设其解为 vtvEtatxvtvEtaxcossinddsincos2222 其中其中a,为任意常数。为任意常数。如果如果0，强迫振动的解仍取上式的形式，则，强迫振动的解

40、仍取上式的形式，则其中的其中的a,将为时间将为时间t的函数，现以的函数，现以a,作为新变量，作为新变量，进行变量转换。进行变量转换。通过类似与自治系统的平均法的计算，即可求通过类似与自治系统的平均法的计算，即可求出非共振情况的近似解析解。基本思路与以前相同出非共振情况的近似解析解。基本思路与以前相同。可参考有关书籍。可参考有关书籍。10.6.3 共振情况的平均法共振情况的平均法 txxvtfvtEtxxfxtxdd,sindd,dd1222 在单自由度非线性系统强迫振动的方程式中，在单自由度非线性系统强迫振动的方程式中，如果如果v与与 相等或二者的差值为与相等或二者的差值为与 同价的量，则为共

41、同价的量，则为共振情况。振情况。在共振情况下实际上是有周期解的，此时干扰在共振情况下实际上是有周期解的，此时干扰力的幅值与力的幅值与 同量级，为此，方程式为同量级，为此，方程式为下面讨论上述方程在共振情况下的渐近解。下面讨论上述方程在共振情况下的渐近解。设设 是是vt的以的以2 为周期的周期函数，为周期的周期函数，可表示为可表示为 txxvtfdd,NNnnvtitxxfetxxvtfdd,dd,n taxcos若无干扰，即若无干扰，即=0，则系统将有简谐振动，则系统将有简谐振动 将其代入上式将其代入上式,并将展成傅氏级数，由于它对并将展成傅氏级数，由于它对vt来说是周期函数，所以其傅氏级数中

42、将有来说是周期函数，所以其傅氏级数中将有 sin(nv+m)t和和cos(nv+m)t，其中，其中n和和m是整数是整数 在干扰力中将有组合频率为在干扰力中将有组合频率为nv+m 的谐波成分。的谐波成分。很明显，任何一个组合频率成分与固有频率相很明显，任何一个组合频率成分与固有频率相接近时，则具有该频率成分的干扰力对振动特性将接近时，则具有该频率成分的干扰力对振动特性将发生很显著的影响。发生很显著的影响。因此在非线性系统中，共振现象不但发生在因此在非线性系统中，共振现象不但发生在 v时，而且也发生在时，而且也发生在nv+m 时，故在单频干扰力时，故在单频干扰力作用下，共振关系可近似取为作用下，共

43、振关系可近似取为 qpv 几种共振情况：几种共振情况：2、nv 时，时，称为超谐共振；称为超谐共振；1、v 时，时，称为主共振；称为主共振；3、v n 时，时，称为亚谐共振称为亚谐共振 由于共振的类别很多，所以在研究共振情况的由于共振的类别很多，所以在研究共振情况的渐进解时，应指明研究哪种共振类型的渐进解。渐进解时，应指明研究哪种共振类型的渐进解。4、v p/q时，称为分数共振时，称为分数共振 qpv pqv以分数共振为例以分数共振为例 tpqvatxtpqvaxsinddcos 称为调谐参数，称为调谐参数，设解为设解为 以以a、作为新变量进行坐标变换，通过类似与作为新变量进行坐标变换，通过类

44、似与自治系统的平均法的计算，即可求出共振情况的近似自治系统的平均法的计算，即可求出共振情况的近似解析解。基本思路与以前相同。可参考有关书籍。解析解。基本思路与以前相同。可参考有关书籍。10.6.4 龙洗自激振动的近似解析解龙洗自激振动的近似解析解 ArrrrNvcvbvavF)sgn(3111*作用于龙洗两耳上的干摩擦力和相对速度的作用于龙洗两耳上的干摩擦力和相对速度的关系可近似表示为：关系可近似表示为：其中其中 a1 为静摩擦系数，为静摩擦系数，b1、c1为实验曲线拟合的为实验曲线拟合的动摩擦粘滑系数。动摩擦粘滑系数。NA为作用于龙洗双耳上为作用于龙洗双耳上A点的正压力，点的正压力，此函数是

45、由三段光滑非线性曲线组成。此函数是由三段光滑非线性曲线组成。)(1II0tTVSvr 1 II1 II0*2)(2)()(VtTVSFtTtT 1 II0*2)(12)(VSFtTX 参考数值计算结果，略去高阶模态，只考虑参考数值计算结果，略去高阶模态，只考虑龙洗的第一阶模态，则相对速度为龙洗的第一阶模态，则相对速度为龙洗的振动方程为龙洗的振动方程为引入新变量引入新变量 由于当摩擦力函数为跳跃的非确定性函数，由于当摩擦力函数为跳跃的非确定性函数，对求解此系统的解析解增加了很大困难，对求解此系统的解析解增加了很大困难，t)()()(XFXX II13010113II101II101II101 a

46、)()()(sgn(2)(VScSbXVScXVSbXVSaNXFA 引入新变量引入新变量 进行无量纲化后，运动微进行无量纲化后，运动微分方程式为分方程式为其中其中)()()(12XFtXtX 1 II0*1 II1 II0*1)(2)(2)(VSFVXVSFXF 运动微分方程式为运动微分方程式为其中其中 利用数值解，相平面图如图所示。利用数值解，相平面图如图所示。由图可知，方程表示成分段光滑非线性数学方由图可知，方程表示成分段光滑非线性数学方程的形式为程的形式为 010311211110101 II31121111010II1 ,)()()()(,)()()()(),()()(SSXcXbX

47、aSSVXXcXbXaSSVXXXX 其中其中 41 II11131 II011121 II20111126)3(2VcNcVScNbVScbNaAAA)cos()(aX)sin()(aX 由于非线性方程的复杂性，不可能直接得到方由于非线性方程的复杂性，不可能直接得到方程的解析解，因此可用平均法来求近似解析解。程的解析解，因此可用平均法来求近似解析解。因为非线性项是小量，故引进小参数因为非线性项是小量，故引进小参数。设方程。设方程的解为的解为 代入运动方程可得标准方程代入运动方程可得标准方程 01043113211211010II143113211211010II1 )(sin)(sin)(s

48、in)()(sin)(sin)(sin)()sin()cos(ddSSacabaaSSVXacabaaSSVXaa 010321121111010II1321121111010II12 )cos()(sin)(sin)sin()()cos()(sin)(sin)sin()()(cosddSSacabaSSVXacabaSSVX ),(),(),(),(221221 yRyRyGyGya )()(dd)()(dd221221yZyZyYyYy 采用采用KB变换变换注意到运动学关系注意到运动学关系 2II120112II101)(2tansinVScbNyVSA 经计算得：经计算得：0103111

49、10102131112211122112131111213111122121 8321 )4sin4(sin321)3cos3(sin121)cos(sin43)2sin2sin)(y2121(21)23)(8321()cos(sin2121SSycyaSS ycybybcyaycyayY 0 )cos(sin41)sin(cos31)cos(sin21)2sin2(sin41)21(2121010010241421113231122121112121 SSSSycybaZ 为求第一次近似的定常解。令为求第一次近似的定常解。令 0)(dd0)(dd11yZyYy 则可得一次稳定的近似定常解。则

50、可得一次稳定的近似定常解。利用定常解的稳定性条件，计算了近似解析利用定常解的稳定性条件，计算了近似解析解的非零解，近似解析解与数值解进行比较，吻合解的非零解，近似解析解与数值解进行比较，吻合良好。通过近似解析解。还得到角频率变化量与搓良好。通过近似解析解。还得到角频率变化量与搓动速度的关系曲线。动速度的关系曲线。振幅振幅y、角频率变化量、角频率变化量Z1与搓动速度的关系曲线与搓动速度的关系曲线 在解决此类分段非线性系统的求近似解时，可在解决此类分段非线性系统的求近似解时，可先通过相平面分析建立分段光滑非线性数学方程，再先通过相平面分析建立分段光滑非线性数学方程，再利用平均法求解。这是得到此类系