COMSOL Multiphysics中2D非线性磁场有限元仿真

1算例介绍

一电磁铁模型截面及几何尺寸如图1所示，铁芯为软铁，磁化曲线(B-H)曲线如图2所示，励磁电流密度J=250 A/cm2。现需分析磁铁内的磁场分布。

图1电磁铁模型截面图(单位cm)

图2铁芯磁化曲线

2 COMSOL Multiphysics仿真步骤

根据磁场计算原理，结合算例特点，在COMSOL Multiphysics中实现仿真。

(1)设定物理场

COMSOL Multiphysics 4.0以上的版本中，在AC/DC模块下自定义有8种应用模式，分别为：静电场(es)、电流(es)、电流-壳(ecs)、磁场(mf)、磁场和电场(mef)、带电粒子追踪(cpt)、电路(cir)、磁场-无电流(mfnc)。其中，“磁场(mef)”是以磁矢势A作为因变量，可应用于：

①已知电流分布的DC线圈；

②电流趋于表面的高频AC线圈；

③任意时变电流下的电场和磁场分布；

根据所要解决的问题的特点——分析磁铁在线圈通电情况下的电磁场分布，选择2维“磁场(mf)”应用模式，稳态求解类型。

(2)建立几何模型

根据图1，在COMSOL Multiphysics中建立等比例的几何模型，如图3所示。

图3几何模型

有限元仿真是针对封闭区域，因此在磁铁外添加空气域，包围磁铁。

由于磁铁的磁导率，因此空气域的外轮廓线可以理想地认为与磁场线迹线重合，并设为磁位的参考点，即

(21)

式中，L为空气外边界。

(3)设置分析条件

①材料属性

本算例中涉及到的材料有空气和磁铁，在软件自带的材料库中选取Air和Soft Iron。

对于磁铁的B-H曲线，在该节点下将已定义的离散B-H曲线表单导入其中即可。

②边界条件

由于磁铁的磁导率，因此空气域的外轮廓线可以理想地认为与磁场线迹线重合，并设为磁位的参考点，即

(21)

式中，L为空气外边界。

为引入磁铁的B-H曲线，除在材料属性节点下导入B-H表单之外，还需在“磁场(mef)”节点下选择“安培定律”，域为“2”，即磁铁区域，在“磁场>本构关系”处将本构关系选择为“H-B曲线”。此时，即表示将材料性质表达为磁通密度B的函数，也符合以磁矢势A作为因变量时的表达，从而避免在本构关系中定义循环变量。设置窗口如下图所示。

图4磁铁本构关系设置

该模型中，线圈中励磁电流密度为J=250 A/cm2，因此，在“磁场(m ef)”节点下，选择两个“外部电流密度”节点，分别用于设置两个线圈的电流密度。根据式(2)，该电流密度因为z 轴方向的电流密度，且两个线圈的电流密度方向应相反。事先在模型树下定义参数J，表达式为“250e4*A/cm^2+”。设置窗口如下图所示。

(a)线圈1

(b)线圈1

图5线圈电流设置

(4)网格划分

网格节点下直接创建三角形网格，结果如下图所示。

图6网格划分结果

(5)求解

选择“稳态”求解模式，直接进行计算。该模型结构比较简单，求解时间为2 s。4后处理

磁通密度如图7所示。

图7面磁通密度分布

由图6可以看出，磁通密度主要集中分布在磁铁上，在转角处磁通密度较大(图中红色区域)；在空气域磁通密度很小。

磁通密度等值分布图如下所示：

图8磁通密度等值线分布

磁通密度方向如下图所示，图中箭头表示磁通密度方向。

图9磁通密度等值线分布

磁力线分布如下图所示：

图10磁势分布

由图10可看出，磁矢势A围绕线圈，在磁铁中形成闭合曲线，图中线的密度形象表示了磁场的强弱，在转角处的线较稠密。

2 结论

本文在COMSOL Multiphysics中实现了对一简单电磁铁模型磁场分析，并有以下结论：

(1) COMSOL Multiphysics中，引入B-H曲线数据时，需指定本构关系中的设置，避免在本构关系中出现循环变量；同时，需另在该域上定义安培定律；

(2) COMSOL Multiphysics中，等值线所表示的和流线所表示的意义不同；等值线是将值相等的点连接而成，而流线只是表示方向，即线上某点切线表示该点的实际方向；图10中用等值线表示磁矢势A z，即是磁力线；

非线性有限元方法及实例分析

非线性有限元方法及实例分析梁军河海大学水利水电工程学院，南京（210098）摘要：对在地下工程稳定性分析中常用的非线性方程组的求解方法进行研究，讨论了非线性计算的迭代收敛准则，并利用非线性有限元方法分析了一个钢棒单轴拉伸的实例。关键词：非线性有限元，方程组求解，实例分析 1引言有限单元法已成为一种强有力的数值解法来解决工程中遇到的大量问题，其应用范围从固体到流体，从静力到动力，从力学问题到非力学问题。有限元的线性分析已经设计工具被广泛采用。但对于绝大多数水利工程中遇到的实际问题如地下洞室等，将其作为非线性问题加以考虑更符合实际情况。根据产生非线性的原因，非线性问题主要有3种类型[1]： 1．材料非线性问题（简称材料非线性或物理非线性） 2．几何非线性问题 3．接触非线性问题（简称接触非线性或边界非线性） 2 非线性方程组的求解在非线性力学中，无论是哪一类非线性问题，经过有限元离散后，它们都归结为求解一个非线性代数方程组[2]： ()()()00 021212211=… …==n n n n δδδψδδδψδδδψΛΛΛ （1.1）其中n δδδ,,,21Λ是未知量，n ψψψ,,,21Λ是n δδδ,,,21Λ的非线性函数，引用矢量记号 []T n δδδδΛ21＝（1.2） []T n ψψψψΛ21= （1.3）上述方程组（1.1）可表示为 ()0=δψ （1.4）可以将它改写为 ()()()0=?≡?≡R K R F δδδδψ （1.5）其中()δK 是一个的矩阵，其元素是矢量的函数，n n ×ij k R 为已知矢量。在位移有限元中，δ代表未知的结点位移，()δF 是等效结点力，R 为等效结点荷载，方程()0=δψ表示结点平衡方程。在线弹性有限元中，线性方程组

