分享：网络划分的一点体会-映射网格划分技巧

carpenter

近来看到不少人问有关网格划分的问题，感觉存在一些概念模糊。写一点个人对划分网格的一点体会，希望能对新手有些帮助，抛砖引玉：

几个概念：
1.网格划分主要分为映射网格划分(mapped)和自由网格划分(free)，这里要说明一点的是，存在一个不正确的概念，即对于二维来说映射就是划成四边形。

2.有人过度夸大四边形单元和六面体的精确性，从而贬低自由网格划分。其实，过度“傻瓜”操作的自由网格划分可能容易产生尖锐或边长比很大的不良单元，经控制的自由网格划分的精度不至于降低到哪里去。要说到缺点，采用自由网格划分所产生的三角形单元或四面体单元将使自由度大大增加，从而造成运算量的增加，这可能算是一个缺点，尤其对配置较低的设备。



自由网格划分较简单，以下介绍几个关于映射网格划分的技巧和方法：

1.对于面来说，使用映射网格划分需要满足以下几个几何拓扑结构的要求条件：
          *面必须有三条或四边线的边界组成
          *一般来说，对边的应划分成相等的个数

    这里可以看出，要使用映射网格划分技术，首先要保证面有三条或四条边界线，要达到这个目的，经常使用的方法有：切割、搭接(参考ASBW,VSBW,LCCAT等，这里不详述)，保证分割或搭接后的面满足几何拓扑要求。

以下举一简例说明：

/PREP7
ET,1,PLANE42

BLC4,,,4,4
CYL4,,,1
ASBA,1,2

几何图形如下所示，很显然：这是一个五边形，不能直接采用映射网格划分技术

[ 本帖最后由 yangjungang2004 于 2007-6-26 22:18 编辑 ]

citnamor

fengting1989 发表于 2013-8-12 10:04
顶一下！

学习了很多,感谢各位分享!

carpenter

接下来就可采取不用的方法，使上述图形满足要求。
1.比如沿对称轴分割，变为两个四边形，再划分。这个方法有兴趣的可试一下。
2.这里介绍搭接方法，并将网格划分 成从缺口到边沿由密到疏的渐变过渡的形式。

以下是命令流和划分 结果：

/PREP7
ET,1,PLANE42

BLC4,,,4,4
CYL4,,,1
ASBA,1,2

LESIZE,5,,,16
LSEL,S,,,9,10
LESIZE,ALL,,,16,0.5
ALLSEL,ALL
LSEL,S,,,2,3
LCCAT,ALL
LESIZE,ALL,,,8
ALLSEL,ALL

AMESH,ALL

carpenter

网格：

carpenter

今天先写到这里吧，以后再写关于3D 体划分的方法，大同小异。

如果有人感兴趣，可考虑一下以下的体划分方法

carpenter

体2

渊源1983

那你能不能帮我看看这个怎么样划分呢？
http://www.simwe.com/forum/thread-728384-1-1.html

carpenter

今晚有时间，顺便看了一下有人问过的问题。仅做了其中一例子，网格分的还可以再均匀些。如果有人感兴趣，以后再继续其它的部分

[ 本帖最后由 carpenter 于 2006-7-6 22:58 编辑 ]

半个苹果

响应，楼主号召！1呵呵  不会编辑log文件呢！！正在练习中！！

渊源1983

关键是如何分的呢？？
刚做例子
旁人也只有在一边滴口水的份儿了 ～

akpp01

能不能给出命令流，让我们这些菜鸟学习！
akpp01@163.com
谢谢！

jz200428

稍微改了一下!偷懒,哈哈~~~

/PREP7
ET,1,PLANE42
ET,2,SOLID45

BLC4,,,4,4
CYL4,,,1
ASBA,1,2

LESIZE,5,,,16
LSEL,S,,,9,10
LESIZE,ALL,,,16,0.5
ALLSEL,ALL
LSEL,S,,,2,3
LCCAT,ALL
LESIZE,ALL,,,8
ALLSEL,ALL

AMESH,ALL
EXTOPT,ESIZE,10
EXTOPT,ACLEAR,1
VEXT,3,,,,,5

