ANSYS经典APDL程序实例

TOC\o"1-1"\h\z\u1．ANSYSSOLID65环向布置钢筋的例子 22.混凝土非线性计算实例（1）-MISO单压 33.混凝土非线性计算实例（2）-MISO约束压 44.混凝土非线性计算实例（3）-KINH滞回 55.混凝土非线性计算实例（4）-KINH压-拉裂 66.混凝土非线性计算实例（5） 77.混凝土非线性计算实例（6） 88.混凝土非线性计算实例（7）-MISO滞回 99.混凝土非线性计算实例（8） 1010.混凝土非线性计算实例（9）-梁平面应力 1111.四层弹簧－质点模型的地震分析 1212.悬臂梁地震分析 2513.用beam54单元描述变截面梁的例子 3714.变截面梁实例 3815.拱桥浇筑过程分析-单元生死应用实例 3816.简支梁实体与预应力钢筋分析实例 3917.

简单的二维焊接分析－单元生死实例 4018.隧道开挖（三维）的命令流 4419.岩土接触分析实例 5320.钢筋混凝土管的动力响应特性分析实例 5721.隧道模拟开挖命令流（入门） 6122.螺栓连接的模拟实现问题 6323.道路的基层、垫层模量与应力之间的关系 6723.滞回分析 8024.模拟某楼层浇注 8125.在面上施加移动的面力 8327.在任意面施加任意方向任意变化的压力 8528.预紧分析 8629.几何非线性+塑性+接触+蠕变 8730.埋设在地下的排水管道 8932.幕墙企业玻璃简化计算 9233.等截面杆单元生死应用实例 10034.梁板建模联系 10135.圆盘转动例子 10136.简单的例子－如何对结构的振动控制分析 10237.模态分析结果的输出实例 10338.火车过桥动态加载实例（部分） 10439.悬索结构的找形和计算的例题 11340.陶瓷杆撞击铝板的例子 11541.求反作用力的APDL命令法 117实例（部分） 11843.路面分层填筑对路基的影响 11844.一个例子（含地震影响，求振兴与频率） 12145.接触面上的压力总和 12346.施加位置函数荷载 12547.非线性分析考虑刚度退化 12548.一个圆形水池的静力分析 126中混凝土模式预应力模拟的算例 12750.悬臂梁受重力作用发生大变形求其固有频率 12851.循环对称结构模态分析 12952.三角平台受谐波载荷作用的结构响应 13053.三角平台受一地震谱激励的应力分布和支反力 13154.三角平台受时程载荷作用的应力分布和变形过程 13255.经典层合板理论 13356．定易圆轨迹的例子 13757．模拟门式刚架施工－单元生死 13758．钢筋混凝土整体式模型例子 13859．在荷载步之间改变材料属性例子 13960．含预应力的特征值屈曲计算 14061．振型叠加计算及工况组合例子 14162．柱子稳定分析算（预应力，特征值屈曲，初始缺陷） 14363．

