在OpenModelica中通过列表选择流体包调用物性函数

你可以通过两种实用方式实现枚举与物性包的关联，自动调用对应气体的物性函数，无需手动修改路径：

方案一：接口包继承+可替换包绑定（推荐，可扩展性强） #

核心思路是先定义统一的物性函数接口，让各气体包继承该接口，再通过枚举绑定对应的实现包。

步骤1：定义通用物性接口包 #

创建包含所有物性函数模板的基础包，确保所有气体包的函数接口一致：

package Project.Fluids.BaseFluid
  function rho_PT
    input Real P "压力";
    input Real T "温度";
    output Real "密度";
  end rho_PT;
  
  function cp_T
    input Real T "温度";
    output Real "定压比热";
  end cp_T;
  // 可添加其他物性函数模板
end BaseFluid;

步骤2：实现各气体的物性包 #

让空气、氮气包继承基础接口，实现各自的物性计算逻辑：

// 空气物性包
package Project.Fluids.Air
  extends Project.Fluids.BaseFluid;
  
  function rho_PT
    input Real P;
    input Real T;
    output Real;
  algorithm
    // 空气密度的具体实现代码
    rho = ...; // 替换成你的计算逻辑
  end rho_PT;
  
  function cp_T
    input Real T;
    output Real;
  algorithm
    // 空气定压比热的具体实现代码
    cp = ...; // 替换成你的计算逻辑
  end cp_T;
end Air;

氮气包结构与空气包完全一致，仅需替换内部函数的实现逻辑。

步骤3：在模型中关联枚举与物性包 #

用replaceable package配合枚举选项，将选择的气体与对应包绑定：

model FluidPropertyCalculator
  type GasType = enumeration(Air, Nitrogen);
  parameter GasType Gas = GasType.Air;
  
  // 根据枚举动态绑定物性包
  replaceable package SelectedFluid = Project.Fluids.Air constrainedby Project.Fluids.BaseFluid
    annotation(choices(
      choice(GasType.Air) = Project.Fluids.Air,
      choice(GasType.Nitrogen) = Project.Fluids.Nitrogen
    ));
  
  input Real P = 101325 "单位：Pa";
  input Real T = 293.15 "单位：K";
  output Real rho "密度";
  output Real cp "定压比热";
equation
  rho = SelectedFluid.rho_PT(P, T);
  cp = SelectedFluid.cp_T(T);
end FluidPropertyCalculator;

切换Gas枚举选项时，SelectedFluid会自动关联到对应物性包，方程无需任何修改。

方案二：if-else条件判断（适合简单场景） #

如果不想维护接口包，可直接在方程中根据枚举值判断调用对应函数：

model SimpleFluidCalculator
  type GasType = enumeration(Air, Nitrogen);
  parameter GasType Gas = GasType.Air;
  
  input Real P = 101325;
  input Real T = 293.15;
  output Real rho;
equation
  rho = if Gas == GasType.Air then Project.Fluids.Air.rho_PT(P, T)
        else Project.Fluids.Nitrogen.rho_PT(P, T);
  // 其他物性函数同理添加判断
end SimpleFluidCalculator;

这种方式直接但扩展性差，新增气体种类时需手动修改所有物性函数的判断逻辑。

注意事项 #

• 方案一中的所有物性包必须严格遵循基础接口的函数参数、返回值定义，否则constrainedby会报错
• 使用Dymola等工具时，切换枚举选项后需重新编译模型，确保加载对应包的函数实现

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