有限元网格划分方法与基本原理

合集下载

第07讲-有限元网格划分的基本原则及技巧

第07讲-有限元网格划分的基本原则及技巧

7-6
网格疏密
• • 网格疏密是指在结构不同部位采用大小不同的网格,这是为了适应计算数据的分 布特点。 在计算数据变化梯度较大的部位(如应力集中处、几何形状、材料、厚度变化的 位置),为了较好地反映数据变化规律,需要采用比较密集的网格。而在计算数 据变化梯度较小的部位,为减小模型规模,则应划分相对稀疏的网格。这样,整 个结构便表现出疏密不同的网格划分形式。—— 网格数量应增加在结构的关键 部位,在次要部位增加网格是不必要的,也是不经济的。 边界上最好要在8个单元以上,至少不少于4个; 分析结果完成后,需要检查以下各项,误差较大的位置要进行细分: 单元应力的连续性,比较相邻单元应力值的差值; 应力偏差:结点上的单元结点应力和结点平均应力的差值的较大值; 当以上差值与其中的最大应力的比值较大时,该位置的网格需要细分。
精度 计算时间 精确解 1 2 O
7-4


P
网格数量
网格数量(续)
在决定网格数量时应考虑分析数据的类型。 实体单元:
• •
1、在静力分析时,如果仅仅是计算结构的变形,网格数量可以少一些。如 果需要计算应力,则在精度要求相同的情况下应取相对较多的网格。 2、在响应计算中,计算应力响应所取的网格数应比计算位移响应多。 3、在计算结构固有动力特性时,若仅仅是计算少数低阶模态,可以选择较 少的网格,如果计算的模态阶次较高,则应选择较多的网格。
左图中(a)、(b)改 变了结构质量的对称分 布,应避免。 (c)是 比较理想的结果。
(a)
7-8
(b)
(c)
单元的形状及评价
• 形状比(长边与短边距离之比) 一般实体单元的长宽比越大,分析误差也越大。 对于板壳单元,评价应力为主时不宜超过1:3,评价位移为主时不宜超过1:5; 对于块体单元,评价应力为主时不宜超过1:2,评价位移为主时不宜超过1:3; 在应力分布几乎没有变化的区域里使用的单元,适当放大也没问题。 倾角(表示单元偏离直角四边形的程度(Angular Deviation)) 四边形的内倾角最好是在45度~135度之间,不要超过15度~165度。 锥度(限于四边形) 用几何偏离(Geometric Deviation)表示四边形单元的变形程度。

Deform网格划分原则及方法

Deform网格划分原则及方法

[原]Deform网格划分原则及方法2009-04-04 23:48引言:划分网格是建立有限元模型的一个重要环节,它要求考虑的问题较多,需要的工作量较大,所划分的网格形式对计算精度和计算规模将产生直接影响。

为建立正确、合理的有限元模型,这里介绍网格划分时的一些基本原则及方法。

关键词:Deform 网格局部细化一、网格划分的原则1 网格数量网格数量的多少将影响计算结果的精度和计算规模的大小。

一般来讲,网格数量增加,计算精度会有所提高,但同时计算规模也会增加,所以在确定网格数量时应权衡两个因数综合考虑。

图1中的曲线1表示结构中的位移随网格数量收敛的一般曲线,曲线2代表计算时间随网格数量的变化。

可以看出,网格较少时增加网格数量可以使计算精度明显提高,而计算时间不会有大的增加。

当网格数量增加到一定程度后,再继续增加网格时精度提高甚微,而计算时间却有大幅度增加。

所以应注意增加网格的经济性。

实际应用时可以比较两种网格划分的计算结果,如果两次计算结果相差较大,可以继续增加网格,相反则停止计算。

图1 位移精度和计算时间随网格数量的变化在决定网格数量时应考虑分析数据的类型。

在静力分析时,如果仅仅是计算结构的变形,网格数量可以少一些。

如果需要计算应力,则在精度要求相同的情况下应取相对较多的网格。

在热分析中,结构内部的温度梯度不大,不需要大量的内部单元,这时可划分较少的网格。

2 网格疏密网格疏密是指在结构不同部位采用大小不同的网格,这是为了适应计算数据的分布特点。

