Difference between revisions of "Simantics System Dynamics"

From End-User Wiki
Jump to: navigation, search
(Arithmetic Operators)
(Components)
Line 20: Line 20:
 
[[File:ComponentTypes.png|right|frame|Component types]]
 
[[File:ComponentTypes.png|right|frame|Component types]]
 
;Auxiliary
 
;Auxiliary
:Auxiliary is the most basic variable you can use. It represents a single value or a mathematical expression. There are different types of auxiliary variables currently in the system: Auxiliary, Parameter, Constant and WithLookup. Auxiliary is the default type. Parameters are single values that the user can change. When only parameters are changed, the model simulates faster, because the system does not have to recompile the model. WithLookup is a special variable that has an expression and a lookup table. The expresssion defines what value is taken from the defined table.
+
:Auxiliary is the most basic variable you can use. It represents a single value or a mathematical expression. There are different types of auxiliary variables currently in the system: Auxiliary, Parameter, Constant, Delay and WithLookup. Auxiliary is the default type. Parameters are single values that the user can change. When only parameters are changed, the model simulates faster, because the system does not have to recompile the model. WithLookup and Delay are special variables. WithLookup has an expression and a lookup table. The expresssion defines what value is taken from the defined table. Delay variable delays an equation a given time with a given delay order.
 
;Dependency
 
;Dependency
 
:Dependency is an arrow that connects two components. It means that the value of the variable from which the arrow starts is used to calculate the value of the variable where the arrow ends.
 
:Dependency is an arrow that connects two components. It means that the value of the variable from which the arrow starts is used to calculate the value of the variable where the arrow ends.

Revision as of 10:37, 25 October 2012

Simantics System Dynamics

Simantics System Dynamics is currently the only open source modelling and simulating tool for Simantics. Simantics System Dynamics is under development and will go through some changes in the future. New features will be added and old ones improved according to the needs of the modellers. New features include module libraries and different analysis tools.

This documentation introduces you to the current version of Simantics System Dynamics. Documentation includes the basic modelling principles and a guide on how to model system dynamics models with Simantics System Dynamics. If you like to get to know the tool better and try modelling and simulating yourself, install the software and try our basic and advanced tutorials!

Introduction to System Dynamics Simulation

System Dynamics

System dynamics is an approach to understanding different organizations, markets and other complex systems and their dynamic behavior. System dynamics modelling in Simantics is a free modelling tool. See installation instructions.

Model

Model structure

System dynamics model is generally understood as the model configuration. In this tool, the model contains also other components. Model can have Operating interfaces that enable the use of the model with simple controls. Operating interfaces are configured separately and when using them, the complex model is hidden from the user. Experiments are the way to simulate the model. In the future, you can have experiments with different configurations, for example different initial values for some parameters. In that way, you don't have to always configure the model for different scenarios. Modules folder contains all the different module types in your model and you can create new module types there.


Components

Most of the components you can use in your models are basic system dynamics components. The modularity of the models introduces two additional components, Modules and Inputs. All the components are explained below.

Component types
Auxiliary
Auxiliary is the most basic variable you can use. It represents a single value or a mathematical expression. There are different types of auxiliary variables currently in the system: Auxiliary, Parameter, Constant, Delay and WithLookup. Auxiliary is the default type. Parameters are single values that the user can change. When only parameters are changed, the model simulates faster, because the system does not have to recompile the model. WithLookup and Delay are special variables. WithLookup has an expression and a lookup table. The expresssion defines what value is taken from the defined table. Delay variable delays an equation a given time with a given delay order.
Dependency
Dependency is an arrow that connects two components. It means that the value of the variable from which the arrow starts is used to calculate the value of the variable where the arrow ends.
Flow
Flow connects clouds, valves and stocks. Flow represents an actual flow of something from stocks or clouds to stocks or clouds. There has to be at least one valve in a flow and the system creates it automatically, if none of the ends of the flow is a valve.
Valve
Valve regulates the rate of a flow. The value of a valve is automatically used in calculating the level of an adjacent stock. Valves behave just like Auxiliary variables but look different and you can connect also flows to them.
Stock
The value of a stock variable is an integral of flows leaving and flows arriving to the variable. The integral is calculated automatically from the valves that are connected to the variable with flow connections. You must give an initial value to the stock. Initial value can be a single value or an equation. You can use values of other variables to calculate the initial value.
Cloud
Cloud is not a variable. It represents a starting or ending point of a flow, if it is not in the scope of the model.
Module
Modules enable structural modeling. Modules are defined just like the basic model configuration, but the module component hides the actual configuration. You can only connect dependency connections into the module and dependency connections from the module must end to Input variables. The interface of the module is defined using input and output variables in the configuration of the module. All variable types can be set as output variables. If a variable is an output variable, its font is bold.
Input
Input variables are the way of getting values from other modules. Inputs look like auxiliary variables except their font is italics. You can set a default value to the input variable in case it is not connected to any variable. Connections are made from the modules properties, when the module is populated. Input doesn't have to be connected with an arrow to a module. If the variable has no connections, it can get values from a higher level in the hierarchy. If an input is connected to an output, the output and its module are shown below the variable using dot notation.


