求助一个非线性的钢混问题

yangerrol

他总是说我约束不够  算不出来  轻大家帮我看下  谢谢


FINI
/CLE
/PREP7
ET,1,SOLID65,,,,,,,1
ET,2,LINK8
MP,EX,1,13585
MP,PRXY,1,0.2
MP,EX,2,2e5
MP,PRXY,2,0.3
r,1,(0.003)**(0.003)*3.14
r,2,(0.004)**(0.004)*3.14
fc=14.3
ft=1.43
tb,concr,1
tbdata,,0.5,0.95,ft,-1
tb,miso,1,,14
tbpt,,0.0002,fc*0.19
tbpt,,0.0004,fc*0.36
tbpt,,0.0006,fc*0.51
tbpt,,0.0008,fc*0.64
tbpt,,0.0010,fc*0.75
tbpt,,0.0012,fc*0.84
tbpt,,0.0014,fc*0.91
tbpt,,0.0016,fc*0.96
tbpt,,0.0018,fc*0.99
tbpt,,0.002,fc
tbpt,,0.0033,fc*0.85
tb,biso,2
tbdata,,300,0
BLC4,,,0.2,0.75,0.2
wpoffs,0.035
wprota,,,90
vsbw,all
wpstyl,defa
wpoffs,0.1
wprota,,,90
vsbw,all
wpstyl,defa
wpoffs,0.165
wprota,,,90
vsbw,all
wpstyl,defa
wpoffs,,,0.035
vsbw,all
wpstyl,defa
wpoffs,,,0.1
vsbw,all
wpstyl,defa
wpoffs,,,0.165
vsbw,all
wpstyl,defa
*do,i,1,10   
wpoffs,,0.035+(i-1)*0.07
wprota,,90
vsbw,all
wpstyl,defa
*enddo
wpoffs,,0.75-0.035
wprota,,90
vsbw,all
wpstyl,defa
elemsize=0.025
lsel,s,loc,x,0.035
lsel,r,loc,z,0.035
LSEL,R,LOC,Y,0.035,0.75-0.035
cm,zj1,line
latt,2,2,2
lesize,all,elemsize
lsel,s,loc,x,0.1
lsel,r,loc,z,0.035
LSEL,R,LOC,Y,0.035,0.75-0.035
cm,zj2,line
latt,2,2,2
lesize,all,elemsize
lsel,s,loc,x,0.165
lsel,r,loc,z,0.035
LSEL,R,LOC,Y,0.035,0.75-0.035
cm,zj3,line
latt,2,2,2
lesize,all,elemsize
lsel,s,loc,x,0.035
lsel,r,loc,z,0.1
LSEL,R,LOC,Y,0.035,0.75-0.035
cm,zj4,line
latt,2,2,2
lesize,all,elemsize
lsel,s,loc,x,0.165
lsel,r,loc,z,0.1
LSEL,R,LOC,Y,0.035,0.75-0.035
cm,zj5,line
latt,2,2,2
lesize,all,elemsize
lsel,s,loc,x,0.035
lsel,r,loc,z,0.165
LSEL,R,LOC,Y,0.035,0.75-0.035
cm,zj6,line
latt,2,2,2
lesize,all,elemsize
lsel,s,loc,x,0.1
lsel,r,loc,z,0.165
LSEL,R,LOC,Y,0.035,0.75-0.035
cm,zj7,line
latt,2,2,2
lesize,all,elemsize
lsel,s,loc,x,0.165
lsel,r,loc,z,0.165
LSEL,R,LOC,Y,0.035,0.75-0.035
cm,zj8,line
latt,2,2,2
lesize,all,elemsize
lsel,s,tan1,y
lsel,r,loc,x,0.035,0.165
lsel,r,loc,z,0.035,0.165
lsel,r,loc,y,0.035,0.75-0.035
cm,gujin,line
latt,2,1,2
lesize,all,elemsize
vatt,1,,1
meshkey,1
esize,elemsize
vmesh,all
allsel,all
nsel,s,loc,y,0
d,all,all

nsel,s,loc,y,0.75
cm,top,node
d,top,ux,0.056
allsel,all

/solu
antype,0
nsubst,100
outres,all,all
autos,on
neqit,50
cnvtol,u,,0.015
solve

ztggzjq

是软件的问题，用ansys12.0输入应力应变曲线时不允许出现下降段，而10.0没有这个问题
tb,miso,1,,10
tbpt,,0.0002,fc*0.19
tbpt,,0.0004,fc*0.36
tbpt,,0.0006,fc*0.51
tbpt,,0.0008,fc*0.64
tbpt,,0.0010,fc*0.75
tbpt,,0.0012,fc*0.84
tbpt,,0.0014,fc*0.91
tbpt,,0.0016,fc*0.96
tbpt,,0.0018,fc*0.99
tbpt,,0.002,fc
可以按照以上运行。