当前位置：文档之家 › ANSYS命令流教程

ANSYS命令流教程

  • 格式：docx
  • 大小：798.28 KB
  • 文档页数：11

引用ANSYS 入门教程(26) - 网格划分高级技术(b)四、单元有效性检查不良的单元形状会导致不准确的结果，然而并没有判别单元形状好坏的通用标准，也就是说一种单元形状对一个分析可能导致不准确的结果，但可能对另一种分析的结果又是可接受的。

在计算过程中，ANSYS可能不出现单元形状警告信息，也可能会出现很多个单元形状警告信息，这都不能说明单元形状就一定会导致准确或不准确的结果，因此单元形状的好坏和结果的准确性完全依赖用户的判断和分析。

1. 单元形状参数限值设置命令：SHPP, Lab, VALUE1, VALUE2ANSYS 单元形状检查是缺省的，但控制单元形状检查的参数可以修改。

⑴Lab = ON：激活单元形状检查。

VALUE1 可取:ANGD：SHELL28 单元角度检查。

ASPECT：单元纵横比检查。

如四边形单元警告限值为20，错误限值为1E6；PARAL：对边平行度检查。

如无中间节点的四边形该项的警告限值为70°，如超过150°则给出错误信息。

MAXANG：最大角度检查。

无中间节点的四边形单元该项警告限值为155°，而其错误限值为179.9°；JACRAT：雅可比率检查。

简单地说，雅可比率表达了“单元”模拟“实际”的计算可靠性，比率越高越不可靠。

如h 单元的警告限值为30，超过30 单元形状就很不理想（与母单元形状相差甚远）。

WARP：歪曲率检查。

对于四边形面单元、壳单元或体单元的面等，当其严重歪曲时造成不好的单元形状，此值越高表示单元歪曲越严重。

也可用ALL 关闭或激活所有选项。

⑵Lab = WARN：仅激活警告模式，对超过错误限制的单元只给出警告信息而不致网格划分失败。

而Lab=ON 则一旦超过错误限制时将导致网格划分失败。

⑶Lab = OFF：完全关闭单元形状检查，可通过设置VALUE1 的值而关闭个别形状检查。

如VALUE1 可取ANGD、ASPECT、PARAL、MAXANG、JACRAT、WARP 及ALL 等。

⑷Lab = STATUS：列表输出当前形状检查限制参数及检查结果情况。

⑸Lab = SUMMARY：列表输出所选择单元的形状检查结果。

⑹Lab = DEFAULT：恢复单元形状检查限值的缺省值。

⑺Lab = OBJECT：是否将单元形状检查结果保存于内存中的控制参数；如VALUE1=1、YES 或ON（缺省）则保存在内存中；如VALUE1=0、NO 或OFF 则不保存在内存中。

⑻Lab = LSTET：检查雅可比率时选择在积分点还是角点取样控制；如VALUE= 1、YES 或ON 则选择积分点；如VALUE1=0、NO 或OFF（缺省）则选择角点取样。

⑼Lab = MODIFY：重新设置一个形状参数检查限值，此时VALUE1 为修改的形状参数限值的数据位置，而VALUE2 则为修改的新限值。

如拟修改纵横比率检查的警告限值，通过SHPP,STATUS 列表可以看出，该数据的位置为1，缺省设置为20.0。

可用SHPP, MODIFY, 1, 1000 将此限值修改为1000。

如拟修改h单元的雅可比率警告限值，通过SHPP,STATUS 列表查得该数据的位置为31，其缺省设置为30.0。

可用SHPP, MODIFY, 31, 100 将此限值修改为100。

使用shpp,defa 将恢复系统的缺省限值设置。

⑽Lab = FLAT：确定显示非零或非常数Z 坐标单元的警告和错误限值。

2. 网格检查命令⑴逐个单元数据完整性检查命令：CHECK, Sele, LevlSele - 拟检查的单元。

如Sele 为空，检查所有单元数据；Levl - 仅当Sele=ESEL 时，其值可取WARN（选择生成警告和错误信息的单元）和ERR（仅检查生成错误信息的单元，这是缺省选项）。

该命令对每个单元的数据完整性和单元形状进行检查，也是在求解之前自动进行的检查。

如单元材料、实常数、约束及单元形状等，然后在输出窗口列出结果。

⑵网格连通性检查命令：MCHECK, Lab其中Lab=ESEL，该选项可不选泽正确的单元，仅选择有问题的单元。

