我的照着书上的例子做了一个模型，怎么会出现这样的情况?
书上的是矩形梁，运算没问题，我改成了T形梁，就出错了。原矩形梁：钢筋混凝土梁
/clear,nostart
/config,nres,5000
/prep7
!定义单元类型
et,1,solid65 !混凝土单元
et,2,pipe59 !钢筋单元
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
能正常运算。
我做的T形梁：
fini
/clear,nostart
/config,nres,5000
/prep7
!定义单元类型
et,1,solid65 !混凝土单元
et,2,pipe59 !钢筋单元
et,3,plane42 !用于拉伸成体单元之用
!定义实常数
r,1,16,7.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,11,500
fill,1,11
ngen,6,11,1,11,1,,-20
ngen,13,8,59,63,1,,-20
ngen,12,1000,1,159,1,,,-100
/view,1,1,1,1
!!!!!!!箍筋，受压钢筋*******!!!!
type,2
real,2
mat,3
!!腹板箍筋
*do,ii,16,49,11
e,ii,ii+11
*enddo
*do,ii,60,140,8
e,ii,ii+8
*enddo
*do,ii,16,17,1
e,ii,ii+1
*enddo
*do,ii,18,51,11
e,ii,ii+11
*enddo
*do,ii,62,142,8
e,ii,ii+8
*enddo
*do,ii,148,149,1
e,ii,ii+1
*enddo
!产生整个模型的箍筋
egen,12,1000,1,54,1
!纵向受压钢筋
*do,ii,16,10016,1000
e,ii,ii+1000
*enddo
*do,ii,18,10018,1000
e,ii,ii+1000
*enddo
!!!!!纵向受拉钢筋************!!!!!!!!!!
type,2
real,1
mat,2
*do,ii,148,10148,1000
e,ii,ii+1000
*enddo
*do,ii,150,10150,1000
e,ii,ii+1000
*enddo
/view,1,1,1,1
/pnum,elem,0
/pnum,node,0
/eshape,1
eplot
!!!**********生成面单元，以便拉伸成体单元!!!!!!
k,1
k,2,0,-100
k,3,150,-100
k,4,150,-340
k,5,350,-340
k,6,350,-100
k,7,500,-100
k,8,500
a,1,2,3,4,5,6,7,8
lsel,s,loc,y,0
lesize,all,,,10
lsel,all
lsel,s,loc,y,-100
lesize,all,,,3
lsel,all
lsel,s,loc,y,-340
lesize,all,,,4
lsel,all
lsel,s,loc,x,0
lsel,a,loc,x,500
lesize,all,,,2
lsel,all
lsel,s,loc,x,150
lsel,a,loc,x,350
lesize,all,,,6
lsel,all
type,3
amesh,all
type,1 !拉伸成混凝土单元
real,3
mat,1
extopt,esize,25
extopt,aclear,1
vext,all,,,,,-1200
/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,-340 !左边界约束
nsel,r,loc,z,-1000
d,all,uy
d,all,ux
nsel,all
asel,s,loc,z,0 !对称面约束
da,all,symm
nsel,all
nsel,s,loc,y,0 !施加外部载荷
nsel,r,loc,z,-400
d,all,dy,-5
allsel
gplot
nlgeom,on
NSUBST,200
OUTRES,ALL,ALL
NEQIT,50
pred,on
CNVTOL,F,,0.05,2,0.5
allsel
SOLVE
FINISH
菜鸟求指导。
求解前检查一下，节点 1382， 1383， 1384 等是否在当前实体集中？
或者，求解前，执行命令：
nsel, inve
esel, inve
然后对节点和单元列表，看是否有节点和/或单元。 htbbzzg 发表于 2014-4-4 08:53
求解前检查一下，节点 1382， 1383， 1384 等是否在当前实体集中？
或者，求解前，执行命令：
nsel,...
反向选择之后什么都没有了。 htbbzzg 发表于 2014-4-4 08:53
求解前检查一下，节点 1382， 1383， 1384 等是否在当前实体集中？
或者，求解前，执行命令：
nsel,...
我找到了，我选的作用约束的位置有问题，不在结点上。谢了~！
页:
[1]