Modeling principles

System dynamics modelling is much more than just mathematical formulas and nice graphs. Models are ways of communicating. There are some basic principles in system dynamics modelling that make the models easier to read and understand. You do not have to apply these principles to simulate models, but using them makes it easier for you to communicate your model to others.

Variable names
Variables names should be nouns, not verbs. The names should be positive: for example it is easier to understand that satisfaction decreases than dissatisfaction rises.
Variable names can and should have multiple words, if they are needed. Unfortunately current version of the tool does not support spaces in variable names. A good practice is to use CamelCase if name contains multiple words. This means that you start every word woth a capital letter (example: ExampleVariable).
Connections
You should try to avoid overlapping dependency arrows. In some situations, however, it is not possible. Dependencies should also form distinctive loops, if there is a loop. It makes it easier to read and understand the model and its behavior.
Graphical annotations
System dynamics contains usually annotations for loops, polarities, delays and so on. Currently those annotations are not supported, but the support will be added later. Annotations are important, again, for communicating the behavior of the model.

Installation Instructions

System dynamics tool is a provided with the Simantics platform.

  1. Download installer from Simantics download page
  2. Install the program to a directory without spaces (limitation inherited from OpenModelica).
  3. Run the application through the launcher generated by the installer
(default: Simantics\Simantics-1.6.0-Sysdyn-32\Simantics-1.6.0-Sysdyn-32)

OpenModelica is used to build and simulate the models. Simantics platform has integrated OpenModelica 1.8.1 for Windows environments. For other versions and other environments you need to install OpenModelica.

Workbench

BasicWorkbench.png


  1. Diagram
    Diagram is the area where you will graphically modify your model. Diagrams are built from elements that can be dragged from Symbols view or populated using shortcut keys.
  2. Model Browser
    Model browser shows the structure of your model and all items related to it.
    Symbols view is stacked with the model browser and used for dragging elements to diagrams.
  3. Properties
    Property view shows the selected variable's properties. Property view has a different layout depending on the type of the selected component. The view can also have different tabs depending on the component type. Basic tabs are with variables are Equation and Additional information.
  4. Trend
    Trend view shows the graphical representation of the values of the selected value over time. For the trend to be shown, a simulation has to be run. The view shows always the results of the latest run, but you can save results of a simulation and show them in the same trend with results from another simulation. Values view and Dependencies are stacked with Trend view and can also be used for analyzing the model.
  5. Experiment controls
    Experiment controls are shown when an experiment is active. Experiment is activated by double clicking an experiment in the model browser. With the experiment control, you can start simulation runs and save simulation results.


If you accidentally close a view, you can reopen them from Window->Show View->Other...

Modelling

Basic modelling

Basic modelling functions enable you to create and configure models. System dynamics modeling is basically pretty simple, so with these instructions you can build small and also very large models. The tricky part is writing all the expressions and adjusting the model so that it actually tells you something.

Basic 1.png
Creating a new model
Start a new model by right-clicking the model browser and selecting New->Model or from the main menu File->New Model.


Basic 2.png
Creating a new module type
Create a new module type by right-clicking on the Modules-folder and selecting New->Module. This creates a new module type that you can populate to your other modules and the model configuration.


