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

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

BOND SYMBOL EXPR OUT DEFIN PARM INIT

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

GS4のコード

BOND SYMBOL EXPR OUT DEFIN PARM INIT
SE1    21
OJ     -21   20    23 
II1    23
ZJ     -20   25    24
OJ     -25   33    34
RE1    24
RE2    34
II2    33
/
BC: Y1    EF24      伝達トルク
BC: Y2    FL23
BC: Y3    FL33
/
SE1  
PA: PX1    20.0            Torque Source
      E=PX1
EOD
II1                        
PA: PI1    0.30            moment of inertia
      L=Z/PI1
EOD
RE1
PA: PR1    1.0             Basic Torque
      IF(T.LT.2.0) THEN
       F1=0.0
       ELSE
       IF(T.LT.3.0) THEN
        F1=T-2.0
       ELSE
        IF(T.LT.6.0) THEN
         F1=1.0
        ELSE
         IF(T.LT.6.2) THEN
          F1=31.0-5.0*T
         ELSE
          F1=0.0
          ENDIF
         ENDIF
        ENDIF
       ENDIF
       R=Z*PR1*F1
EOD
RE2
PA: PR2    0.02            Load Friction
      R=PR2*Z
EOD
II2
PA: PI2    0.4             Moment of Inertia
      L=Z/PI2
EOD
/
 0.0D+00           10.0             0.05              200 
/

20-Simによるコード



I23

parameters
	real i = 0.3;
equations
    state = int(p.e);
    p.f = state / i;

I33
parameters
	real i = 0.4;
equations
    state = int(p.e);
    p.f = state / i;

MR241
equations
	p.e = r * p.f;

R34
parameters
	real r = 0.02;
equations
	p.e = r * p.f;

Se21
parameters
	real effort = 20.0;
variables
	real flow;
equations
	p.e = effort;
	flow = p.f;

Table242
parameters					
	string filename= 'd:\20-sim\Simtec\GS4.txt';
equations					
	output = table (filename, input);

Time243
equations
	output = time;
     

GS4.txt

0.0  0.0
2.0  0.0
3.0  1.0
6.0  1.0
6.2  0.0
10.0 0.0