SIMCOS

SIMCOS (an acronym standing for SIMulation of COntinuous Systems) is a computer language and a development environment for computer simulation. In 1989 it was developed by Slovenian experts led by Borut Zupančič.

The purpose of the language is simulation of dynamic mathematical models of systems, given as set of ordinary differential equations. It is an equation oriented and compiler type of language. Despite its name it can be used for discrete simulation as well. The language suits well to the CSSL'67 standard of simulation languages so portability among other languages conforming to the same standard (e.g. Tutsim, ACSL etc.) is quite simple. It is a DOS based software occasionally it is slightly modified so it can be run under actual versions of Microsoft Windows. Apart from the simulation itself it can also perform parametrisation (a series of simulations with different values of parameters), linearisation of models and optimisation (finding such values of parameters that a criterion function is minimised).

When a simulation scheme must be prepared it must be described in the SIMCOS language. It can be "drawn" (similarly as with an analogue computer) using an enclosed block library graphics tool (it contains basic elements such as integrators, amplifiers, summators, some basic input signals etc.) but more often it is entered as a program using one of text editors, e.g. Edit enclosed with DOS. Whichever form of entry of the model is used, the first phase of simulation reprocesses it into space of states form and rewrites the program into Fortran and prepares files with input parameters. This Fortran program is compiled into an executable file (.EXE) and executed. The executable program reads parameter values from input files, performs the simulation and writes requested calculated values into another file. When it terminates, SIMCOS takes control again and can display results as a graphic plot.

The "heart" of the executable is function INTEG which can solve differential equations using one of several numerical methods. First it reads necessary values (e.g. values of parameters, initial conditions) from files then it calls the function DERIV where the model is actually described as series of functions of its derivatives. The returned values are used at the selected numerical method. Requested calculated results are written into the file and the whole procedure is repeated until the termination condition is satisfied.

Continuous simulation of dead time (its Laplace transform is e − s T {\displaystyle e^{-sT}} ) is not a trivial task and usually we use one of Padé approximations. We will simulate Padé approximation of 2nd order

e − s T = ˙ ( s T ) 2 − 6 s T + 12 ( s T ) 2 + 6 s T + 12 {\displaystyle e^{-sT}{\dot {=}}{\frac {(sT)^{2}-6sT+12}{(sT)^{2}+6sT+12}}}

and 4th order:

e − s T = ˙ ( s T ) 4 − 20 ( s T ) 3 + 180 ( s T ) 2 − 840 s T + 1680 ( s T ) 4 + 20 ( s T ) 3 + 180 ( s T ) 2 + 840 s T + 1680 . {\displaystyle e^{-sT}{\dot {=}}{\frac {(sT)^{4}-20(sT)^{3}+180(sT)^{2}-840sT+1680}{(sT)^{4}+20(sT)^{3}+180(sT)^{2}+840sT+1680}}.}

Input signal is a unit step, communication interval equals 0.01s, length simulation run is 5s, results will be compared with output of built-in discrete function delay (it requires additional array (del in our case) of appropriate size).

y1 is a result of simulation of Padé approximation of 2nd order, y2 is a result of simulation of Padé approximation of 4th order and y3 is result of the discrete function delay.

When transfer functions of both Padé approximation are developed using one of simulation schemes, the model can be described with the following program:

program pade
constant tm=1.0
constant tfin=5
array del(101)
variable t=0.0
u=step(t,0.)

u11d=12/(tm*tm)*u-12/(tm*tm)*y1
u11=integ(u11d,0.)
u21d=u11-u*6/tm-y1*6/tm
u21=integ(u21d,0.)
y1=u21+u
u12d=u*1680/(tm*tm*tm*tm)-y2*1680/(tm*tm*tm*tm)
u12=integ(u12d,0.)
u22d=u12-u*840/(tm*tm*tm)-y2*840/(tm*tm*tm)
u22=integ(u22d,0.)
u32d=u22+u*180/(tm*tm)-y2*180/(tm*tm)
u32=integ(u32d,0.)
u42d=u32-u*20/tm-y2*20*tm
u42=integ(u42d,0.)
y2=u42+u
y3=delay(u,tm,#del,ci)

cinterval ci=0.01
hdr Pade approximation of dead time
prepar y1,y2,y3
output 10,y1,y2,y3
termt(t.ge.tfin)

end

After the simulation run is finished the results can be displayed as plots. It is possible to trace values of plots, select which plots to display, turning on a grid, zoom etc.

• Borut Zupančič's homepage
• LMSC download page (the link to SIMCOS is at the bottom)