请教大家一个问题 有关超弹性 多载荷步

mark98306

我的一个橡胶薄膜
有初始应力
在计算的时候
第一载荷步 我先加一个初应力  求解
第二载荷步  在上面初应力的基础上，使得上表面某些节点发生一定的强迫位移
这样应该不会有问题吧
为什么第二步总是实现不了


没什么提示
就是结果不对
默认是没有第一步计算的初始应力

/filnam,rubber2
/PREP7  
lamda=4.08
r1=75e-3
r0=r1/lamda
rc1=22e-3
t1=0.06e-3
t0=t1*lamda*lamda
rc0=rc1/lamda
ET,1,PLANE182   
KEYOPT,1,3,1
TB,HYPE,1,1,2,MOON  
TBDATA,,0.045e6,0.0412e6,0,,,
RECTNG,0,r1,0,t1
RECTNG,0,rc1,0,t1
allsel,
aovlap,all
lsel,s,loc,y,t1/2
LESIZE,all, , ,3
allsel
lsel,s,loc,x,rc1/2
LESIZE,all, , ,130,
lsel,s,loc,x,rc1+0.00001,r1-0.00001
LESIZE,all, , ,150,
allsel,
AMESH,all
FINISH  
/SOL
ANTYPE,0
solcontrol,on
AUTOTS,1
time,1
NLGEOM,1
NSUBST,1,,  
OUTRES,ALL,ALL  
NEQIT,20
lsel,s,loc,x,0
NSLL,S,1
D,all,UX,0
lsel,s,loc,x,r1
NSLL,S,1
D,all,ux,0
nsel,s,loc,y,0
dsym,symm,y
istress,0.243e8,,0.243e8
ALLSEL,ALL
solve
finish
/solu
PSTRESS,ON
solcontrol,on
AUTOTS,1
NLGEOM,1
NSUBST,2,100,2  
OUTRES,ALL,ALL  
lsel,s,loc,x,0
NSLL,S,1
D,all,UX,0
lsel,s,loc,x,r1
NSLL,S,1
D,all,ux,0
nsel,s,loc,y,0
dsym,symm,y
nsel,s,loc,y,t1
nsel,r,loc,x,0,rc1
d,all,uy,-t1/2
allsel,
solve

mark98306

第一步初始应力已经加上了
想在这个基础上 让上表面节点发生一定的强迫位移

mark98306

/filnam,rubber2
/PREP7  
lamda=4.08
r1=75e-3
r0=r1/lamda
rc1=22e-3
t1=0.06e-3
t0=t1*lamda*lamda
rc0=rc1/lamda
ET,1,PLANE182   
KEYOPT,1,3,1
TB,HYPE,1,1,2,MOON  
TBDATA,,0.045e6,0.0412e6,0,,,
RECTNG,0,r1,0,t1
RECTNG,0,rc1,0,t1
allsel,
aovlap,all
lsel,s,loc,y,t1/2
LESIZE,all, , ,3
allsel
lsel,s,loc,x,rc1/2
LESIZE,all, , ,130,
lsel,s,loc,x,rc1+0.00001,r1-0.00001
LESIZE,all, , ,150,
allsel,
AMESH,all
FINISH  
/SOL
ANTYPE,0
solcontrol,on
AUTOTS,1
time,1
NLGEOM,1
NSUBST,1,,  
OUTRES,ALL,ALL  
NEQIT,20
lsel,s,loc,x,0
NSLL,S,1
D,all,UX,0
lsel,s,loc,x,r1
NSLL,S,1
D,all,ux,0
nsel,s,loc,y,0
dsym,symm,y
istress,0.243e8,,0.243e8
ALLSEL,ALL
solve

!第二步，在上一步初始应力的基础上，让上表面的某些节点发生一定的强迫位移，结果怎么不对呢？
PSTRESS,ON  ！已经打开预应力效应开关了，怎么没有上一步的初始应力呢？
nsel,s,loc,y,t1
nsel,r,loc,x,0,rc1
d,all,uy,-t1/2
allsel,
solve

mark98306

我知道有一种方法可以实现
就是用两个工况
一个工况加载初应力
一个工况加强迫位移
然后两个工况叠加
不过这个方法不太准确

mark98306

/filnam,rubber2
/PREP7  
lamda=4.08
r1=75e-3
r0=r1/lamda
rc1=22e-3
t1=0.06e-3
t0=t1*lamda*lamda
rc0=rc1/lamda
ET,1,PLANE182   
KEYOPT,1,3,1
TB,HYPE,1,1,2,MOON  
TBDATA,,0.045e6,0.0412e6,0,,,
RECTNG,0,r1,0,t1
RECTNG,0,rc1,0,t1
allsel,
aovlap,all
lsel,s,loc,y,t1/2
LESIZE,all, , ,3
allsel
lsel,s,loc,x,rc1/2
LESIZE,all, , ,130,
lsel,s,loc,x,rc1+0.00001,r1-0.00001
LESIZE,all, , ,150,
allsel,
AMESH,all
FINISH  
/SOL
ANTYPE,0
solcontrol,on
AUTOTS,1
time,1
NLGEOM,1
NSUBST,1,,  
OUTRES,ALL,ALL  
NEQIT,20
lsel,s,loc,x,0
NSLL,S,1
D,all,UX,0
lsel,s,loc,x,r1
NSLL,S,1
D,all,ux,0
nsel,s,loc,y,0
dsym,symm,y
istress,0.243e8,,0.243e8
nsel,s,loc,y,t1
nsel,r,loc,x,0,rc1
d,all,uy,-t1/2
allsel,
solve