moduleMConcrete!混凝土模板 14464.混凝土开裂实例 14865．螺栓网格划分 14966．自由液面的土石坝平面渗流分析 15067．导出刚度矩阵 15268．某混凝土拱坝工程施工期及运行期温度场仿真分析 15269．移动温度荷载计算 15670．SHSD用于壳-实体装配实例An 15771．ansys显示－隐式－回弹分析实例 16072．工况组合的经典例子 1681．ANSYSSOLID65环向布置钢筋的例子!一个管道，环向配筋率为1%，纵向配筋率为0.5%，径向配筋率为0.1%!FINISH/CLEAR/PREP7!*!单元属性ET,1,SOLID65!*KEYOPT,1,1,0！：单元主要选项设置KEYOPT,1,5,0！5：NumberoftheKEYOPTtobedefined(KEYOPT(KNUM)).KEYOPT,1,6,0！0：ValueofthisKEYOPTKEYOPT,1,7,1!*!实参数1：不同方向配筋R,1,2,.001,,,2,.01,RMORE,90,,2,.005,90,90,!Addsrealconstantstoaset.!材料属性!混凝土基本材料属性MPTEMP,,,,,,,,！给材料属性定义温度表MPTEMP,1,0MPDATA,EX,1,,30E2！给温度表定义相应的温度数据MPDATA,PRXY,1,,.2!屈服准则TB,MISO,1,1,5,！Activatesadatatablefornonlinearmaterialpropertiesorspecialelementinput.TBTEMP,0！Definesatemperatureforthedatatable.TBPT,,0.0005,15！Definesapointonanonlineardatacurve在应力应变区线上定义一点TBPT,,0.001,21TBPT,,0.0015,24TBPT,,0.002,27TBPT,,0.003,24!破坏准则TB,CONC,1,1,9,TBTEMP,0TBDATA,,.5,.9,3,30,,!抗拉和抗压是10倍的关系TBDATA,,,,1,,,MPTEMP,,,,,,,,MPTEMP,1,0!钢材基本属性MPDATA,EX,2,,200E3MPDATA,PRXY,2,,.27!屈服准则TB,BISO,2,1,2,TBTEMP,0TBDATA,,310,2E3,,,,!管道内径和外径CYL4,0,0,3000,,,,10000CYL4,0,0,2000,,,,10000VSBV,1,2!定义局部柱坐标CSWPLA,11,1,1,1,!Definesalocalcoordinatesystemattheoriginoftheworkingplane.KWPAVE,11!movetokeypoint11WPRO,,-90.000000,VSBW,3WPCSYS,-1,0!Definestheworkingplanelocationbasedonacoordinatesystem.KWPAVE,1WPRO,,,-90.000000VSBW,ALLESIZE,500,0,!注意：设定单元局部坐标VATT,1,1,1,11!*VSWEEP,ALL/DEVICE,VECTOR,1!Displayselementswithshapesdeterminedfromtherealconstantsorsection/REPLO!注意：红色代表最大配筋方向，绿色代表其次，蓝色表示最小配筋方向2.混凝土非线性计算实例（1）-MISO单压

!MISO单压

FINISH

/CLEAR

/PREP7

ET,1,SOLID65

R,1

MP,EX,1,30E3

MP,NUXY,1,0.2TB,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,9BLOCK,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

FINISH3.混凝土非线性计算实例（2）-MISO约束压

!MISO约束压

FINISH

/CLEAR

/PREP7

ANTYPE,STATICET,1,SOLID65R,1MP,EX,1,30E3

MP,NUXY,1,0.2TB,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,9BLOCK,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

FINISH4.混凝土非线性计算实例（3）-KINH滞回

!KINH滞回

FINISH

/CLEAR

/PREP7

ANTYPE,STATICET,1,SOLID65R,1MP,EX,1,30E3

MP,NUXY,1,0.2TB,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.0BLOCK,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,STATICET,1,SOLID65R,1MP,EX,1,30E3

MP,NUXY,1,0.2TB,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.0BLOCK,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,STATICET,1,SOLID65R,1MP,EX,1,30E3

MP,NUXY,1,0.2TB,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.9BLOCK,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,0NSEL,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,STATICET,1,SOLID65R,1MP,EX,1,30E3

MP,NUXY,1,0.2

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,9BLOCK,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,STATICET,1,SOLID65R,1MP,EX,1,30E3

MP,NUXY,1,0.2TB,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.9BLOCK,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,35NSUBST,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,FYADD,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,STATICET,1,SOLID65R,1MP,EX,1,30E3

MP,NUXY,1,0.2TB,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,9BLOCK,0,50,0,50,0,50ESIZE,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.9ET,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,1N,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

/REPLOTNSEL,S,LOC,Z,0

D,ALL,UZ,0NSEL,S,LOC,Y,0