非线性有限元分析

轨道结构的非线性有限元分析姜建华练松良摘要实际轨道结构受载时的力学行为,属于典型的非线性力学问题。钢轨垫层刚度、钢轨抗扭刚度和扣件扣压力的大小是影响轨距扩大的主要因素。根据非线性有限元接触理论,建立了能准确反映扣件、钢轨与垫层的拧紧接触,以及受载车轮与钢轨侧向滑动接触的力学计算模型;并研究计算了不同扣件压力下,由于受载车轮与钢轨侧向滑动接触引起的轨距扩大问题。关键词轮轨关系,扣件压力,非线性弹性力学,有限元分析 1 引言实际工程中常见的非线性问题一般可以归纳为三类:材料非线性、几何非线性以及边界条件非线性。材料非线性问题是由于材料的非线性本构关系所引起的,例如材料的弹塑性变形,材料的屈服和硬化等;几何非线性问题是由于结构的位移或变形相当大,以至必须按照变形后的几何位置来建立平衡方程;边界条件非线性问题是指边界条件随位移变化所引起的非线性问题。通常情况下,我们所遇到的非线性问题多数是上述三类非线性问题的组合[1,2]。实际轨道结构受载时的力学行为,属于典型的非线性力学问题。比如基于轮轨接触的材料非线性、几何非线性及边界条件非线性问题,以及扣件、钢轨、垫层三者间相互作用时所表现的边界条件非线性行为等。所以,机车车辆在轨道结构上行驶时引起的力学现象是相当复杂的。以往在研究轨道各部分应力应变分布规律时,通常采用连续弹性基础梁理论或连续点支承,偶尔简单考虑扣件的作用和弹性垫层的使用。不管用哪一种支承方式建立模型,都由于这样那样的假设而带有一定程度的近似性。所以,如何利用现代力学理论的最新成果以及日益发展的计算机技术,根据轨道结构的具体情况,建立更为完整更为准确的轨道结构计算模型,为轨道设计部门提供更加可靠的设计依据或研究思路,已十分必要。本文提出了用非线性有限元理论研究轮轨系统和轨道结构的思路。作为算例之一,本文将根据非线性有限元理论,建立能准确反映扣件、钢轨与垫层的拧紧接触,以及受载车轮与钢轨侧向滑动接触的力学计算模型。 2 轨道结构的有限元接触模型对于非线性问题,不管是材料非线性、几何非线性,还是边界条件非线性,总是最终归结为求解一组非线性平衡方程及其控制方程。例如用位移作为未知数进行有限元分析时,最后可得到一组平衡方程及其控制方程为 : 图1 轮轨系统的对称性模型简图 [K(u)]{u}={R}(1) (u)= (u)(2)其中:{u}为节点位移列阵;{R}为节点载荷列阵; [K(u)]为总体刚度矩阵; (u)为边界条件。它们 36 姜建华:同济大学工程力学系,副教授、博士,上海200092

有限元仿真技术的发展及其应用

有限元仿真技术的发展及其应用许荣昌孙会朝(技术研发中心) 摘要:介绍了目前常用的大型有限元分析软件的现状与发展,对其各自的优势进行了分析,简述了有限元软件在冶金生产过程中的主要应用领域及其发展趋势,对仿真技术在莱钢的应用进行了展望。关键词:有限元仿真冶金生产发展趋势 0 前言自主创新,方法先行,创新方法是自主创新的根本之源,同时,随着市场竞争的日益激烈,冶金企业的产品设计、工艺优化也由经验试错型向精益研发方向发展,而有限元仿真技术正是这种重要的创新方法。近年来随着计算机运行速度的不断提高,有限元分析在工程设计和分析中得到了越来越广泛的应用,比如,有限元分析在冶金、航空航天、汽车、土木建筑、电子电器、国防军工、船舶、铁道、石化、能源、科学研究等各个领域正在发挥着重要的作用,主要表现在以下几个方面:增加产品和工程的可靠性;在产品的设计阶段发现潜在的问题;经过分析计算,采用优化设计方案,降低原材料成本;缩短产品研发时间;模拟试验方案,减少试验次数,从而减少试验成本。与传统设计相比,利用仿真技术,可以变经验设计为科学设计、变实测手段为仿真手段、变规范标准为分析标准、变传统分析技术为现代的计算机仿真分析技术,从而提高产品质量、缩短新产品开发周期、降低产品整体成本、增强产品系统可靠性,也就是增强创新能力、应变能力和竞争力(如图1、2) 。图1 传统创新产品(工艺优化)设计过程为大循环作者简介:许荣昌(1971-),男,1994年毕业于武汉钢铁学院钢铁冶金专业,博士,高级工程师。主要从事钢铁工艺技术研究工作。图2 现代CA E 创新产品(工艺优化)设计过程为小循环 1 主要有限元分析软件简介目前,根据市场需求相继出现了各种类型的应用软件,其中NASTRAN 、ADI N A 、ANSYS 、 ABAQUS 、MARC 、MAGSOFT 、COS MOS 等功能强大的CAE 软件应用广泛,为实际工程中解决复杂的理论计算提供了非常有力的工具。但是,各种软件均有各自的优势,其应用领域也不尽相同。本文将就有限元的应用范围及当今国际国内C AE 软件的发展趋势做具体的阐述,并对与冶金企业生产过程密切相关的主要有限元软件ANSYS 、AB AQUS 、MARC 的应用领域进行分析。 M SC So ft w are 公司创建于1963年,总部设在美国洛杉矶,M SC M arc 是M SC Soft w are 公司于1999年收购的MARC 公司的产品。MARC 公司始创于1967年,是全球首家非线性有限元软件公司。经过三十余年的发展,MARC 软件得到学术界和工业界的大力推崇和广泛应用,建立了它在全球非线性有限元软件行业的领导者地位。随着M arc 软件功能的不断扩展,软件的应用领域也从开发初期的核电行业迅速扩展到航空、航天、汽车、造船、铁道、石油化工、能源、电子元件、机械制造、材料工程、土木建筑、医疗器材、冶金工艺和家用电器等,成为许多知名公司和研究机构研发新产品和新技术的重要工具。在航空业M SC N astran 软件被美国联邦航空管理局(F AA )认证为领取飞行器适 13

