Modelica® - A Unified Object-Oriented Language for Systems Modeling Language Specification Version 3.5-dev

Appendix C Modelica DAE Representation

In this appendix, the mapping of a Modelica model into an appropriate mathematical description form is discussed.

In a first step, a Modelica translator transforms a hierarchical Modelica simulation model into a “flat” set of Modelica “statements”, consisting of the equation and algorithm sections of all used components by:

  • expanding all class definitions (flattening the inheritance tree) and adding the equations and assignment statements of the expanded classes for every instance of the model.

  • replacing all connect-equations by the corresponding equations of the connection set (see section 9.2).

  • mapping all algorithm sections to equation sets.

  • mapping all when-clauses to equation sets (see section 8.3.5).

As a result of this transformation process, a set of equations is obtained consisting of differential, algebraic and discrete equations of the following form where (v:=[x˙;x;y;t;m;pre(m);p] ):

c:=fc(relation(v)) (C.1a)
m˙:=fm(v,v) (C.1b)
0=fx(v,v) (C.1c)

and where

p Modelica variables declared as parameter or constant, i.e., variables without any time-dependency.
t Modelica variable time, the independent (real) variable.
x(t) Modelica variables of type Real, appearing differentiated.
m(te) Modelica variables of type discrete Real, Boolean, Integer which are unknown. These variables change their value only at event instants te. pre(m) are the values of m immediately before the current event occurred.
y(t) Modelica variables of type Real which do not fall into any other category (= algebraic variables).
c(te) The conditions of all if-expressions generated including when-clauses after conversion, see section 8.3.5).
relation(v) A relation containing variables vi, e.g. v1 > v2, v3 >= 0.

For simplicity, the special cases of the noEvent() operator and of the reinit() operator are not contained in the equations above and are not discussed below.

The generated set of equations is used for simulation and other analysis activities. Simulation means that an initial value problem is solved, i.e., initial values have to be provided for the states x, section 8.6. The equations define a DAE (Differential Algebraic Equations) which may have discontinuities, a variable structure and/or which are controlled by a discrete-event system. Such types of systems are called hybrid DAEs. Simulation is performed in the following way:

  1. 1.

    The DAE (C.1c) is solved by a numerical integration method. In this phase the conditions c of the if- and when-clauses, as well as the discrete variables m are kept constant. Therefore, (C.1c) is a continuous function of continuous variables and the most basic requirement of numerical integrators is fulfilled.

  2. 2.

    During integration, all relations from (C.1a) are monitored. If one of the relations changes its value an event is triggered, i.e., the exact time instant of the change is determined and the integration is halted. As discussed in section 8.5, relations which depend only on time are usually treated in a special way, because this allows to determine the time instant of the next event in advance.

  3. 3.

    At an event instant, (C.1) is a mixed set of algebraic equations which is solved for the Real, Boolean and Integer unknowns.

  4. 4.

    After an event is processed, the integration is restarted with 1.

Note, that both the values of the conditions c as well as the values of m (all discrete Real, Boolean and Integer variables) are only changed at an event instant and that these variables remain constant during continuous integration. At every event instant, new values of the discrete variables m and of new initial values for the states x are determined. The change of discrete variables may characterize a new structure of a DAE where elements of the state vector x are disabled. In other words, the number of state variables, algebraic variables and residue equations of a DAE may change at event instants by disabling the appropriate part of the DAE. For clarity of the equations, this is not explicitly shown by an additional index in (C.1).

At an event instant, including the initial event, the model equations are reinitialized according to the following iteration procedure:

    known  variables: x, t, p
    unkown variables: dx/dt, y, m, pre(m), c
   // pre(m) = value of m before event occured
   loop
      solve (1) for the unknowns, with pre(m) fixed
      if m == pre(m) then break
      pre(m) := m
   end loop

Solving (C.1) for the unknowns is non-trivial, because this set of equations contains not only Real, but also Boolean and Integer unknowns. Usually, in a first step these equations are sorted and in many cases the Boolean and Integer unknowns can be just computed by a forward evaluation sequence. In some cases, there remain systems of equations (e.g. for ideal diodes, Coulomb friction elements) and specialized algorithms have to be used to solve them.

Due to the construction of the equations by ”flattening” a Modelica model, the hybrid DAE (C.1) contains a huge number of sparse equations. Therefore, direct simulation of (C.1) requires sparse matrix methods. However, solving this initial set of equations directly with a numerical method is both unreliable and inefficient. One reason is that many Modelica models, like the mechanical ones, have a DAE index of 2 or 3, i.e., the overall number of states of the model is less than the sum of the states of the sub-components. In such a case, every direct numerical method has the difficulty that the numerical condition becomes worse, if the integrator step size is reduced and that a step size of zero leads to a singularity. Another problem is the handling of idealized elements, such as ideal diodes or Coulomb friction. These elements lead to mixed systems of equations having both Real and Boolean unknowns. Specialized algorithms are needed to solve such systems.

To summarize, symbolic transformation techniques are needed to transform (1) in a set of equations which can be numerically solved reliably. Most important, the algorithm of Pantelides should to be applied to differentiate certain parts of the equations in order to reduce the index. Note, that also explicit integration methods, such as Runge-Kutta algorithms, can be used to solve (C.1c), after the index of (C.1c) has been reduced by the Pantelides algorithm: During continuous integration, the integrator provides x and t. Then, (C.1c) is a linear or nonlinear system of equations to compute the algebraic variables y and the state derivatives dx/dt and the model returns dx/dt to the integrator by solving these systems of equations. Often, (C.1c) is just a linear system of equations in these unknowns, so that the solution is straightforward. This procedure is especially useful for real-time simulation where usually explicit one-step methods are used.