在进行PushOver分析之后就可以进行时程分析了，当然不是一定要先PushOver然后再时程分析，一般是先进行模态分析再进行时程分析，时程分析需要地震波的输入，这里选择了EL Centro波。模型还是弹性的悬臂梁，截面为弹性材料构成，弹性模量取200000，整个分析的单位为N和mm，如下图所示。
EL Centro波可以参考：EL Centro地震波数据
波形如下图所示：
计算分析的命令流如下所示：
|
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 |
#弹性悬臂梁PushOver分析By伊犁的秋天 #基本参数设置 wipe; model basic -ndm 2 -ndf 3; file mkdir EQ; #模型定义 node 1 0 0 node 2 0 1000 #边界条件 fix 1 1 1 1; #定义水平方向质量M=G/g mass 2 200 1.e-9 0. #坐标转换 geomTransf Linear 1; #定义弹性梁柱单元，命令：element elasticBeamColumn eleTag iNode jNode A E Iz transfTag element elasticBeamColumn 1 1 2 10000 200000 8333333 1; #定义记录器 recorder Node -file EQ/node2.out -time -node 2 -dof 1 2 3 disp; recorder Node -file EQ/node1.out -time -node 1 -dof 1 2 3 disp; recorder Drift -file EQ/drift.out -time -iNode 1 -jNode 2 -dof 1 -perpDirn 2 ; #定义竖向荷载 pattern Plain 1 Linear { load 2 0. -2000. 0.; } constraints Plain; numberer Plain; system BandGeneral; test NormDispIncr 1.0e-8 6; algorithm Newton; integrator LoadControl 0.1; analysis Static; analyze 10; loadConst -time 0.0; #定义地震作用 #地面加速度定义（由于地震波文件单位是m/s^2转换成mm/s^2,放大1000倍） set accelSeries "Series -dt 0.02 -filePath ELCENTRO.DAT -factor 1000"; #施加地震作用 pattern UniformExcitation 2 1 -accel $accelSeries; #定义瑞雷阻尼（根据第一振型计算） rayleigh 0. 0. 0. [expr 2*0.02/pow([eigen 1],0.5)]; # 建立分析 constraints Plain; numberer Plain; system BandGeneral; test NormDispIncr 1.0e-8 10; algorithm Newton; integrator Newmark 0.5 0.25 ; analysis Transient; analyze 1560 0.02; |
分析之后，打开EQ文件夹，可以看到分析结果的文件。通过Excel绘制出顶部自由节点的时程曲线：
最后附上计算的tcl文件：点击下载