NSEL,R,LOC,Z,740

D,ALL,UY,0

NSEL,ALLFINISH/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

FINISH11.四层弹簧－质点模型的地震分析

包括频率分析，谱分析和时程分析

作者：陆新征，清华大学土木工程系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

OUTRES,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,,!时程分析最大位移数据文件-0.0014

-0.0108

-0.0101

-0.0088

-0.0095

-0.012

-0.0142

-0.0128

-0.011

-0.0085

-0.0085

-0.finish

/clear12.悬臂梁地震分析!定义单元类型

/prep7

ET,1,BEAM3

R,1,0.09,0.000675,0.3,!0.3*截面

MP,EX,1,2.06e11

MP,NUXY,1,0.3

MP,DAMP,1,0.02!阻尼比0.02

MP,dens,1,7.85e3

!

!定义节点单元和位移约束

*do,i,1,11

n,i,(i-1)*0.5,0,0!5m长,分10段

*enddo

*do,i,1,10

e,i,i+1

*enddo

d,1,ALL,0

!

!定义和读入时程曲线

NT=1000!时程曲线有NT个点

DT=0.01!时间间隔

*dim,ac,,NT

!读入数据,这个数据文件可以用excel等软件来写

/input,tianjin,txt!天津波东西向共10秒

!

/SOLU

NSUBST,1,,,1!1个子步

OUTRES,ALL,1!输出每个子步的结果

ANTYPE,TRANS!时程分析

*do,i,1,NT

ACEL,0,ac,0

TIME,i*DT

solve

*enddo

!

/POST26

NSOL,2,11,U,Y,

PLVAR,2,3,,,下面是文件tianjin.txt，这个文件可以用excel等软件制作或编写很短的程序就能写出来先

贴前500个点

ac(1)=-0.06334598

ac(2)=-0.04417088

ac(3)=-0.02188456

ac(4)=-0.00622243

ac(5)=0.01599961

ac(6)=0.03431334

ac(7)=0.05332774

ac(8)=0.07603510

ac(9)=

ac(10)=0.11568701

ac(11)=0.13706695

ac(12)=0.15446468

ac(13)=0.17966937

ac(14)=0.16881546

ac(15)=0.10252541

ac(16)=0.03081323

ac(17)=-0.02835768

ac(18)=-0.04577148

ac(19)=0.00371867

ac(20)=0.05160179

ac(21)=0.08542017

ac(22)=0.14456214

ac(23)=0.15453540

ac(24)=0.07974091

ac(25)=-0.00708701

ac(26)=-0.08872422

ac(27)=-0.15415934

ac(28)=-0.16546641

ac(29)=-0.15415934

ac(30)=-0.15538390

ac(31)=-0.15533569

ac(32)=-0.15030567

ac(33)=-0.14931896

ac(34)=-0.14984284

ac(35)=-0.14375542

ac(36)=-0.14322831

ac(37)=-0.14130631

ac(38)=-0.07721145

ac(39)=0.02666387

ac(40)=0.06092257

ac(41)=0.03247169

ac(42)=0.01391368

ac(43)=-0.00049497

ac(44)=-0.02224453

ac(45)=-0.04070290

ac(46)=-0.05431446

ac(47)=-0.05096541

ac(48)=-0.03585931

ac(49)=-0.02533004

ac(50)=-0.01427366

ac(51)=-0.00180630

ac(52)=0.00995074

ac(53)=0.02199384

ac(54)=0.03315949

ac(55)=0.04807596

ac(56)=0.05553902

ac(57)=0.03931121

ac(58)=-0.00056568

ac(59)=-0.05451052

ac(60)=-0.07338993

ac(61)=-0.02214490

ac(62)=0.05202605

ac(63)=0.11738081

ac(64)=0.17217417

ac(65)=0.18113498

ac(66)=0.15603314

ac(67)=0.14236373

ac(68)=0.12819293

