显示动力学与隐式样动力学区别
显式与隐式算法区别
显式(explicit)和隐式(implicit)这两个词在有限元分析中大家可能经常看到,特别是涉及到动力学分析时。
但其实广义的说他们分别对应着两种不同的算法:显式算法(explicit method)和隐式算法(implicit method)。
所以不论在动力学或者静力学中都有涉及到。
显式算法:不直接求解切线刚度,不进行平衡迭代,计算速度快,时间步长只需要足够小,一般不存在收敛问题,需要的内存也小。
隐式算法:每一增量步都需要对静态方程进行平衡迭代,且每次迭代需要求解大量的线性方程组,这一特点使之占用大量的资源。
但该算法增量步可以很大,至少比显式算法大的多,实际计算中会受到迭代次数及非线性程度的影响我们都知道有限元分析FEA在计算微分方程(differential equations)时,由于计算本身的局限,比如计算机储存的位数有限,以及方程本身的复杂性,计算机运用的是数值算法(numerical algorithm)来逼近真实解的。
有限元分析中数值算法的基础是欧拉法(Euler method),欧拉法又分为forward Euler method 和backward Euler method,这两种方法被简称为显式法(explicit method)和隐式法(implicit method)。
中心差分法:(动力学分析)用有限差分代替位移对时间的求导,将运动方程中的速度与加速度用位移的某种组合来标示,这样就将常微分方程组的求解问题转化为代数方程组的求解问题,并假设在每个小的时间间隔内满足运动方程。
首先我们来看看这两种算法的区别。
显式算法(explicit method )(forward Euler method )考虑常微分方程:初始条件:设为每一步的时间步长, 在Tn 时刻,. (n=0,1,2,3...),在T(n+1)时刻有:所以在显式算法中,T(n+1)时刻的值由T(n)时刻决定,也就是说当前时刻的值由上一时刻的值决定。
隐式求解与显式求解.
大多数非线性动力学问题一般多是采用显式求解方法,特别是在求解大型结构的瞬时高度非线性问题时,显示求解方法有明显的优越性。
下面先简要对比一下隐式求解法和显示求解法。
动态问题涉及到时间域的数值积分方法问题。
在80年代中期以前,人们基本上采用纽曼法进行时间域的积分。
根据纽曼法,位移、速度和加速度有着如下关系:u(i+1)=u(i)+△t*v(i)[(1—2p)a(i)+2p*a(i+1)] (1)v(i+1)=V(i)+△t[(1-2q)a(i)+2qa(i+1)] (2)上面式子中 u(i+1),u(i)分别为当前时刻和前一时刻的位移,v(i+1)和V(i)为当前时刻和前一时刻的速度,a(i+1)和a(i)为当前时刻和前一时刻的加速度,p和q为两个待定参数,△t为当前时刻与前一时刻的时问差,符号 * 为乘号。
由式(1)和式(2)可知,在纽曼法中任一时刻的位移、速度、加速度都相互关联,这就使得运动方程的求解变成一系列相互关联的非线性方程的求解,这个求解过程必须通过迭代和求解联立方程组才能实现。
这就是通常所说的隐式求解法。
隐式求解法可能遇到两个问题。
一是迭代过程不一定收敛,二是联立方程组可能出现病态而无确定的解。
隐式求解法最大的优点是它具有无条件稳定性,即时间步长可以任意大。
如果采用中心差分法来进行动态问题的时域积分,则有如下位移、速度和加速度关系式:u(i+1)=2u(i)-u(i-1)+a(i)(△t)^2 (3)v (i+1)=[u (i+1)-u (i-1)]/2(△t) (4)式中u(i-1),为i -1时刻的位移。
由式(3)可以看出,当前时刻的位移只与前一时刻的加速度和位移有关,这就意味着当前时刻的位移求解无需迭代过程。
另外,只要将运动过程中的质量矩阵和阻尼矩阵对角化,前一时刻的加速度求解无需解联立方程组,从而使问题大大简化,这就是所谓的显式求解法。
显式求解法的优点是它既没有收敛性问题,也不需要求解联立方程组,其缺点是时间步长受到数值积分稳定性的限制,不能超过系统的临界时间步长。
显式算法与隐式算法的区别
显式算法与隐式算法的区别(总4页) --本页仅作为文档封面,使用时请直接删除即可----内页可以根据需求调整合适字体及大小--显式算法与隐式算法的区别所谓显式和隐式,是指求解方法的不同,即数学上的出发点不一样。
并不是说显式只能求动力学问题,隐式只能求静力学问题,只是求解策略不通。
显式求解是对时间进行差分,不存在迭代和收敛问题,最小时间步取决于最小单元的尺寸。
过多和过小的时间步往往导致求解时间非常漫长,但总能给出一个计算结果。
解题费用非常昂贵。
因此在建模划分网格时要非常注意。
隐式求解和时间无关,采用的是牛顿迭代法(线性问题就直接求解线性代数方程组),因此存在一个迭代收敛问题,不收敛就的不到结果。
两者求解问题所耗时间的长短理论上无法比较。
实际应用中一般感觉来说显式耗时多些。
由于两者解题的出发点,所以一般来说显式用于求解和时间相关的动力学问题。
隐式用来求解和时间无关的静力学问题。
但也不是绝对的。
比如,用隐式求解时,为了克服迭代不收敛,改用显式算,但是要多给点时间,这样虽然克服了不收敛的问题,但是求解的时间费用也是相当客观的。
另外,隐式也可以求解动力学问题。
这是ansys里面的两种求解方法。
大多数非线性动力学问题一般多是采用显式求解方法,特别是在求解大型结构的瞬时高度非线性问题时,显示求解方法有明显的优越性。
下面先简要对比一下隐式求解法和显示求解法。
动态问题涉及到时间域的数值积分方法问题。
在80年代中期以前,人们基本上采用纽曼法进行时间域的积分。
根据纽曼法,位移、速度和加速度有着如下关系:u(i+1)=u(i)+△t*v(i)[(1—2p)a(i)+2p*a(i+1)] (1)v(i+1)=V(i)+△t[(1-2q)a(i)+2qa(i+1)] (2)上面式子中 u(i+1),u(i)分别为当前时刻和前一时刻的位移,v(i+1)和V(i)为当前时刻和前一时刻的速度,a(i+1)和a(i)为当前时刻和前一时刻的加速度,p和q为两个待定参数,△t为当前时刻与前一时刻的时问差,符号 * 为乘号。
显式有限元和隐式有限元
显式算法和隐式算法,有时也称为显式解法和隐式解法,是计算力学中常见的两个概念,但是它们并没有普遍认可的定义,下面只是我的一些个人理解。
一、两种算法的比较1、显式算法基于动力学方程,因此无需迭代;而静态隐式算法基于虚功原理,一般需要迭代计算。
显式算法,最大优点是有较好的稳定性。
动态显式算法采用动力学方程的一些差分格式(如广泛使用的中心差分法、线性加速度法、Newmark法和wilson法等),不用直接求解切线刚度,不需要进行平衡迭代,计算速度快,时间步长只要取的足够小,一般不存在收敛性问题。
因此需要的内存也比隐式算法要少。
并且数值计算过程可以很容易地进行并行计算,程序编制也相对简单。
但显式算法要求质量矩阵为对角矩阵,而且只有在单元积分点计算尽可能少时速度优势才能发挥, 因而往往采用减缩积分方法,容易激发沙漏模式,影响应力和应变的计算精度。
静态显式法基于率形式的平衡方程组与Euler向前差分法,不需要迭代求解。
由于平衡方程式仅在率形式上得到满足,所以得出的结果会慢慢偏离正确值。
为了减少相关误差,必须每步使用很小的增量。
除了欧拉向前差分法外,其它的差分格式都是隐式的方法,需要求解线性方程组。
2、隐式算法隐式算法中,在每一增量步内都需要对静态平衡方程进行迭代求解,并且每次迭代都需要求解大型的线性方程组,这以过程需要占用相当数量的计算资源、磁盘空间和内存。
该算法中的增量步可以比较大,至少可以比显式算法大得多,但是实际运算中上要受到迭代次数及非线性程度的限制,需要取一个合理值。
二、求解时间使用显式方法,计算成本消耗与单元数量成正比,并且大致与最小单元的尺寸成反比,应用隐式方法,经验表明对于许多问题的计算成本大致与自由度数目的平方成正比,因此如果网格是相对均匀的,随着模型尺寸的增长,显式方法表明比隐式方法更加节省计算成本。
三、两种方法的应用范围:a)在求解动力学问题时,将方程在空间上采用有限元法(或其他方法)进行离散后,变为常微分方程组F=M(u)+C(u)+K(u)。
Ansys显示算法和隐式算法知识完全解读
An sys 显示算法和隐式算法知识完全解读这是 ansys 里面的两种求解方法。
大多数非线性动力学问题一般多是采用显式求解方法, 特别是在求解大型结构的瞬时高度非 线性问题时, 显示求解方法有明显的优越性。
下面先简要对比一下隐式求解法和显示求解法。
动态问题涉及到时间域的数值积分方法问题。
在 80 年代中期以前,人们基本上采用纽曼法 进行时间域的积分。
根据纽曼法,位移、速度和加速度有着如下关系:u(i+1)=u(i)+ △ t*v(i)[(1 — 2p)a(i)+2p*a(i+1)]⑴ v(i+1)=V(i)+A t[(1-2q)a(i)+2qa(i+1)](2) 上面式子中u(i+1),u(i)分别为当前时刻和前一时刻的位移, 一时刻的速度,a(i+1)和a(i)为当前时刻和前一时刻的加速度, 为当前时刻与前一时刻的时问差,符号 * 为乘号。
由式 (1)和式 (2)可知,在纽曼法中任一时 刻的位移、 速度、 加速度都相互关联, 这就使得运动方程的求解变成一系列相互关联的非线性方程的求解, 这个求解过程必须通过迭代和求解联立方程组才能实现。
能出现病态而无确定的解。
隐式求解法最大的优点是它具有无条件稳定性, 任意大。
如果采用中心差分法来进行动态问题的时域积分, 则有如下位移、 速度和加速度关系式:u(i+1)=2u(i)-u(i-1)+a(i)( △ t)A 2(3) v(i+1)=[u(i+1)-u(i-1)] / 2(△ t) (4)式中 u(i-1) ,为 i-1 时刻的位移。
由式 (3)可以看出,当前时刻的位移只与前一时刻的加速度 和位移有关, 这就意味着当前时刻的位移求解无需迭代过程。
另外, 只要将运动过程中的质 量矩阵和阻尼矩阵对角化, 前一时刻的加速度求解无需解联立方程组, 从而使问题大大简化,这就是所谓的显式求解法。
显式求解法的优点是它既没有收敛性问题,也不需要求解联立方 程组,其缺点是时间步长受到数值积分稳定性的限制,不能超过系统的临界时间步长。
AbaqusExplicit和AbaqusStandard比较
AbaqusExplicit和AbaqusStandard比较非线性动力学Abaqus/Explicit和Abaqus/Standard比较Abaqus/Standard没有稳定时间增量的限制:完成指定分析,一般需要较少的时间增量。
因为每个增量步都需要求解线性方程组,所以每个增量步的费用非常昂贵。
Abaqus/Explicit有稳定时间增量的限制:完成指定分析,一般需要较多的时间增量。
因为不需要求解线性方程组,所以每个增量步的费用较低。
–所有的计算费用与单元计算相关(形成和组装I )。
计算时:Abaqus/Standard适合于模拟与模型的振动频率相比研究响应周期较长的问题。
由于时间增量的限制,使用显式动力学求解效率很低。
用于具有适度非线性的问题,其中非线性是平滑的(比如,塑性)。
如果是平滑的非线性响应,Abaqus/Standard需要较少次数的迭代找到收敛解。
Abaqus/Explicit适合于模拟高速动力学问题–需要较少的时间增量;隐式动力学效率较低。
对于包含不连续的非线性问题,一般更加可靠。
–接触行为是不连续的,并且包含碰撞。
对于隐式时间积分,接触和碰撞将产生求解问题。
–其它不连续行为包括屈曲和材料失效。
非线性动力学例子显式动力学过程概述显示动力学求解器与隐式求解器,比如Abaqus/Standard,是互为补充的。
从用户的角度出发,隐式与显式方法显著的区别为:显式方法需要小的时间增量。
只与模型的最高自然频率相关。
与载荷类型和载荷持续时间无关。
一般的,增量步的数量级为10,000到1,000,000个增量,但是每个增量步内的计算费用相对较小。
隐式方法对于时间增量没有内在的限制。
增量大小一般由精度和收敛性控制。
同显式方法相比,隐式模拟中的增量步成数量级的少。
然而,因为必须在每个增量步内求解系统总体方程组,隐式方法在每个增量步中的费用比显式方法大很多。
知道两种方法的特征之后,在实际分析过程中就可以为所分析的问题选择合适的方法。
显示算法和隐式算法单点积分算法和全积分算法
先说说显示算法和隐式算法:这是ansys里面的两种求解方法。
大多数非线性动力学问题一般多是采用显式求解方法,特别是在求解大型结构的瞬时高度非线性问题时,显示求解方法有明显的优越性。
下面先简要对比一下隐式求解法和显示求解法。
动态问题涉及到时间域的数值积分方法问题。
在80年代中期以前,人们基本上采用纽曼法进行时间域的积分。
根据纽曼法,位移、速度和加速度有着如下关系:u(i+1)=u(i)+△t*v(i)[(1—2p)a(i)+2p*a(i+1)] (1)v(i+1)=V(i)+△t[(1-2q)a(i)+2qa(i+1)] (2)上面式子中 u(i+1),u(i)分别为当前时刻和前一时刻的位移,v(i+1)和V(i)为当前时刻和前一时刻的速度,a(i+1)和a(i)为当前时刻和前一时刻的加速度,p和q为两个待定参数,△t为当前时刻与前一时刻的时问差,符号* 为乘号。
由式(1)和式(2)可知,在纽曼法中任一时刻的位移、速度、加速度都相互关联,这就使得运动方程的求解变成一系列相互关联的非线性方程的求解,这个求解过程必须通过迭代和求解联立方程组才能实现。
这就是通常所说的隐式求解法。
隐式求解法可能遇到两个问题。
一是迭代过程不一定收敛,二是联立方程组可能出现病态而无确定的解。
隐式求解法最大的优点是它具有无条件稳定性,即时间步长可以任意大。
如果采用中心差分法来进行动态问题的时域积分,则有如下位移、速度和加速度关系式:u(i+1)=2u(i)-u(i-1)+a(i)(△t)^2 (3)v(i+1)=[u(i+1)-u(i-1)]/2(△t) (4)式中u(i-1),为i-1时刻的位移。
由式(3)可以看出,当前时刻的位移只与前一时刻的加速度和位移有关,这就意味着当前时刻的位移求解无需迭代过程。
另外,只要将运动过程中的质量矩阵和阻尼矩阵对角化,前一时刻的加速度求解无需解联立方程组,从而使问题大大简化,这就是所谓的显式求解法。
显示算法与隐式算法
以下内容转自abaqus版面的总结:显式一般用于动态问题的分析, 对于大型问题, 或复杂的接触情况可能需要几百万的增量步的计算, 所用时间可能是几天或更长. 而隐式的增量步长要长得多, 一般用于静态问题的求解.显式算法别explicit method use direct iterative method, which has small cost in eachtime increment but require relatively small increment. Abaqus pre-determinethe time increment based on wave propagation speed and minimum meshsize. This method could be efficient for highly nonlinear and contact problem.For quasi-static problem, properly adjust model parameter as density and totaltime is important to achieve good computation time.standard-隐式算法Implicit method use newton method for iteration, which means high cost foreach time increment but could mean large time increment. Convergencecould be a problem in this case. It could be efficient for linear and some nonlinear problem. More materials, elements and procedures are available in standard.所谓显式和隐式,是指求解方法的不同,即数学上的出发点不一样。
显式算法与隐式算法的区别
显式算法与隐式算法得区别1、显式算法最大优点就是有较好得稳定性。
动态显式算法采用动力学方程得一些差分格式(如广泛使用得中心差分法、线性加速度法、Newmark法与wilson法等),不用直接求解切线刚度,不需要进行平衡迭代,计算速度快,时间步长只要取得足够小,一般不存在收敛性问题。
因此需要得内存也比隐式算法要少。
并且数值计算过程可以很容易地进行并行计算,程序编制也相对简单。
但显式算法要求质量矩阵为对角矩阵,而且只有在单元级计算尽可能少时速度优势才能发挥。
因而往往采用减缩积分方法,容易激发沙漏模式,影响应力与应变得计算精度。
静态显式法基于率形式得平衡方程组与Euler向前差分法,不需要迭代求解。
由于平衡方程式仅在率形式上得到满足,所以得出得结果会慢慢偏离正确值。
为了减少相关误差,必须每步使用很小得增量。
2、隐式算法隐式算法中,在每一增量步内都需要对静态平衡方程进行迭代求解,并且每次迭代都需要求解大型得线性方程组,这以过程需要占用相当数量得计算资源、磁盘空间与内存。
该算法中得增量步可以比较大,至少可以比显式算法大得多,但就是实际运算中上要受到迭代次数及非线性程度得限制,需要取一个合理值。
3、求解时间t使用显式方法,计算成本消耗与单元数量成正比,并且大致与最小单元得尺寸成反比;应用隐式方法,经验表明对于许多问题得计算成本大致与自由度数目得平方成正比;因此如果网格就是相对均匀得,随着模型尺寸得增长,显式方法表明比隐式方法更加节省计算成本。
所谓显式与隐式,就是指求解方法得不同,即数学上得出发点不一样。
并不就是说显式只能求动力学问题,隐式只能求静力学问题,只就是求解策略不通。
显式求解就是对时间进行差分,不存在迭代与收敛问题,最小时间步取决于最小单元得尺寸。
过多与过小得时间步往往导致求解时间非常漫长,但总能给出一个计算结果。
解题费用非常昂贵。
因此在建模划分网格时要非常注意。
隐式求解与时间无关,采用得就是牛顿迭代法(线性问题就直接求解线性代数方程组),因此存在一个迭代收敛问题,不收敛就得不到结果。
显示动力学与隐式样动力学区别
显式求解是对时间进行差分,不存在迭代和收敛问题,最小时间步取决于最小单元的尺寸。
过多和过小的时间步往往导致求解时间非常漫长,但总能给出一个计算结果。
解题费用非常昂贵。
因此在建模划分网格时要非常注意。
隐式求解和时间无关,采用的是牛顿迭代法(线性问题就直接求解线性代数方程组),因此存在一个迭代收敛问题,不收敛就的不到结果。
两者求解问题所耗时间的长短理论上无法比较。
实际应用中一般感觉来说显式耗时多些。
由于两者解题的出发点,所以一般来说显式用于求解和时间相关的动力学问题。
隐式用来求解和时间无关的静力学问题。
但也不是绝对的。
比如,用隐式求解时,为了克服迭代不收敛,改用显式算,但是要多给点时间,这样虽然克服了不收敛的问题,但是求解的时间费用也是相当客观的。
另外,隐式也可以求解动力学问题。
1、显式算法基于动力学方程,因此无需迭代;而静态隐式算法基于虚功原理,一般需要迭代计算2、显式算法最大优点是有较好的稳定性。
动态显式算法采用动力学方程的一些差分格式(如广泛使用的中心差分法、线性加速度法、Newmark法和wilson法等),不用直接求解切线刚度,不需要进行平衡迭代,计算速度快,时间步长只要取的足够小,一般不存在收敛性问题。
因此需要的内存也比隐式算法要少。
并且数值计算过程可以很容易地进行并行计算,程序编制也相对简单。
但显式算法要求质量矩阵为对角矩阵,而且只有在单元级计算尽可能少时速度优势才能发挥, 因而往往采用减缩积分方法,容易激发沙漏模式,影响应力和应变的计算精度。
静态显式法基于率形式的平衡方程组与Euler向前差分法,不需要迭代求解。
由于平衡方程式仅在率形式上得到满足,所以得出的结果会慢慢偏离正确值。
为了减少相关误差,必须每步使用很小的增量。
3、隐式算法隐式算法中,在每一增量步内都需要对静态平衡方程进行迭代求解,并且每次迭代都需要求解大型的线性方程组,这个过程需要占用相当数量的计算资源、磁盘空间和内存。
该算法中的增量步可以比较大,至少可以比显式算法大得多,但是实际运算中上要受到迭代次数及非线性程度的限制,需要取一个合理值。
LSDYNA显示动力学简介
2
显式动力学的应用
应力波 传播 接触 碰撞
制造 工艺 耐撞性 分析
3 隐式动力学简介
cx kx F (t ) mx
对于特定积分参量的线性问题 对于非线性问题
无条件稳定
可采用较大的时间步
线性逼近,需要多次迭代
收敛需要较小时间步
小结
无须进行收敛检查
收敛
时间步
谢谢
LS/DYNA 显式算法
n1/2 i
x
n1/2 i
x t
n
t f * h / c min
P F H Mx
n n n
n
中心差分法
P:为外力和体力的合力 ∆t为时间增量, f为时间步稳定因子 n时刻加速度→n+1/2时刻各个方向速度 →n+1时刻坐标值 F:为应力放散矢量 H:为沙漏阻力 h为单元的特征尺寸,c为局部应力波波速
无须进行迭代 迭代
时间步大小受CFL条件限制
Explicit Dynamics
LS/DYNA 显式动力学
目录
显式动力学简介 显式算法 隐式动力学简介
1
2
显式动力学的应用
3
4隐式算法Biblioteka 小结1 显式动力学简介
x x
n1 i
x x
n i
n1/2 i n i
t
n1/2
n+1时刻的坐标值 n+1/2时刻各个方向速度 时间步CFL条件限制
隐式与显式动力学的区别-弹性动力学有限元基本解法
1.1. 弹性动力学有限元基本解法结构系统的通用运动学方程为:tR KU U C U M =++ (1) 求解该动力学振动响应主要有三类方法:(1)时域法(2)频域法(3)响应谱法时域法又可分为:(1)直接积分法,(2)模态叠加法。
直接积分法又可分为中心差分法(显式),Wilson θ(隐式)法以及Newmark (隐式)法等。
本文介绍中心差分法(显式)与Newmark (隐式)法。
1 中心差分法(显式)假定0,1t ,2t ,…,n t 时刻的节点位移,速度与加速度均为已知,现求解)(t t t n ∆+时刻的 结构响应。
中心差分法对加速度,速度的导数采用中心差分代替,即为:)2(12t t t t t t U U U t U ∆+∆-+-∆= )(21t t t t t U U tU ∆-∆+-∆= (2) 将(2)式代入(1)式后整理得到tt t R U M ˆˆ=∆+ (3) 式(3)中C tM t M ∆+∆=211ˆ2 t t t t t U C tM t U M t K R R ∆-∆-∆-∆--=)211()2(ˆ22 分别称为有效质量矩阵,有效载荷矢量。
R ,M ,C ,K 为结构载荷,质量,阻尼,刚度矩阵。
求解线性方程组(3),即可获得t t ∆+时刻的节点位移向量t t U ∆+,将t t U ∆+代回几何方程与物理方程,可得t t ∆+时刻的单元应力和应变。
中心差分法在求解t t ∆+瞬时的位移t t U ∆+时,只需t t ∆+时刻以前的状态变量t U 和t t U ∆-,然后计算出有效质量矩阵M ˆ,有效载荷矢量tR ˆ,即可求出t t U ∆+,故称此解法为显式算法。
中心差分法,在开始计算时,需要仔细处理。
t =0时,要计算t U ∆,需要知道t U ∆-的值。
因此应该有一个起始技术,因而该算法不是自动起步的。
由于0U ,0U ,0U 是已知的,由t =0时的(2)式可知: 02002U t U t U U t ∆+∆-=∆-中心差分法中时间步长t ∆的选择涉及两个方面的约束:数值算法的稳定性和计算时间。
显式隐式求解
<显式动力学&隐式动力学分析>1、显式算法基于动力学方程,因此无需迭代;而静态隐式算法基于虚功原理,一般需要迭代计算2、显式算法最大优点是有较好的稳定性。
动态显式算法采用动力学方程的一些差分格式,不用直接求解切线刚度,不需要进行平衡迭代,计算速度快,时间步长只要取的足够小,一般不存在收敛性问题。
因此需要的内存也比隐式算法要少。
并且数值计算过程可以很容易地进行并行计算,程序编制也相对简单。
但显式算法要求质量矩阵为对角矩阵,而且只有在单元级计算尽可能少时速度优势才能发挥, 因而往往采用减缩积分方法,容易激发沙漏模式,影响应力和应变的计算精度。
静态显式法基于率形式的平衡方程组与Euler 向前差分法,不需要迭代求解。
由于平衡方程式仅在率形式上得到满足,所以得出的结果会慢慢偏离正确值。
为了减少相关误差,必须每步使用很小的增量。
3、隐式算法隐式算法中,在每一增量步内都需要对静态平衡方程进行迭代求解,并且每次迭代都需要求解大型的线性方程组,这个过程需要占用相当数量的计算资源、磁盘空间和内存。
该算法中的增量步可以比较大,至少可以比显式算法大得多,但是实际运算中上要受到迭代次数及非线性程度的限制,需要取一个合理值。
4、求解时间使用显式方法,计算成本消耗与单元数量成正比,并且大致与最小单元的尺寸成反比; 应用隐式方法,经验表明对于许多问题的计算成本大致与自由度数目的平方成正比;因此如果网格是相对均匀的,随着模型尺寸的增长,显式方法表明比隐式方法更加节省计算成本隐式算法将冲压成型过程的计算作为动态问题来处理后,就涉及到时间域的数值积分方法问题。
在80年代中期以前,人们基本上使用纽曼法进行时间域的积分。
根据纽曼法,位移、速度和加速度有着如下的关系:[][](1)()()(12)()2(1)(1)()(12)()2(1)u i u i t v i p a i pa i v i v i t q a i qa i +=+∆⨯-+++=+∆-++上面式子中, u(i+1)和u(i) 分别为当前时刻和前一时刻的位移, v(i+1)和v(i)为当前时刻和前一时刻的速度, a(i+1)和a(i)为当前时刻和前一时刻的加速度,p 和q 为两个待定参数。
abaqus显式方法和隐式方法
abaqus是一个常用的有限元分析软件,广泛应用于工程领域。
在abaqus中,有两种常用的有限元分析方法,分别是显式方法和隐式方法。
这两种方法在不同的工程场景下有着不同的适用性和特点,本文将对abaqus显式方法和隐式方法进行介绍和比较。
一、abaqus显式方法1. 简介abaqus显式方法又称为动力显式有限元法,适用于求解瞬态非线性动力学问题。
其基本原理是通过将时间步长划分为相对较小的子步长,在每个子步长内计算结构的响应。
显式方法的特点是计算速度快,适用于求解高速碰撞、爆炸等暂态动力学问题。
2. 适用场景abaqus显式方法适用于以下工程场景:- 高速动态加载问题:如车辆碰撞、子弹撞击等。
- 爆炸冲击问题:如爆破工程、军事领域的冲击载荷等。
- 材料动态行为问题:如弹性-塑性材料在高速加载下的动态响应。
3. 优缺点优点:- 计算速度快,适用于大规模复杂结构的动力学仿真。
- 对动态非线性效应的处理能力强。
- 需要选择合适的时间步长和子步长,计算稳定性受到限制。
- 适用范围受到限制,不能用于稳态问题的分析。
二、abaqus隐式方法1. 简介abaqus隐式方法又称为静态隐式有限元法,适用于求解稳态和瞬态非线性静力学问题。
其基本原理是通过建立增量形式的非线性方程组,采用迭代算法求解非线性方程组的平衡解。
隐式方法的特点是计算稳定,适用于求解稳态和缓变动力学问题。
2. 适用场景abaqus隐式方法适用于以下工程领域:- 结构强度分析:如钢结构、混凝土结构的承载能力分析。
- 热-机-结构耦合问题:如热载荷下的构件稳定性分析、变温环境下的材料性能分析。
- 持续动态加载问题:如地震、风载等动态荷载下的结构响应分析。
3. 优缺点优点:- 计算稳定性好,适用于求解长时间的稳态和瞬态非线性问题。
- 对非线性效应的处理能力强。
- 计算速度相对较慢,在处理大规模结构的动力学仿真时耗时较长。
- 算法求解过程复杂,需要对非线性迭代收敛条件进行调整。
显式求解方法和隐式求解方法对比
采用有限元方法开展结构的动力学分析最终归结为求解离散后的常微分方程组tR KU U C U M =++ 。
在时域内求解该方程最常用的方法是直接积分法,而又根据求解过程中是否需要迭代求解线性方程组,将直接积分法分为隐式积分方法和显式积分方法两类。
隐式积分法认为t+Δt时刻系统的状态不仅与t时刻状态有关,且与t+Δt时刻某些量有关。
因此隐式算法是根据t n 及t n-1...时刻体系的物理量值建立关于以t n+1时刻物理量为未知量的线性方程组,通过求解方程组确定t n+1时刻的物理量(常用的方法有线性加速度法、常平均加速度法、Newmark 方法、Wilson-θ法、Houbolt 方法等)。
而显式积分法认为t+Δt时刻系统的状态仅与t时刻状态有关可,因此可由t n 及t n-1...时刻体系的物理量值直接外推t n+1时刻物理量值(如中心差分法),不需要求解线性方程组,实现了时间离散的解耦。
两种算法的比较 :(1)隐式算法隐式算法基于虚功原理,要迭代计算。
隐式算法在每一增量步内都需要对静态平衡方程进行迭代求解,并且每次迭代都需要求解大型的线性方程组,这一过程需要占用相当数量的计算资源、磁盘空间和内存。
理论上在这个算法中的增量步可以很大,但是实际运算中上要受到接触以及摩擦等条件的限制。
随着单元数目的增加,计算时间几乎呈平方次增加。
由于需要矩阵求逆以及精确积分,对内存要求很高。
隐式算法的不利方面就是收敛问题不容易解决,且在开始起皱失稳时,在分叉点处刚度矩阵出现奇异。
(2)显式算法显示算法基于动力学方程,无需迭代,包括动态显式和静态显式算法。
动态显式算法采用动力学方程的中心差分格式,不用直接求解切线刚度,不需要进行平衡迭代,计算速度快,也不存在收敛控制问题。
该算法需要的内存也比隐式算法要少。
数值计算过程可以很容易地进行并行计算,程序编制也相对简单。
它也有一些不利方面。
显式算法要求质量矩阵为对角矩阵,而且只有在单元级计算尽可能少时速度优势才能发挥, 因而往往采用减缩积分方法,容易激发沙漏模式,影响应力和应变的计算精度。
显式表达式
显式表达式
显式表达式,可以理解为呈现出来的表达式,或者是显示出来的表达式,相对于隐性表达式。
显式算法基于动力学方程,因此无需迭代;而静态隐式算法基于虚功原理,一般需要迭代计算。
显式算法最大优点是有较好的稳定性。
动态显式算法采用动力学方程的一些差分格式(如广泛使用的中心差分法、线性加速度法、Newmark法和wilson法等),不用直接求解切线刚度,不需要进行平衡迭代,计算速度快,时间步长只要取的足够小,一般不存在收敛性问题。
因此需要的内存也比隐式算法要少。
并且数值计算过程可以很容易地进行并行计算,程序编制也相对简单。
显式与隐式的区别
使用显式方法,计算成本消耗与单元数量成正比,并且大致与最小单元的尺寸成反比,应用隐式方法,经验表明对于许多问题的计算成本大致与自由度数目的平方成正比,因此如果网格是相对均匀的,随着模型尺寸的增长,显式方法表明比隐式方法更加节省计算成本。
显式算法是建立在i时刻的运动平衡方程,不需要迭代,运算简单但是对步长要求很高,因为其影响精度和稳定性;而隐式算法是建立在i+1时刻的,因此需要迭代,过程复杂些,但是更加精确。