下载文档原格式
1网格数量网格数量的多少将影响计算结果的精度和计算规模的大小。
一般来讲，网格数量增加，计算精度会有所提高，但同时计算规模也会增加，所以在确定网格数量时应权衡两个因数综合考虑。
网格较少时增加网格数量可以使计算精度明显提高，而计算时间不会有大的增加。
当网格数量增加到一定程度后，再继续增加网格时精度提高甚微，而计算时间却有大幅度增加。
所以应注意增加网格的经济性。
实际应用时可以比较两种网格划分的计算结果，如果两次计算结果相差较大，可以继续增加网格，相反则停止计算。
在决定网格数量时应考虑分析数据的类型。
在静力分析时，如果仅仅是计算结构的变形，网格数量可以少一些。
如果需要计算应力，则在精度要求相同的情况下应取相对较多的网格。
同样在响应计算中，计算应力响应所取的网格数应比计算位移响应多。
在计算结构固有动力特性时，若仅仅是计算少数低阶模态，可以选择较少的网格，如果计算的模态阶次较高，则应选择较多的网格。
在热分析中，结构内部的温度梯度不大，不需要大量的内部单元，这时可划分较少的网格。
2网格疏密网格疏密是指在结构不同部位采用大小不同的网格，这是为了适应计算数据的分布特点。
在计算数据变化梯度较大的部位(如应力集中处)，为了较好地反映数据变化规律，需要采用比较密集的网格。
而在计算数据变化梯度较小的部位，为减小模型规模，则应划分相对稀疏的网格。
这样，整个结构便表现出疏密不同的网格划分形式。
采用疏密不同的网格划分，既可以保持相当的计算精度，又可使网格数量减小。
因此，网格数量应增加到结构的关键部位，在次要部位增加网格是不必要的，也是不经济的。
划分疏密不同的网格主要用于应力分析(包括静应力和动应力)，而计算固有特性时则趋于采用较均匀的钢格形式。
这是因为固有频率和振型主要取决于结构质量分布和刚度分布，不存在类似应力集中的现象，采用均匀网格可使结构刚度矩阵和质量矩阵的元素不致相差太大，可减小数值计算误差。
同样，在结构温度场计算中也趋于采用均匀网格。
有限元的网格划分技术对于有限元分析来说，网格划分是其中最关键的一个步骤，网格划分的好坏直接影响到解算的精度和速度。
网格化有三个步骤：定义单元属性（包括实常数）、在几何模型上定义网格属性、划分网格。
定义网格的属性主要是定义单元的外形、大小。
单元大小基本上在线段上定义，可以用线段数目或长度大小来划分，可以在线段建立后立即声明，或整个实体模型完成后逐一声明。
采纳BottOm-UP方式建立模型时，采纳线段建立后立即声明比较便利且不易出错。
例如声明线段数目和大小后，叁制对象时其属性将会一•起夏制，完成上述操作后便可进行网格化命令。
网格化过程也可以逐步进行，即实体模型对象完成到某个阶段就进行网格话，如所得结果满足，则连续建立其他对象并网格化。
网格的划分可以分为自由网格（free meshing）、映射网格（mapped meshing）和扫略网格（SWeeP meshing）等。
一、自由网格划分自由网格划分是自动化程度最高的网格划分技术之一,它在面上可以自动生成三角形或四边形网格，在体上自动生成四周体网格。
通常状况下，可采用ANSYS的智能尺寸掌握技术（SMARTSIZE命令）来自动掌握网格的大小和疏密分布，也可进行人工设置网格的大小（AESIZE、LESIZE、KESIZE、ESIZE等系列命令）并掌握疏密分布以及选择分网算法等（ MOPT 命令）。
对于简单几何模型而言，这种分网方法省时省力，但缺点是单元数量通常会很大，计算效率降低。
同时，由于这种方法对于三维简单模型只能生成四周体单元，为了获得较好的计算精度，建议采纳二次四周体单元（92号单元）。
假如选用的是六面体单元，则此方法自动将六面体单元退化为阶次全都的四周体单元，因此，最好不要选用线性（•阶次）的六面体单元（没有中间节点，比如45号单元）,由于该单元退化后为线性的四周体单元，具有过大的刚度，计算精度较差；假如选用二次的六面体单元（比如95 号单元）,由于其是退化形式，节点数与其六面体原型单元全都，只是有多个节点在同一位置而己,因此,可以采用TCHG命令将模型中的退化形式的四周体单元变化为非退化的四周体单元（如92号单元），削减每个单元的节点数量，提高求解效率。
有限元网格划分摘要:总结近十年有限元网格划分技术发展状况。
首先,研究和分析有限元网格划分的基本原则；其次,对当前典型网格划分方法进行科学地分类,结合实例,系统地分析各种网格划分方法的机理、特点及其适用范围,如映射法、基于栅格法、节点连元法、拓扑分解法、几何分解法和扫描法等；再次,阐述当前网格划分的研究热点,综述六面体网格和曲面网格划分技术；最后,展望有限元网格划分的发展趋势。
关键词:有限元网格划分；映射法；节点连元法;拓扑分解法；几何分解法；扫描法；六面体网格1 引言有限元网格划分是进行有限元数值模拟分析至关重要的一步，它直接影响着后续数值计算分析结果的精确性。
网格划分涉及单元的形状及其拓扑类型、单元类型、网格生成器的选择、网格的密度、单元的编号以及几何体素。
在有限元数值求解中，单元的等效节点力、刚度矩阵、质量矩阵等均用数值积分生成，连续体单元以及壳、板、梁单元的面内均采用高斯（Gauss）积分，而壳、板、梁单元的厚度方向采用辛普生（Simpson）积分。
2 有限元网格划分的基本原则有限元方法的基本思想是将结构离散化,即对连续体进行离散化,利用简化几何单元来近似逼近连续体,然后根据变形协调条件综合求解。
所以有限元网格的划分一方面要考虑对各物体几何形状的准确描述,另一方面也要考虑变形梯度的准确描述。
为正确、合理地建立有限元模型,这里介绍划分网格时应考虑的一些基本原则。
2.1 网格数量网格数量直接影响计算精度和计算时耗,网格数量增加会提高计算精度,但同时计算时耗也会增加。
当网格数量较少时增加网格，计算精度可明显提高,但计算时耗不会有明显增加;当网格数量增加到一定程度后,再继续增加网格时精度提高就很小,而计算时耗却大幅度增加。
所以在确定网格数量时应权衡这两个因素综合考虑。
2.2 网格密度为了适应应力等计算数据的分布特点,在结构不同部位需要采用大小不同的网格。
在孔的附近有集中应力,因此网格需要加密;周边应力梯度相对较小,网格划分较稀。
机械设计中有限元分析的几个关键问题在机械设计中，有限元分析是一种常用的工具和方法。
它可以帮助工程师们对机械结构进行仿真和分析，评估其性能和可靠性，优化设计方案，减少试验成本和开发周期。
在进行有限元分析时，也存在一些关键问题需要注意和解决。
下面将介绍几个常见的有限元分析的关键问题。
1. 网格划分：网格划分是有限元分析的第一步，也是最关键的一步。
合理的网格划分对于结果的准确性和计算效率至关重要。
过于粗糙的网格会导致计算结果不精确，而过于细密的网格则会增加计算量。
需要根据设计要求和边界条件合理划分网格，尽量在重要的应力集中区域和位移较大的区域细化网格，以获得更准确的结果。
2. 材料本构模型：材料本构模型是用来描述材料力学性质的数学模型，对有限元分析结果的准确性和可靠性有重要影响。
选择合适的本构模型需要考虑材料的性质、应变应力关系和加载条件等因素。
常用的本构模型有弹性模型、塑性模型、粘弹性模型等。
在选择本构模型时，需要根据具体应用场景和加载条件进行合理选择，并进行验证和校准。
3. 边界条件：边界条件是有限元分析中非常重要的一个因素。
它直接影响着模型的应力分布和位移结果。
在设置边界条件时，需要根据实际问题的要求进行准确的设置。
一般包括固支边界、强制位移边界、加载边界等。
在实际应用中，边界条件的设置需要考虑结构的约束和外部加载的作用，并进行合理的假设和简化。
4. 模型验证：模型验证是确保有限元分析结果准确性和可靠性的关键环节。
在进行有限元分析前，可以进行一些简化模型或者理论计算，对部分区域或者特定加载情况进行验证。
验证的方法可以包括理论计算、试验验证、实际工程应用等。
验证的目的是检验有限元模型的准确性和可靠性，进一步提高分析结果的精确性。
5. 结果后处理：有限元分析的结果后处理是对分析结果进行展示和进一步分析的过程。
合适的结果后处理可以帮助工程师们更好地理解分析结果，发现问题和优化设计。
常用的结果后处理方法包括应力和位移的分布图、应变云图、动态变化曲线等。
有限元⽅法的发展及应⽤有限元⽅法的发展及应⽤摘要：有限元法是⼀种⾼效能、常⽤的计算⽅法。
有限元法在早期是以变分原理为基础发展起来的，所以它⼴泛地应⽤于以拉普拉斯⽅程和泊松⽅程所描述的各类物理场中。
⾃从1969年以来，某些学者在流体⼒学中应⽤加权余数法中的迦辽⾦法或最⼩⼆乘法等同样获得了有限元⽅程，因⽽有限元法可应⽤于以任何微分⽅程所描述的各类物理场中，⽽不再要求这类物理场和泛函的极值问题有所联系。
基本思想：由解给定的泊松⽅程化为求解泛函的极值问题。
1有限元法介绍1.1有限元法定义有限元法（FEA，Finite Element Analysis）的基本概念是⽤较简单的问题代替复杂问题后再求解。
它是起源于20世纪50年代末60年代初兴起的应⽤数学、现代⼒学及计算机科学相互渗透、综合利⽤的边缘科学。
有限元法的基本思想是将求解域看成是由许多称为有限元的⼩的互连⼦域组成，对每⼀单元假定⼀个合适的(较简单的）近似解，然后推导求解这个域总的满⾜条件(如结构的平衡条件），从⽽得到问题的解。
这个解不是准确解，⽽是近似解，因为实际问题被较简单的问题所代替。
由于⼤多数实际问题难以得到准确解，⽽有限元不仅计算精度⾼，⽽且能适应各种复杂形状，因⽽成为⾏之有效的⼯程分析⼿段。
有限元法最初应⽤在⼯程科学技术中,⽤于模拟并且解决⼯程⼒学、热学、电磁学等物理问题。
1.2有限元法优缺点有限元⽅法是⽬前解决科学和⼯程问题最有效的数值⽅法，与其它数值⽅法相⽐，它具有适⽤于任意⼏何形状和边界条件、材料和⼏何⾮线性问题、容易编程、成熟的⼤型商⽤软件较多等优点。
（1）概念浅显，容易掌握，可以在不同理论层⾯上建⽴起对有限元法的理解，既可以通过⾮常直观的物理解释来理解，也可以建⽴基于严格的数学理论分析。
（2）有很强的适⽤性，应⽤范围极其⼴泛。
它不仅能成功地处理线性弹性⼒学问题、费均质材料、各向异性材料、⾮线性应⽴-应变关系、⼤变形问题、动⼒学问题已及复杂⾮线性边界条件等问题，⽽且随着其基本理论和⽅法的逐步完善和改进，能成功地⽤来求解如热传导、流体⼒学、电磁场等领域的各类线性、⾮线性问题。
一、基本有限元网格概念1．单元概述ﻫ几何体划分网格之前需要确定单元类型.单元类型的选择应该根据分析类型、形状特征、计算数据特点、精度要求和计算的硬件条件等因素综合考虑。
为适应特殊的分析对象和边界条件,一些问题需要采用多种单元进行组合建模。
ﻫ 2.单元分类选择单元首先需要明确单元的类型，在结构有限元分析中主要有以下一些单元类型：平面应力单元、平面应变单元、轴对称实体单元、空间实体单元、板单元、壳单元、轴对称壳单元、杆单元、梁单元、弹簧单元、间隙单元、质量单元、摩擦单元、刚体单元和约束单元等。
根据不同的分类方法,上述单元可以分成以下不同的形式。
ﻫ3。
按照维度进行单元分类根据单元的维数特征，单元可以分为一维单元、二维单元和三维单元。
ﻫ一维单元的网格为一条直线或者曲线。
直线表示由两个节点确定的线性单元。
曲线代表由两个以上的节点确定的高次单元,或者由具有确定形状的线性单元。
杆单元、梁单元和轴对称壳单元属于一维单元,如图1~图3所示。
ﻫ二维单元的网格是一个平面或者曲面，它没有厚度方向的尺寸.这类单元包括平面单元、轴对称实体单元、板单元、壳单元和复合材料壳单元等,如图4所示。
二维单元的形状通常具有三角形和四边形两种，在使用自动网格剖分时,这类单元要求的几何形状是表面模型或者实体模型的边界面。
采用薄壳单元通常具有相当好的计算效率。
ﻫﻫ三维单元的网格具有空间三个方向的尺寸，其形状具有四面体、五面体和六面体,这类单元包括空间实体单元和厚壳单元，如图5所示．在自动网格划分时，它要求的是几何模型是实体模型(厚壳单元是曲面也可以）。
ﻫ4.按照插值函数进行单元分类根据单元插值函数多项式的最高阶数多少，单元可以分为线性单元、二次单元、三次单元和更高次的单元。
线性单元具有线性形式的插值函数,其网格通常只具有角节点而无边节点，网格边界为直线或者平面.这类单元的优点是节点数量少，在精度要求不高或者结果数据梯度不太大的情况下，采用线性单元可以得到较小的模型规模.但是由于单元位移函数是线性的，单元内的位移呈线性变化,而应力是常数,因此会造成单元间的应力不连续,单元边界上存在着应力突变，如图6所示。
有限元方法的发展及应用有限元方法的发展可以追溯到20世纪50年代，当时数学家、工程师和物理学家开始使用有限元方法来解决结构力学问题。
最早的有限元方法是基于简单的三角形或四边形划分网格，通过近似的方式将连续介质离散化为有限数量的元素。
然后，通过求解一个代数方程组来得到数值解。
这种方法由于计算量小、理论基础牢固而得到了广泛应用。
随着计算机科学的发展，有限元方法得到了更广泛的应用。
计算机技术的进步使得复杂的有限元模型能够被处理，并且计算速度得到了大幅提升。
有限元方法的应用也从最初的结构力学问题扩展到了流体力学、热传导、电磁场、生物医学工程等领域。
有限元方法在工程领域具有很大的应用潜力。
在结构工程中，有限元方法可以用于分析房屋、桥梁和建筑物等结构的强度和刚度。
在汽车工程中，有限元方法可以用于分析汽车的碰撞和安全性能。
在航空航天工程中，有限元方法可以用于分析飞机的气动力学特性和结构强度。
在电子工程和电力工程中，有限元方法可以用于分析电路和传输线的电磁场特性。
有限元方法的应用不仅限于工程领域，还涉及到了其他学科的研究。
在生物医学工程中，有限元方法可以用于模拟人体组织的生物力学行为，如骨骼系统、心脏和血管的应力分布等。
在地球科学中，有限元方法可以用于分析地下水流动、地震波传播和岩土工程等问题。
在物理学中，有限元方法可以用于分析电磁场、热传导和量子力学等问题。
总之，有限元方法的发展及其应用已经取得了巨大的成功。
它在工程、力学、物理和地球科学等领域中得到了广泛应用，并为实际工程问题的解决提供了有效的数值方法。
然而，有限元方法的进一步发展仍面临着一些挑战，需要继续改进算法和技术，以满足更加复杂和多样化的工程问题的需求。
有限元网格划分及发展趋势摘要:总结近十年有限元网格划分技术发展状况。
首先,研究和分析有限元网格划分的基本原则；其次,对当前典型网格划分方法进行科学地分类,结合实例,系统地分析各种网格划分方法的机理、特点及其适用范围,如映射法、基于栅格法、节点连元法、拓扑分解法、几何分解法和扫描法等；再次,阐述当前网格划分的研究热点,综述六面体网格和曲面网格划分技术；最后,展望有限元网格划分的发展趋势。
关键词:有限元网格划分；映射法；基于栅格法；节点连元法;拓扑分解法；几何分解法；扫描法；六面体网格1 引言作为有限元走向工程应用枢纽的有限元网格划分,是有限元法的一个非常重要的研究领域,经历了40多年的发展历程。
有限元网格划分算法研究中的某些难点问题始终未能得到真正意义上的解决,它们的解决对工程问题具有重要的现实价值和理论意义。
有限元分析的基本过程可分为三个阶段:有限元模型的建立(即前处理)、有限元解算、结果处理和评定(即后处理)。
根据经验,有限元分析各阶段所用的时间为】【1:40%-45%用于模型的前处理,50%-55%用于后处理,而分析计算只占5%左右;更有文献】【2指出有限元建模占有限元分析一半以上的工作量,甚至高达80%。
因此,有限元分析的前后处理一直都是有限元分析的瓶颈问题,严重地阻碍着有限元分析技术的应用和发展。
许多学者对有限元网格生成方法近30年的研究进行了概括和总结】【4。
近年来,【3,对某些重要分支领域的研究进展方面也做出了贡献】有限元网格生成方法研究有两个显著特点:(1)经历了一个进化过程,一些方法的研究与应用出现停滞,而另外一些方法在不断地深入、完善和发展,成为适应性强、应用范围广泛的通用方法;(2)领域和主题在不断扩展和深入,研究重点由二维平面问题转移到三维曲面和三维实体问题,从三角形、四面体网格自动生成转移到四边形、六面体网格自动生成。
2 有限元网格划分的基本原则有限元方法的基本思想是将结构离散化,即对连续体进行离散化,利用简化几何单元来近似逼近连续体,然后根据变形协调条件综合求解。
所以有限元网格的划分一方面要考虑对各物体几何形状的准确描述,另一方面也要考虑变形梯度的准确描述。
为正确、合理地建立有限元模型,这里介绍划分网格时应考虑的一些基本原则。
2.1 网格数量网格数量直接影响计算精度和计算时耗,网格数量增加会提高计算精度,但同时计算时耗也会增加。
当网格数量较少时增加网格，计算精度可明显提高,但计算时耗不会有明显增加;当网格数量增加到一定程度后,再继续增加网格时精度提高就很小,而计算时耗却大幅度增加。
所以在确定网格数量时应权衡这两个因素综合考虑。
2.2 网格密度为了适应应力等计算数据的分布特点,在结构不同部位需要采用大小不同的网格。
在孔的附近有集中应力,因此网格需要加密;周边应力梯度相对较小,网格划分较稀。
由此反映了疏密不同的网格划分原则:在计算数据变化梯度较大的部位,为了较好地反映数据变化规律,需要采用比较密集的网格;而在计算数据变化梯度较小的部位,为减小模型规模,网格则应相对稀疏。
2.3 单元阶次单元阶次与有限元的计算精度有着密切的关联,单元一般具有线性、二次和三次等形式,其中二次和三次形式的单元称为高阶单元。
高阶单元的曲线或曲面边界能够更好地逼近结构的曲线和曲面边界,且高次插值函数可更高精度地逼近复杂场函数,所以增加单元阶次可提高计算精度。
但增加单元阶次的同时网格的节点数也会随之增加,在网格数量相同的情况下由高阶单元组成的模型规模相对较大,因此在使用时应权衡考虑计算精度和时耗。
2.4 单元形状网格单元形状的好坏对计算精度有着很大的影响,单元形状太差的网格甚至会中止计算。
在网格划分时应保证合理的单元形状,即使只有一个单元形状很差或畸形(畸形网格通常有单元节点交叉编号、节点重合、内角大于180°和单元面积为零等类型)时,也可能给计算结果带来很大的误差,甚至使得计算无法进行下去。
单元形状评价一般有以下几个指标:(1)单元的边长比、面积比或体积比以正三角形、正四面体、正六面体为参考基准,理想单元的边长比为1,线性单元可接受的边长比小于3,二次单元小于10。
(2)扭曲度:单元面内的扭转和面外的翘曲程度。
(3)节点编号:节点编号对于求解过程中总刚矩阵的带宽和波前因数有较大的影响,从而影响计算时耗和存储容量的大小。
因此合理的节点编号有利于刚度矩阵对称、带状分布等求解效率,从而提高计算速度。
2.5 单元协调性单元协调是指单元上的力和力矩能够通过节点传递给相邻单元。
为保证单元协调,必须满足的条件是:(1)一个单元的节点必须同时也是相邻点,而不应是内点或边界点。
(2)相邻单元的共有节点具有相同的自由度性质，例如SHELL和SOLID两种单元,前者为六自由度,而后者为平动三自由度,所以SOLID只能向SHELL传递平动位移,不能传递转动位移。
另外,有相同自由度的单元网格也并非一定协调。
例如,BEAM3和SHELL41都具有三个自由度,但BEAM3具有三个平动自由度,SHELL41有两个平动自由度和一个转动自由度,所以SHELL41只能向BEAM3传递平动位移,却不能传递转动位移。
3 网格生成通用方法回顾有限元网格划分方法难以准确分类,分类方法有很多,可以按产生的单元类型、生成单元的维数、自动化程度等进行分类。
而且目前的网格生成法分类中,对某一具体的网格生成法叫法不一致,如有人将拓扑分解法叫做结构分解法,将节点连元法叫做填充法。
本节仅对众多网格划分方法中的典型方法的算法思想进行回顾和总结。
3.1 映射法映射法的基本思想是实际图形与标准图形的双向映射,具体有三个步骤:(1)根据形体边界的参数方程,利用适当的映射函数,将待划分的物理域映射到参数空间,形成规则参数域;(2)对参数域进行网格划分;(3)将参数空间内单元的网格反向映射到欧氏空间,从而生成实际的网格。
这种网格控制机理有以下几个缺点:(1)映射法不是完全面向几何特征的,所以很难完成自动化,尤其是对于3D区域;(2)网格局部控制能力差;(3)各映射块之间的网格密度相互影响程度很大,改变某一映射块的网格密度,其它映射块的网格都要做相应的调整;(4)对于形状较为复杂的形体适应性差,要求事先根据所要产生的网格类型将目标域分割成一系列可映射的子区域。
子域分解繁琐,所需人工交互多,难以实现网格自动的生成。
3.2 基于栅格法基于栅格法(Grid-based Approach)也叫空间分解法,最早由Thacker等】【5提出。
该算法的基本流程:首先用一组不相交的栅格(cells)覆盖在物体之上,既可在栅格的规则点处布置节点,也可在栅格单元中随机布置节点;再对栅格和物体进行相交检测,保留完全或部分落在目标区域之内的栅格,删除完全落在目标区域之外的栅格;然后对与物体边界相交的栅格进行调整、剪裁、再分解等操作,使其更准确地逼近目标区域;最后对内部栅格和边界栅格进行栅格级的网格剖分,进而得到整个目标区域的有限元网格。
3.3 节点连元法节点连元法一般分为两步:(1)在物体的边界和有效区域内按照网格密度的要求均匀布点;(2)根据一定的准则将这些节点连接成三角形或四面体网格。
最早提出的节点生成法是随机布点法,该方法不能保证布点的均匀性,同时为了确定可用节点,要进行点距检查和多次试探,算法效率低。
近年出现的长方形网格直接布点法、等距水平扫描线直接布点法等方法在一定程度上提高了节点连元法的计算效率,但Delaunay三角剖分方法(简称DT法)的出现才真正使得节点连元法成为目前主流的全自动网格划分方法之一。
3.4 拓扑分解法拓扑分解法(Topology Decomposition Approach),它首先是由英国剑桥大学的Wordenwaber提出来的。
拓扑分解法是从形体的拓扑因素着手进行分割,而不过问元素的具体形状。
首先假设网格顶点全部由目标边界顶点组成,那么可以用一种三角化算法将目标用尽量少的三角形完全分割覆盖。
该方法生成的网格一般相当粗糙,必须与其它方法相结合,通过网格加密等过程,才能生成合适的网格。
拓扑分解法的原理简单,易于处理,但它只从拓扑关系入手,不考虑几何因素,因此难以获得良好形状的单元网格,也不适用于含有曲面的三维形体。
该方法现在已普遍使用于目标初始三角化,用来实现从实体表述到初始三角化表述的自动化转换。
3.5 几何分解法几何分解法(Geometry Decomposition Approach)最大的特点是节点和单元同步生成。
该方法较多地考虑了待分域的几何特征,确保生成质量较好的网格单元。
几何分解法可以进一步分成区域递归细分法、单元移去方法和子域分离方法。
(1)区域递归细分法(Recursive Subdivision Method)由Bykat A 提出[45,46],他将实体分割成多个凸区域,并在分割线上按网格密度要求插入一组新的节点,分割的最终结果为三角形网格。
(2)单元移去方法(Element Removing Method)是从实体的边界向里一次移去一个单元网格,迭代执行下去,推进波前法就是该方法的成功代表。
推进波前法(Advancing Front Technique,AFT)最初是由Lo提出,AFT方法的基本流程是:首先离散待剖分域的边界,形成首尾相连的线段集合或拓扑相容的三角形面片集合,这种离散后的域边界称为前沿;再从前沿开始,依次插入一个节点,并连接生成一个新的单元,更新前沿,将前沿向待剖分域的内部推进。
(3)子域分离法是将问题区域一次移去一块,直到剩下最后一块。
该方法一般可以分为两类:一类是每次分离出一个单元;另一类是每次分离出一个拓扑简单的子域块,然后利用其它划分算法细化子域块。
这种方法一般以层次的形式利用大量的移去操作,并考虑当前的移去操作对以后的移去子块选择的影响。
该方法的每一步移去操作所需的计算量很大,所以该算法在生成细化网格方面计算效率并不高,但是可以利用该算法,将整个区域分解为一个个易于网格化的子区域,其总的效率将会大大提高。
3.6 扫描法扫描法是将离散化的基本单元形体进行旋转、扫描、拉伸等操作,获得高维网格的一种方法。
这种方法难度较低,容易实现,在当今大多数商用CAD软件和有限元前置处理软件中均有这种功能。
但是这种方法只适合于形状简单的三维物体,且主要靠人机交互来实现,自动化程度低。
4 研究热点近年来有限元分析在各种工程领域中得到了广泛的应用,网格划分技术的理论基础已日趋成熟,特别是在二维空间以及三维空间的四面体网格划分方面。
近几年有限元网格划分的研究领域已由二维平面问题转移到三维实体,研究重点已经由三角形(四面体)网格转变为四边形(六面体)网格,注重网格的全自动生成、网格自适应等研究,本章主要总结当前网格划分领域的研究热点。