有限元分析软件比较分析

有限元分析软件有限元分析是对于结构力学分析迅速发展起来的一种现代计算方法。它是50 年代首先在连续体力学领域--飞机结构静、动态特性分析中应用的一种有效的数值分析方法，随后很快广泛的应用于求解热传导、电磁场、流体力学等连续性问题。有限元分析软件目前最流行的有：ANSYS、ADINA、ABAQUS、MSC 四个比较知名比较大的公司，其中ADINA、ABAQUS 在非线性分析方面有较强的能力目前是业内最认可的两款有限元分析软件，ANSYS、MSC 进入中国比较早所以在国内知名度高应用广泛。目前在多物理场耦合方面几大公司都可以做到结构、流体、热的耦合分析，但是除ADINA 以外其它三个必须与别的软件搭配进行迭代分析，唯一能做到真正流固耦合的软件只有ADINA。ANSYS是商业化比较早的一个软件，目前公司收购了很多其他软件在旗下。ABAQUS专注结构分析目前没有流体模块。MSC是比较老的一款软件目前更新速度比较慢。ADINA是在同一体系下开发有结构、流体、热分析的一款软件，功能强大但进入中国时间比较晚市场还没有完全铺开。结构分析能力排名：ABAQUS、ADINA、MSC、ANSYS 流体分析能力排名：ANSYS、ADINA、MSC、ABAQUS 耦合分析能力排名：ADINA、ANSYS、MSC、ABAQUS 性价比排名：最好的是ADINA，其次ABAQUS、再次ANSYS、最后MSC ABAQUS 软件与ANSYS 软件的对比分析： 1．在世界范围内的知名度：两种软件同为国际知名的有限元分析软件，在世界范围内具有各自广泛的用户群。ANSYS 软件在致力于线性分析的用户中具有很好的声誉，它在计算机资源的利用，用户界面开发等方面也做出了较大的贡献。ABAQUS软件则致力于更复杂和深入的工程问题，其强大的非线性分析功能在设计和研究的高端用户群中得到了广泛的认可。由于ANSYS 产品进入中国市场早于ABAQUS，并且在五年前ANSYS 的界面是当时最好的界面之一，所以在中国，ANSYS 软件在用户数量和市场推广度方面要高于ABAQUS。但随着ABAQUS北京办事处的成立，ABAQUS软件的用户数目和市场占有率正在大幅度和稳步提高，并可望在今后的几年内赶上和超过ANSYS。 2．应用领域：ANSYS 软件注重应用领域的拓展，目前已覆盖流体、电磁场和多物理场耦合等十分广泛的研究领域。ABAQUS 则集中于结构力学和相关领域研究，致力于解决该领域的深层次实际问题。 3．性价比：ANSYS 软件由于价格政策灵活，具有多种销售方案，在解决常规的

(完整word版)有限元分析软件的比较

有限元分析软件的比较（购买必看）－转贴随着现代科学技术的发展，人们正在不断建造更为快速的交通工具、更大规模的建筑物、更大跨度的桥梁、更大功率的发电机组和更为精密的机械设备。这一切都要求工程师在设计阶段就能精确地预测出产品和工程的技术性能，需要对结构的静、动力强度以及温度场、流场、电磁场和渗流等技术参数进行分析计算。例如分析计算高层建筑和大跨度桥梁在地震时所受到的影响，看看是否会发生破坏性事故；分析计算核反应堆的温度场，确定传热和冷却系统是否合理；分析涡轮机叶片内的流体动力学参数，以提高其运转效率。这些都可归结为求解物理问题的控制偏微分方程式，这些问题的解析计算往往是不现实的。近年来在计算机技术和数值分析方法支持下发展起来的有限元分析（FEA，Finite Element A nalysis）方法则为解决这些复杂的工程分析计算问题提供了有效的途径。在工程实践中，有限元分析软件与CAD系统的集成应用使设计水平发生了质的飞跃，主要表现在以下几个方面：增加设计功能，减少设计成本；缩短设计和分析的循环周期；增加产品和工程的可靠性；采用优化设计，降低材料的消耗或成本；在产品制造或工程施工前预先发现潜在的问题；模拟各种试验方案，减少试验时间和经费；进行机械事故分析，查找事故原因。在大力推广CAD技术的今天，从自行车到航天飞机，所有的设计制造都离不开有限元分析计算，FEA在工程设计和分析中将得到越来越广泛的重视。国际上早20世纪在50年代末、60年代初就投入大量的人力和物力开发具有强大功能的有限元分析程序。其中最为著名的是由美国国家宇航局（NASA）在1965年委托美国计算科学公司和贝尔航空系统公司开发的NASTRAN有限元分析系统。该系统发展至今已有几十个版本，是目前世界上规模最大、功能最强的有限元分析系统。从那时到现在，世界各地的研究机构和大学也发展了一批规模较小但使用灵活、价格较低的专用或通用有限元分析软件，主要有德国的ASKA、英国的PA FEC、法国的SYSTUS、美国的ABQUS、ADINA、ANSYS、BERSAFE、BOSOR、COSMOS、ELAS、MARC和STARDYNE等公司的产品。以下对一些常用的软件进行一些比较分析： 1. LSTC公司的LS-DYNA系列软件

非线性有限元分析

