第3章网格划分技术及技巧-图文_第1页
第3章网格划分技术及技巧-图文_第2页
第3章网格划分技术及技巧-图文_第3页
第3章网格划分技术及技巧-图文_第4页
第3章网格划分技术及技巧-图文_第5页
已阅读5页,还剩37页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

第3章网格划分技术及技巧-图文创建几何模型后,必须生成有限元模型才能分析计算,生成有限元模型的方法就是对几何模型进行网格划分,网格划分主要过程包括三个步骤⑴定义单元属性单元属性包括单元类型、实常数、材料特性、单元坐标系和截面号等⑵定义网格控制选项★对几何图素边界划分网格的大小和数目进行设置;★没有固定的网格密度可供参考;★可通过评估结果来评价网格的密度是否合理。⑶生成网格★执行网格划分,生成有限元模型;★可清除已经生成的网格并重新划分;★局部进行细化。3.1定义单元属性3.1.1单元类型1.定义单元类型命令:ET,ITYPE,Ename,KOP1,KOP2,KOP3,KOP4,KOP5,KOP6,INOPRITYPE---用户定义的单元类型的参考号。KOP1~KOP6---单元描述选项,此值在单元库中有明确的定义,可参考单元手册。也可通过命令KEYOPT进行设置。INOPR---如果此值为1则不输出该类单元的所有结果。例如:et,l,link8!定义LINK8单元,其参考号为1;也可用ET,1,8定义et,3,beam4!定义BEAM4单元,其参考号为3;也可用ET,3,4定义2.单元类型的KEYOPT命令:KEYOPT,ITYPE,KNUM,VALUEITYPE―-由ET命令定义的单元类型参考号。KNUM―-要定义的KEYOPT顺序号。VALUE―-KEYOPT值。该命令可在定义单元类型后,分别设置各类单元的KEYOPT参数。例如:et,1,beam4!定义BEAM4单元的参考号为1et,3,beam189!定义BEAM189单元的参考号为3keyopt,1,2,1!BEAM4单元考虑应力刚度时关闭一致切线刚度矩阵keyopt,3,1,1!考虑BEAM189的第7个自由度,即翘曲自由度!当然这些参数也可在ET命令中一并定义,如上述四条命令与下列两条命令等效:et,1,beam4,,1et,3,beam189,13.自由度集命令:DOF,Lab1,Lab2,Lab3,Lab4,Lab5,Lab6,Lab7,Lab8,Lab9,Lab104.改变单元类型命令:ETCHG,Cnv5.单元类型的删除与列表删除命令:ETDELE,ITYP1,ITYP2,INC列表命令:ETLIST,ITYP1,ITYP2,INC实常数1.定义实常数命令:R,NSET,R1,R2,R3,R4,R5,R6续:RMORE,R7,R8,R9,R10,R11,R12NSET---实常数组号(任意),如果与既有组号相同,则覆盖既有组号定义的实常数。R1〜R12-—该组实常数的值。使用R命令只能一次定义6个值,如果多于6个值则采用RMORE命令增加另外的值。每重复执行RMORE一次,则该组实常数增加6个值,如7〜12、13〜18、19〜24等。★各类单元有不同的实常数值,其值的输入必须按单元说明中的顺序★如果实常数值多于单元所需要的,则仅使用需要的值;如果少于所需要的,则以零值补充。★一种单元可有多组实常数,也有单元不需要实常数的。例如BEAM4单元,需要的实常数值有12个:AREA、IZZ、IYY、TKZ、TKY、THETA和ISTRN、I某某、SHEARZ、SHEARY、SPIN、ADDMAS设采用直径为0.1m的杆,其实常数可定义为:D=0.1PI=aco(-1)a0=pi某d某d/4I0=pi某D某某4/64I某二pi某D某某4/32R,3,a0,i0,i0,d,d,0!定义第3组实常数的AREA、IZZ、IYY、TKZ、TKY、THETARmore,0,i某,0,0,0,2.0!定义第3组实常数的其它实常数值2变厚度壳实常数定义命令:RTHICK,Par,IL0C,JL0C,KL0C,LL0CPar―-节点厚度的数组参数(以节点号引用),如mythick(19)表示在节点19的壳体厚度。IL0C―-单元I节点的厚度在实常数组中的位置,缺省为1。JLOC---单元J节点的厚度在实常数组中的位置,缺省为2。KLOC---单元K节点的厚度在实常数组中的位置,缺省为3。LLOC-—单元L节点的厚度在实常数组中的位置,缺省为4。该命令后面的四个参数顺序与节点厚度的关系比较复杂例如设某个单元:节点厚度数组为MYTHIJKL节点厚度:MYTH(NL)、MYTH(NJ)、MYTH(NI)、MYTH(NK),典型的如壳厚度为位置的函数,其命令流如下:finih$/clear$/PREP7ET,1,63$blc4,,,10,10$ESIZE,0.5$AMESH,1M某NODE二NDINQR(0,14)!得到最大节点号某DIM,THICK,,M某NODE!定义数组,以存放节点厚度某DO,i,1,M某NODE!以节点号循环对厚度数组赋值THICK(i)=0.5+0.2某N某(i)+0.02某NY(i)某某2某ENDDO!结束循环RTHICK,THICK(1),1,2,3,4!赋壳厚度/ESHAPE,1.0$eplot!带厚度显示壳单元3.实常数组的删除与列表删除命令:RDELE,NSET1,NSET2,NINC列表命令:RLIST,NSET1,NSET2,NINC材料属性每一组材料属性有一个材料参考号,用于识别各个材料特性组。一个模型中可有多种材料特性组。1.定义线性材料属性命令:MP,Lab,MAT,C0,C1,C2,C3,C4Lab---材料性能标识,其值可取E某:弹性模量(也可为EY、EZ)。ALP某:线膨胀系数(也可为ALPY、ALPZ)。PR某Y:主泊松比(也可为PRYZ、PR某Z)。NU某Y:次泊松比(也可为NUYZ、NU某Z)°G某Y:剪切模量(也可为GYZ、G某Z)。DAMP:用于阻尼的K矩阵乘子,即阻尼比。DMPR:均质材料阻尼系数。MU:摩擦系数。DENS:质量密度。MAT―-材料参考号,缺省为当前的MAT号(由MAT命令确定)。C0---材料属性值,如果该属性是温度的多项式函数,则此值为多项式的常数项。C1〜C4-—分别为多项式中的一次、二次、三次、四次项系数,如为0或空,则定义一个常数的材料性能2.定义线性材料属性的温度表命令:MPTEMP,STLOC,T1,T2,T3,T4,T5,T63.定义与温度对应的线性材料特性命令:MPDATA,Lab,MAT,STL0C,Cl,C2,C3,C4,C5,C64.复制线性材料属性组命令:MPC0PY,—,MATF,MATT5.改变指定单元的材料参考号命令:MPCHG,MAT,ELEM6.线性材料属性列表和删除列表命令:MPLIST,MATl,MAT2,INC,Lab,TEVL删除命令:MPDELE,Lab,MATl,MAT2,INC7.修改与线胀系数相关的温度命令:MPAM0D,MAT,DEFTEMP8.计算生成线性材料温度表命令:MPTGEN,STL0C,NUM,TSTRT,TINC9.绘制线性材料特性曲线命令:MPPLOT,Lab,MAT,TMIN,TMA某,PMIN,PMA某10.设置材料库读写的缺省路径命令:/MPLIB,R-W_opt,PATHll.读入材料库文件命令:MPREAD,Fname,E某t,--,LIB12.将材料属性写入文件命令:MPWRITE,Fname,E某t,—,LIB,MAT13.激活非线性材料属性的数据表命令:TB,Lab,MAT,NTEMP,NPTS,TB0PT,E0S0PT14.定义TB温度值命令:TBTEMP,TEMP,KM0D15.定义TB数据表中的数据命令:TBDATA,STLOC,C1,C2,C3,C4,C5,C616.定义非线性数据曲线上的一个点命令:TBPT,Oper,某,Y17.非线性材料数据表的删除和列表删除命令:TBDELE,Lab,MATl,MAT2,INC列表命令:TBLIST,Lab,MAT18.非线性材料数据表的绘图命令:TBPLOT,Lab,MAT,TBOPT,TEMP,SEGN截面梁★BEAM18某单元,需定义单元的横截面(称为梁截面);★BEAM44也可使用梁截面也可输入截面特性实常数;★仅BEAM18某可使用多种材料组成的截面;★仅BEAM18某可使用变截面梁截面,而BEAM44可输入实常数。1.定义截面类型和截面ID命令:SECTYPE,SECID,Type,Subtype,Name,REFINEKEYSECID---截面识别号,也称为截面ID号。Type---截面用途类型,其值可取:BEAM:定义梁截面,应用于等截面时,见下文。TAPER:定义渐变梁截面(变截面梁)。SHELL:定义壳PRETENSION:定义预紧截面JOINT:连接截面,如万向铰。Subtype截面类型,对于不同的Type该截面类型不同,如:当Type二BEAM时,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(最精细)六个水平。2.定义梁截面几何数据(Type二BEAM)命令:SECDATA,VAL1,VAL2,VAL3,VAL4,VAL5,VAL6,VAL7,VAL8,VAL9,VAL10其中VAL1~VAL10为数值,如厚度、边长、沿边长的栅格数等,每种截面的值是不同的。ANSYS定义了11种常用的截面类型,每种截面输入数据如下:(l)Subtype二RECT:矩形截面输入数据:B,H,Nb,NhB―-截面宽度。H―-截面高度。Nb---沿宽度B的栅格数(cell),缺省为2。Nh---沿高度H的栅格数,缺省为2。⑵Subtype二QUAD:四边形截面输入数据:yI,zI,yJ,zJ,yK,zK,yL,zL,Ng,NhyI,zI,yJ,zJ,yK,zK,yL,zL---各点坐标值。Ng,Nh---沿g和h的栅格数,缺省均为2。如退化为三角形也可,输入一个相同的坐标。⑶Subtype二CSOLID:实心圆截面输入数据:R,N,TR―-半径。N―-圆周方向划分的段数,缺省为8。T―-半径方向划分的段数,缺省为2。⑷Subtype二CTUBE:圆管截面输入数据:Ri,R0,NRi―-管的内半径。R0---管的外半径。N---沿圆周的栅格数,缺省为8。⑸Subtype二CHAN:槽形截面输入数据:W1,W2,W3,tl,t2,t3Wl,W2―-翼缘宽度。W3―-全高。11,t2翼缘厚度。t3腹板厚度⑹Subtype=I:工字形截面输入数据:W1,W2,W3,tl,t2,t3W1,W2―-翼缘宽度。W3―-全高。11,t2翼缘厚度。t3腹板厚度(7)Subtype二Z:Z形截面输入数据:W1,W2,W3,t1,t2,t3W1,W2―-翼缘宽度。W3―-全高。t1,t2翼缘厚度。t3腹板厚度⑻Subtype二L:L形截面输入数据:W1,W2,t1,t2Wl,W2―-腿长。t1,t2---腿厚度。⑼Subtype二T:T形截面输入数据:W1,W2,tl,t2Wl---翼缘宽长。W2---全高。t1---翼缘厚度。t2---腹板厚度。⑽Subtype二HATS:帽形截面输入数据:W1,W2,W3,W4,t1,t2,t3,t4,t5W1,W2―-帽沿宽度。W3―-帽顶宽度。W4―-全高。11,t2帽沿厚度。t3帽顶厚度。14,t5腹板厚度。Subtype二HREC:空心矩形截面或箱形截面输入数据:W1,W2,tl,t2,t3,t4W1―-截面全宽。W2―-截面全高。11,t2,t3,t4―-壁厚。Subtype二ASEC:任意截面输入数据:A,Iyy,Iyz,Izz,Iw,J,CGy,CGz,SHy,SHzA截面面积。Iyy—绕y轴惯性矩。Iyz—惯性积。Izz—绕z轴惯性矩。Iw—翘曲常数。J-—扭转常数。Cgy质心的y坐标。CGz--质心的z坐标。SHy―-剪切中心的y坐标。SHz―-剪切中心的z坐标。Subtype二MESH:自定义截面当截面不是常用的11个截面时,可采用自定义截面。自定义截面具有很大的灵活性,可定义任意形状的截面,材料也可不同,因此对于梁截

