[求助]请大家帮我看看这个命令流,为什么我中间的圆柱固体就是加不上角速度

niathjing

我要做的是中间一个旋转的圆柱体,流体以一定的速度经过后,看圆柱体周围流体的流动情况.可是我搞不好怎么把旋转坐标系加到中间圆柱上,让他能以一定的角速度转动.我看帮助里给的两个关于旋转的命令(Cgomga   Cgloc)我加上去了还是没反应,弄了2天了也不转,实在是愚钝啊.请各位大侠帮忙指点一下,给加段命令流.让我圆柱转起来吧.
/prep7
rho = 1000 ! Density
grav = 9.81 ! Gravity
mu = 0.0008 ! Viscosity
R = 0.03 ! Bump radius
L = 2*R
H = 0.075 ! Inlet height
outlen = 5.0*R ! Outlet length
inlen = 5.0*R ! Inlet length
ny1 = 8 ! Number of elements in line 1
ny2 = 2 ! Number of elements in line 2
ny3 = 4 ! Number of elements in line 3
ry2 = 1.0 ! Ratio of element size in line 2
ry3 = 1.0 ! Ratio of element size in line 3
nx1 = 15 ! Number of elements in line 13
nx2 = 15 ! Number of elements in line 21
nx3 = 20 ! Number of elements in line 17
rx1 = 0.5 ! Ratio of element size in line 13
rx2 = 1.0 ! Ratio of element size in line 25
rx3 = 1.5 ! Ratio of element size in line 17
nr1 = 10 ! Number of elements in radial direction near bump
et,1,141
k,1,-inlen,0.0 ! KP 1 - 4 are at the inlet
k,2,-inlen,L
k,3,-inlen,H
k,4,-inlen,-L
K,5,-inlen,-H
l,1,2 ! Lines 1-3 are the inlet line
l,2,3
l,1,4
l,4,5
k,6,-L,0.0 ! KP 5 - 8 are before the bump
k,7,-L,L
k,8,-L,H
k,9,-L,-L
k,10,-L,-H
!
l,6,7
l,7,8
l,6,9
l,9,10
!
k,11,L,0.0 ! KP 9 - 12 are after the bump
k,12,L,L
k,13,L,H
k,14,L,-L
k,15,L,-H
!
l,11,12 ! Lines 7-9 after the bump
l,12,13
l,11,14
l,14,15
!
k,16,outlen,0.0 ! KP 13 - 16 are at the outlet
k,17,outlen,L
k,18,outlen,H
k,19,outlen,-L
k,20,outlen,-H
!
l,16,17 ! Lines 10-12 at the outlet
l,17,18
l,16,19
l,19,20
k,21,
circlr,21,R,,,360,6
l,1,6 ! Lines 13-16 are horizontal lines near inlet
l,2,7
l,3,8
l,4,9
l,5,10
l,8,13
l,7,12
l,9,14
l,10,15
l,13,18
l,12,17
l,11,16
l,14,19
l,15,20
L,7,24
L,6,25
L,9,26
L,14,27
L,11,22
L,12,23
al,1,24,5,23
al,2,25,6,24
al,3,23,7,26
al,4,26,8,27
al,6,28,10,29
al,8,30,12,31
al,10,32,14,33
al,9,33,13,34
al,11,34,15,35
al,12,35,16,36
al,17,41,9,42
al,18,42,29,37
al,19,37,5,38
al,20,38,7,39
al,21,39,30,40
al,22,40,11,41
lsel,s,,,1,15,2
lesize,all,,,ny1
lsel,s,,,23,36
lesize,all,,,nx1
lsel,s,,,2,16,2
lesize,all,,,ny2
lsel,s,,,37,42
lesize,all,,,NR1
lsel,s,,,17,19,2
lsel,a,,,20,22,2
lesize,all,,,ny1
!lsel,s,,,18,21,3
!lesize,all,,,nx1
allsel
amesh,all
! Inlet boundary condition
lsel,s,,,1,3,2
nsll,S,1
d,all,vx,1.0
d,all,vy
!nsel,s,loc,x,-inlen
!nsel,r,loc,y,0.0
!d,all,vx,0.001
! boundary volume fraction (SFL and SFE commands)
sfl,1,vfrc,1,,1
sfl,3,vfrc,1,,1
!!!! The walls
lsel,s,,,17,22
lsel,a,,,26,30,4
lsel,a,,,35
lsel,a,,,24,29,5
lsel,a,,,33
nsll,S,1
d,all,vx,0
d,all,vy,0
allsel
Keyopt,1,3,0
Keyopt,1,4,1
Cgomga,0,-1000,0,
Cgloc,0,1e-010,0,
/solu
!!!! FLOTRAN Input
acel,0.0,0.0,0.0
FLDATA1,SOLU,ALE,1
FLDA,SOLU,VOF,T
FLDA,SOLU,FLOW,T
FLDA,SOLU,TRAN,T
FLDA,SOLU,TURB,T
FLDA,TIME,GLOB,10
FLDA,TIME,NUMB,50
FLDA,TIME,STEP,0.01
FLDA,TIME,APPE,0.1
FLDA,TIME,PRES,1.0E-6
FLDA,ADVM,MOME,SUPG
flda,relx,vx,0.5
flda,relx,vy,0.5
flda,relx,pres,1.0
flda,prot,dens,constant ! Constant density
flda,prot,visc,constant ! Constant viscosity
flda,nomi,dens,rho ! Density set earlier
flda,nomi,visc,mu ! Viscosity set earlier
solve
fini

niathjing

这个是图片

niathjing

怎么没传上去.急!

niathjing

顶一下

niathjing

晕.怎么没人帮忙看啊