Ansys非线性静态实例分析(命令流方式)

2013-08-14 by:CAE分析培训中心来源:仿真在线

Ansys非线性静态实例分析(命令流方式)

问题描述:
    一个周边简支的圆盘,在其中心受到一个冲杆的周期作用。由于冲杆被假定是刚性的,因此在建模时不考虑冲杆,而将圆盘上和冲杆接触的结点的Y 方向上的位移耦合起来。
    由于模型和载荷都是轴对称的,因此用轴对称模型来进行计算。求解通过四个载荷步实现。
问题详细说明:
材料性质:
EX=70000 (杨氏模量)
NUXY=0.325(泊松比)
塑性时的应力—应变关系如下:
           应变               应力
          0.0007857           55
          0.00575             112
0.02925             172
             0.1               241
加载历史:
          时间               载荷
           0                  0
           1                - 6000

2 750
3 -6000
问题描述图:

Ansys非线性静态实例分析(命令流方式)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