八节点六面体有限元网格自动消隐技术
有限元六面体网格的典型生成方法及发展趋势

有限元六面体网格的典型生成方法及发展趋势吕 军,王忠金,王仲仁(哈尔滨工业大学材料科学与工程学院,黑龙江哈尔滨150001)摘 要:工程问题三维有限元仿真的主要困难是模型的建立,而模型的建立需要采用合适的方法来生成高质量的三维有限元网格.以金属塑性成形过程的三维有限元仿真为例,说明了采用六面体单元的必要性.针对典型的有限元六面体网格生成方法,系统地分析了各种方法的实现原理和发展趋势,并探讨了六面体网格生成总的发展趋势.分析结果说明,复杂域内六面体网格全自动生成的实现是全自动网格生成真正走向实用化、通用化必须解决的难题.关键词:数值仿真;有限元法;网格生成;六面体网格中图分类号:TG 302 文献标识码:A 文章编号:036726234(2001)0420485206G eneration of f inite element hexahedral mesh and its trend of developmentL αJ un ,WAN G Zhong 2jin ,WAN G Zhong 2ren(School of Materials Science and Engineering ,Harbin Institute of Technology ,Harbin 150001,China )Abstract :The major difficulty with the 3-D finite element simulation of an engineering problem lies in the construction of models ,which needs the proper generation of 3-D finite element hexahedral mesh of high quality.The necessity to use a hexahedral unit is justified by taking the 3-D finite element simulation of the plastic formation of metals as an example.The theories behind and trends of development of different ways of generating finite element hexahedral meshes are systematically analysed ,and the general trend of development for generation of hexadedral mesh is discussed as well.It is concluded that the full automatic generation of hexahedral mesh in complex domains is the key to the popularization of full automatic genera 2tion of hexahedral mesh.K ey w ords :numerical simulation ;finite element method ;mesh generation ;hexahedral mesh 有限元法是求解工程问题的一种近似数值方法,近年来在工程领域中得到了广泛的应用[1,2].有限元仿真的一个重要步骤是对连续体进行离散化,为使离散出的网格能更精确地逼近连续区域和有限元计算的结果在预定误差范围内,应保证离散化后得到质量较高的网格.为此提出了许多方法来生成有限元网格[35].目前,二维有限元收稿日期:2001-01-20.基金项目:金属精密热加工重点实验室基金资助项目.作者简介:吕 军(1974-),男,博士研究生;王仲仁(1934-),男,教授,博士生导师.网格生成方面已比较成熟,提出了许多行之有效的方法.但在三维有限元网格尤其是六面体网格的生成方面还没有通用的算法,存在许多难点问题需进一步解决.本文论述了采用六面体单元的必要性,并分析了有限元六面体网格的典型生成方法.最后,对六面体网格生成的发展趋势进行了探讨.1 采用六面体单元的必要性在有限元仿真过程中,单元类型的选择对整个有限元仿真的计算效率、自动化程度、计算精度等都将产生重要影响.因此单元类型的选择在各 第33卷 第4期 哈 尔 滨 工 业 大 学 学 报 Vol.33,No.42001年8月 J OURNAL OF HARB IN INSTITU TE OF TECHNOLO GY Aug.,2001个领域的有限元仿真中都占有重要地位.在体积成形刚塑性/刚粘塑性有限元仿真中,单元类型选择的重要性尤为突出.这是因为刚塑性/刚粘塑性有限元仿真有三个突出的特点:(1)塑性成形往往是一个大变形过程,有限元计算中通常需要进行多次网格重划分,而网格重划分需要耗费大量的时间,且每次都会损失一定的精度.(2)塑性变形过程是一个非线性问题,需要进行迭代求解,计算效率问题更为突出.(3)刚塑性/刚粘塑性有限元仿真时必须进行多次工件与模具间的动态接触处理,每一次的处理都会使得工件的有限元模型产生一定的体积损失而影响计算精度.这三种问题的处理都与单元类型的选择密切相关.选择合理的单元类型,就可用较少的网格重划分次数、相同或较少的高斯积分点数来达到较高的计算精度和计算效率,这在有限元仿真中是非常重要的.六面体单元由于变形特性好、计算精度高等优点而在很多三维有限元仿真领域中得到了广泛的应用.在金属体积成形的三维有限元仿真中,要求单元既要有一定的“刚性”(即抗畸变能力)以避免频繁的网格重划分,又要有一定的“柔性”(即良好的变形特性)以准确地仿真变形过程,还必须有较高的计算精度.在体积成形三维有限元仿真中常采用四面体和六面体单元.大量计算结果表明:采用六面体单元进行三维有限元仿真可采用较少的网格重划分次数达到较高的计算精度,故六面体单元是金属体积成形过程三维有限元仿真的首选单元[6,7].2 有限元六面体网格的典型生成方法六面体网格在三维有限元仿真中有四面体网格无法比拟的优越性,但现有的有关三维有限元网格生成方法的文献往往偏重于介绍四面体网格的生成方法,对六面体网格生成方法介绍得很少.实现可靠的、高质量的六面体有限元网格自动生成是三维有限元仿真领域的瓶颈问题.当前,有限元六面体网格的典型生成方法主要有以下几种.2.1 映射单元法映射单元法是三维网格生成中最早使用的方法之一.这种方法先把三维实体交互地分成几个大的20节点六面体区,然后使用形函数映射技术把各个六面体区域映射为很多细小的8节点六面体单元[8].这种方法易于实现,可以生成规整的结构化网格;缺点是当三维实体的表面是十分复杂的自由曲面时,该方法的逼近精度不高,且人工分区十分麻烦、难以实现自动化.近年来,一些研究者采用“整体规划技术(Integer programming technique)”来进行实体的自动分区[9],但该技术很难对复杂形体(如塑性加工中的复杂锻件)进行自动分区.曲面映射是三维映射的特例,采用曲面映射技术可以对几何曲面进行离散化处理[10,11].文献[12]详细研究了基于映射单元法的有限元六面体网格自动生成技术,采用加权因子控制网格生成过程中自然坐标的分割,可以生成密度不同的有限元网格.对原域为单连通凸区域的简单形体及原域为复连通凹区域的复杂形体,该种方法均可生成质量较高的网格,生成的六面体网格如图1所示.图1 映射单元法Fig.1 Mapped element method 映射单元法的发展趋势是:实现简单、规则形状形体的自动分区,提高手工分区的交互性,能方便地进行复杂的三维形体的分区.2.2 基于栅格法这种方法预先产生网格模板,然后将要进行网格化的物体加到其上,并在实体内部尽可能多地填充规则的长方体或正方体网格,在实体的边界上根据实体边界的具体特征更改网格的形状和相互连接关系,使得边界上的六面体单元尽可能地逼近物体的边界形状.文献[13]采用这种方法生成了六面体单元.这种方法能实现网格生成的自动化,网格的生成速度也非常快.其最大弱点是边界单元的质量较差;另一个缺点是所生成的单元尺寸相近,网格密度很难得到控制.1998年发布的MARC/HexMesh模块中采用了基于栅格法,并对这种方法进行了改进,使得初始填充在实体内部的单元尺寸较大、实体边界单元的尺寸较小,这样可以较好地控制网格密度[14].但是,对于复杂三维形体,MARC/ HexMesh模块产生的边界六面体网格的质量仍然不够理想.Tekkaya[15]将改进八叉树法与基于・684・哈 尔 滨 工 业 大 学 学 报 第33卷栅格法相结合来生成边界过渡网格,首先根据工件的边界来区分内部网格和边界网格,然后应用网格细化和均匀化处理来改进边界单元的质量,如图2所示.图2 基于栅格法Fig.2 Grid 2based method 这种方法的发展趋势是:控制实体内部的初始规则网格的尺寸,以控制最终形成的网格的密度;采用网格结构重组(包括拆分和合并单元)和网格优化算法来提高边界单元的质量.2.3 几何变换法这种方法由二维四边形网格经过旋转、扫描、拉伸等几何变换而形成六面体网格,几何变换后删除重节点及四边形、进行单元及节点的重新编号[16].这种方法生成的六面体网格如图3所示;优点是比较容易实现,在当今大多数的大型CAD 软件前置处理中均有此功能.但是,这种方法只适用于形状简单的三维形体,且主要依靠人机交互的方式来实现.图3 几何变换法Fig.3 G eometry transformation method 这种方法的发展趋势是使四边形有限元网格能够以自由曲线为路径进行扫描,尽量减少人机交互的步骤.2.4 改进八叉树法这种方法的基础是三维物体的八叉树表示.所作的改进类似于四叉树法的改进,但三维物体的边界处理更加复杂;它具有改进四叉树法同样的利弊.Y erry 等[17]首先提出并实现了这种方法,他们将物体边界简化为42种可能的模式(18种单平面和24种双平面切割八叉元).这种方法与基于栅格法结合生成“过渡网格”,效果较好[15].著名的有限元分析软件MARC/Auto 2Forge 模块中采用了这种方法,生成的六面体网格如图4所示.图4 改进八叉树法Fig.4 Modified 2octree method 这种方法的发展趋势是:与基于栅格法结合使用来提高过渡网格的质量,并减少仿真过程中的畸形单元,提高形体尖角处单元的质量.2.5 模块拼凑法Yang 等[18,19]把工件分为一定数量的子模块,然后对每一类形状简单的子模块规定一种六面体网格生成方法,整个工件的有限元网格即可由这些子模块内的网格拼凑而成,如图5所示.但是,实际生产中的工件(尤其是模锻件)的形状往往非常复杂,很难对其进行子模块的自动划分,采用专家系统的方法往往也是不可行的.因此,这种方法只能针对形状相对简单和变化较少的工件来生成六面体网格. 这种方法的发展趋势是完善专家系统的知识库,使其能适应更复杂形状工件的子模块自动划分.2.6 单元转换法由于多种四面体网格自动生成算法已经达到实用化的程度,在自动生成四面体网格后,可以把一个直边四节点四面体单元分成四个六面体单元[20,21],这样可以把四面体网格自动地转化为六面体网格,如图6所示.这种方法的缺点是得到的网格是杂乱无章的非结构化六面体网格,网格的质量不高;为了较好地逼近复杂物体的曲面边界,需要生成较多的直边四面体单元,因而也将得到数量极多的六面体单元,这会使得有限元仿真的时间过长.对同时具有内外复杂边界的三维问题(如内部有空洞缺陷的复杂锻件分析),该方法是实现六面体网格自动生成的一种比较有效的方・784・ 第4期 吕 军,等:有限元六面体网格的典型生成方法及发展趋势法.文献[22]对单元转换法进行了改进,将十节点曲边四面体转换为六面体,并采用非线性约束优化算法大幅度提高了六面体网格的单元质量.图5 模块拼凑法Fig.5 Modular method图6 单元转换法Fig.6 Element conversion method 这种方法的发展趋势是减少不必要的四面体单元的数量,采用网格结构重组技术以剔除不必要的单元,采用约束优化算法提高六面体单元的质量.2.7 B 样条曲面拟合插值法这种方法基于三维物体的边界曲面B 样条表示,采用插值拟合曲面来生成六面体网格[23].在几何构形确定的情况下,这种方法即可自动生成六面体网格.通过调整B 样条函数中的参数可以控制网格密度,生成的六面体网格如图7所示.这种方法的优点是边界曲面逼近好,形体的几何表示与网格生成在数学方法上一致;缺点是局部网格的处理比较困难,这是整体域剖分所带来的问题.图7 B 样条曲面拟合插值法Fig.7 B 2spline surface interpolation method 这种方法的发展趋势是采用B 样条曲面和实体造型相结合的方式来描述三维物体,采用模块法来处理物体内部的局部网格.2.8 采用波前法逐层由实体表面向实体内部生成六面体网格(Plastering algorithm)Blacker 和Meyers [24]于1993年提出了这种方法,该方法实际上是二维四边形网格逐层推进生成法[25](Paving algorithm )在三维空间上的拓展.在三维实体内部,各个六面体单元的边与边、面与面之间的相互关系十分复杂,并且只有满足一定条件的实体表面上的节点才能生成完全的六面体网格,故这种方法的实现具有很高的难度.生成的网格如图8所示.该方法生成的六面体网格的单元质量(尤其是边界单元的质量)是所有算法中最好的,但该方法的实现仍需解决一些技术细节上的问题.图8 Plastering 算法Fig.8 Plastering algorithm 这种方法的发展趋势是优化实体表面的布点,避免在向实体内部逐层推进时产生尺寸过小和形状不合理的单元,避免单元间的裂缝.2.9 采用中轴面分解和整体规划技术生成六面体网格这种方法首先将三维实体分解成一定数量的简单子域,然后在每个子域内生成六面体网格[26].在将实体分解成子域过程中采用中轴面(Medial Surface )分解技术,并采用整体规划技术来确定每条边的分割数,进而控制六面体网格的密度[9].根据形体的中轴面可以确定必要的子域,子域可以定义为13种可能类型中的一种[27].中轴面分解方法也可以拓展应用于带有凹边或凹顶点的实体及退化情况,从而可以实现复杂实体(如带有孔、凹角等)的六面体网格生成[28].该方法生成的六面体网格的单元质量很高并且疏密有致,如图9所示. 这种方法的发展趋势是实现复杂形体的全自动中轴面分解;尽可能形成容易网格化的子域;提高边界单元的质量,避免产生形状不好的单元(如・884・哈 尔 滨 工 业 大 学 学 报 第33卷狭长单元).图9 中轴面分解法Fig.9 Medial surface subdivision method3 六面体网格生成的发展趋势有限元六面体网格生成问题近年来成为三维网格生成方法研究的热点和难点,出现了许多种算法,但至今尚未提出一种通用的有限元六面体网格自动生成方法.下列问题将成为六面体网格生成的研究前沿,也将是未来的发展趋势.3.1 开发复杂域六面体网格的全自动生成方法全自动网格生成方法因其高效性、处理复杂情况的能力和便于集成到计算机集成制造系统(CIMS )等优点,已成为网格生成的发展趋势.现有的网格自动生成方法在时效、稳定性和通用性等方面与实用要求都有一些差距,六面体网格生成方面尤为突出,问题的关键在于开发有效、高效的自动生成方法.复杂域的网格生成是全自动网格生成的前提,人们正在研究能在任意复杂域内生成六面体网格的方法,并注重方法的可靠性.复杂域内六面体网格全自动生成的实现是全自动网格生成真正走向实用化、通用化必须解决的难题.3.2 网格密度定义和控制技术的研究三维网格密度定义和控制一直没有行之有效的方法,生成疏密有致的六面体网格并且使密网格和疏网格之间的单元均匀过渡将是六面体网格生成的发展趋势之一.3.3 基于几何造型的六面体网格生成及其集成基于几何造型的网格生成是通往有限元仿真技术集成到计算机集成制造系统的必由之路,目前的网格生成方法大都从造型系统中得到形体描述,但几乎都与造型系统松散结合.实现六面体网格生成与几何造型系统的集成将是未来的发展方向.3.4 六面体网格显示技术及正确性检测六面体网格缺乏有效的显示技术.正因为缺乏直观的显示来验证其正确性,其正确性的检测就显得格外重要;目前在这方面尚缺乏快速有效的方法.六面体网格显示技术的关键是提出正确、可靠和高效的消隐处理算法.4 结 论在三维有限元仿真中采用六面体单元有很多优点,实现可靠、高质量的有限元六面体网格自动生成一直是CAD/CAE 领域内的一个难点,也是制约三维有限元仿真走向实用化的瓶颈问题.解决这个问题的关键是在现有方法的基础上提出稳定、高效和通用的复杂域六面体网格全自动生成方法,以使得六面体网格在三维有限元仿真领域中的应用真正走向实用化.参考文献:[1]GHOUAL I M A ,DUVAU T G.Local analytical de 2sign sensitivity analysis of the forging problem using FEM[J ].Comput Methods Appl Mech Eng ,1998,163:55270.[2]FISH F ,PANDHEERADI M ,BEL SKY V.E fficient solutions schemes for interface problems[J ].Finite El 2ements in Analysis and Design ,1996,22:2672280.[3]施云生,沈国强.基于边界适应的有限元网格自动生成及局部调整技术[J ].锻压技术,1998(4):28230.[4]JOUN M S ,L EE M C.Quadrilateral finite element generation and mesh quality control for metal forming simulation[J ].Int J Num Methods Eng ,1997,40:405924075.[5]JOE B.Tetrahedral mesh generation in polyhedral re 2gions based on convex polyhedron decompositions[J ].Int J Num Methods Eng ,1994,37:6932713.[6]SCHN EIDERS R.A 2grid based algorithm for the gen 2eration of hexahedral element meshes [J ].Eng with Comput ,1996(12):1682177.[7]TEKK A Y A A E ,K AVA K L I S.3-D simulation of metal forming processes with automatic mesh genera 2tion[J ].Steel Res ,1995,66(9):3772383.[8]COO K W A ,OA KES W R.Mapping methods for generating three 2dimensional meshes [J ].Comput in Mech Eng ,1982,8:67272.[9]TAM T ,ARMSTRON G C G.Finite element mesh controlled by integer programing [J ].Int J Num Methods Eng ,1993,36:258122605.[10]VARDHAN I K V R ,PRASAD N S.Mesh genera 2tion for spherical and conical surfaces using transfinite interpolation[J ].Comput &Struct ,1989,32(6):135921362.[11]LAU T S ,LO S H.Finite element mesh generation・984・ 第4期 吕 军,等:有限元六面体网格的典型生成方法及发展趋势over analytical curved surfaces[J].Comput&Struct,1996,59(2):3012309.[12]蒋浩民,刘润广,王忠金,等.基于映射法的三维有限元网格自动划分[J].塑性工程学报,1998,5(3):27231.[13]L EE Y K,Y AN G D Y.A new automatic mesh gen2eration technique and its application to the finite ele2 ment analysis of practical forging process[J].AdvTech Plasticity,1996,1:4092413.[14]News Letter[J].Nippon MARC Analysis Res Corpo2ration Japan,1998,3:9.[15]TEKK A Y A A E.Fully automatic simulation of bulkmetal forming processes[A].Proc NUMIFORM’98[C].Rotterdam:Netherlands,1998.[16]NA GESH K,SRIK AN T A.Automatic mesh genera2tion in2-D and3-D objects[J].Adv Eng S oft2 ware,1989,11(1):19225.[17]YERR Y M A,SHEPHARD M S.Automatic threedimensional mesh generation by the modified2octreetechnique[J].Int J Num Methods Eng,1984,20(11):196521990.[18]Y AN G D Y,Y OON J H,L EE N K.Modularremeshing:a practical method of3-D remeshing inforging of complicated Parts[J].Adv Tech Plasticity,1990,1:1712178.[19]Y OON J H,Y AN G D Y.A three dimensional rigid2plastic finite element analysis of bevel gear forging byusing a remeshing Technique[J].Int J Mech Sci,1990,32(4):2772291.[20]陈 军.虚拟模具制造及金属成形过程三维仿真技术[D].上海:上海交通大学,1996.[21]XIE G,RAMAEKER J A H,Graded mesh genera2tion and transformation[J].Finite Elements in Analy2sis and Design,1994,17:41255.[22]左 旭.集成于CAD系统的汽车零件多工位体积成形三维CAE仿真[D].上海:上海交通大学,1998.[23]王忠金.模锻过程的三维数值模拟及连杆终锻成形规律的研究[D].长春:吉林工业大学,1995.[24]BLACKER T D,MEYERS R J.Seams and wedgesin plastering:a3-D hexahedral mesh generation al2gorithm[J].Eng with Comput,1993,9:83293. [25]BLACKER T D,STEPHENSON M B.Paving:anew approach to automated quadrilateral mesh genera2tion[J].Int J Num Methods Eng,1991,32:8112847.[26]L I T S,MCKEA G R M,ARMSTRON G C G.Hex2ahedral meshing using midpoint subdivision and integerprogramming[J].Comput Methods Appl Mech Eng,1995,124:1772193.[27]PRICE M A,ARMSTRON G C G,SABIN M A.Hexahedral mesh generation by medial surface subdivi2sion:PartⅠ.solids with convex edges[J].Int JNum Methods Eng,1995,38:333523359.[28]PRICE M A,ARMSTRON G C G.Hexahedral meshgeneration by medial surface subdivision:PartⅡsolids with flat and concave edges[J].Int J NumMethods Eng,1997,40:1112136.(责任编辑 王小唯)・94・哈 尔 滨 工 业 大 学 学 报 第33卷。
基于八叉树的六面体网格自动生成算法

