网格划分1 1一、网格种类对应网格,自由网格1.平面2-D结构，对应网格的面积为四边形，相对应边的元素数目一定要相等，网格后元素为四边形，如图(a)及图(c)所示，但也可以强制进行自由网格化(网格化后元素为非四边形)，如图(b)及图(d)。相对应边元素数目相同的对应网格与自由网格示意图相对应边元素数目相同但元素长度不同的对应网格与自由网格示意图2 2相对应边元素数目不同，元素形状不同的自由网格示意图2.自由网格化时则相对应边元素个数不一定相同，如图(e)所示，网格化后元素可为全部三角形，如图(f).或四边形和三角形的混合，如图(e)所示。3.对三角形，进行对应网格化时，则其三边元素数目一定要相等且为偶数，如图(g)、(h)所示。3 34.对3-D结构，对应网格化时，体积必为六面体，相对应的线段元素数目一定要相等如图所示，所形成元素也是六面体。反之，自内网格化时，任何形状的体积皆可，元素形状一定为四面体的三角锥。4 4二、网格化步骤实体模型网格化有限元模型的建立，需要进行3个步骤：1建立、选取元素的数据。2设定网格建立所需的参数。3产生网格。5 5第一个步骤是建立元素的数据这些数据包括元素的种类元素的几何常数元素的材料性质及元素形成时所在的坐标系统，当对象进行网格化后，元素的属性是什么，其命令在前面已有详细介绍。当然，我们可以设定不同种类的元素，相同的元素又可设定不同的几何常数，也可以设定不同的材料特性，以及不同的元素坐标系统。6 6定义X对象网格后元素属性XATT，MAT,REAL,TYPE,ESYS,SECSUM定义X对象网格后元素属性(ATTributes)。MAT,REAL,TYPE,是前面所定义元素材料特性、几何常数、元素类型的号码，与直接建立法相比较，仅将原三行命令合并为一个命令，X代表对象名称。ESYS为建立元素时所在坐标系统号码，通常圆柱形体积可在圆柱坐标系统下建立。系统默认MAT,REAL，TYPE皆为第一组及ESYS为卡式坐标。7 7第二个步骤设定网格建立所需的参数1、元素的数据设定后，即可进行网格划分所需的参数设定.2、最主要在于定义对象边界(即线段)元素的大小与数目。3、网格设定所需要的参数将决定网格的大小、形状，这一个步骤非常重要，将影响分析时的正确性及经济性。4、网格太细也许会得到较好的结果，但并非网格细就是最好的结果，因为网格太密太细，会造成占用大量的分析时间。有时较细的网格与较粗的网格比较起来，较细的网格分析的精确度只增加百分之几，但占用的计算机资源比较粗的网格却是数倍之多，同时较细的网格在复杂的结构中常会造成不同网格划分时连接的因难，这一点不能不特别注意。8 8声明网格时使用何种方法网格化MSHKEY,KEY声明网格时使用自由网格化(KEY=0，系统默认)、对应网格化(KEY=1)或对应网格及自由网格混合(KEY=2)。对应网格化优先考虑，而且SMRTSIZE命令将不起作用。KEY=2时，适用于2-D实体模型，并不适用3-D实体模型，因3-D网格中不允许六面体元素及角推元素共存。如果实体模型不符合对应网格的需求，执行此命令则会产生错误消息。9 9声明网格化时元素的形状MSHAPE，KFY,Dimension声明网格(Meshing)化时元素的形状(SHAPE)。2-D实体模型(Dimension=2D)时采用四边形(KEY=0)或全部为三角形(KEY=1)。虽然对应网格的系统默认元素为四边形和六面体，我们仍可利用该命令，强制元素为三角形和角锥。1010定义所选择的线段进行元素网格化时元素的大小LESIZE，NL1,SIZE，ANGSIZ，NDIV,SPACE,KFORC,LAYER1，LAYER2定义所选择的线段(Line,NLl)进行元素(Element)网格化时元素的大小(SIZE)，元素的大小可用线段的长度(SIZE)或该条线段要分割元素的数目(NDIV)来确定。SPACE为间距比，最后一段长与最先一段长的比值，正值代表以线段的方向为基准，负值代表以中央为基难，系统默认为等间距。NLl=ALL为目前所有的线段。1111定义通过点的所有线段进行元素网格化时元素的大小KESIZE,NPT,SIZE,FACT1,FACT2定义通过点(Keypoint，NPT)的所有线段进行元素(EIement)网格化时元素的大小(SIZE)，不含LESIZE所定义的线段。元素的大小仅能用元素的长度(SIZE)输入。该命令必须成对使用，因为线段基本上含两点。NPT=ALL为通过目前所有点的线段。如仅对一条线段而言，使用该命令两次其结果和直接采用LESIZE一样。1212定义元素网格化时元素的大小ESIZE,SIZE,NDIV该命令以目前所有对象为基准，也就是全部线段(不含LESIZE，KESIZE所定义的线段)，定义元素(Elemet)网格化时元素的大小(SIZE)，元素的大小可以用元素的边长(SIZE)或者线段要分成元素的数目(NDIV)来确定。1313系统默认元素大小DESIZE，MINL，MINH,MXEL,ANGL，ANGH，EDGMIN，EDGMX，ADJF,ADJM该命令为系统默认元素大小，也就是全部线段(不含LESIZE、KESIZE、ESIZE所定义的线段)。系统默认低阶元素(即元素节点与节点间无附加节点，如BEAM3，PLAND2，SOLID45)时，线段最少(MINL)为三段，高级元素(即元素的节点与节点间附加一个节点，如PLANE82，SOLID95)线段最少为两段，线段最多为15段。可修改MINL，MINH,MXEL的值，调整元素的大小，LPLOT无法检查其分割状态。1414自由网格化时网格大小的高级控制SMRTSIZE,SIZVAL，FAC,EXPND，TRANS，ANGL，ANGH，GRATIO，SMHLC，SMANC，MXITR、SPRX该命令用于自由网格化时，在网格大小的高级控制(不合LESIZE、KESIZE、ESIZE所定义的线段)中，系统默认值就不起作用。一般由DESIZE控制元素大小，DESIZE及SMRTSIZE是互相独立的命令，仅能存在一个，执行SMRTSIZE命令后DESIZE自动无效了。该命令可由两个层面输入，第一个是所有参数自行调整，对初学者而言，可暂时不考虑，第二个是利用系统默认的等级输入(SIZVAL)。每一个等级都已设定其后的参数值，共有十级，系统的默认值为第六级，级数越高网格越粗。该命令原理是利用目前实体模型线段长度、曲率自行进行最佳网格化。1515第三步网格的产生当元素属性及网格所需的参数定义后，便可对实体模型进行网格化。线元素由线段网格产生，平面元素由面积网格产生，立体元素由体积网格产生。1616将一组线进行网格化LMESH，NLl，NL2，NINC将一组(NLl，NL2，NINC)线(Lines)进行网格化(MESH)并产生节点及元素，若NLlALL，则NL2及NINC都忽略掉，网格化目前所有线段。1717将一组面积进行网格化AMESH，NAl，NA2，NINC将一组(NAl，NA2，NINC)面积进行网格化(MESH)并产生节点及元素，若NAl=ALL，则NA2及NINC都忽略掉，网格化目前所有面积。1818将一组体积进行网格化VMESH，NVl，NV2，NINC将一组(NVl，NV2，NINC)体积(voLume)进行网格化(MESH)并产生节点及元素，若NVl=ALL，则NV2及NINC都忽略掉，网格化目前所有体积。1919将网格化的一组线段上节点与元素清除LCLEAR，NL1，NL2，NINC将网格化的一组(NLI，NL2，NINC)线段上节点与元素清除(CLEAR)。清除的节点及元素，只与所选的线有关而不管在此线上的节点及元素是否已被筛选出宋，但若一个节点是由很多元素共享，而元素不在所选出来的线上，则不会被删除，节点、元素的号码等都会因节点、元素删除而去掉。2020将网格化的一组面积上节点与元素清除ACLEAR,NA1，NA2，NINC将网格化的一组(NA1，NA2，NINC)面积上节点与元素清除。清除的节点及元素，只与所选的面积有关而不管在此面积上的节点及元素是否已被筛选出来，但若一个节点是由很多元素共享，而元素不在所选出来的面积上，则不会被删除，节点、元素的号码等都会因节点、元素删除而去掉。2121将网格化的一组体积上的节点与元素清除VCLEAR，NVl，NV2，NINC将网格化的一组(NVl，NV2，NINC)体积上的节点与元素清除。清除的节点及元素，只与所选的体积有关而不管在此体积上的节点及元素是否已被筛选出来，但若一个节点是由很多元素共享，而元素不在所选出来的体积上，则不会被删除，节点、元素的号码等都会因节点、元素删除而去掉。2222/PREP7ET,1,PLANE42K,1,-2.5$K,2,6,-2.5$K,3,2.5K,4,6,2.5CSYS,1$K,5,10,-30$K,6,10,30CSYS,0SAVEL,1,2$L,4,3CSYS,1$L,2,4$L,5,6CSYS,0A,1,2,4,3$A,2,5,6,4SAVE网格划分综合练习2323AMESH,ALL此时网格以系统默认DESIZE进行自由网格化，但由于面积符合对应网格的要求，故不需MSHKEY及MSHAPE命令，系统优先采用对应网格四边形元素。元素属性AATT默认为第一组也可省略，以下为改变DESIZE后网格化的效果。2424RESUMEDESIZE,5AMESH,ALLRESUMEDESIZE,10,20AMESH,ALL2525RESUMESMRTSIZE,6AMESH,ALLRESUMESMRTSIZE,8AMESH,ALLRESUMESMRTSIZE,1AMESH,ALL2626RESUME$MSHAPE,1$DESIZE,5$AMESH,ALLRESUME$MSHAPE,1$DESIZE,10,20$AMESH,ALL2727RESUME$MSHAPE,1$SMRTSIZE,6$AMESH,ALLRESUME$MSHAPE,1$SMRTSIZE,8$AMESH,ALLRESUME$MSHAPE,1$SMRTSIZE,1$AMESH,ALL2828RESUME$LESIZE,1,6,6$LESIZE,2,6,3LESIZE,3,5$LESIZE,4,5,2$ESIZE,4SAVEESIZE,4$AMESH,ALLLESIZE，NL1,SIZE，ANGSIZ，NDIV,SPACE,KFORC,LAYER1，LAYER22929ACLEAR,ALL$ESIZE,7$AMESH,ALLRESUME$LESIZE,5,5$AMESH,ALL3030RESUME$ET,1,PLANE2$ET,2,PLANE82ESIZE,4AMESH,1AATT,2AMESH,23131线的划分/UNITS,SI/PREP7ET,1,BEAM3MP,EX,1,207e9R,1,1e-4,2.083e-10,0.005K,1$K,2,0.3L,1,2SaveESIZE,4LATT,1,1,1,0LMESH,1RESUMEESIZE,0.03LATT,1,1,1,0LMESH,13232实体划分/PREP7K,1K,2,1K,3,2,1K,4,-1,1.5K,5,-0.5,2,3K,6,3,1.5,2K,7,4,3,2K,8,-1.5,3,3V,1,2,3,4,5,6,7,8ET,1,SOLID45SAVEVMESH,1RESUMEDESIZE,10,20VMESH,13333RESUMESMRTSIZE,6VMESH,1RESUMEMSHAPE,1DESIZE,5VMESH,1RESUMEMSHAPE,1SMRTSIZE,6VMESH,13434实体模型外力3535实体模型外力1、外力的负载对结构系统而言，不管是直接建立法或间接建立法，力的形式都是相同的。2、直接建立法时，外力是直接作用在节点、平面元素的边或者立体元素的面上。实体模型建立法时，由于节点及元素产生后较难掌握节点、元素号码及其相对应元素边及面的号码，故在元素或节点上很不容易直接施加外力。所以对实体模型而言，施力方式可采用直接建构在实体模型的点、线段、面积上，其相对应效果如同作用于节点、元素边及元素面。当网格化完成后，原先作用在实体模型上的力，可自动转换至相对应的节点及元素上。3、实体模型外力的命令和直接建立法外力的命令非常相似。3636直接法建立外力实体建模后，有些简易的外力，仍然采用直接建立法方式，即完成网格化后，直接施加外力至节点或元素上。集中力施加集中力是作用于一点的力，故较容易掌握，但必须注意施力点要建立一个节点，因为所有与实体模型有关的节点在网格化中，必定会成为一个节点，所以只要知道施力点的号码、网格后再进行检查该点相对应节点号码便可直接施加外力。假设第7点为外力施加点完成网格后检查节点号码如下：KSEL,S,KP,7!选择7点为有效点NSLK!选择有效点上的节点NLIST!显示该节点资料NPLOT!检查该节点3737间接法外力建立间接法的外力直接施加外力于实体模型上，若施加在点、线段及面积上，网格化后都会自动转换到有限元模型上，此转换手续会在/soLu处理器中下达SOLVE命令后自动进行.3838D,NODE,Lab,VALUE,VALUE2,NEND,NINC,Lab2,Lab6NODE:节点号码Lab:相对节点的每一个节点受自由度约束的形式Lab=UX,UY,UZ,ROTX,ROTY,ROTZ(结构力学量)=TEMP(热学量)VALUE:约束值的大小NODE,NEND,NINC:选取节点的范围定义节点自由度约束(有限元模型上加约束)3939与D相对应，在点上定义约束条件DK，KPOI，Lab，VALUE,VALUE2,KEXPND，Lab2，Lab3，Lab4，Lsb5,Lab6该命令与D相对应，定义约束条件(Displaceme)在点上，KPOI为受限制点的号码，VALUE为受约束点的值。Lab-Lab6与D命令相同，可借着KEXPND去扩展定义在不同点间的节点所受的约束.4040在线段上定义约束条件DL，LINE，AREA，Lab，VALUEl，Value2在线段(LINE)上定义约束条件(Displacement)，LINE，AREA为受约束线段及线段所属面积的号码，Lab与D命令相同但增加对称(Lab=SYMM)与反对称(Lab=ASYM)，Value1为受约束的值。4141在面积上定义约束条件DA,AREA，Lab，Value1，Value2在面积上定义约束条件(DispIacement)，AREA为受约束面积的号码，Lab与D命令相同但增加对称(Lab=SYMM)与反对称(Lab=ASYM)，valuel为受约束的值。4242DTRAN将DX(XK，L，A)命令所定义的约束条件转换(TRANsfer)到有限元模块上，当下达SOLVE命令后，该命令为自动执行，使用者也可先转换(在SOLVE指令之前)并检查其结果。4343定义集中力(有限元模型上加约束)F,NODE,Lab,VALUE,VALUE2,NEND,NINCNODE:节点号码Lab:外力的形式Lab=FX,FY,FZ,MX,XY,XZ(结构力学的力,力矩方向),=HEAT(热学的热流量)=AMP,CHRG(电学的电流,电荷)=FLUXVALUE:外力的大小NODE,NEND,NINC:选取节点的范围4444与F相对应，在点上定义集中外力FK,KPOI，Lab，VALUE，VALUE2该命令与F相对应，在点(Keypoint)上定义集中外力(Force)，KPOI为受外力点的号码,VALUE为外力的值。Lab与F命令相向。4545FTRAN将FK命令所定义的外力(Force)转换到有限元模块上，当下达SOLVE命令后,该命令为自动执行，也可先转换(在solve指令之前)并检查其结果。4646定义分布力SFE,ELEM,LKEY,Lab,KVAL,VAL1,VAL2,VAL3,VAL4ELEM:单元号码LKEY:建立单元后,依节点顺序,平面单元四个边,该参数是边号Lab=PRES(表面分布压力)CONV(热力学的对流)VAL1,VAL4:相对应作用于单元边及面上节点的值.4747与SFE相对应，在面积线上定义分布力SFL,LINE，Lab，VALI，VALJ，VAL2I，VAL2J该命令与SFE相对应，在面积线(Line)上定义分布力(surfaceload)作用的方式和大小，应用于2-D的实体模型表面力。LINE为线段的号码，Lab的定义与SFE相同，VALI，VALJ为当初建立线段时点顺序的分布力值。4848与SFE相对应，在体积的面积上定义分布力SFA，AREA，LKEY,Lab，VALUE，VALUE2该命令与SFE相对应，在体积的面积(Area)上定义分布力(surfaceload)作用的方式和大小，应用于3-D的实体模型表面力.AREA为面积的号码，LKEY为当初建立体积时面积的顺序，选择AREA与LKEY其中一个输入。Lab的定义与SFE相同，VALUE为分布力值。4949SFTRAN将SFL及SFA命令所定义的外力转换到有限元模块上，当下达SOLVE命令后，该命令为自动执行，也可先转换并检查其结果。/PSF，Item，Comp，KEY,KSHELL在图形窗口中绘制(plot)分布力(SurFaceload)状态。Item、Comp为外力的类别和方向KEY0为不显示,KEY=1以面型表示,KEYl=2以箭头表示.5050此外还有DKLIST、DKDELE、DLLIST、DLDELE、DALIST、DADELE、FKLIST、FKDELE,SFLLIST、SFLDELE、SFALIST、SFADELE为显示及删除己定义约束条件、已定义外力、己定义分布力，与其他对象的显示和删除相同5151例子：加温度载荷/UNITS,BIN/PREP7ET,1,PLANE55,1MP,KXX,1,50K,1,0.5K,2,1.5K,3,1.5,0.1K,4,0.5,0.1A,1,2,3,4LESIZE,1,5LESIZE,2,1LESIZE,3,5LESIZE,4,1AATT,1,1,1,0AMESH,1FINISH/SOLUDK,1,TEMP,100DK,4,TEMP,100DK,2,TEMP,0DK,3,TEMP,0SOLVEFINISH/POST15252/PREP7ET,1,PLANE42MP,EX,1,30E6K,1$K,2,15K,3,15,2$K,4,15,5K,5,15,8 $K,6,15,10K,7,10$K,8,15,5,1L,1,2$L,2,3$L,5,6L,6,7$L,7,1CIRCLE,4,3,8,5,180NUMMRG,KPAL,ALL$ESIZE,4$AMESH,ALLFINISH/SOLULSEL,S,LINE,5NSLL,S,1SF,ALL,PRES,100,0ALLSELLSEL,S,LINE,6,7NSLL,S,1SF,ALL,PRES,500节点间加均布力ALLSEL/PBC,ALL,1/PSF,PRES,NORM,2NSEL,S,LOC,X,14.99,15.01DSYM,SYMM,XALLSELSOLVEFINISH/POST1NSLL,Type,NKEYSelectsthosenodesassociatedwiththeselectedlines./PBC,Item,-,KEY,MIN,MAX,ABSShowsboundarycondition(BC)symbolsandvaluesondisplays.ALLRepresentsallappropriatelabels.KEY1Showsymbol.5353修改上面讲的例子将此例子在点1,4处加固定约束,点2,3处x方向加集中力100N.线L3上加均布压力10/UNITS,BIN/PREP7ET,1,PLANE42MP,EX,1,2.07E9K,1,0.5$K,2,1.5K,3,1.5,0.1$K,4,0.5,0.1A,1,2,3,4LESIZE,1,5LESIZE,2,1LESIZE,3,5LESIZE,4,1AATT,1,1,1,0AMESH,1FINISH/SOLU.SOLVEFINISH/POST15454Dk,1,all,0Dk,4,all,0Fk,2,fx,100Fk,3,fx,100Sfl,3,pres,105555