首页/文章/ 详情

abaqus中的隐式分析步如何设置参数

25天前浏览357
图片

概述

  隐式分析常常用于长持时、动荷载计算中,在abaqus中的参数设置中,大致可分为两种:自动增量步长和固定增量步长。下面详细介绍这两种计算方式和具体的参数设置,以及合适的场输出设置。

步长的概念

  要解释“步长”的概念,最好从静力计算开始。在静力计算中,离散后的有限元方程为

  理论上来讲,静力计算中,没有时间的概念,只需要对矩阵求一次逆,就能把问题求解出来。但是abaqus计算静力问题也是采用的增量迭代法,abaqus将荷载拆分,然后依次施加在结构上,最终将全部的荷载施加在结构上,每一次的荷载大小和简单的理解为“步长”。

    在动力问题中,

  abaqus采用了更加复杂的增量迭代法,因为动力问题中有了时间,abaqus就把时间离散,在每一个很小的时间间隔内,计算被视为静力计算,但事实上并不是静力计算,仅仅是数学意义上的分段线性近似。同时,在每一个小的时间间隔,计算又被拆分成更小的计算单位。
  那么,步长数值的确定就是一个非常重要的问题,轻则影响计算时间和精度,重则计算不收敛。这就涉及到计算的收敛判断和步长计算,这部分的内容,非常、极度的复杂,abaqus的非线性计算性能能够傲视群雄很大程度上是因为它的收敛与发散判断非常的高效而准确。
  当然,这个帖子并不是探讨abaqus是如何计算收敛与发散的,就算是蚂蚁能把大象吃了我也整不明白这玩意儿。帖子主要将abaqus中关于隐式计算中的参数设置


主页面设置

  在进行主要数据的设置前,abaqus提供了一些前置的参数设置,如下图图片  这里面的Time Period指的是计算的总时长,即荷载作用在结构上的时间,这个务必设置准确,这个参数后面需要用到。
  下面的Nlgeom是几何非线性,这个按钮的打开与否取决于自己的计算涉及几何大变形与否。这个页面其他的设置不太重要,看看即可。

自动增量步长

  自动增量步长常常用于我们对模型能否收敛不是很确定的时候,一般都是把最小增量步长设置的非常的小。
  进入incrementation界面,选择Automatic选项,即自动增量步长,如下图图片  自动增量步长,顾名思义,就是abaqus自动确定计算过程中每一迭代步的增量步长,不需要我们的人工干预,但是需要我们设置一些初始的参数,下面依次介绍这些参数。
  首先是Maximum number of increments,这个参数是我们允许这次计算的最大增量步数,如果超过了这个增量步数量计算还没有到我们在前面设置的计算总时长,abaqus就会直接退出计算,需要注意,这并不是意味着我们的计算失败了,我们只需要把这个数字调大一些就行了。
  initial incrament size是初始的增量步步长,即第一次迭代的步长,一般来说,这个参数最好设置小一点,因为万事开头难,计算机算东西有时候也是这样,一开始就把步长整的很过分,就很难继续计算下去。
  Minimum incrament size是整个计算过程中我们允许的最小增量步长。abaqus在遇到收敛困难的时候会缩小步长,如果一直缩小到我们允许的最小步长计算还是没有收敛,abaqus就会直接退出计算,如果我们把这个参数胡调整的非常的小还是计算不下去,那可能是自己的模型不正确。 其他的参数不是非常的重要,看看即可。

固定增量步长

  固定增量步长常常应用于我们对模型非常有自信,觉得模型一定可以收敛的时候。
  固定增量步长和自动增量步长在同一个页面,如下图所示
图片  这里的参数只有两个需要设置。
  Maximum number of increments是增量步的数量,这个参数是可以计算出来的,前面我们设置了计算总时长,依据下面的increment size就可以确定总的增量步数。假如计算总时长是10s,我们设定的固定增量步长是0.02s,那么总增量步数就是500,这三个数值要对应正确。

场输出设置

  这里为什么突然开始说场输出设置了呢?因为不同的计算方式需要匹配合适的输出频率,主要是因为不同的计算方式有不同的增量步长,如下图所示,abaqus提供了非常多的输出方式图片  这个页面里的内容非常多,本次帖子主要关注Frequency,即输出频率。
  首先是Last increment,这个参数只在每一次分析步的最后一个增量步输出数据,其余计算过程中的增量步数据均不被输出,这种设置不常用。
  Every n increment是经常用的设置,他的意思每n个增量步输出一次数据,在固定增量步长计算中,这个参数常常被设置为1,这样就能得到固定时间间隔的数据,非常有利于后处理,但是自动增量步长就不太适合选择这种设置,因为自动增量步长计算中的步长不是等间隔的,我们后续还需要插值处理,比较麻烦。
  Every spaced time intervals指的是每隔固定时间输出一次数据,这种设置非常适合自动增量步长,无论增量步长是多少,abaqus会自动进行时间插值计算,最终只按照固定的时间间隔输出数据。

