PROFILE (1051094), страница 2
Текст из файла (страница 2)
fldata14,temp,bulk,tem
fldata26,stab,pres,1e-4
fldata26,stab,turb,1e-4
ca=sqrt(temp0*cpheat)-1.0
fldata31,capp,velo,1
ca=2*vel
fldata31,capp,umin,-ca
fldata31,capp,umax, ca
fldata31,capp,vmin,-ca
fldata31,capp,vmax, ca
vvx=vel*cos(alfa*pi/180) !velocity components
vvy=vel*sin(alfa*pi/180) !velocity components
!Setting boundary conditions
alls
/prep7
ddele,all,all
/solu
an_div=60/(n3/2) !angle between 2 nodes in the farfield
csys,1
cmsel,s,ext !node selection at the inlet
nsel,r,loc,y,90+alfa-an_div,270+alfa+an_div
cm,inle,node
cmsel,s,ext !node selection at the outlet
nsel,u,loc,y,90+alfa-an_div,270+alfa+an_div
cm,outl,node
csys,0
cmsel,s,inle
d,all,vx,vvx
d,all,vy,vvy
d,all,enke,kk
d,all,ends,ep
cmsel,s,outl
d,all,pres,0
cmsel,s,wall
d,all,vx,0
d,all,vy,0
alls
/pbc,all,,1
ic,all,vx,0
ic,all,vy,0
ic,all,pres,0
ic,all,ends,0
ic,all,enke,0
ic,all,temp,temp0
solve
/post1
set,1
cmsel,s,wall
intsrf,pres
*get,fpx,intsrf,,item,fx
*get,fpy,intsrf,,item,fy
*get,mpz,intsrf,,item,mz
fx=fpx
fy=fpy
m_025=mpz
pi=3.1415
ffx=fx*cos(alfa*pi/180)+fy*sin(alfa*pi/180)
ffy=-fx*sin(alfa*pi/180)+fy*cos(alfa*pi/180)
c=vel*vel*rho*lc/2
cl=ffy/c
cd=ffx/c
cm=m_025/c
cmsel,s,w_up
*get,num1,node,,count
*dim,pw1,table,num1,1,1
*dim,xx1,table,num1,1,1
*dim,n11,array,num1,1,1
*do,i,1,num1
cm,current,node
*get,mnx,node,,mnloc,x
nd=node(mnx,0,0)
*do,k,1,5
*if,nx(nd),ne,mnx,then
nsel,u,,,nd
nd=node(mnx,0,0)
*endif
n11(i)=nd
*enddo
cmsel,s,current
xx1(i)=nx(nd)+0.25*lc
nsel,u,,,nd
*enddo
cmsel,s,w_do
*get,num2,node,,count
*dim,pw2,table,num1,1,1
*dim,xx2,table,num1,1,1
*dim,n22,array,num1,1,1
*do,i,1,num2,1
cm,current,node
*get,mnx,node,,mnloc,x
nd=node(mnx,0,0)
*do,k,1,5
*if,nx(nd),ne,mnx,then
nsel,u,,,nd
nd=node(mnx,0,0)
*endif
n22(i)=nd
*enddo
cmsel,s,current
xx2(i)=nx(nd)+0.25*lc
nsel,u,,,nd
*enddo
alls
save
===============================================================================================================
Файл VEL
!Macro for drawing velocity filed
/edge,,0
!/device,vector,on
/contour,,,0,,vel*2
/clabel,,on
plnsol,v,sum
/annot,on
/tlabel,-0.5,0.5,CL=
/tlabel,-0.4,0.5,%CL%
/tlabel,-0.5,0.4,CD=
/tlabel,-0.4,0.4,%CD%
/tlabel,-0.5,0.3,CM=
/tlabel,-0.4,0.3,%CM%
/tlabel,-0.5,0.2,ALFA=
/tlabel,-0.3,0.2,%ALFA%
!/annot,off
!/device,vector,off
===============================================================================================================
Файл PRE
!Macro for drawing pressure filed
/edge,,0
/contour,,,-10000,,10000
/clabel,,on
plnsol,pres
/annot,on
/tlabel,-0.5,0.5,CL=
/tlabel,-0.4,0.5,%CL%
/tlabel,-0.5,0.4,CD=
/tlabel,-0.4,0.4,%CD%
/tlabel,-0.5,0.3,CM=
/tlabel,-0.4,0.3,%CM%
/tlabel,-0.5,0.2,ALFA=
/tlabel,-0.3,0.2,%ALFA%
!/annot,off
===============================================================================================================
Файл CP
/post1
*do,i,1,num1
*get,pc,node,n11(i),pcoe
pw1(i)=-pc
*enddo
/yrange,-1.3,3.6
/xrange,0,lc
/axlab,y,"neg pres "
/axlab,x,"chord"
/color,curve,gree,1
*vplot,xx1(1),pw1(1)
/noer
*do,i,1,num2,1
*get,pc,node,n22(i),pcoe
pw2(i)=-pc
*enddo
/color,curve,red ,1
*vplot,xx2(1),pw2(1)
/annot,on
/tlabel,-0.5,0.5,CL=
/tlabel,-0.4,0.5,%CL%
/tlabel,-0.5,0.4,CD=
/tlabel,-0.4,0.4,%CD%
/tlabel,-0.5,0.3,CM=
/tlabel,-0.4,0.3,%CM%
/tlabel,-0.5,0.2,ALFA=
/tlabel,-0.3,0.2,%ALFA%
/annot,off
/erase
alls
5
Представительство CAD FEM GmbH в СНГ: 107497 Москва, Щелковское шоссе, 77, офис 1703Тел.: (095) 468-8175 Тел.\факс: (095) 913-2300 e-mail: cadfem@online.ru http://www.cadfem.ru
![](https://s.studizba.com/z.php?f=/templates/si/i/noavatar.png&w=100&h=100&t=1"){kind=avatar}
