Разделы сайта

Читаемое

Обновления Oct-2018

Промышленность Ижоры -->  Динамика жидкости: уравнения 

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 [ 51 ] 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182

1 SUBROUTINE BODY(N,H,IPR.FMN,B)

3 С CALCULATES BODY AND CONTROL POINT COORDINATES

4 с FOR AN ELLIPSE WITH MINOR SEMI-AXIS, В

6 DIMENSION X(50),Y(50)ДС(50),YC(50),DS(50),FN{50,50)

7 1 .FT{50,50),RHS(50),SDE(50),CI(50),SI(50)

8 COMMON X,Y,XC,YC,DS,FN,FT,RHS,PI.CPI,CI,SI

9 1 ,UINF,VINF,SDE

10 С

11 С BODY POINTS

12 С

13 FAC = SQRTd.O - FMN*FMN)

14 NHLFF = N/2 + 1

15 NHH=:NHLFF1

16 AN = NHLFF - 1

17 DTH = PI/AN

18 С

19 DO 2 I=1,HLFF

20 AI = I - 1

21 TH = AI*DTH

22 TH =r PI - TH

23 X(I) = COS(TH)

24 Yd) = B*SIN(TH)

25 С

26 С PRANDTL-GLAUERT TRANSFORMATION

27 С

28 2 Yd) = Y(I)*FAC

29 С

30 С REFLECT FOR COORDINATES OF LOWER HALF.

31 С

32 DO 3 I=NHH,N

33 X(I)=X(N+2-I)

34 3 Yd) = - Y(N+2-I)

35 X(M)=Xd)

36 Y(M)=Yd)

37 С

38 С PLACE CONTROL POINTS AT THE CENTER OF PANELS.

39 С

40 DO 4 1=1.N

41 XCd) (Xd)+Xd+l))*0.5

42 4 YCd) = (Yd)+Yd+l))*0.5-

43 С

44 С CALCULATE PANEL SPANS,COS AND SINE OF ANGLES.

45 С

46 DO 5 1=1,N

47 SX=Xd+l)-Xd)

48 SY=Yd+l)-Yd)

49 DSd)=SQRT(SX*SX+SY*SY)

50 CId) = (X(I+l)-Xd))/DS(I)

51 5 SI(I) = (Y(I+l)-Yd))/DSd)

52 IFdPR .EQ. 0)RETURN

53 С

54 С OUTPUT ELEMENT PARAMETERS

55 С

56 VRITE(6,6)

57 6 FORMAT{2X,ELEMENT PARAMETERS)

58 DO 7 I = l.N

59 7 WRITE(6,8)I.Xd).Yd) ,XC{I),YC(I),DSd).CI(I),SI(I)

60 8 FORMAT(2X,I=,12 X,Y=,2F8.4, XC.YC= ,2F8.4. SPA

61 1 F8.4. CI.SI= ,2F8.4)

62 RETURN

63 END

Рис. 14.8. Распечатка программы BODY.



1 SUBROUTINE M.ATELM(N,M,IPR)

3 С CALCULATES MATRIX ELEMENTS AND RHS.

5 DIMENSION X(50),Y(50),XC(50),YC(50),DS(50),FN{50,50)

6 1 .FT{50,50),RHS(50),SDE(50),CI(50).SI(50)

7 COMMON X,Y,XC,YC,DS,FN,FT,RHS,PI,CPI,CI,SI

8 1 ,UINF,VINF,SDE

10 DO 2 K=1,N

11 DO I J=1,N

12 IFCK.EQ.J) FN(K.J)=2.*PI

13 IF(K.EQ.J) FT(K.J)=0.0

14 IF(K.EQ.J) GO TO 1

15 DYJ=SI(J)*DS(J)

16 DXJ=CI(J)*DS{J)

17 SPH=DS(J)*0.5

18 XD=XC{K)-XC(J)

19 YD=YC(K)-YC(J)

20 RKJ=SQRT(XD*XD+YD*YD)

21 BKJ=ATAN2{YD,XD)

22 ALJ=:ATAN2(DYJ.DXJ)

23 GKJ=ALJ-BKJ

24 ZIK=RKJ*COS(GKJ)

25 ETK=-RKJ*SIN{GKJ)

26 R1S=({ZIK+SPH)**2)+ETK*ETK

27 R2S=((ZIK-SPH)**2)+ETK*ETK

28 QT AL0G(R1S/R2S)

29 DEN=ZIK*ZIK+ETK*ETK-SPH*SPH

30 GNM=ETK*DS(J)

