第14章结构动力计算续论ppt课件.ppt

第14章 结构动力计算绪论,14-1 多自由度体系的自由振动,14-2 多自由度体系主振型的正交性和主振型矩阵,14-3 多自由度体系的强迫振动,14-4 无限自由度体系的自由振动,14-5 无限自由度体系的自由振动的常微分方 程求解器解法,14-6 近似法求频率,14-7 矩阵位移法求刚架的自振频率,14-8 用求解器求解自振频率和振型,14-9 小结,14-1 多自由度体系的自由振动,1.刚度法,振动方程为,设振动方程解的形式为,将上式代入振动方程，得,若得到非零解，则,展开形式为,（a）,解行列式，得到n个体系的自振频率,令,由此可求出第 i 振型,（b）,式（b）是一组齐次方程，只能确定主振型的形状，但不能位移地确定它的振幅。,振型的标准化,规定某个元素的值，如第一个元素等于1，或者 最大的一个元素等于1,规定主振型满足下式,例14-1 试求图示刚架的自振频率和振型。设横梁的变形忽略不计，层间刚度系数和质量如图所示。,解（1）求自振频率,刚度矩阵和质量矩阵分别为,频率方程为,展开，得,用试算法求得方程的三个根为,因此，三个自振频率为,进一步求得,（2）求振型,令Y311，解得,令Y321，解得,令Y331，解得,将,代入振型方程，得,刚度法振动方程为,故频率方程为,2 柔度法,展开为,相应的振型方程为,例 142 试用柔度法重做例141。,解（1）求自振频率由各层的刚度系数得到各层柔度系数为,柔度矩阵为,频率方程为,展开，得,解得,因此，三个自振频率为,（2）求主振型,任选体系的两个振型,体系的质量矩阵为,则，第一个正交关系为,14-2 多自由度体系主振型的正交性和主振型矩阵,1 主振型的正交性,另一种证明方法,令振型方程中的i分别等于k、l，得,将（a）式两边分别左乘Y（l）T、（b）式两边分别左乘Y（k）T，得,考虑KT=K，MT=M，将（d）式两边转置，得,式（c）-式（d），得,第一个正交关系,将第一个正交关系代入（c），得,对刚度也正交,对于k=l，定义,第k振型的广义质量,第k振型的广义刚度,以Y（k）T前乘下式,得,即,由此得,由广义刚度和质量求自振频率,主振型正交关系的应用,判断主振型的形状特点,第二振型分为两个区，各居结构的两侧，只有这样才能满足正交条件；第三振型分为三区，交替位于结构的不同侧。这样才能符合与第一、第二主振型都彼此正交的条件。,确定位移展开公式中的系数,任意一个位移向量都可按主振型展开,用Y（j）TM前乘上式两边,由正交性，得,由此求得系数为,例 143 验算例141中所求得的主振型的正交性，求出每个主振型相应的广义质量和广义刚度，并求频率,解 由 例141 得知刚度矩阵和质量矩阵分别为,三个主振型分别为,（1）验证对质量矩阵的正交性,同理,（2）验证对刚度矩阵的正交性,同理,（3）求广义质量,同理,（4）求广义刚度,同理,（5）求频率,主振型向量组成的方阵,转置矩阵为,2 主振型矩阵,故,同理,由振型的正交性可知，非对角线上的元素等于零，主对角线上的元素为各振型的广义质量。所以,振动方程为,简谐荷载,若,14-3 多自由度体系的强迫振动,1 n个自由度体系在简谐荷载下的强迫振动,在平稳阶段，各质点也作简谐振动，即,代入振动方程，整理后，得,令,若D00，则,讨论,故,当荷载频率与其中任意一个自振频率相等时,都可能出现共振现象，因此，对n个自由度体系，存在n个共振区。,振动方程,将位移向量按振型分解,代入振动方程，并前乘YT,令FP=YTFP（t）广义荷载向量,振动方程变为,2 多自由度体系在一般荷载下的强迫振动,由于M*、K*都是对角阵，方程已经解偶，即,同理，令,则,振型分解法,由杜哈梅积分，得,初始条件为,代入初始条件，得,例 14-4 已知结构的频率和振型，试求图示结构在突加荷载FP1作用下的位移和弯矩。,解（1）主振型矩阵,（2）建立坐标变化关系,（3）求广义质量,（4）求广义荷载,（5）求正则坐标,（6）求质点位移,质点1的位移时程曲线,实线：,虚线：,（7）求弯矩,振动过程中质点所受的荷载与惯性力之和为,截面1的弯矩为,截面1弯矩时程曲线,实线：,虚线：,只考虑第一振型,（8）讨论,由于第一和第二主振型分量并不是同时达到最大值，因此不能简单地把两分量的最大值相加。,第二主振型分量的影响比第一主振型分量的影响要小的多。,阶次愈高的振型分量的影响愈小，通常可以计算前23个低阶振型的影响，就可以得到满意的结果。,按无限自由度体系计算可以了解近似计算方法的应用范围和精确程度。,将无限自由度体系简化为有限自由度体系进行计算，是不完整的。,对某种类型的结构，直接按无限自由度体系计算也有方便之处。,14-4 无限自由度体系的自由振动,在无限自由度体系的动力计算中，时间和位置坐标都是独立变量。振动方程是偏微分方程。,等截面梁弯曲时的静力平衡方程为,在自由振动时，唯一的荷载就是惯性力，即,因此，等截面梁弯曲时的自由振动方程为,用分离变量法求解，令,代入振动方程，并整理得,左边是x的函数，右边是t的函数。因此，两边都与x、t无关。,故得两个常微分方程,两个方程的解分别为,则，振动方程的解为,C1C4由边界条件确定,例 14-5 试求等截面简支梁的自振频率和主振型。,右边：,振幅曲线简化为,解：边界条件引入振幅曲线,左边：,得：,令系数行列式=0,得,故,这样就得到了无限多个自振频率和对应的振型曲线,14-5 无限自由度体系自由振动的常微分方程求解器解法,等截面两弯曲时的自由振动偏微分方程为,n=1：表示下段结果；n=2：表示上段结果。,令,代入振动方程，得,边界条件,顶部（x=0）：弯矩=0、剪力=0,中部（x=H2）：水平位移、转角、弯矩、剪力都连续,底部（x=H）：水平位移=0、转角=0,将特征值问题转化为标准的非线性ODE问题,首先，利用区域映射技巧作坐标变换,于是有,这个变化将两段区间影射为标准的单位区间0，1,微分方程变为,边界条件变为,顶部（x=0，=0）：自由,中部（x=H2，=1）：连续,底部（x=H，=0）：固定,微分方程已变成常微分方程组特征值问题,利用平凡的ODE技巧和等价的ODE技巧将其转化为标准的非线性ODE问题.,建议平凡的ODE，即,取振型归一化条件为,分段考虑并利用坐标变换，有,利用等价ODE技巧将该积分转化为标准的ODE问题.,这样，就形成了一个标准的非线性常微分方程组，可直接利用标准的ODE求解器的非线性功能求解。,利用COLSYS求解的计算步骤如下：,设要求解前N个特征值,（1）设初始解,（2）对第k（k=1，2，MN振型求正交化的初始解,其中,（4）用COLSYS求解如下的一个一阶线性ODE问题,然后，求出,（5）回到第（2）步作第k+1步求解。,例 14-6 图示变截面柱，计算数据如下：,1（下）段：2（上）段：,其中s为一比例系数。计算s=1.0,0.5,0.1三种情况。,解：前5个自振频率在下表中给出，相应的振型如图所示。,例14-6的自振频率,计算结果表明：（1）当上下段的质量比和刚度比变小（即s变小）时，基本频率变大；但高阶频率不一定如此。（2）在三种情况中，s=0.1时的振型在顶部位移很大（注意上下部的位移比），通常这种现象称为鞭梢效应；当s更小时，鞭梢效应将更严重。,一个无阻尼的弹性体系自由振动时，在任一时刻的总能量（应变能与动能之和）保持不变。,14-6 近似法求自振频率,1 能量法求第一频率瑞利（Rayleigh）法,理论基础：能量守恒原理,例 具有分布质量的等截面梁，自由振动时，位移可表示为,梁的弯曲应变能为,位移表示式对时间微分，得速度表达式为,最大值为,最大值为,梁的动能为,位移和应变能为零，体系的总能量为Tmax,速度和动能为零，体系的总能量为Vmax,由能量守恒原理，可得,由此得到计算频率的公式,若梁上还有集中质量mi，计算公式为,如果Y（x）是第i振型，则得到的就是第i频率的精确解,取某个静荷载下的位移曲线作为Y（x）。,这时，应变能可用荷载作的功来代替，即,频率计算公式为：,取结构自重的变形曲线作为Y（x）。,例 14-7 试求等截面简支梁的第一频率,解（1）将抛物线作为Y（x）。,（2）将均布荷载作用下的位移曲线作为Y（x）。,（3）将正弦曲线作为Y（x）。,（4）讨论。,正弦曲线是第一主振型的精确解，因此由它求得的是第一频率的精确解。根据均布荷载作用下的挠度曲线求得的结果具有很高的精度。,例 14-8 试求图14-15所示楔形悬臂梁的自振频率。设梁的截面宽度b=1，截面高度为直线变化：,解,单位长度质量,截面惯性矩,设位移形状函数为,代入频率计算公式，得,精确解为,误差为3%,理论基础：哈密顿（W.R.Hamilton)原理,在所有可能的运动状态中,精确解使,驻值,得哈密顿泛函,驻值,Y(x)是满足边界条件的任意可能位移函数,2.能量法求最初几个频率瑞利-里兹（Rayleigh-Ritz）法,瑞利-里兹（Rayleigh-Ritz）法的具体步骤:,(1)将体系的自由度折减为n个自由度,位移函数表示为,:n个可能的位移函数;,a:待定系数。,（2）将位移函数代入哈密顿泛函，得,令,得,应用驻值条件,得,写成矩阵形式,令系数行列式为零，即,可求得最初几个自振频率的近似值。,例 14-9 试求等截面悬臂梁的最初几个频率。,设可能位移为,解,其中,（1）第一次近似,得,驻值条件为,令,得,（2）第二次近似解,得,令,则，第一、二频率的近似值,(误差为0.48%）,(误差为58%）,这里第一频率的精度已大为提高。,例 14-10 试用集中质量法去等截面简支梁的自振频率。,解,3 集中质量法,例 14-11 试求框架的最低频率。,解,读者可自行验证，对称振型的频率大于反对称振型的频率,14-7 矩阵位移法求刚架的自振频率,1 单元的泛函,将刚架分成有限个单元，任一单元的哈密顿泛函为,刚架的泛函,根据刚架泛函为驻值的条件，求的非零解，得到刚架频率,可用单元的结点位移表示,单元的结点位移幅值为,杆件的位移幅值函数可表示为,形状函数列阵,其中,单元的刚度矩阵,单元的质量矩阵,对单元泛函叠加，得,将EP改用刚架的结点位移幅值来表示。,2 刚架的泛函,应用驻值条件，得,频率方程为,3 驻值条件和频率方程,频率方程为,精确解为,误差为1.6%,例14-12 试求梁的自振频率。,解（1）对称振型取半边结构作为一个单元，只有一个待定的结点位移。,将半边结构分为两个单元,待定的结点位移幅值为,驻值条件为,令系数行列式为零，求得三个频率及其误差如下：,（2）反对称振型,取半边结构，分成两个单元，得另外三个频率,例 14-13 试用矩阵位移法从做例14-11,解,总刚度矩阵及总质量矩阵,待定的结点位移幅值为,驻值条件为,对称振动时，得,求得,反对称振动时，得,其中,频率方程为,求得,按从大到小的顺序重新排列,14-8 用求解器求解自振频率与振型,对一般平面结构，可以给出振型和频率；对于无限自由度体系，可以给出全部精确解；解出的振型可以用静态、动态两种方式显示。,14-9 小结,讨论了多自由度体系的振动问题，深化了主振型、主振型的正交性、主振型矩阵的概念；对于一般荷载，介绍了主振型叠加法，将多自由度体系的振动问题转化为单自由度体系的计算问题是这个方法的核心；近似计算方法中，能量法是计算自振频率的一种有效的近似方法。,