シミュレーションの演習のコード

GSA1のコードを以下に示します。このコードは BGSP用です。第１行の 

BOND SYMBOL EXPR OUT DEFIN PARM INIT

はコードの一部ですから省略してはいけません。また、BGSPではカラム指定が厳格で1カラムでもずれるとエラーになります。また、許される文字は英数字、コロン(:)、スラッシュ(/)、加減乗算記号(+-*)、半角空白のみであり、制御文字が入るとエラーになります。倍角の空白はエラーを生じますが、ソースをプリントしても気づきにくいので注意が必要です。
方法1と方法2に分けて記述します。方法１はダミーボンドを使用しますが、方法2では使用しません。

GSA1のコード

方法1のコード

BOND SYMBOL EXPR OUT DEFIN PARM INIT
SFP    1
ZJ     -1    2     3
RER    2
DE1    3
DF1    10
OJ     -10   4     5
REB    4
ZJ     -5    6     7
CIB    6
OJ     -7    8     9
REV    8
SE0    9
/
BC: ELRV  EN2
OP: ELFL
FB: X11    EN10
FB: X12    EN7
      OP=X11-X12
EOD
BC: P-PR  EF2
BC: F-OUT FL8
/
SFP
PA: PFPUMP 5.0D-5         FLOW OF PUMP
       F=PFPUMP
EOD
SE0
PA: PPOUT  0.0D-00         PRESSURE OUTSIDE
       E=PPOUT
EOD
RER
PA: PPTHR  2.8D07         THRESHOLD PRESSURE
PA: PREL   4.1D09          RESISTANCE OF RELIEF VALVE
        IF(Z.LT.PPTHR) THEN
          R=0.0D-00
        ELSE
          R=(Z-PPTHR)/PREL
        ENDIF
EOD
REV
PA: PRREST 8.9D09          RESISTANCE OF RESTRICTION
         R=ABS(Z)**0.5D-00/PRREST
EOD
REB
PA: PPLOSS 0.91D09          RESITANCE OF LINE
         R=PPLOSS*Z
EOD
CIB
PA: PVB    3.3D-12         VOLUME/BULK MODULOUS
         C=Z/PVB
EOD
/
0.0D-00            5.0D-00           0.01D-00          500
/

方法2のコード

BOND SYMBOL EXPR OUT DEFIN PARM INIT
SFP    1
ZJ     -1    2     3     4 
RER    2
CIA    3
OJ     -4    5     6
REB    5
ZJ     -6    7     8
CIB    7
OJ     -8    9     10
REV    9
SE0    10
/
BC: ELRV  EN2
OP: ELFL
FB: X11    EN4
FB: X12    EN8
      OP=X11-X12
EOD
BC: P-PR  EF2
BC: F-OUT FL9
/
SFP
PA: PFPUMP 5.0D-05
       F=PFPUMP
EOD
CIA
PA: PAC1   1.0D-13
CIA
       C=Z/PAC1
EOD
SE0
PA: PPOUT  0.0D-00
       E=PPOUT
EOD
RER
PA: PPTHR  2.8D+07
PA: PREL   4.1D+09
        IF(Z.LT.PPTHR)THEN
         R=0.0D-00
        ELSE
          R=(Z-PPTHR)/PREL
        ENDIF
EOD
REV
PA: PRREST 8.9D+09         RESISTANCE OF RESTRICTION
         R=ABS(Z)**0.5/PRREST
EOD
REB
PA: PPLOSS 0.91D+09        RESITANCE OF LINE
         R=Z/PPLOSS
EOD
CIB
PA: PVB    3.3D-12         VOLUME/BULK MODULOUS
         C=Z/PVB
EOD
/
0.0D-00           5.0D-00           0.01D-05          500
/