Configuring a model
Model configuration can be opened by double clicking Configuration in the model browser.
Configuring modules
Configuration of a module type can be opened by double clicking the module type you want to configure. You can also open the configuration of a module from a diagram, when a module has been populated to that diagram, by right-clicking the module and selecting Show Module. When opening modules from diagram, the opened diagram knows to which diagram the module has been populated and can show the connections between the modules. Keep in mind that when making changes to a module, the changes apply to all instances of the module!
Populate variables
You can drag variables to a diagram from symbol view. You can also populate variables using shortcut keys
Populate modules
Modules are populated from the model browser. Just drag the module you want to populate from the Modules folder to a diagram.
Create connections
There are two types of connections: dependencies and flows. Both are created basically the same way. Hold Alt down and click on a variable. Left click starts a dependency, right click starts a flow. Both are ended to another variable with a left click.
Flows can also be started and ended to an empty spot in the diagram. If there is no start or end variable, a cloud will be created. Also if start or end is not a valve, a new valve is created in the middle of the flow.
There are some restrictions on what connections can be made, but don't worry, the user interface won't let you do connections that are not allowed.
Connections between modules
Outside the module
You can connect variables to variables in modules like this:
MODULE -----> INPUT
or
ANY VARIABLE -----> MODULE 
This is just the visual configuration, but you need those connections to really connect variables in the module's properties.
In Inputs -tab you select which variables you connect to inputs inside the module. In Outputs -tab you select which variables you lift from the module to inputs outside it.
Inside the module
Input variables get values from outside the module
Output variables can send their values outside the module. From the variable properties, select Additional Information and press Is Output
Configure variables
Select a single variable from diagram or model browser. The properties of the variables are shown in the equation view and you can modify them.
Export model
To export your model to a file, select your model from the model browser, right-click and from the context menu choose Export->Model. Select the folder where to export your model, give the file a name and press Save. You do not need to export a model to Save it, the model is automatically saved in your database. Export can be used to for example create different versions of the model, create backups or to transport it to another database.
Import model
Right-click on the model browser and select Import->Model. Browse to your .tg file and select open. The model is added to your model browser.

Shortcut and control keys

Shortcut keys for configuring a model on diagram.

Esc Cancel operations (e.g. connection and rename).
Shift + A Hover Auxiliary at the cursor position, populate with left mouse button.
Shift + S Hover Stock at the cursor position, populate with left mouse button.
Shift + C Hover Cloud at the cursor position, populate with left mouse button.
Shift + V Hover Valve at the cursor position, populate with left mouse button.
Shift + I Hover Input at the cursor position, populate with left mouse button.
Alt + left mouse button Start an arrow from a variable. End to another variable by clicking left mouse button.
Alt + right mouse button Start flow from a variable. End by clicking left mouse button.

If a flow is not started or ended on to a variable, a cloud will be created to that end.
If a new flow does not have a valve at either end, a valve will be created in the middle of the flow.

Delete Remove selected variables
Ctrl + left mouse button Select multiple variables
Mouse wheel or
+ and -
Diagram zoom
drag(mouse3) or
drag(shift + any mouse button)
Diagram pan
Ctrl + C Copy selected elements
Ctrl + X Cut selected elements
Ctrl + V Paste copied or cut elements
G Show / hide grid
R Show / hide ruler

Other shortcut keys can be found selecting Window -> Preferences from the main menu. Keys are located in General -> Keys.

Simulation

There are three different ways for simulating a model. Different simulations are represented as different types of experiments in Simantics System Dynamics.

Experiments are created from the context menu of experiments folder in a model.

NewExperiment.png

Double-click on an experiment activates it and displays experiment type specific controls in the main tool bar.

Normal simulation
Experiment is the basic simulation type. It simulates a model from start time to end time and results can be viewed after the simulation has been run.
Game simulation
Game simulations allow simulating an arbitrary number of steps, changing parameter values and continuing simulation with the new values from where it previously ended.
Simulation playback
Simulation playback works essentially as the basic experiment, but it offers some additional visualization options.

Simulation result and model structure analysis

Simulation results can be viewed as trends in trend view or as plain numbers in values view. The views display results of variables that are selected in the workbench. Users can also create custom charts.

Model structure can be examined in the Structure view. There are two options: Dependency structures of individual variables or the complete hierarchy structure of the model.

Trend
Shows you the values of the selected variable(s) graphically over the simulation time. You can modify the trend and zoom it using the context menu (right-click) of the trend.
Trend.png
Values
Shows you the values of the selected variable(s) in a table form. If the value has a star (*) at the end, it means that the value is an approximate. There is no value for that variable for the given time step. Selected values can be copied to other software in csv format.
Values.png
Compare results
You can compare different results of the same model by saving simulation results and displaying the saved results side by side with other results. You can save your results after simulating by clicking the save button Saveresults.png on your experiment controls. The saved results appear to model browser under the active experiment. To show the results on trends and tables, right-click on the result and select Show on charts.
CompareTrend.png
Dependencies
Dependencies display which variables affect and which are affected by the selected variable. The direction and number of steps that are traveled from the selected variable can be configured at the bottom of the view.
Dependencies.png
Model Hierarchy
Model Hierarchy displays the hierarchical structure of the model. The module type where a selected variable is located is highlighted.
ModelHierarchy.png

Custom charts

Custom charts are user-defined displays of simulation result data. They can be used and re-used in various places. The three types of custom charts (line, bar and pie charts) provide model developers means for examining and developing their models and communicating the simulation results to others.

