大神来解答
浏览：65033
我在运行ANSYS土木工程应用实例（第三版）中的混凝土开裂命令流时，发现居然运行不起来求大神指导！命令流如下：!采用均匀网格划分，同时取消两种垫片,DMX=6.415,最大第一主应力263.629，28个子部计算完毕
/clear,nostart
/config,nres,5000
/prep7
!定义单元类型
et,1,solid65 !混凝土单元
et,2,pipe20
!钢筋单元
et,3,plane42
!用于拉伸成体单元之用
!定义实常数
r,1,18,8.99
!纵向受拉钢筋实常数
r,2,8,3.99
!纵向受压钢筋，箍筋实常数
r,3
!混凝土
!定义材料属性
mp,ex,1,24000
!混凝土材料属性
mp,prxy,1,0.2
tb,conc,1,1,9
tbdata,,0.4,1,3,-1
mp,ex,2,2e5
!纵向受拉钢筋
mp,prxy,2,0.3
tb,bkin,2,1,2,1
tbdata,,350
mp,ex,3,2e5
!横向箍筋，受压钢筋材料属性
mp,prxy,3,0.25
tb,bkin,3,1,2,1
tbdata,,200
/pnum,node,1
/pnum,elem,1
!!产生所有的节点**********!!!
n,1
n,9,150
fill,1,9
ngen,11,9,1,9,1,,30
ngen,19,1000,1,99,1,,,-75
/view,1,1,1,1
!!!!!!!箍筋，受压钢筋*******!!!!
type,2
real,2
mat,3
!水平箍筋
*do,ii,11,16,1
e,ii,ii+1
*enddo
*do,ii,83,88,1
e,ii,ii+1
*enddo
!竖直箍筋
*do,ii,11,74,9
e,ii,ii+9
*enddo
*do,ii,17,80,9
e,ii,ii+9
*enddo
!产生整个模型的箍筋
egen,19,1000,1,28,1
!此时最大的单元编号为532
!纵向受压钢筋
*do,ii,83,17083,1000
e,ii,ii+1000
*enddo
*do,ii,89,17089,1000
e,ii,ii+1000
*enddo
!!!!!纵向受拉钢筋************!!!!!!!!!!
type,2
real,1
mat,2
*do,ii,11,17011,1000
e,ii,ii+1000
*enddo
*do,ii,17,17017,1000
e,ii,ii+1000
*enddo
!此时最大单元编号为604
/view,1,1,1,1
/pnum,elem,0
/pnum,node,0
/eshape,1
eplot
!!!**********生成面单元，以便拉伸成体单元!!!!!!
k,1
k,2,150,
k,3,150,300
k,4,,300
a,1,2,3,4
lsel,s,loc,y,0
lsel,a,loc,y,300
lesize,all,,,8
lsel,all
lsel,s,loc,x,0
lsel,a,loc,x,150
lesize,all,,,10
type,3
amesh,all
type,1
!拉伸成混凝土单元
real,3
mat,1
extopt,esize,18
extopt,aclear,1
vext,all,,,,,-1350
/pnum,mat,1
/pnum,node,0
/pnum,elem,0
/view,1,1,1,1
eplot
allsel
nummrg,all
numcmp,all
eplot
/solu
nsel,s,loc,y,-0
!右边界约束
nsel,r,loc,z,-1275
d,all,uy
d,all,ux
asel,s,loc,z,0
!对称面约束
da,all,symm
nsel,all
!施加外部载荷
nsel,s,loc,y,300
nsel,r,loc,z,-600
d,all,uy,-5
allsel
gplot
nlgeom,on
NSUBST,200
OUTRES,ALL,ALL
NEQIT,50
pred,on
CNVTOL,F, ,0.05,2,0.5
allsel
SOLVE
FINISH
/post1
allsel
/DEVICE,VECTOR,1
PLCRACK,0,1
!显示第一开裂位置
PLCRACK,0,2
!显示第二开裂位置
!/DEVICE,VECTOR,0
!PLNSOL,U,Y,0,1
!ESEL,S,TYPE,,2
!PLNSOL,S,1,0,1
/clear,nostart
/config,nres,5000
/prep7
!定义单元类型
et,1,solid65 !混凝土单元
et,2,pipe20
!钢筋单元
et,3,plane42
!用于拉伸成体单元之用
!定义实常数
r,1,18,8.99
!纵向受拉钢筋实常数
r,2,8,3.99
!纵向受压钢筋，箍筋实常数
r,3
!混凝土
!定义材料属性
mp,ex,1,24000
!混凝土材料属性
mp,prxy,1,0.2
tb,conc,1,1,9
tbdata,,0.4,1,3,-1
mp,ex,2,2e5
!纵向受拉钢筋
mp,prxy,2,0.3
tb,bkin,2,1,2,1
tbdata,,350
mp,ex,3,2e5
!横向箍筋，受压钢筋材料属性
mp,prxy,3,0.25
tb,bkin,3,1,2,1
tbdata,,200
/pnum,node,1
/pnum,elem,1
!!产生所有的节点**********!!!
n,1
n,9,150
fill,1,9
ngen,11,9,1,9,1,,30
ngen,19,1000,1,99,1,,,-75
/view,1,1,1,1
!!!!!!!箍筋，受压钢筋*******!!!!
type,2
real,2
mat,3
!水平箍筋
*do,ii,11,16,1
e,ii,ii+1
*enddo
*do,ii,83,88,1
e,ii,ii+1
*enddo
!竖直箍筋
*do,ii,11,74,9
e,ii,ii+9
*enddo
*do,ii,17,80,9
e,ii,ii+9
*enddo
!产生整个模型的箍筋
egen,19,1000,1,28,1
!此时最大的单元编号为532
!纵向受压钢筋
*do,ii,83,17083,1000
e,ii,ii+1000
*enddo
*do,ii,89,17089,1000
e,ii,ii+1000
*enddo
!!!!!纵向受拉钢筋************!!!!!!!!!!
type,2
real,1
mat,2
*do,ii,11,17011,1000
e,ii,ii+1000
*enddo
*do,ii,17,17017,1000
e,ii,ii+1000
*enddo
!此时最大单元编号为604
/view,1,1,1,1
/pnum,elem,0
/pnum,node,0
/eshape,1
eplot
!!!**********生成面单元，以便拉伸成体单元!!!!!!
k,1
k,2,150,
k,3,150,300
k,4,,300
a,1,2,3,4
lsel,s,loc,y,0
lsel,a,loc,y,300
lesize,all,,,8
lsel,all
lsel,s,loc,x,0
lsel,a,loc,x,150
lesize,all,,,10
type,3
amesh,all
type,1
!拉伸成混凝土单元
real,3
mat,1
extopt,esize,18
extopt,aclear,1
vext,all,,,,,-1350
/pnum,mat,1
/pnum,node,0
/pnum,elem,0
/view,1,1,1,1
eplot
allsel
nummrg,all
numcmp,all
eplot
/solu
nsel,s,loc,y,-0
!右边界约束
nsel,r,loc,z,-1275
d,all,uy
d,all,ux
asel,s,loc,z,0
!对称面约束
da,all,symm
nsel,all
!施加外部载荷
nsel,s,loc,y,300
nsel,r,loc,z,-600
d,all,uy,-5
allsel
gplot
nlgeom,on
NSUBST,200
OUTRES,ALL,ALL
NEQIT,50
pred,on
CNVTOL,F, ,0.05,2,0.5
allsel
SOLVE
FINISH
/post1
allsel
/DEVICE,VECTOR,1
PLCRACK,0,1
!显示第一开裂位置
PLCRACK,0,2
!显示第二开裂位置
!/DEVICE,VECTOR,0
!PLNSOL,U,Y,0,1
!ESEL,S,TYPE,,2
!PLNSOL,S,1,0,1
