首页/文章/ 详情

你已经是一个成熟的工程师了,还不知道这些CAE分析要点和步骤吗?

7月前浏览817

计算机辅助工程(computer aided engineering,CAE)可以在设计阶段对结构进行校核、优化,使工程师在产品未生产之前就对设计的经济性、可靠性、安全性进行评估。

在各种CAE方法中,有限元法(finite element method,FEM)在工程领域应用最广,也是技术相对比较成熟的一种方法。一个合格的有限元分析(finite element ana lysis,FEA)工程师,至少应该具备三方面的要求:

(1)坚实的理论基础,主要包括力学理论(结构仿真工程师)和有限元理论;

(2)软件运用经验,能熟练应用常用的有限元软件;

  (3)  工程实践经验,对于各种工程问题能够准确地判断并确定分析方案。

有限元分析前,工程师应首先针对实际工程问题进行判断,依据工况确定问题类型,判定是否需要有限元分析(很多问题用基本力学计算或者查手册就能更快更准确地得到答案)。

其次, 对有限元分析项目进行规划并计算,包括模型简化及计算规模、分析类型和边界条件的确定。

最后, 根据有限元分析结果,提出相应结论和建议,包括分析项目的可靠性、 安全性判定,优化的可能,危险的处理等。

由上可知,有限元分析工程师的工作是对专业知识及实践经验的综合性体现。那么,我们如何判断在一个问题能否采用有限元分析呢?或者采用有限元分析必须满足什么条件呢?

1.CAE分析需满足的条件

一般而言,一个 CAE 分析必须满足下列四要素:

(1)问题应具有清晰的物理概念,需判定所涉及的工程问题的数学问题类型。一个工程问题,按数学一般分为:①稳定场(椭圆)方程,用于描述静平衡、 稳态热等;②扩散(抛物线)方程,用于描述动力学、瞬态热等;③波动(双曲线)方程,用于描述应力波等波动现。

(2)明确问题系统的所需的各项属性。如已知上面的3种控制方程,还需要有初始条件或/和边界条件,才可以得到方程的解析解。因此,问题系统中需要具备基本的自身参数,如弹性模量、泊松比、长度、 截面积等,还要具备系统的外界参数,如力、力矩等。

(3)工程问题的数学表征或抽象。实际工程问题往往存在于大量的数据中,需要抽取或换算得到数学表征参数或建立数学模型。例如,一对齿轮副进行静平衡(静力学)计算,除了知道其模型尺寸、 材料的弹性模量、泊松比、齿轮间的摩擦系数以外,还需要通过计算求出其载荷(力和力矩),以保证各个齿轮加载后整个系统的力平衡和力矩平衡;如果进行瞬态计算,则需要知道齿轮的密度,载荷以转速度形式加载。 

(4)问题在计算机上实现的可行性和高效性。任何有限元分析都基于一定的假设,例如,连续性 是实现有限元计算的必要条件,各向同性、对称性则是实现有限元高效计算的简化手段。 

此外,有限元法是实际工程设计中一种数学模型求解的辅助方法,需为实际工程而服务,主要解决的 是难以被实验验证的工程问题, 因此切忌为数学分析而分析。

那么作为一名有限元分析工程师,在遇到一个需要有限元分析的实际问题时,该如何入手呢?又应该特别注意什么问题呢?

通常分析实际工程问题,进行有限元分析时需从简化模型、确定边界条件、输入合理的参数和采用合适的网格几方面入手。因此,这也是分析需特别注意几点问题。

2.模型简化

在决定需要进行有限元分析后,对分析的模型及其工况在理论和本质上均要有清晰的认识,避免不合理和不切实际的分析。有限元模型必须与分析目的、计算机性能匹配,并不是模型越精确计算精度越高。

有限元分析的模型越精确就意味着模型越复杂,进而要求软件进行更加复杂的矩阵化简求解。这样一来,模型的误差虽然小了,计算误差反而增大,导致最终得不到合适的结果。因此,在建立有限元模型时,应尽量采用尽可能简单的模型,无需保留实物模型的所有细节特征。