Custom charts are created to the Charts folder in model browser.

Creating a new chart

Charts are configured the same way as variables: first select the chart from model browser and then configure it in the property view. Charts can be viewed in the same active trend view as any variable by selecting the chart from model browser. In addition to trend view, charts can be populated directly to the model diagram or to a separate chart panel that can contain multiple charts.

Line chart

Line char shows the time series of a variable. It is capable of displaying multiple multi-dimensional variables at the same time. Simulation playback mode adds a marker to the chart to display the current time in playback simulation. Below is an example chart and its settings.

Line chart with an array variable

General settings of line chart

General
Name The name of the chart. Displayed in e.g. model browser
Title The title of the chart. Displayed in the chart
Type Line chart (default) or area chart
Hide Options to hide chart grid, title or legend
X-axis
Variable A variable that is used as the x-axis variable. The default variable is time and it needs not to be typed in the field.
Label Optional other label for the variable than its name
Min Minimum value of the x-axis
Max Maximum value of the x-axis
Hide Options to hide x-axis lable, line, tick marks and tick labels

Axis and variables settings of line chart

Variable
Variable The full path of the variable. This field has content assistant, which shows the possible variables as you type.
Range If the variable has multiple dimensions, this field allows to select, which dimensions are displayed. It is also possible to sum the variables in the dimensions.
Label Optional label for the variable. If label is empty, the variable name is displayed.
Color Color for the lines of this variable
Line width The width of the lines drawn for this variable

There are also configuration options for y-axis. They are displayed when an axis is selected in the Axis and variables tab.

Y-axis
Label Label of the axis
Min Minimum value in the axis
Max Maximum value in the axis
Color The color of the axis lines, ticks and labels
Hide Options to hide parts of the axis: label, axis line, tick marks, tick labels

Bar chart

Bar charts display variables in a certain point in time. The point can be defined for the complete chart in cart's general settings and individually for each variable in variable's settings. If no time is defined, the last simulation step is shown in game or normal simulation mode. Playback mode displays the current playback time.

General bar chart settings

General
Name The name of the chart. Name is displayed in model browser.
Title The title of the chart. Title is displayed in the chart, if not hidden
Type Normal or stacked. Normal type displays all variables in separate bars. Stacked mode is able to display one-dimensional array variables in stacks.
Time A general time setting for the entire chart. The values at this time are displayed in the chart
Hide Option for hiding grid, title and Legend. Legend is hidden by default since it only contains the information that the current simulation is displayed.

Bar chart in a model diagram

Charts can be added to model diagrams simply by dragging the chart from model browser to diagram.

Settings for a variable in bar chart

Variable
Add / Remove Add or remove variables
Variable The full path of the variable. This field has content assistant, which shows the possible variables as you type.
Range If the variable has multiple dimensions, this field allows to select, which dimensions are displayed. It is also possible to sum the variables in the dimensions.
Label Optional label for the variable. If label is empty, the variable name is displayed
Time The time for which the value of this variable is displayed. If this is empty, a general (or default) value of the time is used.

Settings for axis in bar chart

Domain axis
Label An optional label for x axis
Label angle Axis labels can be displayed in an angle. This value determines the angle in degrees.
Color The color of the axis line, tick lines and labels
Hide Options for hiding axis label, line, tick marks and tick labels
Range axis
Label An optional label for y axis
Min Minimum value of y axis
Max Maximum value of y axis
Color The color of the axis line, tick lines and labels
Hide Options for hiding axis label, line, tick marks and tick labels

Pie chart

Pie charts display proportions of variables at a certain time. The time, similar to bar charts, can be configured for the whole chart and individually for each variable. If no time is configured, the value at the latest time step (or the time of simulation playback) is displayed.

General pie chart settings

General
Name The name of the chart. Displayed in e.g. model browser
Title Title The title of the chart. Displayed in the chart
Time A general time setting for the entire chart. The values at this time are displayed in the chart
Hide Options for hiding title, section labels and legend

Variable settings in pie chart

Variables
Add / Remove Add or remove variables
Variable The full path of the variable. This field has content assistant, which shows the possible variables as you type.
Range If the variable has multiple dimensions, this field allows to select, which dimensions are displayed. It is also possible to sum the variables in the dimensions.
Label Optional label for the variable. If label is empty, the variable name is displayed
Color Color for the sector of this variable
Time The time for which the value of this variable is displayed. If this is empty, a general (or default) value of the time is used.
Exploded If set, the sector(s) of the variable(s) are "exploded" away from the center of the pie chart

Pie chart in chart panel

Chart panel

ChartPanel.png

