显式算法与隐式算法区别

这篇具有很好参考价值的文章主要介绍了显式算法与隐式算法区别。希望对大家有所帮助。如果存在错误或未考虑完全的地方,请大家不吝赐教,您也可以点击"举报违法"按钮提交疑问。

基本概念

        有限元计算微分方程由于计算机本身的从存储限制及方程的复杂性,采用数值方法来逼近真实解。有限元算法基础——欧拉法、拉格朗日算法,以欧拉法为例,主要分为forward Euler Method及backward Euler Method
    Forward Euler Method——显示算法,也即Explicit Method。
    Backward Euler Method——隐式算法,也即Implicit Method。

算法区别

1.概述

  • 在显式算法中, T ( n + 1 ) T(n+1) T(n+1)时刻的值由 T ( n ) T(n) T(n)时刻决定,也就是说当前时刻的值由上一时刻的值决定。
  • 1、求解常微分方程: d y d t = f ( t , y ) \frac{dy}{dt}=f(t,y) dtdy=f(t,y)
    2、初始条件: y ( t 0 ) = y 0 y(t_0)=y_0 y(t0)=y0
    3、设 h h h为每一步的时间步长,在tn时刻: t n = t 0 + n h t_n=t_0+nh tn=t0+nh
    4、 T ( n + 1 ) T(n+1) T(n+1)时刻数值: y n + 1 = y n + h f ( t n , y n ) y_{n+1}=y_n+hf(t_n,y_n) yn+1=yn+hf(tn,yn)
  • 在隐式算法中, T ( n + 1 ) T(n+1) T(n+1)时刻的值不仅由 T ( n ) T(n) T(n)时刻决定,还由当前时刻 T ( n + 1 ) T(n+1) T(n+1)决定;也就是说当前时刻的值由上一时刻和当前时刻的值共同决定;隐式算法往往需要求解二次方程。
  • 1、求解常微分方程: d y d t = f ( t , y ) \frac{dy}{dt}=f(t,y) dtdy=f(t,y)
    2、初始条件: y ( t 0 ) = y 0 y(t_0)=y_0 y(t0)=y0
    3、设 h h h为每一步的时间步长,在tn时刻: t n = t 0 + n h t_n=t_0+nh tn=t0+nh
    4、 T ( n + 1 ) T(n+1) T(n+1)时刻数值: y n + 1 = y n + h f ( t n + 1 , y n + 1 ) y_{n+1}=y_n+hf(t_{n+1},y_{n+1}) yn+1=yn+hf(tn+1,yn+1)

2.收敛性

  • 显示算法不存在收敛性的问题(因为不进行收敛计算),每个时刻的值由上一时刻所确定,但当时间步取得较大时,就会脱离真实值。
  • 隐式算法是无条件收敛的,在隐式算法中,在求解二次方程的同时,会通过牛顿法对每一步进行迭代收敛,直至收敛到指定的偏差。

3.时间步长的依赖性

        显示算法的准确性严重依赖于计算时间步长,要想获得比较准确的结果,时间步需要非常小,但是显式算法是绝对收敛的,对时间步长要求不高。文章来源地址https://www.toymoban.com/news/detail-670404.html

4.计算效率

  • 显示动力学由于是间接求得位移 x x x的,要取得足够的精度,需要取很短的时间步长,所以需要进行很多步计算,但每一步计算需要的时间很短。
  • 隐式动力学由于直接求得位移 x x x,所以不存在收敛性问题,但由于动力学分析往往涉及非线性(几何非线性、材料非线性等),可能每一步都要进行刚度的求逆,非常耗时,但时间步长相对显示算法可以取得很大。

5.适用范围

  • 显式算法:一般用于作用时间短,非线性程度高的动力学分析、准静态分析。如爆炸、碰撞、跌落及倒塌等。
  • 隐式算法:一般用于作用时间较长、非线性程度不高的动力学、静力学分析及准静态分析。

到了这里,关于显式算法与隐式算法区别的文章就介绍完了。如果您还想了解更多内容,请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章,希望大家以后多多支持TOY模板网!

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处: 如若内容造成侵权/违法违规/事实不符,请点击违法举报进行投诉反馈,一经查实,立即删除!

领支付宝红包 赞助服务器费用

