Modelica理想气体混合包无报错但实例化失败问题求助

一、代码错误排查与修复 #

内部错误多由语法不严谨、类型不匹配或未声明变量引发，以下是你的代码问题及修复方案：

核心问题点 #

• 未声明变量：getMolarFractions函数中totalMoles未在protected块声明，导致实例化失败。
• 类型不匹配：getTotalMoles函数protected块中，将数组类型的properties.MM赋值给标量类型的MM，需改为数组声明。
• 向量运算错误：getHeatCapacity中Mi*heatCapacities是矩阵乘法，实际需要点积计算总热容，应使用Modelica.Math.Vectors.dot函数。
• 未定义符号：getHeatCapacity函数的replaceable input gasProperties[n] properties中n未定义，需改为与输入Mi尺寸匹配的[:]。
• 初始值长度不兼容：mixedGas中Mi的初始值{1}仅1个元素，与nS动态长度不匹配，需改为动态填充。

修复后的完整代码 #

package gasLibrary
 
partial package mixedGas
    replaceable String[:] gasNames = {"nonUsablePackage"};
    parameter Integer nS = size(gasNames,1);
    replaceable gasProperties[nS] properties = {gasProperties()};
    replaceable Modelica.Units.SI.MassFraction[nS] Mi = fill(1/nS, nS);
    replaceable gasState state = gasState(P = 101.3*10^3, T = 300);

      function getHeatCapacity
        input Modelica.Units.SI.MassFraction[:] Mi;
        input gasProperties[:] properties;
        output Modelica.Units.SI.HeatCapacity Cv;
      protected
        Modelica.Units.SI.HeatCapacity[:] heatCapacities = {properties.Cv};
      algorithm
        Cv := Modelica.Math.Vectors.dot(Mi, heatCapacities);
      end getHeatCapacity;

      function getMolarFractions
        input Modelica.Units.SI.MassFraction[:] Mi;
        input gasProperties[:] properties;
        output Real[:] Xi "molar fractions";
      protected
        Modelica.Units.SI.MolarMass[:] MM = {properties.MM};
        Real totalMoles;
      algorithm
        totalMoles := sum(Mi / MM);
        Xi := Mi/MM/totalMoles;
      end getMolarFractions;

      function getTotalMoles
        input gasProperties[:] properties;
        input Modelica.Units.SI.MassFraction[:] Mi;
        output Real totalMoles;
      protected
        Modelica.Units.SI.MolarMass[:] MM = {properties.MM};
      algorithm
        totalMoles := sum(Mi/MM);
      end getTotalMoles;
    
      function getSpecificVolume
            input gasProperties[:] properties;
            input Modelica.Units.SI.MassFraction[:] Mi;
            input gasState state;
            output Modelica.Units.SI.SpecificVolume v;
          algorithm
            v := getTotalMoles(properties, Mi)* Modelica.Constants.R * state.T / state.P;
          end getSpecificVolume;
          function getDensity
            input gasProperties[:] properties;
            input Modelica.Units.SI.MassFraction[:] Mi;
            input gasState state;
            output Modelica.Units.SI.Density rho;
          algorithm
            rho := 1/getSpecificVolume(properties, Mi, state);
          end getDensity;
end mixedGas;

  record gasProperties
    constant Modelica.Units.SI.MolarMass MM = 0;
    constant Modelica.Units.SI.HeatCapacity Cv = 0;
  end gasProperties;

  record gasState
  Modelica.Units.SI.Pressure P;
  Modelica.Units.SI.Temperature T;
  end gasState;

  package singleGasProperties
  gasProperties O2(MM = 0.032, Cv = 950);
  gasProperties N2(MM = 0.028, Cv = 1150);
  gasProperties CO2(MM = 0.044, Cv = 800);
  gasProperties fuel(MM=0.114, Cv = 2200);
  end singleGasProperties;

  package fuelCombustion
  extends mixedGas;
  
  redeclare String[:] gasNames = {"O2", "N2", "CO2", "fuel"};
  redeclare gasProperties[nS] properties = {singleGasProperties.O2, singleGasProperties.N2, singleGasProperties.CO2, singleGasProperties.fuel};
  redeclare Modelica.Units.SI.MassFraction[nS] Mi = {0.2, 0.8, 0, 0};
  
  end fuelCombustion;
end gasLibrary;

二、Modelica深入学习资源建议 #

• 《Modelica语言规范》：官方权威文档，详细涵盖语言语法、语义、面向对象特性及建模规则，是理解Modelica底层逻辑的核心资料。
• 《Modelica工程建模与仿真》：进阶实战类书籍，包含多领域复杂建模案例、仿真优化技巧，适合从入门过渡到专业建模。
• Modelica协会官方进阶教程：由协会推出的系列课程，覆盖库开发、多物理场耦合、性能调优等高级主题，系统性强。
• 开源库源码分析：研读Modelica Standard Library（MSL）、ThermoSysPro等成熟开源库的源码，学习专业建模范式和最佳实践。
• Modelica会议论文集：关注Modelica Conference的会议论文，追踪建模技术前沿动态和行业应用案例。

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