Chart panel is able to house multiple charts at the same time. Chart panel is a normal view in the workbench, which means it can be resized and located anywhere in the workbench. Eclipse also allows detaching the tab completely from the workbench, which means that the chart panel can be dragged to another screen.

Chart panel can be either horizontally or vertically oriented. Different orientations help in configuring the workbench perspective.

Chart panel orientation

Charts are added to the panel by dragging from model browser. The order of the charts is easily changed by dragging charts from their header. Charts can be minimized or completely removed from the buttons in their headers.

By default, chart panel is hidden in the upper-left corner of the system dynamics perspective.

Multidimensional variables

Modeling

Models with multidimensional variables look just like any other models. The structure of the models is replicated givin multiple indexes to variables. For users with programming background, notation Variable[] may be familiar. In system dynamic modeling we need to give names to the indexes. Instead of using numbers to define the indexes, like in normal programming, we use enumerations. The steps to create a model with multidimensional variables are as follows (with examples):

  1. Create the model structure
    Model with multidimensional variables looks just like any other model.
    Model with multidimensional variables
  2. Create Enumerations and define the indexes
    Enumerations are created by right-clicking configuration and selecting New->Enumeration.
    New enumeration
    Define enumeration indexes by adding as many indexes as you want. Rename the indexes by selecting them and then clicking on them again.
    Index definitions
  3. Add Enumerations for variables
    Select the variable that you want to be multidimensional. From the Indexes -tab in the property view, move the wanted enumerations to the right. The order of the enumerations does matter.Select enumerations for variables
  4. Define equations for all possible indexes
    Multidimensional variable can be thought as a multidimensional table. Each cell of the table needs to have an expression or a value. A cell cannot have multiple definitions. All cells can be defined in one expression like in the following example. E1 and E2 have both three indexes, so the resulting definition can be {{1,2,3},{4,5,6},{7,8,9}}.
    Define equations for all indexes


Expressions

Values for all cells in the variable matrix can be defined in a single expression.

DefineEquations.png

DefineEquations2.png

In many situations, it is however more clear to define separate expressions for each cell or blocks.

  1. Define range for the expression
    DefineRange.png
  2. Define the expression (dimensions of the defined range and the result of the expression must match!)
    DefineEquations3.png
  3. Create a new expression
    CreateNewExpression.png
  4. Select the new expression
    SelectNewExpression.png

Repeat until there is value for each cell. There must be exactly one value for each cell.

Array slices

Many times it is useful to access a slice of a multidimensional variable. In this chapter, we will use the following example variable:

enumeration E1 = one, two, three
enumeration E2 = eins, zwei, drei

Auxiliary[E1, E2] = {{1,2,3},{4,5,6},{7,8,9}}

The whole variable can be accessed in three different ways:

  1. Auxiliary
  2. Auxiliary[E1, E2]
  3. Auxiliary[:, :]

The following explains different methods for accessing parts of the variable

// Single cell
Auxiliary[one, eins] = 1
// Slices
Auxiliary[one, E2] = {1,2,3}
Auxiliary[E1, zwei] = {{2},{5},{8}}

// In addition to single cells and the whole enumeration range, a subrange of the enumeration can be used
Auxiliary[two : three, E2] = {{4,5,6},{7,8,9}}
Auxiliary[one : two, zwei : drei]  = {{2,3},{5,6}}

The syntax for accessing parts of variables can be used in both expression range definitions and in expressions.

Arithmetic Operators

Arithmetic operations are defined in Modelica. Below are examples of different operations.

Addition and substraction

Addition and substraction are calculated elementwise

{1,2,3} - {0,2,3} = {1,0,0}
{{2,2},{4,4}} + {{8,8},{1,1}} = {{10,10},{5,5}}
{1,2,3} - {1,2} = ERROR! Different array sizes
Division

Division with a scalar

{5,10,15} / 5 = {1,2,3}
5 / {5,10,15} = ERROR! not allowed

Division with an array

{1,2,3} / {1,2,3} = ERROR! not allowed

// Elementwise
{1,2,3} ./ {1,2,3} = {1,1,1}
Multiplication

Multiplication with a scalar

{1,2,3} * 2 = {2,4,6}
5 * {1,2,3} = {5,10,15}

Matrix multiplication with an array

{{3,4},{5,6}} * {1,2} = {11,17}
{{3,4},{5,6}}  * {1,2,3} = ERROR! incompatible array sizes

//Elementwise
{1,2} .* {1,2} = {1,4}

Real[3,2] c = {{1,2},{3,4},{5,6}};  
Real[2,2] d = {{3,4},{5,6}};
Real[2,2] cd;
cd = c[2:3, :] .* d; // Result: {{9,16},{25,36}}
Power