常用作法是:去掉非关键位置 的小孔和槽,用圆孔代替螺纹孔,用直角代替圆角及倒角。 

如果保留实物模型的诸多细微几何特征,会导致分析结果的应力集中(可用应力集中系数K表征应力集中程度,收敛),甚至出现应力奇异状态(表现出网格相关性,不收敛)。如下图所示。

  尺寸变化处应力集中系数K与尺寸相关

 

L形区域的平面问题的转角处应力奇异

对于实体为桁架的模型,整体分析不考虑剪力的影响时可以将实体模型简化为梁模型, 这样总体计算量较小,精度也较高。如果单向剪力较大,则可将实体模型抽壳为壳模型计算。只有需关注桁架局部的详细受力,且双向剪力较大时,才使用实体模型。


 可见, CAE分析必须根据理论和经验判断模型几何细节的相关性,进而确定模型的简化方式。但是,有时一些模型几何细节开始时显得不重要,简化分析后在这些细节处应力较大,则可以在有限元模型中恢复几何细节或采用子模型分析


3.边界条件正确


将实际工程问题转化为力学问题,分析对象的选取、载荷工况和施加载荷的确定、边界条件(位移约束条件)的确定、结构的刚度和质量、载荷传递路径和应力集中等问题的处理是CAE 分析的关键。


有限元分析时,必须按照实际工况的边界条件,且满足有限元平衡方程,才能求解得到正确结果。


在静力学分析中容易出现边界条件不足的问题,容易产生刚体运动。虽然有的软件会自动处理,如ANSYS Workbench 会自动将弱弹簧(Weak Springs)施加到模型可能出现刚体 位移的位置,如下图。但是还是建议工程师设置好充分的约束后,将弱弹 簧设置为 Off

 

ansysworkbench 中弱弹簧处理刚体运动


针对约束不足的正确方法是对模型先进行模态分析,观察是否具有刚体模态(模态分析出的固有频率在 0~1Hz),依据其频率对应的模态形状,进一步分析是否存在刚体运动(单个零件),或者存在零件之间接触不足(组件或部件)。


 接触分析(无摩擦接触、粗糙接触、摩擦接触)涉及迭代计算,如果在接触面体上施加力载荷时,往往难以收敛;改为位移载荷,则相对简单得多。


4.合理的参数设置


进行有限元分析,需要对软件有深刻的认识,做到每输入一个参数都清楚知道这个参数的意义和作用,这不仅仅是需要熟悉软件的界面,更多的是需要理解有限元和力学的理论,有时甚至需要对参数进行一些常识性的辨识。


在进行有限元分析时,合理设置参数至关重要,它直接影响到分析结果的准确性和可靠性。那么如何在有限元分析中合理设置材料参数呢?


在设置参数时,一般需要注意以下几点:


  • 单位制统一确保所有输入的材料参数单位一致,以避免计算错误。很多通用有限元软件中,材料参数的输入没有单位,需要使用者在输入参数时统一单位制,如abaqus、LSDYNA等。如果材料参数单位错误,那么计算结果也必然错误。

  • 材料参数选择通常需要输入的参数包括弹性模量、泊松比和密度。这些参数足以描述材料的基本弹性行为。需要注意的是,采用不同的材料本构模型,即不同目的的分析所需的材料参数是不一样的。比如动力学分析中,必须要有质量或密度参数,传热分析必须要有导热系数等传热参数

  • 泊松比限制自然界几乎所有的材料都为正泊松比材料,金属铝和铜分别为0.133、0.127,典型的聚合物泡沫为0.11~0.14左右,对于不可压缩材料(橡胶类材料)泊松比的数值接近0.5。当泊松比接近0.5时,材料模型会变得不稳定。比如橡胶作为不可压缩材料,具有很强的材料非线性,在计算过程中,往往会因为产生过大的变形导致网格畸变严重(下图),又或者产生体积自锁,而难以收敛,这时需要考虑使用“几乎不可压缩材料”选项。如ansysworkbench中可以使用U-P单元技术解除体积自锁,插入command如下图所示。


   
   

 
   
   
 



