Разделы сайта
Читаемое
Обновления Mar-2024
|
Промышленность Ижоры --> Динамика жидкости: уравнения 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
|
© 2003 - 2024 Prom Izhora
При копировании текстов приветствуется обратная ссылка |