者 对其 进 行 了改进 , 广 泛应 用 于 三 维六 面 体 并 网格 的 自动 生成 。笔者 根据六 面体 网格 生成 方 法 的研 究 现 状 , 择 自动 栅 格法 作 为 本文 的基 选
1 实体 模 型 建 立
三 维六 面体 网格生 成方法 的第 一步 就是 要 建 立实 体模 型 。传 统栅 格法 的初始 网格 要足 够
大, 能包 含整 个 实 体模 型 。栅 格 法 是 用 大 小 均 匀 的 网格 然 后拓展 到包 含整个 实体 模型 的大 包 围盒 内 , 这样 做算 法 简 单 , 算 速度 快 , 算 法 计 但 有 致命 缺 陷 , 就 是 网格 大小 反 映不 出原 实体 那 模 型 的几 何 特征 。传统栅 格法 的初 始 网格也要 足 够大 , 能包 含 整 个 实体 模 型 。但 是 与 之 不 同 的是八叉 树算 法是 用一 种大小 不 同的 网格拓 展 到 包含 整 个 实体 模 型 的大 包 围盒 内 , 样 做 网 这 格 大小 能 反 映原 实 体模 型 的几何 特 征 , 其是 尤
文件作 为输入 数据 。如果 实体 模 型是 由 2组 三
角面片组成 的 , 这 2 三角面片是不允 许有 重 那 组 叠和缝 隙情况 出现 的 。图 1所示 是 一个 由 3 S D
文 件 显 示 的 哑 铃 模 型 的三 角 面 片 , 小 8 3 B 大 . 6k ,
收稿 日期 : 0 1 1 一 8 2 1 - l2
本方 法 。
成 为解决 各种 工程 问题 的重要 方法 。有 限元 法
的基 本概念 是用 较简 单 的问题 代替 复杂 问题 后 再求 解 。从 几何 角度 说就 是用 网格 去逼 近实 体
8节点6面体减缩积分单元的积分点个数