31 QN=2.0*ATAN2(GNM.DEN)

32 UKJ=QT*CI(J)-QN*SI(J)

33 VKJ=QT*SI(J)+QN*CI{J)

34 FN(K,J)=-UKJ*SI(K)+VKJ*CI(K)

35 FT{K,J)=UKJ*CI(K)+VKJ*SI(K)

36 1 CONTINUE

37 RHS(K)=UINF*SI(K)-VINF*CI(K)

38 2 CONTINUE

39 С

40 IFCIPR .LE. DRETURN

41 3 WRITE(6.4)

42 4 FORMAT(2X,MATRIX ELEMENTS = NORMAL VELOCITY COMPONENTS

43 DO 5 K=1,N

44 5 WRITE(6,8) K,{FN(K,J).J=1,N)

45 WRITE(6.6)

46 6 FORMAT(2X,TANGENTIAL VELOCITY COMPONENTS)

47 DO 7 K=1,N

48 7 WRITE(6.8)K.(FT(K,J),J=1,N)

49 8 FORMAT(2X,I5,(10F10.5))

50 WRITE{6,9)

51 9 F0RMAT(2X.RIGHT HAND SIDE)

52 WRITE(6,10)(RHS(K),K=1,N)

53 10 FORMAT(2X,10F10.5)

54 RETURN

55 END

Рис. 14.9. Распечатка программы MATELM.



1 5UBR0UTIME SURVL(N,B,rHN)

3 С CALCULATES VELOCITIES AND PRESSURE AT TE CONTROL POINTS

4 С gEX IS THE EXACT VELOCITY AT THE SURFAOS OF THE ELLIPSE

6 DIMENSION X(50),Y(50)ДС(50)ЛС(50) DS($0),FN(50,5O)

7 1,FT(50,50) ,RHS(50) ,SDE(50).CI(50),SI(50)f

8 COMMON X,Y,XC,YC,DS,FN,rT,RHS,PI,CPI,CI,SI

9 1,UINF,VINF,SDE

10 С

11 TAC = SQRTd. - FMN*F)UI)

12 GAM = 1.4

13 CI = 0.5*(GAM-1.)*FHN*F4M

14 C2 0.5*GAM*FMN*FMM

15 CMP GAN/(GAN-1.)

16 WRITE(6,1) .

17 1 FORMAT(2X,VELOCITY AND PRESSURE AT THE tONTROL POINTS)

18 DO 4 K 1;N

19 QTS 0.0

20 QNS=0.0

21 DO 2 J=1,N

22 QTS=QTS+FT(K,a)*SDE(J)

23 2 QNS=QNS+FN(K,J)*SDE(J)

24 С

25 QNK QNS + VINF*CI(K) - UINF*SI(K)

26 QTK QTS + VINF*SI(K) + UINF*CI(K)

27 UU UINF-QNS*SI(K)+QTS*CI(K)

28 VV VINF+QNS*CI(K)+QTS*SI(K)

29 UU = UU/FAC/FAC

30 VY = VV/FAC

31 PP 1.-UU*UU-YV*VV

32 IF(FMM .GT. 0.05)PP = ((l.+Cl*PP)**CHP-l.)/C2

33 С

34 DUM = B*B*XC{K)

35 DUM - YC(K)*YC(K) + DUM*DUK

36 QEX = (1. + B)*YC(K)/SQRT(DUH)

37 С

38 WRITE(6,3)XC(K),YC(K),gNK,W.UU.VV,PP,QE3C

39 3 F0RMAT(1X/XC,YC=\2F6.3,* QH,QT ,2F€ Л,

40 1 U,V=,2F6.3, P ?,F6.3,* QBX ,F6.3)

41 4 CONTINUE

42 RETURN

43 END

Рис. 14.10. Распечатка программы SURVL.

Ординаты точек (хрур) считываются из файла входных данных, расположенного на устройстве с логическим номером L

Параметры, используемые в подпрограмме PANEL, описаны в табл. 14.2, а на рис. 14.12 приведена типичная выдача результатов расчета обтекания эллипса, длина меньшей полуоси которого ft = 0.5. На рис. 14.12 в каждой контрольной точке приведены значения нормальной QN и тангенциальной QT относительно локального наклона тела составляющих скорости. Нормальная составляющая QN равна нулю. Можно напомнить, что это граничное условие (14.4) используется для определения плотности источников. Значение касательной составляющей QT



1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 [ 51 ] 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182

© 2003 - 2018 Prom Izhora
При копировании текстов приветствуется обратная ссылка