基于SALOME与Python的翼型压力载荷计算

Hey, 看起来你在MEDCoupling处理二维翼型压力载荷时卡壳了，你的核心思路是完全正确的：必须先定位翼型的边缘边界单元，再关联对应的压力场数据来计算载荷。下面我帮你梳理正确的步骤，同时修正代码里的关键问题：

1. 正确提取翼型的边缘边界单元 #

在二维网格里，翼型的边缘是1D线单元，这类单元只会被一个2D面单元共享（内部线单元会被两个面单元共享）。你之前尝试用连通性来筛选，但其实MEDCoupling有更直接的API：

import medcoupling as mc
import numpy as np

# 读取2D网格文件
ailMesh2DMc = mc.ReadUMeshFromFile("naca10.med", "MeshFlow10", 0)

# 直接提取所有1D边界单元（只被单个2D单元共享的线单元）
all_boundary_cells = ailMesh2DMc.getBoundaryCells()

# 筛选出属于翼型的边界单元（如果网格包含远场等其他边界，需要通过坐标范围过滤）
# 先获取边界单元的中心坐标
boundary_centers = all_boundary_cells.computeCellCenterOfMass()
# 假设你的翼型x范围在-1到1，y范围在-0.5到0.5，根据实际网格调整这个区间
wing_ids = boundary_centers.findIdsInRange(
    mc.DataArrayDouble([-1.0, -0.5]), 
    mc.DataArrayDouble([1.0, 0.5])
)
# 得到翼型的边缘网格
ailSkinMc = all_boundary_cells[wing_ids]

2. 正确读取压力场并关联到翼型边缘 #

你的代码里错误地覆盖了读取的压力数组，导致数据丢失。压力场通常定义在2D面单元上，所以我们需要找到翼型边缘所属的父面单元，再提取对应压力：

# 读取存储在med文件中的压力场（ON_CELLS表示压力定义在2D面单元上）
pressure_field = mc.ReadFieldCell(
    "naca10.med", 
    ailMesh2DMc.getName(), 
    0, 
    "Pressure", 
    200, 
    201
)

# 获取每个边界单元对应的父2D面单元
parent_face_cells = ailMesh2DMc.getParentCellsOfBoundary(all_boundary_cells)
# 筛选出翼型边缘对应的父面单元
wing_parent_faces = parent_face_cells[wing_ids]

# 提取这些面单元上的压力值
wing_pressures = pressure_field.getArray()[wing_parent_faces]

# （可选）创建翼型边缘上的压力场，方便后续可视化或处理
wing_press_field = mc.MEDCouplingFieldDouble.New(mc.ON_CELLS, mc.ONE_TIME)
wing_press_field.setMesh(ailSkinMc)
wing_press_field.setArray(wing_pressures)
wing_press_field.setName("Pressure_on_Airfoil")

3. 计算二维压力载荷 #

二维场景下，载荷计算是压力 × 单元长度 × 法向量（默认单位厚度，所以面积等价于长度）：

# 转换压力单位（假设读取的是bar，转成Pa）
wing_pressures *= 1e5

# 获取每个翼型边缘单元的长度
edge_lengths = ailSkinMc.getMeasureField(True).getArray()

# 获取每个边缘单元的法向量（注意方向：如果力的方向反了，给法向量乘-1即可）
edge_normals = ailSkinMc.buildOrthogonalField().getArray()

# 计算每个单元的力向量
cell_force_vectors = wing_pressures * edge_lengths * edge_normals

# 计算总载荷（x和y方向分量）
total_force = cell_force_vectors.accumulate()
print(f"总力X分量: {total_force[0]:.2f} N/m，总力Y分量: {total_force[1]:.2f} N/m")

代码里的关键错误修正 #

• 你之前写的arr=mc.DataArrayDouble(ailMesh2DMc.getNumberOfCells(),1)完全覆盖了读取的压力数组，导致所有压力数据丢失，这是最核心的问题。
• 提取边界单元时，用getBoundaryCells()比手动处理连通性更可靠，MEDCoupling已经封装了这个边界提取逻辑。
• 压力场定义在2D面单元上，边缘单元本身没有压力数据，必须从它所属的父面单元获取。

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