非线性有限元分析 1 概述在科学技术领域，对于许多力学问题和物理问题，人们已经得到了它们所应遵循的基本方程（常微分方程或偏微分方程）和相应的定解条件（边界条件）。但能够用解析方法求出精确解的只是少数方程性质比较简单，并且几何形状相当规则的问题。对于大多数工程实际问题，由于方程的某些特征的非线性性质，或由于求解区域的几何形状比较复杂，则不能得到解析的答案。这类问题的解决通常有两种途径。一是引入简化假设，将方程和几何边界简化为能够处理的情况，从而得到问题在简化状态下的解答。但是这种方法只是在有限的情况下是可行的，因为过多的简化可能导致误差很大甚至是错误的解答。因此人们多年来一直在致力于寻找和发展另一种求解途径和方法——数值解法。特别是五十多年来，随着电子计算机的飞速发展和广泛应用，数值分析方法已成为求解科学技术问题的主要工具。已经发展的数值分析方法可以分为两大类。一类以有限差分法为代表，主要特点是直接求解基本方程和相应定解条件的近似解。其具体解法是将求解区域划分为网格，然后在网格的结点上用差分方程来近似微分方程，当采用较多结点时，近似解的精度可以得到改善。但是当用于求解几何形状复杂的问题时，有限差分法的精度将降低，甚至发生困难。另一类数值分析方法是首先建立和原问题基本方程及相应定解条件相等效的积分提法，然后再建立近似解法并求解。如果原问题的方程具有某些特定的性质，则它的等效积分提法可以归结为某个泛函的变分，相应的近似解法实际上就是求解泛函的驻值问题。诸如里兹法，配点法，最小二乘法，伽辽金法，力矩法等都属于这一类方法。但此类方法也只能局限于几何形状规则的问题，原因在于它们都是在整个求解区域上假设近似函数，因此，对于几何形状复杂的问题，不可能建立合乎要求的近似函数。 1960年，R.W.CLOUGH发表了有限单元法的第一篇文献“The Finite Element Method in Plane Stress Analysis”，这同时也标志着有限单元法（FEM）的问世。有限单元法的基本思想是将连续的求解区域离散为一组有限个，且按一定方式相互联接在一起的单元的组合体。由于单元能按不同的联结方式进行组合，且单元本身又可以有不同形状，因此可以模型化几何形状复杂的求解域。并且可以利用在每一个单元假设的近似函数来分片地表示全求解域上待求的未知场函数，从而使一个连续的无限自由度问题变成离散的有限自由度问题。现已证明，有限单元法是基于变分原理的里兹法的另一种形式，从而使里兹法分析的所有理论基础都适用于有限单元法，确认了有限单元法是处理连续介质问题的一种普遍方法。利用变分原理建立有限元方程和经典里兹法的主要区别是有限单元法假设的近似函数不是在全求解域而是在单元上规定的，而且事先不要求满足任何边界条件，因此可以用来处理很复杂的连续介质问题。在短短四十余年的时间里，有限单元的分析方法已经迅速地发展为适合于使用各种类型计算机解决复杂工程问题的一种相当普及的方法。如今，有限元广泛地应用于各个学科门类，已经成为工程师和科研人员用于解决实际工程问题，进行科学研究不可或缺的有力工具。有限单元法的应用围已由弹性力学平面问题扩展到空间问题，板壳问题，由静力平衡问题扩展到稳定问题，动力问题和波动问题。分析的对象从弹性材料扩展到塑性，粘弹性，粘塑性和复合材料等，从固体

有限元分析在结构分析和计算机仿真中的应用

第20卷增刊重庆交通学院学报2001年11月VoI.20Sup.JOURNAL OF CHONGOING JIAOTONG UNIVERSITY NoV.， ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 2001文章编号：1001-716 （2001）S0-0124-03 有限元分析在结构分析和计算机仿真中的应用" 韩西，钟厉，李博（重庆交通学院结构工程部级重点实验室，重庆400074）摘要：简要论述了有限元分析方法在结构分析和计算机仿真的发展趋势和应用情况. 关键词：有限元分析；结构分析；计算机仿真中图分类号：TU311.41文献标识码：B 自1943年数学家Courant第一次尝试应用定义在三角形区域上的分片连续函数的最小位能原理求解St.Venant扭转问题以来，许多数学家、物理学家和工程师由于各种原因都涉及过有限单元的概念.但由于即使一个小规模的工程问题，用有限元分析都将产生较大的计算工作量.直到1960年后，随着计算机技术的发展，有限元分析这门特别依赖数值计算的学科才真正进入了飞速发展阶段.到目前为止，有限元法已成为最强有力的数值分析方法之一，在固体力学、流体力学、机械工程、土木工程、电气工程等领域得到了广泛的应用.由于其所涉及问题和算法基本上全部来源于工程实际、应用于工程中，其解决工程实际问题的能力愈来愈强.由于计算机技术作为有限元分析的计算平台和应用支撑工具，故有限元分析成为CAE（Computer Aided Engineering，计算机辅助工程）这一学科类的主要研究内容.与此同时，由于有限元分析所建模型具有和实际结构相对应的几何、材料、力学特性，对实际结构具有“真实”的模拟特性，和单纯的几何仿真有着本质的区别，所以可以称之为“真实的仿真”（ReaIity SimuIation），可以想象，其模型和计算的数据量将比单纯的几何仿真要大得多，当前，计算机并行多处理器技术正迅猛发展，如SGI OONU-MA体系使计算能力达到工程应用水平，极大地促进了有限元分析计算的发展［2］. 1现状与发展趋势 1.1现状 1956年，Tuner，CIough等人将刚架位移法推广应用于弹性力学平面问题，并用于分析飞机结构，这是现代有限元法第一次成功的尝试.他们第一次给出了用三角形单元求解平面应力问题的正确解答，其研究工作打开了利用计算机求解复杂平面问题的新局面.1963~1964年，BesseIing、MeIosh 和Jones等证明有限单元法是基于变分原理的Ritz 法的另一种形式，从而使Ritz分析的所有理论基础都适用于有限元法，确认了有限元法是处理连续介质问题的一种普遍方法. 几十年来，有限元法的应用已由弹性力学平面问题扩展到空间问题、板壳问题，由静力平衡问题扩展到稳定问题、动力问题和波动问题；分析的对象从弹性材料扩展到塑性、粘塑性和复合材料等；从固体力学扩展到流体力学、传热学等连续介质力学领域.在工程分析中的作用已从分析比较扩展到优化设计并和CAD（计算机辅助设计）结合越来越紧密. 有限元分析理论的逐步成熟主要经历了60年代的探索发展时期，70~80年代的独立发展专家应用时期和90年代与CAD相辅相成的共同发展、推广使用时期. 有限元分析作为一种强有力的数值分析方法，在结构分析和仿真计算中有着极大的应用价值.目前，结构仿真中的静力分析、动力分析、稳定性计 "收稿日期：2001-03-19 作者简介：韩西（1964-），男，重庆人，工学博士，副教授，主要从事振动工程、结构损伤识别、结构动力及计算机仿真分析研究.