8节点6面体减缩积分单元的积分点个数在有限元分析中,积分点的选择对于结果的精度至关重要。
而对于8节点6面体减缩积分单元来说,积分点的个数更是一个重要的问题。
在本文中,我们将深入探讨这个问题,以便更好地理解积分点个数对于有限元分析的影响。
1. 8节点6面体减缩积分单元的概念我们需要了解什么是8节点6面体减缩积分单元。
这个单元是常用的有限元网格单元之一,用于模拟三维空间内复杂结构的应力、应变等物理量。
它由8个节点和6个面组成,可以较好地适用于各种不规则几何形状的建模。
2. 积分点的作用在有限元分析中,积分点是用来对单元上的物理量进行数值积分的点。
积分点的选择直接影响到数值解的精度和计算效率。
选择合适数量的积分点对于有限元分析的准确性至关重要。
3. 积分点个数的选择对于8节点6面体减缩积分单元来说,积分点个数的选择需要考虑到单元内部的物理量分布情况以及计算成本。
一般来说,可以通过分析单元内部的函数的奇偶性,以及对物理量的精度要求来确定积分点的个数。
4. 影响积分点个数选择的因素在实际工程中,有些因素会影响积分点个数的选择。
如果要模拟的结构具有复杂的几何形状,可能需要更多的积分点来精确地进行积分。
另外,对于非线性分析和接触问题,积分点的选择也会有所不同。
5. 个人观点和理解在实际工程中,我们经常会遇到需要选择积分点个数的情况。
在我看来,应当根据具体的建模情况和分析要求来确定积分点个数,以便在保证精度的前提下尽量减少计算成本。
总结在本文中,我们深入探讨了8节点6面体减缩积分单元的积分点个数选择问题。
通过分析单元的特点、物理量分布情况和计算成本,我们可以更好地选择合适数量的积分点,从而保证数值解的精度和计算的高效性。
在工程实践中,合理选取积分点个数对于有限元分析结果的准确性至关重要。
只有在深入理解了8节点6面体减缩积分单元的特点和积分点个数选择的影响因素后,我们才能更好地应用于实际工程中,并取得更准确的结果。
一种多功能映射网格分方法