在计算数据变化梯度较大的部位(如应力集中处),为了较好地反映数据变化规律,需要采用比较密集的网格。

而在计算数据变化梯度较小的部位,为减小模型规模,则应划分相对稀疏的网格。

这样,整个结构便表现出疏密不同的网格划分形式。

图2是中心带圆孔方板的四分之一模型,其网格反映了疏密不同的划分原则。

小圆孔附近存在应力集中,采用了比较密的网格。

板的四周应力梯度较小,网格分得较稀。

有限元计算原理

有限元计算原理

有限元计算原理
有限元计算原理是一种工程分析的方法,用于求解各种结构及连续体的力学问题。

其基本思想是将结构或连续体分割成有限数量的小单元,然后通过对这些小单元进行计算,再将其组合起来求解整体问题。

这种方法可以将结构或连续体的力学行为分析得非常精确,可以获得结构的应力应变分布、位移分布等信息。

有限元计算的原理可以概括为以下几个步骤:
1. 网格划分:将结构或连续体划分成许多小单元,即有限元,这些小单元通过节点连接起来构成整个结构。

2. 求解力学方程:根据结构或连续体的几何形状和物理特性,建立相应的力学方程组。

通常采用弹性力学理论来描述结构或连续体的力学行为。

3. 边界条件的处理:给定结构或连续体的边界条件,如固支、约束力等,在有限元网格中对应的节点上施加相应的约束。

4. 单元刚度矩阵的组装:通过计算每个小单元的刚度矩阵,将其组装成整个结构或连续体的整体刚度矩阵。

5. 单元荷载向量的组装:根据给定的荷载条件,在每个小单元上计算相应的荷载向量,将其组装成整个结构或连续体的荷载向量。

6. 求解位移和应力:根据组装好的整体刚度矩阵和荷载向量,通过求解线性方程组,得到结构或连续体中每个节点的位移和应力。

7. 后处理:根据求解得到的位移和应力,可以计算出结构或连续体的各种物理量,比如应变、应力、变形等。

通过这种有限元计算的方法,可以对各种复杂的结构或连续体进行力学分析和优化设计。

机械零件有限元分析-5-第四讲-网格划

机械零件有限元分析-5-第四讲-网格划

THANKS
感谢观看
理现象。
均匀性
网格的分布应尽量均匀,以提 高计算精度和稳定性。
局部细化
对于关键区域或需要更高精度 的地方,应进行局部网格细化

边界条件处理
在边界区域,应根据实际情况 处理网格,以避免出现奇异性
和不合理的解。
03
网格划分的方法和技术
结构化网格划分
01
02
03
结构化网格
按照一定的规则和顺序对 有限元模型进行网格划分, 每个网格单元具有相同或 相似的形状和尺寸。
详细描述
对于形状不规则、结构复杂的机械零件,网格划分变得困难,需要采用特殊的有 限元网格划分方法,如自适应网格、非结构化网格等。
实例三:多物理场耦合的网格划分
总结词
多物理场、耦合、复杂度增加
详细描述
对于涉及多个物理场耦合的机械系统,如热-力耦合、流-固耦合等,网格划分变得更加复杂。需要采用多物理场 耦合的有限元网格划分方法,如分区耦合、全局耦合等。
网格划分的重要性和意义
网格划分是有限元分析的关键 环节,它决定了模型的离散精 度和计算规模。
合适的网格划分能够提高计算 精度,降低模型的自由度,从 而减少计算时间和资源消耗。
不合理的网格划分可能导致计 算精度降低,甚至出现数值不 稳定或计算失败的情况。
02
网格划分的基本概念
网格划分的定义
网格划分是将连续的物理模型离散化 为有限个小的单元,每个单元称为网 格或节点。
自适应移动节点
03
根据计算结果动态移动网格节点,以保持网格质量。
05
实例分析
实例一:简单零件的网格划分
总结词
规则、简单、容易划分
详细描述

有限元网格剖分与网格质量判定指标

有限元网格剖分与网格质量判定指标

有限元网格剖分与网格质量判定指标有限元网格剖分与网格质量判定指标一、引言有限元法是一种常用的数值分析方法,广泛应用于工程、力学等领域。