第10章(非线性有限元)分解

公式号、图号等第十章非线性动力有限元法当机械结构受到较大的外载荷，或受到持续时间较短的冲击载荷作用时，结构会产生过大的变形, 以至于必须考虑结构几何大变形对结构整体刚度及固有频率的影响，即所谓的几何非线性影响。另外, 对于多数非线性动力学问题，还需要考虑材料非线性、接触非线性等方面的影响。非线性动力学分析求解的基本方程有如下形式 0=-+P I u M (4.141) 式中，Ku u C I += 为粘性效应项，考虑阻尼、粘塑、粘弹等效应。P 为外部激励。对于考虑各种非线性效应的动力学问题求解，需要对动力学方程进行直接时间积分。即非线性动力有限元分析具有如下特点：（1）问题分析过程需要考虑时间积分效应，不必做模态分析，不必提取固有频率；（2）采用直接积分方法求解非线性动力学方程，需要对时间作积分计算，因此计算量远远大于线性模态动力学方法；（3）非线性动力学分析中可以施加不同类型的载荷，包括结点力、非零位移、单元载荷；（4）在每个时间步上，进行质量、阻尼、及刚度的集成，采用完整矩阵，不涉及质量矩阵的近似；（5）可以同时考虑几何、材料和接触等多种非线性效应。非线性动力有限元分析程序常采用隐式Hilber-Hughes-Taylor 法进行时间积分运算。这种方法适于模拟非线性结构的动态问题，对于冲击、地震等激发的结构动态响应以及一些由于塑性或粘性阻尼造成的能量耗散，隐式算法特别有效。隐式积分方法需要对刚度矩阵求逆计算，并通过多次迭代求解增量步平衡方程。隐式Hilber-Hughes-Taylor 时间积分算法为无条件稳定，对时间步长没有特别的限制。采用子空间法也可以对动力学平衡方程作时间积分运算。子空间法是提取模态分析得到的各阶特征模态，并采用与线性模态动力学分析方法相近的分析方式进行求解。对于带有微小非线性效应的问题，如材料小范围进行入屈服、结点转角不大的情况，子空间法效率比进接积分法要高。此外，非线性动力有限元分析还可以采用显式动态算法，如中心差分法。显式时间积分算法为有条件稳定，其临界稳定时间步长限制了时间步长的大小，与有限元模型最小单元尺寸、材料应力波速等有关。显式时间积分法适于模拟高速冲击、接触等问题。上述方法的选择需要综合考虑计算量、分析问题的规模、单元限制等多方面因素，需要丰富的有限元模拟的理论、经验和实践知识。以下以几何非线性问题和材料非线性问题为例介绍非线性有限元法，其中粘弹粘塑性非线性材料问题的分析是典型的非线性动力有限元的求解思想。 9.1 几何非线性问题的有限元法几何非线性问题一般是指物体经历大的刚体位移和转动，但固连于物体坐标系中的应变分量仍假设为小量, 即大位移小应变情况。

空气弹簧力学性态的非线性有限元模拟仿真

空气弹簧力学性态的非线性有限元模拟仿真兰艳1,2, 蔡海涛1, 王成国2, 刘金朝2 江军3 (1.中南大学应用数学系,长沙410083;2.铁道科学院机车车辆研究所,北京100008 3.湘潭大学数学与计算科学学院,湖南湘潭411105)[摘要] 利用非线性有限元方法,在国内第一次给出了提速客车上空气弹簧的全实体单元模拟仿真方法.运用ABAQ US 软件,给出了详细的力学性态分析过程与结果,通过与实验结果的比较,改进了提速客车空气弹簧设计的方法,并提供了重要的理论依据. 关键词:空气弹簧;非线性;有限元;模拟仿真;全实体单元;ABAQUS 中图分类号:O245;O246;TP319 文献标识码:A 文章编号:10005900(2005)01009004 The Nonlinear Finite Element Imitating of Airspring Mechanics Quality LAN Yan 1,2, CAI Hai tao 1, W ANG Chen guo 2, LI U Jin chao 2, JIANG Jun 3 (1.Mathematics Department of Central South Uni versity,Changsha 410083China; 2.Ins ti tute of Equipment,Chinese Academy of Rail way Seieuces,Beijing 100008China; 3.Mathematics Department of Xiangtan Uni versity,Xiangtan 411105China) Abstract We analyze the non-linear character of the air-spring by using the non-linear FE M,and give an si mu - lation method wi th whole substance unit first in China.In compu tation,we get the detail mechanical property analyzing process and result by ABAQUS software.After comparing with the experiment result,author also i mproves the design method of the air-spring. Key words: air-spring;non-linear;FEM;Si mulate;whole substance unit;AB AQUS 空气弹簧的性能对车辆运行的稳定性有重要的影响.新型空气弹簧的研制对改善和提高列车、高速列车和城市轨道车辆的运行品质,延长零部件的使用寿命都具有极其重要的意义,并对其他相关行业的发展也起着重要作用.从早期的主要依靠试验和实际运用中得到的数据来改进研制,到现在欧美、日本及我国的研究部门运用计算机模拟仿真设计的方法来研制开发空气弹簧,在空气弹簧的有限元计算和控制理论方面已经有了一些成果[4,5,6,9] ,不过,国际上的空气弹簧模拟仿真设计因其强非线性性,目前还没有找到解决空气弹簧的非线性性理论问题的办法.本文所做的工作是在文献[3,10,7]的基础上,对空气弹簧的非线性性给以一定说明,并第一次通过有限元软件AB AQUS [1]将空气弹簧的模拟仿真设计,用实体单元给出其相关的力学性态分析.通过与实验结果的比较,改进了提速客车空气弹簧设计的方法,并提供了重要的理论依据.1 空气弹簧模拟仿真中的非线性问题空气弹簧的性能分析过程中,包含有几何非线性、材料非线性、边界非线性三大非线性问题,即为一三重耦合的强非线性问题.空气弹簧的胶囊和橡胶堆的材料是典型的超弹性材料橡胶,其计算就牵涉到了材料的非线性,本文中空气弹簧的材料模型是依Moongy-Rovlin 模型[8]: W =C 1(I 1-3)+C 2(I 2-3) 给出,其中C 1,C 2的取值依赖于实验数据而得.对此模型,既可运用大变形非线性弹性有限元法中的全拉格朗日模式来求解,也可用更新拉格朗日模式来求解.胶囊和橡胶堆在模拟的过程中变形是很大的,传统有限元计算的小变形情况的理论就不再适用了,因而必须考虑有限元计算中的几何非线性问题.而收稿日期:20031115 基金项目:湖南省教育厅基金资助项目(02C571) 作者简介:蔡海涛(1935-)男,湖南南县人,教授,博士生导师,E-mail:lanyan@https://www.360docs.net/doc/e416065562.html,; 第27卷第1期2005年3月湘潭大学自然科学学报Natural Science Journal of Xiangtan University Vol.27No.1Mar.2005

