2004年第32卷第1期流体机械59文章编号:1005)0329(2004)01)0059)03Fortran与Matlab混合编程实现溴化锂水溶液的物性计算可视化王书中,由世俊,孙贺江,李秋生(天津大学,天津300072)摘要:简要介绍了Fortran与Matlab混合编程的基本原理,通过Fortran与Matlab的接口函数,实现了Matlab对Fortran函数的直接调用,并以溴化锂水溶液物性程序为基础,利用Matlab强大的矩阵计算及图形绘制与编辑功能,实现了溴化锂水溶液物性计算的可视化,为进行溴化锂吸收式制冷的热力计算及分析提供了方便、快捷的图形计算工具。关键词:Fortran;Matlab;溴化锂水溶液;物性计算中图分类号:TB64;TP39文献标识码:AVisualizationofH2O-LiBrSolutionPhysicalPropertiesComputationbyMixedProgrammingwithFortranandMatlabWANGShu-zhong,YOUSh-ijun,SUNHe-jiang,LIQiu-sheng(TianjinUniversity,Tianjin300072,China)Abstract:ThebasicprincipleofmixedprogrammingwithFortranandMatlabisbrieflyintroduced.basedonthephysicalpropertiessubroutinesofH2O-LiBrsolutioncompliedinFortranandusingthepowerfulmatrixcomputationandfiguredrawingandcompilingfuctionofMatlab,thephysicalpropertiescomputationofH2O-LiBrsolutionisvisualized,aconvenientandeffectivetoolforthehea-tdynamicscomputationandanalyseofLiBrabsorptionrefrigerationcycleissupplized.Keywords:Fortran;Matlab;H2O-LiBr;physicalpropertiescomputation1前言运算速度快、Matlab矩阵计算及图形处理功能强的优点,实现了溴化锂水溶液物性计算的可视化。随着溴化锂吸收式制冷机组在技术与应用上本文基于Matlab的科学运算环境,利用Mat-的不断发展,人们对机组的优化设计、变工况条件lab外接程序接口API技术,通过Mex接口函数编下最佳运行参数的动态调整和自动控制等方面提译由Fortran6.5编写的溴化锂水溶液物性计算子出了更高的要求。传统计算要借助于溴化锂水溶程序,生成Matlab环境下可以直接调用的动态链液的物性数据图表,从图表上查出或手算出有关接子程序DLL文件,然后充分利用Matlab强大的物性参数和热力过程特性数据。这种方法繁琐,图形绘制及编辑功能,编程生成1M文件。并在费时,误差较大,已远远不能满足目前的应用需图形用户界面设计的基础上,编写实用且易于操求。为此,能否找到一种可以动态绘图,用鼠标实作的用户界面,实现人机交互绘图、实时显示状态时显示状态点参数,根据已知的的状态参数求其点、自动绘制流程图等一系列功能。他未知的参数,并根据用户需要直接绘制溴化锂吸收式制冷循环流程图的方法,已成为急需解决2实现溴化锂水溶液物性计算可视化的步骤的问题。基于此,本文从满足以上几种功能出发,利用Fortran与Matlab混合编程,充分结合Fortran(1)Fortran环境下溴化锂物性子程序的编写收稿日期:2003)03)0460FLUIDMACHINERYVol132,No11,2004溴化锂水溶液的物性实验数据经过多次改进,并运行结果如下:不断扩充其使用范围,由McNeely给出的溴化锂>>Mexhlibr.f%水溶液的平衡方程是常见的形式[1];文献[2]则在编译前面已经编写完的Fortran-Mex文件总结前人实验数据的基础上,将溴化锂水溶液的hlibr.f使用范围从原来的浓度范围0~65%,温度范围0>>h=hlibr(40,50)%~120e扩充到0~75%,0~190e;文献[3]则提计算浓度N=40%,温度T=50的焓值供了国内溴化锂水溶液的物性数据图表。本文的h=106.6005%溴化锂水溶液物性子程序的编制是基于文献[1]焓值的计算结果的方程形式,并利用文献[2]的方程与数据将其范