PyMAPDL风轮机建模求助：创建3个120度间隔叶片关键点遇阻

问题原因 #

你当前代码的核心问题是工作平面旋转是叠加执行的：每次循环调用mapdl.wprota(thzx=angle)时，是在当前已旋转的工作平面基础上继续旋转，而非基于全局坐标系旋转。比如第二次循环时，实际旋转角度是0°+120°=120°，第三次则是120°+240°=360°（回到初始方向），最终导致三个叶片的关键点方向重复，出现建模异常。

另外，每次循环都调用mapdl.kwpave定位工作平面到机舱关键点，这一步本身没问题，但叠加旋转的逻辑错误直接导致了后续关键点位置错误。

修正方案 #

推荐使用局部坐标系来创建叶片关键点，这种方式无需依赖工作平面的旋转，定位更精准：

blade_radius = 50000
blade_count = 3

# 获取机舱关键点坐标（假设nacelle_keypoints[-1]是机舱的安装中心点）
nacelle_kp = nacelle_keypoints[-1]
kp_coords = mapdl.kpt(nacelle_kp)  # 获取该关键点的全局坐标

# 在机舱中心点创建局部柱坐标系（Z轴为全局Z轴，用于旋转定位叶片）
local_cs = mapdl.local(11, 1, kp_coords[0], kp_coords[1], kp_coords[2])

blade_tips_keypoints = []
blade_lines = []

for i in range(blade_count):
    angle = i * 360 / blade_count
    # 在局部柱坐标系下创建叶片尖端关键点（r=blade_radius, theta=angle, z=0）
    tip_kp = mapdl.k(0, blade_radius, angle, 0, csys=11)
    blade_tips_keypoints.append(tip_kp)
    
    # 创建叶片线条
    blade_lines.append(mapdl.l(p1=nacelle_kp, p2=tip_kp))

# 恢复全局坐标系（可选，不影响后续操作）
mapdl.csys(0)

如果坚持使用工作平面的方式，需要每次循环后重置工作平面到全局状态：

blade_radius = 50000
blade_count = 3

blade_tips_keypoints = []
blade_lines = []

for i in range(blade_count):
    angle = i * 360 / blade_count
    
    # 重置工作平面到全局坐标系
    mapdl.wplane(0)
    # 移动工作平面到机舱关键点
    mapdl.kwpave(p1=nacelle_keypoints[-1])
    # 基于全局坐标系旋转工作平面
    mapdl.wprota(thzx=angle)

    blade_tips_keypoints.append(mapdl.k(0, 0, blade_radius))

for i in range(blade_count):
    blade_lines.append(mapdl.l(p1=nacelle_keypoints[-1], p2=blade_tips_keypoints[i]))

说明 #

• 局部坐标系方法更稳定，尤其适合需要精准角度定位的场景，避免工作平面操作的叠加干扰。
• 两种方案都确保三个叶片的尖端关键点分别位于0°、120°、240°的方位上，与中心点（机舱关键点）连接后形成正确的三叶结构。

内容的提问来源于stack exchange，提问作者Nicolás Olabarría