Fortran与Matlab混合编程实现溴化锂水溶液的物性计算可.pdf

2004年第32卷第1期流体机械59

文章编号:1005)0329(2004)01)0059)03

Fortran与Matlab混合编程实现溴化锂
水溶液的物性计算可视化

王书中,由世俊,孙贺江,李秋生
(天津大学,天津300072)

摘要:简要介绍了Fortran与Matlab混合编程的基本原理,通过Fortran与Matlab的接口函数,实现了Matlab对Fortran
函数的直接调用,并以溴化锂水溶液物性程序为基础,利用Matlab强大的矩阵计算及图形绘制与编辑功能,实现了溴化
锂水溶液物性计算的可视化,为进行溴化锂吸收式制冷的热力计算及分析提供了方便、快捷的图形计算工具。
关键词:Fortran;Matlab;溴化锂水溶液;物性计算
中图分类号:TB64;TP39文献标识码:A

VisualizationofH2O-LiBrSolutionPhysicalPropertiesComputation
byMixedProgrammingwithFortranandMatlab

WANGShu-zhong,YOUSh-ijun,SUNHe-jiang,LIQiu-sheng
(TianjinUniversity,Tianjin300072,China)

Abstract:ThebasicprincipleofmixedprogrammingwithFortranandMatlabisbrieflyintroduced.basedonthephysicalproperties

subroutinesofH2O-LiBrsolutioncompliedinFortranandusingthepowerfulmatrixcomputationandfiguredrawingandcompilingfuction

ofMatlab,thephysicalpropertiescomputationofH2O-LiBrsolutionisvisualized,aconvenientandeffectivetoolforthehea-tdynamics
computationandanalyseofLiBrabsorptionrefrigerationcycleissupplized.

Keywords:Fortran;Matlab;H2O-LiBr;physicalpropertiescomputation

1前言运算速度快、Matlab矩阵计算及图形处理功能强
的优点,实现了溴化锂水溶液物性计算的可视化。
随着溴化锂吸收式制冷机组在技术与应用上本文基于Matlab的科学运算环境,利用Mat-
的不断发展,人们对机组的优化设计、变工况条件lab外接程序接口API技术,通过Mex接口函数编
下最佳运行参数的动态调整和自动控制等方面提译由Fortran6.5编写的溴化锂水溶液物性计算子
出了更高的要求。传统计算要借助于溴化锂水溶程序,生成Matlab环境下可以直接调用的动态链
液的物性数据图表,从图表上查出或手算出有关接子程序DLL文件,然后充分利用Matlab强大的
物性参数和热力过程特性数据。这种方法繁琐,图形绘制及编辑功能,编程生成1M文件。并在
费时,误差较大,已远远不能满足目前的应用需图形用户界面设计的基础上,编写实用且易于操
求。为此,能否找到一种可以动态绘图,用鼠标实作的用户界面,实现人机交互绘图、实时显示状态
时显示状态点参数,根据已知的的状态参数求其点、自动绘制流程图等一系列功能。
他未知的参数,并根据用户需要直接绘制溴化锂
吸收式制冷循环流程图的方法,已成为急需解决2实现溴化锂水溶液物性计算可视化的步骤
的问题。基于此,本文从满足以上几种功能出发,
利用Fortran与Matlab混合编程,充分结合Fortran(1)Fortran环境下溴化锂物性子程序的编写

收稿日期:2003)03)04
60FLUIDMACHINERYVol132,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]的方程与数据将其范