Elementwise power of array variable

{1,2} .^ 2 = {1,4}

Builtin Modelica functions

Modelica has some builtin functions that help using multidimensional variables. This chapter introduces some of the builtin functions.

Dimension and size functions

We will use the same example variable as previously.

enumeration E1 = one, two, three
enumeration E2 = eins, zwei, drei

Auxiliary[E1, E2] = {{1,2,3},{4,5,6},{7,8,9}}
ndims(A)

The number of dimensions in array A

ndims(Auxiliary) = 2
size(A, i)

The size of dimension i in array A

size(Auxiliary, 2) = 3
size(A)

A vector of length ndims(A) containing the dimension sizes of A

size(Auxiliary) = {3,3}


Construction functions

In addition to normal array constructing, a special construction functions can be used.

zeros(n1,n2,n3,...)

An array full of zeros with dimensions n1 x n2 x n2 x ...

zeros(2, 2) = {{0,0}, {0,0}}
ones(n1,n2,n3,...)

An array full of ones with dimensions n1 x n2 x n2 x ...

ones(2, 2) = {{1,1}, {1,1}}
fill(s,n1,n2,n3)

Like zeros() and ones(), but with user defined value (s) for array elements.

fill(3,2,2) = {{3,3}, {3,3}}
identity(n)

Creates an n x n integer identity matrix with ones on the diagonal and all other elements zero.

identity(3) =
  1 0 0
  0 1 0
  0 0 1
diagonal(v)

Constructs a square matrix with elements of vector v on the diagonal and all other elements zero.

diagonal({1,2,3}) =
  1 0 0
  0 2 0
  0 0 3
linspace(x1,x2,n)

Constracts a Real vector from x1 to x2 with n equally spaced elements.

linspace(2,5,4) = {2,3,4,5}
Reduction functions

Reduction functions reduce arrays to scalars.

min(A)

Returns the minimum value in array A.

Real A = {{1,2},{3,4}}
min(A) = 1
max(A)

Returns the maximum value in array A.

Real A = {{1,2},{3,4}}
max(A) = 4
sum(A)

Returns the sum of values in array A.

Real A = {{1,2},{3,4}}
sum(A) = 10    // 1 + 2 + 3 + 4
product(A)

Returns the product of values in array A.

Real A = {{1,2},{3,4}}
product(A) = 24   // 1 * 2 * 3 * 4


Using functions with iterators

Functions min(A), max(A), sum(A) and product(A) reduce arrays to scalars. When you use multidimensional variables, you will most probably like to reduce less dimensions. This can be achieved using iterators with reduction functions. The result is constructed as an array, if curly brackets {} are used to enclose the expression.

enumeration E1 = one, two, three
enumeration E2 = eins, zwei, drei
Auxiliary[E1, E2] = {{1,2,3},{4,5,6},{7,8,9}}

AuxiliarySum[E1] = {sum( Auxiliary[ i , E2 ] ) for i in E1} // Result: {6, 15, 24}

/* Same as:
{sum(Auxiliary[one, E2]), sum(Auxiliary[two, E2]), sum(Auxiliary[three, E2])}
{sum({1,2,3}), sum({4,5,6}), sum({7,8,9})}
*/

One expression can have multiple iterators. The Modelica specification defines the following format, but it is NOT yet supported in OpenModelica:

{sum(Array[ i, j, E3]) for i in E1, j in E2}  // Dimensions reduced from 3 to 2

To use multiple iterators, use the following format (NOTE the reversed order of iterators!):

{{sum(Array[ i, j, E3]) for j in E2} for i in E1}  // Dimensions reduced from 3 to 2

The range doesn't have to be the whole enumeration. Subranges can also be used.

{sum( Auxiliary[ i , eins : zwei ] ) for i in E1.one : E1.two} // Result: {3, 9}

/* Same as
{sum(Auxiliary[one, eins : zwei]), sum(Auxiliary[two, eins : zwei])}
{sum({1,2}), sum({4,5})}
*/

Simulation results

Multidimensional variables provide multiple results for the same variable. One result for each index. The trend view clutters very quickly when you add dimensions to the variables.

Trend view for a 3x3 variable

The clutter can be reduced by selecting which enumeration indexes are shown in charts. Select an enumeration and tick the indexes that you want to show. The same settings apply to each variable that uses the enumeration. This way you can follow an interesting index throughout the model.

Select indexes to show in charts


Array Variables in Modules

Replaceable enumeration inside a module

You can also use array variables inside modules. Enumeration need to be defined separately for each module type and added to all necessary variables, also inputs and outputs (see Modeling).