ac(69)=0.10743007

ac(70)=0.09416241

ac(71)=0.07652363

ac(72)=0.05719426

ac(73)=0.04319058

ac(74)=0.02332124

ac(75)=0.00492073

ac(76)=-0.00940756

ac(77)=-0.04927159

ac(78)=-0.13574918

ac(79)=-0.23619831

ac(80)=-0.28188944

ac(81)=-0.26245400

ac(82)=-0.23472306

ac(83)=-0.21773353

ac(84)=-0.19507116

ac(85)=-0.17273341

ac(86)=-0.15239482

ac(87)=-0.12869753

ac(88)=-0.10692225

ac(89)=-0.08414096

ac(90)=-0.06152039

ac(91)=-0.05068257

ac(92)=-0.04799239

ac(93)=-0.04467870

ac(94)=-0.04122037

ac(95)=-0.03772668

ac(96)=-0.03407871

ac(97)=-0.03062681

ac(98)=-0.02680529

ac(99)=-0.02322804

ac(100)=-0.01947723

ac(101)=-0.01553678

ac(102)=-0.01203667

ac(103)=-0.00798695

ac(104)=-0.00334905

ac(105)=-0.00462182

ac(106)=-0.03004507

ac(107)=-0.08054765

ac(108)=-0.10697689

ac(109)=-0.07439592

ac(110)=-0.01944508

ac(111)=0.02653852

ac(112)=0.07484268

ac(113)=0.12548025

ac(114)=0.17680563

ac(115)=0.21375450

ac(116)=0.18098070

ac(117)=0.08950845

ac(118)=0.03323984

ac(119)=0.05083042

ac(120)=0.11876930

ac(121)=0.18590787

ac(122)=0.17401262

ac(123)=0.08684078

ac(124)=-0.00345190

ac(125)=-0.09242682

ac(126)=-0.17069249

ac(127)=-0.17872764

ac(128)=-0.13418716

ac(129)=-0.09391815

ac(130)=-0.05721033

ac(131)=-0.01592890

ac(132)=0.02479650

ac(133)=0.06137897

ac(134)=0.08479986

ac(135)=0.10140371

ac(136)=0.12139840

ac(137)=0.14102989

ac(138)=0.16063246

ac(139)=0.18017076

ac(140)=0.20074397

ac(141)=0.21970373

ac(142)=0.21152073

ac(143)=0.16972826

ac(144)=0.12305363

ac(145)=0.08278464

ac(146)=0.03976118

ac(147)=-0.00771375

ac(148)=-0.04257670

ac(149)=-0.03893517

ac(150)=-0.00746948

ac(151)=0.01970221

ac(152)=0.04440550

ac(153)=0.07499374

ac(154)=0.09928563

ac(155)=0.10520273

ac(156)=0.09103835

ac(157)=0.07317137

ac(158)=0.06160396

ac(159)=0.04555935

ac(160)=0.02963688

ac(161)=0.01710525

ac(162)=-0.00092565

ac(163)=-0.01085068

ac(164)=-0.01018537

ac(165)=-0.01261841

ac(166)=-0.01169919

ac(167)=-0.01216202

ac(168)=-0.01623424

ac(169)=-0.01104031

ac(170)=-0.02329875

ac(171)=-0.07698646

ac(172)=-0.15369652

ac(173)=-0.22725999

ac(174)=-0.25317177

ac(175)=-0.22106008

ac(176)=-0.18460616

ac(177)=-0.15669201

ac(178)=-0.12168124

ac(179)=-0.09077158

ac(180)=-0.05768922

ac(181)=-0.02212883

ac(182)=0.00782303

ac(183)=0.04577148

ac(184)=0.08010089

ac(185)=0.11013953

ac(186)=0.15554783

ac(187)=0.17380692

ac(188)=0.12869753

ac(189)=0.03899945

ac(190)=-0.06785852

