请教关于回弹问题的分析

forsetii

一张板经过上下模具压制成型后想考察撤掉模具后板本身的回弹量，请问这个过程该怎么做，ansys的静态分析可以完成吗？

下面是我做的压制过程命令流

/prep7

k,1,-300,150
k,2,225,150
k,3,750,150
k,4,-300,180
k,5,750,180
k,6,125,230

k,7,125,230
k,8,125,430
k,9,325,230
k,10,325,430
k,11,225,180

LARC,       9,       7,      11
LSTR,       7,       8  
LSTR,       8,      10  
LSTR,      10,       9  
LSTR,       4,       5  
LSTR,       5,       3  
LSTR,       3,       1  
LSTR,       1,       4  
FLST,2,4,4  
FITEM,2,3   
FITEM,2,2   
FITEM,2,4   
FITEM,2,1   
AL,P51X
FLST,2,4,4  
FITEM,2,8   
FITEM,2,5   
FITEM,2,7   
FITEM,2,6   
AL,P51X
!*  
VOFFST,1,100, ,

APLOT   
!*  
VOFFST,2,100, ,

SAVE






CYL4,0,0,150,0,0,180,-100   

CYL4,450,0,150,0,0,180,-100

et,1,solid95
mptemp,1,293,473,873,1273,1473,1773
mpdata,ex,1,1,2.1e11,1.95e11,1.55e11,0.45e11,0.01e11,0.0005e11  
mpdate,nuxy,1,1,0.3,0.3,0.3,0.3,0.3,0.3
mpdate,alpx,1,1,1.23e-5,1.23e-5,1.23e-5,1.23e-5,1.23e-5,1.23e-5
tb,biso,1,6,2   
tbtemp,293,1
tbdata,1,530e6,2.1e10   
tbtemp,473,2
tbdata,1,258e6,1.95e10  
tbtemp,873,3
tbdata,1,137e6,1.55e10  
tbtemp,1273,4   
tbdata,1,30e6,0.45e10   
tbtemp,1473,5   
tbdata,1,10e6,0.01e10   
tbtemp,1773,6   
tbdata,1,1e6,0.0005e10  
mp,dens,1,7815  
mp,kxx,1,43
mp,c,1,580  
mp,mu,1,0.3
MP,MU,1,0.3
MAT,1   
R,3
REAL,3  
et,3,targe170   
et,4,conta174   
KEYOPT,4,9,0
KEYOPT,4,10,1   

esize,25
mat,1
vsweep,1,1,3,1

esize,25
mat,1
vsweep,2,2,8,1

esize,25
mat,1
vsweep,3,13,14,1

esize,25
mat,1
vsweep,4,19,18,1

asel,s,,,6
nsla,s,1
cm,mubiao1,node

asel,s,,,15
nsla,s,1
cm,mubiao2,node

asel,s,,,20
nsla,s,1
cm,mubiao3,node

asel,s,,,10
nsla,s,1
cm,jiechu1,node

asel,s,,,12
nsla,s,1
cm,jiechu2,node

allsel,all
cmsel,s,mubiao1
type,3
mat,1
real,3
esurf
allsel,all
cmsel,s,jiechu1
type,4
mat,1
real,3
esurf
r,3,0,0,0.01
keyopt,4,5,1
keyopt,4,10,1
allsel,all

allsel,all
cmsel,s,mubiao2
type,3
mat,1
real,4
esurf
allsel,all
cmsel,s,jiechu2
type,4
mat,1
real,4
esurf
r,4,0,0,0.01
keyopt,4,5,1
keyopt,4,10,1
allsel,all

allsel,all
cmsel,s,mubiao3
type,3
mat,1
real,5
esurf
allsel,all
cmsel,s,jiechu2
type,4
mat,1
real,5
esurf
r,5,0,0,0.01
keyopt,4,5,1
keyopt,4,10,1
allsel,all

/solu
cmsel,s,mubiao1
d,all,ux,0
d,all,uy,-20
asel,s,,,21
nsla,s,1
d,all,ux,0
d,all,uy,0
d,all,uz,0

asel,s,,,22
nsla,s,1
d,all,ux,0
d,all,uy,0
d,all,uz,0

asel,s,,,17
nsla,s,1
d,all,ux,0
d,all,uy,0
d,all,uz,0

asel,s,,,16
nsla,s,1
d,all,ux,0
d,all,uy,0
d,all,uz,0

asel,s,,,15
nsla,s,1
d,all,ux,0
d,all,uy,0
d,all,uz,0

asel,s,,,20
nsla,s,1
d,all,ux,0
d,all,uy,0
d,all,uz,0

allsel,all
time,2
nsubst,20
autots,on
nropt,full,,off
nropt,unsym,,off
neqit
lnsrch,auto
nlgeom,on
outres,all,all
solve

pfsr

重新启动分析，读入db文件，upgeom使结构configuration回到上个分析结束是的状态，从rst文件读入上个分析结束是壳单元的厚度和应力，用rimport命令。save，solve，OK.

pfsr

我在上一个回复中答错了，开始以为你用的壳单元，而且rimport是一个跟DYNA相关的命令。对于你这个问题，再加一步load step即可，如下。

/solu
antype,0,rest,1,9,continue     ！接着开始分析

cmsel,s,mubiao1
d,all,uy,20                         ！移开上面的模板
allsel,all
da,9,all,0.0                         ！约束变形板的两端，或你可以在合适的地方加约束。
da,11,all,0.0                       ！必须要有约束，不然板就会发生刚体位移

time,1.0
nsubst,20
nlgeom,on
solve
save

/post1
set,2,last
plnsol,eppl,eqv
andata,0.5,,,,,1            ！查看动画

forsetii

pfsr ,您好,谢谢您的回复,给了我很大帮助.
一点疑问,面9和11在给了全约束之后板本身的回弹行为会受到影响,无法考察真实的回弹量和应力.只给面9或者面11约束的办法我也试过,但是板出现了弯曲.是否可以通过给出板的密度,通过自身的重力与下面两个柱体表面接触,作为y向的约束呢?

另外您提到的会出现刚体位移的情况如果在ls-dyna里面来做是否可以得到解决?我是个初学者,对这些缺少基本的概念,给您添麻烦了.

pfsr

我试了一下，这个分析考虑自重就没问题了。但有的分析自重小回弹力大，还是会出现刚体位移的。
你可以试一下施加对称约束，或是在不同的面或线上施加ux,uy,uz约束，使回弹尽量不受约束的影响。还可以试一下只把上面的模板移开一点，使之在回弹的时候再一次发生接触。
如果不加约束ls-dyna也会出现刚体位移，但如果把时间限制的很小，刚体位移就会很小，对分析影响不大。或者回弹的刚体位移就是想要的分析结果。