CHECK 命令对单个单元进行检查，而MCHECK 则根据单元的连接方式检查网格潜在的问题，如单元的交叠等。

其检查内容主要有：①方向：当两个面单元共线时，检查每个单元的节点顺序是否与其法线方向一致；②体：当两个体单元共面时，检查每个单元的完整体符号是否一致；③封闭面：检查形成封闭面的单元外表面，以防网格中出现“裂缝”；④网格空洞：当环绕内部空腔的单元面数量很少时，有可能出现遗漏的单元从而形成空洞。

五、网格修改如果对生成的网格不满意，可用下列方法改变网格：①重新设置单元尺寸，并划分网格（只有GUI 才可）；②清除网格，重新设置单元尺寸，并划分网格；③细化局部网格；④改进网格（仅实用于四面体网格）。

1. 清除网格关键点网格清除命令：KCLEAR, NP1, NP2, NINC线网格清除命令：LCLEAR, NL1, NL2, NINC面网格清除命令：ACLEAR, NA1, NA2, NINC体网格清除命令：VCLEAR, NV1, NV2, NINCNX1,NX2,NINC - 为图素范围和编号增量，NX1 可取ALL 或组件名。

该系列命令用于清除既有网格，并可重新对线设置单元网格划分数目或尺寸，然后再重新对几何模型进行网格划分。

2. 细化局部网格节点附近细化命令：NREFINE, NN1, NN2, NINC, LEVEL, DEPTH, POST, RETAIN 单元附近细化命令：EREFINE, NE1, NE2, NINC, LEVEL, DEPTH, POST, RETAIN 关键点附近细化命令：KREFINE, NP1, NP2, NINC, LEVEL, DEPTH, POST, RETAIN 线附近细化命令：LREFINE, NL1, NL2, NINC, LEVEL, DEPTH, POST, RETAIN面附近细化命令：AREFINE, NA1, NA2, NINC, LEVEL, DEPTH, POST, RETAIN Nx1,Nx2,NINC - 图素编号范围与编号增量。

LEVEL - 细化等级，其取值范围1（缺省）～5，值越高网格越密。

当LEVEL=1 时，则采用单元边长的1/2 进行细化生成新的单元。

DEPTH - 从所选图素向外根据单元数设置网格细化的深度，缺省为1。

POST - 单元细化时质量处理控制参数。

如POST=SMOOTH，进行光滑处理，且可能会改变节点位置；如POST=CLEAN（缺省），进行光滑处理，可能会删除存在的单元而重新细分，且节点位置也会改变；如POST=OFF 则不进行任何处理，即节点位置不变也不删除重分。

RETAIN - 所有单元都是四边形网格，在细化时，如RETAIN=ON（缺省）则细化网格也为四边形网格，而不管单元质量如何；如RETAIN=OFF 则允许用三角形网格，以保证网格质量。

但对于下列情况则不能细化：①含有初始条件的节点、耦合节点、约束方程的节点等；②含有边界条件、荷载的节点或单元；③六面体单元、楔形单元和金字塔单元不能细化。

3.4 网格划分实例- 基本模型的网格划分1. 圆圆面的网格划分一般可将圆切分为四等份或八等份，实现映射网格划分。

! EX3.12 圆的网格划分finish $ /clear $ /prep7et,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 ! 对上述线设置网格划分个数为8lsel,r,length,,r0 ! 选择上述线中长度为半径的线lesize,all,,,8,0.1,1 ! 设置这些线的网格划分数和间隔比amesh,all $ allsel ! 圆面B 划分网格2. 圆环圆环面的网格划分与圆面类似，但因由4 条边组成，可更加方便地对网格进行控制。

下面取1/4 圆环面进行单元划分。

! EX3.13 圆环的网格划分finish $ /clear $ /prep7et,1,plane82 $ r0=10 ! 定义单元类型和圆半径参数cyl4,,,r0/3,,r0,90 $ cyl4,2*r0,,r0/10,,r0,90 ! 创建两个1/4 环面asel,s,loc,x,-r0,r0 ! 选择环面Alsla,s$lesize,all,,,8 ! 选择环面A 的所有线，定义网分数lsel,r,length,,r0*2/3 $ lesize,all,,,3,,1 ! 选择径向线，网分数修改为3mshape,0,2d $ mshkey,1 $ amesh,all ! 定义单元形状、划分类型、划分单元ALLSEL $ asel,s,loc,x,2*r0,4*r0 ! 选择环面Blesize,5,,,12 $ lesize,7,,,6 ! 定义外周线和内周线网分数分别为12 和6lsel,s,length,,r0*9/10 $ lesize,all,,,7 ! 选择径向线，网分数为7amesh,all ! 划分环面B 的单元网格3. 圆柱柱体的网格划分方法与圆面类似，空心柱体的网格划分方法同环面类似，而柱面则可直接划分网格。

