### 如何使用Python绘制压焓图 为了实现这一目标，`CoolProp`库是一个非常合适的选择。该库提供了丰富的功能用于处理制冷剂和其他流体的热物理性质数据[^2]。 下面展示一段利用`CoolProp`以及`matplotlib`来创建R134a压焓图（Pressure-Enthalpy Diagram）的示例代码： ```python import CoolProp.CoolProp as CP import numpy as np import matplotlib.pyplot as plt fluid = 'R134a' # 获取饱和线的数据点 p_min, p_max = 1e4, 40e5 # 压力范围 Pa n_points = 100 pressures = np.linspace(p_min, p_max, n_points) h_liq = [] h_vap = [] for p in pressures: h_liq.append(CP.PropsSI('H', 'P', p, 'Q', 0, fluid)) h_vap.append(CP.PropsSI('H', 'P', p, 'Q', 1, fluid)) plt.figure(figsize=(8, 6)) plt.plot(h_liq, pressures / 1e5, label='Saturation Liquid') plt.plot(h_vap, pressures / 1e5, label='Saturation Vapor') # 添加等温线作为辅助线条 T_list = [-10 + i * 5 for i in range(15)] # 温度列表 Celsius colors = ['gray'] * len(T_list) for T, c in zip(T_list, colors): try: p_T = CP.PropsSI(['P'], 'T', 273.15 + T, 'Q', [0], fluid)[0] h_T = CP.PropsSI('H', 'T', 273.15 + T, 'P', p_T, fluid) plt.axvline(x=h_T, color=c, linestyle='--', alpha=0.5) except Exception as e: continue plt.title(f'Pressure-Enthalpy Diagram of {fluid}') plt.xlabel('Specific Enthalpy (J/kg)') plt.ylabel('Pressure (bar)') plt.legend() plt.grid(True) plt.show() ``` 这段程序首先定义了所需的工作介质——在这里选择了常用的制冷剂R134a；接着分别计算不同压力下的液相和气相熵值，并将其绘制成两条曲线表示饱和液体线与饱和蒸汽线之间的区域即为两相区。此外还加入了若干条代表特定温度条件下的垂直虚线以便更好地解释图表中的信息[^1]。