呼,论坛再一次好 不热闹。今天尝试了一下之前做过的东西,额,关于固有频率的,权当砖头抛上来吧。
/prep7
et,1,solid45
et,2,fluid30
et,3,fluid30,,1
mp,ex,1,2.1e11
mp,prxy,1,0.3
mp,dens,1,7800
mp,dens,2,789
mp,sonc,2,1400 !
mp,dens,3,789
mp,sonc,3,1400
!mp,mu,1,,0.5 !不影响结果
!mp,mu,2,,0.5
*SET,r,0.01
*SET,t,0.001
*SET,l,0.2
!geometry
cylind,r-t,r,0,l,-90,90
cm,pipe,volu !
cylind,0,r-t,0,l,-90,90
vsel,s,,,pipe !
vsel,invet
cm,water,volu
nummrg,all !重合的线只留下一个,保证结合面节点重合。vital
!mesh the pipe
vsel,s,,,pipe
vatt,1,,1
mshape,0,2d
mshkey,1
LESIZE,6, , ,2
LESIZE,8, , ,2
LESIZE,7, , ,20
LESIZE,11, , ,100
vsweep,all
!mesh the water
vsel,s,,,water
R,3,1e5,
vatt,2,3,3 !分配材料、实常数、单元
mshape,0,2d
mshkey,1
LESIZE,7, , ,20
LESIZE,11, , ,100
vsweep,all
! ACEL,0,9.8,0, !重力加速度
!
!dof
asel,s,loc,z,0
da,all,uz,0
da,all,ux,0
da,all,uy,0
da,all,ROTX
da,all,ROTY
!to change the dof,try to set next...
asel,s,loc,z,l
da,all,uz,0
da,all,ux,0
da,all,uy,0
da,all,rotx,0
da,all,roty,0
!
asel,s,,,5,6,1
da,all,symm
asel,s,,,10,11
da,all,symm
sfa,4,1,pres,1e8 !此处面的正方向和空管相反,故为1e8。
!to note SFI nodes
csys,1
nsel,s,loc,x,r-t-0.00025,r-t+0.00025
esln,s !
esel,r,type,,3 !基于上一步选择单元类型为3
emodif,all,type,2 !
sf,all,fsi !加载流固耦合fluid/ansys-field surface-
allsel,all
finish
/solu
solve
finish
/SOL
antype,2
PSTRES,1 !预应力打开
modopt,unsym,10
mxpand,10,,,0
solve |