在有限元方法中,对于复杂的几何体,需要将其分割成多个简单的几何单元,称为有限元。

而有限元的形状和尺寸对计算结果的精度和稳定性有重要影响。

因此,有限元网格剖分和网格质量判定指标的选择和优化是提高有限元方法计算精度和效率的关键。

二、有限元网格剖分的基本原则和方法有限元网格剖分的基本原则是要确保网格足够细密,以捕捉几何体的细节和特征。

一般来说,有限元网格剖分可以分为以下几个步骤:1. 几何体建模:根据实际问题建立几何体模型,可以使用CAD软件进行建模。

2. 离散化:将几何体分割成简单的几何单元,如三角形、四边形或六面体等。

3. 网格生成:根据几何单元的尺寸和形状要求生成网格。

一般可采用三角形剖分算法或四边形剖分算法进行网格生成。

4. 网格平滑:对生成的网格进行平滑处理,以提高网格的质量。

三、网格质量判定指标网格质量判定指标是用来评价和衡量网格质量好坏的指标。

一个好的网格是指网格单元形状较正、网格单元之间大小相近、网格单元的边界规则等。

常用的网格质量判定指标包括:1. 网格单元形状度:用于评价网格单元的形状正交性和变形。

常用的形状度指标有内角度、调和平均内角度和狄利克雷三角形剖分等。

2. 网格单元尺寸误差:用于评价网格单元尺寸与理想尺寸之间的差异。

常用的尺寸误差指标有网格单元长度标准差、最大和最小网格单元尺寸比等。

3. 网格单元的四边形度:用于评价四边形网格的形状规则性。

常用的四边形度指标有圆度、直角度和Skewness等。

四、网格质量优化方法为了改善有限元网格质量,可以采用以下方法:1. 网格加密:通过将大尺寸网格单元划分为小尺寸网格单元,提高网格的细密度。

2. 网格平滑:通过对矩阵约束或拉普拉斯平滑等方法对网格进行平滑处理,改善网格单元的形状。

3. 网格优化:通过对网格单元的拓扑结构和形状进行优化,提高网格的质量。

电磁场计算中的有限元方法教程

电磁场计算中的有限元方法教程

电磁场计算中的有限元方法教程引言电磁场计算是电磁学领域中重要的研究内容之一,广泛应用于电气工程、通信工程、电子技术等领域。

而有限元方法(Finite Element Method,简称FEM)是一种常用的数值计算技术,可以解决电磁场计算中的复杂问题。

本文将介绍有限元方法在电磁场计算中的基本原理、步骤和应用。

一、有限元方法简介有限元方法是一种通过将待求解区域划分成有限数量的小单元,利用单元上的近似函数构造整个区域上的解的数值计算方法。

有限元方法的基本思想是在每个小单元内近似解以建立一个代数方程组,通过将这些方程组联立得到整个区域上的解。

有限元方法具有处理复杂几何形状、边界条件变化和非线性问题的优势,因此被广泛应用于工程和科学计算中。

二、电磁场方程建立在电磁场计算中,关键是建立合适的电磁场方程。

常见的电磁场方程包括静电场方程、恒定磁场方程、麦克斯韦方程等。

根据具体情况选择适用的方程,并根据材料的性质和边界条件确定相应的方程形式。

三、有限元网格划分有限元方法需要将计算区域划分为有限数量的小单元。

在电磁场计算中,通常采用三角形或四边形单元来进行划分,这取决于计算区域的几何形状和分辨率要求。

划分过程需要考虑电场变化的特点和计算精度的需求,合理划分网格对精确计算电磁场起着重要的作用。

四、有限元方程的建立有限元网格划分完成后,需要建立相应的有限元方程组。

以求解静电场问题为例,我们可以利用能量最小原理、偏微分方程等方法建立有限元方程组。

有限元方程组的建立需要考虑电场的连续性、边界条件和材料特性等。

五、有限元方程求解有限元方程组的求解是求解电磁场分布的核心任务。

根据具体的方程形式和计算区域的几何形状,可以采用直接法、迭代法、近似法等方法来求解方程。

在电磁场计算中,常用的求解算法包括高斯消元法、迭代法、有限元法和有限差分法等。

