刚塑性有限元数值模拟中产生误差的原因及改进方法
- 格式:docx
- 大小:93.52 KB
- 文档页数:8
有限元仿真误差有限元仿真是一种常用的工程分析方法,它通过将物体划分为小的元素,对每个元素的特性进行计算,然后将这些元素组合成大的模型,以模拟物体在所受力的情况下的应力和位移。
然而,尽管有限元仿真成为了工程领域中广泛使用的工具,但是在实际应用中,由于多种因素的影响,仿真结果可能会产生误差。
以下是一些可能会导致仿真误差的因素。
第一个因素是模型几何形状的简化。
在有限元仿真中,模型的几何形状可能需要进行简化以便进行运算。
例如,将曲面转化成平面或是将立体形状简化成二维形状,并将每个形状分为小元素以便进行计算。
但是,这些简化方式可能会丢失真实模型的细节和复杂性,从而影响仿真结果的精确度。
第二个因素是材料的物理性质的误差。
在进行有限元仿真时,需要求得物体的材料性质,如弹性模量和泊松比等。
这些物理量通常从实验结果中估计而来,但是实验结果往往受到多种因素干扰,如温度、湿度和应力等。
这些因素会影响材料行为的精确度,从而影响仿真结果的准确性。
第三个因素是边界条件的误差。
在有限元仿真中,外部力和边界条件被认为是恒定的,并被认为是物体的一部分。
但在实际场景中,这些条件可能不是恒定的,它们可能会在时间和空间上变化。
例如,在风力作用的情况下,风速和方向可能会随时间而变化,而在地震等自然灾害中,由于地震波在空间中传递,物体的边界条件可能会随着时间和空间的变化而发生变化。
这些外部因素的变化可能会影响物体的行为,并导致仿真结果的误差。
第四个因素是数值计算误差。
由于有限元计算是基于数学算法进行的,所以数值计算误差也可能导致仿真结果的误差。
常见的数值计算误差包括离散化误差、舍入误差和迭代误差等。
在有限元仿真中,如果不采用正确的方法进行模型构建、材料属性计算、边界条件设定和数值计算,那么可能会导致仿真误差的产生。
因此,为了减小仿真误差的影响,需要采取以下方法:首先,通过对模型进行更精细的几何形状建模和材料特性测量以提高数据的准确性和精度,尽量减小模型简化造成的误差。
试验和有限元的误差全文共四篇示例,供您参考第一篇示例:试验和有限元分析是工程领域常用的两种方法,它们常常用于预测和分析结构在不同载荷条件下的响应。
无论是试验还是有限元分析,都存在着误差,因此了解和评估这些误差是非常重要的。
本文将探讨试验和有限元分析中的误差,以及如何有效地管理和减小这些误差。
让我们来看看试验中存在的误差。
试验通常涉及到测量物理量,如应力、应变、位移等。
由于测量设备的精度、环境条件、人为操作等因素,测量结果往往会存在一定的误差。
测量设备的刻度可能不够精确,环境温度和湿度可能会影响到测量结果的准确性,操作人员的技术水平也会对测量结果产生影响。
试验中还可能会出现一些偶然误差,如设备故障、实验样品的缺陷等。
这些偶然误差在一定程度上会影响试验结果的准确性。
对于试验中可能存在的误差,我们需要采取相应的措施来减小这些误差的影响。
比如说,可以通过校准测量设备、控制实验环境、提高操作技术来减小误差,并且在试验结果分析时考虑到可能的误差范围,以便更准确地评估结构的响应。
与试验不同,有限元分析是一种数值计算方法,它通过将结构分割成有限个小单元,利用数学方程对这些小单元进行求解,从而得到结构的响应。
有限元分析中也存在着误差。
有限元分析中的误差可以来自模型的简化。
由于实际结构往往非常复杂,我们在进行有限元建模时往往需要对结构进行简化,例如忽略一些小的细节,这样会导致模型与实际结构存在一定的差异,从而引入误差。
有限元分析中的误差还可能来自数值计算的方法和参数选择。
数值计算方法的选取、边界条件的处理、网格划分的精度等因素都会对有限元分析结果的精度产生影响。
在进行有限元分析时,需要认真选择合适的数值计算方法,合理处理边界条件,以及进行网格收敛性分析,以减小这些误差的影响。
有限元分析中还可能存在由于数值计算误差引起的问题。
使用有限元方法进行求解时,使用的数值积分、迭代收敛条件等都可能会引入数值计算误差,从而影响到结果的准确性。
有限元模型如何查错【一】How to find errors In finite-element models在建立有限元模型的过程 中很容易出错,如果你知 道如何查错,修正这些错 误将会变得很简单PAUL KUROWSKI President ACOM Consulting London, Ontario, Canada BARNA SZABO Professor of Mechanics Washington University St. Louis, Mo.翻译文献出自:MACHINE DESIGN SEPTEMBER 25,1997感谢: 冲令狐,walter2003 编辑sugar,erin 两位网友又见飘整理 August 27, 2006中华钢结构论坛 【正文部分】:有限元分析的第一步就是建立被分析对象的数学模型,这要求我们思索建模 的理论基础如弹性理论,板的 Reissner 理论,塑性变形理论等,和考虑问题的其 它信息如几何描述、材料特性,约束和荷载等等。
分析的目的就是由这些条件, 计算得到精确解 u _ EX 并同时得到位移 u _ EX 的应力函数 F ( u _ EX ) 如 Von Mises 应力等。
应力函数 F ( u _ EX ) 仅仅依赖于数学 模型的定义,而与求解该数学问题的数值近似计算方法无关;同时应力函数F ( u _ EX ) 也不依赖于网格划分、网格类型和单元尺寸。
函数 F ( u _ EX ) 与模型实体物理性质之间的差异,被称为“模型错误” 。
下一步就是使用有限元方法去找到精确解 u _ EX 的近似值 u _ FE 。
这个过程 包括选择网格划分和构件类型,如对二维板用八节点(矩形)单元,依此类推。
网格划分&单元定义被称为有限元的离散化。
离散化产生的误差可以被定义为:e=F (u _ EX ) − F (u _ FE ) F (u _ EX )大部分的分析应该把这个误差控制在 10%以内。
试验和有限元的误差试验是科学研究中常用的一种方法,通过实际操作来验证理论的正确性。
然而,试验结果往往会受到各种误差的影响,这些误差可能来自实验设备的精度、环境条件的变化以及操作人员的技术水平等因素。
对于试验结果的误差,科学家们一直在努力去理解和控制。
有限元方法是一种数值计算方法,广泛应用于工程领域中的结构分析和流体力学等问题的求解。
然而,即使使用有限元方法进行计算,得到的结果也不是完全准确的,会存在一定的误差。
这些误差可能来自离散化网格的选择、数值积分的近似以及模型参数的估计等方面。
试验和有限元方法都有其自身的优势和局限性。
试验能够直接观测到真实系统的行为,具有很高的可信度。
然而,试验过程中的误差会影响结果的准确性和可重复性。
而有限元方法则能够通过数学模型对系统进行建模,具有较高的计算效率和灵活性。
但由于模型的简化和假设,有限元计算结果与真实系统可能存在一定的差异。
为了减小试验和有限元计算的误差,科学家们采取了一系列的措施。
在试验中,他们会提高实验设备的精度和稳定性,控制环境条件的变化,并进行多次重复实验来获取可靠的结果。
在有限元计算中,他们会优化网格的选择,采用更精确的数值积分方法,并进行灵敏度分析和误差估计等工作。
尽管如此,试验和有限元方法的误差仍然无法完全消除。
在实际工程应用中,科学家和工程师们需要根据具体情况,综合考虑试验结果和有限元计算结果,做出合理的判断和决策。
同时,他们也在不断探索新的方法和技术,以提高试验和有限元计算的准确性和可靠性。
试验和有限元方法是科学研究和工程实践中常用的手段,但它们都存在一定的误差。
我们需要认识到误差的存在,并采取相应的措施来减小误差的影响。
只有这样,我们才能更好地理解和解决实际问题,推动科学技术的发展。
本周热点:有限元计算误差的影响因素有限元作为一种数值计算方法,它的计算结果一般与真实解存在误差,影响这些误差的因素有那些?如何减小误差?何种情况下不存在误差(不考虑由于计算机本身的计算误差)?我发表一下个人的一些想法,请各位指正,有限元仿真的结果基本上和真实解都会存在误差的,可从多个方面来说。
1.就是在有限元模拟的时候,我们都要对模型进行一些简化,这一定或多或少影响计算精度的;2.有限元求解的时候,由于各个项目的差异,我们定义各种参数(和实际的一定有差异)例如滑动摩擦系数的值等等,这也会影响理论公式的计算精度;3.建立有限元模型的时候网格的划分,熟练人员和不熟练人员的网格划分有很大差别,这更是影响着求解的计算精度;4.有限元求解本身就是近似计算,它用近似模型替代实际模型,所以计算的最终结果一定和实际存在着一定的差别;5.即使有限元的计算结果正好等于实际值,但是有的实际解在实际中根本没办法测量或者说即使测量了由于采取的手段的诧异,它的结果也不一定非常的精确,这样来说实际的解本身也存在误差;¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥我就考虑到这么多,请各位多多指正。
至于减小误差,我个人认为,这是一种经验的积累,随着我们资历的加深,对分析所采用的各种手段(采用什么样的网格?材料模型?各种参数控制?等等的一些)理解的更加透彻,计算精度一定会更加的精确,由于有限元算法和程序不是我们这些CAE操作人员所能决定的,所以这里对软件本身就不给意见了。
有限元计算中的误差来源与处理方法在有限元计算中,用户最关心自己计算的结果是否准确与合理,那么今天就和大家谈谈有限元计算的误差来源,按照误差来源类型主要分为两类:1、有限元理论假设引入的误差有限元这种数值计算方法,为了实现对现实问题的计算,引入一些力学假设,即(1)连续介质假设,认为计算模型是理想连续,没有孔洞,即位移具有单值性,但是实际产品在加工中难以避免会引入初始的孔洞缺陷,例如铸造件。
(2)材料均匀性假设,认为模型的材料参数不随空间变化为变化,是均匀的,但是实际产品在工艺处理过程中,例如淬火,都会使得材料的力学性能发生改变,并不能严格满足材料均匀性。
有限元理论与真实物理世界的差别通过一个系数来进行折中,也就是有限元计算得到一个基准结果,再通过实验进行对比,计算得到二者的转换系数。
2.有限元计算过程的误差有限元的计算过程主要包括前处理,求解和后处理三个阶段,有限元计算过程的误差主要发生在前处理和求解阶段。
第一个误差来源,即模型简化,发生在前处理阶段模型简化的是否恰当直接影响的计算结果,一般情况数值计算模型与产品的加工CAD模型还是有区别的,数值计算模型只要求把模型的主要特征反映出来,而舍去一下细致的特征,从而平衡了求解精度与计算效率。
第二个误差来源,即材料参数,就是反映计算模型的材料参数,如弹性模量,泊松比等,但是如前所述,即使找到材料参数也不能完全反应由于工艺造成的材料非均匀性,如果输入了错误的材料参数,则计算得到结果没有价值。
第三个误差来源,即工况对接,所谓工况对接,就是在软件中设置的位移约束和载荷与模型实际工况的对应情况,软件只提供了有限的位移约束和载荷类型,而分析模型可能收到的约束和载荷非常多,如果二者对应不正确,则计算结果也没有价值。
第四个误差来源,即网格划分,这个大家容易理解,也是目前很多文献提到的比较多的一个误差,要想将该误差降到最低,必须通过合理的网格加密得到网格无关解。
3.总结与处理有限元理论假设引入的误差是有限元理论与生俱来的,无法避免只能通过一些实验数据进行修正,这也说明要想有限元对工程设计产生作用,必须有与实验对比的过程。
有限元计算误差的影响因素1.网格划分网格划分是有限元方法中最关键的一步,网格的划分对计算结果具有很大的影响。
当网格划分不够细致时,会导致网格近似真实物理结构的能力较差,从而引入较大的误差。
而当网格划分过于细致时,会增加计算量,造成不必要的计算误差。
因此,网格划分需要根据具体问题的特点进行合理选择。
2.材料参数有限元方法在计算中需要使用材料的本构模型和材料的物理性质等参数。
如果这些参数的值与真实材料参数相差较大,就会引入较大的误差。
因此,确定准确的材料参数对于减小有限元计算误差非常重要。
3.边界条件边界条件是指在计算区域内界面及周边所给出的条件。
边界条件的选择和给定不准确都会对计算结果产生很大影响。
合理选择边界条件是保证计算结果准确性的关键。
4.计算方法和算法不同的有限元计算方法和算法对计算结果的准确性也有影响。
例如高阶元素和低阶元素、隐式算法和显式算法等的选择都会对计算误差产生影响。
5.近似假设有限元方法在对实际问题进行数值计算时,通常要对问题进行简化和近似处理。
这些简化和近似假设可能会导致误差的产生。
因此,在进行有限元计算时需要对问题的简化和近似假设进行合理的评估。
6.数值积分在有限元分析中,求解离散形式的形式方程通常需要进行数值积分。
数值积分是将连续函数在一个有限区间中近似表示为离散点的加权和。
数值积分的精度和稳定性会直接影响到计算结果的准确性。
7.迭代收敛有限元求解器通常会使用迭代算法来求解非线性和时间依赖问题。
迭代算法的收敛速度和稳定性对计算误差也会有一定影响。
8.舍入误差总结起来,有限元计算误差的影响因素包括网格划分、材料参数、边界条件、计算方法和算法、近似假设、数值积分、迭代收敛和舍入误差等。
在进行有限元计算前,需要认真评估这些影响因素,并采取相应的措施来减小计算误差,以获得准确可靠的计算结果。
刚塑性有限元数值模拟中产生误差的原因及改进方法1 引言塑性加工过程的有限元数值模拟,可以获得金属变形的详细规律,如网格变形、速度场、应力和应变场的分布规律,以及载荷-行程曲线。
通过对模拟结果的可视化分析,可以在现有的模具设计上预测金属的流动规律,包括缺陷的产生(如角部充不满、折叠、回流和断裂等)。
利用得到的力边界条件对模具进行结构分析,从而改进模具设计,提高模具设计的合理性和模具的使用寿命,减少模具重新试制的次数。
在制造技术飞速发展、市场竞争日益加剧的今天,塑性加工过程的计算机模拟可在模具虚拟设计、制造阶段就能充分检验模具设计的合理性,减少新产品模具的开发研制时间,对用户需求做出快速响应,提高市场竞争能力。
由此可见,金属成型过程的有限元模拟已是模具计算机集成制造系统中必不可少的模具设计检验环节。
金属成形工艺分体积成形和板料成形两大类,相应地,用于分析其流动规律的有限元法也分为两类,即:刚塑性、刚粘塑性有限元和弹塑性有限元。
体积成形中的挤压成形和锻造成形在实际生产中应用很广,中外学者在这方面进行了很多研究,其中二维模拟技术已相当成熟,三维模拟是目前的世界研究热点。
刚塑性、刚粘塑性有限元模拟能否对模具设计的合理性做出可靠校验,取决于模拟的精度和效率。
作者结合从事二维塑性有限元模拟的经验和当前的三维塑性有限元模拟系统开发的实践,对刚塑性、刚粘塑性有限元模拟过程中产生误差的原因进行了全面的详细分析,并提出相应的解决方法,同时以具体实例说明。
2 刚塑性、刚粘塑性有限元模拟中产生误差的原因及改进方法2.1 刚塑性有限元法求解的数学基础刚塑性有限元法是假设材料具有刚塑性的特点,把实际的加工过程定义为边值问题,从刚塑性材料的变分原理或上界定理出发,接有限元模式把能耗率表示为节点速度的非线性函数,利用数学上的最优化原理,在给定变形体某些表面的力边界条件和速度边界条件的情况下,求满足平衡方程、本构方程和体积不变条件的速度场和应力场。
刚塑性有限元数值模拟中产生误差的原因及改进方法作者:陈军邓清张卫刚1 引言塑性加工过程的有限元数值模拟,可以获得金属变形的详细规律,如网格变形、速度场、应力和应变场的分布规律,以及载荷-行程曲线。
通过对模拟结果的可视化分析,可以在现有的模具设计上预测金属的流动规律,包括缺陷的产生(如角部充不满、折叠、回流和断裂等)。
利用得到的力边界条件对模具进行结构分析,从而改进模具设计,提高模具设计的合理性和模具的使用寿命,减少模具重新试制的次数。
在制造技术飞速发展、市场竞争日益加剧的今天,塑性加工过程的计算机模拟可在模具虚拟设计、制造阶段就能充分检验模具设计的合理性,减少新产品模具的开发研制时间,对用户需求做出快速响应,提高市场竞争能力。
由此可见,金属成型过程的有限元模拟已是模具计算机集成制造系统中必不可少的模具设计检验环节。
金属成形工艺分体积成形和板料成形两大类,相应地,用于分析其流动规律的有限元法也分为两类,即:刚塑性、刚粘塑性有限元和弹塑性有限元。
体积成形中的挤压成形和锻造成形在实际生产中应用很广,中外学者在这方面进行了很多研究,其中二维模拟技术已相当成熟,三维模拟是目前的世界研究热点。
刚塑性、刚粘塑性有限元模拟能否对模具设计的合理性做出可靠校验,取决于模拟的精度和效率。
作者结合从事二维塑性有限元模拟的经验和当前的三维塑性有限元模拟系统开发的实践,对刚塑性、刚粘塑性有限元模拟过程中产生误差的原因进行了全面的详细分析,并提出相应的解决方法,同时以具体实例说明。
2 刚塑性、刚粘塑性有限元模拟中产生误差的原因及改进方法2.1 刚塑性有限元法求解的数学基础刚塑性有限元法是假设材料具有刚塑性的特点,把实际的加工过程定义为边值问题,从刚塑性材料的变分原理或上界定理出发,接有限元模式把能耗率表示为节点速度的非线性函数,利用数学上的最优化原理,在给定变形体某些表面的力边界条件和速度边界条件的情况下,求满足平衡方程、本构方程和体积不变条件的速度场和应力场。
速度场的真实解使以动可容速度场建立的能量泛函取极小值。
但所得到的塑性力学的微分方程组一般不能用解析法求解,常采用数值解近似,而采用数值解,则会出现各种误差。
误差取决于所用的数值方法。
下述处理方式易引起系统误差。
2.1.1时间和空间的离散化刚塑性有限元分析的对象是一个非线性变化过程,即材料应力-应变关系的非线性和几何边界条件的非线性。
解决这一问题可以采用线性小变形拟合非线性大变形,如图1所示。
每一个小变形过程的选取须足够小,同时兼顾逼近的精度和效率。
对于刚塑性材料来说,每个加载步长△S,即(△t·Vdz)应小于某一规定值(坯料当前高度的1.0%)[1]。
作者认为,三维模拟的位移加载步长不应超过边界单元最小边长的1/4,以减缓接触边界非平面性的程度,更好的模拟金属的流动规律。
另一方面,所分析的变形材料是一个空间连续体,而有限元法的思想是把无限的连续用有限的连续近似,即用网格离散变形体。
变形场量在单元内连接,这时将产生离散误差,一般地,有限元网格划分得越细,引起的离散误差越小。
但是刚塑性有限元分析必须同时兼顾精度和效率,所以单元不可能过于细化,采用局部网格细分可以满足要求。
如图2所示是作者对方坯反挤工艺三维刚塑性有限元模拟时采用局部网格细分后的网格变形。
文献[2~4]中采用自适应网格离散变形体。
但对于一定的变形过程,即使网格划分得再细,误差仍然存在,这就是形函数误差,也叫做插值误差,是指单元的形函数反映真实变形规律的能力。
所以高阶单元反映材料的变形能力较好,但计算效率偏低;如果同时考虑计算精度、效率以及可视化问题,二维问题宜采用四节点四面体等参元,三维问题宜采用八节点六面体等参元。
三维问题考虑到网格再划分的方便,可采用十节点四面体等参元[5],但可视性较差,文献[6]中采用多种单元类型的处理方法,可以减小一定的插值误差,但系统实现繁杂。
2.1.2解析式的数值计算对于刚塑性有限元分析,单元的刚度矩阵由下式求出:式中[B]——形函数矩阵上式无法直接积分求出,必须进行插值,最有效的方法是高斯插值。
二维问题可采取形心和另外四个高斯积分点插值,三维问题采取形心和八个高斯积分点插值,误差极小。
例如三维问题采用8个高斯积分点插值,上模速度若设为1.0,与27个高斯积分点相比,节点速度的绝对误差小于10-5。
与单刚的求法相似,摩擦力边界条件对单刚和载荷列阵的贡献可采用辛普森插值的方法,二维问题的线性边界取5个积分点,三维问题的类似平面边界中的两个局部坐标方向各取5个积分点,即可满足精度。
在采用计算机进行数值计算时,还会产生截断误差和舍入误差,误差的积累可能产生病态的线性方程,无法获得精确解。
应根据采用的计算机硬件和编程语言,合理选择数值精度。
2.2 刚塑性有限元法求解的力学基础刚塑性有限元求解时,假定材料各向同性,而且体积不变,即忽略弹性变形。
实际上,各种原材料由于生产方式的不同和材料成分的不均匀,多呈现一定程度的各相异性,一般在分析轧制工艺和板材成形时考虑材料的各向异性,而对于其它体积成形问题多不考虑,这种做法仅产生较小的偶然误差。
另一方面,因为在刚塑性有限元分析时,处处体积不变的条件不易满足,常采用三种方法实现近似满足,即Lagrangian乘子法、罚函数法和泊松比接近于0.5法。
第一种方法对力的求解最精确,但线性方程组的求解量大,后一种方法的求解的精度取决于材料的体积可压缩率与实际的接近程度,常用于分析特殊成型过程,如粉末烧结成形。
而罚函数法由于其求解效率高而应用最广。
但这种方法是用惩罚因子与单元平均等效应变速率的乘积近似作为单元的静水应力,这也是产生误差的主要原因。
当模具的塑性有限元法是根据变形能量的泛函进行变形场量的求解,其实质上是一种稍微精确的上限法,上限法所引起的误差是存在的。
刚塑性有限元模拟时,一般不考虑体积力和惯性力,前者引起的误差足够小;对于惯性力,当材料高速成形时,能量泛函中需考虑惯性力做功,此时材料在模腔内的充填规律与低速成形时不大相同,尤其对于刚粘塑性材料[8]。
温度条件的简化也是产生误差的原因。
对于冷挤压问题,从一开始材料接触模具边界的成分相当多,剧烈的摩擦容易引起局部温度的升高,变形载荷和应力分布也发生相应的改变。
因此耦合模拟是必要的。
对于热锻成形过程,尤其对于多工位成形,温度变化相当大,也需要耦合模拟。
刚塑性有限元模拟时,材料的应力-应变函数关系(多采用各相同性强化模式,而极少采用随动强化模式)是一项极其重要的初始条件。
这一具体的函数关系式必须要由专门的实验确定,它的精确性决定了对实际成形过程模拟的近似程度。
文献[9]中的实验表明,变形力和应力的值具有与应力-应变曲线同等的精度,流动特性对应力-应变曲线的偏差则不敏感,该文还尤其强调了应力-应变曲线斜率的精确性要求。
2.3 物理模型在模拟过程中的技术处理刚塑性有限元理论在具体的系统实现时有诸多技术问题。
这些问题的处理直接决定了模拟系统的精度和效率。
若处理不恰当,则会产生误差。
产生误差的原因主要有以下几方面:1)初始速度场的生成刚塑性有限元模拟开始要生成精度合理的初始速度场。
生成方法有多种[10],实践证明,适用于二维和三维任意边界以及速率敏感材料的应属有限元法线性化本构关系和直接迭代生成初始速度场。
三维问题中由触节点的局部坐标方向和坯料放置不当而生成的速度场可能无法用于加载迭代运算,尤其在多工步成形时2)摩擦力边界条件的施加在模具的作用下发生塑性变形的金属与模具表面之间存在着剧烈的摩擦,这在数值模拟过程中,是一项很重要的边界条件。
该边界条件的简化直接影响模拟系统的可靠性。
目前有多种摩擦力的数学模型,应用最为广泛的是C.C.Chen和S.Kobayashi 提出的反正切函数模型[11],可以有效地处理各种情况的边界,尤其是具有分流点的属流动,其模型表达式为:由于上式中α的取值极小,所以一般边界自由节点一旦与模具边界接触,摩擦力几乎是最大值(m·k)。
这可能与实际存在误差,较完善的方法是附加一项考虑相对滑动速度的修正系数。
另外,三维模拟时若采用六面体八节点等参单元,单元与模具的接触部分是四个节点组成的直纹面(即四个节点可能不在同一个平面内),构造与该直纹面对应的虚拟接触平面(关键是其法线方向)至关重要。
3)收敛准则速度场迭代是否收敛常采用的判定准则有三种:节点相对速度误差泛数收敛,节点相对力误差泛数收敛和一阶功率泛函收敛。
如果衰减因子β选取合理(收敛因子的自动调整是决定计算效率最重要的因素),则第一个条件最先满足。
应当指出的是,收敛精度不能取的太低,否则容易引起以后加载步速度场迭代的发散。
同时,对边界接触模具节点的脱模法向力的判断也不准确。
4)刚性区的处理刚塑性有限元法将变形体的弹性变形区视为刚性区,这种区分的准则是极限应变速率ε0。
当单元的平均等效速率小于极限应变速率时,则认为该单元为刚塑性区。
刚性区变形功率的处理有以下两种方法:对于不同的锻压工艺过程的模拟,极限应变速率的取值有所不同。
一般对于变形成形过程,初始变形时刚性区较大,边界约束较小,极限应变速率的取值稍大一点,为10-2~10-3,稍后可逐渐减小极限应变速率的值;而对于挤压变形过程,边界约束较多,极限应变速率的取值可以稍小一点,为10-3~10-4。
如果极限应变速率的取值不合适,则速度场迭代难以收敛,甚至得到完全不同的速度场,模拟的经验就显得特别重要。
文献[12]中提出一种新的不必区分刚塑性区的处理方法,可以减少刚塑性区的区分而引起的误差。
5)动态边界的自动处理由于刚塑性有限元模拟的是金属的整个变形过程(非稳态成形过程),因此网格形状的动态刷新是模拟得以进行下去的保证,动态边界的自动处理正是实现这一功能的前提条件。
这一处理技术主要包括:①边界自由节点接触模具时间步长增一的确定;②时间加载步长增量确定后的位置刷新;③边界触模节点位置的调整;④边界触模节点的脱模判断。
时间加载步长增量若以下面的式子确定:式中△ti——边界自由节点接触模具的时间加载步长增量△t1——保证速度场迭代不发散的时间加载步长增量△t2——保证体积损失不超过1.0%的时间加载步长增量则可以保证较高的位置刷新精度,但计算效率低,有时在模拟模具的局部形状比较特殊的成型过程时,时间加载步长增量经常为零,出现迭代过程中的死点。
为了避免这种情况,可以采用下式确定时间加载步长增量:这样即可以保证刷新的精度,又可以提高模拟的效率。
作者在进行连杆锻造过程滚挤预成型工步的三维刚塑性有限元模拟时,在初始的加载步内,时间加载步长增量的确定先采用式(5),随后再采用式(6),实践证明效率较高。
图3为连杆多工位锻造成型中的预成形滚挤工艺三维仿真的变形网格图。