[物理]结构动力学.ppt

《结构动力学基础》 《结构动力学基础》目录 绪论 体系的运动方程建立 单自由度体系的振动 多自由度体系的振动 频率和振型的实用计算方法 随机振动初步 结构地震反应分析 结构的振动控制 一、绪论 1.1 阪神地震录像 1.2 动力荷载及其分类 1.3 结构动力学的研究内容和任务 1.4 结构动力分析中体系的自由度 1.5 结构的动力特性 1.6 建立结构运动方程的一般方法 四、 多自由度体系的振动 多自由度无阻尼自由振动 振型的正交性 多自由度的受迫振动 杆系结构有限元动力分析 多自由度时程分析方法 结论与讨论 五、实用计算方法 能量法求基本频率 迭代法求频率、振型 结论与讨论 4.4 杆系结构有限元动力分析 4.4.2 几点说明 1）以单元上无荷载作用，仅产生单位位移的形函数作为单元位移场，这是常用的一种近似处理。 2）结构一致质量矩阵和结构刚度矩阵非零元素分布一样。 3）Clough教授曾经指出，对于框架结构，将杆件一半质量集中在杆端，用集中质量法计算不仅在处理后可减少未知数个数（自由度），而且往往精度更好。 4）当采用集中质量法时，[M]中相应转动自由度的对角线元素（转动惯量）为零，假设位移编码将转动自由度集中在最后编，则无阻尼运动方程分块形式为 [M1][ü]+[K11][u]+[K12][?]=[R1] [K21][u]+[K22][?]=[R2] 由此消去[?]，可得只有线位移自由度的方程。 4.4 杆系结构有限元动力分析 4.4.2 几点说明 5）如果分析时用集中质量法且不考虑轴向变形，则集装后最终质量矩阵是每层质量对角排列的形式。这是目前杆系模型的常用计算方案。 6）对于上述杆系模型的计算程序，质量矩阵很简单。但是集装形成刚度矩阵时，要做4）中所述的“静力缩聚”。当[R2]=[0]时，[K1]=[K11]-[K12][K22]-1[K21]， 运动方程为 [M1][ü]+[K1][u]=[R1] (43) 自由度数等于框架的层数。 7）本节基本原理是对杆系结构进行说明的，象计算结构力学力里一样，思路、方法也可用于其他位移有限元动力分析。 8）程序Vibra可用来计算杆系结构的自振特性等等，请大家使用。 4.5 多自由度时程分析方法 4.5.1 多自由度的线加速度法 在3.3节介绍了单自由度线加速度法，从运动方程的相似性 mü+cú+ku=P(t) [M]{ü}+[C]{ú}+[K]{u}={P(t)} 显然在[0,?t]时间间隔内假设加速度线性变化,则将3.3节m,c,k,P换成[M]、[C]、[K]、{P(t)}，即可得到多自由度线加速度法的等效刚度和等效荷载。 数值积分能做线性、非线性时程分析，对非正交阻尼矩阵也能求解。重要、高层结构要用时程分析。 4.5.2 多自由度的Wilson-?法 线加速度法要求?t小于系统最短周期的1/10，当自由度很多时频率将很高周期很短，这一要求使计算很费时间。而且进一步数学分析表明它是条件稳定的。 4.5 多自由度时程分析方法 Wilson提出，假设??[0,??t]加速度线性变化,仿线加速度法进行推导,可得 [K]*= a0[M]+a1[C] +[K] (44) {P(t+??t)}*={P(t)}+? ({P(t+?t)}-{P(t)})+ +[M](a0{u(t)}+a2{ú(t)}+2{ü(t)})+ +[C](a1{u(t)}+2{ú(t)}+a3{ü(t)}) (45) [K]*{u(t+??t)}={P(t+??t)}* (46) 由式(46)可解出{u(t+??t)}，进一步可以求的t+?t时刻的状态向量。 4.5.3 Wilson-?法的步骤 1）形成系统[M]、[C]、[K]； 2）确定初始状态向量{u(0)}、{ú(0)}、{ü(0)}； 3）确定? (一般为1.4)和?t;按以下公式计算常数 4.5 多自由度时程分析方法 a0=6/(??t)2; a1=3/(??t); a2=2a1; a3=??t/2; a4=a0/?; a5=-a2/? ; a6=1-3/? ; a7=?t/2; a8=?t2/6 (47) 4）按式(44)计算等效刚度； 5）对等效刚度进行LDLT 分解，获得D和L； 6）按式(45)计算等效荷载； 7）用线性方程组的LDLT法解{u(t+??t)}