六、计算结果的后处理在得到有限元方法计算的电磁场分布结果后,需要进行相应的后处理,进行数据分析和可视化。

有限元结合格子boltzmann方法

有限元结合格子boltzmann方法

有限元结合格子boltzmann方法随着计算机技术的飞速发展,数值模拟方法在工程领域中的应用越来越广泛。

有限元法(FEM)和格子Boltzmann方法(LBM)作为两种常见的数值方法,各自具有独特的优势。

将这两种方法相结合,可以充分发挥它们在计算流体力学、材料科学等领域的潜力。

本文将简要介绍有限元结合格子Boltzmann方法的基本原理及其在工程中的应用。

一、有限元法与格子Boltzmann方法简介1.有限元法(FEM)有限元法是一种将连续域问题转化为离散问题求解的数值方法。

它通过将复杂的几何形状划分成简单的单元(如三角形或四边形),在每个单元内采用插值函数近似求解偏微分方程,从而实现整个域上的问题求解。

2.格子Boltzmann方法(LBM)格子Boltzmann方法是一种基于微观粒子的动力学行为的宏观现象模拟方法。

它通过离散化的Boltzmann方程,在格子网络上模拟粒子的碰撞和传播过程,从而得到宏观物理量(如速度、密度等)。

二、有限元结合格子Boltzmann方法的基本原理有限元结合格子Boltzmann方法的主要思想是将FEM的高精度与LBM 的微观模拟相结合,以解决复杂的流体力学问题。

具体步骤如下:1.划分网格:在计算域内同时采用有限元和格子Boltzmann方法进行网格划分,其中有限元网格主要用于求解宏观物理量,而格子Boltzmann网格则用于模拟微观粒子的运动。

2.确定边界条件:根据实际问题,为有限元和格子Boltzmann方法设置相应的边界条件。

3.求解宏观物理量:利用有限元法求解宏观物理量,如速度、压力等。

4.更新微观粒子分布函数:在格子Boltzmann网格上,根据微观粒子的碰撞和传播过程,更新粒子的分布函数。

5.反向映射:将格子Boltzmann方法得到的微观粒子信息映射到有限元网格上,更新宏观物理量。

6.迭代求解:重复步骤3-5,直至满足收敛条件。

三、有限元结合格子Boltzmann方法在工程中的应用有限元结合格子Boltzmann方法在工程领域具有广泛的应用前景,以下列举几个典型应用:1.计算流体力学:结合FEM的高精度和LBM的微观模拟,可以更准确地预测复杂流场中的流动现象。

三维有限元模型

三维有限元模型

三维有限元模型一、引言三维有限元模型是一种数学计算方法,用于分析和解决复杂的结构问题。

它可以将实际结构转化为由许多小单元组成的离散化模型,并通过数学方程求解每个单元的应力、应变等物理量,最终得出整个结构的响应。

本文将介绍三维有限元模型的基本原理、建模方法和求解过程。

二、三维有限元模型基本原理1. 有限元法基本思想有限元法是一种数值计算方法,它将一个连续的物理问题转化为由许多小单元组成的离散化问题,在每个小单元上建立数学模型,并通过求解代数方程组来得到整个系统的响应。

在三维有限元模型中,通常采用四面体或六面体等简单形状的单元进行离散化。

2. 三维有限元模型建立过程(1)几何建模:根据实际结构进行几何建模,包括确定结构尺寸、形状等。

(2)网格划分:将几何模型划分为许多小单元,并确定每个单元节点坐标。

(3)材料参数:根据实际材料性质确定每个单元的杨氏模量、泊松比等物理参数。

(4)载荷边界条件:根据实际工况确定结构所受载荷和边界条件。

(5)约束边界条件:根据实际结构确定约束边界条件,如支座、铰链等。

(6)求解:将以上信息输入计算机中,通过数学方法求解每个单元的应力、应变等物理量,并得出整个结构的响应。

三、三维有限元模型建模方法1. 网格划分方法三维有限元模型的网格划分可以采用手动或自动方式进行。

手动划分需要经验丰富的工程师进行,通常用于简单结构;自动划分则是利用计算机软件进行,可以快速生成复杂结构的网格。