主流CAE有限元分析软件的比较

随着现代科学技术的发展，人们正在不断建造更为快速的交通工具、更大规模的建筑物、更大跨度的桥梁、更大功率的发电机组和更为精密的机械设备。这一切都要求工程师在设计阶段就能精确地预测出产品和工程的技术性能，需要对结构的静、动力强度以及温度场、流场、电磁场和渗流等技术参数进行分析计算。例如分析计算高层建筑和大跨度桥梁在地震时所受到的影响，看看是否会发生破坏性事故；分析计算核反应堆的温度场，确定传热和冷却系统是否合理；分析涡轮机叶片内的流体动力学参数，以提高其运转效率。这些都可归结为求解物理问题的控制偏微分方程式，这些问题的解析计算往往是不现实的。近年来在计算机技术和数值分析方法支持下发展起来的有限元分析（FEA，Finite Element Analysis）方法则为解决这些复杂的工程分析计算问题提供了有效的途径。在工程实践中，有限元分析软件与CAD系统的集成应用使设计水平发生了质的飞跃，主要表现在以下几个方面：增加设计功能，减少设计成本；缩短设计和分析的循环周期；增加产品和工程的可靠性；采用优化设计，降低材料的消耗或成本；在产品制造或工程施工前预先发现潜在的问题；模拟各种试验方案，减少试验时间和经费；进行机械事故分析，查找事故原因。在大力推广CAD技术的今天，从自行车到航天飞机，所有的设计制造都离不开有限元分析计算，FEA 在工程设计和分析中将得到越来越广泛的重视。国际上早20世纪在50年代末、60年代初就投入大量的人力和物力开发具有强大功能的有限元分析程序。其中最为著名的是由美国国家宇航局（NASA）在1965年委托美国计算科学公司和贝尔航空系统公司开发的NASTRAN有限元分析系统。该系统发展至今已有几十个版本，是目前世界上规模最大、功能最强的有限元分析系统。从那时到现在，世界各地的研究机构和大学也发展了一批规模较小但使用灵活、价格较低的专用或通用有限元分析软件，主要有德国的ASKA、英国的PAFEC、法国的SYSTUS、美国的ABQUS、ADINA、ANSYS、BERSAFE、BOSOR、COSMOS、ELAS、MARC和STARDYNE等公司的产品。以下对一些常用的软件进行一些比较分析： 1. LSTC公司的LS-DYNA系列软件 LS-DYNA是一个通用显式非线性动力分析有限元程序，最初是1976年在美国劳伦斯利弗莫尔国家实验室(Lawrence Livermore National Lab.)由J.O.Hallquist 主持开发完成的，主要目的是为核武器的弹头设计提供分析工具，后经多次扩充和改进，计算功能更为强大。此软件受到美国能源部的大力资助以及世界十余家著名数值模拟软件公司（如ANSYS、MSC.software、ETA等）的加盟，极大地加强了其的前后处理能力和通用性，在全世界范围内得到了广泛的使用。在软件的广告中声称可以求解各种三维非线性结构的高速碰撞、爆炸和金属成型等接触非线性、冲击载荷非线性和材料非线性问题。即使是这样一个被人们所称道的数值模拟软件，实际上仍在诸多不足，特别是在爆炸冲击方面，功能相对较弱，其欧拉混合单元中目前最多只能容许三种物质，边界处理很粗糙，在拉格朗日——欧拉结合方面不如DYTRAN灵活。虽然提供了十余种岩土介质模型，但每种模型都有不足，缺少基本材料数据和依据，让用户难于选择和使用。2. MSC.software公司的DYTRAN软件当前另一个可以计算侵彻与爆炸的商业通用软件是MSC.Software Corporation ( MSC公司) 的MSC.DYTR AN程序。该程序在是在LS-DYNA3D的框架下，在程序中增加荷兰PISCES INTERNATIONAL公司开发的PICSES的高级流体动力学和流体——结构相互作用功能，还在PISCES的欧拉模式算法基础上，开发了物质流动算法和流固耦合算法。在同类软件中，其高度非线性、流—固耦合方面有独特之处。MSC.DYTR AN的算法基本上可以概况为：MSC.DYTRAN采用基于Lagrange格式的有限单元方法（FEM）模拟结构的变形和应力，用基于纯Euler格式的有限体积方法（FVM）描述材料（包括气体和液体）流动，对通过流体与固体界面传递相互作用的流体—结构耦合分析，采用基于混合的Lagrange格式和纯Euler 格式的有限单元与有限体积技术，完成全耦合的流体-结构相互作用模拟。MSC.DYTRAN用有限体积法跟踪

