Pause JModelica et Passe Entrées incrémentielles lors de la simulation

Question

Salut Modelica Communauté,

Je voudrais lancer deux modèles en parallèle JModelica mais je ne suis pas sûr de savoir comment passer des variables entre les modèles. Un modèle est un modèle python et l'autre est un modèle EnergyPlusToFMU.

Les exemples de la documentation JModelica contiennent les entrées de période de simulation complètes définies avant la simulation du modèle. Je ne comprends pas comment on pourrait configurer un modèle qui fait une pause pour les entrées, ce qui est une caractéristique clé des UFA et de la co-simulation.

Quelqu'un peut-il me fournir un exemple ou un morceau de code qui montre comment cela pourrait être implémenté dans JModelica? Est-ce que je mets la commande de simulation dans une boucle? Si oui, comment puis-je gérer les périodes d'échauffement et l'initialisation sans perdre de données au cours des périodes précédentes?

Merci pour votre temps,

Justin

Answer 1

réponse tardive, mais dans le cas où il est repris par d'autres ...

Vous pouvez en effet mettre la simulation en boucle, vous venez besoin de garder une trace de l'état de votre système, de sorte que vous pouvez le ré-initier à chaque itération. Prenons l'exemple suivant:

Ts = 100 
x_k = x_0 

for k in range(100): 
    # Do whatever you need to get your input here 
    u_k = ... 

    FMU.reset() 
    FMU.set(x_k.keys(), x_k.values()) 

    sim_res = FMU.simulate(
     start_time=k*Ts, 
     final_time=(k+1)*Ts, 
     input=u_k 
    ) 

    x_k = get_state(sim_res) 

Maintenant, j'ai écrit une petite fonction pour saisir l'état, x_k, du système:

# Get state names and their values at given index 
def get_state(fmu, results, index): 
    # Identify states as variables with a _start_ value 
    identifier = "_start_" 
    keys = fmu.get_model_variables(filter=identifier + "*").keys() 
    # Now, loop through all states, get their value and put it in x 
    x = {} 
    for name in keys: 
     x[name] = results[name[len(identifier):]][index] 
    # Return state 
    return x 

Cela repose sur la configuration "state_initial_equations": True option de compilation.