Dymola中无法单独修改数组节点初始焓值的问题求助

问题背景 #

我搭建了一个简单流体模型：包含源、网格化管道和汇。使用Dymola 2024x时，其在初始化变量中选择部分节点的焓h，但初始值不足以找到解，导致求解器初始化仿真失败。

我已知一个更优的猜测值，通过在模型文本中添加h(each start=my_better_value)可完成初始化，但无法仅为Dymola提示的特定节点设置值。例如，在变量浏览器的initial values中仅显示my_module.h[2]，尝试修改该值时，命令面板出现报错：

A heterogenous array in array of records my_module.h is not supported.
Warning: Undeclared variable or command: my_module.h

修改未生效。我确信过去曾独立修改过数组值，想了解该问题原因，以及除each start命令外的初始值设置方法。

编辑 - 实际示例 #

以下是基于ThermoSysPro库的基础测试用例：

model BOP
  ThermoSysPro.Fluid.BoundaryConditions.SourceQ sourceQ(
    Q0=10,
    T0=500);

  ThermoSysPro.Fluid.BoundaryConditions.SinkP sinkP;

  ThermoSysPro.Fluid.HeatExchangers.DynamicOnePhaseFlowPipe
    dynamicOnePhaseFlowPipe(
    Ns=1,
    h(each start=2.e6));
equation 
  connect(dynamicOnePhaseFlowPipe.C1, sourceQ.C);
  connect(dynamicOnePhaseFlowPipe.C2, sinkP.C);
end BOP;

问题原因 #

1. 记录数组类型限制：报错中的"A heterogenous array in array of records"表明my_module.h是记录类型数组，而非普通数值数组。Dymola对记录数组的单个元素修改支持有限，因为记录包含多个属性，直接修改下标元素会触发类型不兼容校验失败。
2. 变量路径解析问题："Undeclared variable or command"说明Dymola无法在当前上下文正确识别my_module.h[2]的完整路径，可能是该变量是模型内部隐式生成的变量，未在顶层组件中公开声明，导致变量浏览器的路径映射不准确。

替代初始值设置方法 #

除了each start批量设置外，可采用以下几种针对性方案：

1. 模型文本中指定单个元素初始值 #

若明确需要修改的节点下标，可直接在组件声明中为单个元素设置start值，语法如下：

ThermoSysPro.Fluid.HeatExchangers.DynamicOnePhaseFlowPipe
  dynamicOnePhaseFlowPipe(
    Ns=1,
    h[2].start=2.2e6); // 仅为第2个节点设置初始焓值

注：多维数组需写出完整下标（如h[1,2].start），具体取决于模型中数组的维度定义。

2. Dymola脚本命令赋值 #

通过Dymola命令行脚本操作，先确认变量完整路径，再赋值：

// 列出所有初始化变量，获取目标变量的准确路径
listVariables();
// 为特定节点设置初始值（需替换为实际变量路径）
setInitial("dynamicOnePhaseFlowPipe.h[2]", 2.2e6);

若h是记录类型，需指定记录内的具体属性（如h[2].value，需参考ThermoSysPro库的变量定义）。

3. 添加初始化方程约束 #

在模型的initial equation块中添加强制初始化约束，优先级最高，能确保数值生效：

model BOP
  ThermoSysPro.Fluid.BoundaryConditions.SourceQ sourceQ(
    Q0=10,
    T0=500);

  ThermoSysPro.Fluid.BoundaryConditions.SinkP sinkP;

  ThermoSysPro.Fluid.HeatExchangers.DynamicOnePhaseFlowPipe
    dynamicOnePhaseFlowPipe(
    Ns=1);

initial equation
  // 强制指定第2个节点的焓初始值
  dynamicOnePhaseFlowPipe.h[2] = 2.2e6;
equation 
  connect(dynamicOnePhaseFlowPipe.C1, sourceQ.C);
  connect(dynamicOnePhaseFlowPipe.C2, sinkP.C);
end BOP;

注：需确保该约束与其他初始化条件（如源边界的温度设置）无冲突，否则会导致初始化方程组无解。

4. 初始化窗口高级设置 #

在Dymola的「Simulation Setup」→「Initialization」选项卡中，切换到「Advanced」模式，找到目标变量的完整路径，手动输入初始值。这种方式可绕过变量浏览器的数组类型限制，直接对单个元素赋值。

内容的提问来源于stack exchange，提问作者Songio