ANSYS结构有限元分析中的网格划分技术及其应用实例
- 格式:doc
- 大小:76.50 KB
- 文档页数:9
【分享】复杂几何模型的系列网格划分技术众所周知,对于有限元分析来说,网格划分是其中最关键的一个步骤,网格划分的好坏直接影响到解算的精度和速度。
在ANSYS中,大家知道,网格划分有三个步骤:定义单元属性(包括实常数)、在几何模型上定义网格属性、划分网格。
在这里,我们仅对网格划分这个步骤所涉及到的一些问题,尤其是与复杂模型相关的一些问题作简要阐述。
一、自由网格划分自由网格划分是自动化程度最高的网格划分技术之一,它在面上(平面、曲面)可以自动生成三角形或四边形网格,在体上自动生成四面体网格。
通常情况下,可利用ANSYS的智能尺寸控制技术(SMARTSIZE命令)来自动控制网格的大小和疏密分布,也可进行人工设置网格的大小(AESIZE、LESIZE、KESIZE、ESIZE等系列命令)并控制疏密分布以及选择分网算法等(MOPT命令)。
对于复杂几何模型而言,这种分网方法省时省力,但缺点是单元数量通常会很大,计算效率降低。
同时,由于这种方法对于三维复杂模型只能生成四面体单元,为了获得较好的计算精度,建议采用二次四面体单元(92号单元)。
如果选用的是六面体单元,则此方法自动将六面体单元退化为阶次一致的四面体单元,因此,最好不要选用线性的六面体单元(没有中间节点,比如45号单元),因为该单元退化后为线性的四面体单元,具有过刚的刚度,计算精度较差;如果选用二次的六面体单元(比如95号单元),由于其是退化形式,节点数与其六面体原型单元一致,只是有多个节点在同一位置而已,因此,可以利用TCHG命令将模型中的退化形式的四面体单元变化为非退化的四面体单元,减少每个单元的节点数量,提高求解效率。
在有些情况下,必须要用六面体单元的退化形式来进行自由网格划分,比如,在进行混合网格划分(后面详述)时,只有用六面体单元才能形成金字塔过渡单元。
对于计算流体力学和考虑集肤效应的电磁场分析而言,自由网格划分中的层网格功能(由LESIZE命令的LAYER1和LAYER2域控制)是非常有用的。
ansys有限元求解基本方法-回复ANSYS有限元求解基本方法有限元法(Finite Element Method,FEM)是一种高效且广泛应用于工程领域的数值分析方法。
它将复杂的实际结构问题转化为计算机理解的离散化网格,然后应用数值方法对这个网格进行计算,以求解结构的行为和性能。
ANSYS是目前应用最广泛的有限元软件之一,本文将介绍ANSYS 中有限元求解的基本方法。
1. 网格划分(Meshing):在使用ANSYS进行有限元分析之前,首先需要将复杂的实际结构转化为有限元网格。
网格划分是有限元分析的第一步,它直接影响到后续的求解精度和计算效率。
ANSYS提供了多种网格划分工具,包括自动划分和手动划分。
自动划分是指ANSYS根据用户设定的参数自动生成网格,手动划分是指用户手动绘制网格。
2. 定义材料和边界条件:在进行有限元分析之前,需要定义材料的力学性质和边界条件。
材料的力学性质包括弹性模量、泊松比、屈服强度等,边界条件包括约束和外载荷等。
ANSYS提供了简单易用的界面,可以方便地输入这些参数。
3. 选择求解器和求解方法:在网格划分和参数定义完成后,需要选择适当的求解器和求解方法。
ANSYS提供了多种求解器和求解算法,用于求解不同类型的问题。
选择合适的求解器和求解方法可以提高计算效率和求解精度。
4. 求解并后处理:在进行有限元求解之前,可以进行预处理操作,如自适应网格划分、模型简化等,以提高求解效率。
然后,通过点击求解按钮,ANSYS将自动进行有限元求解。
求解完成后,可以进行后处理操作,如显示位移、应力、应变等结果,以及生成图形和报表等。
ANSYS在有限元求解过程中还提供了许多高级功能,如非线性分析、动力学分析、热传导分析等。
这些高级功能可以进一步扩展ANSYS的应用范围和分析能力。
综上所述,ANSYS有限元求解的基本方法包括网格划分、定义材料和边界条件、选择求解器和求解方法、求解并后处理。
通过这些步骤,可以对复杂的实际结构进行准确、可靠的分析和设计,为工程实践提供重要的支持和指导。
ANSYS模态分析教程及实例讲解ANSYS是一款常用的有限元分析软件,可以用于执行结构分析、热分析、流体分析等多种工程分析。
模态分析是其中的一项重要功能,用于计算和分析结构的固有振动特性,包括固有频率、振型和振动模态,可以帮助工程师了解和优化结构的动态响应。
以下是一份ANSYS模态分析教程及实例讲解,包含了基本步骤和常用命令,帮助读者快速上手模态分析。
1.创建模型:首先需要创建模型,在ANSYS界面中构建出待分析的结构模型,包括几何形状、材料属性和边界条件等。
可以使用ANSYS的建模工具,也可以导入外部CAD模型。
2.网格划分:在模型创建完毕后,需要进行网格划分,将结构划分为小的单元,使用ANSYS的网格划分功能生成有限元网格。
网格划分的细腻程度会影响分析结果的准确性和计算时间,需要根据分析需要进行合理选择。
3.设置材料属性:在模型和网格创建完毕后,需要设置材料属性,包括弹性模量、密度和材料类型等。
可以通过ANSYS的材料库选择已有的材料属性,也可以自定义材料属性。
4.定义边界条件:在模型、网格和材料属性设置完毕后,需要定义结构的边界条件,包括约束和加载条件。
约束条件是指结构受限的自由度,例如固定支撑或限制位移;加载条件是指施加到结构上的载荷,例如重力或外部力。
5.运行模态分析:完成前面几个步骤后,就可以执行模态分析了。
在ANSYS中,可以使用MODAL命令来进行模态分析。
MODAL命令需要指定求解器和控制选项,例如求解的模态数量、频率范围和收敛准则等。
6.分析结果:模态分析完成后,ANSYS会输出结构的振动特性,包括固有频率、振型和振动模态。
可以使用POST命令查看和分析分析结果,例如绘制振动模态或振动模态的频率响应。
下面是一个实际的案例,将使用ANSYS执行模态分析并分析分析结果。
案例:矩形板的模态分析1.创建模型:在ANSYS界面中创建一个矩形板结构模型,包括矩形板的几何形状和材料属性等。
刘丽贤,马国鹭,赵登峰:基于APDL的ANSYS网格划分及应用拳DIM,LA,ARRAY,8,51将实体所属线号放入LA二维数组中’kSET,LA(1,1),15,16,4,3,19,54,24,521给二维LA数组赋值水SET,LA(1,5),2,1,5,6,12,8,7,11木DIM,LB,ARRAY,51将对不同线段划分的段数放入LB一维数组中LB(1)=20,30,10,20,40堆DO,AB,l,51用双重循环按照设定的段数划分实体所属线木DO,C。
1,8LESIZE,LA(C,AB),,,LB(AB),,,,,l!调用线号数组LA并通过调用段数数组LB对其设置划分的段数拳ENDDO!结束嵌套循环木ENDDO!结束外部循环木DIM,LD,ARRAY,9,2LD(1,1)=1,13,24,4,9,7,16,10,141将要被划分网格的实体编号按照一定顺序放入二维数组LD中LD(1,2)=2,26,21,27,25,23,22,15,17,kDO,LLD,l,21用双重循环划分实体MAT.LID!给被划分的实体赋材料属性REAL。
LED!给被划分的实体赋实常数MSHAPE.0.2DMSHKEY,1,IcDO。
U正,1,9ASEL,S,LD(LLE,LED)!调用LD数组选中将被划分的实体AMESH,LD(UJE,LLD)!调用LD数组选中划分的实体宰ENDDO!结束嵌套循环木ENDDO!结束外部循环以上是对该型电视机的CRT网格划分的命令流,以相同方式划分电视机壳体肋板和前后外壳等部件。
APDL划分网格b自由划分f。
g梧图2实体模型网格划分图2是用APDL通过以上方式和自由划分网格对长虹SF21366型电视机的实体模型的网格划分。
由单元信息表1可反映出通过APDL划分的电视机有限元模型的网格质量较好。
表1单元信息对比表3结论ANSYS软件经过几十年的发展,日趋成熟。
它不但具有良好的数据库管理和强大的前后处理功能,而且还时刻追踪先进的计算方法和计算机技术,不断提高分析精度和扩展自身的开放性,并提供良好的二次开发功能。
0 前言利用ANSYS分析钢筋混凝土结构时,其有限元模型主要有分离式和整体式两种模型。
这里结合钢筋混凝土材料的工作特性,从模型建立到非线性计算再到结果分析的全过程讲述了利用ANSYS进行钢筋混凝土结构分析的方法与技巧,并以钢筋混凝土简支梁为例,采用分离式有限元模型,说明其具体应用。
1 单元选取与材料性质1. 1 混凝土单元ANSYS中提供了上百种计算单元类型,其中Solid65单元是专门用于模拟混凝土材料的三维实体单元。
该单元是八节点六面体单元,每个节点具有三个方向的自由度( UX , UY , UZ) 。
在普通八节点线弹性单元Solid45 的基础上,该单元增加了针对于混凝土的材性参数和组合式钢筋模型,可以综合考虑包括塑性和徐变引起的材料非线性、大位移引起的几何非线性、混凝土开裂和压碎引起的非线性等多种混凝土的材料特性。
使用Solid65 单元时,一般需要为其提供如下数据:1)、实常数(Real Constants) :定义弥散在混凝土中的最多三种钢筋的材料属性,配筋率和配筋角度。
对于墙板等配筋较密集且均匀的构件,一般使用这种整体式钢筋混凝土模型。
如果采用分离式配筋,那么此处则不需要填写钢筋实常数。
2)、材料模型(Material Model) :在输入钢筋和混凝土的非线性材料属性之前,首先必须定义钢筋和混凝土材料在线弹性阶段分析所需的基本材料信息,如:弹性模量,泊松比和密度。
3)、数据表(Data Table) :利用数据表进一步定义钢筋和混凝土的本构关系。
对于钢筋材料,一般只需要给定一个应力应变关系的数据表就可以了,譬如双折线等强硬化(bilinear isotropic hardening)或随动硬化模型( kinematic hardening plasticity)等。
而对于混凝土模型,除需要定义混凝土的本构关系外,还需要定义混凝土材料的破坏准则。
在ANSYS中,常用于定义混凝土本构关系的模型有:1)多线性等效强化模型(Multilinear isotropic hardening plas2ticity ,MISO模型),MISO模型可包括20条不同温度曲线,每条曲线可以有最多100个不同的应力-应变点;2)多线性随动强化模型(Multilinear kinematic hardening plas2ticity ,MKIN 模型),MKIN 模型最多允许5个应力-应变数据点;3)Drucker2Prager plasticity(DP)模型。
基于ANSYS软件的1+6钢丝绳网格划分策略及仿真基于ANSYS软件的1+6钢丝绳网格划分策略及仿真【摘要】钢丝绳作为一种常用的工程材料,在结构、航空、汽车等领域得到了广泛应用。
钢丝绳的应力分布与网格划分密切相关,因此,本文基于ANSYS软件对1+6钢丝绳进行了网格划分策略的研究,并仿真分析了其力学性能。
【关键词】ANSYS;1+6钢丝绳;网格划分;仿真1. 引言钢丝绳是由多股金属丝绳或其他材料绞合而成,其具有高强度、耐腐蚀性能好等优点,在结构工程、航天工程等领域有着广泛应用。
在工程实践中,为了准确计算钢丝绳的力学性能,常常需要使用有限元分析软件进行仿真。
2. ANSYS软件的基本原理ANSYS是一种常用的有限元分析软件,可以对各种结构进行力学仿真,因此非常适合用于钢丝绳的网格划分和力学性能仿真。
ANSYS软件具有强大的前后处理能力,可以实现复杂结构的建模、网格划分以及结果的后处理。
3. 1+6钢丝绳的网格划分策略在对1+6钢丝绳进行网格划分时,需要考虑到以下几个因素:(1)钢丝绳的几何形状:钢丝绳由一根主绳和六根辅绳组成,每根绳线的直径、长度等参数都需要考虑进去。
(2)网格密度:钢丝绳的网格划分密度对仿真结果的准确性有重要影响,因此需要根据预期的结果精度进行合理的网格划分。
(3)约束条件:钢丝绳通常会受到一定的约束条件限制,比如固定端、自由端等,这些条件也需要在网格划分过程中考虑进去。
在基于ANSYS软件进行网格划分时,可以采用以下步骤:(1)建立几何模型:根据实际钢丝绳的几何形状,使用ANSYS软件中的建模工具创建几何模型,包括主绳和辅绳的位置、长度等参数。
(2)网格划分:根据上一步建立的几何模型,使用ANSYS软件中的网格划分工具对钢丝绳进行网格划分。
在网格划分时需要注意细节部分的网格密度,可以使用不同大小的网格对主绳和辅绳进行划分。
(3)加载边界条件:根据实际情况,设置钢丝绳的边界条件,比如约束条件等。
ANSYS案例——20例ANSYS经典实例】针对【典型例题】3.3.7(1)的模型,即如图3-19所示的框架结构,其顶端受均布力作用,用有限元方法分析该结构的位移。
结构中各个截面的参数都为:113.010PaE=,746.510mI-=,426.810mA-=,相应的有限元分析模型见图3-20。
在ANSYS平台上,完成相应的力学分析。
图3-19框架结构受一均布力作用图3-20单元划分、节点位移及节点上的外载解答对该问题进行有限元分析的过程如下。
1.基于图形界面的交互式操作(tepbytep)(1)进入ANSYS(设定工作目录和工作文件)程序→ANSYS→ANSYSInteractive→Workingdirectory(设置工作目录)→Initialjobname(设置工作文件名):beam3→Run→OK(2)设置计算类型(3)选择单元类型(4)定义材料参数ANSYSMainMenu:Preproceor→MaterialProp→MaterialModel→Struc tural→Linear→Elatic→Iotropic:E某:3e11(弹性模量)→OK→鼠标点击该窗口右上角的“”来关闭该窗口(5)定义实常数以确定平面问题的厚度ANSYSMainMenu:Preproceor→RealContant…→Add/Edit/Delete→Add→Type1Beam3→OK→RealContantSetNo:1(第1号实常数),Cro-ectionalarea:6.8e-4(梁的横截面积)→OK→Cloe(6)生成几何模型生成节点ANSYSMainMenu:Preproceor→Modeling→Creat→Node→InActiveCS→Nodenumber1→某:0,Y:0.96,Z:0→Apply→Nodenumber2→某:1.44,Y:0.96,Z:0→Apply→Nodenumber3→某:0,Y:0,Z:0→Apply→Nodenumber4→某:1.44,Y:0,Z:0→OK生成单元ANSYSMainMenu:Preproceor→Modeling→Create→Element→AutoNum bered→ThruNode→选择节点1,2(生成单元1)→apply→选择节点1,3(生成单元2)→apply→选择节点2,4(生成单元3)→OK(7)模型施加约束和外载左边加某方向的受力ANSYSMainMenu:Solution→DefineLoad→Apply→Structural→Force/Moment→OnNode→选择节点1→apply→Directionofforce:F某→VALUE:3000→OK→上方施加Y方向的均布载荷ANSYSMainMenu:Solution→DefineLoad→Apply→Structural→Preure→OnBeam→选取单元1(节点1和节点2之间)→apply→VALI:4167→VALJ:4167→OK左、右下角节点加约束(8)分析计算(9)结果显示(10)退出系统(11)计算结果的验证与MATLAB支反力计算结果一致。
ANSYS有限元分析中的网格划分有限元分析中的网格划分好坏直接关系到模型计算的准确性。
本文简述了网格划分应用的基本理论,并以ANSYS限元分析中的网格划分为实例对象,详细讲述了网格划分基本理论及其在工程中的实际应用,具有一定的指导意义。
作者: 张洪才关键字: CAE ANSYS 网格划分有限元1 引言ANSYS有限元网格划分是进行数值模拟分析至关重要的一步,它直接影响着后续数值计算分析结果的精确性。
网格划分涉及单元的形状及其拓扑类型、单元类型、网格生成器的选择、网格的密度、单元的编号以及几何体素。
从几何表达上讲,梁和杆是相同的,从物理和数值求解上讲则是有区别的。
同理,平面应力和平面应变情况设计的单元求解方程也不相同。
在有限元数值求解中,单元的等效节点力、刚度矩阵、质量矩阵等均用数值积分生成,连续体单元以及壳、板、梁单元的面内均采用高斯(Gauss)积分,而壳、板、梁单元的厚度方向采用辛普生(Simpson)积分。
辛普生积分点的间隔是一定的,沿厚度分成奇数积分点。
由于不同单元的刚度矩阵不同,采用数值积分的求解方式不同,因此实际应用中,一定要采用合理的单元来模拟求解。
2 ANSYS网格划分的指导思想ANSYS网格划分的指导思想是首先进行总体模型规划,包括物理模型的构造、单元类型的选择、网格密度的确定等多方面的内容。
在网格划分和初步求解时,做到先简单后复杂,先粗后精,2D单元和3D单元合理搭配使用。
为提高求解的效率要充分利用重复与对称等特征,由于工程结构一般具有重复对称或轴对称、镜象对称等特点,采用子结构或对称模型可以提高求解的效率和精度。
利用轴对称或子结构时要注意场合,如在进行模态分析、屈曲分析整体求解时,则应采用整体模型,同时选择合理的起点并设置合理的坐标系,可以提高求解的精度和效率,例如,轴对称场合多采用柱坐标系。
有限元分析的精度和效率与单元的密度和几何形状有着密切的关系,按照相应的误差准则和网格疏密程度,避免网格的畸形。
一、前言有限元网格划分是进行有限元数值模拟分析至关重要的一步,它直接影响着后续数值计算分析结果的精确性。
网格划分涉及单元的形状及其拓扑类型、单元类型、网格生成器的选择、网格的密度、单元的编号以及几何体素。
从几何表达上讲,梁和杆是相同的,从物理和数值求解上讲则是有区别的。
同理,平面应力和平面应变情况设计的单元求解方程也不相同。
在有限元数值求解中,单元的等效节点力、刚度矩阵、质量矩阵等均用数值积分生成,连续体单元以及壳、板、梁单元的面内均采用高斯(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、ADINA等)的建模功能都不是很强,尤其是在面对包含复杂空间曲面的产品结构时表现出明显的不足,同时不利于建立后续的单元网格划分模型。
因此,利用现有CAD 平台(如CATIA、UGII、PRO/E)完成网格划分工作,或借助专业网格划分软件HyperMesh、AIEnviroment等来完成任务是比较好的方法。
下面分别以包含大量空间自由曲面的汽车覆盖件产品和宇航业中常用的大型整体网格筋壳体为对象,简述有限元网格划分的基本原理方法和应用。
图1 IGES文件输出的图素信息分页二、有限元网格划分方法与基本原理1.有限元网格划分的指导思想有限元网格划分的指导思想是首先进行总体模型规划,包括物理模型的构造、单元类型的选择、网格密度的确定等多方面的内容。
在网格划分和初步求解时,做到先简单后复杂,先粗后精,2D单元和3D单元合理搭配使用。
为提高求解的效率要充分利用重复与对称等特征,由于工程结构一般具有重复对称或轴对称、镜象对称等特点,采用子结构或对称模型可以提高求解的效率和精度。
利用轴对称或子结构时要注意场合,如在进行模态分析、屈曲分析整体求解时,则应采用整体模型,同时选择合理的起点并设置合理的坐标系,可以提高求解的精度和效率,例如,轴对称场合多采用柱坐标系。
有限元分析的精度和效率与单元的密度和几何形状有着密切的关系,按照相应的误差准则和网格疏密程度,避免网格的畸形。
在网格重划分过程中常采用曲率控制、单元尺寸与数量控制、穿透控制等控制准则。
在选用单元时要注意剪力自锁、沙漏和网格扭曲、不可压缩材料的体积自锁等问题。
典型有限元软件平台都提供网格映射划分和自由适应划分的策略。
映射划分(Mapped/IsoMesh)用于曲线、曲面、实体的网格划分方法,可使用三角形、四边形、四面体、五面体和六面体,通过指定单元边长、网格数量等参数对网格进行严格控制,映射划分只用于规则的几何图素,对于裁剪曲面或者空间自由曲面等复杂几何体则难以控制。
自由网格划分(Free/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单元,Beam3具备平动方向的三个自由度,而SHELL41包括两个平动自由度(UX/UY)和一个旋转自由度(RTOTZ),因此SHELL41只能传递BEAM3的平动位移,不能传递旋转方向的值。
(3)ANSYS中三维梁单元与三维壳单元具有相同的六个自由度:壳单元旋转自由度与平面旋转刚度相关,为虚拟刚度,不是真实的自由度,同时,要注意三维梁单元与壳单元出现不匹配的问题。