Matlab调用Refprop查询物性参数

Matlab调用Refprop的方法教程 共2个教程，refprop查物性参数 教程一 Matlab调用物性计算是仿真工程中常见的问题。在64位环境下的调用稍有不同，在使用中，需要掌握物性函数必要的使用技巧。((所用到的文件，见V1附件，下载解压 就可以直接运行） 一、matlab 调用必备的程序工具 1、matlab 2、refprop（9或9.1）安装包均可 3、refpropm.m 和rp_proto64.m(这个文件是为了转化dll的数据型使其在64位平台上可使用 4、REFPRP64_thunk_pcwin64.dll，这个动态链接库文件是为了把数据从32位转化为64位 二、调用方法 首先要安装refprop软件 安装好refprop后，找到refprop的安装目录，找到其中的refprop64.dll文件 建立一个matlab工作文件夹，如work文件夹 在work文件夹下放入如下文件夹/文件 refprop文件夹 refprop文件夹下必须包括fluids文件夹（refprop安装目录下），refprop64.dll,REFPRP64_thunk_pcwin64.dll refpropm.m rp_proto64.m 自己编写的运行脚本 总的来说，也就是现在需要有一个work文件夹，work中包含了refprop文件夹和refpropm.m rp_proto64.m 文件。refprop文件夹中需要有refprop64.dll,REFPRP64_thunk_pcwin64.dll和fluids文件夹，然后在work文件夹中建立主函数执行脚本，运行物性函数。具体物性函数的使用参照下文。 5.将matlab的工作文件夹调整为work文件夹 （或者不想调，可以直接在work文件夹下建立脚本，运行时候直接会跳出是否将当前文件夹设置为工作文件夹，选是就行） （原理说明部分，可跳过） 三、物性函数使用方法 在matlab调用refprop中，任何物性计算都是通过refpropm子函数实现的，也就是说，refprop将详细的调用工作集成在了refpropm子函数中了。这样大大方便基础的物性运算，但对与更底层的计算例如物质相互作用参数等进行了隐藏了。 必要的函数使用说明在refpropm.m的说明部分都已经包含了，使用说明部分： % Examples: % 1) P = refpropm('P','T',373.15,'Q',0,'water') gives % Vapor pressure of water at 373.15 K in [kPa] % % 2) [S Cp] = refpropm('SC','T',373.15,'Q',1,'water') gives % Entropy and Cp of saturated steam at 373.15 K % % 3) D = refpropm('D','T',323.15,'P',1e2,'water','ammonia',[0.9 0.1]) % Density of a 10% ammonia/water solution at 100 kPa and 323.15 K. % % 4) [x y] = refpropm('X','P',5e2,'Q',0.4,'R134a','R32',[0.8, 0.2]) % Temperature as well as gas and liquid compositions for a mixture % of two refrigerants at a certain pressure and quality. % Note that, when 'X' is requested, two variables must be sent, the % first contains the liquid phase composition and the second % the vapor phase composition. % % 5) T=refpropm('T','C',0,' ',0,'water') % Critical temperature % % 6) T=refpropm('T','M',0,' ',0,'r410a.mix') % Maximum temperature that can be used to call properties. % Shows how to call a predefined mixture. 教程二 Matlab 调用Refprop REFPROP（REference Fluid PROPerties）是一款国际权威工质物性计算软件，该物性计算软件能够被MATLAB和Fortran调用计算。 example: ① P = refpropm('P','T',373.15,'Q',0,'water') 其中第一个元素P为函数的返回值；第二个元素温度T为输入值为373.15K，第三个为干度Q,为0。因此，此式的含义就是温度为373.15K时，水的压力值 P为1.014179966562579e+02kPa ② [S Cp] =refpropm('SC','T',373.15,'Q',1,'water') 此式是为了求取温度为373.15K时的饱和水蒸气的定压比热容和焓值 ③ D =refpropm('D','T',323.15,'P',1e2,'water','ammonia',[0.9 0.1]) 此式是为了求取10%的water和90%的ammonia的密度； ④ [x y] = refpropm('X','P',5e2,'Q',0.4,'R134a','R32',[0.8, 0.2]) 此式是为了求取在特定压力和质量下两种混合制冷剂的气相的质量分数和液相的质量分数。 x = 0.857407902025834 0.142592097974166 y = 0.713888145117291 0.286111854882709 ⑤ T=refpropm('T','C',0,' ',0,'water') 此式是求取水的临界温度 ⑥ D_viscosity=refpropm('V','T',273.15,'P',101.325,'air.mix') 此式是为了求取温度为273.15K，压力为101.325kPa时空气的运动黏度为1.725784151112505e-05。 其中，在查取如R410A时，其有两种选择，一种是假设为纯流体，一种是混合物。两种计算结果如下： 前者： 密度为：12.2849 后者 T_average_ref=refpropm(‘D’,‘T’,245.7935,‘P’,300,‘r410a.ppf’)