8 结构的动力计算
8.1一般概念 一、结构的动力荷载及分类 动力荷载,是指荷载的大小、方向、位置随时间迅速变化的荷载;它使结构质量产生不容忽视的加速度,使结构发生明显的振动,即在平衡位置附近往返运动。 静力荷载,是指荷载的大小、方向、位置不随时间变化的荷载; 同时考虑其对结构的影响来看,如果荷载变化极其缓慢,使结构质量产生的加速度可以忽略不计时,仍属于静力荷载 动力荷载分类:周期荷载、冲击荷载、突加荷载、随机荷载
P(t ) P t t P(t ) P tr t (1)周期荷载:随时间周期性变化的荷载 (2)冲击荷载:作用于结构上的荷载值在很短的时间内急剧增大或减小的荷载 P(t ) t P tr (3)突加荷载:在瞬间内将全部重量加于结构或移去的荷载
二、动力计算的内容和研究方法 P(t ) t P(t ) t (4)随机荷载:不能用确定的函数表示,非确定性的荷载 首先要确定动力计算简图,明确动力荷载的性质和规律,然后进行分析。无论是确定结构的动力特性,或是计算动力反应,都是从研究结构质量的运动规律入手,把质点的位移作为基本未知量,建立体系的运动方程,进行分析。
动静法是根据达朗贝尔(d’Alembert)原理,设想将惯性力I(t)加于振动体系的质点上,则任一瞬时体系中的实有各力与惯性力处于平衡状态 动力特性,是指结构的固有的振动频率,基本振动形式(主振型)和阻尼特性等。这些是结构自身的固有特性,与外部作用因素无关。 动力反应,是指动力荷载作用下,结构产生的内力、位移、速度、加速度等。不仅与荷载的大小、方向、作用位置及其变化规律有关,即是时间的函数;还与结构的动力特性有关。 与静力计算的对比:两者都是建立平衡方程,但动力计算,利用动静法,建立的是形式上的平衡方程。力系中包含了惯性力,考虑的是瞬间平衡,荷载、内力都是时间的函数。建立的平衡方程是微分方程 动静法是根据达朗贝尔(d’Alembert)原理,设想将惯性力I(t)加于振动体系的质点上,则任一瞬时体系中的实有各力与惯性力处于平衡状态
8.2动力计算简图和动力自由度 单自由度体系 动力计算中要引入惯性力,因此计算简图要考虑质量的分布 一个动力体系运动过程中确定其任一时刻全部质量位置所需的独立几何参数的数目,为该体系的动力自由度 实际结构的质量都是连续分布的,是无限自由度体系,选取动力计算简图是,常将无限自由度体系化为有限自由度体系。 EI L EI L y(x,t) 单自由度体系
EI L 多自由度体系
8.3单自由度体系的自由振动 一、运动微分方程 (一)动力平衡方程法(刚度法) EI L 体系在没有外部动力荷载作用,而由初始位移(y 0)和初始速度(v 0)引起的振动,叫做自由振动 一、运动微分方程 根据动静法,建立质点的运动方程,可采用两种方式 (一)动力平衡方程法(刚度法) 取质点隔离体为研究对象,质点所受各力保持平衡
m y jw m y(t) 建立运动方程时考虑质点所受的力有: (1)重力 W 为静力荷载 (2)弹性恢复力 与位移成正比,方向与位移指向相反。k为刚度系数,其意义是使质点沿位移方向产生的单位位移时所需的在质点上所加的力 (3)阻尼力 与质点的速度成正比,方向与速度相反。c为粘滞阻尼系数 (4)惯性力 其大小为质点质量与质点加速度之积,方向与加速度方向相反 m y jw y(t) m 可写出平衡方程:
(二)位移方程法(柔度法) m f(柔度系数) y(t) 因 ,可得出质点振动的运动微分方程: 因 ,可得出质点振动的运动微分方程: 当动力位移由质点的静力平衡位置算起时,可不考虑质点的重力 (二)位移方程法(柔度法) m y(t) f(柔度系数) 按动静法,体系的动力位移可看为是由于惯性力和阻尼力静力作用所引起的 可得方程: 柔度系数 f 和刚度系数 k 有如下关系:
二、无阻尼自由振动 (一)运动微分方程解 m 则两种方法所得方程可写成统一形式 令 单自由度体系无阻尼自由振动的运动微分方程: 它是二阶常系数线性齐次微分方程,其通解为: 常数C1,C2由初始条件确定 静力平衡位置 m 设 t=0 质点位移方程:
可知,自由振动由两部分组成:一部分是由初始位移 y引起,按余弦规律振动;另一部分是初始速度 v 引起,按正弦规律振动。 令 可得: 振幅 表示合成运动仍为简谐运动,其中A和φ为: 初相位 y t y t y t y -y T T -A A T
(二)自振周期与频率 m 由运动方程可知自由振动是简谐周期运动 周期 频率 1)结构的周期、频率只与结构自身的质量、刚度(柔度)系数有关,与外 因无关,是结构自身的固有的特性,称为固有周期、固有频率; 2)结构的频率与质量的平方根成反比,与结构刚度系数的平方根成正比; 3)结构的固有周期和频率是结构动力性能的重要标志。 例1 计算结构的频率和周期(EI为常数) m h h f
(三)质点的振动规律 例2 计算结构的频率和周期 1 m h EI 质点的位移、速度、加速度和惯性力分别为: 例2 计算结构的频率和周期 1 m EI h (三)质点的振动规律 质点的位移、速度、加速度和惯性力分别为: (1)当体系处在平衡位置时,加速度及惯性力为零,而速度最大 达到振幅位置时,速度为零,而位移、加速度及惯性力同时达到最大 (2)位移、加速度和惯性力同步变化,利用这一性质,可在质点振幅位置建立运动方程,所得运动方程是代数方程而不是微分方程 (3)弹性力指向永与位移方向相反,而惯性力永与位移方向相同
例: 求图示梁频率 EI=∞ 2a a 此梁为一个自由度体系,振动达到幅值时,两质点的振幅为 A1 A2,惯性力幅值为 由平衡方程 即
二、阻尼对自由振动的影响 实验证明,振动中的结构,不仅产生与变形成比例的弹性内力,还产生非弹性的内力。 事实上,由于非弹性力的存在,自由振动会衰减直到停止,非弹性力起着减小振幅的作用,使振动衰减,最后停止振动。 有阻尼的自由振动运动微分方程: 当 ξ ≥1 时,体系的运动为非振动状态; 当 ξ < 1时,称为低阻尼,体系呈振动运动。 阻尼比 当 ξ < 1,及初始条件 方程可解为: 或 其中
yk yk+1 y t 第k个振幅为 经过一周期相邻两振幅比值 可通过试验的方法测算阻尼比 阻尼比 ξ 越大,振幅衰减越快;振幅按等比级数递减 t 可通过试验的方法测算阻尼比 由自由振动试验曲线量测出任何相邻的两个振幅,算出对数衰减率δ,则可求得阻尼比 可用相隔 i 个周期的两个振幅计算δ,可提高ξ 精确度
8.4单自由度体系的受迫振动 一、运动微分方程 m m m 体系在振动过程中有动力荷载P(t)或支座运动等外部干扰作用时,其振动称为受迫(或强迫)振动。 m m y(t) 由质点的平衡可得: 单自由度体系强迫振动的振动微分方程 或: 若体系的动力荷载不在质点上作用 m y(t)
二、简谐荷载下的无阻尼受迫振动 m 由位移方程可得: y(t) 即 1 f 11 或 1 f 1P 上式中: 设单自由度体系在质点上作用简谐荷载为: 不考虑阻尼,振动微分方程: 扰力幅值 荷载频率
(一)质点的位移方程 齐次解: 特解: 将特解及其二阶导数代入振动微分方程中可确定: 质点位移方程为齐次解和特解之和: 设零初始条件,即 频比,扰力频率与自振频率比 质点位移方程为齐次解和特解之和: 设零初始条件,即 最后 质点的位移方程 由上式可以看出,振动是由两部分合成的;式右第一项是按荷载频率θ的振动。第二项是按自振频率ω的振动。后一部分是由荷载作用引起的称为伴生自由振动。实际上由于存在阻尼,伴生自由振动在短时间内即行消失,最后剩下的仅按荷载频率变化的振动,称为纯强迫振动。在振动开始两种振动共存阶段,称作过渡阶段,以后的纯强迫振动称为平稳阶段或稳态强迫震动。
(二)稳态强迫振动的动力反应 μ β 3 2 1 质点的位移方程: 振幅 当β →0时, μ →1,荷载变化得很慢,可当作静荷载处理 振动频率和荷载频率相同,二者完全同步 质点的位移方程: 振幅 扰力幅值产生的静位移 动力系数 动力系数:最大动力位移与相应静力位移的比值,是衡量动力反应大小的重要指标 1 2 3 β μ 当β →0时, μ →1,荷载变化得很慢,可当作静荷载处理 当0< β <1时, μ >1,并且随β的增大而增大 当β =1时, μ = ∞。即当荷载频率接近于自振频率时,振幅会无限增大。称为“共振” 当β >1时, μ的绝对值随β的增大而减小,且为负值,质点的位移和扰力的指向相反
(1)单自由度体系 (2)受简谐荷载作用 (3)荷载位于质点上 体系满足前面的条件时,不仅荷载、质点位移和惯性力同步变化,而且各种的动力反应的动力系数 μ相同,故可利用 μ求体系的动力反应。 当求结构的最大动力反应时,可用乘以动力系数的扰力幅值P代替扰力幅值P与惯性力幅值 作用,用静力方法计算。 例1:已知m=300kg,EI=90×105N.m2 ,k=48EI/l3 ,P=20kN,θ=80s-1 ,求梁中点的位移幅值及最大动力弯矩。 l/2 EI m 解:1)求ω 2)求β 3)求 ymax, Mmax
例2: 图示简支梁跨中有一集中质量m,支座A处受动力矩Msinθt 作用,求质点的动位移和A的动转角的幅值。 EI l/2 m Msinθt B A 解:体系的动力荷载Msinθt 不是作用在质点,因而不能直接利用 μ求动位移,可由建立体系的振动方程来求解。 1 1)设惯性力和动力荷载分别为单位力和单位力偶作用在体系上,并作出相应的弯矩图 M1 M2 M1 l/4 运用图乘法可得: 1 M2 1 2)质点的动位移是惯性力 I (t) 和动力荷载共同作用下产生的,按叠加原理表示为 式中: 自振频率 将柔度系数代入上式,并整理得 等效荷载幅值
由质点位移方程可得,受迫振动的稳态解为: 质点的动位移幅值为 其中 为动荷载 幅值M所引起的质点静位移yjp,μ为动力系数。 3) 支座A 处的动转角也是由惯性力I (t)和动力荷载共同作用下产生的,按叠加原理表示为: 将y(t)求二阶导数代入上式,可得: 将柔度系数和 代入可得:
支座处动转角幅值为 其中 为动荷载 幅值M所引起的静转角, 为动力系数。 上例表明,动荷载不作用在质量上时,质点的位移的动力系数和支座处动转角的动力系数是不同的,即体系不能用统一的动力系数表示
三、阻尼对受迫振动的影响 振动微分方程: 齐次解: 特解: 全解: 上式右第一大项是按自振频率ωr的振动,由于存在阻尼,这部分很快消失。余下的第二大项是按扰力频率θ的纯强迫振动 频率比: 纯强迫振动位移方程: 如令 与无阻尼的质点位移方程相比,多了相位差φ,质点随仍为简谐振动,但与荷载不同步,位移变化之后于荷载变化。
μ β 1 2 3 振幅A可写为: 动力系数: 阻尼使动力系数减小,动力系数μ与频率比β及阻尼比ξ有关: 1)μ曲线随阻尼比增大而趋于平缓,在β =1附近μ值降低比较快。 2)当β =1时,有阻尼情况 μ≠∞ ξ=0 1 2 β 3 μ ξ=0.2 此时动力系数为: ξ=0.3 ξ=0.5 而最大动力系数不是在β =1处,而是在β值略小于1处 由 ξ=1.0 得
四、在任意动力荷载作用下受迫振动 P (一)瞬时冲量的反应 P(t) Pdt P t dt 任意动力荷载作用下体系的动力反应,可用瞬时冲量作用的反应推导。 P (一)瞬时冲量的反应 体系在静止状态突然作用荷载P,停留时间dt。则瞬时冲量Pdt引起的振动可视为由初始条件引起的自由振动。 由动量定理可得: P(t) t P Pdt t=dt时的速度: dt时间内的平均速度: dt t=dt的位移:
P(t) P(τ)dτ P(τ) t τ dτ 当 t >dt后,体系的振动相当于以d v 和d y为初始条件的自由振动,即取: 上式是在t=0时的作用瞬时冲量的反应,如果在 t =τ时作用瞬时冲量P(τ)d τ P(t) t P(τ)dτ 则有t >τ时的位移反应 : P(τ) τ dτ
(二)任意荷载作用的反应 P(t) P(τ )d τ t dτ τ t 杜哈梅(Duhamel)积分 P(τ )d τ 上式为初始处于静止状态的无阻尼单自由度体系受任意动力荷载作用下计算质点位移的一般公式。 dτ τ t 初始位移y0和初始速度v0不为零在任意荷载作用下的位移公式: 考虑阻尼,则杜哈梅积分为
P(t) t P td 例 单自由度体系在零初始条件下,质点上受到图示的短时突加荷载作用,求质点的位移反应。 解: 荷载表达式为 例 单自由度体系在零初始条件下,质点上受到图示的短时突加荷载作用,求质点的位移反应。 P(t) t 解: 荷载表达式为 P td 计算分两个阶段 1)当 时,将荷载表达式代入杜哈梅积分,得 突加荷载位移反应
2)当 时,为由位移 和速度 引起的自由振动,由杜哈梅积分可得 3)最大反应分析。若将上述两个阶段的位移反应表达式中的频率换以周期来表示则分别为 时,最大位移发生在第一阶段。当 时,有最大位移,动力系数为2 时,最大位移发生在第二阶段。当 时,有最大位移,动力系数为
8.5两个自由度体系的自由振动 一、运动微分方程 m2 I2 (一)位移方程法(柔度法) m1 I1 y2(t) y1(t) f 21 在自由振动过程中任意时刻 t,质量m1、m2的位移y1(t)、y2(t)可看做体系在当时惯性力I1、I2作用下的静力位移。 f 21 f 22 1 f 12 可得两个自由度体系自由振动微分方程: f 11 1
k21 k22 k11 k12 (二)动力平衡方程法(刚度法) m1 m2 m2 1 I2 S2 1 m1 I1 S1 y2t) y1(t) 沿两个自由度方向截取包括相应质量的两个隔离体,各隔离体上作用相应的弹性力和惯性力,建立平衡方程。从而的振动微分方程 两个自由度体系自由振动微分方程
二、固有振动及固有频率、主振型的确定 (一)体系的固有振动 体系中各质点按同频率、同相位的简谐振动,称为固有振动,又叫同步振动。其振动频率叫结构的固有频率。 n个自由度体系有n个固有频率。与每个固有频率相应,体系的振动有一定的振动形式,称作结构的固有振型,又叫主振型,或简称振型。 结构的固有频率和主振型是结构的重要动力特性,固有振动或自由振动分析的目的就是确定结构的固有频率和主振型 对于两个自由度体系的固有振动,微分方程式的解为: 将上式代入两个自由度体系的振动微分方程中,可得:
(二)体系的固有频率 或 固有振动基本方程 固有振动基本方程为关于振幅的A1 、A2齐次方程,振幅不会全为零,则其系数行列式必为零,可得: 频率方程 或
频率方程展开,并考虑 f12=f21 ,且令 ,可得: 或 上式可求得两个实根,为: 或 可解得两个固有频率ω1 和ω2 。其中数值较小的为ω1 ,称为第一频率或基本频率,数值较大的为ω2 ,称作第二频率
(三)体系的主振型 由固有振动 可得 体系振动过程中,质点的位移大小虽然不断变化,但两个质量位置的之比始终相等,即振动形式是确定不变的,故可用振幅之比表示体系的主振型。 两个自由度体系有两个固有频率,每个固有频率与一个特定的振动形式即主振型相对应。已求出ω1 ω2 后,可利用振动基本方程来求出各振幅比值,从而确定主振型。 由 由
柔度系数的副系数f12=f21及刚度系数的副系数k12=k21是有正负号的,因此,求出的振幅比及相对振幅值也可有正号或负号。 关于正负号的规则是,在计算开始时,先规定位移y1(t)及y2(t)的正方向,求柔度系数的单位力及求刚度系数的单位位移均按位移yi(t)的正向施加,则求得的柔度系数fij或刚度系数kij指向与yi(t)正向相同者为正,反之为负。 (四)主振型的正交性 m1 m2 主振动是简谐自由振动,位移和惯性力同时达到幅值。 第一主振型可看做相应惯性力幅值ω12m1A 1(1)和ω12m2A 2(1)作用所产生的静力位。 第二主振型可看作相应惯性力幅值ω22m1A 1(2)和ω22m2A 2(2) 作用所产生的静力位移。
m1 m2 m1 m2 第一主振型 第二主振型 对两个状态应用虚功互等定理,有 整理得: 因 ,则存在: 主振型的正交条件
m 例1 求图示结构自振频率和振型。 A l/2 解:体系为静定结构,有两个自由度 1)求柔度系数,由图乘法和弹簧内力虚功计算,得 B EI k=3EI /l3 A B E D 例1 求图示结构自振频率和振型。 解:体系为静定结构,有两个自由度 1)求柔度系数,由图乘法和弹簧内力虚功计算,得 1 2)求自振频率,代入两个自由度体系自由振动的频率方程,得 M1 l /2 1/2 令: 1 则: M2 展开并解得: l /4 1/2
可得相应的频率: 3)求振型并绘振型图 当λ1= 27.083 时, 当λ2= 2.917 时, 1 1 1.553 0.644 第一振型 第二振型
8.6一般多自由度体系的自由振动 mn 一、运动微分方程 (一)柔度法方程 mi m2 m1 yn(t) yi(t) y2(t) yn(t) y1(t) (一)柔度法方程 n个自由度体系自由振动的某一个时刻质点mi产生的位移为yi(t),此位移可看作由各质点加上惯性力 静力作用产生的,可得n个自由度体系自由振动微分方程
写成矩阵形式: 其中 加速度列阵 位移列阵 质量矩阵 当仅考虑质量的移动惯性作用时,质量矩阵为n阶对角方阵,其元素mi为沿 yi方向的质量 柔度矩阵 柔度矩阵为n阶对称方阵,fij=fji。fij的意义是沿体系的位移yi方向加单位力所产生的沿yi方向的位移
(二)刚度法方程 仿照两个自由度体系的自由振动微分方程,可写出n个自由度体系自由振动微分方程 写为矩阵方程: 其中: 刚度矩阵 刚度矩阵为n阶对称方阵,kij=kji。kij的意义是使体系沿yj方向产生单位位移、沿其他位移方向的位移为零,所引起的附加约束i方向的力
二、固有振动特性的确定 (一)基本方程 设n个自由度体系自由振动仍为简谐振动,各质点位移为: 即 振幅列阵 或 其中[E]为n阶单位矩阵
(二)固有频率 (三)主振型 基本方程是振幅的齐次方程,振幅不会全为零,因此有 或 上两式即为频率方程。展开行列式为ω2的n次代数方程,即可得n个固有频率ωi 将个频率由小到大排列,ω1、ω2…… ωi,称为频率谱,其中ω1最小称为第一频率或基本频率。由此可知,n个自由度体系有n个固有频率 (三)主振型 与每一个固有频率对应都有一个主振型。由自由振动基本方程不能唯一确定主振型的各振幅值,只能确定主振型的形状。为了使表示主振型的相对振幅有确定的值,需要对定某种条件,如此确定的主振型称为规准化主振型,一种简单的表示方法是规定振幅列阵种某个元素值为1,如令第一个元素为1。这样的主振型叫做归一化主振型。
(四)主振型的正交性 如令振幅列阵中第一个元素A1(i),则主振型列阵为: 采用规准化主振型,自由振动基本方程可写为 或 对第i主振型和第j主振型,有 用{φ}T(j)前乘(a)式和 {φ}T(i)前乘(b)式
式(c)转置后其值不变,为 [K]和[M]皆是对称方阵,有 上式可改为 式(e)减式(d),整理后得 当ωi≠ωj时,得到 主振型的第一正交条件 将第一正交条件代入式(d)可得到 主振型的第二正交条件
8.7多自由度体系在简谐荷载下的 受迫振动 P1(t) P2(t) Pi(t) Pn(t) m1 m2 mi mn 一、运动微分方程 P1 yi(t) y2(t) yn(t) y1(t) 设n个自由度体系中在各质点上沿质点运动方向作用有简谐荷载 1 任一质点位移可表示为: fii f2i fni f1i 其中 P1 P2 Pi ΔiP Δ2P ΔnP Δ1P ΔiP为各扰力幅值在yi方向产生的静力位移