9 famous examples of concrete in ANSYS 
  by Yulee@smth
--------------------------------------------------------------------------------
  标 题: 大牛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 
发信人: 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 
发信人: 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 
发信人: 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 
发信人: 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 
发信人: 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 
发信人: 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 
发信人: 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,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 
发信人: windoe (风娃娃◎为了论文,冲啊!!!!), 信区: Structure 
  标 题: 大牛yulee的混凝土非线性计算实例(9) Ansys 
  发信站: BBS 水木清华站 (Thu Dec 4 19:35:38 2003), 转信 
!梁平面应力 
  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