如何用Modelica面向对象方法实现空气组分质量分数计算？

原代码核心问题 #

1. 变量关联断裂：Constituent模型内的M_total是独立变量，未与全局计算的空气总摩尔质量关联，导致w = y*M/M_total方程因M_total无定义无法求解。
2. 实例独立性错误：Example_1_1中定义的constituents与Composition内部的constituents是完全独立的数组，Composition计算的M_total无法传递给外部组分实例。
3. 循环范围错误：sum(constituents[i].M*constituents[i].y for i in size(constituents))中，size(constituents)返回的是数组维度（如(3,)），无法遍历所有组分，应使用i in 1:nc。
4. 顶层模型无输出关联：Example_1_1未建立变量传递关系，也未定义输出变量，无法获取计算结果。

修正后的实现代码 #

1. 组分模型（Constituent） #

将M_total改为输入变量，明确从外部接收总摩尔质量，同时将w设为输出变量方便结果导出：

model Constituent
  parameter Real y "摩尔分数";
  parameter Real M "摩尔质量";
  input Real M_total "空气总摩尔质量";
  output Real w "质量分数";
equation
  w = y * M / M_total;
end Constituent;

2. 组分管理模型（Composition） #

封装组分数组与总摩尔质量计算逻辑，将计算得到的M_total传递给每个组分实例：

model Composition
  parameter Integer nc = 3;
  parameter Real y[nc] = {0.781, 0.209, 0.01};
  parameter Real M[nc] = {28, 32, 39.9};
  
  Constituent[nc] constituents(y = y, M = M, M_total = M_total);
  output Real M_total "空气总摩尔质量";
  output Real w[nc] = constituents.w "各组分质量分数";
equation
  M_total = sum(constituents[i].y * constituents[i].M for i in 1:nc);
end Composition;

3. 顶层示例模型（Example_1_1） #

直接实例化管理模型，无需额外方程：

model Example_1_1
  Composition air;
end Example_1_1;

运行验证 #

在OMEdit中运行Example_1_1后，可在仿真结果面板查看：

• air.M_total：约28.96
• air.w：{0.754, 0.231, 0.014}，与预期结果一致

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