ac(191)=-0.12427820

ac(192)=-0.09443562

ac(193)=-0.04283704

ac(194)=-0.00661454

ac(195)=0.03005471

ac(196)=0.07215894

ac(197)=0.11737439

ac(198)=0.12984817

ac(199)=0.08094940

ac(200)=0.00863940

ac(201)=-0.03943013

ac(202)=-0.03732813

ac(203)=0.00579817

ac(204)=0.05437874

ac(205)=0.10425458

ac(206)=0.12601379

ac(207)=0.05491228

ac(208)=-0.07981163

ac(209)=-0.17312552

ac(210)=-0.16219451

ac(211)=-0.08033552

ac(212)=-0.00402401

ac(213)=0.06996374

ac(214)=0.14205196

ac(215)=0.12564416

ac(216)=0.02699170

ac(217)=-0.07388489

ac(218)=-0.16207558

ac(219)=-0.19582647

ac(220)=-0.15537426

ac(221)=-0.11322825

ac(222)=-0.08221896

ac(223)=-0.04049720

ac(224)=-0.00473110

ac(225)=0.03313699

ac(226)=0.06339419

ac(227)=0.07094402

ac(228)=0.08264322

ac(229)=0.09398243

ac(230)=0.10272148

ac(231)=0.11706583

ac(232)=0.08520804

ac(233)=-0.00966790

ac(234)=-0.12072986

ac(235)=-0.20040973

ac(236)=-0.19588110

ac(237)=-0.13388182

ac(238)=-0.08552302

ac(239)=-0.03347447

ac(240)=0.02360408

ac(241)=0.06844670

ac(242)=0.13263476

ac(243)=0.13517709

ac(244)=-0.00161667

ac(245)=-0.12267759

ac(246)=-0.10150655

ac(247)=-0.02301591

ac(248)=0.03731206

ac(249)=0.10196938

ac(250)=0.18120569

ac(251)=0.21782672

ac(252)=0.16884118

ac(253)=0.09337819

ac(254)=0.04204959

ac(255)=-0.01785091

ac(256)=-0.08301926

ac(257)=-0.11096233

ac(258)=-0.10820787

ac(259)=-0.10117551

ac(260)=-0.09735720

ac(261)=-0.09355818

ac(262)=-0.08819069

ac(263)=-0.08428238

ac(264)=-0.07970878

ac(265)=-0.07474948

ac(266)=-0.07102116

ac(267)=-0.06531620

ac(268)=-0.06104149

ac(269)=-0.05688249

ac(270)=-0.05052508

ac(271)=-0.04616360

ac(272)=-0.04679998

ac(273)=-0.08248573

ac(274)=-0.14751908

ac(275)=-0.15090993

ac(276)=-0.08145723

ac(277)=-0.01117530

ac(278)=0.06283816

ac(279)=0.11612411

ac(280)=0.07451484

ac(281)=-0.02657709

ac(282)=-0.12738618

ac(283)=-0.20841272

ac(284)=-0.22112435

ac(285)=-0.17711097

ac(286)=-0.14081776

ac(287)=-0.10858393

ac(288)=-0.06811886

ac(289)=-0.03359339

ac(290)=0.00408186

ac(291)=0.04104038

ac(292)=0.05570615

ac(293)=0.05431446

ac(294)=0.05417304

ac(295)=0.05539117

ac(296)=0.05399305

ac(297)=0.05511476

ac(298)=0.05607256

ac(299)=0.05353987

ac(300)=0.05654502

ac(301)=0.05659966

ac(302)=0.05304490

ac(303)=0.06132433

ac(304)=0.04616360

ac(305)=-0.01192739

ac(306)=-0.07660720

ac(307)=-0.14828083

ac(308)=-0.18884873

ac(309)=-0.12763689

ac(310)=-0.03275131

ac(311)=0.03611965

ac(312)=0.12820899

ac(313)=0.18819949

