rumoca 0.7.28

Modelica compiler written in RUST
Documentation 
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

// Simplified versions of Modelica.Blocks.Nonlinear blocks for balance testing
package NonlinearBlocks
  // SISO interface base
  block SISO
    input Real u;
    output Real y;
  end SISO;

  // Limiter: extends SISO, adds 1 protected variable (simplifiedExpr)
  // Should be: 2 unknowns (y, simplifiedExpr), 2 equations
  block Limiter
    extends SISO;
    parameter Real uMax = 1;
    parameter Real uMin = -uMax;
    parameter Boolean strict = false;
    Real simplifiedExpr;
  equation
    simplifiedExpr = u;
    // simplified binding
    if strict then
      y = smooth(0, noEvent(if u > uMax then uMax else if u < uMin then uMin else u));
    else
      y = smooth(0, if u > uMax then uMax else if u < uMin then uMin else u);
    end if;
  end Limiter;

  // VariableLimiter: extends SISO, adds 2 inputs and 1 protected
  // Should be: 2 unknowns (y, simplifiedExpr), 2 equations
  block VariableLimiter
    extends SISO;
    input Real limit1;
    input Real limit2;
    parameter Boolean strict = false;
    Real simplifiedExpr;
  equation
    simplifiedExpr = u;
    if strict then
      y = smooth(0, noEvent(if u > limit1 then limit1 else if u < limit2 then limit2 else u));
    else
      y = smooth(0, if u > limit1 then limit1 else if u < limit2 then limit2 else u);
    end if;
  end VariableLimiter;

  // SlewRateLimiter: extends SISO, has state y and protected val
  // Should be: 2 unknowns (y state, val algebraic), 2 equations (der(y)=..., val=...)
  block SlewRateLimiter
    extends SISO;
    parameter Real Rising = 1;
    parameter Real Falling = -Rising;
    parameter Real Td = 0.001;
    parameter Boolean strict = false;
    Real val = (u - y) / Td; // binding equation
  equation
    if strict then
      der(y) = smooth(1, if noEvent(val < Falling) then Falling else if noEvent(val > Rising) then Rising else val);
    else
      der(y) = if val < Falling then Falling else if val > Rising then Rising else val;
    end if;
  end SlewRateLimiter;

  // DeadZone: extends SISO
  // Should be: 1 unknown (y), 1 equation
  block DeadZone
    extends SISO;
    parameter Real uMax = 1;
    parameter Real uMin = -uMax;
  equation
    y = smooth(0, if u > uMax then u - uMax else if u < uMin then u - uMin else 0);
  end DeadZone;
end NonlinearBlocks;