第3章网格划分剖析

《第3章网格划分剖析》由会员分享，可在线阅读，更多相关《第3章网格划分剖析（118页珍藏版）》请在金锄头文库上搜索。

1、授课教师： 叶 茂 副研究员 工程北320,第3章 网格划分,目 录 3.1 定义单元属性 单元类型 /实常数/材料属性 /梁截面/设置几何模型的单元属性 3.2 网格划分控制 单元形状控制及网格类型选择/单元尺寸控制 内部网格划分控制/划分网格 3.3 网格划分高级技术 面映射网格划分/体映射网格划分 扫掠生成体网格/单元有效性检查/网格修改 3.4 网格划分实例 基本模型的网格划分/复杂面模型的网格划分 复杂体模型的网格划分,创建几何模型后，必须生成有限元模型才能分析计算，生成有限元模型的方法就是对几何模型进行网格划分，网格划分主要过程包括三个步骤： 定义单元属性 单元属性包括单元类型、实

2、常数、材料特性、单元坐标系和截面号等。 定义网格控制选项 对几何图素边界划分网格的大小和数目进行设置； 没有固定的网格密度可供参考； 可通过评估结果来评价网格的密度是否合理。 生成网格 执行网格划分，生成有限元模型； 可清除已经生成的网格并重新划分； 局部进行细化。,3.1 定义单元属性单元类型,定义单元类型 命令：ET, ITYPE, Ename, KOP1, KOP2, KOP3, KOP4, KOP5, KOP6, INOPR ITYPE-用户定义的单元类型的参考号。 Ename-ANSYS单元库中给定的单元名或编号，它由一个类别前缀和惟一的编号组成，类别前缀可以省略，而仅使用单元编号。

3、 KOP1KOP6-单元描述选项，此值在单元库中有明确的定义，可参考单元手册。也可通过命令KEYOPT进行设置。 INOPR-如果此值为1则不输出该类单元的所有结果。 例如： et,1,link8 !定义LINK8单元，其参考号为1；也可用ET,1,8定义 et,3,beam4 !定义BEAM4单元，其参考号为3；也可用ET,3,4定义,2. 单元类型的KEYOPT 命令：KEYOPT,ITYPE,KNUM,VALUE ITYPE-由ET命令定义的单元类型参考号。 KNUM-要定义的KEYOPT顺序号。 VALUE-KEYOPT值。 该命令可在定义单元类型后，分别设置各类单元的KEYOPT参数

4、。例如： et,1,beam4 !定义BEAM4单元的参考号为1 et,3,beam189 !定义BEAM189单元的参考号为3 keyopt,1,2,1 !BEAM4单元考虑应力刚度时关闭一致切线刚度矩阵 keyopt,3,1,1 !考虑BEAM189的第7个自由度，即翘曲自由度 !当然这些参数也可在ET命令中一并定义，如上述四条命令与下列两条命令等效： et,1,beam4,1 et,3,beam189,1,3.1 定义单元属性单元类型,3. 自由度集 命令：DOF,Lab1,Lab2,Lab3,Lab4,Lab5,Lab6,Lab7,Lab8,Lab9,Lab10 4. 改变单元类型 命

5、令：ETCHG,Cnv 5. 单元类型的删除与列表 删除命令：ETDELE,ITYP1,ITYP2,INC 列表命令：ETLIST,ITYP1,ITYP2,INC,3.1 定义单元属性单元类型,3.1 定义单元属性实常数,1. 定义实常数 命令：R,NSET,R1,R2,R3,R4,R5,R6 续：RMORE,R7,R8,R9,R10,R11,R12 NSET-实常数组号（任意），如果与既有组号相同，则覆盖既有组号定义的实常数。 R1R12-该组实常数的值。 使用R命令只能一次定义6个值，如果多于6个值则采用RMORE命令增加另外的值。每重复执行RMORE一次，则该组实常数增加6个值，如712