2. 材料模型在三维有限元模型中,通常采用线性弹性模型来描述材料行为。

这种模型假设材料是各向同性的,并且满足胡克定律。

如果需要考虑非线性效应,则需要采用非线性材料模型。

3. 载荷和边界条件在三维有限元模型中,载荷和边界条件是建模的重要组成部分。

载荷可以是静载荷、动载荷或温度载荷等,边界条件可以是支座、铰链等。

四、三维有限元模型求解过程1. 单元刚度矩阵单元刚度矩阵是计算每个单元应力和应变的关键。

它由每个单元的杨氏模量、泊松比和几何信息确定。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

结构有限元分析中的网格划分技术及其应用实例结构有限元分析中的网格划分是否直接关系到解算的效果。

本文简述了网格划分应用的基本理论,并以空间自由曲面覆盖件和大型整体网络钢筋壳体产品的有限元分析中的网格划分为实例对象,详细讲述了空间自由和三维实体的网格划分基本理论及其在工程中的实际应用,非常具有现实意义和借鉴价值。

一、前言有限元网格划分是进行有限元数值模拟分析至关重要的一步,它直接影响着后续数值计算分析结果的精确性。

网格划分涉及单元的形状及其拓扑类型、单元类型、网格生成器的选择、网格的密度、单元的编号以及几何体素。

从几何表达上讲,梁和杆是相同的,从物理和数值求解上讲则是有区别的。

同理,平面应力和平面应变情况设计的单元求解方程也不相同。

在有限元数值求解中,单元的等效节点力、刚度矩阵、质量矩阵等均用数值积分生成,连续体单元以及壳、板、梁单元的面内均采用高斯(Gauss)积分,而壳、板、梁单元的厚度方向采用辛普生(Simpson)积分。

辛普生积分点的间隔是一定的,沿厚度分成奇数积分点。

由于不同单元的刚度矩阵不同,采用数值积分的求解方式不同,因此实际应用中,一定要采用合理的单元来模拟求解。

CAD软件中流行的实体建模包括基于特征的参数化建模和空间自由曲面混合造型两种方法。

Pro/E和SoildWorks是特征参数化造型的代表,而 CATIA与Unigraphics等则将特征参数化和空间自由曲面混合造型有机的结合起来。

现有CAD软件对表面形态的表示法已经大大超过了CAE软件,因此,在将CAD实体模型导入CAE软件的过程中,必须将CAD模型中其他表示法的表面形态转换到CAE软件的表示法上,转换精度的高低取决于接口程序的好坏。

在转换过程中,程序需要解决好几何图形(曲线与曲面的空间位置)和拓扑关系(各图形数据的逻辑关系)两个关键问题。

其中几何图形的传递相对容易实现,而图形间的拓扑关系容易出现传递失败的情况。

数据传递面临的一个重大挑战是,将导入CAE程序的CAD模型改造成适合有限元分析的网格模型。

在很多情况下,导入CAE程序的模型可能包含许多设计细节,如细小的孔、狭窄的槽,甚至是建模过程中形成的小曲面等。

这些细节往往不是基于结构的考虑,保留这些细节,单元数量势必增加,甚至会掩盖问题的主要矛盾,对分析结果造成负面影响。

CAD模型的“完整性”问题是困扰网格剖分的障碍之一。

对于同一接口程序,数据传递的品质取决于CAD模型的精度。

部分CAD模型对制造检测来说具备足够的精度,但对有限元网格剖分来说却不能满足要求。

值得庆幸的是,这种问题通常可通过CAD软件的“完整性检查”来修正。

改造模型可取的办法是回到CAD系统中按照分析的要求修改模型。

一方面检查模型的完整性,另一方面剔除对分析无用的细节特征。

但在很多情况下,这种“回归”很难实现,模型的改造只有依靠 CAE软件自身。

CAE中最直接的办法是依靠软件具有的“重构”功能,即剔除细部特征、缝补面和将小面“融入”大曲面等。

有些专用接口在模型传递过程中甚至允许自动完成这种工作,并且通过网格剖分器检验模型的“完整性”,如发现“完整性”不能满足要求,接口程序可自动进行“完整性”修复。

