Ansys非線性靜態(tài)實(shí)例分析(命令流方式)

2013-08-14 by:CAE分析培訓(xùn)中心來源:仿真在線

Ansys非線性靜態(tài)實(shí)例分析(命令流方式)

問題描述:
    一個(gè)周邊簡支的圓盤,在其中心受到一個(gè)沖桿的周期作用。由于沖桿被假定是剛性的,因此在建模時(shí)不考慮沖桿,而將圓盤上和沖桿接觸的結(jié)點(diǎn)的Y 方向上的位移耦合起來。
    由于模型和載荷都是軸對(duì)稱的,因此用軸對(duì)稱模型來進(jìn)行計(jì)算。求解通過四個(gè)載荷步實(shí)現(xiàn)。
問題詳細(xì)說明:
材料性質(zhì):
EX=70000 (楊氏模量)
NUXY=0.325(泊松比)
塑性時(shí)的應(yīng)力—應(yīng)變關(guān)系如下:
           應(yīng)變               應(yīng)力
          0.0007857           55
          0.00575             112
0.02925             172
             0.1               241
加載歷史:
          時(shí)間               載荷
           0                  0
           1                - 6000

2 750
3 -6000
問題描述圖:

Ansys非線性靜態(tài)實(shí)例分析(命令流方式)ansys分析圖片1

Ansys命令流:

fini
/cle
/title,circular plate loaded by a circular punch - kinematic hardening
rpl=65
rpu=5
h=6.5
exx=70000
sts1=55      !yield stress
stn1=sts1/exx
sts2=112
stn2=0.00575
sts3=172
stn3=0.02925
sts4=241
stn4=0.1

nex=15          ! No. of elements along the radius
net=2            ! No. of elements in the plate's thickness

/prep7
et,1,42,,, 1
! define node for convenient postprocessing
n,1,rpl,0
n,2,0,0
n,3,0,h/2

!define geometry
k,1,, -(h/2)
k,2,rpu,-(h/2)
k,3,rpl, -(h/2)
kgen,2,1,3,1,,h,,3
nex1=nint(0.8*net)
nex2=nex -nex1

l,1,2
l,4,5
l,2,3
l,5,6

a,1,2,5,4
a,2,3,6,5

lesize,1,,,nex1
lesi ze,2,,,nex1
lesize,3,,,nex2,2.5
lesize,4,,,nex2,2.5

esize,(h/net)
amesh,all
nummrg,node
nodes
nsel,s,loc,x,0,rpu
nsel,r,loc,y,(h/2)
cp,1,uy,all
nsel,all
esel,all
fini

/solu
nsel,s,loc,x,rpl
nsel,r,loc,y,0
d,all,uy
nsel,s,loc,x,0
d,all,ux
nsel,all
outres,all,all
fini

/prep7
mp,ex,1,exx
mp,nuxy,1,0.325

tb,mkin,1
tbtemp,,strain
tbdata,,stn1,stn2,stn3,stn4
tbtemp,,
tbdata,,sts1,sts2,sts3,sts4
fini

/solu
pred,on
outres,all,all
nsubst,1
time,1e -6
f,3,fy,0
lswrite

autots,on
nsubst,10
time,1
f,3,fy, -6000
lswri te

time,2
f,3,fy,750
lswrite

time,3
f,3,fy, -6000
lswrite

lssolve,1,4
fini

/post1
set,2
/dscal,1,1
pldisp,2
fini

/post26
nsol,2,2,u,y,uy2
rforce,3,1,f,y,ry1
add,2,2,,,uy2,,, -1
/grid,1
/axlab,x,deflection [mm]
/axlab,y,force [n]
xval,2
plvar,3
prvar,2,3
fini