6、、1318、1924等。 各类单元有不同的实常数值，其值的输入必须按单元说明中的顺序； 如果实常数值多于单元所需要的，则仅使用需要的值；如果少于所需要的，则以零值补充。 一种单元可有多组实常数，也有单元不需要实常数的。,3.1 定义单元属性实常数,例如BEAM4单元，需要的实常数值有12个： AREA、IZZ、IYY、TKZ、TKY、THETA和 ISTRN、IXX、SHEARZ、SHEARY、SPIN、ADDMAS 设采用直径为0.1m的圆杆，其实常数可定义为： D=0.1 PI=acos(-1) a0=pi*d*d/4 I0=pi*D*4/64 IX=pi*D*4/32 R,3,a0,i0

7、,i0,d,d,0 !定义第3组实常数的AREA、IZZ、IYY、TKZ、TKY、THETA Rmore,0,ix,0,0,0,2.0 !定义第3组实常数的其它实常数值,3.1 定义单元属性实常数,2. 变厚度壳实常数定义 命令：RTHICK,Par,ILOC,JLOC,KLOC,LLOC Par-节点厚度的数组参数（以节点号引用），如mythick(19)表示在节点19的壳体厚度。 ILOC-单元I节点的厚度在实常数组中的位置，缺省为1。 JLOC-单元J节点的厚度在实常数组中的位置，缺省为2。 KLOC-单元K节点的厚度在实常数组中的位置，缺省为3。 LLOC-单元L节点的厚度在实常数组中

8、的位置，缺省为4。 该命令后面的四个参数顺序与节点厚度的关系比较复杂。 3. 实常数组的删除与列表 删除命令：RDELE, NSET1,NSET2,NINC 列表命令：RLIST, NSET1,NSET2,NINC 其中NSET1,NSET2,NINC-实常数组编号范围和编号增量，缺省时NSET2等于NSET1且NINC=1。NSET1也可为ALL 。,3.1 定义单元属性材料属性,每一组材料属性有一个材料参考号，用于识别各个材料特性组。一个模型中可有多种材料特性组。 1. 定义线性材料属性 命令：MP,Lab,MAT,C0,C1,C2,C3,C4 Lab-材料性能标识，其值可取： EX：弹性

9、模量（也可为EY、EZ）。 ALPX：线膨胀系数（也可为ALPY、ALPZ ）。 PRXY：主泊松比（也可为PRYZ、PRXZ）。 NUXY：次泊松比（也可为NUYZ、NUXZ）。 GXY：剪切模量（也可为GYZ、GXZ）。 DAMP：用于阻尼的K矩阵乘子，即阻尼比。 DMPR：均质材料阻尼系数。 MU：摩擦系数。 DENS：质量密度。 MAT-材料参考号，缺省为当前的MAT号（由MAT命令确定）。 C0-材料属性值，如果该属性是温度的多项式函数，则此值为多项式的常数项。 C1C4-分别为多项式中的一次、二次、三次、四次项系数，如为0或空，则定义一个常数的材料性能,3.1 定义单元属性材料属性

10、,2. 定义线性材料属性的温度表 命令：MPTEMP,STLOC,T1,T2,T3,T4,T5,T6 3. 定义与温度对应的线性材料特性 命令：MPDATA,Lab,MAT,STLOC,C1,C2,C3,C4,C5,C6 4. 复制线性材料属性组 命令：MPCOPY,-,MATF,MATT 5. 改变指定单元的材料参考号 命令：MPCHG,MAT,ELEM 6. 线性材料属性列表和删除 列表命令：MPLIST,MAT1,MAT2,INC,Lab,TEVL 删除命令：MPDELE,Lab,MAT1,MAT2,INC 7. 修改与线胀系数相关的温度 命令：MPAMOD,MAT,DEFTEMP 8.

11、 计算生成线性材料温度表 命令：MPTGEN,STLOC,NUM,TSTRT,TINC 9. 绘制线性材料特性曲线 命令：MPPLOT,Lab,MAT,TMIN,TMAX,PMIN,PMAX,3.1 定义单元属性材料属性,10. 设置材料库读写的缺省路径 命令：/MPLIB,R-W_opt,PATH 11. 读入材料库文件 命令：MPREAD,Fname,Ext,-,LIB 12. 将材料属性写入文件 命令：MPWRITE,Fname,Ext,-,LIB,MAT 13. 激活非线性材料属性的数据表 命令：TB,Lab,MAT,NTEMP,NPTS,TBOPT,EOSOPT 14. 定义TB温度