当几何模型距 CAE分析的要求相差太大时,还可利用CAE程序的造型功能修正几何模型。

“布尔运算”是切除细节和修理非完整特征的有效工具之一。

目前数据传递一般可通过专用数据接口,CAE程序可与CAD程序“交流”后生成与CAE 程序兼容的数据格式。

另一种方式是通过标准图形格式如IGES、 SAT和ParaSolid传递。

现有的CAD平台与通用有限元平台一般通过IGES、STL、Step、Parasolid等格式来数据交换,早期 IGES接口应用比较广泛,但由于该标准本身的不严格性,导致多数复杂模型的传递以失败告终,如图1所示为某汽车覆盖件在UGII中以IGES格式输出时产生的信息,可以看出其包含大量有限元分析不必要的几何信息。

而SAT与ParaSolid标准较为严格,被多数CAD程序采用。

由于典型通用有限元软件(如MSC.PATRAN、MSC.MARC、ANSYS、ABAQUS、AD INA等)的建模功能都不是很强,尤其是在面对包含复杂空间曲面的产品结构时表现出明显的不足,同时不利于建立后续的单元网格划分模型。

因此,利用现有CAD平台(如CATIA、UGII、PRO/E)完成网格划分工作,或借助专业网格划分软件HyperMesh、AIEnviroment等来完成任务是比较好的方法。

下面分别以包含大量空间自由曲面的汽车覆盖件产品和宇航业中常用的大型整体网格筋壳体为对象,简述有限元网格划分的基本原理方法和应用。

图1 IGES文件输出的图素信息分页有限元网格划分方法与基本原理1.有限元网格划分的指导思想有限元网格划分的指导思想是首先进行总体模型规划,包括物理模型的构造、单元类型的选择、网格密度的确定等多方面的内容。

在网格划分和初步求解时,做到先简单后复杂,先粗后精,2D单元和3D单元合理搭配使用。

为提高求解的效率要充分利用重复与对称等特征,由于工程结构一般具有重复对称或轴对称、镜象对称等特点,采用子结构或对称模型可以提高求解的效率和精度。

利用轴对称或子结构时要注意场合,如在进行模态分析、屈曲分析整体求解时,则应采用整体模型,同时选择合理的起点并设置合理的坐标系,可以提高求解的精度和效率,例如,轴对称场合多采用柱坐标系。

有限元分析的精度和效率与单元的密度和几何形状有着密切的关系,按照相应的误差准则和网格疏密程度,避免网格的畸形。

在网格重划分过程中常采用曲率控制、单元尺寸与数量控制、穿透控制等控制准则。

在选用单元时要注意剪力自锁、沙漏和网格扭曲、不可压缩材料的体积自锁等问题。

典型有限元软件平台都提供网格映射划分和自由适应划分的策略。

映射划分(Mapped/ IsoMesh)用于曲线、曲面、实体的网格划分方法,可使用三角形、四边形、四面体、五面体和六面体,通过指定单元边长、网格数量等参数对网格进行严格控制,映射划分只用于规则的几何图素,对于裁剪曲面或者空间自由曲面等复杂几何体则难以控制。

自由网格划分(F ree/Paver)用于空间自由曲面和复杂实体,采用三角形、四边形、四面体进行划分,采用网格数量、边长及曲率来控制网格的质量。

例如,在MSC.MARC中,其转换(Convert)用法是几何模型转换为网格模型,点转换为节点,曲线转换为线单元,面转换为三角形、四边形等。

网格自动划分(AutoMesh)则是在任意曲面上生成三角形或者四边形,对任意几何体生成四面体或者六面体。

网格重划分(Remesh)是在每一步计算过程中,检查各单元法向来判定各区域的曲率变化情况,在曲率较大变形剧烈的区域单元,进行网格加密重新划分,如此循环直到满足网格单元的曲率要求为止。

网格重划分的思想是通过网格加密的方法来提高分析的精度和效率。

网格自适应划分(Adaptive Refinement)的思想是在计算步中,升高不满足分析条件的低阶单元的阶次来提高分析的精度和效率,应用比较广泛。

自适应网格划分必须采用适当的单元,在保证单元阶次的基础上,原本已形成的单元刚度矩阵等特性保持不变,才能同时提高精度和效率。

