《ANSYS经典APDL程序实例.doc》由会员分享，可在线阅读，更多相关《ANSYS经典APDL程序实例.doc（169页珍藏版）》请在三一文库上搜索。
1、1ANSYS SOLID65环向布置钢筋的例子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.用beam 54单元描述变截面梁的例子3714.变截面梁实例3815.拱桥浇筑过程分析-单元生
2、死应用实例3816.简支梁实体与预应力钢筋分析实例3917.简单的二维焊接分析单元生死实例4018.隧道开挖（三维）的命令流4419.岩土接触分析实例5320.钢筋混凝土管的动力响应特性分析实例5721.隧道模拟开挖命令流（入门）6122.螺栓连接的模拟实现问题6323.道路的基层、垫层模量与应力之间的关系6723.滞回分析8024.模拟某楼层浇注8125.在面上施加移动的面力8327.在任意面施加任意方向任意变化的压力8528.预紧分析8629.几何非线性+塑性+接触+蠕变8730.埋设在地下的排水管道8932.幕墙企业玻璃简化计算9233.等截面杆单元生死应用实例10034.梁板建模联系1
3、0135.圆盘转动例子10136.简单的例子如何对结构的振动控制分析10237.模态分析结果的输出实例10338.火车过桥动态加载实例（部分）10439.悬索结构的找形和计算的例题11340.陶瓷杆撞击铝板的例子11541.求反作用力的APDL命令法11742.LS-DYNA实例（部分）11843.路面分层填筑对路基的影响11844.一个例子（含地震影响，求振兴与频率）12145.接触面上的压力总和12346.施加位置函数荷载12547.非线性分析考虑刚度退化12548.一个圆形水池的静力分析12649.ANSYS中混凝土模式预应力模拟的算例12750.悬臂梁受重力作用发生大变形求其固有频率1
4、2851.循环对称结构模态分析12952.三角平台受谐波载荷作用的结构响应13053.三角平台受一地震谱激励的应力分布和支反力13154.三角平台受时程载荷作用的应力分布和变形过程13255.经典层合板理论13356定易圆轨迹的例子13757模拟门式刚架施工单元生死13758钢筋混凝土整体式模型例子13859在荷载步之间改变材料属性例子13960含预应力的特征值屈曲计算14061振型叠加计算及工况组合例子14162柱子稳定分析算（预应力，特征值屈曲，初始缺陷）14363module MConcrete !混凝土模板14464.混凝土开裂实例14865螺栓网格划分14966自由液面的土石坝平面渗
5、流分析15067导出刚度矩阵15268某混凝土拱坝工程施工期及运行期温度场仿真分析15269移动温度荷载计算15670SHSD用于壳-实体装配实例An15771ansys显示隐式回弹分析实例16072工况组合的经典例子16816911ANSYS SOLID65环向布置钢筋的例子! 一个管道，环向配筋率为1%，纵向配筋率为0.5%，径向配筋率为0.1%! FINISH/CLEAR/PREP7 !* ! 单元属性ET,1,SOLID65!* KEYOPT,1,1,0！：单元主要选项设置KEYOPT,1,5,0！5：Number of the KEYOPT to be defined (KEYOPT
6、(KNUM).KEYOPT,1,6,0！0：Value of this KEYOPTKEYOPT,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, ！Activates a data table for nonline
7、ar material properties or special element input.TBTEMP,0！Defines a temperature for the data table.TBPT,0.0005,15！Defines a point on a nonlinear data curve在应力应变区线上定义一点TBPT,0.001,21TBPT,0.0015,24 TBPT,0.002,27TBPT,0.003,24!破坏准则 TB,CONC,1,1,9, TBTEMP,0TBDATA,.5,.9,3,30,!抗拉和抗压是10 倍的关系TBDATA,1, MPTEMP, M
8、PTEMP,1,0 !钢材基本属性MPDATA,EX,2,200E3MPDATA,PRXY,2,.27 !屈服准则TB,BISO,2,1,2, TBTEMP,0TBDATA,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,-
9、90.000000, VSBW, 3 WPCSYS,-1,0 !Defines the working plane location based on a coordinate system.KWPAVE, 1 WPRO,-90.000000 VSBW,ALLESIZE,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、：红色代表最大配筋方向，绿色代表其次，蓝色表示最小配筋方向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,
11、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 MSHK
12、EY,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
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约束压 !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
14、, 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,1600
15、E-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 N
16、SUBST,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）- KIN
17、H滞回 !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 TBP
18、T,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,
19、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 /A
20、XLAB,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
21、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 /REP
22、LOT /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 /CLE
23、AR /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,1800
24、E-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
25、 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 E
26、T,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.
27、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
28、,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
29、,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
30、, 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 N
31、SEL,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 T
32、IME,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
33、,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,MK
34、IN,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 /SO
35、LU 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
36、 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
37、,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
38、,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 NSU
39、BST,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.02ET,1,BEAM188KEYOPT,1,7,1KEYOPT,1,8,1ET,2,MASS21,4MP,EX,1,210E9 MP,NUXY,1,0.3 MP
40、,DENS,1,7800 MPDATA,DAMP,1,DAMPRATIOSECTYPE,3,BEAM,ISECDATA,0.2,0.2,0.5,0.01,0.01,0.008R,1,1.6E2,1.6E2R,2,1.2E2,1.2E2K,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,10000L,1,2L,2,3L,3,4L,4,5LSEL,1,4LATT,1,0,1,6,3 ALLSLESIZE,ALL, , ,3 , , , , ,1 LMESH,ALL NUMMRG,ALLNUMCMP,ALLTYPE,2REAL,1E,
41、2TYPE,2REAL,2E,8E,14E,20*SET,NT,1001*SET,DT,0.02 *DIM,AC,NT *VREAD,AC(1),RECORD,TXT(F8.3)/SOLU D,1,ALL,ANTYPE,0SOLVEFINISH/SOLU !模态分析ANTYPE,2MODOPT,SUBSP,8 MXPAND,8, , ,1 SOLVEFINI!得到自振频率1*GET,FREQ1,MODE,1,FREQ/ESHAPE,1.0 /SOLUANTYPE,SPECTR !谱分析SPOPT,SPRSMDAMP,1,DAMPRATIO, , , , , , SED,1SVTYPE,2,1F
42、REQ,0.167,0.25,0.333,0.44444,0.5,0.667,1,1.25,1.667SV,0.154350625,0.191590625,0.210210625,0.224175625,0.250716714,0.329514922,0.484352764,0.598723486,0.786897371FREQ,2.222,10,100000SV,1.034212766,1.034212766,0.3528!按7度多遇地震，地震影响系数0.08 第一组III场地TG0.45SRSS,0.01,DISPSOLVEFINISH/POST1/INP,MCOMNSORT,U,SUM *GET,DMAX_SPE,SORT, ,MAX !反应谱最大变形!开始地震计算FINISH/CONFIG,NRES,20000/SOLU ANTYPE,TRANS TRNOPT,FULL ALPHAD,2*DAMPRATIO*FREQ1*2*3.1415926BETAD,2*DAMPRATIO/(FREQ1*2*3.1415926)*DO,I,1,100