Abaqus科普
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2025-01-11
最近编辑:25天前
外太空土豆儿
博士 我们穷极一生,究竟在追寻什么?
获赞 26粉丝 8文章 70课程 0
点赞
收藏
作者推荐

基于非凸多边形单元的光滑有限元法——任意网格构建和热传导分析

一、研究背景和难点多边形单元因其在域离散化过程中的高度灵活性、出色的网格细化及自适应性能力,在数值算法领域取得了显著的进展。然而,非凸单元的应用对数值算法的稳健性构成了严峻的挑战。针对这一难题,本文采用光滑有限元方法(SmoothedFiniteElementMethod,S-FEM)来求解采用非凸单元离散的二维热传导问题,为处理非凸单元的数值算法提供了一种新的视角和解决方案。二、非凸多边形网格构建及问题域离散为了验证S-FEM在处理任意凹网格的鲁棒性,本研究通过图像采点构建了多种具有高度复杂性的凹多边形单元。下图所示即为依据图像信息所构造的一个具代表性的“狗形”凹多边形单元示例,该示例直观展示S-FEM在应对复杂多边形单元挑战时的应用潜力。将上述图像中抽象所得的网格嵌入至标准化的坐标系内,随后在该正方形单元内部构建支撑单元,以确保对整个正方形区域实现一个既无重叠又无遗漏的离散化处理。在此基础上,将此正方形区域视为一个基单元,采用逆坐标映射技术,将其应用于整个问题域的离散化过程中。这一方法允许我们采用任意复杂且非凹的多边形单元来对问题域进行离散,可充分验证数值算法的鲁棒性,灵活性和适用范围。三、光滑有限元算法针对任意复杂形态的多边形单元,S-FEM采用耳切技术来构造光滑积分域。并基于光滑梯度技术,实现了将单元内部的面积分转化为光滑域边界上的线积分。这策略有效地避免了雅可比矩阵的计算,从而显著增强了S-FEM在处理网格畸变问题时的鲁棒性。如上图所示,对于含有n个节点的任意多边形单元,通过耳切技术,可以生成n-2个光滑积分域。根据下述公式,可以验证上述光滑积分域满足S-FEM在求解热传导问题时所需的稳定性条件。 表1光滑有限元求解热问题的稳定性条件维度最少的光滑域数量1D2D3D四、数值算例及验证为了验证本文所提出方法的有效性和适用性,我们设计了基于均匀介质及多材料情形的数值算例。下图展示了一个均匀介质热传导问题的具体算例,其中采用了具有独特几何特征的“狗形”单元和“兔形”单元来对问题域进行离散化处理,充分检验S-FEM在处理复杂几何网格时的性能和准确性。下图呈现了一个涉及多种材料的热传导问题数值算例,旨在验证所提方法在自适应网格的鲁棒性。为此,采用不同网格密度,多种类型的缝合多边形单元对问题域进行离散化处理。五、结论S-FEM不仅在线性单元求解中展现出极高的精度和收敛性,而且在非凸多边形单元的应用上也表现出良好的适用性,本研究不仅扩展了S-FEM的应用范围,使其能够更有效、更灵活的离散问题域。基于梯度光滑技术,S-FEM在处理非凸网格时展现出自稳定性,无需引入任何额外的稳定项。这一特性使得S-FEM在保持与传统有限元方法相同格式的同时,能够更稳健地应对网格畸变和非凸性带来的挑战。六、论文引用格式:WuSW,YangR,NiuRP,CaoLX,LiuX,WanDT,LiuGR,Anovelformulationforheatconductionusingnon-convexmeshesbasedonsmoothedfiniteelementmethod.InternationalJournalofHeatandMassTransfer,2025,237:126310.WuSW,GR.Liu,JiangC,LiuX,LiuK,WanDT,YueJH,Arbitrarypolygonmeshforelasticandelastoplasticanalysisofsolidsusingsmoothedfiniteelementmethod.ComputerMethodsinAppliedMechanicsandEngineering,2023,405:115874七、作者信息吴少伟,博士,现为长沙理工大学汽车与机械工程学院讲师,主要从事光滑有限元理论和算法,非线性计算固体力学、非线性材料本构模型、核-热-力耦合的多物理场软件开发等研究。主持国家自然科学青年基金,湖南省自然科学青年基金2项,以核心成员参与国家自然科学基金重点项目。以第一作者或通讯作者在计算力学权威期刊CMAME、IJNME、CM、IJMS、AMMOD、IJHMT等杂志发表SCI论文9篇。担任多个期刊审稿人。邮箱:wusw@csust.edu.cnResearchgate:https://www.researchgate.net/profile/Wu-Shaowei-3/research来源:有限元先生

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