此外,在材料参数设置方面,材料参数设置与采用的材料本构模型相关。本构模型(Constitutive Model)是一种数学模型(如混凝土的本构如下图),用于描述材料如何响应不同的物理加载和应力条件。它是材料力学和结构分析领域的重要概念,帮助工程师和研究人员理解和预测材料的行为。


 


本构模型在有限元分析中广泛使用。有限元分析是一种数值方法,用于模拟复杂结构和系统的行为。模拟采用的本构模型决定了材料的行为,因此在仿真分析中材料本构模型中材料参数的合理性非常重要。在实际工程中一般材料参数需要进行标定,以确保仿真模拟的精度。


5.适用的网格划分


目前,很多仿真软件都拥有非常智能的自主划分网格能力,整个计算甚至颠覆了常用有限元软件的流程,不用把划分网格作为一个必要的操作。但采用软件自主划分网格,往往结果是:重要部分(应力集中区、接触区、大变形区)的结果不准确或不收敛; 不重要部分网格过细,浪费了计算时间。


网格的划分往往需要实践经验,当然可以参考软件提供的网格质量进行评估一般情况下,网格划分遵循以下原则:


1. 圆形模型采用古钱币的切分划法,中间正方形的边长为圆形的半径,圆周上至少需 要 40 等分;


2.圆环模型采用多体划分;


3.模型厚度方向至少需要 3 层单元;缩减积分时厚度方向至少划分 4 层单元;


4.接触比较难收敛时,可以采用无中间节点的网格形式(单元类型发生变化)。


 对于一个模型最理想的网格形状是什么?答案就是整个模型都是大小一致的正方形和立方体,当然由于模型存在斜角、圆弧等要素,这个要求很难满足,因此,需要根据形状拓扑的规律尽量满足上述条件。 网格的数量到底多少合适呢?


正确的方法是将网格尺寸定义为参数变量,利用 优化分析方法,确定其敏感度。如果网格细化到一定程度,应力结果不会有大的偏差,就说明网格密度合适。  

 


来源:薛定谔的Cube
WorkbenchAbaqus静力学非线性通用理论材料单元技术控制ANSYS
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2025-12-31
最近编辑:7月前
巴郡撸猫人
硕士 签名征集中
获赞 15粉丝 25文章 73课程 0
点赞
收藏
作者推荐

一看就会,一学就废? 如何用ANSYS Workbench中的梁单元进行分析?