一种多功能映射网格分方法张少杰;陈江;李鹏;凌涛;刘浩吾【摘要】在有限元当中网格(包括影射网格和自由网格)划分对于计算结果的影响至关重要.提出一种多功能映射网格划分方法,该方法以平面的网格划分为基础,曲面和某些特殊实体的网格可以通过其投影面或基准面的网格划分数据自动生成.算例表明,该方法结构简单,适用范围广,通过合理规划可以得到质量良好的网格.【期刊名称】《现代电子技术》【年(卷),期】2008(031)008【总页数】3页(P70-72)【关键词】有限元;网格生成;四边形;六面体【作者】张少杰;陈江;李鹏;凌涛;刘浩吾【作者单位】四川大学,水电学院,四川,成都,610065;四川大学,水电学院,四川,成都,610065;四川大学建筑与环境学院,四川,成都,610065;四川大学建筑与环境学院,四川,成都,610065;四川大学,水电学院,四川,成都,610065【正文语种】中文【中图分类】TP391有限元法是当今最有效的数值计算方法之一,其应用范围涉及航空航天、土木工程、机械制造等众多领域,随着计算机技术的发展,其优越性更加突出。
一个完整的有限元计算包括前处理、求解、后处理3个部分,其中前处理约占整个计算周期的70%,求解占10%,后处理占20%,方便、快捷的前处理方式必然缩短整个计算周期以提高效率。
前处理一般包括建立模型、定义材料几何参数、网格划分、定义边界条件等,就网格划分而言,大体可分为自由网格划分和映射网格划分,目前自由网格划分的方法主要有:Delaunay三角化方法[1]、前沿生成法[2]、Paving单元法[3,4]等。
这些方法有些已用在商品化的有限元软件中,发挥了较大作用。
而有限元计算中普遍推崇的还是映射网格划分,本文提出一种多功能的映射网格划分方法,并用程序实现平面、曲面、特殊实体的一体化映射网格划分。
1 算法简介本节首先介绍平面的映射网格划分方法,然后给出曲面及特殊实体的处理方式。
无厚度六面体 8 结点单元 刚度阵

无厚度六面体 8 结点单元刚度阵1. 引言在有限元分析中,无厚度六面体 8 结点单元是一种常用的元素类型,用于模拟薄膜结构或薄壁结构的行为。
对于这种元素类型,刚度阵是非常重要的参数,它描述了该元素在有限元模拟中的应力-应变关系。
本文将首先介绍无厚度六面体 8 结点单元的基本原理,然后深入探讨其刚度阵的计算方法和意义,最后共享个人对这一主题的理解和观点。
2. 无厚度六面体 8 结点单元的基本原理无厚度六面体 8 结点单元是一种二维的有限元单元,常用于模拟薄膜或薄壁结构的行为。
它具有8个结点,每个结点有两个自由度(通常是位移),因此总共有16个刚度分量。
该单元的节点布置和形状函数的选取对于模拟结果的准确性至关重要。
在进行有限元分析时,需要建立该单元的几何属性、材料参数和边界条件等,以便对结构进行合理的模拟。
3. 无厚度六面体 8 结点单元刚度阵的计算方法刚度阵描述了有限元单元在受力情况下的应力-应变关系。
对于无厚度六面体 8 结点单元,其刚度阵的计算涉及到单元的几何属性、材料参数和形状函数等。
一般来说,刚度阵的计算可以通过数学推导和积分计算得到。
在实际工程中,通常会借助有限元软件来自动计算刚度阵,以减少人工计算的错误和工作量。
4. 无厚度六面体 8 结点单元刚度阵的意义刚度阵是描述有限元单元行为的重要参数,它可以用于计算单元在受力情况下的应力、变形和应变等信息。
通过分析刚度阵,可以了解单元在受力情况下的响应情况,进而评估结构的稳定性和安全性。
刚度阵还可以用于拟合材料模型、进行结构优化和预测结构的疲劳寿命等工作。
5. 个人观点和理解我对无厚度六面体 8 结点单元刚度阵的理解是,它是描述该单元行为的重要参数,能够直观地反映出结构在受力情况下的应力-应变关系。
在工程实践中,我们可以通过分析刚度阵来评估结构的性能,并基于此进行设计优化和决策制定。
在实际工程中,我认为熟练掌握刚度阵的计算方法和分析技巧是非常重要的,能够有效提高工程设计的准确性和效率。
一种大型三维有限元网格的显示和高效消隐方法

一种大型三维有限元网格的显示和高效消隐方法
陈良玉;杨文通
【期刊名称】《计算机辅助设计与图形学学报》
【年(卷),期】1997(009)002
【摘要】本文提出了一种面向AutoCAD的大型三限元网格的显示和高效消隐方法,可实现8-20个节点三维正确显示,并大大节省消隐的运算时间。
【总页数】4页(P164-167)
【作者】陈良玉;杨文通
【作者单位】东北大学机械工程学院;东北大学机械工程学院
【正文语种】中文
【中图分类】TP391.41
【相关文献】
1.利用OpenGL实三维有限元网格图的消隐显示 [J], 徐全生;冯艳君
2.适用于三维有限元网格显示的逐步消隐法 [J], 詹梅;刘郁丽
3.一种三维有限元网格的快速消隐算法 [J], 杨小辉;方宗德;杨青
4.基于B—Rep的三维有限元网格高效消隐技术 [J], 王志强;洪嘉振
5.三维有限元网格自动剖分及其高效消隐技术研究 [J], 袁海平;付玉华;杨国春因版权原因,仅展示原文概要,查看原文内容请购买。
三维有限元网格自动剖分及其高效消隐技术研究

