我如何使用 'fmi2GetDirectionalDerivative'？

Question

我正在尝试从 FMU 获取模型雅可比矩阵，根据以下文献，我可以使用 fmi2GetDirectionalDerivative 来执行此操作，但我不确定我需要做什么。

我的问题是：

我可以在 Dymola 或 MATLAB 中调用这个函数吗？
屏幕截图或示例会很有帮助。

https://ep.liu.se/ecp/132/091/ecp17132831.pdf

Answer 1

我不熟悉在 MATLAB 中调用 DLL 函数，但这是 Python 中的示例。 FMPy (https://github.com/CATIA-Systems/FMPy) 在 python.

中有运行个 FMU 的包装器

我已经针对我在此处编写的一个简单模型对此进行了测试 ()。在这种情况下，已知数是状态或输入的值参考，未知数是导数或输出的值参考。

当通过 Dymola 作为 Model Exchange FMU 而不是 Co-Simulation FMU 导出时，我已成功提取雅可比行列式。

def get_jacobian(fmu, vr_knowns, vr_unknowns):
    """
    populates jacobian from list of knowns and unknowns
    can be only called after the current sim time and inputs are set
    """
    jacobian = []
    try:
        for vr_known in vr_knowns:
            for vr_unknown in vr_unknowns:
                jacobian.extend(
                    fmu.getDirectionalDerivative(
                        vUnknown_ref=[vr_unknown],
                        vKnown_ref=[vr_known],
                        dvKnown=[1.0]
                    ))
        print_status(f'Jacobian Elements: {jacobian}')
    except Exception as e:
        print("[ERROR] cannot compute jacobian at current timestep")
        print(f"[ERROR] {e}")

我使用此代码片段使用 FMPy 收集状态和导数的值引用：

# get FMU model description object
model_description = fmpy.read_model_description(
    os.path.join(fmu_path, fmu_filename)
)

# collect the value references
vrs = {}
for variable in model_description.modelVariables:
    vrs[variable.name] = variable.valueReference

# collect list of states and derivatives
states = []
derivatives = []
for derivative in model_description.derivatives:
    derivatives.append(derivative.variable.name)
    states.append(re.findall('^der\((.*)\)$',derivative.variable.name)[0])

# collect the value references for states and derivatives
vr_states = [vrs[x] for x in states]
vr_derivatives = [vrs[x] for x in derivatives]

Answer 2

三菱电机在最近的 NAM Modelica Conference 2020 上发表了一篇可能相关的论文。

我如何使用 'fmi2GetDirectionalDerivative'？

How could I use 'fmi2GetDirectionalDerivative'?

modelica

dymola

fmi