ex3-3.txt

《ANSYS工程分析进阶实例》-王呼佳-源代码 看文件吧.不多说了.

!得到重力作用下的应力场
fini
/cle

/filname,try01,1
/prep7
et,1,42
mp,ex,1,2e11
mp,nuxy,1,0.3
mp,dens,1,7800
blc4,,,1,2
esize,0.1
amesh,all
fini

/solu
nsel,s,loc,y,
d,all,uy
nsel,s,loc,x,
nsel,a,loc,x,1,
d,all,ux
alls
acel,,10,
iswrite,1
solv
fini

/post1
plnsol,s,eqv
PLNSOL, U,SUM
PLNSOL, U,y
fini

!===================
!调入初应力文件法 
fini
/cle

/filname,try02,1
/prep7
et,1,42
mp,ex,1,2e11
mp,nuxy,1,0.3
mp,dens,1,7800
blc4,,,1,2
esize,0.1
amesh,all
fini

!全部结点全约束，读初应力文件try01.ist
!计算节点反力就是重力场作用下的等效节点力
/solu
d,all,all
isfile,read,try01,ist,,2
isfile,list
outres,all,all
alls
solve

/post1
set,last
PLNSOL, U,SUM
PLNSOL, U,y
plnsol,s,eqv
finish

!施加等效节点力
/filname,try03,1
/solu
lsclear,all
nsel,s,loc,y,
d,all,uy
nsel,s,loc,x,
nsel,a,loc,x,1,
d,all,ux
alls
ldread,reac,1,1,,,'try02','rst',
solv

SFL,3,PRES,1e5
alls
solve
fini

/post1
set,last
PLNSOL, U,SUM
PLNSOL, U,y
plnsol,s,eqv
finish
!＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝
!常规方法
!重力+压力
fini
/cle

/filname,try04,1
/prep7
et,1,42
mp,ex,1,2e11
mp,nuxy,1,0.3
mp,dens,1,7800
blc4,,,1,2
esize,0.1
amesh,all
fini

/solu
nsel,s,loc,y,
d,all,uy
nsel,s,loc,x,
nsel,a,loc,x,1,
d,all,ux
alls
acel,,10,
SFL,3,PRES,1e5
solv
fini

/post1
plnsol,s,eqv
PLNSOL, U,SUM
PLNSOL, U,y
fini