哪位大神帮忙看下，这个CFPR钢管混凝土命令流得出的位移荷载曲线为折线
浏览：4878
finish
/clear
/PREP7
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!定义材料属性
et,1,181
r,1,0.504
mp,ex,1,2.3e5
MP,PRXY,1,0.3
tb,bkin,1,1,2,1
tbdata,,3550
et,2,65
mp,ex,2,3.0e4 !混凝土弹模
mp,prxy,2,0.2 !混凝土泊松比
tb,conc,2 !定义非线性材料特性
tbdata,,0.5,0.9,2.55,-1,, !张开系数，闭合，单轴抗拉强度，单轴抗压
tb,miso,2,1,12 !多线性等向强化
tbpt,,0.0001,3.0 !应变应力关系
tbpt,,0.0003,7.076
tbpt,,0.0005,11.156
tbpt,,0.0007,14.726
tbpt,,0.0009,17.786
tbpt,,0.0011,20.336
tbpt,,0.0013,22.376
tbpt,,0.0015,23.906
tbpt,,0.0017,24.926
tbpt,,0.002,25.5
tbpt,,0.0022,25.5
tbpt,,0.0023,25.5
et,3,45
MP,EX,3,2.06E5
MP,PRXY,3,0.3
TB,bkin,3,1,2, 1 !双线性随动强化模型
TBDATA,,345
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!设置参数区
R1=60
R2=63
bb=4/3
zzz=600
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!建立模型
csys,1
k,1, R1/bb,0
k,2, R1/bb,90
k,3, R1/bb,180
k,4, R1/bb,270
k,5, R1,0
k,6, R1,90
k,7, R1,180
k,8, R1,270
k,9, R2,0
k,10, R2,90
k,11, R2,180
k,12, R2,270
csys,0
l,1,2
l,2,3
l,3,4
l,4,1
csys,1
l,5,6
l,6,7
l,7,8
l,8,5
l,9,10
l,10,11
l,11,12
l,12,9
csys,0
l,1,5
l,2,6
l,3,7
l,4,8
l,5,9
l,6,10
l,7,11
l,8,12
al,1,2,3,4
al,1,13,5,14
al,2,14,6,15
al,3,15,7,16
al,4,16,8,13
al,17,9,18,5
al,18,10,19,6
al,19,11,20,7
al,20,12,17,8
!
alls
csys,1
lsel,s,loc,y,45
lplot
lesize,all,,,4,
!
alls
csys,1
lsel,s,loc,y,135
lplot
lesize,all,,,4,
!
alls
csys,1
lsel,s,loc,y,225
lplot
lesize,all,,,4,
!
alls
csys,1
lsel,s,loc,y,315
lplot
lesize,all,,,4,
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
LSEL,S,LOC,Y,0
LSEL,A,LOC,Y,90
LSEL,A,LOC,Y,180
LSEL,A,LOC,Y,270
LPLOT
LSEL,U,LOC,X,R1,R2
LPLOT
lesize,all,,,2,
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
LSEL,S,LOC,Y,0
LSEL,A,LOC,Y,90
LSEL,A,LOC,Y,180
LSEL,A,LOC,Y,270
LPLOT
LSEL,U,LOC,X,0,R1
LPLOT
lesize,all,,,1
!
alls
csys,1
asel,s,loc,x,0,r1-5
aplot
type, 1
amesh ,all
alls
asel,u,loc,x,0,r1-5
aplot
type,1
amesh,all
alls
type,2
mat,2
real,2
asel,s,loc,x,0,r1-5
extopt,esize,zzz/20
EXTOPT,ACLEAR,1
EXTOPT,ATTR,0,0,0
csys,0
VEXT,all, , ,0,0,zzz,,,,
eplot
alls
csys,1
type,3
mat,3
real,3
asel,u,loc,x,0,r1-5
asel,u,loc,z,1,zzz
extopt,esize,zzz/20
EXTOPT,ACLEAR,1
EXTOPT,ATTR,0,0,0
csys,0
VEXT,all, , ,0,0,zzz,,,,
alls
eplot
!接触单元设置
et,4,170
et,5,173
keyopt,5,9,0
keyopt,5,12,2
r,5
real,5
mp,mu,2,0.2
csys,1
vsel,s,type,,3
asel,r,loc,x,r1
aplot
cm,target,area
type,4
!设置目标单元
nsla,s,1
esurf,all
esel,s,type,,4
eplot
/psyms,esys,1
allsel
vsel,s,type,,2
asel,r,loc,x,r1
cm,contact,area
type,5
nsla,s,1
esurf,all
esel,s,type,,5
eplot
/psyms,esys,1
allsel
gplot
csys,1
type,1
mat,1
real,1
*do,i,0,15
*do,z,0,zzz/20-1
node1=node(r2,i*90/4,z*20)
node2=node(r2,i*90/4,(z+1)*20)
node4=node(r2,(i+1)*90/4,z*20)
node3=node(r2,(i+1)*90/4,(z+1)*20)
e,node1,node2,node3,node4
*enddo
*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!! 求解
NSEL,S,LOC,Z,0
NPLOT
D,all, , , , , ,ALL, , , , ,
ALLSEL,ALL
EPLOT
!————————
/solu
antype,static
outres,basic,all !输出选项
outpr,basic,all !窗口输出，这个没什么影响
autots,on !自动计算时间步长
nlgeom,on !打开是大变形，非线性
outres,all,all !输出所有子步结果
time,1
allsel
NSEL,S,LOC,Z,zzz
NPLOT
D,all, ,-25, , , ,UZ, , , , ,
nsubst,200 ,50 ,200 !荷载子步
neqit,20
eqslv,-1
ncnv,1 !终止分析选项，终止分析且退出程序
lnsrch,on !线性搜索，使计算稳定
CNVTOL,F, ,0.05,2, ,
allsel
solve
/POST26
numvar,10
csys,1
nsel,s,loc,z,zzz
*get,ncount,node,,count !获取被选择节点的数量定义到变量ncount中
nnow=0
nsol,6,ndnext(nnow),u,z,Zdisplancement
add,6,6,,,,,,-1
*do,i,1,ncount
nnow=ndnext(nnow)
rforce,8,nnow,f,z,zforce
add,7,7,8,,,,,,1
*enddo
xvar,6
/axlab,x,deflection(mm)
/axlab,y,totle load(N)
prod,7,,,,,,,-1
plvar,7
/clear
/PREP7
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!定义材料属性
et,1,181
r,1,0.504
mp,ex,1,2.3e5
MP,PRXY,1,0.3
tb,bkin,1,1,2,1
tbdata,,3550
et,2,65
mp,ex,2,3.0e4 !混凝土弹模
mp,prxy,2,0.2 !混凝土泊松比
tb,conc,2 !定义非线性材料特性
tbdata,,0.5,0.9,2.55,-1,, !张开系数，闭合，单轴抗拉强度，单轴抗压
tb,miso,2,1,12 !多线性等向强化
tbpt,,0.0001,3.0 !应变应力关系
tbpt,,0.0003,7.076
tbpt,,0.0005,11.156
tbpt,,0.0007,14.726
tbpt,,0.0009,17.786
tbpt,,0.0011,20.336
tbpt,,0.0013,22.376
tbpt,,0.0015,23.906
tbpt,,0.0017,24.926
tbpt,,0.002,25.5
tbpt,,0.0022,25.5
tbpt,,0.0023,25.5
et,3,45
MP,EX,3,2.06E5
MP,PRXY,3,0.3
TB,bkin,3,1,2, 1 !双线性随动强化模型
TBDATA,,345
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!设置参数区
R1=60
R2=63
bb=4/3
zzz=600
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!建立模型
csys,1
k,1, R1/bb,0
k,2, R1/bb,90
k,3, R1/bb,180
k,4, R1/bb,270
k,5, R1,0
k,6, R1,90
k,7, R1,180
k,8, R1,270
k,9, R2,0
k,10, R2,90
k,11, R2,180
k,12, R2,270
csys,0
l,1,2
l,2,3
l,3,4
l,4,1
csys,1
l,5,6
l,6,7
l,7,8
l,8,5
l,9,10
l,10,11
l,11,12
l,12,9
csys,0
l,1,5
l,2,6
l,3,7
l,4,8
l,5,9
l,6,10
l,7,11
l,8,12
al,1,2,3,4
al,1,13,5,14
al,2,14,6,15
al,3,15,7,16
al,4,16,8,13
al,17,9,18,5
al,18,10,19,6
al,19,11,20,7
al,20,12,17,8
!
alls
csys,1
lsel,s,loc,y,45
lplot
lesize,all,,,4,
!
alls
csys,1
lsel,s,loc,y,135
lplot
lesize,all,,,4,
!
alls
csys,1
lsel,s,loc,y,225
lplot
lesize,all,,,4,
!
alls
csys,1
lsel,s,loc,y,315
lplot
lesize,all,,,4,
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
LSEL,S,LOC,Y,0
LSEL,A,LOC,Y,90
LSEL,A,LOC,Y,180
LSEL,A,LOC,Y,270
LPLOT
LSEL,U,LOC,X,R1,R2
LPLOT
lesize,all,,,2,
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
LSEL,S,LOC,Y,0
LSEL,A,LOC,Y,90
LSEL,A,LOC,Y,180
LSEL,A,LOC,Y,270
LPLOT
LSEL,U,LOC,X,0,R1
LPLOT
lesize,all,,,1
!
alls
csys,1
asel,s,loc,x,0,r1-5
aplot
type, 1
amesh ,all
alls
asel,u,loc,x,0,r1-5
aplot
type,1
amesh,all
alls
type,2
mat,2
real,2
asel,s,loc,x,0,r1-5
extopt,esize,zzz/20
EXTOPT,ACLEAR,1
EXTOPT,ATTR,0,0,0
csys,0
VEXT,all, , ,0,0,zzz,,,,
eplot
alls
csys,1
type,3
mat,3
real,3
asel,u,loc,x,0,r1-5
asel,u,loc,z,1,zzz
extopt,esize,zzz/20
EXTOPT,ACLEAR,1
EXTOPT,ATTR,0,0,0
csys,0
VEXT,all, , ,0,0,zzz,,,,
alls
eplot
!接触单元设置
et,4,170
et,5,173
keyopt,5,9,0
keyopt,5,12,2
r,5
real,5
mp,mu,2,0.2
csys,1
vsel,s,type,,3
asel,r,loc,x,r1
aplot
cm,target,area
type,4
!设置目标单元
nsla,s,1
esurf,all
esel,s,type,,4
eplot
/psyms,esys,1
allsel
vsel,s,type,,2
asel,r,loc,x,r1
cm,contact,area
type,5
nsla,s,1
esurf,all
esel,s,type,,5
eplot
/psyms,esys,1
allsel
gplot
csys,1
type,1
mat,1
real,1
*do,i,0,15
*do,z,0,zzz/20-1
node1=node(r2,i*90/4,z*20)
node2=node(r2,i*90/4,(z+1)*20)
node4=node(r2,(i+1)*90/4,z*20)
node3=node(r2,(i+1)*90/4,(z+1)*20)
e,node1,node2,node3,node4
*enddo
*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!! 求解
NSEL,S,LOC,Z,0
NPLOT
D,all, , , , , ,ALL, , , , ,
ALLSEL,ALL
EPLOT
!————————
/solu
antype,static
outres,basic,all !输出选项
outpr,basic,all !窗口输出，这个没什么影响
autots,on !自动计算时间步长
nlgeom,on !打开是大变形，非线性
outres,all,all !输出所有子步结果
time,1
allsel
NSEL,S,LOC,Z,zzz
NPLOT
D,all, ,-25, , , ,UZ, , , , ,
nsubst,200 ,50 ,200 !荷载子步
neqit,20
eqslv,-1
ncnv,1 !终止分析选项，终止分析且退出程序
lnsrch,on !线性搜索，使计算稳定
CNVTOL,F, ,0.05,2, ,
allsel
solve
/POST26
numvar,10
csys,1
nsel,s,loc,z,zzz
*get,ncount,node,,count !获取被选择节点的数量定义到变量ncount中
nnow=0
nsol,6,ndnext(nnow),u,z,Zdisplancement
add,6,6,,,,,,-1
*do,i,1,ncount
nnow=ndnext(nnow)
rforce,8,nnow,f,z,zforce
add,7,7,8,,,,,,1
*enddo
xvar,6
/axlab,x,deflection(mm)
/axlab,y,totle load(N)
prod,7,,,,,,,-1
plvar,7