阶谱单元(Hierachical Element)充分发挥了自适应网格划分的优点,在计算中通过不断增加初始单元的边上的节点数,从而使单元插值函数的阶次在前一阶的基础上不断增加,通过引入新增节点的插值函数来提高求解的精度和效率。

例如,三节点三角形单元升为六节点三角形单元,四节点四边形单元升阶为8节点四边形单元,四节点四面体单元升阶为8节点、10节点、20节点四面体。

2.有限元网格划分的基本方法有限元网格划分方法有两种,对于简单的结构多采用直接建立单元模型的网格直接生成法,当对象比较复杂时,多通过几何自动生成法来完成,即在几何元素描述的物理基础上自动离散成有限单元。

有限元单元可以按几何维数划分为一维、二维和三维单元,而在实际应用中采用拓扑结构单元,包括常用的质量单元、弹簧元、杆与梁管单元、平面三角形单元、平面四边形单元、膜单元、等参单元、壳单元和三维实体单元。

有限元网格划分,对于二维平面、三维曲面和三维实体网格有以下几种划分方法:(1)覆盖法:基于四边形的网格划分,要求网格划分的平面或曲面必须是完整裁减曲面,该曲面边界必须是裁减曲线;(2)前沿法:通过把曲面等参变换到二维空间进行网格划分,然后映射到三维空间曲面上,把曲面划分成完全的四边形单元或三角形单元;(3)Delaunay三角形法:主要用于由至少一条封闭曲线所围成的单连通域或多连通域内生成三角形单元,趋向于等边三角形。

充分考虑了几何形状中细微的几何特征,并在微小特征处划分成较细的单元,在不需要密网格处,采用稀疏单元网格。

(4)转换扩展法:针对曲面几何形状比较规则的几何区域进行网格划分,其网格生成速度快,网格质量高。

由节点扩展为线单元,从线单元生成平面二维单元,从二维单元生成三维单元。

它不仅仅用于三维网格的生成,同时可进行一维、二维网格和几何体的生成,包括移动、镜像、拉伸、旋转、扫描三维实体的扩展方式、扩展系数和扩展方向。

3.网格质量的评估单元的质量和数量对求解结果和求解过程影响较大,如果结构单元全部由等边三角形、正方形、正四面体、立方六面体等单元构成,则求解精度可接近实际值,但由于这种理想情况在实际工程结构中很难做到。

因此根据模型的不同特征,设计不同形状种类的网格,有助于改善网格的质量和求解精度。

单元质量评价一般可采用以下几个指标:(1)单元的边长比、面积比或体积比以正三角形、正四面体、正六面体为参考基准。

理想单元的边长比为1,可接受单元的边长比的范围线性单元长宽比小于3,二次单元小于10。

对于同形态的单元,线性单元对边长比的敏感性较高阶单元高,非线性比线性分析更敏感。

(2)扭曲度:单元面内的扭转和面外的翘曲程度。

(3)疏密过渡:网格的疏密主要表现为应力梯度方向和横向过渡情况,应力集中的情况应妥善处理,而对于分析影响较小的局部特征应分析其情况,如外圆角的影响比内圆角的影响小的多。

(4)节点编号排布:节点编号对于求解过程中的总体刚度矩阵的元素分布、分析耗时、内存及空间有一定的影响。

合理的节点、单元编号有助于利用刚度矩阵对称、带状分布、稀疏矩阵等方法提高求解效率,同时要注意消除重复的节点和单元。

4.装配结构中单元的协调(1)自由度不同的单元不协调:例如,ANSYS中SHELL63、BEAM4和SOLID45三种单元,前二者均包含六个自由度,而Solid45只包含三个平动自由度,因此后者只传递前二者的平动位移,不传递R旋转方向的位移。

(2)有相同自由度的单元不总是协调的:例如,ANSYS中BEAM3和SHELL41单元,Bea m3具备平动方向的三个自由度,而SHELL41包括两个平动自由度(UX/UY)和一个旋转自由度(RTOTZ),因此SHELL41只能传递BEAM3的平动位移,不能传递旋转方向的值。

相关文档
最新文档