例如：! EX3.14 圆柱面和圆柱体的网格划分! 圆柱面的网格划分finish $ /clear $ /prep7r0=10 $ h0=50 $ et,1,shell63 ! 定义参数和单元类型cyl4,,,r0 $ adele,1 $ cm,l1cm,line ! 创建面，删除面保留线，定义组件Licmk,50 $ k,51,,,h0 $ l,50,51 $ adrag,l1cm,,,,,,5 ! 创建拖拉路径并拖拉线创建柱面(L1cm - 组件名)lsel,s,loc,z,0 $ lesize,all,,,6 ! 每条圆周线网格划分数为6lsel,s,length,,h0 $ lesize,all,,,8 ! 每条柱面侧线网格划分数为8 mshape,0,2d $ mshkey,1 $ amesh,all ! 定义单元形状、网格划分类型、划分网格!圆柱体finish $ /clear $ /prep7r0=10 $ h0=50 $ et,1,solid95 ! 定义参数和单元类型cyl4,,,r0,,,,h0 ! 创建圆柱体wprota,,90 $ vsbw,all ! 切分圆柱体wprota,,,90 $ vsbw,all ! 再切分圆柱体mshape,0,3d $ mshkey,1 ! 定义单元形状、网格划分类型lsel,s,loc,z,0 $ lesize,all,,,6 ! 每条圆周线网格划分数为6lsel,s,length,,h0 $ lesize,all,,,8 ! 每条柱面侧线网格划分数为8vmesh,all4. 圆锥和圆台圆锥体的网格划分应以1/4 圆锥进行，然后利用对称命令创建其余体和网格。

合集下载

ANSYS APDL命令流详解-4创建几何模型

ANSYS APDL命令流详解-4创建几何模型

。关键点必须按顺时针或逆时针顺序输入，同时按右手规则确 定面的正法线方向。
当关键点数≥4时，应该保证所有关键点位于同一平面或曲面 内，即在当前坐标系下有一相同的坐标值，如Z相同，则该面位 于XY平面内。
★如果相邻两关键点已经存在线（直线或曲线），则创建面时 使用该线，该线形状与当前坐标系无关；
★如果存在多条线，则采用其中最短的线（直线）。 ★如果相邻关键点没有线，则创建面时边的形状决定当前坐标 系，如在直角坐标系下生成直线边，而在柱坐标系下生成曲线 边。但是一旦由这些关键点创建了面，再改变当前坐标系也不 能改变面的形状了。
NL1~NL6---将要拖拉的线号，也可为ALL或元件名，线必须是连续的。 NLP1~NLP6---路径线的编号，也必须是连续的。也可为元件名。
★用ADRAG创建的面， 其线和关键点号由系统自 动定义
★相邻面共用线、相邻线 共用关键点。
★拖拉线与拖拉路径不一 定相交，拖拉线仅仅将路 径作为方向和参考长度， 该命令在创建复杂曲面时 较为方便。
2. 通过线创建面 命令：AL, L1, L2, L3, L4, L5, L6, L7, L8, L9, L10 其中L1～L10为线编号，最少要3条线，当采用输入线号时最
多10条线。生成面的正法线方向按右手规则由L1的方向确定。 当L1为负值时则表示面的正法线方向相反。L1可为ALL、P或元 件名，当L1=ALL时面的法线由L2定义面的法线方向，当L2为空 时则默认为最小编号的线，且此时线数不受限制。
K,50,,,2
!创建关键点及线
K,51,,1,4$L,1,50$L,50,51
ADRAG,11,12,,,,,path1 !沿路径PATH1拖拉线L11和L12创建面
4. 线绕轴旋转生成弧面
（完整版）ANSYS最常用命令流+中文注释（超级大全）

（完整版）ANSYS最常用命令流+中文注释（超级大全）

（完整版）ANSYS最常用命令流+中文注释（超级大全）ANSYS最常用命令流+中文注释VSBV, NV1, NV2, SEPO, KEEP1, KEEP2 —Subtracts volumes from volumes，用于2个solid相减操作，最终目的是要nv1-nv2=?通过后面的参数设置，可以得到很多种情况：sepo项是2个体的边界情况，当缺省的时候，是表示2个体相减后，其边界是公用的，当为sepo的时候，表示相减后，2个体有各自的独立边界。

