ansys中混凝土结构的计算

(1)转贴一个例题，来自中国有限元联盟论坛 -运行没有问题/title, fixed - fixed concrete beam example /prep7 et,1,65 mp,ex,1,3e7 ! steel rebar (units are pounds, inches) MP,PRXY, 1,0.3mp,ex,2,1e6, ! concrete MP,PRXY, 2,0.3mp,dens,2,.00025 tb,concr,2 tbdata,1,.3,.5,200,4000 ! shear coeffs, tensile and compress strength r,1,1,.03,0,0 ! mat 1 (steel), 3 percent reinforcement in x dir r,2,1,.01,0,0 ! mat 1 (steel), 1 percent r,3,1,.04,0,0 ! mat 1 (steel), 4 percent block,100,5,5 block,100,5,10,5 block,100,10,15,5 !vovlap,all NUMMRG,KP, , , ,LOW numcmp,volu esize,5 mat,2 ! concrete material real,1 ! rebar vmesh,1 real,2 vmesh,2 real,3 vmesh,3 nsel,s,loc,x d,all,all nsel,s,loc,x,100 d,all,all nsel,all fini /solu nsel,s,loc,y,15 sf,all,pres,100 nsel,all OUTRES,ALL,ALL, nsub,10 solve fini /post1 /DEVICE,VECTOR,1 set,last plcrack fini (2)! 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 , !材料属性!混凝土基本材料属性MPTEMP, MPTEMP,1,0 MPDATA,EX,1,30E2MPDATA,PRXY,1,.2 !屈服准则 TB,MISO,1,1,5, !多线性等向强化，材料号为1，数据的温度数为1，对某一给定温度数据的点数为5。TBTEMP,0 !为材料表定义温度值。TBPT,0.0005,15TBPT,0.001,21TBPT,0.0015,24 TBPT,0.002,27TBPT,0.003,24!破坏准则 TB,CONC,1,1,9, !材料号为1，数据的温度数为1，对某一给定温度数据的点数为9。TBTEMP,0 !为材料表定义温度值。TBDATA,.5,.9,3,30, !给当前数据表定义数据TBDATA,1, MPTEMP, MPTEMP,1,0 !钢材基本属性MPDATA,EX,2,200E3MPDATA,PRXY,2,.27 !屈服准则TB,BISO,2,1,2, !双线性等向强化，材料号为2，数据的温度数为1，对某一给定温度数据的点数为2。TBTEMP,0 !为材料表定义温度值为0度。TBDATA,310,2E3, !给当前数据表定义数据。! 管道内径和外径CYL4,0,0,3000, , , ,10000 !生成内径为2000，外径为3000，长度为10000的管道CYL4,0,0,2000, , , ,10000 VSBV, 1, 2 !体相减，在体1中减去体2! 定义局部柱坐标CSWPLA,11,1,1,1, !通过当前定义的工作平面的原点为中心建立局部坐标系，第一个1表示为柱坐标。后 !两个1表示比率，一般不用KWPAVE, 11 !移动工作平面到关键点11WPRO,-90., !旋转工作平面VSBW, 3 !由工作平面剖分体3WPCSYS,-1,0 !通过当前坐标系的R-0来定义工作平面，-1表示图形窗口编号，0为坐标系的编号。KWPAVE, 1 !移动工作平面到关键点1WPRO,-90. !旋转工作平面VSBW,ALL !工作平面剖分体ESIZE,500,0, !设置单元尺寸大小! 注意：设定单元局部坐标VATT, 1, 1, 1, 11 !给体分配属性，材料号为1，实常数为1，TYPE号为1，坐标系号为11。!* VSWEEP,ALL !划分网格，体扫掠划分/DEVICE,VECTOR,1 !绘制等值线/ESHAPE,1.0 /REPLO ! 注意：红色代表最大配筋方向，绿色代表其次，蓝色表示最小配筋方向(3)/clear beam189单元的定性命令，能运行的通/filname,beam189test,1/PREP7 ET,1,BEAM189 !定义单元类型及其KEYOPT选项KEYOPT,1,7,2KEYOPT,1,8,3KEYOPT,1,9,3SECTYPE, 1, BEAM, RECT, , 0 !定义第一种截面SECOFFSET, CENT SECDATA,300,500,0,0,0,0,0,0,0,0 SECTYPE, 2, BEAM, RECT, , 0 !定义第二种截面 SECOFFSET, CENT SECDATA,500,500,0,0,0,0,0,0,0,0 MPTEMP, !定义第一种材料弹性模量，泊松比，密度MPTEMP,1,0 MPDATA,EX,1,2.06e5 MPDATA,PRXY,1,0.3 MPTEMP, MPTEMP,1,0 MPDATA,DENS,1,7.85e-9K,1,0,0,0, !定义关键点KEYPOINTK,2,6000,0,0, K,3,6000,4000,0,K,4,0,4000,0, K,5,0,0,-5000, K,6,6000,0,-5000, K,7,6000,4000,-5000,K,8,0,4000,-5000, LSTR, 1, 2 !通过KEYPOINT定义线LINE，梁LSTR, 2, 3 LSTR, 3, 4 LSTR, 4, 1 LSTR, 1, 5 !通过KEYPOINT定义线LINE，柱LSTR, 2, 6 LSTR, 3, 7 LSTR, 4, 8 LSEL,S,1,2,1,!给LINE赋属性材料，截面，单元类型LATT,1, ,1, ,6, ,1LSEL,S,3,4,1,LATT,1, ,1, ,8, ,1LSEL,S,5,LATT,1, ,1, ,8, ,2LSEL,S,6,LATT,1, ,1, ,7, ,2LSEL,S,7,LATT,1, ,1, ,6, ,2LSEL,S,8,LATT,1, ,1, ,5, ,2ALLSEL,ALLLESIZE,ALL,1000, , , , , , ,1 !指定划分单元长度LMESH,ALL!划分单元/soluNSEL,S,LOC,Z,-5000 !按位置选择需要施加约束的节点D,ALL, , , , , ,ALL, , , , , !柱底施加固定约束 ACEL,0,0,9810, !施加重力荷载ESEL,S,SEC,1 !选择全部梁SFBEAM,ALL,1,PRES,-500,-500, , , , ,!对全部梁施加压力PRESSUREALLSEL,ALLSOLVE!求解/post1!绘制位移结果PLNSOL,U,Z,0,1 /ESHAPE,1!绘制应力结果(4)finish 带有三个圆孔的平面支座的分析，可以运行的通/clear,start/filname,ex2_2/title,2d plane stress bracket/prep7blc4,0,0,80,100cyl4,80,50,50cyl4,0,20,20cyl4,0,80,20blc4,-20,20,20,60aadd,allcyl4,80,50,30cyl4,0,20,10cyl4,0,80,10asba,6,1asba,4,3asba,1,2et,1,plane82keyopt,1,3,3r,1,20mp,ex,1,mp,prxy,1,0.3aesize,all,7amesh,allfinish/soluantype,0dl,7,all,0dl,8,all,0dl,9,all,0dl,10,all,0dl,11,all,0dl,12,all,0dl,13,all,0dl,14,all,0fk,9,fy,-1000solvefinish/post1plnsol,s,eqvancntr,10,0.15finish(5)/title,buckling of a bar 杆件的屈曲分析，可以运行的通/prep7et,1,plane42,3r,1,0.01mp,ex,1,2.1e11mp,prxy,1,0.3n,1n,21,1fillngen,2,21,1,21,0.01e,22,23,2,1egen,20,1,1finish/solud,1,all,22,21f,21,fy,-0.5,42,21pstrss,onsolvefinish/soluantype,bucklebucopt,subsp,1mxpand,1solvefinishFinish在接头固定处定义预紧，能运行的通/clear/prep7 /title,Sample application of PSMESH!标题，PSMESH运用示例et,1,92mp,ex,1,1e7 MP,PRXY, 1,0.3 mp,alpx,1,1.3e-5mp,ex,2,3e7 MP,PRXY, 2,0.3 mp,alpx,2,8.4e-6tref,70 /foc,-.09,.34,.42 /dist,.99/ang,-55.8 /view,.39,-.87,.31/pnum,volu,1 /num,1cylind,0.5, -0.25,0, 0,180cylind,0.5, 1,1.25, 0,180cylind,0.25, 0,1, 0,180wpoff,.05 wpstyle,0cylind,0.35,1, 0,0.75, 0,180wpoff,-.1 wpstyle,0cylind,0.35,1, 0.75,1, 0,180vglue,all numc,all vplotmat,1 smrt,off vmesh,4,5mat,2 vmesh,1,3/pnum,mat,1 eplotpsmesh,example,volu,1,0,z,0.5,n CM,lines,LINE /dist,1.1 cmplot/solueqslve,pcg,1e-8asel,s,loc,y da,all,symm asel,alldk,1,ux dk,12,ux dk,1,uzsload,1,9,tiny,force,100,1,2/title,Sample application of PSMESH - preload only solve !Finally, we construct the actual solution of interest. We want to!know what happens to the preload in the bolt, and the stress field around!it, when the assembly temperature rises to 150 F.!Both the preload and the stresses increase because, for a uniform!temperature rise, there is greater thermal expansion in the aluminum plates!than in the steel bolt. Any method for applying preload that did not!allow the load to change would be unable to predict this result. !/post1plnsol,s,z/soluantype,restart tunif,150 /title,Sample application of PSMESH - uniform 150 deg solve /post1 plnsol,s,z(6)发信人: windoe (风娃娃为了论文，冲啊！), 信区: Structure 标 题: 大牛yulee的混凝土非线性计算实例（1） Ansys 发信站: BBS 水木清华站 (Thu Dec 4 19:19:34 2003), 转信 ，运行没有问题!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 (7)发信人: windoe (风娃娃为了论文，冲啊！), 信区: Structure 标 题: 大牛yulee的混凝土非线性计算实例（2） Ansys 发信站: BBS 水木清华站 (Thu Dec 4 19:20:18 2003), 转信 ，运行没有问题!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 (8)发信人: windoe (风娃娃为了论文，冲啊！), 信区: Structure 标 题: 大牛yulee的混凝土非线性计算实例（3） Ansys 发信站: BBS 水木清华站 (Thu Dec 4 19:22:00 2003), 转信 ，运行有问题!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 (9)发信人: windoe (风娃娃为了论文，冲啊！), 信区: Structure 标 题: 大牛yulee的混凝土非线性计算实例（4） Ansys 发信站: BBS 水木清华站 (Thu Dec 4 19:24:06 2003), 转信 ，运行没有问题!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 (10)发信人: windoe (风娃娃为了论文，冲啊！), 信区: Structure 标 题: 大牛yulee的混凝土非线性计算实例（5） Ansys 发信站: BBS 水木清华站 (Thu Dec 4 19:26:10 2003), 转信 ，运行有问题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 (11)发信人: windoe (风娃娃为了论文，冲啊！), 信区: Structure 标 题: 大牛yulee的混凝土非线性计算实例（6） Ansys 发信站: BBS 水木清华站 (Thu Dec 4 19:27:21 2003), 转信 ，运行有问题!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 (12)发信人: windoe (风娃娃为了论文，冲啊！), 信区: Structure 标 题: 大牛yulee的混凝土非线性计算实例（7） Ansys 发信站: BBS 水木清华站 (Thu Dec 4 19:28:29 2003), 转信 ，运行没有问题!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 (13)发信人: windoe (风娃娃为了论文，冲啊！), 信区: Structure 标 题: 大牛yulee的混凝土非线性计算实例（8） Ansys 发信站: BBS 水木清华站 (Thu Dec 4 19:31:21 2003), 转信 ，运行有问题!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