结构动力学大作业详解.docVIP

结构动力学大作业 班 级 土木卓越1201班 学 号 U座机电话号码3 姓 名 陈祥磊 指导老师 叶 昆 2014.12.30 结构动力学大作业 ——SDOF体系在任意荷载作用下的动力响应 一、结构参数 计算结构为右图所示的 1、 2、 3、结构参数中；。 确定各阶频率和振型 多自由度体系自由振动时的运动方程为 ...... 写成矩阵形式即为 假设此方程的解答为，带入到运动方程中得到振动方程 此方程要有非零解必须满足频率方程，可解得各阶主频率 再根据 可求出结构的主振型。在主振型中，通常将最后一个位移值设定为1，只要在程序中加入下列语句： MDOF.YMatrix :,i MDOF.YMatrix :,i /MDOF.YMatrix MDOF.ND,i 运行程序之后得到如下结果： 1、各阶频率和周期 W1 12.座机电话号码 T1 0.座机电话号码3 W2 37.1座机电话号码 T2 0.1座机电话号码 W3 58.座机电话号码 T3 0.10座机电话号码 W4 75.座机电话号码 T4 0.0座机电话号码 W5 85.座机电话号码 T5 0.0座机电话号码 各阶阵型 Y1 Y2 Y3 Y4 Y5 0.座机电话号码7 -0.座机电话号码6 1.座机电话号码 -1.座机电话号码6 1.座机电话号码7 0.座机电话号码9 -1.0座机电话号码 0.座机电话号码 1.座机电话号码9 -3.座机电话号码 0.座机电话号码8 -0.座机电话号码4 -1.座机电话号码 0.座机电话号码8 3.座机电话号码2 0.座机电话号码7 0.座机电话号码8 -0.座机电话号码 -1.座机电话号码6 -2.座机电话号码 1 1 1 1 1 由主振型可以求得广义质量矩阵和广义刚度矩阵，均为对角矩阵。程序如下所示： MDOF.MGMatrix zeros MDOF.ND,MDOF.ND ; for i 1:MDOF.ND MDOF.MGMatrix i,i transpose MDOF.YMatrix :,i *MDOF.MMatrix*MDOF.YMatrix :,i ; end MDOF.KGMatrix zeros MDOF.ND,MDOF.ND ; for i 1:MDOF.ND MDOF.KGMatrix i,i transpose MDOF.YMatrix :,i *MDOF.KMatrix*MDOF.YMatrix :,i ; end 三、使用能量法计算近似的一阶频率； 使用能量法求解一阶频率就是瑞利法，Rayleigh法主要用于求的近似解，原理是能量守恒定律：；最终求解的表达式为： 因为采用了近似解答，相当于在精确解条件下添加了约束，提高了结构的整体刚度，所以第一频率会大于精确值。在实际计算时，假设的位移幅值函数决定最终计算的精确度，越接近实际的位移模式，Rayleigh法的计算结果越精确。假设的位移幅值函数必须满足以下几点： ①必须满足运动边界条件，即几何和自然边界条件； ②所设位移幅值函数必须与实际振型形状大致相近。 通常可取结构在耨个静位移作用下的弹性曲线作为的近似表达式，此式应变能可用相应荷载所做的功来代替。 计算时采用列表计算和编程计算两种方法互为佐证以证明结果的正确性，计算结果如下： 1、列表计算法 层数 质量m MiXi MiXi2 W 1 1000000 座机电话号码 2000*E6 0.0245 0.0245 24500 600.25 2 1000000 座机电话号码 2000*E6 0.0196 0.0441 44100 1944.81 3 1000000 座机电话号码 2000*E6 0.0147 0.0588 58800 3457.44 4 1000000 1座机电话号码 2000*E6 0.0098 0.0686 68600 4705.96 5 1000000 座机电话号码 2000*E6 0.0049 0.0735 73500 5402.25 汇总 269500 16110.71 12.80 SDOF_Time_History_Analysis之后，得到结构的位移与时间关系图 放大之后可得细节图 矩阵迭代法 矩阵迭代法是采用逐步逼近的方法确定结构的频率和振型。体系做自由振动时各质点的位移幅值为 这是一个迭代关系式，通过假定最初的带入计算，可以求出，继而不断求得，当与足够接近时即可将作为第一振型。其运算过程如下： ①首先假定一个振型，通常假定为代入上式右边，进行求解后即可得到和主振型的第一次近似值； ②然后以第一次近似值代入上式进行计算，得到和主振型的第二次近似值，如此下去，直到最后两次的计