keep1与keep2是询问相减后，保留哪个体？当第一个为keep时，保留nv1,都缺省的时候，操作结果最终只有一个体，比如：vsbv,1,2,sepo,,keep,表示执行1-2的操作，结果是保留体2，体1被删除，还有一个1-2的结果体，现在一共是2个体（即1-2与2），且都各自有自己的边界。

如vsbv,1,2,,keep,,则为1-2后，剩下体1和体1-2，且2个体在边界处公用。

同理，将v换成a 及l是对面和线进行减操作！mp,lab, mat, co, c1,…….c4 定义材料号及特性lab: 待定义的特性项目（ex,alpx,reft,prxy,nuxy,gxy,mu,dens）ex: 弹性模量nuxy: 小泊松比alpx: 热膨胀系数reft: 参考温度reft: 参考温度prxy: 主泊松比gxy: 剪切模量mu: 摩擦系数dens: 质量密度mat: 材料编号（缺省为当前材料号）co: 材料特性值，或材料之特性，温度曲线中的常数项c1-c4: 材料的特性-温度曲线中1次项，2次项，3次项，4次项的系数定义DP材料：首先要定义EX和泊松比：MP，EX，MA T，……MP，NUXY，MAT，……定义DP材料单元表（这里不考虑温度）：TB，DP，MA T进入单元表并编辑添加单元表：TBDATA，1，CTBDATA，2，ψTBDATA，3，……如定义：EX=1E8，NUXY=0.3，C=27，ψ=45的命令如下：MP，EX，1，1E8MP，NUXY，1，0.3TB，DP，1TBDATA，1，27TBDATA，2，45这里要注意的是，在前处理的最初，要将角度单位转化到“度”，即命令：*afun,degVSEL, Type, Item, Comp, VMIN, VMAX, VINC, KSWP Type，是选择的方式，有选择（s）,补选（a），不选(u)，全选(all)、反选（inv）等,其余方式不常用Item, Comp 是选取的原则以及下面的子项如volu 就是根据实体编号选择，loc 就是根据坐标选取，它的comp就可以是实体的某方向坐标！其余还有材料类型、实常数等MIN, VMAX, VINC，这个就不必说了吧！,例：vsel,s,volu,,14vsel,a,volu,,17,23,2上面的命令选中了实体编号为14，17，19，21，23的五个实体VDELE, NV1, NV2, NINC, KSWP: 删除未分网格的体nv1:初始体号nv2:最终的体号ninc:体号之间的间隔kswp=0:只删除体kswp=1:删除体及组成关键点,线面如果nv1=all,则nv2,ninc不起作用其后面常常跟着一条显示命令VPLO,或aplo,nplo,这个湿没有参数的命令，输入后直接回车，就可以显示刚刚选择了的体、面或节点，很实用的哦！Nsel, type, item, comp, vmin, vmax, vinc, kabs 选择一组节点为下一步做准备Type: S: 选择一组新节点（缺省）R: 在当前组中再选择A: 再选一组附加于当前组U: 在当前组中不选一部分All: 恢复为选中所有None: 全不选Inve: 反向选择Stat: 显示当前选择状态Item: loc: 坐标node: 节点号Comp: 分量Vmin,vmax,vinc: ITEM范围Kabs: “0” 使用正负号“1”仅用绝对值下面是单元生死第一个载荷步中命令输入示例：!第一个载荷步TIME,... !设定时间值（静力分析选项）NLGEOM,ON !打开大位移效果NROPT,FULL !设定牛顿－拉夫森选项ESTIF,... !设定非缺省缩减因子（可选）ESEL,... !选择在本载荷步中将不激活的单元EKILL,... !不激活选择的单元ESEL,S,LIVE !选择所有活动单元NSLE,S !选择所有活动结点NSEL,INVE !选择所有非活动结点（不与活动单元相连的结点）D,ALL,ALL,0 !约束所有不活动的结点自由度（可选）NSEL,ALL !选择所有结点ESEL,ALL !选择所有单元D,... !施加合适的约束F,... !施加合适的活动结点自由度载荷SF,... !施加合适的单元载荷BF,... !施加合适的体载荷SA VESOLVE请参阅TIME,NLGEOM,NROPT,ESTIF,ESEL,EKILL,NSLE,NSEL,D, F,SF和BF命令得到更详细的解释。

ANSYS APDL命令流详解-10网格划分高级技术

