ANSYS土木工程经典实例命令流大全.doc

______________________________________________________________________________________________________________
l 1．ANSYS SOLID65环向布置钢筋的例子 3
l 2.混凝土非线性计算实例（1）- MISO单压 5
l 3.混凝土非线性计算实例（2）- MISO约束压 6
l 4.混凝土非线性计算实例（3）- KINH滞回 9
l 5.混凝土非线性计算实例（4）- KINH压-拉裂 11
l 6.混凝土非线性计算实例（5） 12
l 7.混凝土非线性计算实例（6） 14
l 8.混凝土非线性计算实例（7）- MISO滞回 16
l 9.混凝土非线性计算实例（8） 18
l 10.混凝土非线性计算实例（9）-梁平面应力 20
l 11.四层弹簧－质点模型的地震分析 22
l 12.悬臂梁地震分析 48
l 13.用beam 54单元描述变截面梁的例子 72
l 14.变截面梁实例 73
l 15.拱桥浇筑过程分析-单元生死应用实例 74
l 16.简支梁实体与预应力钢筋分析实例 75
l 17. 简单的二维焊接分析－单元生死实例 77
l 18.隧道开挖（三维）的命令流 84
l 19.岩土接触分析实例 101
l 20.钢筋混凝土管的动力响应特性分析实例 109
l 21.隧道模拟开挖命令流（入门） 116
l 22.螺栓连接的模拟实现问题 119
l 23.道路的基层、垫层模量与应力之间的关系 129
l 23.滞回分析 151
l 24.模拟某楼层浇注 153
l 25.在面上施加移动的面力 155
l 27.在任意面施加任意方向任意变化的压力 159
l 28.预紧分析 160
l 29.几何非线性+塑性+接触+蠕变 162
l 30.埋设在地下的排水管道 167
l 32.幕墙企业玻璃简化计算 172
l 33.等截面杆单元生死应用实例 188
l 34.梁板建模联系 189
l 36.简单的例子－如何对结构的振动控制分析 192
l 37.模态分析结果的输出实例 194
l 38.火车过桥动态加载实例（部分） 196
l 39.悬索结构的找形和计算的例题 213
l 40.陶瓷杆撞击铝板的例子 218
l 41.求反作用力的APDL命令法 221
l 42.LS-DYNA实例（部分） 222
l 43.路面分层填筑对路基的影响 223
l 44.一个例子（含地震影响，求振兴与频率） 227
l 45.接触面上的压力总和 231
l 46.施加位置函数荷载 235
l 47.非线性分析考虑刚度退化 236
l 48.一个圆形水池的静力分析 237
l 49.ANSYS中混凝土模式预应力模拟的算例 238
l 50.悬臂梁受重力作用发生大变形求其固有频率 240
l 51.循环对称结构模态分析 242
l 52.三角平台受谐波载荷作用的结构响应 244
l 53.三角平台受一地震谱激励的应力分布和支反力 246
l 54.三角平台受时程载荷作用的应力分布和变形过程 248
l 55.经典层合板理论 250
l 56．定易圆轨迹的例子 257
l 57．模拟门式刚架施工－单元生死 257
l 58．钢筋混凝土整体式模型例子 260
l 59．在荷载步之间改变材料属性例子 262
l 60．含预应力的特征值屈曲计算 263
l 61．振型叠加计算及工况组合例子 265
l 62．柱子稳定分析算（预应力，特征值屈曲，初始缺陷） 268
l 63． module MConcrete !混凝土模板 271
l 64.混凝土开裂实例 279
l 65．螺栓网格划分 280
l 66．自由液面的土石坝平面渗流分析 281
l 67．导出刚度矩阵 285
l 68．某混凝土拱坝工程施工期及运行期温度场仿真分析 286
l 69．移动温度荷载计算 293
l 70．SHSD用于壳-实体装配实例An 295
l 71．ansys显示－隐式－回弹分析实例 299
l 72．工况组合的经典例子 314

1．ANSYS 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 real 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,,,,,,,,
MPTEMP,1,0
!钢材基本属性
MPDATA,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,
VSBW, 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
! 注意：红色代表最大配筋方向，绿色代表其次，蓝色表示最小配筋方向

2.混凝土非线性计算实例（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,,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,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
MSHKEY,1
ESIZE,10
VMESH,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,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约束压

!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
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
!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
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
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

4.混凝土非线性计算实例（3）- 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.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
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
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,STRAIN
/AXLAB,Y,STRESS
XVAR,4
PLVAR,5
FINISH

5.混凝土非线性计算实例（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,,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,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
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,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.混凝土非线性计算实例（5）

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,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-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,UY,ALL
NSEL,ALL
FINISH
/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
/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,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
!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
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
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

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
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
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,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
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
LSWRITE,2
TIME,60
NSUBST,25,0,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单压
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,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,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
FINISH
/VIEW,1,1,1,1
/REPLOT
/SOLU
OUTRES,ALL,1
TIME,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
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
TBPT,,2800E-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,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
NSEL,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
LSWRITE,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
ET,1,BEAM188
KEYOPT,1,7,1
KEYOPT,1,8,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, , ,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,
ANTYPE,0
SOLVE
FINISH
/SOLU !模态分析
ANTYPE,2
MODOPT,SUBSP,8
MXPAND,8, , ,1
SOLVE
FINI
!得到自振频率1
*GET,FREQ1,MODE,1,FREQ
/ESHAPE,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.3528
!按7度多遇地震，地震影响系数0.08 第一组III场地TG＝0.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
TIME,I*0.02
OUT