三维有限元网格自动剖分及其高效消隐技术研究[摘要]对三维多材质任意边界条件结构有限元网格的生成技术与消隐技术进行了探讨,建立了一套行之有效的三维有限元自动剖分通用算法和远远优于常规消隐方法的高效消隐技术,并编写了相应的通用计算程序。
应用结果表明,其能够高效快速简便地达到理想的剖分网格和消隐效果,确保了计算数据的正确可靠,大大提高了计算效率。
[关键词]三维有限元,自动剖分,等参单元,超单元,消隐技术1 引言常规有限元分析计算中前处理工作繁琐而艰巨,且随着数据输入量的增大极易出错。
对于复杂结构,网格生成的工作已成为有限元数值模拟工作的一个瓶颈[1,2],如何简化前处理工作、实现网格自动生成对于推广有限元分析至关重要。
目前已有许多成熟的方法用于规则的或部分非规则的几何图形的剖分,但功能强大的自动剖分算法成果多捆绑于大型的有限元软件中,这与有限元计算程序的广泛应用极不协调。
同时,立体消隐图的绘制是前后处理技术的一个重要部分[3],为便于检查剖分情形及显示计算结果,需要对模型的剖分网格进行消隐处理,使绘出的网格图具有较强的立体感。
目前比较常用的消隐算法有画家算法[4]、Z- buffer算法、扫描线算法、Warnock算法等,各种消隐算法各具特点,适用于不同的情况。
但对于动辄达十几万单元的网格处理,都是一个沉重的负担。
经研究,本文提出了一种实用的三维有限元自动剖分通用算法,并同时在原有消隐技术的基础上建立了远远优于常规消隐方法的超大单元面高效消隐判断技术,大大提高了常规有限元前处理和消隐计算的效率。
2 空间六面体超单元变换剖分原理根据结构的几何外形、内部材料类型以及网格的布置要求,将复杂三维几何结构划分为若干子区域,这些子区域称为超单元[5]。
采用等参变换的方法,求出超单元内各网格结点的整体坐标,此即为超单元变换的基本原理。
图1 空间六面体8结点等参单元Fig.13D 8-node Isoparametric Element对于图1所示的空间六面体8结点等参单元,其单元内任一点的整体坐标变换公式为3 相邻超单元边界重复结点的消除有限元网格划分要求具有不间断和不重复的连接单元和结点编号,因此相邻超单元边界重复结点必须进行消除处理,并统一编号。
8节点6面体减缩积分单元的积分点个数

8节点6面体减缩积分单元的积分点个数这篇文章将探讨8节点6面体减缩积分单元的积分点个数。
我们将从简单的概念开始,逐步深入,以便全面地理解该主题。
1. 介绍8节点6面体减缩积分单元是有限元法中常用的积分单元之一。
它由8个节点和6个面组成,可以用于模拟各种结构的力学行为。
在使用该单元进行数值计算时,我们需要确定积分点的个数,以提高计算的准确性和效率。
2. 理论背景在有限元法中,积分点的个数对计算结果的精度有着重要的影响。
通常情况下,更多的积分点可以提高计算的准确性,但同时也会增加计算的复杂性和耗时。
确定适当的积分点个数是一个关键的问题。
3. 积分点个数的选择为了确定8节点6面体减缩积分单元的积分点个数,我们可以考虑以下几个因素:a. 问题的性质:不同的问题对积分点个数的要求不同。
对于像静力学问题这样的简单情况,较少的积分点个数可能已经足够;而对于动力学或热传导等更复杂的问题,可能需要更多的积分点。
b. 精度要求:如果需要更高的计算精度,我们可能需要增加积分点的个数。
c. 计算资源:考虑到计算资源的限制,我们需要在计算精度和计算效率之间进行权衡。
通常情况下,我们希望尽可能减少计算时间,同时保持结果的准确性。
4. 个人观点和理解在对8节点6面体减缩积分单元的积分点个数进行选择时,我倾向于综合考虑问题的性质、精度要求和计算资源等因素。
我认为在一般情况下,使用8个积分点已经可以获得合理的计算结果。
如果问题比较复杂或精度要求较高,可以考虑增加积分点的个数。
5. 总结与回顾在本文中,我们讨论了8节点6面体减缩积分单元的积分点个数选择的问题。
我们通过综合考虑问题的性质、精度要求和计算资源等因素,得出了使用8个积分点的合理选择。
然而,根据具体情况,也可以考虑增加积分点的个数以提高计算精度。
8节点6面体减缩积分单元的积分点个数是一个关键的选择,它直接影响着数值计算的精度和效率。
在实际应用中,我们需要考虑多种因素来确定合适的积分点个数。
六面体单元生成方法及相关技术

