Д6853 Ширяев ЮН Расчет горизонтального кожухотруб конденсатора

( 0 4

( . ..* + ,) *, 90

«GKTkond.exe». - -

Fortran 90.

-

.

!> 6

parameter (pi=3.14) parameter (g=9.81) real q0,tn,tg1,tg2,wg real dn,dv,dl

real dr,tr,sr

real lg,ng,bg,prg real ra,roa,la,na real en

real no

real taum,tgs,reg,grg,nug,nul,nut,gamma,alg,fv,fo,fn,k1 real sg,sv,lo,b,taua,a,ko,qo,f,l,ntr,m1,d,ltr,nsr,m

print*,'% [!]' read*,tg1

print*,'% [!]' read*,tg2

23

print*,'! [ / ]' read*,wg

!>

taum=(tg2-tg1)/log((tn-tg1)/(tn-tg2)) tgs=tn-taum

!- %&!

print*,'- %&! :' print*,'" ',tgs print*,'%'

read*,lg

print*,'. ' read*,ng

print*,'. ,' read*,bg

print*,' "' read*,prg

!- %&!

print*,'- %&! :' print*,'" ',tn print*,'%' read*,la

print*,'. ' read*,na

print*,' ' read*,ra

print*,'" + ' read*,roa

!-

2 print*,'- :' print*,'1 - ;'

print*,'2 - ;' print*,'3 - .'

24

read*,tip

!-

print*,'- :' print*,'1 - (384)'

print*,'2 - (45)' read*,i

if (i==1) ltr=384 if (i==2) ltr=45

3 print*,'- :' print*,'$ ' read*,dn

print*,'- ' read*,dv

print*,'% ' read*,dl

!-

if (tip==3) then

print*,'- . ' read*,beta

end if

! Re reg=wg*dv/ng ! Nu

if (reg>=10000) nug=0.021*reg**0.8*prg**0.43 if (reg<=2300) then

grg=g*(dv**3)*bg*0.5*taum/(ng**2) if ((grg*prg)<800000)

nug=1.55*(reg*prg/100)**(1/3)*(1+0.01*(reg/100)**(2/3))

if ((grg*prg)>=800000) nug=0.15*reg**0.33*prg**0.43*grg**0.1 end if

if ((reg>2300).and.(reg<10000)) then grg=g*(dv**3)*bg*0.5*taum/(ng**2)

25

if ((grg*prg)<800000) nul=1.55*(2300*prg/100)**(1/3)*(1+0.01*(2300/100)**(2/3))

if ((grg*prg)>=800000) nul=0.15*2300**0.33*prg**0.43*grg**0.1 nut=0.021*10000**0.8*prg**0.43

gamma=1-exp(1-reg/2300) nug=gamma*nut+(1-gamma)*nul

end if

! .

alg=nug*lg/dv

!-

4 if (tip==2) then

print*,'- :' print*,' '

read*,dr

print*,'% + ' read*,tr

print*,'C ' read*,sr

print*,'% ' read*,lr

end if

! +

fv=pi*dv

fo=pi*dn

if (tip==1) fn=fo

if (tip==2) fn=(pi*dn*(sr-0.5*(dn+dv))/sr)+pi/sr*((dr**2- dn**2)/2+dr*0.5*(dn+dv))

if (tip==3) fn=fo*beta

! .

k1=1/(1/alg*fo/fv+0.5*(dn-dv)/ltr*fo/(0.5*(dn+dv))+dl*fo/fv) !-

5 print*,'- :'

26

print*,'" ,' read*,sg print*,'" ,' read*,sv

print*,' ' read*,lo

! .

b=0.728*(g*ra*roa*la**3/na/dn)**0.25 !

t1=0.1*taum

t2=0.9*taum

20 taua=(t2+t1)/2 if (tip==1) then

alagl=b*taua**(-0.25) qa=b*taua**0.75 fot=q0/qa ntr=fot/pi/dn/lo ntr=aint(ntr)+1

nsr=0.92*ntr**0.5*(sg/sv)**0.5

nsr=aint(nsr)+1 a=b*nsr**(-1/6)/k1

end if

if (tip==2) then alagl=b*taua**(-0.25) m=(2*alagl/lr/tr)**0.5 h1=0.5*(dr-dn)*(1+0.35*log(dr/dn))

thmh=(exp(m*h1)-exp(-m*h1))/(exp(m*h1)+exp(-m*h1)) e=thmh/m/h1

hpr=0.785*(dr**2-dn**2)/dr fve=pi*(dr**2-dn**2)/2/sr fgo=pi*dn*(1-tr/sr)+pi*dr*tr/sr

27

psi=1.1*(fve/fo)*e**0.75+(dn/hpr)*fgo/fo

qa=b*psi*taua**0.75

fot=q0/qa

ntr=fot/pi/dn/lo

ntr=aint(ntr)+1

nsr=0.92*ntr**0.5*(sg/sv)**0.5

nsr=aint(nsr)+1 a=b*nsr**(-1/6)*psi/k1

end if

if (tip==3) then alagl=b*taua**(-0.25) qa=b*beta*taua**0.75 fot=q0/qa ntr=fot/pi/dn/lo ntr=aint(ntr)+1

nsr=0.92*ntr**0.5*(sg/sv)**0.5

nsr=aint(nsr)+1 a=b*beta*nsr**(-1/6)/k1

end if

if ((taum-(taua+a*taua**0.75))>0.1) then t1=taua

goto 20 end if

if ((taum-(taua+a*taua**0.75))<-0.1) then t2=taua

goto 20 end if

ala=a*k1*taua**(-0.25)

! .

ko=1/(1/ala+1/k1)

!

qo=ko*taum

28

! +

f=q0/qo

!

l=f/pi/dn

! + ( ) ntr=L/lo

ntr=aint(ntr)+1

!D ,

m1=((ntr-0.25)/0.75)**0.5 m1=anint(m1)

!

d=m1*sg

!-

print*,'.. . ',ko*fo/fv print*,'.. . ',ko*fo/fn print*,' '

print*,'" + ( )-',f print*,'" + ( .)--',f*fv/fo print*,'" + ( .)-',f*fn/fo

if (tip==2) print*,'5 - : ' read*,i

if (i==1) goto 1 if (i==2) goto 2 if (i==3) goto 3 if (i==4) goto 5 if (i==5) goto 4 end

30

% !

31

3 4 1 .

' .

#

!"# "$ % # % & " '

* . .

.0 1

/ 2 1

. / 1 / .,

. 1 / . . 2

1 0

@.#. @

! $.7.

!

$.-. #

____________________________________________________________________

" 25.12.2009. & 60×84 1/16

. . . 1,86. ". . 2,0. .- . . 1,75

% 250 .. : H C 86

_____________________________________________________________________

!" #$"%. 191002, !-", . @, 9 77 !" #$"%. 191002, !-", . @, 9

32