ac(314)=0.14638130

ac(315)=0.08587656

ac(316)=0.04687712

ac(317)=-0.00992503

ac(318)=-0.05112932

ac(319)=-0.04594182

ac(320)=-0.02787236

ac(321)=-0.01266341

ac(322)=0.00448683

ac(323)=0.02053465

ac(324)=0.03622250

ac(325)=0.05326346

ac(326)=0.07118507

ac(327)=0.08044479

ac(328)=0.05676358

ac(329)=0.00339084

ac(330)=-0.04208173

ac(331)=-0.08426631

ac(332)=-0.14127417

ac(333)=-0.15549961

ac(334)=-0.09603622

ac(335)=-0.02469365

ac(336)=0.03537398

ac(337)=0.10094088

ac(338)=0.16460185

ac(339)=0.23567121

ac(340)=0.28361538

ac(341)=0.24446167

ac(342)=0.16147134

ac(343)=0.09640583

ac(344)=0.02374550

ac(345)=-0.05336631

ac(346)=-0.11976565

ac(347)=-0.20272385

ac(348)=-0.26170835

ac(349)=-0.23775071

ac(350)=-0.17382941

ac(351)=-0.11820362

ac(352)=-0.06053688

ac(353)=-0.00005464

ac(354)=0.05626861

ac(355)=0.11858930

ac(356)=0.17472613

ac(357)=0.18699101

ac(358)=0.15499821

ac(359)=0.11930282

ac(360)=0.09198007

ac(361)=0.05843489

ac(362)=0.02561930

ac(363)=-0.00373474

ac(364)=-0.04124286

ac(365)=-0.06626756

ac(366)=-0.05113896

ac(367)=-0.01100174

ac(368)=0.02835125

ac(369)=0.06458018

ac(370)=0.10447314

ac(371)=0.14322188

ac(372)=0.15820906

ac(373)=0.14118737

ac(374)=0.11825825

ac(375)=0.10317144

ac(376)=0.08392241

ac(377)=0.06351954

ac(378)=0.04701854

ac(379)=0.02607570

ac(380)=0.00612922

ac(381)=-0.01090853

ac(382)=-0.03450940

ac(383)=-0.04898876

ac(384)=-0.03926621

ac(385)=-0.01904332

ac(386)=0.00154596

ac(387)=0.02102962

ac(388)=0.04026257

ac(389)=0.06230462

ac(390)=0.08158258

ac(391)=0.10218793

ac(392)=0.12321112

ac(393)=0.11375535

ac(394)=0.07273105

ac(395)=0.03080037

ac(396)=-0.01007609

ac(397)=-0.05136395

ac(398)=-0.09258431

ac(399)=-0.13515459

ac(400)=-0.14030030

ac(401)=-0.09506557

ac(402)=-0.04646250

ac(403)=-0.00435506

ac(404)=0.04341557

ac(405)=0.08796249

ac(406)=0.13636307

ac(407)=0.17954402

ac(408)=0.19444442

ac(409)=0.19544077

ac(410)=0.19704138

ac(411)=0.19950657

ac(412)=0.20282348

ac(413)=0.20244101

ac(414)=0.20844488

ac(415)=0.20146714

ac(416)=0.14091419

ac(417)=0.06138540

ac(418)=-0.00770411

ac(419)=-0.09293786

ac(420)=-0.14660308

ac(421)=-0.11154407

ac(422)=-0.04191139

ac(423)=0.01233558

ac(424)=0.06892559

ac(425)=0.13579096

ac(426)=0.18228883

ac(427)=0.16931044

ac(428)=0.10556270

ac(429)=0.04134571

ac(430)=-0.01331586

ac(431)=-0.07976341

ac(432)=-0.13345756

ac(433)=-0.12389251

ac(434)=-0.05997443

ac(435)=0.00207307

ac(436)=0.05719426

ac(437)=0.12773974