各位小伙伴们大家好呀~,今天咋们来盘一盘ANSYSWorkbench中的梁单元。并根据一个简单的案例分析,说明ANSYSWorkbench中关于梁单元的一些使用要点。首先需要说明的是,本文的主要内容依旧是来源于周炬、苏金英老师所著《ANSYSWorkbench有限元分析实例详解》一书中的案例,经过笔者整理,以供学习参考。那么话不多说,我们开始吧~1梁单元概述小伙伴们肯定知道,当我们的结构杆件是长条形状时,直接采用实体模型进行分析,必将消耗大量计算机资源。所以一般在结构分析时我们会采用梁单元进行分析处理。这个问题我们在:如何快速将3D杆系模型转化为CAE梁单元线体分析模型?这篇文章告诉你正确的方法!这片篇文章讲过,此处不在赘述。在工作中,一般当结构长度对横截面的比率超过10:1时,沿长度方向的应力将会成为主要分析对象,且横截面始终保持不变时,采用梁单元分析将会极大提高我们的效率。通常梁单元可用于分析主要受侧向或横向载荷的结构,如建筑桁架、桥梁、螺栓等。在ANSYSWorkbench中默认为铁木辛柯(Timoshenko)梁单元,即Beam188和Beam189。当然ANSYSWorkbench中还包括许多其他梁单元,如beam3(2D)和beam4(3D)为代表的欧拉梁。本文以Beam188和Beam189为例。铁木辛柯梁是能考虑剪切变形的梁。具体地说,它的位移和截面转角是独立插值的,而不是由位移的导数来求得。当对于梁的高度(或者直径)较大,高跨比增加时,必须考虑剪切变形和转动惯量的梁。相对的,不需要考虑的就是欧拉-伯努力梁。简单地说,铁木辛柯梁是考虑剪切的梁,欧拉梁就是材料力学里讲的梁,忽略剪切作用。欧拉梁基于平截面面假定,弯曲是主要变形,忽略剪切变形的影响,其计算公式通过平衡微分方程得到,而非变形协调方程。因此,在ANSYSWorkbench中默认的铁木辛柯(Timoshenko)梁单元,适用于描述短梁、层合梁以及波长接近厚度的高频激励时梁的表现。因此,ANSYSWorkbench中的梁单元可计算弯曲、轴向、扭转和横向剪切变形,适用性较好。那么ANSYSWorkbench中Beam188和Beam189梁单元有什么区别呢?Beam188采用线性多项式作为形函数,Beam189采用二次多项式作为形函数。当在ANSYSWorkbench(以2020软件为例)软件中的Mesh设置中,Mesh—Elementorder设置为linener(无中节点的线性单元)时,即为Beam188;Mesh-Elementorder设置为quadratic(有中节点的二次单元)时,即为Beam189。如图1所示。图1Beam188和Beam189单元及选择我们知道,有限元对单元特性的描述包括单元形状、节点数目、自由度和形函数等。表1为Beam188和Beam189单元的对比。在ANSYSWorkbench2020中默认为线性单元,即beam188。一般来说,线性单元需要更多的网格数才能达到二次单元的精度。选用二次单元可提高计算精度,这是因为二次单元的曲线或曲面边界能够更好地逼近结构的曲线和曲面边界,且二次插值函数可更高精度地逼近复杂场函数,所以当结构形状不规则、应力分布或变形很复杂时可以选用高阶单元。但高阶单元的节点数较多,在网格数量相同的情况下由高阶单元组成的模型规模要大得多,计算内存消耗也多,因此,在使用时应权衡考虑计算精度和时间。接下来,我们就以一个案例讲一下在ANSYSWorkbench2020中如何采用梁单元分析吧!2梁单元分析案例2.1草图绘制首先,我们在ANSYSWorkbench2020主界面中拖入staticstructural,如下图2所示。打开几何(geometry)下的DM模块。并在草绘界面绘制一边长为30mm、40mm、50mm的三角形,如图3所示。图2选择静力结构分析模块图3绘制skecth草图2.2生成线体和定义截面在concept菜单下,利用“linefromskecth”创建线体,并定义一个梁截面,赋值给线体。分别如下图4~图6所示。需要注意,当各线体的需赋予不同的截面时,不同界截面的线体利用skecth草图绘制时,应绘制在不同的草图上。图4创建线体图5创建截面图6赋予截面2.3调整截面方向和位置在定义截面后,有时我们会发现截面方向不正确,这事我们需要进行调整。首先,在工具选项卡上,切换为线体选择。如图7所示。图7线体选择前,切换选择模式然后,选中需要调整截面方向的线体。再按照详细设置面板的参数进行旋转方向和角度的调整。如下图8所示。最终调整为如图9所示的状态。需要注意,存在多个线体时,最后需用FormNewPart将三根梁合并为一个部件(part)。否则在进入model模块后,就需要进行Connections设置。因为FormNewPart操作,会将体素合并为一个整体,相当于杆件之间为刚接。图8截面方向调整图9最终状态2.4材料单元等设置关闭几何建模,打开model模块。在Geometry中可以看到模型的基本尺寸、质量等参数,其中ElementControl的设置非常重要,在梁模型中一般选择ProgramControlled。图10ElementControl设置点击LineBody,可以在Assignment处设置材料,本例采用默认的StructuralSteel;此外,还可以在Color处设置模型颜色。OffsetType处可选择梁偏移模式,分别为Centroid(默认)、ShearCenter(剪切中心)、Origin(原点)、UserDefined(自定义,注意偏移的X、Y距离不是基于整体坐标系,而是基于单个梁模型的坐标系)。其他设置还可根据需要调整。如下图11所示。图11LineBody下的一些设置值得注意的是,modeltype下可以选择线体的类型,对应了不同的单元类型。如下图12所示。图12modeltype下的设置2.5网格划分与梁单元类型选择Mesh可以不设置,计算时ANSYSWorkbench会自动对模型划分网格,即定义模型为Beam188单元。当然也可手动更改为189单元,如下图13所示。图13更改单元类型2.6边界条件设置梁的边界条件设置如下,并在StaticStructural中设置如下选项:①AnalysisSettings中的StepControls不需要设置,全部默认;WeakSprings(弱弹簧)、LargeDeflection(大变形)、InertiaRelief(减小惯性)均设置Off。这是因为本例为小变形的静力学分析,没有收敛计算。②如图14所示对模型施加边界条件,将选择过滤为选择点,对左上角点施加SimplySupported和FixedRotation约束。SimplySupported施加在梁或壳体的边缘或者端点上,约束平移,但是允许旋转;FixedRotation施加在梁或壳体的表面、边缘或者端点上,约束旋转,但是允许平移。因为梁为6自由度模型,本例如果不施加FixedRotation约束将无法计算。对于梁模型,同时对一点施加SimplySupported和FixedRotation可等效为FixedSupport。③对右上角点施加沿X方向的Force:点击DefineBy选择Components,在XComponent处输入100N。④对左下角点施加Displacement:点击DefineBy选择Components,在XComponent选择Free,YComponent和ZComponent输入0mm。图14约束和荷载施加2.7求解及后处理完成以上设置即可进行求解。求解后,点击BeamTool可以看到DirectStress、MinimumCombinedStress、MaximumCombinedStress等结果。其中:①直接应力(DirectStress)为轴向应力;②最小弯曲应力为四个弯曲应力的最小值;③最大弯曲应力与最小弯曲应力相反;④最小组合应力(MinimumCombinedStress)是直接应力和最小弯曲应力的线性组合,即N+Mmin;⑤最大组合应力(MaximumCombinedStress)是直接应力和最大弯曲应力的线性组合,即N+Mmax;。图15BeamTool下的应力选项应力结果如下:图16BeamTool下的应力结果在ANSYSWorkbench中也可以在DesignModeler中自定义梁截面。对于自定义截面一定不能忽视截面斜角、圆角等细微参数,否则会有较大的计算误差。且如果截面为自定义,后处理里就没有BeamTool,只能在点选对应结果。本案例中,在beamresult下,可以查看梁轴力(AxialForce)、弯矩(BendingMoment)、扭矩(TorsionalMoment)、剪力(ShearForce)、剪力-弯矩图(ShearMomentDiagram)。对于某一横梁,在查看其弯矩剪力图前,需定义路径。以斜边梁为例,定义在其上的路径。建立一个Path,将结果映射到Path上。右键Model→Insert→ConstructionGeometry→Path,然后在Detailsofpath中将pathtype切换为edge,如下所示。图17path定义再点击ShearMomentDiagram(总剪力弯矩图),选择上面定义的路径,可以看到斜边梁各结果图,如图18所示。图18剪力-弯矩图熟悉ANSYS经典界面的小伙伴常在后处理中提取smisc和nmisc数据,ANSYSWorkbench可以通过UserDefinedResult来提取该数据库,并可以进行数学处理。右键点击Solution,在Insert菜单下选择UserDefinedResult,如图19所示。图19用户自定义结果在设置面板中的Expression处可分别定义为smisc1、smisc2、smisc5、smisc6、smisc7、smisc11、smisc12、smisc13、sqrt(ux^2+uy^2+uz^2)。如下图以smisc2为例。图20自定义结果参见Beam188单元表可知,以上自定义输入对应求解分别代表:单元的轴向应力、弯矩、剪应力、轴向应变、剪应变、截面积和总变形(可与TotalDeformation结果对比,两者一致)。需要注意UserDefinedResult后处理输出结果没有单位,读者必须了解各输出类型的物理意义,再按照选择单位制读取结果数据。好啦~,以上就是今天分享的全部内容了。关于梁单元分析的一些其他内容,比如如何分析变截面梁、如何设置铰接、如何作为螺栓分析,咱们后续接着盘~来源:薛定谔的Cube

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