Source code for pyfmi.examples.fmu_with_input_function

#!/usr/bin/env python 
# -*- coding: utf-8 -*-

# Copyright (C) 2014 Modelon AB
# This program is free software: you can redistribute it and/or modify
# it under the terms of the GNU Lesser General Public License as published by
# the Free Software Foundation, version 3 of the License.
# This program is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# GNU Lesser General Public License for more details.
# You should have received a copy of the GNU Lesser General Public License
# along with this program. If not, see <>.

import os as O
import numpy as N
import pylab as p

from pyfmi import load_fmu

curr_dir = O.path.dirname(O.path.abspath(__file__));
path_to_fmus = O.path.join(curr_dir, 'files', 'FMUs')
path_to_fmus_me1 = O.path.join(path_to_fmus,"ME1.0")
path_to_fmus_cs1 = O.path.join(path_to_fmus,"CS1.0")

[docs]def run_demo(with_plots=True): """ Demonstrates how to simulate an FMU with an input function. See also """ fmu_name = O.path.join(path_to_fmus_me1,'SecondOrder.fmu') # Create input object input_object = ('u', N.cos) # Load the dynamic library and XML data model = load_fmu(fmu_name) # Simulate res = model.simulate(final_time=30, input=input_object, options={'ncp':3000}) x1_sim = res['x1'] x2_sim = res['x2'] u_sim = res['u'] time_sim = res['time'] assert N.abs('x1') - (-1.646485144)) < 1e-3 assert N.abs('x2')*1.e1 - (-7.30591626709)) < 1e-3 assert N.abs('u')*1.e1 - (1.54251449888)) < 1e-3 if with_plots: fig = p.figure() p.subplot(2,1,1) p.plot(time_sim, x1_sim, time_sim, x2_sim) p.subplot(2,1,2) p.plot(time_sim, u_sim)
if __name__=="__main__": run_demo()