固有特性近似计算.pptx

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,#,第四章,多自由度系统固有特性的近似计算,多自由度系统固有频率与主振型的计算是振动分析中的最重要和最基本的计算。,当自由度数较大时，求解计算的工作量较大；,工程上可能只需求出较低阶的若干固有频率与主振型。,固有特性的近似计算方法：,瑞利能量法；,邓克利法；,李兹法；,矩阵迭代法；,子空间迭代法；,(,略,),传递矩阵法。,4.1,瑞利（,Rayleigh,）能量法,先假设系统的振型，再用能量法计算其固有频率。,设有一,n,自由度系统，其质量矩阵为,M,刚度矩阵为,K,，它的动能与势能为：,系统作某阶主振动时，最大动能与最大势能：,根据机械能守恒定律,T,max,=,U,max,，可求：,R,I,(A),称为第一瑞利商。,当,A,分别取为系统的各阶主振型,A,(i),时，即可求出各阶固有频率：,注意：由于,A,是假设的振型，因此求出的各固有频率只能是估计值；由于高阶振型很难做出合理的假设，故一般只能估算最低阶固有频率。,结论：若假设的,A,接近于第一阶主振型,A,(1),，由瑞利能量法计算出来的,w,2,确实接近于第一阶固有频率的平方值,w,1,2,，而且比实际值大（此即所谓的上限估值）。,4.1,瑞利（,Rayleigh,）能量法,结论,：若假设的,A,接近于第一阶主振型,A,(1),，由瑞利能量法计算出来的,w,2,确实接近于第一阶固有频率的平方值,w,1,2,，而且比实际值大（此即所谓的上限估值）。,证明如下：,对于,n,自由度系统存在,n,个特征值,w,i,2,，对应有,n,个特征矢量,A,N,(i),(,设已经正则化），它们是相互线性独立的，可以构成,n,维线性空间的一个完备基。任何一个假设振型,A,都可以表示为,n,个正则振型,A,N,(i),的线性组合：,其中,c,1,，,c,2,，,，,c,n,为比例因子，表示相应主振型在假设振型中所占比例的大小。,若假设的,A,接近于第一阶主振型,A,(1),，,则有：,4.1,瑞利（,Rayleigh,）能量法,4.1,瑞利（,Rayleigh,）能量法,假设,上限估值,结论,：若假设的,A,接近于第一阶主振型,A,(1),，由瑞利能量法计算出来的,w,2,确实接近于第一阶固有频率的平方值,w,1,2,，而且比实际值大（此即所谓的上限估值）。,4.1,瑞利（,Rayleigh,）能量法,结论,：若假设的,A,接近于第一阶主振型,A,(1),，由瑞利能量法计算出来的,w,2,确实接近于第一阶固有频率的平方值,w,1,2,，而且比实际值大（此即所谓的上限估值）。,将,得：,代入,瑞利能量法也可用于由柔度矩阵,d,建立系统运动方程的情况：,两次导数代入上式，得：,又因：,4.1,瑞利（,Rayleigh,）能量法,同样可以证明对于同一假设振型,A,总存在有：,即用第二瑞利商算出的固有频率比用第一瑞利商算出的更接近真值。,例,4.1,在如图（例,3.5,）所示的三自由度系统，,求此系统的固有频率和主振型。,解：前面已经算出系统的质量矩阵和刚度矩阵：,柔度矩阵为：,4.1,瑞利（,Rayleigh,）能量法,不妨先粗糙地取假设振型,A,=1 1 1,T,由此假设振型可求得：,若取,A,=1 2 3,T,可求得：,若取,A,=3 5 6,T,可求得：,4.1,瑞利（,Rayleigh,）能量法,第一个假设振型,A,=1 1 1,T,：相当在质量,m,1,上沿坐标方向作用一单位力时三个质量的静态位移的相对值。,第二个假设振型,A,=1 2 3,T,：相当在质量,m,3,上沿坐标方向作用一单位力时三个质量的静态位移的相对值。,第三个假设振型,A,=3 5 6,T,：相当在各质量上沿坐标方向同时作用一单位力时三个质量的静态位移的相对值。,实际振型,A,(1),=0.455 0.801 1.000,T,实际,第三种假设振型与实际最接近，其第二瑞利商与第一阶固有频率的平方最接近。,4.2,邓克利（,Dunkerley,）法,瑞利能量法给出了系统最低阶固有频率的上限估值。,邓克利法将给出系统最低阶固有频率的下限估值。,n,自由度系统的自由振动位移方程：,设解为：,将其代入振动位移方程，并以,w,2,除全式，得振型方程：其特征方程为：,当系统的质量矩阵,M,为对角阵时：,特征方程变为：,即：,展开得：,4.2,邓克利（,Dunkerley,）法,特征方程：,其刚度为：,考虑到系统的固有频率：,则近似地只保留第一项,1/,w,1,2,。,等式左边：,等式右边：,是第,i,个质量处作用单位力时，系统在该处的柔度系数。,设想系统只有一个质量,m,i,存在，则成为单自由度系统，,设这种假想的系统的固有频率为,W,i,，则有：,故，有：,4.2,邓克利（,Dunkerley,）法,特征方程展开后得,将各阶固有频率代入特征方程并相加得,系统最低阶固有频率的平方值,w,1,2,的倒数，近似地等于各质量,m,i,（,i=12,，,，,n,）单独存在时所得各固有频率平方值,W,i,2,的倒数的和。,由于上式的左边舍去了一些正数项，由此求出,1/,w,1,2,的比实际值偏大，即,w,1,2,比实际值偏小，因此，用邓克利法将给出系统最低阶固有频率的下限估值。,上式右边所有各项的和实际上是特征方程,第二个矩阵的迹，故邓克利法又叫迹法。,邓克利法也适用质量矩阵为非对角阵的情况，只是对应矩阵迹的计算复杂些。,4.2,邓克利（,Dunkerley,）法,这里：,n,阶方阵主对角线上元素之和称为矩阵的迹。,例,4.2,用邓克利法估算例,3.5,中系统的第一阶固有频率。,解：在例,4.1,中已经求出：,则：,其迹为：,故：实际,矩阵迹的表示符号,例,4.3,已知一均匀悬臂梁的第一阶固有频率,式中,EJ,为梁的抗弯刚度，,M,为梁的质量，,l,为梁长。,若在梁的自由端安装一质量为,m,的激振器，试用邓克利法估算这时系统的一阶固有频率。并分别计算,m,为,M,/20,、,M,/10,、,M,/5,和,M,/2,时的一阶固有频率，说明激振器质量对均匀梁的一阶固有频率的影响。,解：,只考虑梁的质量：,悬臂梁端点的柔度为：,只考虑激振器的质量：,第一阶固有频率：,将,m,/,M,的各值代入上式，可得下表：,表中的误差是相对于,W,1,的。,加激振器后，梁的固有频率有明显的下降，只有当激振器的质量小于梁的重量的,1/20时，才可以忽略其质量对梁的固有频率的影响。另外，考虑到邓克利法给出的是下限估值，实际误差比表中给出的要小些。,激振器质量对均匀梁的一阶固有频率的影响：,例,4.4,设如图所示的二自由度系统，其中,k,1,=ck,，,c,为常数。用迹法估算其基频（一阶固有频率），并将结果与准确值 作比较。,解：,分别只考虑一个质量的假想系统如图所示：,串联,并联,式中，假设模态 为,n,阶列阵，可预先选定；,李兹法不是直接给出假设振型，而是把它表示为,s,个（,1s=n,）独立的假设模态（即假设振型）的线性组合。,4.3,李兹（,Ritz,）法,李兹法是瑞利法的改进。,用李兹法不仅可以算出系统的基频，还可算出系统的前几阶频率和模态。,李兹法将对近似振型做出合理的假设，从而使算出的基频值进一步下降。,瑞利法算出的基频的精度取决于假设的振型对第一阶主振型的近似程度，而且得到的基频总是精确值的上限。,S,根据需要而定，一般小于,n,；,C,j,为待定系数。,将,假设振型,代入第一瑞利商,得：,由于由瑞利商得出的是系统的一阶固有频率的上限估值，因此，待定系数,c,j,的选择应使上式给出的固有频率为最小。,则：,此式对任意待定系数,c,j,的偏导数都应等于零。,4.3,李兹（,Ritz,）法,得：,由于：,且：,4.3,李兹（,Ritz,）法,阶：,1*s,阶：,s*n,阶：,1*n,得出：,整理成矩阵方程得,：,由于,4.3,李兹（,Ritz,）法,4.3,李兹（,Ritz,）法,由,确定的前,s,阶固有频率以及,的逼真程度与开始选定的,s,个假设模态,有关。,4.3,李兹（,Ritz,）法,4.3,李兹（,Ritz,）法,4.3,李兹（,Ritz,）法,若将,代入第二瑞利商,经过上述同样的处理过程，可将原问题转化为下述特征值问题：,4.3,李兹（,Ritz,）法,将各段质量等分两半，分别置于其两端部，各质量之间由刚度为,k,的弹簧相连，,k,由每段杆的拉压刚度确定为：,例,4.5,图为一等直杆，杆长,l,，截面积为,A,，密度为,r,，试用聚缩质量的方法将其离散为有限自由度系统，并用李兹法求杆作纵向振动的第一阶固有频率及主振型的近似解。,解：将杆分成四段，每段质量：,通过聚缩质量将等直杆简化成四自由度系统。,固定端处的质量不参与振动。,此离散系统的质量矩阵、刚度矩阵和柔度矩阵分别为：,因只求第一阶固有频率和主振型，故选取两个假设模态：,广义质量矩阵：,广义刚度矩阵：,特征方程：,特征方程：,将,代入,特征矢量,同理，根据,确定特征方程：,特征矢量,主振型,等直杆的精确解：,4.4,矩阵迭代法,n,自由度无阻尼系统的自由振动,作用力方程,:,设方程的解为：,代入,作用力方程,，消去,得系统的主振型方程：,n,自由度无阻尼系统的自由振动,位移方程,:,设方程的解为：,代入振动,位移,方程，消去,得系统的主振型方程：,4.4,矩阵迭代法,系统的主振型方程：,矩阵迭代法：从假设主振型出发，对上两式进行矩阵迭代运算，从而算出系统的固有频率和主振型。,设动力矩阵,分析,迭代计算：,4.4,矩阵迭代法,第一阶固有频率和主振型的矩阵迭代运算方法：,（,1,）任意选取一个经过归一化的假设振型,用动力矩阵,D,前乘它，并对矩阵相乘结果再进行归一化,得一新振型，即,为新振型矢量归一化后的系数。,（,2,）如果,就再从,A,1,开始，重复上述步骤，得：,为新振型矢量归一化后的系数。,4.4,矩阵迭代法,主振型方程：,（,3,）如果,继续重复上述步骤。,经过,k,次矩阵乘法运算后，得到：,在规定的有效位数内，当发现,就停止运算。,证明如下：,对于,n,自由度系统存在,n,个特征值,w,i,2,，对应有,n,个特征矢量,A,N,(i),(,设已经正则化），它们是相互线性独立的，可以构成,n,维线性空间的一个完备基。任何一个假设振型,A,都可以表示为,n,个正则振型,A,N,(i),的线性组合：,4.4,矩阵迭代法,设初始假设振型,A,0,表示为系统主振型,A,(i),的线性组合：,假设振型,A,0,经过第一次迭代后：,4.4,矩阵迭代法,经过第二次迭代后：,4.4,矩阵迭代法,4.4,矩阵迭代法,例,4.6,用矩阵迭代法求例,3.5,中系统的第一阶固有频率和主振型。,解：,系统的质量矩阵和柔度矩阵分别为：,系统的动力矩阵为：,系统的动力矩阵为：,取,A,0,=1 1 1,T,进行第一次迭代：,进行第二次迭代：,继续重复迭代运算，得：,这时，有：,第一阶主振型为：,在例,3.5,中通过解频率方程求出的三个固有频率和主振型为：,例,4.6,用矩阵迭代法求系统的第一阶固有频率和主振型。,例,4.1,中通过瑞利能量法，,在假设振型,A=1 1 1,T,，求得：,系统的主振型方程：,4.4,矩阵迭代法,在矩阵迭代法的证明中，,设初始假设振型,A,0,表示为系统主振型,A,(i),的线性组合：,4.4,矩阵迭代法,假设振型,用,前乘上式，根据主振型的正交性，可得：,第,i,阶主质量,取：,即可在假设振型,A,0,中清除前,s,阶主振型分量，使迭代结果收敛于第,s+1,阶固有频率和主振型。,4.4,矩阵迭代法,S,阶清型矩阵（,n*n,阶）,4.4,矩阵迭代法,动力矩阵：,清型变换后的动力矩阵：,4.4,矩阵迭代法,需求第二阶固有频率和主振型时，每次迭代前须乘以动力矩阵：,需求第三阶固有频率和主振型时，每次迭代前须乘以动力矩阵：,动力矩阵的递推公式为：,此种方法适用于正定系统。,4.4,矩阵迭代法,半正定系统：,因不存在柔度矩阵，故不能用上述方法求其固有频率和主振型。,半正定系统的主振型方程：,变成：,再改为：,4.4,矩阵迭代法,系统的主振型方程,例,4.7,用矩阵迭代法求例,3.5,中系统的第二阶和第三阶固有频率及主振型。（继续例,4.6,）,解：在例,4.6,中已求得系统的第一阶固有频率的平方值和主振型为：,系统的第一阶主质量为：,系统的第二阶固有频率和主振型：,动力矩阵：,在例,4.6,已求出,选取第二阶初始振型为：,第一次迭代，得：,第二次迭代，得：,迭代公式,如此继续下去，得：,有,系统的第二阶主振型为：,系统的第二阶固有频率的平方为：,系统的第三阶固有频率和主振型：,第二阶主质量为：,因,选取第三阶初始振型为：,各次迭代运算的结果为：,系统的第三阶主振型为：,系统的第三阶固有频率的平方为：,例,3.5,的结果：,与例,3.5,的结果相比：频率值相同，振型稍有不同，是计算的舍入误差造成的。另外，由于第三假设振型较接近实际振型，迭代次数少。,第四章,多自由度系统固有特性的近似计算,固有特性的近似计算方法：,瑞利能量法；,邓克利法；,李兹法；,矩阵迭代法；,子空间迭代法；（略）,传递矩阵法。,4.6,传递矩阵法,4.6.1,轴的扭转振动,如图为一链式轴盘扭转振动系统，各轴段的质量不计，其扭转刚度分别为,k,1,、,k,2,、,、,k,n,，各圆盘的弹性变形不计，其转动惯量分别为,I,1,、,I,2,、,、,I,n,。,分析第,i,个圆盘：,链式轴盘扭转振动系统，用传递矩阵法计算其固有频率和主振型。,约定：（按右手螺旋规则），,各圆盘的转角矢朝右为正；,作用于圆盘右端面的扭矩矢朝右为正；,作用于圆盘左端面的扭矩矢朝左为正。,将每个端面的转角和扭矩组成一个状态矢量,Z,，反映该处的运动与受力状态。,第,i,个圆盘右端面的状态矢量,第,i,个圆盘左端面的状态矢量,第,i,个圆盘分离体的动力学方程：,设轴的扭转振动为简谐振动：,代入动力学方程，得：,圆盘左右两端面转角相等（忽略圆盘的弹性）：,第,i,个圆盘右端面的状态矢量,第,i,个圆盘左端面的状态矢量,4.6,传递矩阵法,分析第,i,段轴：,4.6,传递矩阵法,第,i,段轴,左右端状态传递关系,第,i,个圆盘,左右端面状态传递关系,C,i,称为第,i,个传递矩阵,4.6,传递矩阵法点矩阵,场矩阵,4.6,传递矩阵法,在上图的链式轴盘扭转自由振动系统中的边界条件：,第,i,个圆盘右端面的状态矢量,4.6,传递矩阵法,边界条件：,自由端,:,固定端：,例,4.9,用传递矩阵法求图示三圆盘扭振系统的固有频率和主振型。,第一个圆盘左右端的点矩阵,第二个传递矩阵,第三个传递矩阵,系统为刚体转动,主链系统：,支链系统：,分支结构的轴盘扭转振动系统，用传递矩阵法计算其固有频率和主振型。,转动惯量为 的三个圆盘所在的,A,轴。,转动惯量为 的圆盘所在的,B,轴。,4.6,传递矩阵法,4.6,传递矩阵法,传递矩阵,C,i,4.6,传递矩阵法,建立了,1,点左右两端面的状态矢量的关系,46,传递矩阵法,4.6,传递矩阵法,4.6.2,梁的弯曲振动,将梁结构简化成带有若干集中质量的弹性梁。,4.6,传递矩阵法,4.6,传递矩阵法,4.6,传递矩阵法,4.6,传递矩阵法,且：,有,：,4.6,传递矩阵法,质量,m,i,梁,l,i,4.6,传递矩阵法,4.6,传递矩阵法,4.6,传递矩阵法,4.6,传递矩阵法,边界条件：,自由端,:,边界条件：,固定端：,第四章,多自由度系统固有特性的近似计算,1,。瑞利能量法,第一瑞利商,R,I,(A),：,第二瑞利商,R,II,(A),：,结论,：若假设的,A,接近于第一阶主振型,A,(1),，由瑞利能量法计算出来的,w,2,确实接近于第一阶固有频率的平方值,w,1,2,，而且比实际值大（此即所谓的上限估值）。,2,。邓克利法,瑞利能量法给出了系统最低阶固有频率的上限估值。,邓克利法将给出系统最低阶固有频率的下限估值。,振动微分方程：,特征方程：,振型方程：,特征方程化为：,系统最低阶固有频率的平方值,w,1,2,的倒数，近似地等于各质量,m,i,（,i=1,，,2,，,，,n,）单独存在时所得各固有频率平方值,W,i,2,的倒数的和。,上式右边所有各项的和实际上是特征方程中的第二个矩阵的迹，故邓克利法又叫迹法。,这里：,n,阶方阵主对角线上元素之和称为矩阵的迹。,式中，假设模态 为,n,阶列阵，可预先选定；,李兹法不是直接给出假设振型，而是把它表示为,s,个（,1s=n,）独立的假设模态（即假设振型）的线性组合。,3,。李兹（,Ritz,）法,李兹法是瑞利法的改进。,用李兹法不仅可以算出系统的基频，还可算出系统的前几阶频率和模态。,瑞利法算出的基频的精度取决于假设的振型对第一阶主振型的近似程度，而且得到的基频总是精确值的上限。,S,根据需要而定，一般小于,n,；,C,j,为待定系数。,将,广义刚度矩阵,广义质量矩阵,代入第一瑞利商,得,s,阶振型方程：,若将,代入第二瑞利商,经过上述同样的处理过程，可将原问题转化为下述特征值问题：,4.3,李兹（,Ritz,）法,系统的主振型方程：,矩阵迭代法：从假设主振型出发，对上两式进行矩阵迭代运算，从而算出系统的固有频率和主振型。,设动力矩阵,4,。矩阵迭代法,n,自由度无阻尼系统的自由振动,作用力方程,:,n,自由度无阻尼系统的自由振动,位移方程,:,第一阶固有频率和主振型的矩阵迭代运算方法：,（,1,）任意选取一个经过归一化的假设振型,用动力矩阵,D,前乘它，并对矩阵相乘结果再进行归一化,得一新振型，即,为新振型矢量归一化后的系数。,（,2,）如果,就再从,A,1,开始，重复上述步骤，得：,为新振型矢量归一化后的系数。,主振型方程：,（,3,）如果,继续重复上述步骤。,经过,k,次矩阵乘法运算后，得到：,在规定的有效位数内，当发现,就停止运算。,在假设振型,A,0,中清除前,s,阶主振型分量，可使迭代结果收敛于第,s+1,阶固有频率和主振型。,S,阶清型矩阵（,n*n,阶）,清型变换后的动力矩阵：,需求第二阶固有频率和主振型时，每次迭代前须乘以动力矩阵：,需求第三阶固有频率和主振型时，每次迭代前须乘以动力矩阵：,动力矩阵的递推公式为：,此种方法适用于正定系统。,4.4,矩阵迭代法,半正定系统：,因不存在柔度矩阵，故不能用上述方法求其固有频率和主振型。,半正定系统的主振型方程：,变成：,再改为：,4.4,矩阵迭代法,系统的主振型方程,6,传递矩阵法,6.1,轴的扭转振动,C,i,称为第,i,个传递矩阵,边界条件：,自由端,:,固定端：,6.2,梁的弯曲振动,