上海交通大学博士学位论文
面体单元)、ICEM CFD(ICEM CFD Engineering,四边形和六面体单元)、QulkGrid(Perspective Edge,四边形单元)、TrueGrid(XYZ Scientific Applications,Inc.,四边形和六面体单元)、VGM (NASA,Langley & Lockhead,四边形和六面体单元)等。
( ) Nξ2
(1+ηηi )(1+ ζζi ),
i
=
9,11, 17, 19
( ) Ni(ξ,η,
ζ)=
1 4
(1+ ξξi )1−η 2
(1+ ζζi ),
i
= 10,12, 18, 20
(3-7)
( ) Ni(ξ,η,
ζ
)
=
1 4
(1+ ξξi
)(1 + ηηi
(3-4)
α = xη2 + yη2, β = xξ xη + yξ yη , γ = xξ2 + yξ2
和代数法相比,偏微分方程法的计算较为复杂,需要求解偏微分方程,右 端强迫函数的选取 也不方便,可以通过强迫函数的选取得到希望的映射网格。
·21·
上海交通大学博士学位论文
3.2.3 超单元映射法
第3章 六面体单元生成方法及相关技术
3.1 引言
有限元法是求解偏微分方程描述的连续体问题的一种近似工程方法。为克服实际连续体问题 难于处理的问题,它将分析区域离散化,将偏微分方程转化为线性方程组,采用数值计算方法求 出连续体问题的近似解。用有限元法进行工程分析的主要过程包括三个阶段:(1)有限元模型的 建立和数据准备;(2)用软件分析计算;(3)分析结果的判断和评定。迅速而合理地划分有限元 网格是完成有限元分析的前提和保证。
有限元网格剖分方法概述

有限元网格剖分方法概述在采用有限元法进行结构分析时,首先必须对结构进行离散,形成有限元网格,并给出与此网格相应的各种信息,如单元信息、节点坐标、材料信息、约束信息和荷载信息等等,是一项十分复杂、艰巨的工作。
如果采用人工方法离散对象和处理计算结果,势必费力、费时且极易出错,尤其当分析模型复杂时,采用人工方法甚至很难进行,这将严重影响高级有限元分析程序的推广和使用。
因此,开展自动离散对象及结果的计算机可视化显示的研究是一项重要而紧迫的任务。
有限元网格生成技术发展到现在, 已经出现了大量的不同实现方法,列举如下:映射法映射法是一种半自动网格生成方法,根据映射函数的不同,主要可分为超限映射和等参映射。
因前一种映射在几何逼近精度上比后一种高,故被广泛采用。
映射法的基本思想是:在简单区域内采用某种映射函数构造简单区域的边界点和内点,并按某种规则连接结点构成网格单元。
也就是根据形体边界的参数方程,利用映射函数,把参数空间内单元正方形或单元三角形(对于三维问题是单元立方体或单元四面体)的网格映射到欧氏空间,从而生成实际的网格。
这种方法的主要步骤是,首先人为地把分析域分成一个个简单可映射的子域,每个子域为三角形或四边形,然后根据网格密度的需要,定义每个子域边界上的节点数,再根据这些信息,利用映射函数划分网格。
这种网格控制机理有以下几个缺点:(1)它不是完全面向几何特征的,很难完成自动化,尤其是对于3D区域。
(2)它是通过低维点来生成高维单元。
例如,在2D问题中,先定义映射边界上的点数,然后形成平面单元。
这对于单元的定位,尤其是对于远离映射边界的单元的定位,是十分困难的,使得对局部的控制能力下降。
(3)各映射块之间的网格密度相互影响程度很大。
也就是说,改变某一映射块的网格密度,其它各映射块的网格都要做相应的调整。
其优点是:由于概念明确,方法简单,单元性能较好,对规则均一的区域,适用性很强,因此得到了较大的发展,并在一些商用软件如ANSYS等得到应用。
一种ECVT有限元自动剖分方法

一种ECVT有限元自动剖分方法邹星星;赵进创;傅文利;马增强【摘要】ECT三维传感器场域剖分是当前ECT直接三维成像研究的关键问题,提出一种三维有限元自动剖分方法并得到结果验证.采用四面体作为剖分单元,利用Delaunay四面体剖分算法反复迭代实现三维实体自动剖分;为清晰显示剖分图像,结合画家算法对网格进行消隐处理.仿真结果表明,随着迭代次数增加,剖分精度提高,其剖分时长也会相应增加.根据需求剖分精度设置迭代次数,能有效的剖分传感器场域,实现了三维场域自动剖分,为计算三维敏感场从而提高ECT三维成像精度奠定基础.%As three-dimensional sensor field mesh of ECT is a key issue in direct 3D ECT currently,a method of 3D FEM self-organizing mesh generation is proposed.The method adopts tetrahedron as split unit,and uses Delaunay triangulation algorithm for auto meshing 3D solid.Finally,thepainter algorithm is used to process grid blanking.Simulation results show that with the increase in the number of iterations,mesh accuracy is improved while more time is spent.Setting iteration according to the precision of mesh can effectively mesh the sensor field,achieving the three-dimensional auto meshing,which provides basis for calculating the three-dimensional sensitivity field in order to improve the accuracy of three-dimensional imaging of ECT.【期刊名称】《石家庄铁道大学学报(自然科学版)》【年(卷),期】2017(030)003【总页数】5页(P106-110)【关键词】电容层析体积成像;有限元;自动剖分【作者】邹星星;赵进创;傅文利;马增强【作者单位】广西大学计算机与电子信息学院,广西南宁530004;广西大学计算机与电子信息学院,广西南宁530004;广西大学计算机与电子信息学院,广西南宁530004;石家庄铁道大学电气与电子工程学院,河北石家庄050043【正文语种】中文【中图分类】TP391电容层析成像技术(Electrical Capacitance Tomography)基于奥地利科学家J.Radon于1917年提出的图像重建理论(即Radon变换及其逆变换),通过放置在被测管道周围的传感器阵列获取被测场域不同观测角度下的多个电容测量值(即ECT正问题),以此为投影数据,采用适当的图像重建算法重建被测区域内介质分布(即ECT逆问题)。
八节点六面体单元的变形重构

八节点六面体单元的变形重构八节点六面体单元是一种常用的有限元分析单元,广泛应用于工程领域。
其变形重构是指在八节点六面体单元的基础上进行形状优化,以改善其力学性能和应变分布。
本文将介绍八节点六面体单元的基本原理和应用,探讨变形重构的方法和效果。
1、八节点六面体单元的基本原理八节点六面体单元是一种具有六个面的立方体单元,通过节点和单元之间的连接关系来模拟结构的变形。
它具有较高的自由度,能够较准确地描述结构的力学行为。
在有限元分析中,常用于模拟复杂结构的变形和应力分布。
2、八节点六面体单元的应用八节点六面体单元广泛应用于航空航天、建筑、汽车等工程领域。
例如,在飞机机翼结构的有限元模型中,八节点六面体单元可以用来模拟机翼的复杂变形和应力分布,从而对机翼的强度和刚度进行分析和优化。
3、变形重构的方法为了进一步提高八节点六面体单元的性能,可以进行变形重构。
变形重构的方法主要包括形状优化和拓扑优化。
形状优化通过调整单元的形状来改善其力学性能,如减小应力集中、增加强度等。
拓扑优化则通过改变单元之间的连接关系来优化结构的刚度和轻量化。
4、形状优化形状优化是指通过调整单元的形状来改善其力学性能。
可以通过改变单元的体积、几何形状等参数来实现。
常用的形状优化方法包括有限元逆设计、拟变形法等。
这些方法通过迭代计算,找到最优的形状参数,以达到力学性能的最佳化。
5、拓扑优化拓扑优化是指改变单元之间的连接关系来优化结构的刚度和轻量化。
拓扑优化方法主要包括有限元拓扑优化、拓扑递归算法等。
通过这些方法,可以调整单元的连接关系,使得结构在满足强度和刚度要求的同时,尽量减少材料的使用量。
6、变形重构的效果通过对八节点六面体单元进行变形重构,可以显著改善其力学性能和应变分布。
形状优化可以减小应力集中,提高结构的强度和刚度;拓扑优化可以实现结构的轻量化,降低结构的质量和成本。
变形重构后的八节点六面体单元模型可以更准确地模拟结构的变形和应力分布,为工程设计提供更可靠的基础。
有限元四面体及六面体单元

空间问题有限元分析
(4-109)
(4-110)
基本概念4节点四面体单元应变场的表达
空间问题有限元分析
基本概念
4节点四面体单元应力场的表达
(4-111)
基本概念
空间问题有限元分析
4节点四面体
5.单元的刚度矩阵及节点等效载荷矩阵
6.单元刚度方程
(4-112)
(4-113)
(4-104)
(4-105)
基本概念
空间问题有限元分析
4节点四面体
2.单元位移场的表达
将式(9-3)代入节点条件(9-4)中,可求取待定系数(ai,bi,ci),i=0,1,2,3。在求得待定系数后,可重写式(9-3)为
(4-106)
(4-107)
空间问题有限元分析
基本概念
4节点四面体单元应变场的表达
空间问题有限元分析
(4-102)
(4-103)
基本概念4节点四面体4节点四面体单元几何和节点描述
基本概念
空间问题有限元分析
4节点四面体
2.单元位移场的表达
该单元有4个节点,单元的节点位移有12个自由度(DOF)。因此每个方向的位移场可以设定4个待定系数,根据节点个数以及确定位移模式的基本原则(从低阶到高阶的完备性、唯一确定性),选取该单元的位移模式为
基本概念
4节点四面体4节点四面体单元的常系数应变和应力
空间问题有限元分析
基本概念8节点正六面体1. 单元的几何和节点描述 该单元为由8节点组成的正六面体单元(hexahedron element),每个节点有3个位移(即3个自由度),单元的节点及节点位移如图所示(4-115) (4-116)
复杂机械零件的六面体有限元网格生成方法_曹华军

1/ 2
1/ 2
1/ 2
m )T C ( )dV ve ( e 1
(6)
式中,m 为整体结构的单元数量。 对于某确定问题的有限元计算,可以采用总体 误差指标 和单元误差指标 e 。两个误差满足如下 精度准则
1
有限元网格的误差估算理论
应力改进值 泛函数变分为
v
1 )T N i dV 0 ( i e 2 i 1, 2, , ne
2
(2)
2.1
模型关键区域识别与计算域规划
模型关键区域识别
对 ne 个方程组进行求解即可得到节点的应力 恢复值,从而计算出整体计算域内的应力恢复值。 位移近似解与精确解的误差 eu 定义为
( R) =
( Ri ) i 1
m
i=1,2,...,m
m为结构单元应力恢复使应力在计算域内具有与位移相同的连续性是单元内部应力恢复后的改进值其求解问题可描述为ne在针对具体有限元模型进行仿真计算时基于有限元网格误差估算理论可分析关注区域的误差特性及其导致误差的因素进而优化网格划分策略获得高效高质量网格
第 50 卷第 15 期 2014 年 8 月
机
械
工 程
Hexahedral Mesh Generation Method for Complex Mechanical Structure
CAO Huajun SHU Linsen XU Lei LI Hao
(The State Key Laboratory of Mechanical Transmission, Chongqing University, Chongqing 400030)
(a) #形策略
有限元网格的拓扑分析

n mb rn f en d s drd cn eb d d f ewh l t n s t x u e go t o e e u igt a wit o t o esi e s i h n a h n h h f mar . i Ke r s c mp tr p l a o ; ii lm e t t o ;o oo ya ay i; s ywo d : o u e pi t n f t ee n h d tp lg n lss meh a ci n e me
JANG o gx n W AN ig r n I Xin — i, P n —o g
( I h I ies y Na e a g i g i 3 0 9 C ia Na c a g v ri , n h n a x 0 2 , hn ) l l Un t Jn 3
Ab ta t sr c :M e h g n rt n a d n d —u eig o t z t n ae k y tc nq e n f i s e e ai n o en mb rn p miai r e e h iu si i t o i o ne ee e t n l ss I ii lme t n l ss teb d d ftewh l t n s ti e u e lm n ay i.nfnt ee n ay i,h a wi t o oesi esmarxi rd c d a e a n h h f s b p miig ten d u eig t ee yt en e e m oy s a ei sv da dtec c lt n y o t zn o en mb r , rb e d dme r p c s a e a u ai i h n h h n h l o i t srd c d T p l gc l eain hp o eee n s e sa ay e . eu eu to s mei e u e . o oo ia lt s i ft lme tmeh si n lz d Th s f l h d r o h me frc n tu t g ee e - o e tp lg ar d a smb ig te wh l o oo y mar y o o s c n lm n n d o oo y m ti a se l oe tp lg t x ae r i t xn n h i
《有限元基础教程》_【MATLAB算例】4.8.2(1) 基于8节点六面体单元的空间块体分析(Hexahedral3D8Node)

【MATLAB 算例】4.8.2(1) 基于8节点六面体单元的空间块体分析(Hexahedral3D8Node)如图4-23所示的一个空间块体,在右端部受两个集中力F 作用,其中的参数为:105110Pa,=0.25,=0.2m,=110N E t F μ=⨯⨯。
基于MATLAB 平台,用一个空间8节点六面体单元计算各个节点位移、支座反力以及单元的应力。
(a) 问题描述 (b) 有限元分析模型图4-23 右端部受集中力作用的空间块体解答:对该问题进行有限元分析的过程如下。
(1)结构的离散化与编号将结构离散为一个8节点六面体单元,节点编号如图4-23(b)所示,节点的几何坐标见表4-13。
表4-13 节点的坐标节点 节点坐标/m x y z 1 0.2 0 0 2 0.2 0.8 0 3 0 0.8 0 4 0 0 0 5 0.2 0 0.6 6 0.2 0.8 0.6 7 0 0.8 0.6 80.6节点位移列阵T 111222888(241)[]e u v w u v w u v w ⨯=q (4-194)总的节点载荷列阵T 111222888(241)[]e x y z x y z x y z P P P P P P P P P ⨯=P (4-195)其中,节点外载567110z z P P F N ==-=-⨯;支反力为11x x P R =,11y y P R =,11z z P R =,44x x P R =,44y y P R =,44z z P R =,55x x P R =,55y y P R =,55z z P R =,88x x P R =,88y y P R =,88z z P R =;其余节点载荷分量为零。
(2)计算单元的刚度矩阵(以国际标准单位)首先在MA TLAB 环境下,输入弹性模量E 和泊松比NU ,然后针对题中单元节点坐标,调用函数Hexahedral3D8Node_Stiffness ,就可以得到单元的刚度矩阵k1 (24×24)。
有限元网格划分方法与基本原理

结构有限元分析中的网格划分技术及其应用实例结构有限元分析中的网格划分是否直接关系到解算的效果。
本文简述了网格划分应用的基本理论,并以空间自由曲面覆盖件和大型整体网络钢筋壳体产品的有限元分析中的网格划分为实例对象,详细讲述了空间自由和三维实体的网格划分基本理论及其在工程中的实际应用,非常具有现实意义和借鉴价值。
一、前言有限元网格划分是进行有限元数值模拟分析至关重要的一步,它直接影响着后续数值计算分析结果的精确性。
网格划分涉及单元的形状及其拓扑类型、单元类型、网格生成器的选择、网格的密度、单元的编号以及几何体素。
从几何表达上讲,梁和杆是相同的,从物理和数值求解上讲则是有区别的。
同理,平面应力和平面应变情况设计的单元求解方程也不相同。
在有限元数值求解中,单元的等效节点力、刚度矩阵、质量矩阵等均用数值积分生成,连续体单元以及壳、板、梁单元的面内均采用高斯(Gauss)积分,而壳、板、梁单元的厚度方向采用辛普生(Simpson)积分。
辛普生积分点的间隔是一定的,沿厚度分成奇数积分点。
由于不同单元的刚度矩阵不同,采用数值积分的求解方式不同,因此实际应用中,一定要采用合理的单元来模拟求解。
CAD软件中流行的实体建模包括基于特征的参数化建模和空间自由曲面混合造型两种方法。
Pro/E和S oildWorks是特征参数化造型的代表,而 CATIA与Unigraphics等则将特征参数化和空间自由曲面混合造型有机的结合起来。
现有CAD软件对表面形态的表示法已经大大超过了CAE软件,因此,在将CAD实体模型导入CAE软件的过程中,必须将CAD模型中其他表示法的表面形态转换到CAE软件的表示法上,转换精度的高低取决于接口程序的好坏。
在转换过程中,程序需要解决好几何图形(曲线与曲面的空间位置)和拓扑关系(各图形数据的逻辑关系)两个关键问题。
其中几何图形的传递相对容易实现,而图形间的拓扑关系容易出现传递失败的情况。
数据传递面临的一个重大挑战是,将导入CAE程序的CAD模型改造成适合有限元分析的网格模型。
高精度多用途八结点六面体单元

高精度多用途八结点六面体单元
鹿晓阳;刘玉文
【期刊名称】《机械工程学报》
【年(卷),期】1989(25)2
【总页数】8页(P67-74)
【关键词】八结点六面体;非协调等参元;有限元
【作者】鹿晓阳;刘玉文
【作者单位】
【正文语种】中文
【中图分类】TB115
【相关文献】
1.空间八节点六面体单元中绘制等值面的两种算法 [J], 刘亚东;彭颖红;李从心
2.空间八节点六面体单元中绘制等值面的两种算法 [J], 刘亚东;彭颖红;李从心
3.八结点任意六面体的杂交/混合元 [J], 姚敬之
4.20结点六面体单元边中点的自动形成和优化 [J], 姚敬之
5.基于六面体八节点单元的三维有限元数值模拟表面等值线的生成算法 [J], 李昱明;周飞;郭新贵;彭颖红
因版权原因,仅展示原文概要,查看原文内容请购买。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
通过 &()*%&’ 直接打开该 ’>? 文件, 然后利用 用该系统将三维有限元网格数据以 H’ 表面的数据格式生成 ’>? 文件。 从而提高了大量数据的显示能力, 并且省去了消隐的复杂 &()*%&’ 提供的消隐功能对三维有限 元 网 格 进 行 自 动 消 隐 , 编程工作。 关键词 三维有限元网格
!顺序应与 (IP , ; # IP:B=.I T !"U=^@IU=T , PFK@XYG ) !#! S ( ) " (IP , ; # IP:B=.I T ("U=^@IU=T , PFK@XYH)
图$ 任意六面体 图! 多边形网络 (IP , IP:B=.I T %"U= ; &’U=T )
b"
!"""#$$ 计算机工程与应用
图.
六面体网格表示的未经消隐的工件实体
$
引言
有限元仿真技术目前已广泛应用于金属塑性加工领域, 它
!
描述三维网格的图形交换文件格式
可以提供变形体在成形过程中的各种物理量场, 分析其变化规 律, 预测成形缺陷等。 这对指导工业生产具有十分重要的意义。 但是, 有限元数值计算的结果是一些大量的数据, 必须对这些 数据进行处理, 并将其可视化后才能进行有效的分析。 科学计算可视化是一门融计算机图形学、 图象处理和计算 机视觉综合应用于计算科学的学科, 该学科近年来已取得很大 发展。在科学计算可视化中, 图形的消隐技术占有举足轻重的 地位, 同时它也是有限元网格可视化处理的关键技术之一。因 此, 很有必要研究有限元仿真可视化处理中的自动消隐技术。 目 前 有 许 多 %&’ 软 件 可 以 直 接 进 行 可 视 化 和 自 动 消 隐 , 如 &()*%&’ 、 +,--、 %&.-&、 -’/&0 、 12*3/ 等 , 这 些 软 件 都 具 有 二次开发工具。选择一种先进的图形显示系统, 在其基础上进 行开发, 方便而高效地来显示有限元仿真结果是实现有限元后 处理的捷径。文章选取使用最为广泛的 &()*%&’ 软件作 为 开 发平台。 八节点六面体三维有限元网格通常用该六面体的八个节 点 及 附 加 的 面 的 信 息 来 描 述 。 在 &()*%&’ 中 可 以 用 简 单 的 但用这种方法产 0%4 文件的 15676-8/ 格式实现线框的显示。 (9:;< ) 进行消 生的网格不能直接用 &()*%&’ 提 供 的 消 隐 命 令 隐, 图形显示不直观; 而 且 由 于 0%4 是 &()*%&’ 的 批 处 理 文 件, 对于显示具有大量数据的有限元网格很费时间。 在 &()*= 消隐命令适用于具有面或体 属 性 的 实 体 , 而它提供的 %&’ 中, 开放的 ’>? 文件可以包含图形的所有 信 息 @$A。 另 一 方 面 , ’>? 文件存储结构灵活,只要将图形信息按一定格式写入文件, 就 可直接显示。因 此 , ’>? 是 图 形 构 造 和 图 形 显 示 的 一 种 快 速 、 高效的方法,建立有限元网格在 ’>? 文件中的具有面或体属 性的实体存储 ’>? 格式是实施有限元网格可视化的关键。
&()*%&’ 4$B 支 持 文 本 格 式 和 二 进 制 的 ’>? 文 件 @!A。 其 中文本格式的 ’>? 文件具有操作方便、 易于识别等特点, 可以
直接用高级编程语言生成。
!#$
’>? 文件的一般结构
(,45+1) , 每个 组 占 两 行 , ’>? 文件的最小组成单位是组 , 是一个整数, 第二行是组值 第一行是组码 (,45+1 %5’/ ) , 它的类型取决于组码的数值。 (,45+1 C&6+/ ) (0/%.-58 ) , 每个节由多个组构成。 ’>? 文 件 分 为 多 个 节 典型的 ’>? 文件由标题节、 类节、 表节、 块节、 实体节、 对象节、 /5? 标 志 组 成 。 最 简 单 的 ’>? 文 件 仅 包 含 实 体 节 和 /5? 标 志。
B=. B ; 9.:-;. <=/>? EPFK@ ;
参见 $#! S S 定义三维网格实体, (IP , ; IP:B=.I T "U2P/@G@B=?U=T ) (IP , ; IP:B=.I T 7U="U=T ) (IP , ; IP:B=.I T &&U= $U=T ) (IP , ; IP:B=.I T $"U="#"U=T ) (IP , ; IP:B=.I T !"U="#"U=T ) (IP , ; IP:B=.I T ("U="#"U=T ) (IP , IP:B=.I T %"U= ; $&U=T ) (IP , ; IP:B=.I T %$U= &U=T ) (IP , ; IP:B=.I T %!U= (U=T ) 参见 !#! S S 定义三维网格实体的各个节点, (B4" ; BV$7 ; BWW ) I/:
八节点六面体有限元网格自动消隐技术
李昱明 周 飞 郭新贵 彭颖红 上海交通大学塑性成形工程系 (上海 !"""H" )
/3KL:MNOP<QK*RS(TM:UB#O)L#Q<)#UQ
摘 要 文章基于图形消隐技术, 利用 % 语言, 开发了一套有限元网格与 &()*%&’ 之间进行图形交换的软件系统, 并利
格式的研究, 将三维有限元网格数据以三维网格的数据形式生 利 用 ,-./0,) 的 aB>? 命 令 自 动 对 其 进 行 消 隐 , 成 )*+ 文 件 , 大大改善了有限元分析的后处理能力, 提高了后处理程序的开 发效率。该方法已用于作者开发的有限元分析软件中, 实践证 明, 该方法十分有效。 同时, 它也表明自动消隐技术在有限元网 格可视化过程中起着十分重要的作用。 (收稿日期: $bbb 年 $! 月)
$
一致
C:?KL ;
由于六面体沿棱边展开有多种方案, 因此展开图有 1324 &3( 、 53’ 等不同形式。 ,-./0,) 对这些网格的节点按列处理。 这样就得到了一个节点序列来进行存储。以图 ! 为例, 节点序 列 为 6& , !, !, 5, $, $, 7, ’, $, %, (, !, %, %, &, 7, 7, 58。 当 然 方 案 有 很多种,可以从中选取数据量较少的一种来表达凸三维六面 体。
构成三维网格实体的点序列
( 任意凸三维六面体的 )*+ 描述 (#$ 任意凸三维六面体的 )*+ 描述
六 面 体 在 ,-./0,) 图 形 数 据 库 中 用 () 表 面 来 表 示 , 而 ,-./0,) 现 有 的 三 维 实 体 作 图 命 令 不 能 表 达 不 规 则 的 三 维 多 文 面体, 因而不能用 ,-./0,) 的批处理命令作图。通过研究, 章提出了一种可描述 任 意 三 维 六 面 体 (规 则 或 不 规 则 ) 的 )*+ 格式。 这种格式可用于八节点六面体的直接消隐, 将其推广, 还 可以用于任意多面体。 任意凸三维六面体具有八个节点。这些节点四四组合, 组 成六个曲面或平面实体 (以下统称曲面) , 这六个三维曲面封闭 形成凸六面体。 ,-./0,) 用一个三维多边形网格来表达这些曲 面。 这个网格由 1 和 2 两个交叉方向的网格线形成, 其中每一 个网格表达一个曲面实体, 网格线为六面体的棱边。 可得到如图 ! 例如, 将如图 $ 所示的六面体沿棱边展开, 所示的网格。不同的展开方式得到的网格不同。图中实线表示 棱边。 由于网格的规则性, 用虚线将缺口补齐, 得到一个规则的 网格图, 图 ! 是一个 &3( 的节点示意图。
元网格的 )*+ 文件。在 ,-./0,) ‘$’ 中, 用 )*+M2 命令直接 打开该文件, 将显示出如图 ( 所示的六面体网格表示的工件实 体。然后使用 aM)O 命令进行消隐得到如图 ’ 的消隐后的工件 实体。
5
结论
文章通过对三维八节点六面体单元形体表示的 )*+ 文件
其中, <?@?A 结 构 描 述 有 限 元 单 元 的 结 构 体 , <=/>? 结 构 描述一个单元中的各个节点。
!A9#65-#: [LO<; *Q )P< X2LSP:UO P:;:QX )<UPQ*M*XW, L OWO)<K :O ;<Z<M*S<; )* )2LQOV<2 )P< H’ ?/J J<OP :Q)* )P< &(= )*%&’ OWO)<K TW % MLQX(LX< )P2*(XP ’>? V*2KL)\ ]P<Q )P< H’ ?/J K<OP ;L)L L2< O)*2<; :Q L ’>?, :) ULQ T< L()*3 P:;;<Q TW (O:QX 9-’/ V(QU):*Q *V &()*%&’# B,0C$6(9: H’ ?/J, ’>?, &()*3P:;< , C:O(LM:^L):*Q
列以 0/G/8’ 结束。
$" , !" , H" I"
组码 $" 、 组码 H" 表示多义线 !" 组值为 " , 的高度。 多义线标志, 若组值为 $E , 表示三维多边形 网格。