相关文章

  • 利用Matlab实现Abaqus有限元分析预处理器——从inp文件导入模型信息到自编Matlab有限元程序前处理

    大家好,我非常荣幸有这个机会和大家分享我在有限元分析方面的一些经验。我希望我的这篇文章能够帮助你在有限元分析领域获得更多的知识,提升你的技能。这篇文章主要介绍了如何使用Matlab实现一个有限元分析的预处理器,这个预处理器能够从Abaqus的 *.inp 文件中导入模

    2024年02月10日
    浏览(48)
  • 有限元(FEM)基本知识速阅

    1 什么是有限元 2. 固体力学的偏微分方程 密度 位移 3.本构方程 由弹性模型和泊松比就能确定 本构矩阵 进而确定应力和应变的关系 将含有9各变量的微分方程 变为 u v w 三个待求函数的 微分方程 需要进一步加入边界条件 才能求微分方程 4.边界条件 5.CAD模型与微分方程的关系

    2024年02月15日
    浏览(39)
  • python实现简单的二维有限元计算

    有限元算法依据常见的有限元法教材,简单复现悬臂梁在重力作用下的形变(为了变形更明显,重力大小扩大了10倍),还没来得及写注释。【卧槽快跑,没注释!】 节点是随机函数撒的点,完全没有优化; meshpy库中的Delauny优化算法计算得到三角单元; pygame实现图形绘制,

    2024年02月06日
    浏览(45)
  • 有限元三角形单元的等效节点力

    写等几何程序的时候用到等效节点力,之前没有好好理解等效节点力这一块,现在补充学习一下。首先是三角形单元的等效节点力: 可以结合之前的博客《平面问题有限元》一起进行理解。 可以查看知网文章 :《关于面积坐标在三角形单元上的积分》         下面是推导

    2024年02月04日
    浏览(50)
  • 143基于matlab的2D平面桁架有限元分析

    基于matlab的2D平面桁架有限元分析,可以改变材料参数,输出平面结构外形,各桁架应力,位移及作用力。可查看节点力,程序已调通,可直接运行。 143 matlab 平面桁架 有限元分析 桁架应力 (xiaohongshu.com)

    2024年01月25日
    浏览(44)
  • Matlab-杆单元整体刚度矩阵组装(有限元基础-曾攀)

    一维杆单元的组装: 二维杆单元组装 三维情况下以此类推。

    2024年02月16日
    浏览(42)
  • ANSYS APDL 输出有限元模型刚度矩阵和质量矩阵

    APDL输出刚度矩阵和质量矩阵的命令流代码,后附matlab处理代码 根据以上代码定义,将输出刚度矩阵到‘matkMMF.txt’,质量矩阵到‘matmMMF.txt’,如下图 其中每一行的前两个数值代表在矩阵中的行和列号,第三个数值为在该位置的矩阵元素数值。如‘1 1 7.536000000000000E-01’代表

    2024年02月08日
    浏览(44)
  • 【小呆的力学笔记】非线性有限元的初步认识【二】

    1.2 有限元分析的数学原理 书接上回,我们已经回顾了线性有限元分析的理论基础——线弹性力学的内容,虽然有限元方法本质上等效于弹性力学方程的数值解法,但是有限元方法毕竟不是差分法,它并不直接离散控制方程,而是通过等效于控制方程的另一种形式来实现的。

    2024年02月01日
    浏览(85)
  • 【小呆的力学笔记】非线性有限元的初步认识【一】

    有限元方法是仿真模拟的重要方法,特别是在结构仿真中占据主导地位。从本质上来说,有限元方法和有限差分法等方法一样,都是一种求解物理方程的数值计算方法。计算机技术的迅猛发展使得数值仿真应用越来越广泛,有限元方法作为重要仿真工具内容也越来越丰富。

    2023年04月08日
    浏览(42)
  • Unity-数字孪生-实时显示大型结构体有限元分析的应变云图

    分享针对大型结构体数字孪生的技术思路给到大家,给大家提供一些灵感,促进交流。 1、这段时间优化和改进了Abaqus模型转Unity的流程,完成从超50W结点的大型装配体Abaqus inp文件自动提取各部件的结点、面片信息,并在Unity当中进行重构。 2、成功走通Unity-MySQL实时显示应变

    2024年03月22日
    浏览(104)

觉得文章有用就打赏一下文章作者

支付宝扫一扫打赏

博客赞助

微信扫一扫打赏

请作者喝杯咖啡吧~博客赞助

支付宝扫一扫领取红包,优惠每天领

二维码1

领取红包

二维码2

领红包