注册 登录  
 加关注
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

htbbzzg的博客

与朋友们分享 CAE 资料和经验

 
 
 

日志

 
 

对 1/4 圆扇形划分指定形式的网格  

2013-07-01 09:39:47|  分类: 问与答 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |

问题:

为了分析的需要,要求对一个 1/4 圆扇形划分指定形式的网格,扇形的中心区域划分为三角形,其余部份为四边形,如下图所示,请问如何做到?

                     对 1/4 圆扇形划分指定形式的网格 - htbbzzg - htbbzzg的博客
回答:

对一个扇形划分指定形状的网格,可以有多种方法,比如:先对扇形的一个边界半径划分线网格 (使用 mesh200 单元),然后将该半径连同网格一起绕圆心旋转,扫描生成所需的面网格。

    这里提供另一种方法。不采用扫描方法,而是通过对扇形区域的适当切割然后对不同区域分别划分不同形状的面单元,同样可以得到指定形状的网格。

 命令流如下,其中可以调整的参数如下:

      rminrmax 扇形区域中,三角形单元区域的半径和扇形的最大半径;

      m - 沿扇形周向划分单元个数;

      n - 沿扇形径向划分四边形单元个数;(三角形单元区域径向仅一个单元)。

fini

/clear

/filn,shanxing_mesh 

!*

WPSTYLE,,,,,,,,0

/PREP7 

!*

ET,1,plane183   !  如果划分 2D 网格,可以使用 183 号单元;

                !  如果要扫描生成 3D 实体单元,可以改用 Mesh200 单元。

!*  材料性能定义,是任意输入的。

MPTEMP,,,,,,,, 

MPTEMP,1,0 

MPDATA,EX,1,,200000

MPDATA,PRXY,1,,0.3 

!*

!  /PNUM,KP,1    !  打开 KP 点号显示

!*

rmin=2           !  三角形单元区域的半径

rmax=30          !  四边形单元区域的最大半径,最小半径为 rmin

m=12             !  生成扇形时,自动将扇形沿周向分为 m 段,

       !  每一段周向只划分一个单元

n=12             !  四边形单元区域的径向线段划分单元个数

                 三角形单元区域的径向线段只划分一个单元

!*

k,1, , ,

k,2, , , 10      !  KP 1 2 用于线旋转生成扇形面

k,3,rmin, , ,  

k,4,rmax, , , 

l,1,3            !  kp 13 连线为三角形单元区域的边界线,

        ! 12 点旋转生成三角形单元区域

l,3,4            !  kp 34 连线为四边形单元区域的边界线,

        ! 12 点旋转生成四边形单元区域

!*

rmin=4           !  修改 rmin rmax,后面用来选择四边形单元区域的半径线

rmax=18

!*

/AUTO,1

/TRIAD,LBOT       !  将坐标系移到屏幕左下方

!*

FLST,2,1,4,ORDE,1  

FITEM,2,1  

FLST,8,2,3 

FITEM,8,1  

FITEM,8,2  

AROTAT,P51X, , , , , ,P51X, ,90,m,    !  绘制 90 扇形,分为 m

!*

FLST,2,1,4,ORDE,1  

FITEM,2,2 

FLST,8,2,3 

FITEM,8,1  

FITEM,8,2  

AROTAT,P51X, , , , , ,P51X, ,90,m,    !  绘制 90 扇形,分为 m

!*

TYPE, 1  

MAT,  1

REAL,  

ESYS, 0  

SECNUM,

!*

/VIEW,1,,,1

/ANG,1 

GPLOT  

!*

LESIZE,ALL, , ,1, ,1, , ,1, 

     !  先设置所有线段只划分一个单元,后面修改四边形区域径向线段的单元分段数

!*

csys,1                         !  转换到圆柱坐标系

LSEL,S,LOC,X,rmin,rmax  

       !  选择四边形单元区域所有径向线段 - 半径介于 rmin rmax 之间

!*

LESIZE,ALL, , ,n, ,1, , ,1,    !  设置四边形单元区域所有径向线段划分单元数为 n

LPLOT  

!*

MSHKEY,0

!*

FLST,5,m,5,ORDE,2           

   !   对三角形区域使用 Free 方式划分三角形网格,面 1m,每个面只有一个单元

FITEM,5,1  

FITEM,5,-m  

CM,_Y,AREA 

ASEL, , , ,P51X

CM,_Y1,AREA

CHKMSH,'AREA'  

CMSEL,S,_Y 

!*

AMESH,_Y1  

!*

CMDELE,_Y  

CMDELE,_Y1 

CMDELE,_Y2 

!*

!*  以下对四边形区域使用 MAP 方式划分四边形网格,

!   面 m+1 2*m,每个面 n 个单元:周向 1 个,径向 n 个。

FLST,5,m,5,ORDE,2  

FITEM,5,m+1  

FITEM,5,-2*m  

CM,_Y,AREA 

ASEL, , , ,P51X

CM,_Y1,AREA

CHKMSH,'AREA'  

CMSEL,S,_Y 

!*

MSHKEY,1

AMESH,_Y1  

MSHKEY,0

!*

CMDELE,_Y  

CMDELE,_Y1 

CMDELE,_Y2 

!*

save

 

示例:rmin=2rmax=30m=12n=12   时所生成的网格:

                      对 1/4 圆扇形划分指定形式的网格 - htbbzzg - htbbzzg的博客

 

  

  评论这张
 
阅读(1448)| 评论(7)
推荐 转载

历史上的今天

在LOFTER的更多文章

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2016