ANSYS APDL命令流详解-10网格划分高级技术

1. 体扫掠器的激活 命令：VSWEEP,VNUM,SRCA,TRGA,LSMO VNUM---体的编号，还可取ALL、P及元件名或组件名。 SRCA---源面编号。如果该源面尚未划分网格，则系统自动对
其划分网格然后再扫掠。如果不指定SRCA则由系统 自动确定源面。 TRGA---目标面编号，即SRCA面的对面。如果不指定该面 号，系统自动确定目标面。 LSMO---在扫掠时线光滑处理控制参数。
!定义单元类型，创建六面体 !所有线均划分4个分段 !线7定义12个分段
mshape,0,3d$mshkey,1 !单元形状和划分类型定义
vmesh,all
!划分网格
3.3网格划分高级技术- 扫掠生成体网格
对于3D几何体，除采用自由网格划分和映射网格划分外，还可采用“扫 掠(sweep)网格划分”，体扫掠网格划分就是从源面（如边界面）网格扫掠 整个体生成体单元。如果源面网格由四边形网格组成，则扫掠生成的均为 六面体单元；如果源面网格由三角形网格组成，则扫掠生成的均为楔形体 单元；如果源面网格由四边形和三角形网格组成，则扫掠生成六面体和楔 形体单元。
3.3网格划分高级技术-体映射网格划分
2. 连接面和面加运算
对面进行加运算（AADD命令）或连接（ACCAT命令），以 减少面数从而达到体映射网格划分的条件。
连接面时，如果连接面有边界线，线也必须连接在一起，并 且必须先连接面，再连接线。但是如果相连接的两个面都由四条 边组成，线的连接操作会自动进行。当删除连接面时并不自动删 除相关的连接线，应用LDELE命令删除连接线。
3.3网格划分高级技术-面映射网格划分
4. 简化面映射网格划分AMAP 由于上述操作较为复杂，由AMAP命令采用指定的关键点作
为角点，不需要MSHKEY命令参数，自动地进行面的网格划 分
ANSYS命令流学习笔记workbench中命令流的一些应用

ANSYS命令流学习笔记workbench中命令流的一些应用

!ANSYS命令流学习笔记5workbench中命令流的一些应用学习重点：1. 定义单元类型2. 使用各向异性材料时，定义其单元为圆柱坐标系3. 有角度吊装时，定义吊装约束4. workbench、APDL的联合仿真案例如下：如下图模型，四个顶点通过杆件连接，进行吊装时的有限元分析。

1. 建立模型，设定必要的坐标系。

分成两个solid，内部的圆柱，剩余的矩形部分。

建立两个坐标系，分别用于指定各向异性材料的属性、吊装的固定点。

下图，建立圆柱坐标系，编号100，用于指定各向异性材料。

下图，建立直角坐标系，编号12，用于指定吊装固定点。

2. 建立named selection，方便在命令流中选择必要的元素。

下图，将四个吊装点，中间的圆柱，分别定义为任何名称，必须是英文才能用于APDL命令中。

3. 定义边界条件，施加重力加速度，在static structural 下插入command(APDL) ，内容如下/prep7allscmsel,s,c1,elem !选择c1单元所有节点,既圆柱体的所有单元emodif,all,esys,100 !其坐标系转换为100坐标系，因为缠绕的各向异性材料必须在圆柱坐标系下定义单元的坐标。

!完成对各向异性材料的坐标系设定。

et,10,10 !定义编号为10的，link10单元r,10,0.01 !定义编号为10的实常数0.01，用于定义link10单元的截面积0.01mm^2*get,nmax,node,,num,max !获取node的最大数值，储存在nmax名称的变量里csys,12 !调用csys12坐标系n,nmax+1, !csys坐标原点建立node，后面会将其固定，既吊装的固定点mat,1type,10 !选取编号10的单元类型real,10 !选取编号10的实常数cmsel,s,kk1,node !选择kk1点，kk1已经定义为named selection*get,k1,node,,num,max !获取已选节点中的节点数最大值，既kk1的节点编号，取值为k1cmsel,s,kk2,node ! k2,k3,k4方法类型k1*get,k2,node,,num,maxcmsel,s,kk3,node*get,k3,node,,num,maxcmsel,s,kk4,node*get,k4,node,,num,maxalls !全选e,nmax+1,k1e,nmax+1,k2e,nmax+1,k3e,nmax+1,k4 !建立四个link单元d,nmax+1,allalls/sol此外，下图中单位必须保持一致，不然计算很容易出问题。