发布于:2009-12-23 18:08:23
来自:建筑结构/ANSYS
[复制转发]
我在做钢筋混凝土梁计算时,出现我定义的569号单元初始应变大于1,有点不明白!求助高手!谢谢!命令流如下:/clear,nostart
/config,nres,5000
/prep7
!定义单元类型
et,1,solid65 !混凝土单元
et,2,link8 !钢筋单元
et,3,plane42 !用于拉伸成体单元之用
!定义实常数
r,1,18,8.99 !纵向受拉钢筋实常数
r,2,0.03 !纵向受压钢筋,箍筋实常数
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
全部回复(1 )
只看楼主 我来说两句 抢板凳