ANSYSMesh技巧.doc

网格划分的控制主要考虑以下三个因素：（1）单元形状（elementshape）（2）中节点的设置(midsidenodeplacement)（3）单元尺寸（elementsize） 现在分别加以说明： 单元形状：对于2d的面的划分，可以采用三角形单元或者四边形单元。对于3D的体的划分，要么采用六面体单元，要么采用四面体单元。二者的混合使用一般不推荐使用。若采用(transitionalpyrmidelement)过渡的金字塔单元，可以采用二者的混合使用。 中间节点设置的控制（controlingplacementofmidsidenodes）ANSYS默认情况下，将具有中节点的单元的中节点设置在边界线上或边界的面上。 单元尺寸的设置 （1）对于采用free方式的smartelementsizing(smrtsize)方法：该方法具有如下优点：首先计算面或体中线的单元边的尺寸；其次，若采用四边形单元，所有边的划分为偶数。smrtsize控制方法：basic，简单的设置划分等级(level),1(finemesh)10(coarsemesh)。Advacedcontrol,可以控制划分的质量，使网格尽可能的满足要求。 （2）对于采用mapped方式的默认的单元尺寸（defaultelementsize）通过命令：desize来修改采用mapped方式在每一条线上划分份数。对于大型模型，首先查看模型的划分是很有必要的：例如Et,1,45Mshape,0Mshkey,1Lesize,allLplot改变单元尺寸：Desize.Lesize,all,1Lplot （3）局部网格控制I.esize(整体尺寸控制)，可以采用面或体中最短线之间的距离（一般划分23个单元），来控制整体单元尺寸II.kesize（指定点控制）III.lesize（指定线控制）。上述方法可以联合使用，若指定相互冲突，ANSYS根据以下优先级来确定划分的数量：Lesizekesizeesize （4）内部网格的控制（interiormeshcontrol）以上所述均为通过边界来控制单元的尺寸,也可以通过内部的面或体来控制单元的尺寸,通过mopt来实现。I.mopt,expand：作用从边界的细分到内部的粗分。II.mopt,trans,value:实现从边界上小的单元到内部的单元尺寸的过渡（与(1)联合应用）III.mopt,aorder,on：在多个体划分网格时，为了保证小的面上的单元数和相对较高的质量，采用此命令，保证首先在较小的面上划分网格。问题1一空心半球体，其内半径为1500mm，外半径为2000mm。现要沿其径向粘贴CFRP(Ex=235000Mpa, Ey=2000Mpa, u=0.3), 沿环向粘贴GFRP(Ex=150000Mpa, Ey=1000Mpa, u=0.3)。球体采用solid45单元，CFRP和GFRP均采用shell 41单元，CFRP与GFRP与球体完全粘结(即节点重合)，且CFRP与GFRP间隔粘贴，其水平投影如图所示。如何建立其有限元模型，请各位附上各自的命令流。要求：尽量用mapped mesh. 沿环向分成36个单元，沿厚度方向分成2个单元，沿径向分成5个单元。并求其承载能力在粘贴CFRP和GFRP前后改变多少？条件没给全：1、cfrp和gfrp的各向异性的材性没给全，只好随便猜了。2、cfrp和gfrp的厚度没给，只好随便取了0.1mm。3、半球体的材性没给，随便取了一个值。4、既然是弹性分析，何来“承载能力”之说？模型建好了，命令流如下，还没仔细检查，请指正。fini/clear,start/prep7/uis,msgpop,3 !不显示任何警告!定义局部坐标系LOCAL,12,2,0,2,0, ,90, ,1,1,!定义单元、材性et,1,41et,2,45r,1,0.0001 !假设FRP厚0.1mmmp,ey,1,235e9 !CFRP的材性(注意此处为ey)mp,ex,1,2e9mp,ez,1,2e9mp,prxy,1,0.03mp,pryz,1,0.3mp,prxz,1,0.03MP,GXY,1,100E9MP,GYZ,1,100E9MP,GXZ,1,100E9!1.0-prxy*2*ey/ex-pryz*2*ez/ey-prxz*2*ez/ex-2.0*prxy*pryz*prxz*ez/ex>0mp,ex,2,150e9 !GFRP的材性mp,ey,2,1e9mp,ez,2,1e9mp,prxy,2,0.3mp,pryz,2,0.3mp,prxz,2,0.3MP,GXY,2,50E9MP,GYZ,2,50E9MP,GXZ,2,50E9mp,ex,3,200e9 !假设半球的材性mp,prxy,3,0.3!建半球剖面并分网*do,i,1,10CYL4,0,0,1.5,(i-1)*9,1.75,i*9CYL4,0,0,1.75,(i-1)*9,2,i*9*enddoaglue,allLESIZE,all, ,1, , , , ,0amesh,all!旋转成半球体EXTOPT,ESIZE,1,0EXTOPT,ACLEAR,1TYPE,2MAT,3VROTAT,all, , , , , ,75, 78,360,36CM,qiu,volu !定义一个叫“球”的组!建GFRPcsys,0ASEL,s,LOC,y,sin(18/180*3.14159)*2,sin(36/180*3.14159)*2ASEL,a,LOC,y,sin(54/180*3.14159)*2,sin(72/180*3.14159)*2csys,2ASEL,r,LOC,x,2,3CM,gfrp,area !定义一个叫“GFRP”的组csys,0TYPE,1MAT,2real,1amesh,all!建CFRPallsel,allCSYS,12,ASEL,s,LOC,y,0,10*do,i,3,36,2ASEL,a,LOC,y,(i-1)*10,i*10*enddocsys,2ASEL,r,LOC,x,2,3csys,0CM,temp,AREAAGEN,2,all, , ,0,0,0, ,1CMSE,U,tempCM,cfrp,area !定义一个叫“CFRP”的组TYPE,1MAT,1real,1LESIZE,all, ,1, , , , ,0amesh,all!耦合allsel,allcpintf,all/solu!加约束ASEL,s,LOC,y,0,0da,all,uy!加荷载csys,2ASEL,s,LOC,x,0,1.5csys,0SFA,all,1,PRES,1000 !假设荷载为1000牛/m2/PSF,PRES,NORM,2,0,1allsel,all!solve问题二设计思路：其实ansys的复制生成功能很强,如果能找到规律可以使问题大大简化 版主能不能鼓励一下,给个分!别让我的积分总是零!哈哈 /PREP7 B1=5 !X方向 B2=5 !Y方向 B3=3 !圆柱在z方向的定位 B4=20!长方体的总长度 CL1=8!大圆柱的高度 CL2=5!小圆柱的高度 R1=2 !小圆柱的半径 R2=3 !大圆柱的半径 ET,1,SOLID45 K,1,0,0,0 K,2,B1,0,0 K,3,B1,B2,0 K,4,0,B2,0 K,14,0,0,-B3 K,15,0,0,-(B3+CL1) K,16,0,0,-(B3+CL1+CL2) K,17,0,0,-B4 L,1,14 L,14,15 L,15,16 L,16,17 CSYS,1 K,5,R1,0,0 K,6,R1,90,0 K,7,R1,45,0 K,8,R2,0,0 K,9,R2,90,0 K,10,R2,45,0 K,11,R1*0.5,0,0 K,12,R1*0.5,90,0 K,13,R1*0.5,45,0 L,6,7 L,7,5 L,9,10 L,10,8 CSYS,0 A,1,12,13,11 A,11,5,7,13 A,13,7,6,12 A,5,8,10,7 A,7,10,9,6 A,8,2,3,10 A,10,3,4,9 AGEN,2,6,7,1,0,0,-B3 AGEN,2,4,7,1,0,0,-(B3+CL1) AGEN,2,1,7,1,0,0,-(B3+CL1+CL2) !生成第一段体长度B3 VDRAG,1,2,3,4,5,6,1 VDRAG,7,1 !生成大圆柱所在的体CL1 VDRAG,8,9,2 !生成小圆柱所在的体CL3 VDRAG,10,11,12,13,3 !生成另外异端的体 VDRAG,14,15,16,17,18,19,4 VDRAG,20,4 CSYS,1 VGEN,2,1,20,1,0,90,0 ESIZE,4 VMESH,ALL问题五 ，圆球顶柱身 /PREP7 K,1,1,0,0, K,2,0,1,0, K,3,0,0.9,0, K,4,.9,0,0,K,5,0,0,0, K,6,1,-1.5,0, K,7,0.9,-1.5,0,LARC,1,2,5,1, ! line 1LARC,4,3,5,.9, ! line 2 L, 4, 7 ! line 3 L, 7, 6 ! line 4 L, 6, 1 ! line 5 L, 2, 3 ! line 6 L, 4, 1 ! line 7 AL,1,6,2,7 AL,3,4,5,7 APLOT FLST,2,2,5,ORDE,2 FITEM,2,1 FITEM,2,-2 FLST,8,2,3 FITEM,8,3 FITEM,8,5 VROTAT,P51X, , , , , ,P51X, ,360, ,ET,1,MESH200 !* ET,2,SOLID95 !* KEYOPT,1,1,7 KEYOPT,1,2,0 ASEL,S,LOC,Y,0 FLST,5,4,5,ORDE,4 FITEM,5,5 FITEM,5,13 FITEM,5,21 FITEM,5,29 CM,_Y,AREA ASEL, , , ,P51XCM,_Y1,AREACHKMSH,AREA CMSEL,S,_Y !* AMESH,_Y1 !* CMDELE,_Y CMDELE,_Y1 CMDELE,_Y2 !* ALLSEL,ALL FLST,5,8,6,ORDE,2 FITEM,5,1 FITEM,5,-8 CM,_Y,VOLU VSEL, , , ,P51XCM,_Y1,VOLUCHKMSH,VOLU CMSEL,S,_Y !* VSWEEP,_Y1 !* CMDELE,_Y CMDELE,_Y1 CMDELE,_Y2问题六finish/clear/filname,rt !* r=37.3/prep7et,1,63r,1,0.1mp,ex ,1,10e6mp,nuxy,1,0.3mp,dens,1,0.1/386.1pcirc,0,rcsys,1k,5,28.3,11.25local,11,0,27.75622,5.521056,0,11.25csys,11k,6,2.25k,7,-2.25k,8,7.08545,3.45k,9,7.08545,-3.45k,10,8.5,1k,11,8.5,-1k,12,20,1k,13,20,-1al,1,2,3,4a,5,6,8,10,12,13,11,9,7aovlap,1,2adele,3,1csys,1agen,16,4,22.5aovlap,all/pnum,area,1/replotcsys,1k,200,24.2857,88.5k,201,24.2857,49.5k,202,24.2857,10.5k,203,27.15,88.5k,204,27.15,49.5k,205,27.15,10.5csys,0k,206,24.4610,1.5k,207,27.3835,1.5l,200,201l,201,202l,202,206l,206,207l,207,205l,205,204l,204,203l,203,200lfillt,4,20,0.5lfillt,51,20,0.5al,69,87,2,3,60,4,51,96,105,20arsym,x,5arsym,y,5,7,2aovlap,27,5,7,8,9adele,5adele,7,9,1adele,1,4adele,6adele,10adele,14adele,18,26nummrg,kp/replot!* !* numcmp,all !* /PNUM,KP,1 !* ALLSEL,BELOW,AREA ASEL,INVE LSEL,INVE KSEL,INVE !* LDELE,all, ,1 KDELE,all numcmp,all !* ALLSEL,ALL !* !* wpro,90.000000, ASBW, 1 numcmp,all !* wpro,90.000000 ASBW, 1 !* FLST,2,2,5,ORDE,2 FITEM,2,2 FITEM,2,4 ADELE,P51X, , ,1 ADELE, 3 !* numcmp,all !* KWPAVE, 46LSBW, 44 LSBW, 48 numcmp,all !* KWPAVE, 33FLST,2,4,4,ORDE,4 FITEM,2,1 FITEM,2,40 FITEM,2,44 FITEM,2,50 LSBW,P51X numcmp,all !* KWPAVE, 36wpro,-25.00000 !* FLST,2,2,4,ORDE,2 FITEM,2,44 FITEM,2,53 LSBW,P51X numcmp,all !* KWPAVE, 25FLST,2,3,4,ORDE,3 FITEM,2,50 FITEM,2,54 FITEM,2,56 LSBW,P51X !* KWPAVE, 27wpro,-20.000000 LSBW, 42 numcmp,all !* KWPAVE, 16FLST,2,3,4,ORDE,3 FITEM,2,1 FITEM,2,38 FITEM,2,49 LSBW,P51X numcmp,all !* FLST,2,2,4,ORDE,2 FITEM,2,31 FITEM,2,-32LCOMB,P51X, ,0 FLST,2,2,4,ORDE,2 FITEM,2,23 FITEM,2,-24LCOMB,P51X, ,0 FLST,2,2,4,ORDE,2 FITEM,2,11 FITEM,2,-12LCOMB,P51X, ,0 FLST,2,2,4,ORDE,2 FITEM,2,47 FITEM,2,52 LCOMB,P51X, ,0 FLST,2,2,4,ORDE,2 FITEM,2,46 FITEM,2,56 LCOMB,P51X, ,0 numcmp,all !* KWPAVE, 17wpro,-25.000000 FLST,2,3,4,ORDE,3 FITEM,2,37 FITEM,2,56 FITEM,2,58 LSBW,P51X numcmp,all !* KWPAVE, 4FLST,2,3,4,ORDE,2 FITEM,2,58 FITEM,2,-60LSBW,P51X numcmp,all !* WPCSYS,-1,0!* KWPAVE, 40wpro,90.000000, LSBW, 37 numcmp,all !* KWPAVE, 43LSBW, 65 numcmp,all !* KWPAVE, 31!* a,2,67,43,47 a,66,67,43,40 a,50,31,30,28 a,29,35,34,51 a,32,43,52,33 a,33,52,55,24 !* a,51,34,23,53 a,53,23,22,20 a,34,33,24,23 a,24,55,44,25 !* ASBW, 1 numcmp,all !* a,47,69,31,50 a,69,43,32,31 a,40,66,46,48 a,40,48,49,54 a,21,27,26,56 a,56,26,15,57 a,57,15,14,10 a,11,19,18,60 a,26,25,16,15 a,25,44,59,16 a,16,59,62,17 !* a,60,18,6,63 a,63,6,8,12a,18,17,4,6a,17,62,65,4 a,4,65,45,3a,3,45,38,5!* KWPAVE, 7LSBW, 14 numcmp,all !* wpro,-90.000000 LSBW, 38 !* KWPAVE, 39LSBW, 101 numcmp,all !* KWPAVE, 37LSBW, 38 !* wpro,-40.000000 KWPAVE, 36LSBW, 102 numcmp,all !* KWPAVE, 42LSBW, 38 numcmp,all !* KWPAVE, 64LSBW, 103 numcmp,all !* a,64,76,75,42 a,42,75,74,36 a,36,74,73,37 a,37,73,72,39 !* wpro,80.000000 KWPAVE, 38LSBW, 100 numcmp,all !* a,38,39,72,77 KWPAVE, 5LSBW, 110 numcmp,all !* a,5,38,77,78 a,5,78,71,7a,7,71,1,70a,9,7,70,13!* KWPAVE, 41!* wpro,-60.000000 LSBW, 103 numcmp,all !* FLST,2,3,4,ORDE,3 FITEM,2,51 FITEM,2,55 FITEM,2,61 LCOMB,P51X, ,0 numcmp,all !* a,54,49,77,41 a,41,77,74,62 !* !* AATT, 1, 1, 1, 0 !* !* /auto,1 aplot !* !* LESIZE,ALL, , ,2, ,1, , ,1,!* FLST,5,44,4,ORDE,36FITEM,5,3 FITEM,5,7 FITEM,5,9 FITEM,5,11 FITEM,5,15 FITEM,5,20 FITEM,5,-22FITEM,5,27 FITEM,5,-29FITEM,5,39 FITEM,5,-40FITEM,5,42 FITEM,5,46 FITEM,5,-49FITEM,5,51 FITEM,5,-53FITEM,5,55 FITEM,5,59 FITEM,5,-60FITEM,5,70 FITEM,5,-71FITEM,5,77 FITEM,5,-78FITEM,5,83 FITEM,5,86 FITEM,5,89 FITEM,5,91 FITEM,5,95 FITEM,5,102FITEM,5,-106 FITEM,5,108FITEM,5,-109 FITEM,5,111FITEM,5,-112 FITEM,5,114FITEM,5,-115 CM,_Y,LINE LSEL, , , ,P51XCM,_Y1,LINECMSEL,_Y !* LESIZE,_Y1, , ,6, , , , ,1 !* FLST,5,4,4,ORDE,4 FITEM,5,61 FITEM,5,80 FITEM,5,-81FITEM,5,116CM,_Y,LINE LSEL, , , ,P51XCM,_Y1,LINECMSEL,_Y !* LESIZE,_Y1, , ,8, , , , ,1 !* !* MSHKEY,1 AMESH,all MSHKEY,0 !* !* ARSYM,X,all, , ,0,0 ARSYM,Y,all, , ,0,0 !* nummrg,node,0.01,0.01,low nummrg,kp,0.01,0.01,low !* !* /auto,1 eplot !* save问题七fini/clear,start/filname,davidswpi/Tital, 3D model!* /prep7et,1,142et,2,92k,1k,2,0.025k,3,0.025,-0.084k,4,-0.084k,5,0.065k,6,0.065,-0.084k,7,0.050,-0.0031k,8,0.150,-0.0031k,9,0.150,-0.042k,10,0.050,-0.042a,2,5,6,3a,7,8,9,10vrotat,1,1,4vrotat,2,9,10vsba,7,4vsba,8,4vsba,5,19vsba,6,19vdel,10,13,1,1vglud,all!*KWPAVE, 9wpro,90.000000,FLST,2,4,6,ORDE,4FITEM,2,2FITEM,2,-3FITEM,2,6FITEM,2,10VSBW,P51X!*numcmp,allFLST,2,10,6,ORDE,6FITEM,2,1FITEM,2,-4FITEM,2,6FITEM,2,-8FITEM,2,10FITEM,2,-12VDELE,P51X, , ,1numcmp,all!*VSEL,S, , , 2ALLSEL,BELOW,VOLU!*KWPAVE, 11!*wpro,90.000000CYL4, , ,0.015, , , ,-0.1 ! 小圆柱半径 0.015ASBA, 8, 15VDELE, 3, , ,1 numcmp,all !* l,16,12l,17,13!* AL,22,13,25,26 VSBA, 2, 13KL,7,0.5, KWPAVE, 18 wprota,-45, VSBW, 4 alls其实，这个体的划分并不是那么简单。大家可以看一下davidswpi的原图，其上面的实心圆柱体的直径比下面空心圆柱体的小园直径大，也就是说如果在空心圆柱体上切出一个与实心圆柱体直径一样大的实体的话，这个体就比较不规则了，要想用规整网格划分以及保证节点重合的话，就不那么容易了，有空偶也试试。nickel1 说得不错，前面几位的模形与 davidswpi 有一定差别，即突出的小圆柱直径小于空心圆柱的孔径，这样比较容易划分。 我做了一个结果，是将 davidswpi 的模型进行分割后得到的。结果见后面的图，命令流如下： fini /clear,start /filname,davidswpi /Tital,3D model !* /prep7et,1,142et,2,92k,1k,2,0.025k,3,0.025,-0.084k,4,-0.084k,5,0.065k,6,0.065,-0.084k,7,0.050,-0.0031k,8,0.150,-0.0031k,9,0.150,-0.042k,10,0.050,-0.042a,2,5,6,3a,7,8,9,10vrotat,1,1,4vrotat,2,9,10vsba,7,4vsba,8,4vsba,5,19vsba,6,19vdel,10,13,1,1vglud,all !* !*以下新加内容 !* KWPAVE, 9wpro,90.000000, FLST,2,4,6,ORDE,4 FITEM,2,2 FITEM,2,-3 FITEM,2,6 FITEM,2,10 VSBW,P51X !* numcmp,all !* !* FLST,2,10,6,ORDE,6 FITEM,2,1 FITEM,2,-4 FITEM,2,6 FITEM,2,-8 FITEM,2,10 FITEM,2,-12VDELE,P51X, , ,1 numcmp,all !* VSEL,S, , , 2ALLSEL,BELOW,VOLU !* KWPAVE, 11!* wpro,90.000000 CYL4, , ,0.02, , , ,-0.1 !小圆柱半径 0.02 可以改变 !* ASBA, 8, 15 VDELE, 3, , ,1 numcmp,all !* l,16,12 l,17,13 !* AL,22,13,25,26VSBA, 2, 13 !* VSEL,S, , , 4ALLSEL,BELOW,VOLU !* KWPLAN,-1, 7, 17, 8 numcmp,all !* FLST,2,3,3,ORDE,3 FITEM,2,7 FITEM,2,-8 FITEM,2,17 KWPAVE,P51XASBW, 10 numcmp,all !* KWPLAN,-1, 8, 7, 13 ASBW, 8 numcmp,all !* AL,27,26,28,6VSBA, 3, 18 !* allsel,all nummrg,kp,1e-4,1e-4,low numcmp,all !* k,18, k,19,1!* FLST,2,1,5,ORDE,1 FITEM,2,6 FLST,8,2,3 FITEM,8,18 FITEM,8,19 VROTAT,P51X, , , , , ,P51X, ,90,1, !* !* LESIZE,ALL,0.02, , , ,1, , ,1, !* MSHAPE,0,3dMSHKEY,1 VMESH,all MSHKEY,0 !* VSYMM,Z,all, , , ,0,0 !* WPCSYS,-1,0KWPAVE, 10CSWPLA,11,0,1,1, !* CSYS,11, VSYMM,Y,all, , , ,0,0 CSYS,0, !* allsel,all nummrg,node,1e-5,1e-5,low nummrg,kp,1e-4,1e-4,low numcmp,all !* WPSTYLE,0 !* /auto eplot !* save问题：建模命令： /prep7 BLC4,0,0,10,10,10 SPH4,0,0,5 VSBV, 1, 2 分块思路i此模型不难！ 分块如下：然后直接vsweep就ok了！方法一：/prep7!正方体一角减去圆球的映射网格，关键在拓扑上形成5面体 blc4,1,1,1 /VIEW,1,1,1 /rep sph4,0.5 vsbv,1,2/pnum,area,1 /rep asel,2,6,2 aplot accat,all !共顶点的三个面连接使在拓扑上形成5面体 alls /pnum,line,1 /rep lccat,8,12!使面只有4条边，才能映射分网 lccat,5,10!使面只有4条边，才能映射分网 lccat,2,3 !使面只有4条边，才能映射分网 mshape,0,3dmshk,1 et,1,45 vmesh,all /view,-1,-1,-1 /rep方法二：*afun,deg /prep7 et,1,mesh200 et,2,solid95 keyopt,1,1,7 keyopt,1,2,0 block,0,1,0,1,0,1 sph4,0,0,0.3 vsbv,1,2,sepo,delete,delete