When defining a module, modeler may not want to restrict the size of the array variable. In many cases the same module structure could be used for both large and small arrays. For example if the module is a project, a project may have three or even twenty phases. In this case the enumerations inside modules need to be set as replaceable.



In this example, the original enumeration that the modeler used had two indexes. We are using an instance of this module, but we need to use array variables that have three indexes instead of two. We are using two variables and an instance of the module.

Module configuration for replacing enumerations inside the module


In the parent configuration, we have used an enumeration with three indexes in the two variables.

Enumeration in the parent configuration that will replace the enumeration in the module


The replacement can be defined in the properties of the module instance. When the module instance is selected, a table with all the replaceable enumerations is shown. By clickin on the cell next to the enumeration, a drop-down box is shown with all the enumerations in the parent module. If an enumeration is selected, it will replace the enumeration inside the module during simulation. The replacement will not, however, show elsewhere in the model.


Enumeration in the parent configuration that will replace the enumeration in the module


When creating replaceable enumerations in modules, modelers need to be very careful. Direct references to single enumeration elements are not allowed, since the enumeration elements will change, if the enumeration is replaced for simulation!

Functions

Modelica provides a convenient way to use functions in your models. You can create your own functions, import complete function libraries and share function libraries to be used in all of your models.

Creating functions

NewFunction.png

You can create new functions to the Functions -folder in your model or to function library folders. Right-click on the folder and select New->Function.

Functions that can be found from Functions -folder can be used in your variable definitions.

Functions are defined using Modelica language. The variables used in the function are defined in the declaration section. Function needs an output and an arbitrary number of inputs. Modelica specification enables use of multiple outputs, but this feature is not supported. The inputs are given in the same order as they are used in calling the function.

Algorithm section defines the actual functionality of the function. In algorithm sections you must use assignments ":=" instead of just plain "=". All the assignments are calculated in the order they are written.

Below is an example of a simple function. For more examples, see the built-in functions provided with the tool and Modelica specifications.


ExampleFunction.png

Function libraries

NewFunctionLibrary.png

Functions can be collected into libraries. This is a good way of organizing your functions. If function is located in a library, you need to append the library name to the function call (e.g. ExampleFunctionLibrary.ExampleFunction(arg1, arg2)). With libraries, you can also have functions with same names (e.g. library1.func(arg) and library2.func(arg)).

There are three types of function libraries: normal function library, built-in function library and shared function library. Normal function libraries can be created to a model and are available only in that model. Built-in libraries are available in all models and calls to built-in functions should not include the library name. Shared functions are available to all models in your workspace, but you need to enable them to each model individually.

Create a new function library by right-clicking on Functions folder, Shared functions folder or other module library.

If you create a shared function library, the library can be added to other models or removed from current model by selecting the Shared Functions -folder and using the properties view.

SharedFunctions.png

External functions

Modelica allows you to use external functions that are programmed using C language. Below is a simple example of using a function that returns the sum of two arguments.

double exampleCFunction(double x, double y)
{
  double res;
  res = x + y;
  return res;
}


ImportCFunctionLibrary.png

You need to create an Object file (.o) of your code, import it to a function and define the function to use the library and function that you created. If you don’t have a C-compiler, you can use the one provided by the Simantics installation (e.g. simantics_home/plugins/org.simantics.openmodelica.win32_x/MinGW/bin/mingw32_gcc.exe -g -O -c ExampleCFunctionCode.c).

ExampleExternalFunction.png


You can use the external function like any other function in your variables. Using equation ExampleExternalFunction(Stock1, 10) in an auxiliary variable (Auxiliary2) gives the following result.

ExternalFunctionResult.png


In addition to .o -files, you can import .c and .h files. By including these to the function, the source code can be reviewed and later by you or others that use the function.


Small external C function can also be written directly to the function code in Simantics Sysdyn:

