shimering 发表于 2008-1-16 10:10

一个圆柱加流体,不太会,把命令流放上,大家帮我看看吧

做论文用,大家帮帮忙,就是在水中的一个圆柱体,大家帮我看看吧:-(


FINISH
/CLEAR
/FILENAME,CUBEBOX
/UNITS,SI
/TITLE,CUBEBOX
/PREP7
!------------------------------------------------
ET,24,24
ET,63,63
MP,EX,1,6.9E10
MP,DENS,1,2700
MP,NUXY,1,0.33
r,6,0.006
r,7,0.009             !板厚
r,8,0.012

!==================================================
!------------------------------------------------
! start to generate k points
Z1=0.18288
*do,i,1,48,1
k,i,0.6096*cos(i*7.5*3.14159265/180),0.6096*sin(i*7.5*3.14159265/180),0
*enddo
kgen,21,1,48,1,0,0,0.18288,100
*do,i,0,19,1
a,1+i*100,101+i*100,102+i*100,2+i*100
*repe,47,1,1,1,1
*enddo
k,50,0,0,0
k,51,0,0,0.18288*20
a,48,148,101,1
*repe,20,100,100,100,100
ASEL,ALL
CM,OUHEMIAN1,AREA
a,50,1,2,3
*repe,23,0,2,2,2
a,50,47,48,1
a,51,2000+1,2000+2,2000+3
*repe,23,0,2,2,2
a,51,2000+47,2000+48,2000+1
ASEL,S,LOC,Z,0
AADD,ALL
ALLS
asel,S,loc,z,20*z1
AADD,ALL
ALLS
ASEL,S,LOC,Z,0
ASEL,A,LOC,Z,20*Z1
CM,OUHEMIAN2,AREA
*do,i,10001,10048,1            !建流体边界点
k,i,1*cos((i-10000)*7.5*3.14159265/180),1*sin((i-10000)*7.5*3.14159265/180),0
*enddo
kgen,21,10001,10048,1,0,0,0.18288,100

*do,i,0,19,1
v,1+100*i,2+100*i,10002+100*i,10001+100*i,101+100*i,102+100*i,10102+100*i,10101+100*i
*repe,47,1,1,1,1,1,1,1,1
v,48+100*i,1+100*i,10001+100*i,10048+100*i,148+100*i,101+100*i,10101+100*i,10148+100*i
*enddo

*do,i,110001,110048,1            !建流体外边界点
k,i,11*cos((i-110000)*7.5*3.14159265/180),11*sin((i-110000)*7.5*3.14159265/180),0
*enddo
kgen,21,110001,110048,1,0,0,0.18288,100


*do,i,0,19,1
v,10001+100*i,10002+100*i,110002+100*i,110001+100*i,10101+100*i,10102+100*i,110102+100*i,110101+100*i
*repe,47,1,1,1,1,1,1,1,1
v,10048+100*i,10001+100*i,110001+100*i,110048+100*i,10148+100*i,10101+100*i,110101+100*i,110148+100*i
*enddo


numstr,line,30000
k,300060,0,0,20*z1
k,300061,0,0,20*z1+1-0.6096
l,300060,300061
k,300062,0,0,0
k,300063,0,0,-1+0.6096
l,300062,300063
asel,s,loc,z,0
vdrag,all,,,,,,30001
asel,s,loc,z,20*z1
vdrag,all,,,,,,30000
k,300064,0,0,20*z1+10
k,300065,0,0,-10
l,300061,300064
l,300063,300065
asel,s,loc,z,20*z1+1-0.6096
vdrag,all,,,,,,30770
asel,s,loc,z,-1+0.6096
vdrag,all,,,,,,30771



NUMMRG,ALL
!*************************MESH***********************
alls
!!!!!!!!!!!!!!!流体
et,30,fluid30         !acoustic elements contact with structure interface
et,31,fluid30,,1,0      !acoustic element in the interior
et,45,solid45
et,130,fluid130
!!the reference pressure in water
r,1,1e-6
mp,ex,1,2.1e11
mp,dens,1,7850
mp,nuxy,1,0.3
mp,dens,2,1000
mp,sonc,2,1500

K,1000001,0,0,1000000
K,100002,0,0,2*Z1
K,100004,0,0,4*Z1
K,100006,0,0,6*Z1
K,100008,0,0,8*Z1
K,100010,0,0,10*Z1
K,100012,0,0,12*Z1
K,100014,0,0,14*Z1
K,100016,0,0,16*Z1
K,100018,0,0,18*Z1
MSHKEY,0
ESIZE,1

asel,s,,,ouhemian1
TYPE,63
MAT,1
REAL,6
AMESH,ALL
ALLS

SMRTSIZE,10
asel,s,,,ouhemian2
TYPE,63
MAT,1
REAL,6
AMESH,ALL
ALLS
NUMMRG,ALL

alls
SMRTSIZE,OFF
CSYS,1                  !mesh流体的时候总是有问题,
esize,,1
mshape,1,3d
TYPE,30
REAL,30
MAT,2
esize,,1
vsel,S,LOC,X,0,1
VSEL,R,LOC,Z,-1+0.6096,20*Z1+0.6096
vmesh,all

vsel,all
vsel,r,loc,z,20*z1+1-0.6096,20*z1+11
vadd,all
esize,,1
mshape,1,3d
TYPE,30
REAL,30
MAT,2
esize,,1
vmesh,all

vsel,all
vsel,r,loc,z,-1+0.6096,-13
vadd,all
esize,,1
mshape,1,3d
TYPE,30
REAL,30
MAT,2
esize,,1
vmesh,all
vsel,all
vsel,r,loc,x,1,11
vadd,all
esize,,1
mshape,1,3d
TYPE,30
REAL,30
MAT,2
esize,,1
vmesh,all
*go,:p
vsel,s,loc,x,0,1
vsel,r,loc,z,-1+0.6096,20*z1+1-0.6096
nslv,s,1            !直接eslv???
esln                !直接eslv???
esel,inve
esel,u,type,,63
TYPE,31
REAL,30
MAT,2
EMODIF,ALL
alls            
asel,s,,,ouhemian1
asel,a,,,ouhemian2
nsla,s,1
sf,all,fsi
asel,s,loc,z,20*z1+10
nsla,s,1
D,all,PRES,0         ! SET PRESSURE AT OUTER RADIUS TO ZERO
alls

finish
/SOLU      
alls            
ANTYPE,MODAL   
MODOPT,unsym,10,0
MXPAND,10,0
LUMPM,off   
SOLVE            
FINISH

:p
页: [1]
查看完整版本: 一个圆柱加流体,不太会,把命令流放上,大家帮我看看吧