ac(438)=0.16605459

ac(439)=0.12641875

ac(440)=0.05727461

ac(441)=0.00607458

ac(442)=-0.04870592

ac(443)=-0.11544916

ac(444)=-0.15751803

ac(445)=-0.15890650

ac(446)=-0.14199734

ac(447)=-0.13028850

ac(448)=-0.11048986

ac(449)=-0.03620642

ac(450)=0.06333955

ac(451)=0.07738501

ac(452)=0.00523250

ac(453)=-0.07066761

ac(454)=-0.14915183

ac(455)=-0.19869019

ac(456)=-0.14114238

ac(457)=-0.03322377

ac(458)=0.05093326

ac(459)=0.14850903

ac(460)=0.23628511

ac(461)=0.21461909

ac(462)=0.07730787

ac(463)=-0.07347670

ac(464)=-0.12770760

ac(465)=-0.07153220

ac(466)=0.01057748

ac(467)=0.06796780

ac(468)=0.13055526

ac(469)=0.20931588

ac(470)=0.25367317

ac(471)=0.23495126

ac(472)=0.19182494

ac(473)=0.15627740

ac(474)=0.12165873

ac(475)=0.08412489

ac(476)=0.04547257

ac(477)=0.00920508

ac(478)=-0.02748667

ac(479)=-0.07277282

ac(480)=-0.09333961

ac(481)=-0.05521761

ac(482)=0.00742449

ac(483)=0.05976873

ac(484)=0.11199726

ac(485)=0.17241201

ac(486)=0.22250640

ac(487)=0.22121114

ac(488)=0.16146491

ac(489)=0.09371567

ac(490)=0.03784560

ac(491)=-0.02506327

ac(492)=-0.09071694

ac(493)=-0.14740339

ac(494)=-0.18046002

ac(495)=-0.18352303

ac(496)=-0.18492757

ac(497)=-0.19405872

ac(498)=-0.19680676

ac(499)=-0.20158605

ac(500)=-0.20909412

再贴后500个点

ac(501)=13.用beam54单元描述变截面梁的例子

作者：清华大学土木系，陆新征

/PREP7!*

MPTEMP,,,,,,,,

MPTEMP,1,0

MPDATA,EX,1,,200e3

MPDATA,PRXY,1,,.2

N,1,0,,,,,,

N,2,5,,,,,,

N,3,10,,,,,,

TYPE,1

MAT,1

REAL,1!*

e,1,2

REAL,2

e,3,2

/eshap,1

EPLOT14.变截面梁实例/prep7

et,1,beam188

keyopt,1,8,2

keyopt,1,9,2

mp,ex,1,3.5e10

mp,dens,1,2549

mp,prxy,1,0.1667

sectype,1,beam,rect!definecrosssectionatfirstendpoint

secdata,.0001,0.5

sectype,2,beam,rect!definecrosssectionatfarend

secdata,3,0.5

sectype,3,taper!newSectionIDfortaperedbeamanalysis

secdata,1,0.0,0.0!section1atlocation(0,0,0)

secdata,2,0.0,,-20.0!section2atlocation(0,-20,0)

k,1,-2.7,0,-27.85

k,2,-2.7,0,-19

k,3,-2.7,0,-13

lstr,1,2

lstr,2,3

lsel,s,,,1

latt,1,,1,,,,1

lesize,all,,,1

lmesh,all

lsel,all

lsel,s,,,2

latt,1,,1,,,,3

lesize,all,,,1

lmesh,all

/view,,1,2,3

/auto

/eshape,1

/rep15.拱桥浇筑过程分析-单元生死应用实例/prep7

et,1,BEAM3

mp,ex,1,2.1e11

mp,prxy,1,.3

mp,dens,1,7800

r,1,25.48e-4,250500E-8,1.5

K,1,,,,

K,2,10,5,,

K,3,20,0,,

LARC,1,3,