1、 1ANSYS 分析 实例集 Whtao1998 汇集 21ANSYS SOLID65 环向布置钢筋的例子 .4 2.混凝土非线性计算实例（1）- MISO 单压 6 3.混凝土非线性计算实例（2）- MISO 约束压 7 4.混凝土非线性计算实例（3）- KINH 滞回 10 5.混凝土非线性计算实例（4） - KINH 压-拉裂 12 6.混凝土非线性计算实例（5） .13 7.混凝土非线性计算实例（6） .15 8.混凝土非线性计算实例（7）- MISO 滞回 17 9.混凝土非线性计算实例（8） .19 10.混凝土非线性计算实例（9）-梁平面应力 .21 11.四层弹簧质点模型的地震
2、分析 .23 12.悬臂梁地震分析 .49 13.用beam 54 单元描述变截面梁的例子 .73 14.变截面梁实例 .74 15.拱桥浇筑过程分析 -单元生死应用实例 75 16.简支梁实体与预应力钢筋分析实例 76 17. 简单的二维焊接分析单元生死实例 .78 18.隧道开挖（三维）的命令流 85 19.岩土接触分析实例 102 20.钢筋混凝土管的动力响应特性分析实例 110 21.隧道模拟开挖命令流（入门） 117 22.螺栓连接的模拟实现问题 120 23.道路的基层、垫层模量与应力之间的关系 130 23.滞回分析 152 24.模拟某楼层浇注 .154 25.在面上施加移动的
3、面力 .156 27.在任意面施加任意方向任意变化的压力 .160 28.预紧分析 .161 29.几何非线性 +塑性 +接触 +蠕变 163 30.埋设在地下的排水管道 .168 32.幕墙企业玻璃简化计算 .173 33.等截面杆单元生死应用实例 .189 34.梁板建模联系 .190 36.简单的例子如何对结构的振动控制分析 .193 37.模态分析结果的输出实例 .195 38.火车过桥动态加载实例（部分） .197 39.悬索结构的找形和计算的例题 .214 40.陶瓷杆撞击铝板的例子 .219 41.求反作用力的 APDL 命令法 222 42.LS-DYNA 实例（部分） 223
4、 43.路面分层填筑对路基的影响 .224 44.一个例子（含地震影响，求振兴与频率） .228 45.接触面上的压力总和 .232 46.施加位置函数荷载 .236 347.非线性分析考虑刚度退化 .237 48.一个圆形水池的静力分析 .238 49.ANSYS 中混凝土模式预应力模拟的算例 239 50.悬臂梁受重力作用发生大变形求其固有频率 .241 51.循环对称结构模态分析 .243 52.三角平台受谐波载荷作用的结构响应 .245 53.三角平台受一地震谱激励的应力分布和支反力 .247 54.三角平台受时程载荷作用的应力分布和变形过程 .249 55.经典层合板理论 .251
5、56定易圆轨迹的例子 .258 57模拟门式刚架施工单元生死 .258 58钢筋混凝土整体式模型例子 .261 59在荷载步之间改变材料属性例子 .263 60含预应力的特征值屈曲计算 .264 61振型叠加计算及工况组合例子 .266 62柱子稳定分析算（预应力，特征值屈曲，初始缺陷） .269 63 module MConcrete !混凝土模板 272 64.混凝土开裂实例 .280 65螺栓网格划分 .281 66自由液面的土石坝平面渗流分析 .282 67导出刚度矩阵 .286 68某混凝土拱坝工程施工期及运行期温度场仿真分析 .287 69移动温度荷载计算 .294 70SHSD用
6、于壳-实体装配实例 An .296 71 ansys 显示隐式回弹分析实例 .300 72工况组合的经典例子 .315 41ANSYS SOLID65环向布置钢筋的例子 ! 一个管道，环向配筋率为1%，纵向配筋率为0.5%，径向配筋率为0.1% ! FINISH /CLEAR /PREP7 !* ! 单元属性 ET,1,SOLID65 !* KEYOPT,1,1,0 KEYOPT,1,5,0 KEYOPT,1,6,0 KEYOPT,1,7,1 !* !实参数1：不同方向配筋 R,1,2,.001, , ,2, .01, RMORE, 90, ,2,.005 ,90 ,90 , !Adds re
7、al constants to a set. !材料属性 !混凝土基本材料属性 MPTEMP, MPTEMP,1,0 MPDATA,EX,1,30E2 MPDATA,PRXY,1,.2 !屈服准则 TB,MISO,1,1,5, TBTEMP,0 TBPT,0.0005,15 TBPT,0.001,21 TBPT,0.0015,24 TBPT,0.002,27 TBPT,0.003,24 !破坏准则 TB,CONC,1,1,9, TBTEMP,0 TBDATA,.5,.9,3,30, !抗拉和抗压是10 倍的关系 TBDATA,1, MPTEMP, 5MPTEMP,1,0 !钢材基本属性 MPD
8、ATA,EX,2,200E3 MPDATA,PRXY,2,.27 !屈服准则 TB,BISO,2,1,2, TBTEMP,0 TBDATA,310,2E3, ! 管道内径和外径 CYL4,0,0,3000, , , ,10000 CYL4,0,0,2000, , , ,10000 VSBV, 1, 2 ! 定义局部柱坐标 CSWPLA,11,1,1,1, ! Defines a local coordinate system at the origin of the working plane. KWPAVE, 11 !move to keypoint11 WPRO,-90.000000, V
9、SBW, 3 WPCSYS,-1,0 !Defines the working plane location based on a coordinate system. KWPAVE, 1 WPRO,-90.000000 VSBW,ALL ESIZE,500,0, ! 注意：设定单元局部坐标 VATT, 1, 1, 1, 11 !* VSWEEP,ALL /DEVICE,VECTOR,1 /ESHAPE,1.0 !Displays elements with shapes determined from the real constants or section /REPLO ! 注意：红色代
10、表最大配筋方向，绿色代表其次，蓝色表示最小配筋方向 62.混凝土非线性计算实例（1）- MISO 单压 !MISO单压 FINISH /CLEAR /PREP7 ET,1,SOLID65 R,1 MP,EX,1,30E3 MP,NUXY,1,0.2 TB,CONCR,1,1 TBTEMP,0 TBDATA,1,0.3,0.5,2.5,-1 TB,MISO,1,1,15 TBPT, 100E-6, 3.0 TBPT, 300E-6, 8.3 TBPT, 600E-6,14.6 TBPT, 900E-6,19.1 TBPT,1100E-6,21.0 TBPT,1250E-6,22.0 TBPT,1
11、400E-6,22.6 TBPT,1550E-6,22.8 TBPT,1650E-6,22.7 TBPT,1800E-6,22.3 TBPT,2000E-6,21.4 TBPT,2800E-6,16.8 TBPT,3200E-6,14.7 TBPT,3800E-6,12.3 TBPT,4600E-6,9.9 !TB,MKIN,1 !TBTEMP,STRAIN !TBDATA,600E-6,1100E-6,1600E-6,3000E-6,4500E-6 !TBTEMP,0 !TBDATA,15,21.5,23,16,9 BLOCK,0,50,0,50,0,50 MSHAPE,0,3D MSHKE
12、Y,1 ESIZE,10 7VMESH,ALL NSEL,S,LOC,X,0 D,ALL,UX,0 NSEL,S,LOC,Y,0 D,ALL,UY,0 NSEL,S,LOC,Z,0 D,ALL,UZ,0 NSEL,S,LOC,Y,50 CP,1,UY,ALL NSEL,ALL FINISH /VIEW,1,1,1,1 /REPLOT /SOLU ANTYPE,STATIC AUTOTS,OFF OUTRES,ALL,1 TIME,50 NSUBST,50 D,1,UY,-.25 NSEL,ALL SOLVE FINISH /POST26 NSOL,2,1,U,Y,UY RFORCE,3,1,F
13、,Y,FY ADD,4,2,STRAIN,-1/50 ADD,5,3,STRESS,-1/2500 /AXLAB,X,STRAIN /AXLAB,Y,STRESS XVAR,4 PLVAR,5 FINISH 3.混凝土非线性计算实例（2）- MISO 约束压 8!MISO约束压 FINISH /CLEAR /PREP7 ANTYPE,STATIC ET,1,SOLID65 R,1 MP,EX,1,30E3 MP,NUXY,1,0.2 TB,CONCR,1,1 TBTEMP,0 TBDATA,1,0.3,0.5,2.5,-1 TB,MISO,1,1,15 TBPT, 100E-6, 3.0 TB
14、PT, 300E-6, 8.3 TBPT, 600E-6,14.6 TBPT, 900E-6,19.1 TBPT,1100E-6,21.0 TBPT,1250E-6,22.0 TBPT,1400E-6,22.6 TBPT,1550E-6,22.8 TBPT,1650E-6,22.7 TBPT,1800E-6,22.3 TBPT,2000E-6,21.4 TBPT,2800E-6,16.8 TBPT,3200E-6,14.7 TBPT,3800E-6,12.3 TBPT,4600E-6,9.9 !TB,MKIN,1 !TBTEMP,STRAIN !TBDATA,600E-6,1100E-6,16
15、00E-6,3000E-6,4500E-6 !TBTEMP,0 !TBDATA,15,21.5,23,16,9 BLOCK,0,50,0,50,0,50 ESIZE,10 VMESH,ALL MSHAPE,0,3D NSEL,S,LOC,X,0 D,ALL,UX,0 9NSEL,S,LOC,Y,0 D,ALL,UY,0 NSEL,S,LOC,Z,0 D,ALL,UZ,0 NSEL,S,LOC,Y,50 CP,1,UY,ALL D,ALL,UX,0 D,ALL,UZ,0 NSEL,ALL FINISH /VIEW,1,1,1,1 /REPLOT /SOLU OUTRES,ALL,1 TIME,2
16、0 NSUBST,20,0,20 D,1,UY,-.04 LSWRITE,1 TIME,30 NSUBST,20,0,20 D,1,UY,-.06 LSWRITE,2 TIME,40 NSUBST,100 D,1,UY,-.1 LSWRITE,3 LSSOLVE,1,3 FINISH /POST26 NSOL,2,1,U,Y,UY RFORCE,3,1,F,Y,FY ADD,4,2,STRAIN,-1/50 ADD,5,3,STRESS,-1/2500 /AXLAB,X,STRAIN /AXLAB,Y,STRESS XVAR,4 PLVAR,5 FINISH 104.混凝土非线性计算实例（3）
17、- KINH 滞回 !KINH滞回 FINISH /CLEAR /PREP7 ANTYPE,STATIC ET,1,SOLID65 R,1 MP,EX,1,30E3 MP,NUXY,1,0.2 TB,CONCR,1,1 TBTEMP,0 TBDATA,1,0.3,0.5,2.5,-1 TB,KINH,1,1,10 TBPT, 150E-6, 4.5 TBPT, 600E-6,14.8 TBPT,1000E-6,20.25 TBPT,1300E-6,22.3 TBPT,1480E-6,22.8 TBPT,1620E-6,22.8 TBPT,1800E-6,22.3 TBPT,2000E-6,21
18、.4 TBPT,3500E-6,12.8 TBPT,5000E-6, 9.0 BLOCK,0,50,0,50,0,50 ESIZE,10 VMESH,ALL MSHAPE,0,3D NSEL,S,LOC,X,0 D,ALL,UX,0 NSEL,S,LOC,Y,0 D,ALL,UY,0 NSEL,S,LOC,Z,0 D,ALL,UZ,0 NSEL,S,LOC,Y,50 CP,1,UY,ALL 11NSEL,ALL FINISH /VIEW,1,1,1,1 /REPLOT /SOLU OUTRES,ALL,ALL TIME,20 NSUBST,20,0,20 D,1,UY,-.04 LSWRITE
19、,1 TIME,35 NSUBST,15,0,15 D,1,UY,-.01 LSWRITE,2 TIME,90 NSUBST,55,0,55 D,1,UY,-.12 LSWRITE,3 TIME,125 NSUBST,35,0,35 D,1,UY,-.05 LSWRITE,4 TIME,175 NSUBST,50,0,50 D,1,UY,-.15 LSWRITE,5 LSSOLVE,1,5 FINISH /POST26 NSOL,2,1,U,Y,UY RFORCE,3,1,F,Y,FY ADD,4,2,STRAIN,-1/50 ADD,5,3,STRESS,-1/2500 /AXLAB,X,S
20、TRAIN /AXLAB,Y,STRESS XVAR,4 PLVAR,5 FINISH 125.混凝土非线性计算实例（4） - KINH 压-拉裂 !KINH压拉裂 压到峰值，泄载、反向加载到拉裂 FINISH /CLEAR /PREP7 ANTYPE,STATIC ET,1,SOLID65 R,1 MP,EX,1,30E3 MP,NUXY,1,0.2 TB,CONCR,1,1 TBTEMP,0 TBDATA,1,0.3,0.5,2.5,-1 TB,KINH,1,1,10 TBPT, 150E-6, 4.5 TBPT, 600E-6,14.8 TBPT,1000E-6,20.25 TBPT,1
21、300E-6,22.3 TBPT,1480E-6,22.8 TBPT,1620E-6,22.8 TBPT,1800E-6,22.3 TBPT,2000E-6,21.4 TBPT,3500E-6,12.8 TBPT,5000E-6, 9.0 BLOCK,0,50,0,50,0,50 ESIZE,10 VMESH,ALL MSHAPE,0,3D NSEL,S,LOC,X,0 D,ALL,UX,0 NSEL,S,LOC,Y,0 D,ALL,UY,0 NSEL,S,LOC,Z,0 D,ALL,UZ,0 13NSEL,S,LOC,Y,50 CP,1,UY,ALL NSEL,ALL FINISH /VIE
22、W,1,1,1,1 /REPLOT /SOLU OUTRES,ALL,ALL TIME,20 NSUBST,20,0,20 D,1,UY,-.04 LSWRITE,1 TIME,45 NSUBST,25,0,25 D,1,UY,.01 LSWRITE,2 LSSOLVE,1,2 FINISH /POST26 NSOL,2,1,U,Y,UY RFORCE,3,1,F,Y,FY ADD,4,2,STRAIN,-1/50 ADD,5,3,STRESS,-1/2500 /AXLAB,X,STRAIN /AXLAB,Y,STRESS XVAR,4 PLVAR,5 FINISH 6.混凝土非线性计算实例（
23、5） FINISH /CLEAR /PREP7 ANTYPE,STATIC 14ET,1,SOLID65 R,1 MP,EX,1,30E3 MP,NUXY,1,0.2 TB,CONCR,1,1 TBTEMP,0 TBDATA,1,0.3,0.5,2.5,-1 TB,KINH,1,1,18 TBPT, 100E-6, 3.0 TBPT, 300E-6, 8.3 TBPT, 600E-6,14.6 TBPT, 900E-6,19.1 TBPT,1100E-6,21.0 TBPT,1250E-6,22.0 TBPT,1400E-6,22.6 TBPT,1550E-6,22.8 TBPT,1650E-
24、6,22.7 TBPT,1800E-6,22.3 TBPT,2000E-6,21.4 TBPT,2400E-6,19.1 TBPT,2800E-6,16.8 TBPT,3200E-6,14.7 TBPT,3600E-6,13.0 TBPT,4100E-6,11.3 TBPT,4600E-6,9.9 BLOCK,0,50,0,50,0,50 ESIZE,10 VMESH,ALL MSHAPE,0,3D NSEL,S,LOC,X,0 D,ALL,UX,0 NSEL,S,LOC,Y,0 D,ALL,UY,0 NSEL,S,LOC,Z,0 D,ALL,UZ,0 NSEL,S,LOC,Y,50 CP,1
25、,UY,ALL NSEL,ALL FINISH 15/VIEW,1,1,1,1 /REPLOT /SOLU OUTRES,ALL,ALL TIME,50 NSUBST,50,0,50 D,1,UY,-.25 SOLVE FINISH /POST26 NSOL,2,1,U,Y,UY RFORCE,3,1,F,Y,FY ADD,4,2,STRAIN,-1/50 ADD,5,3,STRESS,-1/2500 /AXLAB,X,STRAIN /AXLAB,Y,STRESS XVAR,4 PLVAR,5 FINISH 7.混凝土非线性计算实例（6） !MISO约束压 FINISH /CLEAR /PRE
26、P7 ANTYPE,STATIC ET,1,SOLID65 R,1 MP,EX,1,30E3 MP,NUXY,1,0.2 TB,CONCR,1,1 TBTEMP,0 TBDATA,1,0.3,0.5,2.5,-1 TB,KINH,1,1,15 TBPT, 100E-6, 3.0 16TBPT, 300E-6, 8.3 TBPT, 600E-6,14.6 TBPT, 900E-6,19.1 TBPT,1100E-6,21.0 TBPT,1250E-6,22.0 TBPT,1400E-6,22.6 TBPT,1550E-6,22.8 TBPT,1650E-6,22.7 TBPT,1800E-6,2
27、2.3 TBPT,2000E-6,21.4 TBPT,2800E-6,16.8 TBPT,3200E-6,14.7 TBPT,3800E-6,12.3 TBPT,4600E-6,9.9 !TB,MKIN,1 !TBTEMP,STRAIN !TBDATA,600E-6,1100E-6,1600E-6,3000E-6,4500E-6 !TBTEMP,0 !TBDATA,15,21.5,23,16,9 BLOCK,0,50,0,50,0,50 ESIZE,10 VMESH,ALL MSHAPE,0,3D NSEL,S,LOC,X,0 D,ALL,UX,0 NSEL,S,LOC,Y,0 D,ALL,U
28、Y,0 NSEL,S,LOC,Z,0 D,ALL,UZ,0 NSEL,S,LOC,Y,50 CP,1,UY,ALL D,ALL,UX,0 D,ALL,UZ,0 NSEL,ALL FINISH /VIEW,1,1,1,1 /REPLOT /SOLU OUTRES,ALL,1 TIME,20 NSUBST,20,0,20 17D,1,UY,-.04 LSWRITE,1 TIME,30 NSUBST,20,0,20 D,1,UY,-.06 LSWRITE,2 TIME,40 NSUBST,100 D,1,UY,-.1 LSWRITE,3 LSSOLVE,1,3 FINISH /POST26 NSOL
29、,2,1,U,Y,UY RFORCE,3,1,F,Y,FY ADD,4,2,STRAIN,-1/50 ADD,5,3,STRESS,-1/2500 /AXLAB,X,STRAIN /AXLAB,Y,STRESS XVAR,4 PLVAR,5 FINISH 8.混凝土非线性计算实例（7）- MISO 滞回 !MISO滞回 FINISH /CLEAR /PREP7 ANTYPE,STATIC ET,1,SOLID65 R,1 MP,EX,1,30E3 MP,NUXY,1,0.2 TB,CONCR,1,1 TBTEMP,0 TBDATA,1,0.3,0.5,2.5,-1 18TB,MISO,1,1,
30、15 TBPT, 100E-6, 3.0 TBPT, 300E-6, 8.3 TBPT, 600E-6,14.6 TBPT, 900E-6,19.1 TBPT,1100E-6,21.0 TBPT,1250E-6,22.0 TBPT,1400E-6,22.6 TBPT,1550E-6,22.8 TBPT,1650E-6,22.7 TBPT,1800E-6,22.3 TBPT,2000E-6,21.4 TBPT,2800E-6,16.8 TBPT,3200E-6,14.7 TBPT,3800E-6,12.3 TBPT,4600E-6,9.9 BLOCK,0,50,0,50,0,50 ESIZE,1
31、0 VMESH,ALL MSHAPE,0,3D NSEL,S,LOC,X,0 D,ALL,UX,0 NSEL,S,LOC,Y,0 D,ALL,UY,0 NSEL,S,LOC,Z,0 D,ALL,UZ,0 NSEL,S,LOC,Y,50 CP,1,UY,ALL NSEL,ALL FINISH /VIEW,1,1,1,1 /REPLOT /SOLU OUTRES,ALL,ALL TIME,20 NSUBST,20,0,20 D,1,UY,-.04 LSWRITE,1 TIME,35 NSUBST,15,0,15 D,1,UY,-.01 19LSWRITE,2 TIME,60 NSUBST,25,0
32、,25 D,1,UY,-.06 LSWRITE,3 TIME,80 NSUBST,20,0,20 D,1,UY,-.02 LSWRITE,4 TIME,110 NSUBST,30,0,30 D,1,UY,-.08 LSWRITE,5 LSSOLVE,1,5 FINISH /POST26 NSOL,2,1,U,Y,UY RFORCE,3,1,F,Y,FY ADD,4,2,STRAIN,-1/50 ADD,5,3,STRESS,-1/2500 /AXLAB,X,STRAIN /AXLAB,Y,STRESS XVAR,4 PLVAR,5 FINISH 9.混凝土非线性计算实例（8） !MISO单压
33、FINISH /CLEAR /PREP7 ANTYPE,STATIC ET,1,SOLID65 R,1 20MP,EX,1,30E3 MP,NUXY,1,0.2 TB,CONCR,1,1 TBTEMP,0 TBDATA,1,0.3,0.5,2.5,-1 TB,MISO,1,1,10 TBPT, 150E-6, 4.5 TBPT, 600E-6,14.8 TBPT,1000E-6,20.25 TBPT,1300E-6,22.3 TBPT,1480E-6,22.8 TBPT,1620E-6,22.8 TBPT,1800E-6,22.3 TBPT,2000E-6,21.4 TBPT,3500E-6,
34、12.8 TBPT,5000E-6, 9.0 !TB,MKIN,1 !TBTEMP,STRAIN !TBDATA,600E-6,1100E-6,1600E-6,3000E-6,4500E-6 !TBTEMP,0 !TBDATA,15,21.5,23,16,9 BLOCK,0,50,0,50,0,50 ESIZE,10 VMESH,ALL MSHAPE,0,3D NSEL,S,LOC,X,0 D,ALL,UX,0 NSEL,S,LOC,Y,0 D,ALL,UY,0 NSEL,S,LOC,Z,0 D,ALL,UZ,0 NSEL,S,LOC,Y,50 CP,1,UY,ALL NSEL,ALL FIN
35、ISH /VIEW,1,1,1,1 /REPLOT /SOLU OUTRES,ALL,1 21TIME,50 NSUBST,100,0,100 D,1,UY,-.25 SOLVE FINISH /POST26 NSOL,2,1,U,Y,UY RFORCE,3,1,F,Y,FY ADD,4,2,STRAIN,-1/50 ADD,5,3,STRESS,-1/2500 /AXLAB,X,STRAIN /AXLAB,Y,STRESS XVAR,4 PLVAR,5 FINISH 10.混凝土非线性计算实例（9）-梁平面应力 !梁平面应力 FINISH /CLEAR /PREP7 ET,1,SOLID65
36、 R,1 MP,EX,1,30E3 MP,NUXY,1,0.2 TB,CONCR,1,1 TBTEMP,0 TBDATA,1,0.3,0.5,2.5,-1 TB,MISO,1,1,15 TBPT, 100E-6, 3.0 TBPT, 300E-6, 8.3 TBPT, 600E-6,14.6 TBPT, 900E-6,19.1 TBPT,1100E-6,21.0 TBPT,1250E-6,22.0 TBPT,1400E-6,22.6 TBPT,1550E-6,22.8 TBPT,1650E-6,22.7 TBPT,1800E-6,22.3 TBPT,2000E-6,21.4 22TBPT,28
37、00E-6,16.8 TBPT,3200E-6,14.7 TBPT,3800E-6,12.3 TBPT,4600E-6, 9.9 ET,2,LINK8 R,2,30 R,3,6 MP,EX,2,200E3 MP,NUXY,2,0.3 TB,BKIN,2,1 TBTEMP,0 TBDATA,300,0 /ESHAPE,1 /TRIAD,LBOT /VIEW,1,1,1,1 N,1,0,0,0 N,2,20,0,0 NGEN,11,2,1,2,1,0,20 NGEN,43,22,1,22,1,0,0,20 E,1,2,4,3,23,24,26,25 EGEN,10,2,1 TYPE,2 REAL,
38、2 MAT,2 E,3,25 E,4,26 REAL,3 E,19,41 E,20,42 EGEN,42,22,1,14 E,47,63 E,48,64 EGEN,10,88,589,590 /NUMBER,1 /PNUM,TYPE,1 /REPLOT NSEL,S,LOC,Z,0 D,ALL,UZ,0 23NSEL,S,LOC,Y,0 NSEL,R,LOC,Z,740 D,ALL,UY,0 NSEL,ALL FINISH /SOLU ANTYPE,STATIC NSEL,S,LOC,Y,200 NSEL,R,LOC,Z,200 D,ALL,UY,-2 TIME,200 NSUBST,10 L
39、SWRITE,1 D,ALL,UY,-0.21 TIME,21 NSUBST,100 LSWRITE,2 D,ALL,UY,-1 TIME,100 NSUBST,100 LSWRITE,3 NSEL,ALL AUTOTS,1 OUTRES,ALL,1 PRED,ON LSSOLVE,1,3 FINISH 11.四层弹簧质点模型的地震分析 包括频率分析，谱分析和时程分析 作者：陆新征，清华大学土木工程系 FINI /CLEAR /UNITS,SI /PREP7 DAMPRATIO=0.02 !振型阻尼系数都为0.02 24ET,1,BEAM188 KEYOPT,1,7,1 KEYOPT,1,8,
40、1 ET,2,MASS21,4 MP,EX,1,210E9 MP,NUXY,1,0.3 MP,DENS,1,7800 MPDATA,DAMP,1,DAMPRATIO SECTYPE,3,BEAM,I SECDATA,0.2,0.2,0.5,0.01,0.01,0.008 R,1,1.6E2,1.6E2 R,2,1.2E2,1.2E2 K,1,0,0,0 K,2,0,3,0 K,3,0,6,0 K,4,0,9,0 K,5,0,12,0 K,6,0,0,10000 L,1,2 L,2,3 L,3,4 L,4,5 LSEL,1,4 LATT,1,0,1,6,3 ALLS LESIZE,ALL, ,
41、,3 , , , , ,1 LMESH,ALL NUMMRG,ALL NUMCMP,ALL TYPE,2 REAL,1 E,2 TYPE,2 REAL,2 E,8 E,14 E,20 *SET,NT,1001 *SET,DT,0.02 *DIM,AC,NT *VREAD,AC(1),RECORD,TXT (F8.3) /SOLU D,1,ALL, 25ANTYPE,0 SOLVE FINISH /SOLU !模态分析 ANTYPE,2 MODOPT,SUBSP,8 MXPAND,8, , ,1 SOLVE FINI !得到自振频率1 *GET,FREQ1,MODE,1,FREQ /ESHAPE
42、,1.0 /SOLU ANTYPE,SPECTR !谱分析 SPOPT,SPRS MDAMP,1,DAMPRATIO, , , , , , SED,1 SVTYPE,2,1 FREQ,0.167,0.25,0.333,0.44444,0.5,0.667,1,1.25,1.667 SV,0.154350625,0.191590625,0.210210625,0.224175625,0.250716714,0.329514922,0.484352764,0.598723486,0.786897371 FREQ,2.222,10,100000 SV,1.034212766,1.034212766,0
43、.3528 !按7度多遇地震，地震影响系数0.08 第一组III场地TG0.45 SRSS,0.01,DISP SOLVE FINISH /POST1 /INP,MCOM NSORT,U,SUM *GET,DMAX_SPE,SORT, ,MAX !反应谱最大变形 !开始地震计算 FINISH /CONFIG,NRES,20000 /SOLU ANTYPE,TRANS TRNOPT,FULL ALPHAD,2*DAMPRATIO*FREQ1*2*3.1415926 BETAD,2*DAMPRATIO/(FREQ1*2*3.1415926) *DO,I,1,1001 ACEL,AC(I),0,0
44、TIME,I*0.02 26OUTRES,ALL,ALL SOLVE *ENDDO FINISH /POST26 FILE NSOL,2,20,U,X, PLVAR,2, , , , , , , , , , STORE,MERGE XVAR,1 PLVAR,2, *GET,DMAX_TRANS,VARI,2,EXTREM,VMAX, , ! 时程分析最大位移 数据文件 record.txt -0.0014 -0.0108 -0.0101 -0.0088 -0.0095 -0.012 -0.0142 -0.0128 -0.011 -0.0085 -0.0085 -0.0131 -0.0176 -
45、0.0194 -0.0162 -0.0144 -0.0108 -0.0082 -0.0042 -0.0066 -0.0131 -0.019 -0.0196 -0.0066 0.003 0.0141 -0.0049 27-0.0128 -0.0144 -0.0203 -0.026 -0.0325 -0.0306 -0.0172 -0.0197 -0.0163 -0.0164 -0.0067 0.0025 0.015 0.0236 0.0252 0.0336 0.0463 0.0492 0.0419 0.0359 0.0271 0.0235 0.0339 0.0412 0.053 0.0639 0
46、.0732 0.0652 0.0599 0.04 0.04 0.0063 -0.0515 -0.0787 -0.0603 -0.0484 -0.025 -0.0059 0.0134 0.0308 0.0499 0.071 0.0995 0.1219 280.1529 0.1449 0.1155 0.0935 0.0892 0.0926 0.0839 0.0901 0.0993 0.1209 0.0328 -0.1475 -0.2066 -0.1989 -0.2034 -0.1816 -0.1725 -0.1752 -0.1753 -0.1805 -0.163 -0.1347 -0.1087 -
47、0.0782 -0.0429 -0.0017 0.036 0.0785 0.1164 0.1598 0.196 0.2412 0.2729 0.3036 0.32 0.3417 0.2821 0.2324 -0.1198 -0.2373 -0.164 -0.1865 -0.1095 -0.0753 29-0.0173 0.0113 0.0533 0.0895 0.1186 0.1757 0.0576 -0.2631 -0.1547 -0.1729 -0.1012 -0.0579 0.0237 -0.067 -0.198 -0.1641 -0.1685 -0.1481 -0.1231 -0.10
48、01 -0.0751 -0.0523 -0.0271 -0.0044 0.0188 -0.0095 -0.0433 -0.0838 -0.0951 -0.0716 -0.0599 -0.0334 -0.0108 0.0185 0.042 0.0673 -0.0097 -0.0372 -0.004 0.0011 0.0344 0.0565 0.0883 0.113 300.1363 0.0219 0.0241 0.0683 0.0689 0.1318 0.1353 0.204 -0.0931 -0.1308 -0.0692 -0.0546 0.0072 0.0675 -0.1067 -0.1488 -0.1071 -0.1162 -0.0762 -0.0559 -0.0215 -0.0126 -0.0674 -0.0324 -0.0337 -0.0109 0.0017 0.0299 0.0488 0.0608 0.0222 -0.0032 -0.0245 0.0077 0.0211 0.0568 0.0826 0.1206 0.1478 0.1737 0.0421 0.0029 0.0259 0.0293