cj_hmily

D:\ansys work\1.JPGl=10
m=10
h1=5
h2=8
h3=8
r1=1
r2=1.5
x1=-l/2+r1
y1=-m/2
x2=0
y2=-r2
*afun,deg
theta=atan((y2-y1)/(x2-x1))

/prep7

blc5,0,0,l,m,h1
cylind,0,r2,0,h2,0,360  
cyl4,-l/2,-m/2,0,0,r1,90,h3

vovlap,all

wprota,,,90
vsbw,all,,dele  
wprota,,90  
vsbw,all,,dele

vsel,r,loc,x,-l/2,0
vsel,r,loc,y,-m/2,0
wpave,x1,y1,0  
wprota,,-theta
vsbw,all,,dele

wpave,-l/2,-m/2+r1
wprota,,2*theta+270
vsbw,all,,dele  

vsel,inve

wpcsys,1,0
wprota,45,90
vsbw,all,,dele  
wprota,,,90
vsbw,all,,dele  

et,1,45

mshkey,2
mshape,0,3d
esize,0.5
allsel,all
vmesh,all

cj_hmily

D:\ansys work\2.JPGl=10
m=10
h1=5
h2=8
h3=8
r1=1
r2=1.5
*afun,deg
theta=atan(l/m)

/prep7

blc5,0,0,l,m,h1
cylind,0,r2,0,h2,0,360  
cyl4,-l/2,-m/2,0,0,r1,90,h3

vovlap,all

wprota,,,90
vsbw,all,,dele  
wprota,,90  
vsbw,all,,dele

wpave,l/2,0,0  
wprota,,-theta-90  
vsbw,all,,dele
wprota,,-2*(90-theta)   
vsbw,all,,dele

wpave,-l/2,0,0
wpave,-l/2,0,0  
vsbw,all,,dele  
wprota,,2*(90-theta)
vsbw,all,,dele  

et,1,45

mshkey,2
mshape,0,3d
esize,0.5
allsel,all
vmesh,all

cj_hmily

l=10
m=10
h1=5
h2=8
h3=8
r1=1
r2=1.5
*afun,deg
theta=atan(l/m)

/prep7

blc5,0,0,l,m,h1
cylind,0,r2,0,h2,0,360  
cyl4,-l/2,-m/2,0,0,r1,90,h3

vovlap,all

wprota,,90  
vsbw,all,,dele  
wprota,,,90
vsbw,all,,dele  

accat,69,37  
accat,14,62
accat,66,15
accat,34,38
accat,73,70

et,1,45

mshkey,2
mshape,0,3d
esize,0.5
allsel,all
vmesh,all

D:\ansys work\3.JPG

cj_hmily

D:\ansys work\4.JPGl=10
m=10
h=10
r=3

/prep7  

blc4,0,0,l,m,h  
sphere,0,r,0,90
vsbv,1,2,,dele,dele

cskp,11,0,5,4,1
wpcsys,1,11
vsbw,all,,dele  

cskp,12,0,4,7,1
wpcsys,1,12
vsbw,all,,dele  

cskp,13,0,5,7,1
wpcsys,1,13
vsbw,all,,del

cskp,14,0,4,7,5
wpcsys,1,14
vsbw,all,,dele

et,1,45
mshape,0,3d
mshkey,1
esize,1
vmesh,all

inkey

这贴子收藏,我最怕划分网格了!

inkey

to cj_hmily
给说明一下对应的参数!
cskp,11,0,5,4,1
wpcsys,1,11
vsbw,all,,dele  

cskp,12,0,4,7,1
wpcsys,1,12
vsbw,all,,dele  

cskp,13,0,5,7,1
wpcsys,1,13
vsbw,all,,del

cskp,14,0,4,7,5
wpcsys,1,14
vsbw,all,,dele

et,1,45
mshape,0,3d
mshkey,1
esize,1
vmesh,all

ljk0002

不知这个划的作么样？能否指导指导，划网格要有耐心。

ljk0002

不知这个划的作么样？能否指导指导，划网格要有耐心。
继续

2004woody

真的学了很多。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。