单点重启分析的问题

terryy11

在压杆顶端，分两次施加两个力，一个是3×P，一个是P。施加3×P后solve$save一次，然后resume$antype,,rest，再施加P，再solve。按理说第二步求解之后，得到的力和位移肯定比第一步的大，但是读取第一步结果，plnsol,s,eqv后显示DMX=0.004367,SMN=852.598,SMX=.364E+09,而读取第二步结果时，结果是：
DMX=.867E-3,SMN=21.832,SMX=.714E+08,反而还比第一步只有3×P时小了，也就是这两步荷载是分别计算，没有叠加在一起？那resume$antype,,rest起的是什么作用？

附上求解部分命令流



finish
!/clear
/filname,add1
/pmeth,on,0
/solu



allsel
wpcsys,-1,0
csys,4
nsel,s,loc,z,0
nsel,r,loc,x,0
d,all,ux
d,all,uy
d,all,uz
nsel,s,loc,z,h
nsel,r,loc,x,0
d,all,ux
d,all,uy
nsel,s,loc,z,h
nsel,r,loc,x,0
*get,nn1,node,0,count
f,all,fz,-5*p/nn1  allsel
ANTYPE,STATIC
rescontrol,,all,1
NLGEOM,ON
AUTOTS,on
allsel
                           !No.1     !No.1     !No.1     !No.1     !No.1
outres,all,all
solve
save

finish
/clear
/filname,add1
/pmeth,on,0
/solu
resume
antype,,rest
NLGEOM,ON
allsel
nsel,s,loc,z,h
nsel,r,loc,x,0
*get,nn1,node,0,count
f,all,fz,-p/nn1
allsel

AUTOTS,on
allsel
outres,all,all
solve

jg_meng

操作没问题，是你理解的问题，先施加3p的载荷，结构有了变形，重启动后换成1p的载荷，结构变形必然没有1p时候的大（弹性变形部分要回弹），ansys默认是替换载荷，如果想看累加效果，可以通过【DCUM】设置，或者做载荷工况然后进行工况组合查看结果

terryy11

感谢斑竹回答！
我按书上的单点启动例子试验了下：
（1）重启前后两个step分别在同点加同方向，同大小力P，分10步加（nsubst,10），每个子步（time增量=0.1）产生位移50，到time=1.0时位移为500。在step1的时候系统会计算，到了step2，所有substep下的边行都与step1一样，也就是step2的力没有进行计算
（2）重启前后两个step分别在同点加反方向，同大小力P，分10步加（nsubst,10）。step1计算与上面相同，step2的时候,time=1.1时，产生的位移是反方向的100，最后的位移是与step1反方向的500.
        据此推算，如果没有FCUM,ALL，那同点加的两个step的荷载，到了step2系统会计算前后step那个节点上力的差值：在（1）的时候，P1-P2=0,因此step2没有进行计算，所有子步位移都与step1一样，没变化。在（2）的时候，P1-P2=-2P，因此到了step2，每个子步的力不是0.1×P，而是0.1×2P，因此产生反方向的位移才是100，而不是50.不知这样理解是否正确？
（3）加上FCUM,ALL以后，同点加同方向同大小的两个力也能叠加了，结果相当于直接施加2P。
（4）又用shell181单元建立了一个一端固定一段自由的板模型，这时候用上FCUM也无法让同点的两个力叠加，这是为什么？麻烦再帮助看看是哪里问题，谢谢！

附上两个命令流：
1.书中的命令流
finish
/clear
/filname,rtest1
l0=1000
b0=10
h0=20
/prep7

dofsel,all
fcum,add

k,1
k,2,l0
l,1,2
et,1,beam3
mp,ex,1,2.0e5
mp,prxy,1,0.3
r,1,b0*h0,b0*h0*h0*h0/12,h0
lesize,all,,,10
lmesh,all
d,1,all



/solu
antype,0
!nlgeom,on
outres,all,all
autots,off
nsubst,10
f,2,fy,2000
solve
save
finish


!/clear  
/filname,rtest1
/solu
resume
antype,,rest
!f,7,fy,4440
f,2,fy,2000
nsubst,20
solve
save
finish

2.SHELL181命令流
finish
/clear
/filname,be test
/prep7
h=1
b=0.5
t=0.006
p=500


et,1,shell181
ET,7,MESH200
keyopt,7,1,7

UIMP,1,EX, , ,2.06e11,   
UIMP,1,DENS, , ,7850,  
UIMP,1,NUXY, , ,0.3,
mp,mu,1,0.5  

TB,BKIN,1,1,2,1
TBDATA,,460E6,0,,,,   
r,1,t,t,t,t

k,1,-b/2,0,0
k,2,b/2,0,0
k,3,b/2,0,h
k,4,-b/2,0,h
l,1,2
l,2,3
l,3,4
l,4,1
al,all


allsel
aatt,1,1,1

lsel,all
lesize,all,,,10


mopt,QMESH,MAIN
!mopt,trans,2
MSHAPE,0,2D
MSHKEY,2
amesh,all
eplot

allsel
nummrg,all

nsel,all
allsel
dofsel,all
fcum,add

dcum,add

!  No.1     No.1     No.1     No.1     No.1     No.1     No.1     
wpcsys,-1,0
csys,4
nsel,s,loc,z,0
d,all,all

/solu

antype,static
nsel,s,loc,z,h
f,all,fy,2*p/10

allsel

NLGEOM,ON

autots,off
nsubst,10
outres,all,all
solve
save

!   No.2     No.2     No.2     No.2     No.2     No.2     No.2     No.2     No.2     
finish
!/clear,start
/filname,be test
/solu
resume
antype,,rest,,,
wpcsys,-1,0
csys,4


nsel,s,loc,z,h
f,all,fy,-2*p/10
NLGEOM,ON

autots,off
nsubst,10
outres,all,all
allsel
solve
finish

523402750

这个nsubst与autos也有关，继续努力