第1页共17页三维结构的结构化网格划分structure对于三维结构,只有模型区域满足以下条件,才能被划分为结构化网格:①没有孔洞、孤立的面、孤立的边、孤立的点;②面和边上的弧度值应该小于90°;③三维区域内的所有面必须要保证可以运用二维结构化网格划分方法;④保证区域内的每个顶点属于三条边;⑤必须保证至少有四个面(如果包含虚拟拓扑,必须仅包含六条边);⑥各面之间要尽可能地接近90°,如果面之间的角大于150°,就应该对它进行分割;⑦若三维区域不是立方体,每个面只能包含一个小面,若三维区域是立方体,每个面可以包含一些小面,但每个小面仅有四条边,且面被划分为规则的网格形状。三维结构的扫略网格划分Sweep(扫略网格划分技术)对于三维结构,只有模型区域满足以下条件,才能被划分为扫略网格:①连接起始面和目标面的的每个面(称为连接面)只能包含一个小面,且不能含有孤立的边或点;②目标面必须仅包含一个小面,且没有孤立的边或点;③若起始面包含两个及两个以上的小面,则这些小面间的角度应该接近180°;④每个连接面应由四条边组成,边之间的角度应接近90°;⑤每个连接面与起始面、目标面之间的角度应接近90°;⑥如果旋转体区域与旋转轴相交,就不能使用扫略网格划分技术;⑦如果被划分区域的一条或多条边位于旋转轴上,ABAQUS/CAE不能用六面体或楔形单元对该区域进行扫略网格划分,而必须选择Hex-dominated形状的单元;⑧当扫略路径是一条封闭的样条曲线时,该样条曲线必须被分割为两段或更多。第2页共17页总结(1)对于不能采用结构化技术(Structured)和扫略技术(Sweep)进行网格划分的复杂结构,用户可以运用Partition工具将其分割成形状较为简单的区域,并对这些区域进行结构化或扫略网格划分。如果模型不容易分割或分割过程过于繁杂,用户可以选用自由网格划分技术(Free)。本章将通过例7-2详细介绍这两种操作过程。(2)采用映射网格划分(Mappedmeshing)能得到高质量的网格,但ABAQUS/CAE不能直接采用映射网格划分技术,只能通过Usemappedmeshingwhereappropriate选项让程序选择映射网格划分的区域。在以下几种情况下,用户可以选择该项进行映射网格划分:2D+Quad/Quad-dominated+Free+Advancingfront(采用自由网格划分技术和进阶算法,对二维结构划分四边形或四边形占优的单元)、2D+Tri+Free、3D+Hex/Hex-dominated+Sweep+Advancingfront、3D+Tet+Free。(3)中轴算法(Medialaxis)和进阶算法(Advancingfront)是主要的ABAQUS网格划分算法,有四种单元形状(ElementShape)和网格划分技术(Technique)的组合能选用这两种算法:2D+Quad+Free和3D+Hex+Sweep默认选择中轴算法,2D+Quad-dominated+Free和3D+Hex-dominated+Sweep默认选择进阶算法。对于不同的模型,用户应该比较这两种算法,得到合适的网格。下面总结一些获得高质量网格的参数设置。1.尽量采用结构化(Structured)或扫略(Sweep)网格划分技术对三维实体模型划分六面体单元。如果单元扭曲较小,建议选用计算精度和效率都高的非协调模式单元;否则选用二次六面体单元。2.若采用扫略技术划分网格,中轴(Medialaxis)算法和进阶(Advancingfront)算法的选择没有统一的标准,用户需要针对实际模型进行尝试。一般情况下,选择中轴算法包含的Minimizethemeshtransition项或进阶算法包含第3页共17页的Usemappedmeshingwhereappropriate项,可以提高网格质量。3.若复杂模型的分割过程过于耗时,用户可以选用二次四面体单元划分网格。建议读者选择Usemappedtrimeshingonboundingfaceswhereappropriate项,如前所述,ABAQUS/CAE会对形状简单的面选用映射网格划分,通常可以提高网格质量(参见例7-2)。另外,若模型的网格密度足够且重点分析区域位于边界,用户可以选择Increasesizeofinteriorelements项来增加内部单元的尺寸,提高计算效率。4.网格密度是协调计算精度和计算效率的重要参数,但合适的网格密度往往需要根据具体模型而定。一般情况下,用户可以在重点分析区域和应力集中区域加密种子,其他区域可以设置相对较稀疏的种子;如果需要控制一些边界区域的节点位置,可以在设置边种子时进行约束(参见例7-1)。线性静力学分析实例线性静力学问题很容易求解,往往用户更关心的是计算精度和求解效率,希望在获得较高精度的前提下尽量缩短计算时间,特别是大型模型。这主要取决于网格的划分,包括种子的设置、网格控制和单元类型的选取。如第7章中的介绍,应尽量选用精度和效率都较高的二次四边形/六面体单元,在主要的分析部位设置较密的种子;若主要分析部位的网格没有大的扭曲,使用非协调单元(如CPS4I、C3D8I)的性价比很高。对于复杂模型,可以采用分割模型的方法划分二次四边形/六面体单元;有时分割过程过于繁琐,用户可以采用精度较高的二次三角形/四面体单元进行网格划分具有下列特点的问题才适于进行线性瞬态动力学分析:1、系统应该是线性的:线性材料特性,无接触条件,无非线性几何效应。2、响应应该只受较少的频率支配。当响应中各频率成分增加时,例如撞击和冲击问题,振型叠加技术的有效性将大大降低。3、载荷的主要频率应在所提取的频率范围内,以确保对载荷的描述足够精确。4、由于任何突然加载所产生的初始加速度应该能用特征模态精确描述。5、系统的阻尼不能过大。第4页共17页接触1、在接触分析模拟中一般最好在那些将会构成从面的模型部分使用一阶单元,使用二阶单元可能会出现问题,这是由接触算法决定的。其中线性减缩积分单元(C3D8R)和非协调单元(C3D8I)都适合进行接触分析,而修正的二阶四面体单元(C3D10M)是为了应用于复杂的接触模拟问题而设计的,在模型复杂的接触分析中推荐使用,但是计算时间也大大增加。2、网格密度的考虑一般来讲,网格较粗的面作为主面,网格较细的面作为从面。3、稳定的接触关系的建立在做接触分析的相关案例时,一般都先定义一个分析步(或者几个)施加很小的载荷,让各个零部件之间的接触关系平稳地建立起来,之后再建立一个或多个分析步把全部载荷施加到求解模型上,这样虽然增加了更多的分析步,但是减少了收敛的困难,提高了求解效率。特别是在过盈接触分析的案例中,更应该注意多分析步的应用。4、接触主从面的定义在定义接触对中的主从表面时,一般选择刚度较大的面作为主面,这里的刚度不仅指材料的特性,还要考虑分析模型结构的刚度。解析面或者由刚性单元构成的面必须作为主面。主从表面的定义体现在网格密度上,从面的网格应该比主面更加细密,如果主从表面的网格密度大致相等,那么柔性材料的表面应该作为从面。无论是主面还是从面,发生接触的部位应尽量圆滑过渡,不要有尖锐的特征存在。一对接触面的法向方向应该相反,即主面和从面在几何位置上没有发生重叠。这在接触分析中要特别注意,否则可能会导致计算不收敛。第5页共17页分割作用对部件进行分割的目的一般来说有三个方面,第一是将复杂的几何部件分割成相对较简单的区域,以便于高质量的四边形或六面体单元的施加;第二是便于部件局部或者整体区域的网格密度控制;第三是便于对部件分割的不同区域施加不同类型的单元。网格划分的部分可以放在定义部件装配后面完成,这里作者放到定义部件装配前面完成,主要是因为一般来说,如果各个部件先定义装配,之后再定义网格划分,则在后续模型网格划分时为了得到较好的网格质量都会进行模型分割,而模型分割会导致各个模型的一些装配参考发生变化,最终可能导致各个部件的装配关系失败。轴对称特点:1、几何形状轴对称2、边界约束条件对称3、载荷对称简化原则:1、三维对称模型简化为平面此时截面属性section是solid,而不是shell2、二维对称平面模型简化为线模型3、二维轴对称线模型简化为点模型4、模型必须建立在轴对称右侧(Y轴右侧)5、必须在轴对称处的分割边界(如果存在)上施加对称边界条件/反对称边界条件smmetry/antisymmetry/encastrey作用6、单元选择轴对称单元类型7、尽量施加分布载荷,不应施加集中载荷,否则容易引起结果不对称第6页共17页(变为原体:view---ODBDisplayOptions----Sweep/Extrude)1、不同分析步作用取消与作用,:------loadmanager------propagated变为inactive,取消该分析步在另一分析步中的作用。改变继承属性。2、在环境栏中Module后面选择Mesh,进入Mesh模块,在环境栏中Object后面选择Part,可以发现部件的颜色为橙色,说明部件不能使用当前的单元形状(六面体)设置进行网格划分,必须改变单元形状或者对部件进行剖分,使之能使用当前的单元形状进行网格划分。执行Mesh→Controls命令,弹出MeshControls对话框,如图11-14所示,ElementShape栏中选择单元形状为Tet(四面体),其他接受默认设置,单击OK按钮,图形窗口中的模型变为粉色,说明能使用四面体单元对模型进行自由网格划分。3、Merge/cutinstance合并几何体在(part中生成):merge--geometry,options--suppressoriginalinstances,intersectingboundaries---remove合并后可以把原来instance的分开几何体删除以及part中的部件删除。4、动画:animate----timehistory5、参考点的选取:tools---referencepoint----选点--完成6、interaction里面的constraints是工件的约束(工件约束为刚体),assembly里面的constraints是对部件与部件之间位置关系进行一定的约束。部件约束为刚体:Interaction----createconstraints---type:rigidbody-----Continue---regiontype:body---edite---选择实列---referencepoint--edite---选择约束点。7、tools---partition---type:edge、face、cell区别Face----创建分割面Edge---创建分割边Cell----创建分割体8、速度的定义:load---createpredefinedfield---mechanical---velocity(再第7页共17页选择参照点,一般选取该物体质心)10、边界条件的确定:11、刚体只发生平移和转动,不发生变形,可以用一个点的情况控制整个刚体情况,为刚体的参考点,创立刚体部件,需定义刚体参考点,刚体上所以载荷和边界条件都施加在该参考点上。12、初始分析步的边界条件只能为0,载荷也只能在后续分析步建立。13、分析步中存在大位移、大转动、初始应力、几何钢化、突然翻转----打开几何非线性。14、焊接定义;若关心的是整个结构的变形,或者焊缝处不是危险部位,不需要对旱区精确建模。(1)面与面的焊接:绑定约束tie(2)点与面:运动耦合约束kinematiccoupling(3)点与点:可以用运动耦合约束kinematiccoupling,但是把constraintregiontype改为nodeset,不是surface15、stress---sstrain--pe\peeq\pemag\ledisplace--uForce==rf\cfcontact---cstress\cdisp16、三维实体尽量采用6面体