通用显式非线性有限元程序：LS-DYNA

通用显式非线性有限元程序：LS-DYNA LS-DYNA 是世界上最著名的通用显式非线性有限元分析程序，能够模拟真实世界的各种复杂问题，特别适合求解各种二维、三维非线性结构的碰撞、金属成型等非线性动力冲击问题，同时可以求解传热、流体及流固耦合问题。在工程应用领域被广泛认可为最佳的分析软件包。与实验的无数次对比证实了其计算的可靠性。 LS-DYNA 是功能齐全的几何非线性（大位移、大转动和大应变）、材料非线性（140多种材料动态模型）和接触非线性（50多种）软件。它以Lagrange 算法为主，兼有ALE 和Euler 算法；以显式求解为主，兼有隐式求解功能；以结构分析为主，兼有热分析、流体-结构耦合功能；以非线性动力分析为主，兼有静力分析功能（如动力分析前的预应力计算和薄板冲压成型后的回弹计算）;是通用的结构分析非线性有限元程序。特色功能 ? 显式求解为主，兼有隐式算法，适合于求解高度非线性问题； ? 具有多种求解算法，以Lagrange 算法为主，兼有ALE、Euler 算法、SPH （Smoothed Particle Hydrodynamics）光顺质点流体动力算法和边界元法BEM（Boundary Element Method）； ? 具有160多种材料模型，是材料模型非常丰富的有限元软件； ? 具有50多种接触类型，是接触类型非常齐全的有限元软件； ? 极好的并行计算能力，包括分布式并行算法（MPP）和共享内存式并行（SMP）； ? 良好的自适应网格剖分技术，包括自适应网格细分和粗化； ? 行业化的专用功能：如针对汽车行业的安全带单元、滑环、预紧器、牵引器、传感器、加速计、气囊等。客户价值 ? 拥有显式和隐式算法，各向异性材料模型，使得板成型、回弹、预应力计算等，可以连续求解； ? 多种控制选项和用户子程序使得用户在定义和分析问题时有很大的灵活性； ? MPP 版本大幅度减少计算时间，计算效率随计算机数目增多而显著提高； ? 与大多数的CAD/CAE 软件集成并有接口。广州有道科技培训中心 h t t p ://w w w .020f e a .c o m

有限元分析软件优势对比

