关于仿真正方形截面定常风荷载升力阻力系数问题

lorrykey

最近在用ansys flotran 计算正方形截面定常风荷载三分力系数随风攻角的变化图
但是总是不收敛，哪位达人有这方面的资料或是实例可供参考一下，谢谢了
这是我计算用的命令流，欢迎大家讨论！

finish
/clear,all
/filname,aoa_v5,1
err1=1e-4
err2=1e-4
width=0.25
length=0.2
radius=0.01
mesh1=0.1
mesh2=0.001
mesh3=0.0005
velo=5
degree=1
pi=3.1415926535897932384626433832795
aoa=3.1415927*degree/180
ux=velo*cos(aoa)
uy=velo*sin(aoa)
/prep7
et,1,fluid141
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,1,,1.225
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,VISC,1,,1.79e-5  
rectng,-2.5,2.5,-2,2
k,100,-length/2,-width/2
k,101,length/2,-width/2
k,102,length/2,width/2
k,103,-length/2,width/2
l,100,101
l,101,102
l,102,103
l,103,100
lfillt,8,5,radius,,
lfillt,5,6,radius,,
lfillt,6,7,radius,,
lfillt,7,8,radius,,
al,5,6,7,8,9,10,11,12
asba,1,2
save
lsel,,,,1,4
lesize,all,mesh1
lsel,,,,5,8
lesize,all,mesh2
lsel,,,,9,12
lesize,all,mesh3
MSHAPE,0,2D
mshkey,0
amesh,all
allsel,all
finish
/solu
dl,4,,vx,ux,1
dl,4,,vy,uy,1
!dl,1,,vx,0,1
!dl,1,,vy,0,1
!dl,3,,vx,0,1
!dl,3,,vy,0,1
dl,1,,pres,0,1
dl,2,,pres,0,1
dl,3,,pres,0,1
dl,5,,vx,0,1
dl,5,,vy,0,1
dl,6,,vx,0,1
dl,6,,vy,0,1
dl,7,,vx,0,1
dl,7,,vy,0,1
dl,8,,vx,0,1
dl,8,,vy,0,1
dl,9,,vx,0,1
dl,9,,vy,0,1
dl,10,,vx,0,1
dl,10,,vy,0,1
dl,11,,vx,0,1
dl,11,,vy,0,1
dl,12,,vx,0,1
dl,12,,vy,0,1
!* property of the fluid
!*  
FLDATA12,PROP,DENS,4
FLDATA12,PROP,VISC,4
!*  
FLDATA7,PROT,DENS,AIR-SI
FLDATA8,NOMI,DENS,-1
FLDATA7,PROT,VISC,AIR-SI
FLDATA8,NOMI,VISC,-1
FLDATA7,PROT,COND,CONSTANT  
FLDATA8,NOMI,COND,-1,   
FLDATA7,PROT,SPHT,CONSTANT  
FLDATA8,NOMI,SPHT,-1,
!*
fldata2,iter,exec,20000
fldata3,term,vx,err2
fldata3,term,vy,err2
fldata3,term,pres,err1
fldata3,term,temp,err2
fldata5,outp,sumf,50
fldata34,mir,mome,0.6
FLDATA25,RELX,PRES,0.8
!
solve
finish
/post1
set,last
lsel,s,line,,5,12
nsll,s,1
nplot,0
intsrf,pres
*get,fx,intsrf,0,pres,fx
*get,fy,intsrf,0,pres,fy
*get,mz,intsrf,0,pres,mz
fd=fx*cos(aoa)+fy*sin(aoa)
fl=-fx*sin(aoa)+fy*cos(aoa)
rho=1.225
cd=2*fd/velo/velo/rho/width
cl=2*fl/velo/velo/rho/width
cz=2*mz/velo/velo/rho/width/length
*cfopen,resu_v5_presBC,txt,,append
*cfwrite,t,degree,fx,fy,fd,fl,cd,cl,cz
*cfclos