Top / AA_Experiment / Rocket / SampleProgram_f

AA_Experiment/Rocket/SampleProgram_f

_ 圧力比 pe/pc の数値解をニュートン法で求めるサンプルプログラム (Fortran77言語版)

下記の内容で、拡張子".f"または".f77"のファイルを作成し、Fortran77コンパイラでコンパイルし、実行する。

ソースファイルをダウンロード → filenewton.f

C ***********************************************
C * 航空宇宙工学実験「ロケット推進力の測定」    *
C * 与えられたノズル開口比 Ae/At から、         *
C * 圧力比 pe/pc の数値解をニュートン法で求める *
C ***********************************************
C
      PROGRAM NEWTON
C
      REAL PI,EPS,G,DT,DE,AT,AE,AR,PR,F,DF,DX
      INTEGER I
      I=0
      PI=3.14159265
      EPS=1.0E-10
      G=1.402
      DT=2.5
      DE=3.3
      AT=DT*DT*PI/4.0
      AE=DE*DE*PI/4.0
      AR=DE*DE/DT/DT
      WRITE(*,*) 'スロート面積 At =',AT
      WRITE(*,*) '  出口面積 Ae =',AE
      WRITE(*,*) '開口比  Ae/At =',AR
      PR=0.1
      DX=1.0
      WRITE(*,*)
      WRITE(*,*) 'i   pe/pc       f           df'
      DO WHILE (ABS(DX).GT.EPS)
        F=PR**(2.0/G)-PR**(1.0+1.0/G)
     *    -((G-1.0)/2.0*(2.0/(G+1.0))**((G+1.0)/(G-1.0)))/AR/AR
        DF=(2.0/G)*PR**(2.0/G-1.0)-(1.0+1.0/G)*PR**(1.0/G)
        WRITE(*,'(I2,3F12.8)') I,PR,F,DF
        DX=-F/DF
        PR=PR+DX
        I=I+1
      END DO
      WRITE(*,*)
      WRITE(*,*) '圧力比 pe/pc = ',PR
      END
 
Attach file: filenewton.f 874 download [Information]
 
Link: AA_Experiment/Rocket(3465d)
Last-modified: 2005-04-19 (Tue) 12:15:02 (4596d)