几款常用有限元软件的比较 (2013-03-08 16:12:16) 转载▼ 分类：CAE 标签： ansys 经常有朋友问到常用几款有限元软件的优劣问题，笔者也很关注，也用过一些常用软件，下面首先转发网上一些朋友对于常用有限元软件的评价，最后表明自己的观点。下面是网上一些朋友的评价。 ----------------------------------------------------------------------------------------------- 目前流行的CAE分析软件主要有NASTRAN、ADINA 、ANSYS、ABAQUS、MARC、MAGSOFT、COSMOS等。以下为对这些常用的软件进行的比较和评价： LSTC公司的LS-DYNA系列软件。 LSDYNA长于冲击、接触等非线性动力分析。LS-DYNA是一个通用显式非线性动力分析有限元程序，最初是1976年在美国劳伦斯利弗莫尔国家实验室由J.O.Hallquist主持开发完成的，主要目的是为核武器的弹头设计提供分析工具，后经多次扩充和改进，计算功能更为强大。虽然该软件声称可以求解各种三维非线性结构的高速碰撞、爆炸和金属成型等接触非线性、冲击载荷非线性和材料非线性问题，但实际上它在爆炸冲击方面，功能相对较弱，其欧拉混合单元中目前最多只能容许三种物质，边界处理很粗糙，在拉格朗日——欧拉结合方面不如DYTRAN灵活。 MSC.software公司的DYTRAN软件在同类软件中，DYTRAN在高度非线性、流固耦合方面有独特之处。MSC.DYTRAN程序是在LS-DYNA3D的框架下，在程序中增加荷兰PISCES;INTERNATIONAL公司开发的PICSES的高级流体动力学和流体结构相互作用功能，还在PISCES的欧拉模式算法基础上，开发了物质流动算法和流固耦合算法发展而来的。但是，由于MSC.DYTRAN是一个混合物，在继承了 LS-DYNA3D与PISCES优点的同时，也继承了其不足。首先，材料模型不丰富，对于岩土类处理尤其差，虽然提供了用户材料模型接口，但由于程序本身的缺陷，难于将反映材料特性的模型加上去；其次，没有二维计算功能，轴对称问题也只能按三维问题处理，使计算量大幅度增加；在处理冲击问题的接触算法上远不如当前版的LS-DYNA3D全面。 HKS公司的ABAQUS软件 ABAQUS是一套先进的通用有限元系统，属于高端CAE软件。它长于非线性有限元分析,可以分析复杂的固体力学和结构力学系统，特别是能够驾驭非常庞大的复杂问题和模拟高度非线性问题。ABAQUS不但可以做单一零件的力学和多物理场的分析，同时还可以做系统级的分析和研究，其系统级分析的特点相对于其他分析软件来说是独一无二的。需要指出的是，ABAQUS对爆炸与冲击过程的模拟相对不如DYTRAN和LS-DYNA3D。 ADINA ADINA是近年来发展最快的有限元软件，它独创有许多特殊解法, 如劲度稳定法(Stiffness Stabilization),自动步进法(Automatic Time Stepping),外力-变位同步控制法(Load-Displacement Control)以及BFGS梯度矩阵更新法,使得复杂的非线性问题(如接触,

有限元分析的典型应用领域

有限元分析的典型应用领域 6-4：对该问题进行有限元分析的过程如下。 (1)进入ANSYS（设定工作目录和工作文件）程序→ANSYS→ANSYS Interactive →Working directory（设置工作目录）→Initial jobname（设置工作文件名）：Press →Run →OK (2)设置分析特性 ANSYSMain Menu：Preferences…→Structural →OK (3)定义单元类型 ANSYSMain Menu：Preprocessor →Element Type →Add/Edit/Delete... →Add…→Solid: Quad 4node 42 →OK（返回到Element Types窗口）→Options…→K3：Plane Strs w/thk（带厚度的平面应力问题）→OK →Close (4)定义材料参数 ANSYSMain Menu：Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic：EX：2.1e11（弹性模量），PRXY：0.3（泊松比）→OK →点击该窗口右上角的“×”来关闭该窗口 (5)定义实常数以确定平面问题的厚度 ANSYSMain Menu：Preprocessor →Real Constants…→Add/Edit/Delete →Add →Type 1 PLANE42 →OK →Real Constant Set No：1（第1号实常数），THK：3.4（平面问题的厚度）→OK →Close （6）生成几何模型生成上拱形梁 ANSYSMain Menu：Preprocessor →Modeling →Create →Keypoints→In Active CS →NPT Keypoint number：1，X，Y，Z Location in active CS：-4.5，8.5→Apply →同样输入后5个特征点坐标（坐标分别为（-2.25，8.5），（2.25，8.5），（4.5，8.5），（0，13），（0，10.75））→OK →Lines →Lines →Straight Line 用鼠标分别连接特征点1，2和3，4生成直线→OK→Arcs →By End KPs & Rad →用鼠标点击特征点2，3 →OK →用鼠标点击特征点6 →OK →RAD Radius of the arc：2.25→Apply （出现Warning对话框，点Close关闭）→用鼠标点击特征点1，4 →OK →用鼠标点击特征点5 →OK →RAD Radius of the arc：4.5→OK（出现Warning对话框，点Close关闭）→Areas →Arbitrary →By Lines →用鼠标点击刚生成的线→OK 生成下拱形梁 ANSYSMain Menu：Preprocessor →Modeling →Create →Keypoints→In Active CS →NPT Keypoint number：7，X，Y，Z Location in active CS：-4.5，-8.5→Apply →同样输入后5个特征点坐标（坐标分别为（-2.25，-8.5），（2.25，-8.5），（4.5，-8.5），（0，-13），（0，-10.75）→OK →Lines→Lines →Straight Line →用鼠标分别连接特征点7，8和9，10生成直线→OK →Arcs →By End KPs & Rad →用鼠标点击特征点8，9 →OK用鼠标点击特征点12 →OK →RAD Radius of the arc：2.25→Apply （出现Warning对话框，点Close关闭）→用鼠标点击特征点7，10 →OK →用鼠标点击特征点11 →OK →RAD Radius of the arc：4.5→OK（出现Warning对话框，点Close关闭）→Areas →Arbitrary →By Lines →用鼠标点击刚生成的线→OK 生成两根立柱 ANSYSMain Menu：Preprocessor →Modeling →Create →Areas →Rectangle →By 2 Corners →WP X：-4.5，WP Y：-8.5，Width：2.25，Height：17→Apply →WP X：2.25，WP Y：-8.5，Width：2.25，Height：17→OK 粘结所有面

(完整版)国内外主要有限元分析软件比较

有限元分析是对于结构力学分析迅速发展起来的一种现代计算方法。它是50年代首先在连续体力学领域--飞机结构静、动态特性分析中应用的一种有效的数值分析方法，随后很快广泛的应用于求解热传导、电磁场、流体力学等连续性问题。有限元分析软件目前最流行的有：ANSYS、ADINA、ABAQUS、MSC四个比较知名比较大的公司。常见软件有限元分析软件目前最流行的有：ANSYS、ADINA、ABAQUS、MSC四个比较知名比较大的公司，其中ADINA、ABAQUS在非线性分析方面有较强的能力目前是业内最认可的两款有限元分析软件，ANSYS、MSC进入中国比较早所以在国内知名度高应用广泛。目前在多物理场耦合方面几大公司都可以做到结构、流体、热的耦合分析，但是除ADINA以外其它三个必须与别的软件搭配进行迭代分析，唯一能做到真正流固耦合的软件只有ADINA。软件对比 ANSYS是商业化比较早的一个软件，目前公司收购了很多其他软件在旗下。ABAQUS专注结构分析目前没有流体模块。MSC是比较老的一款软件目前更新速度比较慢。ADINA是在同一体系下开发有结构、流体、热分析的一款软件，功能强大但进入中国时间比较晚市场还没有完全铺开。结构分析能力排名：1、ABAQUS、ADINA、MSC、ANSYS 流体分析能力排名：1、ANSYS、ADINA、MSC、ABAQUS 耦合分析能力排名：1、ADINA、ANSYS、MSC、ABAQUS 性价比排名：最好的是ADINA，其次ABAQUS、再次ANSYS、最后MSC ABAQUS软件与ANSYS软件的对比分析 1．在世界范围内的知名度两种软件同为国际知名的有限元分析软件，在世界范围内具有各自广泛的用户群。ANSYS 软件在致力于线性分析的用户中具有很好的声誉，它在计算机资源的利用，用户界面开发等方面也做出了较大的贡献。ABAQUS软件则致力于更复杂和深入的工程问题，其强大的非线性分析功能在设计和研究的高端用户群中得到了广泛的认可。由于ANSYS产品进入中国市场早于ABAQUS，并且在五年前ANSYS的界面是当时最好的界面之一，所以在中国，ANSYS软件在用户数量和市场推广度方面要高于ABAQUS。但随着ABAQUS北京办事处的成立，ABAQUS软件的用户数目和市场占有率正在大幅度和稳步提高，并可望在今后的几年内赶上和超过ANSYS。 2．应用领域