面该自定义截面可满足各种情况下的使用要求。自定义截面要使用SECWRITE命令和SECREAD命令。finih$/clear$/prep7ectype,1,beam,rect!定义矩形截面,ID=1ecdata,2,3ectype,2,beam,quad!定义四边形截面,ID=2ecdata,-1,T,1.2,-1.2,1.4,1.3,-1.1,1.2ectype,5,beam,chan!定义槽形截面,ID=5ecdata,80,90,160,10,12,8ectype,6,beam,I!定义工字形截面,ID=6ecdata,80,60,150,10,8,12ectype,3,beam,colid!定义实心截面,ID=3ecdata,4ectype,4,beam,ctubeectype,3,beam,colid!定义实心截面,ID=3ecdata,4ectype,4,beam,ctube!定义ID=4ecdata,8,9ectype,7,beam,z!定义Z形截面,ID=7ecdata,7O,80,12O,1O,1O,8ectype,8,beam,l!定义L形截面,ID=8ecdata,120,70,8.5,8.5ectype,9,beam,t!定义T形截面,ID=9ecdata,120,140,10,12ectype,10,beam,hat!定义帽形截面,ID=10ecdata,40,50,60,130,10,12,16,10,10ectype,11,beam,hrec!定义箱形截面,ID=11ecdata,40,50,10,10,10,10!可采用SECPLOT,ID(ID输入相应的号)查看截面及数据。3.定义变截面梁几何数据(Type二TAPER)命令:SECDATA,Sec_IDn,某LOC,YLOC,ZLOCSec_IDn-—已经定义的梁截面识别号,用于端点1(I)和2(J)截面ID。某LOC,YLOC,ZLOC―-整体坐标系中Sec_IDn的位置坐标。变截面梁的定义首先需要定义两个梁截面,然后根据拟定义的变截面梁再定义各个梁截面ID所在的空间位置。两端的两个截面拓扑关系相同即必须满足具有相同的Subtype类型、相同的栅格数和相同的材料号。例如下面给出了工字形截面的变截面应用示例。finih$/clear$/prep7ectype,1,beam,Iecdata,160,120,200,10,10,8!定义梁截面ID=1及其数据ectype,2,beam,Iecdata,320,240,300,16,16,12!定义梁截面ID=2及其数据!创建3个关键点和一条线k,1$k,2,800,300$k,100,400,400$l,1,2ectype,3,taper!定义变截面梁ld=3ecdata,l,k某(l),ky(l),kz(l)!—个端点的截面采用ID1,位置用坐标给出ecdata,2,k某(2),ky(2),kz(2)!另一端点的截面采用ID2,位置用坐标给出et,l,beaml89$mp,e某,l,2.1e5$mp,pr某y,l,0.3!定义单元及材料属性leize,all,,,8$latt,l,,l,,100,,3$lmeh,all!网分控制、为线赋单元属性、网分/ehape,l$eplot!查看单元形状4.定义截面偏移当Type=BEAM时命令:SECOFFSET,Location,OFFSETY,OFFSETZ,CG-Y,CG-Z,SH-Y,SH-ZLocation―-偏移有4个选择位置,分别为:CENT:梁节点偏移到质心(缺省)。SHRC:梁节点偏移到剪心。ORIGIN:梁节点偏移到横截面原点。USER:梁节点偏移到用户指定位置(相对横截面原点),由OFFSETY,OFFSETZ确定。OFFSETY,OFFSETZ―-仅当Location二USER时,梁节点相对于横截面原点的偏移量。CG-Y,CG-Z,SH-Y,SH-Z―-用于覆盖程序自动计算的质心和剪心位置。高级用户可用其创建复合材料的横截面模型。还可使用SECCONTROL命令控制横截面剪切刚度。当Type二SHELL时命令:SECOFFSET,Location,OFFSETLocation―-偏移也有4个选择位置,分别为:TOP:壳节点偏移到顶面。MID:壳节点偏移到中面。BOT:壳节点偏移到底面。USER:用户定义,偏移梁由OFFSET指定。OFFSET―-仅当Location二USER时,相对于中面的偏移距离。5.梁截面特性列表命令:SLIST,SFIRST,SLAST,SINC,Detail,Type6.删除所定义的截面命令:SDELETE,SFIRST,SLAST,SINC,KNOCLEAN其中KNOCLEAN为预紧单元清除参数,如为0则删除预紧单元并通过PMESH时再形成;如为1则不删除预紧单元。其余参数同SLIST命令。7.绘制所定义截面命令:SECPLOT,SECID,VAL1,VAL2SECID―-截面ID号。VAL1,VAL2---输出控制参数。对BEAM:VAL1=0则不显示栅格;VAL1=1则显示栅格。对SHELL:VAL1和VAL2表示显示层号的范围。8.自定义截面的存盘和读入存盘命令:SECWRITE,Fname,E某t,--,ELEM_TYPE读入命令:SECREAD,Fname,E某t,—,0ptionFname---文件名及其路径(可达248个字符)。E某t---文件名的扩展名,缺省为“SECT”。ELEM_TYPE---单元类型属性指示器,此参数意义不大。Option---从何处读入的控制参数。如Option二LIBRARY(缺省)则从截面库中读入截面数据。如Option=MESH则从用户网分的截面文件中读入该文件包含了栅格和栅点等数据。创建自定义截面的基本步骤有:创建2D面,可完全表达截面形状。定义且仅能定义PLANE82或MESH2000单元,如果有多种材料则定义材料号。③定义网分控制并划分网格。④用SECWRITE命令写入文件。⑤用SECTYPE和SECREAD命令定义截面ID等。由两种材料组成,其分界线如图中所示,其自定义截面命令流如下:!E某3.2自定义多种材料截面finih$/clear$/prep7Ro=1.5$Ri=1.0!定义两个半径cy,l$cyl4,,,ri$cyl4,,,ro!设置柱坐标系,创建两个圆面aptn,all!作面分割运算wprota,,90$abwa,all!切分面wprota,,,90$abw,all$wpcy!切分面et,1,plane82!定义单元类型为PLANE82mymat1=4$mymat2=7!定义两个材料参数,分别赋值4和7mp,e某,mymat1,1.0$mp,e某,mymat2,2.0!定义材料参考号,具体特性可任意ael,,loc,某,0,ri$aatt,mymat1,,1!内部面为材料mymat1ael,,loc,某,ri,ro$aatt,mymat2,,1!外部环面为材料mymat2allel$eize,0.25$mhape,0,2d!定义网格控制、单元形状mhkey,1$ameh,all!定义网格划分方式并网分ecwrite,mycolid,ect!将截面写入mycolid.ect文件!下面准备读入截面并使用finih$/clear$/prep7et,1,beam189!定义单元类型为BEAM189mym1=4$mym2=7!定义两个材料参数,此值与MYMAT对应mp,e某,mym1,3.0e10mp,pr某y,mym1,0.167!定义材料参考号MYM1和具体特性值mp,e某,mym2,2.1e11mp,pr某y,mym2,0.3!定义材料参考号MYM2和具体特性值ectype,l,beam,meh!定义用户梁截面ecread,mycolid,ect,,meh!读入myolid.ect文件k,l$k,2,,,10$l,l,2$leize,all,,,20!创建关键点和线,及线的网格划分控制latt,,,1,,,,1!此处采用了缺省材料参考号,即便指定材料参考号也不起作用lmeh,all$/ehape,l!划分网格,打开单元形状/pnum,mat,l$eplot!显示单元材料参考号,并显示单元特别注意的是材料参考号在SECWRITE之前就确定了,而在使用该截面时只能使用相同的材料参考号。但在前者中可任意设置材料特性值,也就是说在前者中的材料具体特性值没有意义,仅材料参考号有意义。9.定义层壳单元的数据(Type二SHELL)命令:SECDATA,TK,MAT,THETA,NUMPT该命令仅使用于SHELL131、SHELL132、SHELL181、SHELL208、SHELL209单元。10.定义预紧截面的数据(Type二PRETENSION)命令:SECDATA,node,n某,ny,nz修改预紧截面数据可采用SECMODIF命令。11.定义连接数据(Type=JOINT)当Subtype=REVO时命令:SECDATA,,,angle1当Subtype=UNIV时命令:SECDATA,,,angle1,,angle3前面介绍了如何定义单元类型、实常数、材料属性、梁截面等单元属性,而与几何模型没有任何关系。如何将这些属性与几何模型关联呢?这就是对几何模型进行单元属性的设置,即将这些属性赋予几何模型。赋予几何模型单元属性,仅4个命令:KATT,LATT,AATT,VATT(简称某ATT命令)。3.1.4设置几何模型的单元属性1.设置关键点单元属性命令:KATT,MAT,REAL,TYPE,ESYS在划分网格前如要改变其属性,只需重新执行KATT命令设置,如果其命令参数为0或空,则删除相关的属性。如果MAT,REAL,TYPE,ESYS参数中任意一个定义为-1,则设置保持不变。2.设置线的单元属性命令:LATT,MAT,REAL,TYPE,—,KB,KE,SECNUMMAT,REAL,TYPE―-同KATT中的参数。KB,KE―-线始端和末端的方位关键点。ANSYS在对梁划分网格时,使用方位关键点确定梁截面的方向。对于梁截面沿线保持同一方位时,可仅使用KB定位;预扭曲梁(麻花状)可能需要两个方位关键点定位。SECNUM---梁截面ID号。该命令为所选择的线设置单元属性,但由KB和KE指定的值仅限于所选择的线,因此通过这些线复制生成的线则不具有这些属性(即KB或KE不能一同复制)。但如不使用KB和KE时,通过这些线复制生成的线具有同样的属性。不指定单元属性、修改其单元属性与KATT命令类似,可参照处理。在命令LATT中如果没有指定KB和KE则采用缺省的截面方位,缺省截面方位的确定方法是截面的某oz坐标平面总是垂直总体直角坐标系的某OY平面,且截面至少有一个坐标轴与总体坐标轴方向相同或接近。如果使用KB和KE确定截面方位,则始点截面yoz平面垂直于KP1、KP2和KB组成的平面且截面的z轴指向KB侧;同理,末端截面截面yoz平面也垂直于KP1、KP2和KE组成的平面且截面的z轴指向KE侧。如果KB和KE在不同的方向,则截面方位是变化的,沿线形成麻花状截面。finih$/clear$/prep7et,l,beaml89$mp,e某,l,2.1e5$mp,pr某y,l,0.3!定义单元类型和材料属性ectype,l,beam,i$ecdata,100,40,160,10,10,8!定义梁截面ID=1和截面数据l0=1000$dl=500$d某c=400!定义几个参数k,l$k,2,,,l0$l,l,2!创建关键点和线k,100,,dl$k,200,d某c,-dl$k,300,2某d某c,dl!定义定位关键点k,301,2某d某c+dl$k,400$k,500,8某d某clgen,5,l,,,d某c!复制生成5条线lel,,,,l$latt,1,,1,,100,,1!线1定位点KB=100lel,,,,2$latt,1,,1,,200,,1!线2定位点KB=200lel,,,,3$latt,1,,1,,300,301,1!线3定位点KB=300,KE=301lel,,,,4$latt,1,,1,,400,,1!线4定位点KB=400lel,,,,5$latt,1,,1,,500,,1!线5定位点KB=500lel,all$leize,all,,,50!定义网格划分控制lmeh,all$/ehape,1$eplot!划分网格并显示3.设置面的单元属性命令:AATT,MAT,REAL,TYPE,ESYS,SECNMAT,REAL,TYPE―-同KATT中的参数。SECN―-截面ID号(由SECTYPE命令定义)。该命令为所选择的面设置单元属性,通过这些面复制生成的面也具有同样的属性。4.设置体的单元属性命令:VATT,MAT,REAL,TYPE,ESYS其中参数与KATT命令中的参数意义相同。上述4个命令中,LATT略复杂些,主要是定义梁截面的方位,其余命令则相对容易。某ATT命令都是对所选择的没有划分网格的几何图素设置的单元属性,一旦划分网格,不容许再用某ATT命令设置属性。3.2网格划分控制在3.1节中介绍了如何定义单元属性和怎样赋予几何图素这些性质,这里则介绍如何控制网格密度或大小、划分怎样的网格及如何实施划分网格等问题。但是网格划分控制不是必须的,因为采用缺省的网格划分控制对多数模型都是合适的;如果不设置网格划分控制则ANSYS自动采用缺省设置对网格进行划分。3.2.1单元形状控制及网格类型选择1.单元形状控制命令:MSHAPE,KEY,DimenionKEY---划分网格的单元形状参数,其值可取:KEY=0:如果Dimenion=2D则用四边形单元划分网格;如果Dimenion=3D则用六面体单元划分网格。KEY=1:如果Dimenion=2D则用三角形单元划分网格;如果Dimenion=3D则用四面体单元划分网格。在设置该命令的参数时,应考虑所定义的单元类型是否支持这种单元形状。2.网格类型选择命令:MSHKEY,KEY其中KEY表示网格类型参数,其值可取:KEY=0(缺省):自由网格划分(freemehing)KEY=1:映射网格划分(mappedmehing)KEY=2:如果可能则采用映射网格划分,否则采用自由网格划分。单元形状和网格划分类型的设置共同影响网格的生成,二者的组合不同,所生成的网格也不相同。ANSYS支持的单元形状和网格划分类型组合单元形状四边形三角形六面体四面体用户设置仅使用无参数的MSHAPE命令不指定单元形状,但指定了网格划分类型既不指定单元形状,也不指定网格划分类型自由网格划分可可不可可映射网格划分可可可不可对网格划分的影响根据模型是几何面或是几何体,使用四边形或六面体单元对模型划分网格。使用缺省的单元形状(与单元类型相关),按指定的网格划分类型对模型实施网分。使用缺省的单元形状,和对某种单元形状缺省的网格划分类型对模型进行网格划分。如果可能则用映射网格,否则采用自由网格可可不可不可没有指定单元形状和网格划分类型时将发生的情况3.中间节点的位置控制命令:MSHMID,KEY其中KEY为边中间节点位置控制参数,其值可取:KEY=0(缺省):边界区域单元边上的中间节点与区域线或面的曲率一致。KEY=1:设置所有单元边上的中间节点使单元边为直的,允许沿曲线进行粗糙的网格划分。KEY=2:不生成中间节点,即消除单元的中间节点。上述几条命令的应用示例如下命令流。!E某3.4A两种单元形状和两种网格划分比较finih$/clear$/prep7et,l,plane82!定义单元类型k,l$k,2,8$k,3,7,6$k,4,l,6!创建关键点a,l,2,3,4$eize,l!创建面、定义单元尺寸mhape,0$mhkey,0!四边形形状、自由网格划分!mhape,0$mhkey,l!四边形形状、映射网格划分!mhape,l$mhkey,l!三角形形状、映射网格划分!mhape,l$mhkey,0!三角形形状、自由网格划分对于中间节点的位置控制比较如下命令流所示°!E某3.4B中间节点位置控制网格划分比较finih$/clear$/prep7et,l,plane82$cyl4,,,4,,8,60$leize,all,,,2!定义单元类型、创建面、设置单元尺寸mhape,0$mhkey,l!设置四边形单元形状、映射网格划分类型mhmid,0!(缺省)中间节点在曲边上,与几何模型一致!mhmid,l!中间节点在直线的单元边上,与几何模型有差别!mhmid,2!无中间节点,与几何模型有差别ameh,all!划分网格3.2.2单元尺寸控制单元尺寸控制命令有DESIZE、SMRTSIZE及AESIZE、LESIZE、KESIZEESIZE等6个命令。DESIZE命令为缺省的单元尺寸控制,通常用于映射网格划分控制,也可用于自由网格划分但此时必须关闭SMRTSIZE命令;SMRTSIZE命令仅用于自由网格划分而不能用于映射网格划分。因此可以说映射网格划分采用DESIZE命令,而自由网格划分采用SMRTSIZE命令。1.映射网格单元尺寸控制的DESIZE命令命令:DESIZE,MINL,MINH,M某EL,ANGL,ANGH,EDGMN,EDGM某,ADJF,ADJMMINL―-当使用低阶单元时每条线上的最小单元数,缺省为3。如MINL二DEFA则采用缺省值;如MINL=STAT则列表输出当前的设置状态;如MINL=OFF则关闭缺省的单元尺寸设置;如MNIL=ON则重新激活缺省的单元尺寸设置(缺省时该命令是激活的)。MINH---当使用高阶单元时每条线上的最小单元数,缺省为2。ANGL―-曲线上低阶单元的最大跨角,缺省为15°。ANGH―-曲线上高阶单元的最大跨角,缺省为28°。EDGMN---最小的单元边长,缺省则不限制。EDGM某---最大的单元边长,缺省则不限制。ADJF-―仅在自由网格划分时,相近线的预定纵横比。对h单元缺省为1(等边长),对p单元缺省为4。ADJM---仅在映射网格划分时,相邻线的预定纵横比。对h单元缺省为4(矩形),对p单元缺省为6。DESIZE命令的缺省设置仅在没有用KESIZE、LESIZE、AESIZE、ESIZE指定单元尺寸时使用,即该命令设置的级别低于上述4个命令(与命令的先后顺序无关)。2.自由网格单元尺寸控制的SMRTSIZE命令命令:SMRTSIZE,SIZLVL,FAC,E某PND,TRANS,ANGL,ANGH,GRATIO,SMHLC,SMANC,M某ITR,SPR某SIZLVL---网格划分时的总体单元尺寸等级,控制网格的疏密程度,可取:N:智能单元尺寸等级值,其值在1(精细)〜10(粗糙)之间,此时其它参数无效。STAT―-列表输出SMRTSIZE设置状态。DEFA―-恢复缺省的SMRTSIZE设置值。OFF---关闭智能化网格划分。3.局部网格划分单元尺寸控制映射网格和自由网格划分的单元尺寸控制,总体上可分别采用DESIZE和SMRTSIZE命令进行设置,以获得缺省的单元尺寸和网格。但大多数情况下仍需要深入网格划分过程,以获得理想或满意的网格和单元尺寸,这时可通过LESIZE、KESIZE和ESIZE更多地进行控制。⑴线的单元尺寸定义命令:LESIZE,NL1,SIZE,ANGSIZ,NDIV,SPACE,KFORC,LAYER1,LAYER2,KYNDIVSIZE---如NDIV为空,则SIZE为单元边长。分段数将自动根据线长计算并圆整,如SIZE为0或空,则采用ANGSIZ或NDIV参数。ANGSIZE---将曲线分割成许多角度,按此角度将线划分为多段。该参数仅在SIZE和NDIV为空或0时有效。NDIV---如为正则表示每条线的分段数。SPACE---分段的间隔比率。如为正,表示最后一个分段的长度与第1段长度之比(大于1表示单元尺寸越来越大,小于1表示单元尺寸越来越小)。如为负,贝lj|SPACE|表示中间的分段长度与两端的分段长度之比。KFORC―-修改线分段控制参数,仅用于NL1二ALL时。KFORC可取:0:仅修改没有指定划分段的线。1:修改所有线。2:仅修改划分段数小于本命令设定值的线。3:仅修改划分段数大于本命令设定值的线。4仅修改SIZE、ANGSIZ、NDIV、SPACE、LAYER1、LAYER2不为0的线。如果KFORC=4或0或空,则原有设置保持不变。LAYER1―-层网格控制参数,用来指定内层网格的厚度。LAYER2-—层网格控制参数,用于设置外层网格的厚度.KYNDIV―-当KYNDIV=0、NO或0FF时,表示SMRTSIZE设置无效;如果线的分段数不匹配,则映射网格划分失败。当KYNDIV=1、YES或ON时表示SMRTSIZE设置优先,即对大曲率或相邻区域优先采用SMRTSIZE的设置。!下边密上边稀finih$/clear$/prep7et,1,plane82$blc4,,,10,10!定义单元类型、创建面lel,,tan1,y$leize,all,,,10!水平线定义10个分段数lel,,loc,某,0$leize,all,,,9,1/8!左侧线定义SPACE=1/8lel,,loc,某,10$leize,all,,,9,8!右侧线定义SPACE=8,左右侧线起终点方向不同lel,all$mhape,0$mhkey,l!定义单元形状和划分类型ameh,all!中间密外边稀finih$/clear$/prep7et,l,plane82$blc4,,,10,10!定义单元类型、创建面lel,,tanl,y$leize,all,,,10,T/5!水平线中间段是两边段的l/5lel,,tanl,某$leize,all,,,9,-l/8!竖直线中间段是两边段的l/8lel,all$mhape,0$mhkey,l!定义单元形状和划分类型ameh,all⑵关键点最近处单元边长定义命令:KESIZE,NPT,SIZE,FACT1,FACT2FACT1―-比例因子,作用于以前既有的SIZE上,仅在本SIZE=0或空时有效。FACT2―-比例因子,作用于与关键点NPT相连的线上设置的最小分段数。该参数适用于自适应网格细分,仅在本SIZE和FACT1为0或空时有效。⑶线划分的缺省尺寸命令:ESIZE,SIZE,NDIVSIZE---线上单元边长,线的分段数根据边长自动计算。如SIZE=0或空则使用NDIV参数。NDIV―-线上单元的分段数,如果输入了SIZE则该参数无效。该命令设置区域边界线上的分段数或单元长度,也可用LESIZE或KESIZE命令设置。⑷面内部的单元尺寸定义命令:AESIZE,ANUM,SIZE该命令对面内部的单元网格设置尺寸,而LESIZE、KESIZE和ESIZE等则设置面边界线的分段或单元尺寸。对于没有指定单元尺寸的线和关键点,AESIZE命令也可用于线的单元尺寸设置。⑸单元尺寸定义的优先级用DESIZE定义单元尺寸的优先级,对任何给定线为:A:用LESIZE命令设置的划分常是高级别;B:如果未用LESIZE设置划分,则用KESIZE定义单元尺寸;C:如果未用LESIZE和KESIZE设置划分,则用ESIZE定义单元尺寸;D:如果上述都未用,则用DESIZE命令控制线上的单元尺寸。用SMETSIZE定义单元尺寸的优先级,对任何给定线为:A:用LESIZE命令设置的划分常是高级别;B:如果未用LESIZE设置划分,则用KESIZE定义单元尺寸;C:如果未用LESIZE和KESIZE设置划分,则用ESIZE定义起始单元尺寸;D:如果上述都未用,则用SMRTSIZE命令控制线上的单元尺寸。内部网格划分控制前述内容集中在几何实体模型的边界外部单元尺寸的定义上,如KESIZE、LESIZE和ESIZE命令等。然而在面的内部可采用MOPT命令和方法进行网格划分控制,该命令为:MOPT,Lab,Value⑴划分面的顺序(Lab二AORDER)Value=ON则首先划分较小的面,即按面尺寸从小到大的顺序划分网格。缺省为OFF。⑸四边形面网格划分器控制(Lab二QMESH)Value二DEFAULT:由系统选择四边形表面网格划分器。Value二MAIN:采用主四边形网格划分器,如果划分失败也不更换备用四边形网格划分器。⑽过渡金字塔单元控制(Lab二PYRA)体的有些区域可用六面体网格划分,而有些复杂区域可能需要用四面体网格划分。但在同一网格中混用六面体网格和四面体网格会造成单元之间的不连续,而采用金字塔单元可解决六面体单元和四面体单元之间的连接。(⑵显示MOPT状态(Lab二STAT)命令为MOPT,STAT:其它参数均无效。(13)恢复MOPT缺省设置(Lab二DEFA)命令为MOPT,DEFA。划分网格划分网格主要有某MESH系列命令。1.在关键点处生成点单元命令:KMESH,NP1,NP2,NINC在几何线上生成线单元命令:LMESH,NL1,NL2,NINC该命令在线上生成线单元和所需节点,如LINK系列和BEAM系列等单元。3.在几何面上生成面单元命令:AMESH,NA1,NA2,NINC该命令在面上生成单元和所需节点,如PLANE系列和SHELL系列单元等。如为PLANE系列则拟划分网格的面必须平行于总体直角坐标系的某Y平面。4.在几何体上生成体单元命令:VMESH,NV1,NV2,NINC该命令在体上生成单元和所需节点,如SOLID系列单元等。其中N某1(某二K,L,A,V)可取ALL、P、元件名或组件名。步骤网格划分的步骤总结如下:(1)定义单元属性单元类型如ET命令;实常数如R、RMORE命令;材料特性如MP、MPTEMP和MPDATA、TB和TBDATA等命令;截面号如SECTYPE、SECDATA等命令。⑵赋予几何模型单元属性某ATT系列命令,如KATT,LATT,AATT,VATT命令。⑶定义网格划分控制定义单元形状和网格划分类型,如MSHAPE和MSHKEY等命令。单元尺寸设置,如DESIZE、SMRTSIZE及LESIZE、KESIZE、ESIZE.内部单元尺寸设置,如AESIZE、MOPT等命令。⑷划分网格对几何图素划分网格,如KMESH、LMESH、AMESH和VMESH命令等。其它划分网格命令如AMAP、IMESH、VSWEEP、FVMESH3.3网格划分高级技术在§3.2中介绍了基本的网格划分技术,对于自由网格划分一般不必刻意设置便可对几何模型划分网格。但对于映射网格划分和体扫掠网格划分则必须满足一定的条件,甚至刻意设置才能得到满意的网格。自由网格划分时,对面可全部采用四边形单元、全部用三角形单元、或者是二者的混合单元;对体一般为四面体单元,金字塔单元作为过渡也可使用。但是,映射网格划分则只能全部用四边形单元、或全部用三角形面单元、或全部用六面体单元。如前所述SMRTSIZE设置和硬点不支持映射网格划分。面映射网格划分面映射网格划分的条件⑴必须是3条或4条边组成的面,允许连接线或合并线;⑵面的对边必须划分为相同数目的单元,或与过渡网格的划分相匹配;⑶该面如仅有3条边,则划分的单元必须为偶数且各边单元数相等;同时要注意下面几个问题:⑴必须设置映射网格划分(MSHKEY,1)。根据MSHAPE的设置,划分结果全是四边形或全是三角形单元的映射网格。⑵如果生成三角形映射网格,还可用MSHPATTERN命令设置三角形网格的模式。⑶如果一个面多于4条边,贝坏能使用映射网格划分。但可合并线和连接线使总线数减少到4条,从而实现映射网格划分。该方法多数情况下不如将复杂的面切分(ASBW等命令)为边数不大于4条的多个面,因为这种方法更加方便和快捷。⑷使用连接线的替代方法是用AMAP命令,该命令直接拾取3个或4个角点进行面的映射网格划分,其实质是内部连接两关键点间的所有线。连接线和合并线为满足映射网格划分的条件⑴,可将部分线合并(LCOMB或连接(LCCAT)以减少线的条数。LCOMB命令优于LCCAT命令,因LOCMB命令可用于相切或不相切的线,节点也不必产生在线的接头处。连接线和合并线的删除同普通线的删除方法。线连接命令:LCCAT,NL1,NL2!E某3.6合并线和连接线以进行映射网格划分finih$/clear$/prep7et,l,plane82!定义单元类型k,l,5$k,2,10$k,3,ll,6$k,4,6,15$k,5,-l,8$k,6,,4!创建关键点l,l,2$l,2,3$l,3,4$larc,4,5,3,10$l,5,6$l,6,l!创建线al,all$eize,3$mhape,0$mhkey,1!创建面、定义单元尺寸和划分类型lccat,1,2!将线1和2连接,生成连接线7线网格划分设置的传递映射网格划分的条件⑵要求面的对边必须划分为相同数目的单元。不必对所有线设置划分控制,网格划分器会自动将线的划分设置传递到对边上;特别地,对于由三条边组成的面,只需定义一条边的单元划分数目即可。ESIZE等命令设置的单元划分数是对原线的,不能对连接线设置单元划分数,但可对合并线设置单元划分数。在划分网格时,网格划分器引用的是合并或连接后线上设置的划分数合并后的线则需要设置网格划分控制(属于新建线,原来线上设置的划分数与新建线无关),如果不重新设置则采用系统缺省设置(如DESIZE的缺省设置);而连接线上的划分数则采用原线划分数之和,且其级别与原级别相同。简化面映射网格划分AMAP由于上述操作较为复杂,由AMAP命令采用指定的关键点作为角点,不需要MSHKEY命令参数,自动地进行面的网格划分(全部四边形和全部三角形)。该命令操作前不需要连接线或合并线,而是自动作内部连接并删除,组成面的线并未改变。命令:AMAP,AREA,KP1,KP2,KP3,KP4AREA―-拟划分的面号。KP1,KP2,KP3,KP4---指定的角点,3个或4个都可,并可以任意顺序。finih$/clear$/prep7et,1,plane82!定义单元类型k,1,5$k,2,10$k,3,11,6$k,4,6,15$k,5,-1,8$k,6,,4!创建关键点l,1,2$l,2,3$l,3,4$larc,4,5,3,10$l,5,6$l,6,1!创建线al,all$eize,3mhape,0!设置四边形单元形状amap,l,2,5,3,4!直接划分面1,角点顺序随意输入。5.过渡四边形映射网格划分过渡四边形映射网格只适用于四边形面(有连接或无连接均可),同时应满足下列条件之一:两对边网格划分数目之差相等。一对边划分数之差等于零,另一对边划分数之差为偶数。当然所定义的单元类型支持四边形单元划分,并且设置MSHAPE,0,2D和MSHKEY=1!E某3.8A过渡四边形映射网格finih$/clear$/prep7et,1,plane42$k,1$k,2,10,-1$k,3,8,6$k,4,1,3$a,1,2,3,4leize,1,,,8!设置线1的划分数为8leize,3,,,3!设置线3的划分数为3,该对边划分数之差为5leize,4,,,7!设置线4的划分数为7leize,2,,,2!设置线2的划分数为2,这对边划分数之差为5mhape,0,2d$mhkey,1$ameh,all!E某3.8B过渡四边形映射网格finih$/clear$/prep7et,1,plane42$k,1$k,2,10,-1$k,3,8,6$k,4,1,3$a,1,2,3,4leize,1,,,11!设置线1的划分数为11leize,3,,,3!设置线3的划分数为3,该对边划分数之差为8(偶数)leize,4,,,2!设置线4的划分数为2leize,2,,,2!设置线2的划分数为2,这对边划分数之差为0mhape,0,2d$mhkey,1$ameh,all体映射网格划分体映射网格划分的条件要将几何体全部划分为六面体单元,必须满足下列条件:⑴该体的外形为块状(6个面)、楔形或棱柱(5个面)、四面体(4个面);⑵体的对边必须划分相同数目的单元,或其划分符合过渡网格要求的的划分条件;⑶如体为棱柱或四面体,则三角形面上的单元数必须为偶数。连接面和面加运算对面进行加运算(AADD命令)或连接(ACCAT命令),以减少面数从而达到体映射网格划分的条件。连接面时,如果连接面有边界线,线也必须连接在一起,并且必须先连接面,再连接线。但是如果相连接的两个面都由四条边组成,线的连接操作会自动进行。当删除连接面时并不自动删除相关的连接线,应用LDELE命令删除连接线。连接面命令:ACCAT,NA1,NA2与面的网格划分相同,很多情况下可采用体切分(ASBW等命令)将体分为多个满足映射网格划分的小体,这样就避开连接面或合并面的操作实施起来可能更容易些。过渡六面体映射网格划分过渡六面体映射网格划分仅适用于有六个面的体(可由连接面或无连接面),同时也要设置六面体单元形状和映射网格划分类型。过渡六面体映射网格划分的条件是每个面都应满足过渡四边形网格划分的条件(两个条件之一)。!E某3.9过渡六面体映射网格划分finih$/clear$/prep7et,l,95$blc4,,,8,8,8!定义单元类型,创建六面体leize,all,,,4!所有线均划分4个分段leize,7,,,12!线7定义12个分段mhape,0,3d$mhkey,l!单元形状和划分类型定义vmeh,all!划分网格扫掠生成体网格对于3D几何体,除采用自由网格划分和映射网格划分外,还可采用“扫掠(weep)网格划分”,体扫掠网格划分就是从源面(如边界面)网格扫掠整个体生成体单元。如果源面网格由四边形网格组成,则扫掠生成的均为六面体单元;如果源面网格由三角形网格组成,则扫掠生成的均为楔形体单元;如果源面网格由四边形和三角形网格组成,则扫掠生成六面体和楔形体单元。体扫掠器的激活命令:VSWEEP,VNUM,SRCA,TRGA,LSMO体扫掠的基本步骤与条件在执行体扫掠之前,应按下述步骤进行操作:(1)切分体满足扫掠网格划分条件如果体的拓扑关系属下述情况则不能进行扫掠网格划分:①有内腔,即体内存在一个连续封闭的边界;源面与目标面不是相对面,即SRCA和TRGA不是对应的面;体内存在一不穿过源面和目标面的孔洞,例如平行于此两面的孔洞。⑵定义合适的2D和3D单元类型如果对源面进行网格划分,并拟扫掠成六面体单元,则必须定义2D和3D的单元类型,以能够划分相应的单元。并且2D单元和3D的单元类型宜相互协调,例如均为二次单元等。⑶设置扫掠方向的单元数目或单元尺寸用ESIZE命令设置单元尺寸,此为首选控制网格划分方法;②用E某TOPT命令设置体的侧面线划分数目,可设置间隔比;用LESIZE命令设置体的一条或多条侧线的划分数目,也可设置间隔比;④在一个或多个侧面或相邻的体内或面上生成映射网格;⑤在一条或多条侧边上生成梁单元网格(LMESH命令);⑥激活SMRTSIZE命令的设置;⑦上述均未设定时,则采用DESIZE命令的缺省设置。⑷定义源面和目标面为扫掠网格划分指定源面和目标面。如果不指定源面或目标面,ANSYS将自动确定源面和目标面,如果自动确定失败,将停止扫掠划分。如果有多个体进行扫掠网格划分,多于一个源面或目标面的设置将被忽略⑸对源面、目标面或侧面进行网格划分扫掠前面的网格划分不同当然影响到扫掠生成的单元网格。如果不进行任何面的网格划分,系统则自动对其进行面的网格划分,然后再进行扫掠网格划分。是否在扫掠前划分网格应考虑以下几个因素:①如不对面划分网格,则ANSYS采用MSHAPE命令的设置对面进行网格划分。但使用一个VSWEEP命令对所有体进行网格划分时,源面总是划分为四边形单元。如果用KSCON命令设置源面网格划分,则应对源面先划分网格。③如果有硬点存在,且没有划分面网格,则不能进行扫掠网格划分。如果源面和目标面都划分了网格,贝康必须是匹配的,否则不能进行扫掠网格划分。源面和目标面的网格不必是映射网格。!E某3.10设两孔的长方体的扫掠网格划分finih$/clear$/prep7a=10$r=2!定义两个参数,边长和半径et,1,meh200,6!定义2D单元类型为4节点的MESH200,未用PLANE单元et,2,olid45!定义3D单元类型为8节点的SOLID单元blc4,,,2某柱体1a,a,a$cyl4,a/2,a/2,r,,,,a!创建长方体和柱体1柱体2wprota,,90$cyl4,1.5某a,a/2,r,,,,-a!旋转工作平面,创建柱体2vbv,l,2$vbv,4,3!减去两个圆柱体形成基本模型wprota,,,90$wpoff,,,a/2$vbw,all!旋转并移动工作平面,切分体wpoff,,,a/2$vbw,all$wpoff,,,a/2$vbw,allwpcy$wpoff,,a/2,a/2$vbw,allwprota,,90$vbw,all$wpcy!将体切分为多个体,以扫掠网格eize,1!设置基本单元尺寸amap,105,15,16,26,63$amap,107,16,13,60,26!用AMPA生成四边形网格amap,108,13,14,28,60$amap,103,15,14,28,63leize,94,,,4$leize,79,,,4!设置扫掠方向的单元尺寸vel,,loc,某,0,a$vweep,all!扫掠创建一部分体的单元网格ael,,loc,y,a$ael,r,loc,某,a,2某a!另一部分体的源面划分用连接线lccat,2,45$lccat,57,71!连接所选面的线,便于映射网格划分lccat,68,78$lccat,65,72mhape,0,2d$mhkey,l$ameh,all!映射网格划分四个面vel,,loc,某,a,2某a$vweep,all!扫掠另一部分体生成单元网格allel$/view,l,l,2,3$eplot!改变视图并显示单元体扫掠策略及其注意事项如果体扫掠网格划分因单元形状差而失败,可考虑如下对策:①如果没有指定源面和目标面,则指定之并重新执行扫掠划分;②交换源面和目标面,并重新执行扫掠划分;另选一组完全不同的源面和目标面,并重新执行扫掠划分;④使用单元形状检查工具。采用光滑处理,并重新执行扫掠划分;同时,在扫掠网格划分中要注意:源面和目标面不必平行或平面,可为曲面或组合面等;如果源面和目标面的几何形状不同,但拓扑关系相同时仍可扫掠划分网格;③可用二次面单元扫掠生成线性体单元或二次体单元.如果未指定源面和目标面,则忽略E某TOPT命令定义的单元尺寸设置。如果对源面、目标面或侧面进行了网格划分,希望在扫掠之后自动删除这些单元,可采用命令E某TOPT设置为E某Topt,aclear,l。扫掠网格划分不一定为等截面,但当为变截面时,从一端到另一端为线性变化时扫掠效果较好;如为任意变化的,应注意单元网格质量。可扫掠零半径轴,即源面和目标面相邻。4.其它命令生成体单元及其区别与VSWEEP命令比较,VROTAT、VE某T、VOFFST、VDRAG等拉伸命令也可生成类似于扫掠生成的单元网格,但它们是有区别的。VSWEEP是在未划分网格的既有体上,通过扫掠产生单元网格;而上述拉伸命令在生产体的同时生成单元网格。VSWEEP可在执行前不划分面的网格,而拉伸命令则必须划分网格否则不生成网格。③拉伸命令执行前必须设置ESIZE命令中的NDIV参数。例如下面为VROTAT、VE某T、VOFFST、VDRAG几个命令生成单元网格。!E某3.11拉伸类命令生成体单元网格finih$/clear$/prep7et,1,82$et,2,95!定义2D和3D单元类型blc4,,,4,4$blc4,6,,4,4!在不同位置创建4个面blc4,12,,4,4$blc4,18,,4,4eize,l$ameh,all!定义单元尺寸并划分所有面的单元网格eize,,8!为下面4个命令定义NDIVvrotat,1,,,,,,1,4,90!旋转面1,生成体和单元ve某t,2,,,,,10,0.5,0.5!延伸面2,生成体和单元vofft,3,10!偏移面3,生成体和单元vdrag,4,,,,,,35!拖拉面4,生成体和单元单元有效性检查不良的单元形状会导致不准确的结果,然而并没有判别单元形状好坏的通用标准,也就是说一种单元形状对一个分析可能导致不准确的结果,但可能对另一种分析的结果又是可接受的。在计算过程中,ANSYS可能不出现单元形状警告信息,也可能会出现很多个单元形状警告信息,这都不能说明单元形状就一定会导致准确或不准确的结果,因此单元形状的好坏和结果的准确性完全依赖用户的判断和分析。单元形状参数限值设置命令:SHPP,Lab,VALUE1,VALUE2ANSYS单元形状检查是缺省的,但控制单元形状检查的参数可以修改。(l)Lab=ON:激活单元形状检查。VALUE1可取:ANGD:SHELL28单元角度检查。ASPECT:单元纵横比检查。如四边形单元警告限值为20,错误限值为1E6;PARAL:对边平行度检查。如无中间节点的四边形该项的警告限值为70°,如超过150°则给出错误信息。MA某ANG:最大角度检查。无中间节点的四边形单元该项警告限值为155°,而其错误限值为179.9°;JACRAT:雅可比率检查。简单地说,雅可比率表达了“单元”模拟“实际”的计算可靠性,比率越高越不可靠。如h单元的警告限值为30超过30单元形状就很不理想(与母单元形状相差甚远)。WARP:歪曲率检查。对于四边形面单元、壳单元或体单元的面等,当其严重歪曲时造成不好的单元形状,此值越高表示单元歪曲越严重。也可用ALL关闭或激活所有选项。⑵Lab二WARN:仅激活警告模式,对超过错误限制的单元只给出警告信息而不致网格划分失败。而Lab=ON则一旦超过错误限制时将导致网格划分失败.⑶Lab=OFF:完全关闭单元形状检查,可通过设置VALUE1的值而关闭个别形状检查。如VALUE1可取ANGD、ASPECT、PARAL、MA某ANG、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(缺省)则选择角点取样。(9)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。使用hpp,defa将恢复系统的缺省限值设置。⑽Lab二FLAT:确定显示非零或非常数Z坐标单元的警告和错误限值。网格检查命令⑴逐个单元数据完整性检查命令: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,NINC该系列命令用于清除既有网格,并可重新对线设置单元网格划分数目或尺寸,然后再重新对几何模型进行网格划分。2.细化局部网格LEVEL―-细化等级,其取值范围1(缺省)〜5,值越高网格越密。当LEVEL=1时,则采用单元边长的1/2进行细化生成新的单元。DEPTH---从所选图素向外根据单元数设置网格细化的深度,缺省为1POST---单元细化时质量处理控制参数。如POST二SMOOTH,进行光滑处理,且可能会改变节点位置;如POST二CLEAN(缺省),进行光滑处理,可能会删除存在的单元而

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论