有限元分析中圆、圆柱面及圆柱体的网格划分.docxVIP

有限元分析中圆、圆柱面以及圆柱体的网格划分 简介：有限元分析中网格划分质量决定分析准确性，分析用时，甚至分析对错，掌握经典的几何体的划分是学习有限元的必经之路，本文对圆、圆柱体和圆柱面的网格划分方法给与简介，并给出ANSYS LS-DYNA的例题代码。 关键词：有限元分析；ANSYS；LS-DYNA；网格划分；圆柱体网格划分；圆柱面网格划分 在网上找到ANSYS的圆、圆柱面以及柱划分方法，做了一点修改，改为ANSYS LS-DYNA的划分方法，进行发布。 圆 圆的划分思路是先将圆切分为四份，然后进行划分，划分结果如 REF _Re\h 图1所示： 图 SEQ 图 \* ARABIC 1 圆的网格划分结果 代码如下： finish $ /clear $ /prep7 et,1,plane82 $ r0=10 ! 定义单元类型和圆半径参数 cyl4,,,r0 $ cyl4,3*r0,,,,r0 ! 创建两个圆面 A 和 B，拟分别进行不同的网格划分 wprota,,90 $ asbw,all ! 将圆面水平切分 wprota,,,90 $ asbw,all ! 将圆面 A 竖向切分 wpoff,,,3*r0 $ asbw,all ! 移动工作平面，将圆面 B 竖向切分 wpcsys,-1 ! 工作平面复位但不改变视图方向 asel,s,loc,x,-r0,r0 ! 选择圆面 A 的所有面 lsla,s ! 选择与圆面 A 相关的所有线 lesize,all,,,8 ! 对上述线设置网格划分个数为 8（三条边时相等且为偶数） mshape,0,2d $ mshkey,1 ! 设置四边形单元、映射网格划分 amesh,all ! 圆面 A 划分网格 asel,s,loc,x,2*r0,4*r0 ! 选择圆面B的所有面 lsla,s ! 选择与圆面 B 相关的所有线 lesize,all,,,8 ! 对上述线设置网格划分个数为 8 lsel,r,length,,r0 ! 选择上述线中长度为半径的线 lesize,all,,,8,0.1,1 ! 设置这些线的网格划分数和间隔比 amesh,all $ allsel ! 圆面 B 划分网格 圆柱面 圆柱面的划分结果如 REF _Re\h 图2所示 图 SEQ 图 \* ARABIC 2 圆柱面的网格划分结果 finish /clear /prep7 r0=10 !定义圆半径 h0=50 !定义圆的高度 et,1,shell163 !定义单元类型 cyl4,,,r0 adele,1 !删除面 cm,l1cm,line !将几何元素分组形成组元 k,50 !定义关键点 k,51,,,h0 l,50,51 !利用关键点定义线 adrag,l1cm,,,,,,5 !沿线拉伸成面 lsel,s,loc,z,0 !选择Z=0线 lesize,all,,,6 !对线指定网格尺寸 lsel,s,length,,h0 !选择线 lesize,all,,,8 mshape,0,2d !指定划分单元的形状 mshkey,1 !指定映射网格划分 amesh,all !在面中划分节点或线单元 圆柱体 圆柱用结构化网格划分的思路是将圆柱切分为四份，如 REF _Re\h 图3所示，然后在进行划分，划分效果如 REF _Re\h 图4所