结构力学第14章.ppt

第14章 结构动力计算续论 多自由度体系： 多层房屋的侧向振动 不等高排架的振动 块形基础的平面振动 梁上有几个集中质量的振动 求解方法： 刚度法——建立力的平衡方程（位移法） 柔度法——建立位移协调方程（力法） 两个自由度体系——推广到n个自由度体系。 特性（与单自由度区别）： 固有频率：2个——n个； 主振型：2个——n个； 耦合：各自由度的运动相互影响； 不同坐标写方程式（刚度、柔度法） 矩阵形式及运算 刚度法与柔度法实质相同，形式不同。 根据结构的形式适当采用刚度法或柔度法求系数 2、按柔度法求解 ——频率方程 解为n个正实根λi,即1/ωi2（i=1,2, … n）; 得到n个自振频率：ω1，ω2，…ωn ， ——按从小到大顺序排列， 称为第一、第二…第n频率 ——总称为结构自振频谱 将n个自振频率中的任意一个ωk代入特解： n个自由度结构——n个自振频率： ——相应有n个主振动和主振型——特解； 一般解——主振动的线性组合： 一般情况，各质点的振动 ——n个主振动分量叠加而成； ——各主振动的振幅Ak{Y}(k)和初相角αk 取决于初始条件： n个质点的n个初始位移和n个初始速度 ——确定n个Ak和n个αk。 自振频率和振型——结构固有动力特性； ——主要任务 振幅和初相角——由初始条件确定 ? ? 计算步骤 ? 确定积分步长 ??t ? 确定初始位移、初始速度和初始加速度 ? 形成质量矩阵、阻尼矩阵和刚度矩阵； ? 计算拟静力刚度矩阵 ? 计算t+??t时刻拟静力荷载 ? 计算t+??t时刻位移 ? 计算t+??t时刻速度和加速度； ? 重复以上计算直至整个历程结束。 ? ? Wilson 方法是一种无条件稳定的算法。 ? Wilson 方法会使系统增加“算法阻尼”，从而使振幅减小，周期加长。 ? 从泰勒级数展开角度分析，Wilson 方法的截断误差只有?t2量级。 计算说明 ? 经验表明，步长过大实际上还会出现发散现象，称超越现象。 3 Newmark-? 法 ? 美国伊利诺伊斯大学的纽马克教授对线性加速度法进行了修正。 1/基本假定：在[t, ?t]时间间隔内速度和位移变化 引入阻尼假设 第j 阶振型的广义阻尼系数 第j 阶振型的广义阻尼比 有关阻尼矩阵的补充说明 运动方程中引入了阻尼矩阵C，其元素cij的物理意义为：第j 个位移方向有单位速度(其他质量位移方向的速度为零)所引起的第i 个位移方向的阻尼力，称为阻尼影响系数(damplng influence coefficient)。然而，在处理实际问题时，要直接确定阻尼影响系数是十分困难的。 在振型叠加法对运动方程进行处理时，为了使方程解耦，假设体系的阻尼矩阵对振型满足正交性条件,并引入了广义阻尼系数。虽然如此，但如何确定阻尼影响系数仍是一个尚未解决的问题。 下面介绍C 的一种构成方法。 引入Rayleigh(瑞雷)比例阻尼如下 也即认为阻尼和系统质量、刚度成正比，? 与? 可用振型正交性由阻尼比?i，?j 和频率?i，?j 确定。 §14-6 近似法求自振频率 1、能量法求第一频率——Rayleigh法 根据能量守恒定律，当不考虑阻尼自由振动时，振动体系在任何时刻的动能T和应变能U 之和应等于常数。 根据简谐振动的特点可知：在体系通过静力平衡位置的瞬间，速度最大（动能具有最大值），动位移为零（应变能为零）；当体系达到最大振幅的瞬间（变形能最大），速度为零（动能为零）。对这两个特定时刻，根据能量守恒定律得： Umax=Tmax ω 求Umax ，Tmax 求频率 如梁上还有集中质量mi， 位移幅值 . Yi为集中质量mi处的位移幅值。 假设位移幅值函数Y(x)必须注意以下几点： 1、必须满足运动边界条件： （铰支端：Y=0；固定端：Y=0，Y′=0） 尽量满足弯矩边界条件，以减小误差。剪力边界条件可不计。 2、所设位移幅值函数应与实际振型形状大致接近；如正好与第 n 主振型相似，则可求的ωn的准确解。但主振型通常是未知的，只能假定一近似的振型曲线，得到频率的近似值。由于假定高频率的振型困难，计算高频率误差较大。故 Rayleigh法主要用于求ω1的近似解。 3、相应于第一频率所设的振型曲线，应当是结构比较容易出现的变形形式。曲率小，拐点少。 4、通常可取结构在某个静荷载q（x）（如自重）作用下的弹性曲线作为Y（x）的近似表达式。此时应变能可用相应荷载q（x）所作的功来代替，即 2）假设均布荷载q作用下的挠度曲线作为Y(x) 例1 试求等截面简支梁的第一频率。