PyMAPDL创建圆柱壳模型出现壁面不规则问题求助

问题背景 #

我是Python与MAPDL新手，正在用PyMAPDL编写代码创建圆柱壳模型以研究其力学性能，但生成的实体模型壁面存在不规则现象，该问题已影响网格划分。使用的代码如下：

import numpy as np
from ansys.mapdl.core import launch_mapdl
mapdl = launch_mapdl()
mapdl.clear()
mapdl.prep7()


# 创建实心圆柱
inn_rad=50.0e-2
out_rad=55.0e-2
external_pressure=175.0e6
start_surface=0.0
end_surface=30.0e-2
theta1=0.0
theta2=360.0
Sigma_y=350.0e6
FactorOfSafetey=3
# 设定单位
mapdl.units('SI')
secnum_cyl = mapdl.cyl4(0.0, 0.0, inn_rad, theta1, out_rad, theta2, (end_surface-start_surface))
mapdl.vplot()

尝试过mapdl.cylind(inn_rad, out_rad, start_surface, end_surface, theta1, theta2)命令，结果依旧存在不规则现象。

解决方法 #

1. 优先排查显示精度问题 #

这种壁面不规则大概率是显示渲染问题，而非模型几何本身存在缺陷。可以通过以下方式提升显示精度：

• 调用vplot时指定更高的质量参数：

  mapdl.vplot(quality=100, show_edges=True, edge_color='black')
  

• 在MAPDL中修改全局显示分辨率：

  mapdl.run("/REPLOT,RESOLU,HIGH")
  mapdl.run("/VIEW,1,1,1")  # 切换到等轴测视图，减少透视畸变
  

2. 优化几何创建方式 #

避免依赖单一命令生成几何体，改用分步创建的方式减少参数歧义：

# 定义关键点
mapdl.k(1, inn_rad, 0, start_surface)
mapdl.k(2, out_rad, 0, start_surface)
mapdl.k(3, out_rad, 0, end_surface)
mapdl.k(4, inn_rad, 0, end_surface)

# 创建圆弧与直线
mapdl.larc(1, 2, 0, inn_rad)  # 内圆弧
mapdl.larc(4, 3, 0, out_rad)  # 外圆弧
mapdl.lstr(1, 4)  # 左侧直线
mapdl.lstr(2, 3)  # 右侧直线

# 创建矩形面
mapdl.al(1, 3, 2, 4)

# 旋转面生成圆柱壳体
mapdl.vext(1, theta=360)

3. 几何修复与检查 #

如果模型确实存在几何拓扑问题，执行以下命令修复：

mapdl.run("VGLUE,ALL")  # 粘合所有体，消除潜在的拓扑间隙
mapdl.run("VCHECK,ALL") # 检查体的几何有效性，输出问题报告

4. 网格划分前的预处理 #

网格划分时的参数设置也可能放大显示问题，建议提前指定单元尺寸：

mapdl.esize(0.01)  # 设置单元尺寸为0.01m，可根据模型大小调整
mapdl.vmesh('ALL')  # 对所有体进行网格划分

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