如何简便读取翼型坐标文件并在Python中绘图及几何计算？

方法一：用Pandas快速读取与处理 #

Pandas能自动跳过注释行、转换数据类型，直接生成可操作的DataFrame，是最省心的通用方案。

1. 依赖安装（未安装时执行） #

pip install pandas matplotlib

2. 读取翼型文件 #

Selig格式的翼型文件开头为注释（翼型名称、说明），后续是x/y坐标列。用read_csv可一键处理：

import pandas as pd

# 读取文件，跳过#开头的注释行，指定列名
airfoil_df = pd.read_csv(
    "b29tip-il.dat",
    comment="#",
    header=None,
    names=["x", "y"],
    delim_whitespace=True
)

# 查看数据样例
print(airfoil_df.head())

得到的airfoil_df中x、y均为数值类型，无需手动转换字符串。

3. 绘图与几何计算 #

翼型绘图（Matplotlib）

import matplotlib.pyplot as plt

plt.figure(figsize=(8, 4))
plt.plot(airfoil_df["x"], airfoil_df["y"], label="B29 Tip Airfoil")
plt.axis("equal")  # 保证坐标比例正确，避免翼型变形
plt.xlabel("x")
plt.ylabel("y")
plt.title("Airfoil Shape")
plt.legend()
plt.show()

几何计算示例

# 计算弦长
chord_length = airfoil_df["x"].max() - airfoil_df["x"].min()

# 分离上下表面，计算厚度分布
upper_surface = airfoil_df[airfoil_df["x"].isin(airfoil_df["x"].unique()[::-1])]
lower_surface = airfoil_df.drop_duplicates(subset="x", keep="last")
thickness = upper_surface["y"].values - lower_surface["y"].values
max_thickness = thickness.max()

print(f"弦长: {chord_length:.4f}")
print(f"最大厚度: {max_thickness:.4f} (占弦长的 {max_thickness/chord_length*100:.2f}%)")

方法二：使用专业翼型处理库 #

如果需要插值、网格生成、气动参数计算等专业操作，推荐用airfoils库：

1. 安装库 #

pip install airfoils

2. 读取与操作翼型 #

from airfoils import Airfoil

# 读取Selig格式文件
airfoil = Airfoil.from_selig("b29tip-il.dat")

# 获取数值类型的坐标数组
x, y = airfoil.coords

# 一键绘图
airfoil.plot()

# 计算最大厚度
max_thickness = max(airfoil.thickness())
print(f"最大厚度: {max_thickness:.4f}")

该库封装了大量翼型专属功能，省去手动编写逻辑的麻烦。

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