12、值 命令：TBTEMP,TEMP,KMOD 15. 定义TB数据表中的数据 命令：TBDATA,STLOC,C1,C2,C3,C4,C5,C6 16. 定义非线性数据曲线上的一个点 命令：TBPT,Oper,X,Y 17. 非线性材料数据表的删除和列表 删除命令：TBDELE,Lab,MAT1,MAT2,INC 列表命令：TBLIST,Lab,MAT 18. 非线性材料数据表的绘图 命令：TBPLOT,Lab,MAT,TBOPT,TEMP,SEGN,3.1 定义单元属性 梁截面,BEAM18x单元，需定义单元的横截面（称为梁截面）； BEAM44也可使用梁截面也可输入截面特性实常数； 仅BEA

13、M18x可使用多种材料组成的截面； 仅BEAM18x可使用变截面梁截面，而BEAM44可输入实常数。 1. 定义截面类型和截面ID 命令：SECTYPE,SECID,Type,Subtype,Name,REFINEKEY SECID-截面识别号，也称为截面ID号。 Type-截面用途类型，其值可取： BEAM：定义梁截面，应用于等截面时，见下文。 TAPER：定义渐变梁截面（变截面梁）。 SHELL：定义壳 PRETENSION：定义预紧截面 JOINT：连接截面，如万向铰。,3.1 定义单元属性 梁截面,Subtype-截面类型，对于不同的Type该截面类型不同，如： 当Type=BEAM时

14、，Subtype可取： RECT：矩形截面； QUAD：四边形截面； CSOLID：实心圆形截面； CTUBE：圆管截面； CHAN：槽形截面； I：工字形截面； Z：Z形截面； L：L形截面； T：T形截面； HATS：帽形截面； HREC：空心矩形或箱形 ASEC：任意截面； MESH：自定义截面 当Type=JOINT （有刚度可大角度旋转）时，Subtype可取： UNIV：万向铰； REVO：销铰或单向铰 Name-8个字符的截面名，字符可包含字母和数字。 REFINEKEY-设置薄壁梁截面网格的精细水平，有0（缺省）5（最精细）六个水平。,3.1 定义单元属性 梁截面,2. 定义梁

15、截面几何数据（Type=BEAM） 命令：SECDATA,VAL1,VAL2,VAL3,VAL4,VAL5,VAL6,VAL7,VAL8,VAL9,VAL10 其中VAL1VAL10为数值，如厚度、边长、沿边长的栅格数等， 每种截面的值是不同的。 ANSYS定义了11种常用的截面类型，每种截面输入数据如下： Subtype=RECT：矩形截面 输入数据：B,H,Nb,Nh B-截面宽度。 H-截面高度。 Nb-沿宽度B的栅格数（cell），缺省为2。 Nh-沿高度H的栅格数，缺省为2。,3.1 定义单元属性 梁截面,Subtype=QUAD：四边形截面 输入数据：yI,zI,yJ,zJ,yK,

16、zK,yL,zL,Ng,Nh yI,zI,yJ,zJ,yK,zK,yL,zL-各点坐标值。 Ng,Nh-沿g和h的栅格数，缺省均为2。 如退化为三角形也可，输入一个相同的坐标。,3.1 定义单元属性 梁截面,Subtype=CSOLID：实心圆截面 输入数据：R,N,T R-半径。 N-圆周方向划分的段数，缺省为8。 T-半径方向划分的段数，缺省为2。,3.1 定义单元属性 梁截面,Subtype=CHAN：槽形截面 输入数据：W1,W2,W3,t1,t2,t3 W1,W2-翼缘宽度。 W3-全高。 t1,t2-翼缘厚度。 t3-腹板厚度,3.1 定义单元属性 梁截面,Subtype=I：工字形截面 输入数据：W1,W2,W3,t