Example2.Analysis.PFVariants.AutomationFunctions

Provides functions for automation of simulation and linearization using the models of the upper layer package

Information

This package contains functions that automate the simulation, plotting and linearization of the models under Example2.Analysis.PFVariants.

Extends from Modelica.Icons.FunctionsPackage (Icon for packages containing functions).

Package Content

Name	Description
simulate_and_plot_ramping_pf	Simulates and plots the "Ramping_pf" model
simulate_and_plot_inputs_pf	Simulates and plots the "RampingRandomLoadAndInput_pf" which includes random load and the designed input signal.
linearize_KchangeLinRamping_pf	Function to linearize the model "KchangeLinRamping_pf" at any point in time.

Example2.Analysis.PFVariants.AutomationFunctions.simulate_and_plot_ramping_pf

Simulates and plots the "Ramping_pf" model

Information

Usage:

1. In the Package Browser, right click on the function and select "Call function...". This will open the function's window.



2. In the new function window, click on "Execute", if successful, messages/results are displayed in the command window and plots will appear in Simulation tab.
3. Go back to the function's own window and click on "Close".

Sample Output:

Running the function will result in the following plot.

Extends from Modelica.Icons.Function (Icon for functions).

Example2.Analysis.PFVariants.AutomationFunctions.simulate_and_plot_inputs_pf

Simulates and plots the "RampingRandomLoadAndInput_pf" which includes random load and the designed input signal.

Information

Usage:

1. In the Package Browser, right click on the function and select "Call function...". This will open the function's window.



2. Enter the desired number of cores to use in "corenum". The number of cores to use depends on your computer's resources. , if your computer has more than 12 cores, you can enter the desired number of cores to use in the function's main window. Note that it will take some trial an error to find a trade off on the number of cores to use in your system as per Amdahl's law. See the note below the sample output for more details.
3. In the new function window, click on "Execute", if successful, messages/results are displayed in the command window and plots will appear in Simulation tab.
4. Go back to the function's own window and click on "Close".

Sample output:

Running the function will result in the following plot.

Note:

Please observe that running this function involves simulating the Example2.Analysis.PFVariants.RampingRandomLoadAndInput_pf model, which takes subtantially more time compared to example Example2.Analysis.PFVariants.Ramping_pf, the reason is that both random load and probing signals are added to the model which results in an additional computation burden.

This function calls the Example2.Utilities.SetupSolverSettings.On function that configures several flags to improve simulation speed, however, even by enabling these features, the simulation is time consuming.

For example, running the model using a computer with an 13th Gen Intel(R) Core(TM) i9-13900K, 3.00 GHz, processor and 128 GB DDR4 RAM, setup to use 12 cores, gave the following metrics:

Integration terminated successfully at T = 4050

CPU-time for integration : 1.17e+03 seconds

CPU-time for one grid interval : 234 milliseconds

CPU-time for initialization : 0.005 seconds

Number of result points : 959331

Number of grid points : 5001

Number of accepted steps : 9172936

Number of f-evaluations (dynamics) : 18295816

Number of crossing function evaluations : 9668298

Number of Jacobian-evaluations : 4006190

Number of model time events : 476996

Number of input time events : 0

Number of state events : 2744

Number of step events : 0

Minimum integration stepsize : 2.78e-10

Maximum integration stepsize : 0.00922

Maximum integration order : 5

Looking at the first two lines of the metrics shows that running the model for 67.5 min (4050 sec.) took 19.5 min (1.17e+03 sec.).

Naturally, using a computer with a weaker processor and lower memory will lead to longer run times.

If there is no interest to take into account the random load and probing signal, it is recommended that you run instead the function Example2.Analysis.PFVariants.AutomationFunctions.simulate_and_plot_ramping_pf, which takes substantially less time to execute.

Extends from Modelica.Icons.Function (Icon for functions).

Inputs

Example2.Analysis.PFVariants.AutomationFunctions.linearize_KchangeLinRamping_pf

Function to linearize the model "KchangeLinRamping_pf" at any point in time.

Information

Usage:

1. In the Package Browser, right click on the function and select "Call function...". This will open the function's window.



2. Enter the time for linearization, Example2.Analysis.PFVariants.KchangeLinRamping_pf will be simulated until that point in time and linearized.
3. In the new function window, click on "Execute", if successful, messages/results are displayed in the command window and plots will appear in Simulation tab.
4. Go back to the function's own window and click on "Close".

Sample Output:

Executing the function will produce the following output in the "Commands" window.

Example2.Analysis.PFVariants.AutomationFunctions.linearize_KchangeLinRamping_pf();

Setting up things...

DAE Mode is turned off.

Global optimization is disabled.

Sparse options disabled.

Number of cores reset to default:1

Starting simulation and linearization...

Declaring Modelica_LinearSystems2.StateSpace ss 

System not printed since too large (only dimensions):

   ss.A[52,52]

   ss.B[52,2]

   ss.C[1,52]

   ss.D[1,2]

The linear model is in your working directory.

Open the dslin.mat file following the instructions in the documentation.

You can also analyze the linear model using Dymola, for example typing in the commands window

  Modelica_LinearSystems2.StateSpace.Plot.polesAndZeros(ss, zeros=false, print=true);

  will give you the polte plot shown in the Simulation tab and printed below.

Must translate function Modelica_LinearSystems2.Math.Matrices.LAPACK.dgeevx_eigenValues.

Compilation of function Modelica_LinearSystems2.Math.Matrices.LAPACK.dgeevx_eigenValues successful.

Finished 

Eigenvalues

The following pole plot is generated in the Simulation tab:

Post-processing in MATLAB:

The obtained linear model can be used in any other environment. The linear model is available in the file, dslin.mat, that will appear under your Dymola working directory. It can be loaded in MATLAB using the function:

[A,B,C,D,xName,uName,yName] = tloadlin('dslin.mat')

To use this function, add to the MATLAB workspace the directory and sub-directories under: C:\Program Files\Dymola 2024x\Mfiles

Extends from Modelica.Icons.Function (Icon for functions).

Inputs