function ExampleExternalFunction
input Real x;
input Real y;
output Real z;
external "C" z=exampleCFunction(x,y) annotation(Include="
double exampleCFunction(double x, double y)
{
  double res;
  res = x + y;
  return res;
}");
end ExampleExternalFunction;

Modelica functions

Modelica has built-in functions that can be used anywhere and are not visible in function libraries. This section covers a large number of those functions. For functions related to array variables, see see builtin functions for arrays.

abs(x) Returns the absolute value of x. Expanded into "(if x >= 0 then x else -x)".
acos(x) Inverse cosine.
asin(x) Inverse sine.
atan(x) Inverse tangent.
atan2(x1,x2) four quadrant inverse tangent.
cat(n,A,B,...) General concatenation function that concatenates arrays A,B,... along the nth dimension.
ceil(x) Returns the smallest integer not less than x, the closest integer above x.
cos(x) Cosine.
cosh(x) Hyperbolic cosine.
cross(x,y) Returns the 3-vector cross product of the 3-vectors x and y.
delay(expr, delayTime) Returns the value of expr at the time time-delayTime. The value of expr is returned when time <= time.start + delayTime.
der(x) Time derivative of x. X must be have continuous-time variability.
div(x, y) Returns the algebraic quotient x/y with any fractional part discarted. E.g. div(10,3) = 3.
edge(b) Returns true when the value of the boolean expression b changes. Expanded into (b and not pre(b)).
exp(x) Exponential, base e.
floor(x) Returns the largest integer not greater than x, the closest integer below x.
initial() Returns true at the beginning of the simulation.
integer(x) Returns the largest integer not greater than x as an integer.
log(x) Natural logarithm. (base e, x > 0)
log10(x) Base 10 logarithm. (x > 0)
mod(x, y) Returns the integer modulus of x/y: mod(x,y) = x - floor(x/y) * y.
noEvent(expr) noEvent around an expression causes the expression to NOT generate event.

Important:

If you want a condition to be checked only on time steps, use noEvent.
(e.g. if noEvent(value >= 1) then ... else ...)
pre(x) Returns the preceding value of y from time event that has occured before current time.
rem(x, y) Returns the integer remainder of x/y: rem(x,y) = x - div(x,y) * y.
sample(start, interval) Returns true and triggers time events at times start + i * interval (i=0,1,...).
sign(x) Returns -1 if x is negative, 1 if x is positive. Expanded into "(if x > 0 then 1 else if x < 0 then -1 else 0)".
sin(x) Sine.
sinh(x) Hyperbolic sine.
sqrt(x) Square root of x. The value of x must be greater or equal to 0 or an assertion error occurs.
tan(x) Tangent.
tanh(x) Hyperbolic tangent.
terminal() Returns true at the end of a successful simulation.

Spreadsheets

Spreadsheets allow storing and maintaining values in a familiar format. The current spreadsheet implementation is experimental and does not provide much functionality. Some basic tasks can, however, be performed. You cannot type directly to a cell, you need to first copy some value or text into it. Then it can be modified.

Variable values

Spreadsheets are an easy way to import and manage parameter values. Sheet shells can be referenced to in models with the syntax SheetName(CellReference). Below is an example of getting a value for Auxiliary. The text "Auxiliary" is not necessary in the spreadsheet, but it helps maintaining the sheet.

Sheet reference in an expressionValue of a variable in sheet

You can also refer to sheets with array variables. Reference is simple with one and two dimensional array variables. Below is an example of a sheet with values for two different variables. The names of the indexes are not mandatory, but help reading and maintaining the sheet.

Spreadsheet for multidimensional variables

Reference to a sheet range creates an array variable.

Array variable ArrayTest1[2] with sheet reference

ArrayTest2 has two dimensions. Reference to a larger range creates a two-dimensional variable. By default, the indexes are ordered as in this example. However, if you wish to present the indexes in a different order in the spreadsheet, you can transpose the sheet reference (e.g. transpose(ArrayTestSheet(E2:G3)))

Array variable ArrayTest2[2,3] with sheet reference

History data

Spreadsheets can be used to store and display external history data e.g. from some statistics. The history data should be arranged as columns or rows and include a time variable. Reference variables need to be given the exactly same name as their counterparts in model.

Spreadsheet with history data

History data is used by creating a History dataset to an experiment. The dataset is activated and deactivated in charts by double-clicking the dataset. History data reference is set in the properties of the history dataset.

History data settings

History data properties
Name Name of the history dataset
Sheet Combo box for selecting a sheet where the history data is located
Orientation Columns or rows. How the data is arranged in the sheet.
Range start The upper-left corner cell of the history data range
Range end The lower-right corner cell of the history data range
Time variable The name of the time variable in the history data, if it is different than "time"
Variables in range This text field displays the found variables in the given sheet and range

If the history dataset is activated, history data is displayed whenever it is available for a variable that is displayed in a chart.

History data in two charts


Simantics | System Dynamics Modelling


Preferences-system.svg

Tutorial: Basic System Dynamics Modelling
System dynamics modelling in Simantics is a free modelling tool that is included into the basic installation. This tutorial introduces the basic features of the system dynamics modelling tool.

Preferences-system.svg

Tutorial: Advanced System Dynamics Modelling
This tutorial introduces the more advanced features of the system dynamics modelling tool: Modules and Operating interfaces