有限元网格划分方法与基本原理
第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网格划分原则及方法2009-04-04 23:48引言:划分网格是建立有限元模型的一个重要环节,它要求考虑的问题较多,需要的工作量较大,所划分的网格形式对计算精度和计算规模将产生直接影响。
为建立正确、合理的有限元模型,这里介绍网格划分时的一些基本原则及方法。
关键词:Deform 网格局部细化一、网格划分的原则1 网格数量网格数量的多少将影响计算结果的精度和计算规模的大小。
一般来讲,网格数量增加,计算精度会有所提高,但同时计算规模也会增加,所以在确定网格数量时应权衡两个因数综合考虑。
图1中的曲线1表示结构中的位移随网格数量收敛的一般曲线,曲线2代表计算时间随网格数量的变化。
可以看出,网格较少时增加网格数量可以使计算精度明显提高,而计算时间不会有大的增加。
当网格数量增加到一定程度后,再继续增加网格时精度提高甚微,而计算时间却有大幅度增加。
所以应注意增加网格的经济性。
实际应用时可以比较两种网格划分的计算结果,如果两次计算结果相差较大,可以继续增加网格,相反则停止计算。
图1 位移精度和计算时间随网格数量的变化在决定网格数量时应考虑分析数据的类型。
在静力分析时,如果仅仅是计算结构的变形,网格数量可以少一些。
如果需要计算应力,则在精度要求相同的情况下应取相对较多的网格。
在热分析中,结构内部的温度梯度不大,不需要大量的内部单元,这时可划分较少的网格。
2 网格疏密网格疏密是指在结构不同部位采用大小不同的网格,这是为了适应计算数据的分布特点。
在计算数据变化梯度较大的部位(如应力集中处),为了较好地反映数据变化规律,需要采用比较密集的网格。
而在计算数据变化梯度较小的部位,为减小模型规模,则应划分相对稀疏的网格。
这样,整个结构便表现出疏密不同的网格划分形式。
图2是中心带圆孔方板的四分之一模型,其网格反映了疏密不同的划分原则。
小圆孔附近存在应力集中,采用了比较密的网格。
板的四周应力梯度较小,网格分得较稀。
有限元计算原理
有限元计算原理
有限元计算原理是一种工程分析的方法,用于求解各种结构及连续体的力学问题。
其基本思想是将结构或连续体分割成有限数量的小单元,然后通过对这些小单元进行计算,再将其组合起来求解整体问题。
这种方法可以将结构或连续体的力学行为分析得非常精确,可以获得结构的应力应变分布、位移分布等信息。
有限元计算的原理可以概括为以下几个步骤:
1. 网格划分:将结构或连续体划分成许多小单元,即有限元,这些小单元通过节点连接起来构成整个结构。
2. 求解力学方程:根据结构或连续体的几何形状和物理特性,建立相应的力学方程组。
通常采用弹性力学理论来描述结构或连续体的力学行为。
3. 边界条件的处理:给定结构或连续体的边界条件,如固支、约束力等,在有限元网格中对应的节点上施加相应的约束。
4. 单元刚度矩阵的组装:通过计算每个小单元的刚度矩阵,将其组装成整个结构或连续体的整体刚度矩阵。
5. 单元荷载向量的组装:根据给定的荷载条件,在每个小单元上计算相应的荷载向量,将其组装成整个结构或连续体的荷载向量。
6. 求解位移和应力:根据组装好的整体刚度矩阵和荷载向量,通过求解线性方程组,得到结构或连续体中每个节点的位移和应力。
7. 后处理:根据求解得到的位移和应力,可以计算出结构或连续体的各种物理量,比如应变、应力、变形等。
通过这种有限元计算的方法,可以对各种复杂的结构或连续体进行力学分析和优化设计。
机械零件有限元分析-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方法随着计算机技术的飞速发展,数值模拟方法在工程领域中的应用越来越广泛。
有限元法(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. 单元刚度矩阵单元刚度矩阵是计算每个单元应力和应变的关键。
它由每个单元的杨氏模量、泊松比和几何信息确定。
ANSYS有限元分析中的网格划分
ANSYS有限元分析中的网格划分有限元分析中的网格划分好坏直接关系到模型计算的准确性。
本文简述了网格划分应用的基本理论,并以ANSYS限元分析中的网格划分为实例对象,详细讲述了网格划分基本理论及其在工程中的实际应用,具有一定的指导意义。
作者: 张洪才关键字: CAE ANSYS 网格划分有限元1 引言ANSYS有限元网格划分是进行数值模拟分析至关重要的一步,它直接影响着后续数值计算分析结果的精确性。
网格划分涉及单元的形状及其拓扑类型、单元类型、网格生成器的选择、网格的密度、单元的编号以及几何体素。
从几何表达上讲,梁和杆是相同的,从物理和数值求解上讲则是有区别的。
同理,平面应力和平面应变情况设计的单元求解方程也不相同。
在有限元数值求解中,单元的等效节点力、刚度矩阵、质量矩阵等均用数值积分生成,连续体单元以及壳、板、梁单元的面内均采用高斯(Gauss)积分,而壳、板、梁单元的厚度方向采用辛普生(Simpson)积分。
辛普生积分点的间隔是一定的,沿厚度分成奇数积分点。
由于不同单元的刚度矩阵不同,采用数值积分的求解方式不同,因此实际应用中,一定要采用合理的单元来模拟求解。
2 ANSYS网格划分的指导思想ANSYS网格划分的指导思想是首先进行总体模型规划,包括物理模型的构造、单元类型的选择、网格密度的确定等多方面的内容。
在网格划分和初步求解时,做到先简单后复杂,先粗后精,2D单元和3D单元合理搭配使用。
为提高求解的效率要充分利用重复与对称等特征,由于工程结构一般具有重复对称或轴对称、镜象对称等特点,采用子结构或对称模型可以提高求解的效率和精度。
利用轴对称或子结构时要注意场合,如在进行模态分析、屈曲分析整体求解时,则应采用整体模型,同时选择合理的起点并设置合理的坐标系,可以提高求解的精度和效率,例如,轴对称场合多采用柱坐标系。
有限元分析的精度和效率与单元的密度和几何形状有着密切的关系,按照相应的误差准则和网格疏密程度,避免网格的畸形。
有限元网格剖分原理
有限元网格剖分原理: 1. 引言有限元法是求解复杂工程问题的一种近似数值解法,现已广泛应用到力学、热学、电磁学等各个学科,主要分析工作环境下物体的线性和非线性静动态特性等性能。
有限元法求解问题的基本过程主要包括:分析对象的离散化?有限元求解?计算结果的处理三部分。
曾经有人做过统计:三个阶段所用的时间分别占总时间的40%~50%、5%及50%~55%。
也就是说,当利用有限元分析对象时,主要时间是用于对象的离散及结果的处理。
如果采用人工方法离散对象和处理计算结果,势必费力、费时且极易出错,尤其当分析模型复杂时,采用人工方法甚至很难进行,这将严重影响高级有限元分析程序的推广和使用。
因此,开展自动离散对象及结果的计算机可视化显示的研究是一项重要而紧迫的任务。
可喜的是,随着计算机及计算技术的飞速发展,出现了开发对象的自动离散及有限元分析结果的计算机可视化显示的热潮,使有限元分析的“瓶颈”现象得以逐步解决,对象的离散从手工到半自动到全自动,从简单对象的单维单一网格到复杂对象的多维多种网格单元,从单材料到多种材料,从单纯的离散到自适应离散,从对象的性能校核到自动自适应动态设计/分析,这些重大发展使有限元分析摆脱了仅为性能校核工具的原始阶段,计算结果的计算机可视化显示从简单的应力、位移和温度等场的静动态显示、彩色调色显示一跃成为对受载对象可能出现缺陷(裂纹等)的位置、形状、大小及其可能波及区域的显示等,这种从抽象数据到计算机形象化显示的飞跃是现在甚至将来计算机集成设计/分析的重要组成部分。
2. 有限元分析对网格剖分的要求有限元网格生成就是将工作环境下的物体离散成简单单元的过程,常用的简单单元包括:一维杆元及集中质量元、二维三角形、四边形元和三维四面体元、五面体元和六面体元。
他们的边界形状主要有直线型、曲线型和曲面型。
对于边界为曲线(面)型的单元,有限元分析要求各边或面上有若干点,这样,既可保证单元的形状,同时,又可提高求解精度、准确性及加快收敛速度。
有限元网格划分
本文讨论了有限元网格的重要概念,包括单元的分类、有限元误差的分类与影响因素;并讨论分析结果的收敛性控制方法,并由实例说明了网格质量及收敛性对取得准确分析结果的重要性。
同时讨论了一些重要网格控制的建议及其他网格设定的说明。
一、基本有限元网格概念1.单元概述几何体划分网格之前需要确定单元类型。
单元类型的选择应该根据分析类型、形状特征、计算数据特点、精度要求和计算的硬件条件等因素综合考虑。
为适应特殊的分析对象和边界条件,一些问题需要采用多种单元进行组合建模。
2.单元分类选择单元首先需要明确单元的类型,在结构有限元分析中主要有以下一些单元类型:平面应力单元、平面应变单元、轴对称实体单元、空间实体单元、板单元、壳单元、轴对称壳单元、杆单元、梁单元、弹簧单元、间隙单元、质量单元、摩擦单元、刚体单元和约束单元等。
根据不同的分类方法,上述单元可以分成以下不同的形式。
3.按照维度进行单元分类根据单元的维数特征,单元可以分为一维单元、二维单元和三维单元。
一维单元的网格为一条直线或者曲线。
直线表示由两个节点确定的线性单元。
曲线代表由两个以上的节点确定的高次单元,或者由具有确定形状的线性单元。
杆单元、梁单元和轴对称壳单元属于一维单元,如图1~图3所示。
二维单元的网格是一个平面或者曲面,它没有厚度方向的尺寸。
这类单元包括平面单元、轴对称实体单元、板单元、壳单元和复合材料壳单元等,如图4所示。
二维单元的形状通常具有三角形和四边形两种,在使用自动网格剖分时,这类单元要求的几何形状是表面模型或者实体模型的边界面。
采用薄壳单元通常具有相当好的计算效率。
三维单元的网格具有空间三个方向的尺寸,其形状具有四面体、五面体和六面体,这类单元包括空间实体单元和厚壳单元,如图5所示。
在自动网格划分时,它要求的是几何模型是实体模型(厚壳单元是曲面也可以)。
4.按照插值函数进行单元分类根据单元插值函数多项式的最高阶数多少,单元可以分为线性单元、二次单元、三次单元和更高次的单元。
《2024年有限元网格剖分与网格质量判定指标》范文
《有限元网格剖分与网格质量判定指标》篇一一、引言有限元法是一种广泛应用于工程和科学计算中的数值分析方法。
其核心步骤之一是进行网格剖分,即将求解域划分为一系列小的、相互连接的子域或元素。
网格的质量直接影响到有限元分析的准确性和效率。
因此,本文将重点讨论有限元网格剖分的方法以及网格质量的判定指标。
二、有限元网格剖分1. 网格剖分的基本原则有限元网格剖分应遵循以下基本原则:一是尽可能保持单元的规则性,如六面体单元;二是确保网格的连续性和兼容性;三是考虑网格的适应性,以适应求解域的几何形状和边界条件;四是尽可能减少单元的数量,以节省计算资源。
2. 常见的网格剖分方法(1)自动剖分法:利用计算机程序自动进行网格剖分,如基于Delaunay三角化的剖分方法。
(2)映射法:将求解域映射到参数空间进行剖分,再映射回原空间得到网格。
(3)手动剖分法:根据求解域的几何形状和边界条件,手动进行网格剖分。
三、网格质量判定指标1. 单元形态指标(1)扭曲度(Skewness):用于衡量单元的形状与理想形状的偏差程度,扭曲度越大,单元的形状越不规则,影响计算的精度和效率。
(2)内角分布:单元的内角应尽可能接近标准值(如四边形单元为90度),内角分布的均匀性可以反映单元的规则性。
(3)面积/体积变化率:用于衡量单元尺寸变化对整体网格的影响,变化率越小,网格质量越好。
2. 连接性指标(1)节点连接数:每个节点的连接单元数应适中,过多或过少的连接都可能导致计算误差。
(2)相邻单元的协调性:相邻单元在公共边界上应具有良好的协调性,避免出现不连续或重复的单元边界。
3. 整体性指标(1)网格均匀性:整体网格的尺寸和密度应保持均匀,避免出现过大或过小的单元。
(2)边界拟合度:网格应尽可能贴合求解域的边界,提高边界条件的准确性。
四、结论有限元网格剖分是有限元法的重要步骤之一,而网格质量直接影响到有限元分析的准确性和效率。
本文介绍了有限元网格剖分的基本原则和常见方法,以及网格质量的判定指标。
有限元网格划分方法与基本原理
结构有限元分析中的网格划分技术及其应用实例结构有限元分析中的网格划分是否直接关系到解算的效果。
本文简述了网格划分应用的基本理论,并以空间自由曲面覆盖件和大型整体网络钢筋壳体产品的有限元分析中的网格划分为实例对象,详细讲述了空间自由和三维实体的网格划分基本理论及其在工程中的实际应用,非常具有现实意义和借鉴价值。
一、前言有限元网格划分是进行有限元数值模拟分析至关重要的一步,它直接影响着后续数值计算分析结果的精确性。
网格划分涉及单元的形状及其拓扑类型、单元类型、网格生成器的选择、网格的密度、单元的编号以及几何体素。
从几何表达上讲,梁和杆是相同的,从物理和数值求解上讲则是有区别的。
同理,平面应力和平面应变情况设计的单元求解方程也不相同。
在有限元数值求解中,单元的等效节点力、刚度矩阵、质量矩阵等均用数值积分生成,连续体单元以及壳、板、梁单元的面内均采用高斯(Gauss)积分,而壳、板、梁单元的厚度方向采用辛普生(Simpson)积分。
辛普生积分点的间隔是一定的,沿厚度分成奇数积分点。
由于不同单元的刚度矩阵不同,采用数值积分的求解方式不同,因此实际应用中,一定要采用合理的单元来模拟求解。
CAD软件中流行的实体建模包括基于特征的参数化建模和空间自由曲面混合造型两种方法。
Pro/E和S oildWorks是特征参数化造型的代表,而 CATIA与Unigraphics等则将特征参数化和空间自由曲面混合造型有机的结合起来。
现有CAD软件对表面形态的表示法已经大大超过了CAE软件,因此,在将CAD实体模型导入CAE软件的过程中,必须将CAD模型中其他表示法的表面形态转换到CAE软件的表示法上,转换精度的高低取决于接口程序的好坏。
在转换过程中,程序需要解决好几何图形(曲线与曲面的空间位置)和拓扑关系(各图形数据的逻辑关系)两个关键问题。
其中几何图形的传递相对容易实现,而图形间的拓扑关系容易出现传递失败的情况。
数据传递面临的一个重大挑战是,将导入CAE程序的CAD模型改造成适合有限元分析的网格模型。
有限元法基本原理及应用 尹飞鸿 课件
有限元法基本原理及应用尹飞鸿课件
有限元法基本原理及应用
有限元法是一种数值计算方法,用于求解复杂结构的物理问题。
它可通过将物理系统分割成许多小的有限元素来近似描述系统行为,并根据元素之间的关系和物理方程求解系统的响应。
有限元法的基本原理是建立数学模型,将连续体划分为多个离散的有限元素。
每个有限元素代表了原问题的一个小区域,具有一定的属性和形状。
通过将元素的局部行为进行组装,可以重建出整个物理系统的行为。
有限元法的应用非常广泛,涵盖了许多工程领域。
在结构力学中,有限元法可用于分析和优化建筑、航空航天器、机械设备等的力学性能。
在流体力学中,有限元法可用于模拟流体流动、传热和传质等问题。
在电磁学中,有限元法可用于计算电磁场和电磁波的分布。
有限元法的应用过程包括模型建立、划分网格、选取适当的数值计算方法以及求解和后处理结果等步骤。
模型建立是指将物理问题转化为数学描述,包括确定几何形状、材料性质和加载条件等。
划分网格是将物理模型分割成有限元素,通过合适的网格划分可以提高计算效率和精度。
数值计算方法是选择适当的数值算法来求解离散化的模型方程。
求解和后处理结果是对模拟结果进行分析和可视化展示。
总之,有限元法基于分割和离散化的思想,是一种强大的数值计算方法。
通过应用有限元法,我们可以更好地理解和解决复杂的物理问题,提高工程设计的效率和可靠性。
有限元网格划分原理
有限元网格划分原理
有限元网格划分原理是一种用于计算领域离散化的数值方法。
它将连续的领域划分为有限数量的小单元,每个小单元称为有限元。
这些有限元可被视为数学模型中的局部区域,其内部的物理过程可以被近似为线性或非线性的形式。
有限元网格划分原理的目标是将整个领域划分为足够多的有限元,以便能够准确地描述所研究的问题。
划分时需要考虑几何形状、边界条件、计算资源等因素,以获取一个高效且准确的离散模型。
通常,将整个领域划分为小单元可以更好地逼近真实物理过程,并提供对系统行为的详细理解。
在有限元网格划分过程中,首先确定领域的几何形状和边界条件。
然后,选择适当的离散化方法,将领域划分为小单元,如三角形、四边形、六边形或四面体。
每个小单元内的变量以形函数的形式进行逼近,形函数可根据问题的特点进行选择。
一旦完成网格划分,就可以在每个有限元中设置数学方程,在整个领域上建立一个代数系统。
该系统由一系列线性或非线性方程组组成,其中每个方程对应于一个小单元。
通过求解这些方程,可以获得在整个领域中变量的近似解。
有限元网格划分原理的核心思想是将复杂问题转化为简单的局部问题,并通过将领域划分为小单元来近似描述整个系统。
通过调整网格大小和形状,可以调整计算精度和效率。
因此,有限元网格划分原理是计算力学、结构力学、流体力学等领域中常用的数值方法之一。
有限元方法基本原理
有限元方法基本原理有限元方法(Finite Element Method, FEM)是一种数值计算方法,主要用于求解偏微分方程的数值解。
它最早由Courant、Bubnov和Galerkin等人在20世纪50年代提出,并在以后的几十年中得到了广泛的发展和应用。
有限元方法的基本原理是将要求解的区域分割成若干个小的子区域,通常称为有限元,每个有限元内部的物理量可以用一个简单的数学表达式来表示。
然后,通过在有限元之间建立连续性条件,将整个问题转化为一组代数方程,进而得到数值解。
有限元方法的基本步骤包括:建立有限元模型、离散化、建立代数方程、求解代数方程和后处理。
下面将详细介绍每个步骤的具体内容。
第一步,建立有限元模型。
该步骤主要是对要求解的问题进行数学建模,包括选择适当的坐标系、定义物理量和约束条件等。
通常,物理问题可以通过连续介质假设,将其离散化为一组小的有限元。
第二步,离散化。
将要求解的区域划分为有限个小的子区域,通常称为有限元。
常见的有限元形状包括三角形、四边形和六面体等。
有限元的选择通常是根据问题的几何形状和物理条件来确定的。
第三步,建立代数方程。
有限元方法的核心是建立代数方程,用于描述物理问题在离散点上的数值解。
代数方程通常是通过施加适当的数学形式和边界条件来建立的。
建立代数方程的基本思想是使用一组试验函数来近似描述有限元内部的解。
通常采用Galerkin方法,即在离散点上进行加权残差积分,使得残差的加权平均为零。
第四步,求解代数方程。
一旦代数方程建立完成,就可以使用数值方法求解这组代数方程。
常见的求解方法包括直接法和迭代法等。
直接法适用于方程较小的情况,而迭代法适用于方程较大的情况。
常见的迭代法有Jacobi迭代法、Gauss-Seidel迭代法和共轭梯度法等。
第五步,后处理。
求解代数方程后,需要对结果进行后处理和分析。
后处理包括计算和显示物理量、绘制图形以及进行误差估计等。
通过后处理,可以对模型进行验证,并对结果进行解释和解释。
有限元分析原理与步骤
有限元分析原理与步骤
有限元分析是一种数值计算方法,用于解决工程结构的力学问题。
它将任意复杂的结构分割成为若干个简单的子结构,通过数学模型和计算机软件进行力学分析。
有限元分析的步骤如下:
1. 建立几何模型:根据实际结构的几何形状,使用CAD软件
或者手工绘图等方式建立三维或二维模型。
2. 网格划分:将结构模型划分成若干个小单元,如三角形、四边形或六边形等,这些小单元构成了有限元网格。
3. 选择适当的元素类型:根据结构的特性选择合适的元素类型,如杆件元、梁单元、板单元等。
4. 建立整体刚度矩阵:根据每个小单元的几何形状和材料性质,计算每个小单元的刚度矩阵,将其组装成整个结构的刚度矩阵。
5. 施加边界条件:确定结构的边界条件,如固定支座、约束等。
6. 施加荷载:施加力、压力、温度等荷载条件。
7. 求解方程:通过求解结构的刚度方程,得到结构的位移、应力、应变等结果。
8. 后处理结果:根据求解得到的结果,进行结果的可视化及分
析。
通过以上步骤,有限元分析可以提供结构的力学性能分析,如应力、应变、变形等,为工程设计和优化提供参考依据。
有限元计算单元网格划分
有限元语言及编译器(Finite Element Language And it’s Compiler,以下简称FELAC) 是中国科学院数学与系统科学研究院梁国平研究院于1983年开始研发的通用有限元软件平 台,是具有国际独创性的有限元计算软件,是PFEPG系列软件三十年成果(1983年—2013 年)的总结与提升,有限元语言语法比PFEPG更加简练,更加灵活,功能更加强大。目前 已发展到2.0版本。其核心采用元件化思想来实现有限元计算的基本工序,采用有限元语 言来书写程序的代码,为各领域,各类型的有限元问题求解提供了一个极其有力的工具。 FELAC可以在数天甚至数小时内完成通常需要一个月甚至数月才能完成的编程面或体单元过渡
1)从小单元到大单元过渡时,应使同一节点所连接的单元不致相差太大,避免突然过渡现象。 通常用计算结果调整,保证同一节点所连接的单元精度值至少在0.1以下。单元精度值根据单元内节 点应力与节点平均应力的误差计算。
2)难于过度处最好使用过渡单元,过渡单元的使用要比用同一单元勉强过渡的计算结果要好。 例如:对于复杂体结构间的过渡,最好使用“金子塔”单元过渡。 (5)面或体转接部位的单元 几何模型圆角过渡处的单元划分,根据弧长对应的圆心角和半径确定,对于半径为3mm左右、 圆心角大于90度的转接弧长,通常至少要划分3~4个单元。 (6)高应力区的单元 对高应力区,要进行网格细分应力稳定性计算。即采用多次局部网格细分并进行计算,当前、 后两次计算结果满足所需的精度要求时(通常要求小于0.03)确定网格。 总之,几何模型网格划分时,要在单元类型、单元形态、单元大小、单元过渡和局部应力稳定 等方面下功夫,才能满足工程上的精度要求,达到预期的结果。
有限元法的基本步骤
有限元法的基本步骤有限元法是一种用于求解较为复杂的实际工程问题的数值分析方法。
它将一个连续的物体或系统划分为许多小的单元,然后通过建立在这些单元上的数学方程来模拟和求解实际问题。
在这篇文章中,我们将探讨有限元法的基本步骤,并深入讨论其原理和应用。
1. 确定问题的边界和几何形状在使用有限元法求解实际问题之前,需要先确定问题的边界和几何形状。
通常情况下,问题的边界需要定义为固定边界或自由边界,以便在数学模型中进行处理。
问题的几何形状也需要被建模和描述,这样才能得到准确的计算结果。
2. 划分网格划分网格是有限元法中非常重要的一步。
网格划分是将问题的几何形状划分为一系列小的单元。
这些小单元称为有限元,它们可以是三角形、四边形或其他形状。
网格的划分需要根据问题的几何形状和求解精度来确定,并且需要保证各个有限元之间具有充分的连续性和相互联系,以确保模拟结果的准确性和可靠性。
3. 建立数学模型和方程在确定问题的边界和划分网格之后,下一步是建立与物理现象相关的数学模型和方程。
根据问题的具体情况,可以使用不同类型的方程,如静力学方程、热传导方程、流体力学方程等。
这些方程将物理现象转化为数学表达式,并可以通过有限元法进行求解。
4. 应用边界条件在建立数学模型和方程之后,需要应用边界条件。
边界条件可以是物体的固定边界条件,如固定端或自由端;也可以是物体的外部边界条件,如外力、温度等。
边界条件的正确应用对于求解实际问题非常重要,它们将影响模拟结果的准确性和可靠性。
5. 求解数学方程一旦建立了数学模型、划分网格并应用了边界条件,下一步就是使用数值方法求解数学方程。
有限元法将整个问题转化为一个求解代数方程组的问题,并通过迭代方法求解。
求解过程中需要根据初始条件和边界条件进行迭代计算,直到得到收敛的解。
通过以上的基本步骤,我们可以使用有限元法对复杂的实际工程问题进行数值求解。
有限元法的优点在于可以模拟各种不同的物理现象,并且可以对复杂的几何形状进行建模和求解。
有限元法基本原理
有限元法基本原理
有限元法是解决偏微分方程数值分析问题的重要方法,它根据力
学原理将构件表示成一系列有限个拉普拉斯单元,采用有限个有限量
节点在某种元素的基质上建立的模型来近似求解构件的本构关系。
它
将复杂的本构关系准确地还原为有限数量的有限单元,以此分析不同
物理状态下物体受力和变形机制,可用于弹性、非线性动力学分析及
多物理场耦合场景等复杂问题的分析。
有限元法由三部分组成:网格划分、体积单元的本构建立及节点
的采样,它将整个物体划分成几种封闭的体积单元,选取合适的节点
对每一种单元进行采样,并为各种单元类型形成有适用的本构关系方程,以串联每个构件的局部分析结果。
首先,在网格划分方面,有限元法可以通过不同的体积单元划分、节点采样及本构关系来处理复杂的问题,如曲面、孔洞等,形成封闭
的有限元网格,随后,对复杂的本构关系准确地还原为有限个有限单元,即针对每一种单元类型的形变量,采取合适的节点、布点一系列
的坐标。
最后,有限元法利用耦合方程作为求解强度和变形问题的基础,
在此基础上,有限元法可以应用于多物理场、非线性动力学分析及其
他复杂的物理状态场景。
另外,它还可以帮助测量构件受力和变形机制,使得构件能正确适应环境变化。
由于有限元法处理方法较为简单,而且力学原理深入,因此,它已在工程计算中得到广泛的应用,有效
提高了模型的准确性和精确度,为进一步探索物理现象带来了巨大的
方便。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 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的平动位移,不能传递旋转方向的值。