1. 问题描述
复杂几何体积的实体单元网格划分，如果网格数值设置不合适，有可能在语法完全正常的情况下报错，无法运行。
如下：
QQ截图20161126111406.jpg-26kB
下面以一个半圆环形柱体的网格划分为例，说明这类问题出现的情况及避免的办法。
image_1b2g46fb61jmmvtm8s91ng01d1jl.png-103.6kB
2. 技术背景
ANSYS有限元网格划分技术。
3. 解决方案
实体单元网格划分的一般流程是：
- 几何建模
-
lesize命令，划分轮廓线； -
vsweep命令，扫掠生成有限元网格；
轮廓线网格密度的合理设置，可以避免生成有限单元时发生错误。
4. 实施示例
4.1 问题的重现
执行如下代码：
FINISH $ /CLEAR
/PREP7
!------------------------------单元材料-----------------------------------------
ET,1,SOLID185 ! 钢材
MP,EX,1,20.5E10
MP,PRXY,1,0.28
MP,DENS,1,7850
!------------------------------钢管几何建模-------------------------------------
CYL4,0,0,0.146,0,0.15,180,3 ! 半圆环
! 钢管网格划分
LESIZE,1,,,4 ! 径向厚度划分-V1
LESIZE,3,,,4
LESIZE,6,,,4
LESIZE,8,,,4
LESIZE,2,,,18 ! 圆弧划分
LESIZE,4,,,18
LESIZE,5,,,18
LESIZE,7,,,18
*DO,ILOOP,9,12 ! 轴线划分
LESIZE, ILOOP,,, 30
*ENDDO
VATT,1,,1 ! 材料单元附加
VSWEEP,1
以上代码意味着，对于一个半圆柱环；
-
轴线划分30段
image_1b2g49uaf1sqr1r1ssp91skr1lm12.png-18.9kB - 内外圆弧划分18段
-
厚度方向划分4段
image_1b2g4b0n85581s3v5i0kd9169a1f.png-11.6kB
出现错误：
image_1b2g4e70b8u05c113a4k92jrp1s.png-11.6kB
4.2 问题分析
通过对比确认，发现这一段代码本身语法没有问题。问题在于：网格线的划分线段长度比例过于悬殊。以上划分中网格尺度为：
- 轴向：3000/30 = 100 mm
- 厚度：4/4 = 1 mm
- 圆弧：300*3.14/2/18 = 26
轴向网格尺度为厚度网格100倍，比例过于悬殊，导致ANSYS无法自动划分网格。要解决这一问题，需要轴向网格尺度，使之比例协调。
4.3 网格参数修改
轴向网格代码修改为：
*DO,ILOOP,9,12 ! 轴线划分
LESIZE, ILOOP,,, 100
*ENDDO
运行结果如下：
三维网格.png-112.7kB
横截面网格.png-9.7kB
5. 常见问题
从经验来看，一般对于较为规则的形体，自动网格划分不会出现问题。对于本例的半圆环柱体，应当格外注意。