非线性有限元

非线性有限元介绍1.为什么使用FEA解决有限元问题（1）理解设计的意图。

有限元分析（FEA）是研究不同力学设计的有力工具。

（2）降低产品成本和开发周期。

1) FEA通过以下方式降低产品成本和开发周期；2) 在模具制造前识别成型问题；3) 使模具制造返工成本最低；4) 提前识别设计中缺陷减小样机的成本；5) 使用最少的材料；（3）获得结果的唯一办法。

1) FEA可以用来预测产品在极端工况下的性能，这些在实验中无法复现；2) 能在设计阶段提前考虑这些工况。

（4）很多工况在设计阶段无法预料。

2.收敛定义（1）在有限元中收敛有多重意义：1) 网格收敛；2) 时间积分精度；3) 非线性程序收敛；4) 求解精度；（2）网格收敛1) 增加模型单元数量会使仿真解趋于解析解。

网格收敛对线性和非线性问题都适用；在Abaqus中使用H网格自适应技术；2)进一步加密网格时，结果变化很小或不变时，认为网格达到收敛。

3)网格收敛规则的例外：网格奇异解；材料损伤累计在模型特定区域的局部问题；在Abaqus中使用H自适应网格技术；Abaqus提供特殊技术来减小网格依赖性，解决材料软化局部影响。

4)Abaqus提供评估网格收敛工具。

在打印输出文件（.dat）和结果文件（.fil）输出的节点（SJP）处应变跳变。

5)Abaqus后处理云图设置。

应力云图；不连续云图。

6)自适应网格误差设计。

见Abaqus/Standard 自适应网格课程介绍。

（3）瞬态问题时间积分精度。

对于具有物理时间尺度的瞬态问题，Abaqus提供用户定义参数，以控制对相关方程的积分精度：半增量容差。

1)评估当前增量中途点的最大不平衡力；2)一个增量中允许的最高温度变化；3)增量中根据开始和结束速率条件计算出的蠕变应变增量的最大容差。

（4）非线性求解收敛：见后续。

（5）求解精度。

获得精确解需要满足以下条件：准确解需要工程经验来建立合适有限元模型：材料、载荷、边界和求解程序。

非线性有限元法综述摘要：本文针对非线性有限元法进行综述，分别从UL列式及TL列式、CR列式、几何精确梁、壳理论三个方面介绍其分析思路和发展动态，旨在为相关学者提供一些思路参考。

关键词：几何非线性；UL列式；TL列式；CR列式；几何精确梁、壳理论1引言几何非线性是由于位置改变引起了结构非线性响应。

进行结构几何非线性分析，实质上就是要得到结构真实的变形与受力情况。

有限元方法是进行结构几何非线性分析的最成熟的方法，也是应用最广泛的分析方法.2非线性有限元法研究思路非线性有限元法主要指UL列式法、TL列式法、CR列式法和几何精确梁、壳理论等，它们有着基本相同的思路，即利用虚功原理建立平衡方程。

方程中充分考虑了非线性因素对结构应变和应力的影响，也就是将线性应变和非线性应变都代入到表达式中，然后确定单元的本构关系并选取合适的形函数，导出单元对应的弹性刚度矩阵和几何刚度矩阵，再选取合适的增量-迭代算法进行求解，由此就完成了结构的整个几何非线性分析求解过程。

非线性有限元法将结构的变形过程划分为三个主要阶段：C0状态、C1状态和C2状态，如图1所示。

图1 单元的变形C0状态是单元的初始状态，C1状态是单元受力变形后上一次处于平衡的状态；C2状态是单元的当前状态，也就是所求的状态。

2.1UL法和TL法研究思路UL法和TL法为几何非线性问题提供了新的分析思路。

这两种方法本质上没有很大区别，但是方程建立的参考状态有所不同。

完全拉格朗日法(TL法)是以结构变形前C0状态为参考建立平衡方程的，考虑结构从C0状态到C2状态之间的变形；而更新的拉格朗日法(UL法)以结构变形后C1状态为参考建立平衡方程的[2]，考虑结构从C1状态到C2状态之间的变形。

两种拉格朗日法的主要形式如下：（1）TL列式（2）UL列式从上面两式可以看出：TL法和UL法的另一个不同是TL法的增量平衡方程中考虑了初位移矩阵的影响，而UL法则忽略了其影响，只考虑了弹性刚度矩阵和初应力矩阵的影响。

非线性有限元方法非线性有限元方法是大量应用于工程领域的计算方法，它主要用于求解复杂结构的力学问题，例如材料的变形、破坏和变形控制等。

与线性有限元方法不同，非线性有限元方法考虑因为载荷和边界条件的非线性导致问题的非线性本质，以及材料的非线性行为。

在这篇文章中，我们将讨论非线性有限元方法，包括其应用、工作原理以及其在工程领域中的重要性等内容。

首先，我们来研究一下非线性有限元方法的应用。

非线性有限元方法在许多方面都有应用。

其中最重要的领域是结构力学，包括建筑、航空航天、汽车等领域。

由于这些结构需要承受复杂的载荷，因此非线性有限元方法可以很好地模拟这些结构的行为，预测它们的性能和寿命。

此外，非线性有限元方法还可以应用于材料力学研究中，例如破碎、断裂和塑性变形等方面。

其次，我们来了解一下非线性有限元方法的工作原理。

与线性有限元方法类似，非线性有限元方法通过将结构分成小块进行离散，然后在每个小块中进行力学分析，最后将分析结果合并为整个结构的行为。

但是，与线性有限元方法不同的是，非线性有限元方法考虑到材料的非线性行为，采用迭代的方法计算结构的响应。

通常，在每一次迭代中，我们都将结构的当前状态作为一个初始猜测，然后求解出该状态下的切应力和位移场。

然后我们将这个位移场的结果代入底部，从而更新结构的状态。

如果解决方案收敛，则完成计算，否则就将新的状态再次代入求解。

这种方法的本质是将非线性问题转化为一系列线性问题的求解，通过迭代求解来逼近非线性问题的解。

最后，我们来讨论一下非线性有限元方法在工程领域中的重要性。

非线性有限元方法已成为现代工程设计和分析的不可或缺的工具。

它允许工程师们模拟和预测各种工程机构的行为，以及设计和优化各种结构。

例如，它可以帮助我们了解在不同载荷下建筑和桥梁行为的变化，预测材料的破坏和失效，以及优化汽车和飞机的结构以提高其性能。

总之，非线性有限元方法是一种复杂但十分有用的计算方法，它可以模拟各种结构的行为并预测其性能和寿命。

如何利用非线性有限元法进行力学分析非线性有限元法是一种用于数值分析问题的计算方法，其主要应用于力学分析领域。

这种方法在于其对于复杂结构的建模能力和高精度数值计算能力而备受推崇。

在本文中，将介绍如何对力学问题进行分析，以及如何应用非线性有限元法对力学分析进行模拟。

1. 引言力学分析整体上分为两种类型：静力学分析和动力学分析。

静力学分析研究对于物体的力和静止条件进行研究，其中力一般会造成物体的运动。

而动力学分析则研究运动物体的变化，特别是再一定条件下物体的振动问题等。

因为力学分析问题具有很高的复杂性，很多时候需要使用非线性有限元法来得到更准确的结果。

下面我们将详细介绍使用非线性有限元法进行力学分析的方法和流程。

2. 有限元法简介有限元法是一种现代数值计算方法，它将大工程结构分割为小的有限元。

在每个有限元内，结构的物理性质可以被认为是常量。

(具体内容可以自己百度)3. 如何利用非线性有限元法进行力学分析使用非线性有限元法进行力学分析的核心是将宏观问题转变为微观问题来进行模拟计算。

其中需要注意下面几点：3.1 确定力学分析的类型根据要进行分析的结构本身的性质和应用场景，可能涉及到静力学分析或者动力学分析。

其中静力学分析的计算主要涉及到结构在平衡状态下的情况，而动力学分析主要涉及到结构在某种条件下的运动和振动情况。

因此，在进行力学分析之前需要确定其类型，以便进行后续的计算。

3.2 建立结构模型根据具体情况，需要对结构进行建模。

建模可以通过一定的工具软件实现，或者手工建立结构模型。

模型的建立需要考虑到其复杂性和具体的应用场景。

构建好结构模型之后，需要对其进行精细化剖分得到单元网格，并进行编号。

3.3 确定边界条件在进行力学分析时，还需要考虑结构的边界条件。

边界条件可以通过指定某些点的坐标或者某些角度的变化来确定。

因此，在进行计算时需要根据具体情况设定边界条件，以便进行后续的计算。

3.4 进行数值模拟计算运用有限元法的基本原理，将每个单元的机械性质进行计算，根据力学分析的情况，可以得到结构节点的位移、应变和应力等参数。

非线性有限元之非线性求解方法平衡回顾✧静态平衡是内力I和外载P力量平衡；✧在非线性问题中，模型的内力I可以是以下量的非线性函数；✧在非线性问题中，模型的外力P也可以是某些量的非线性函数，如位移u和时间t。

非线性求解方法1.已知一个分析，知道结构总载荷和初始刚度，目的是找到最后的位移。

线性分析中，一次计算就能求解出最终位移；非线性问题中不可能，因为结构刚度随着结构变形而改变。

2.求解这类非线性问题需要的是一种增量\迭代技术，获得的解是非线性问题准确的近似。

这些方程通常没有精确解。

3.Abaqus使用迭代求解该方程：使用牛顿拉普森方法求解近似解，使误差最小。

4.Abaqus用法：1）载荷历史被拆解为一系列的分析步；每个分析步拆解为一系列增量步；用户为初始时间增量猜测一个值；Abaqus使用自动增量算法确定其他的增量步。

在每个增量步结束时，Abaqus根据载荷与时间关系计算当前负载大小2）使用牛顿拉普森程序迭代求解每个增量结束时的解；根据收敛容差判断牛顿拉普森程序的收敛；如果迭代不收敛，减少增量步的大小；然后使用小增量步重新进行计算。

5.分析步、增量步、迭代步1）分析步仿真载荷历程含有一个或多个分析步。

2）增量步是分析步的一部分；在静态问题中，总载荷被分成很小的增量步。

以便可以沿着非线性路径求解。

3）迭代步迭代步是增量步中寻找平衡解得一次计算尝试。

5.牛顿拉普森方法Abaqus/Standard 基于牛顿拉普森方法的增量迭代求解技术，该方法是无条件稳定（任何大小的增量步都可以）。

增量步大小影响动态分析精度，每个增量步通常要求多次迭代才能满足收敛要求，每个分析步通常有多个增量步，牛顿拉普森定义了一个残差为0位移曲线。

6.牛顿拉普森方法基础。

平衡是u的非线性方程，牛顿拉普森迭代求解在Cu 处的线性方程，Cu是位移u的修正量。

7.残差定义为了得到线性方程组，重写一下平衡方程，R（u）是u的残差。

这个残差表示的是位移u处不平衡力。

第三部分非线性分析第一章非线性有限元概述1.1非线性行为1、 非线性结构的基本特征是结构刚度随载荷的改变而变化。

如果绘制一个非线 性结构的载荷一位移曲线，则 力与位移的关系是非线性函数。

2、 引起结构非线性的原因：a 几何非线性：大应变，大位移，大旋转 （例如钓鱼竿的变形）b 材料非线性：塑性，超弹性，粘弹性，蠕变c 状态改变非线性：接触，单元死活3、 非线性行为一一分析方法特点A 不能使用叠加原理！B 结构响应与路径有关，也就是说加载的顺序可能是重要的。

C 结构响应与施加的载荷可能不成比例。

1.2非线性分析的应用1、 一些典型的非线性分析的应用包括： 非线性屈曲失稳分析金属成形研究碰撞与冲击分析制造过程分析（装配、部件接触等）材料非线性分析 （塑性材料、聚合物）2、 橡胶底密封：一个包含几何非线性（大应变与大变形），材料非线性（橡胶）， 及状态非线性（接触）的例子。

2.1非线性方程组的解法1、求解一个结构的平衡问题通常等于求解结构的总位能的驻值 问题。

结构总位能n ： 口 "3弋门心 2、 增量法：就是将荷载分成一系列的荷载增量，即 ANSYS 中的荷载步或荷载子 步。

A 要点：在每一个荷载增量求解完成后，继续进行下一个荷载增量之前， 刚度矩阵以反映结构刚度的变化。

B 增量法的优点：可以追踪结构变形历程，这对于材料或几何非线性（特别是 极限值屈曲分析）十分有用。

C 增量法的缺点：随着荷载步增量的增加而产生积累误差，导致荷载-位移曲 线飘移。

D 对飘移进行平衡修正，可以大大提高增量法的精度。

应用最广的就是在每一 级载荷增量上用Newton-Raphsor 或其变形的迭代法。

3、 迭代法：割线刚度法：收敛性差，因此很少应用切线刚度法Newto n-Ra phsor 迭代法：切向刚度法中 2.2 Newto n-Ra phsor 迭代法 1、 优点：对于一致的切向刚度矩阵有 二次收敛速度。

一、材料非线性问题的有限单元法1．1 引言以前各章所讨论的均是线性问题。

线弹性力学基本方程的特点是1．几何方程的应变和位移的关系是线性的。

2．物性方程的应力和应变的关系是线性的。

3．建立于变形前状态的平衡方程也是线性的。

但是在很多重要的实际问题中，上述线性关系不能保持。

例如在结构的形状有不连续变化（如缺口、裂纹等）的部位存在应力集中，当外载荷到达一定数值时该部位首先进入塑性，这时在该部位线弹性的应力应变关系不再适用，虽然结构的其他大部分区域仍保持弹性。

又如长期处于高温条件下工作的结构，将发生蠕变变形，即在载荷或应力保持不变的情况下，变形或应变仍随着时间的进展而继续增长，这也不是线弹性的物性方程所能描述的。

上述现象都属于材料非线性范畴内所要研究的问题。

工程实际中还存在另一类所谓几何非线性问题。

例如板壳的大挠度问题，材料锻压成型过程的大应变问题等，这时需要采用非线性的应变和位移关系，平衡方程也必须建立于变形后的状态以考虑变形对平衡的影响。

由于非线性问题的复杂性，利用解析方法能够得到的解答是很有限的。

随着有限单元法在线性分析中的成功应用，它在非线性分析中的应用也取得了很大的进展，已经获得了很多不同类型实际问题的求解方案。

材料非线性问题的处理相对比较简单，不需要重新列出整个问题的表达格式，只要将材料本构关系线性化，就可将线性问题的表达格式推广用于非线性分析。

一般说，通过试探和迭代的过程求解一系列线性问题，如果在最后阶段，材料的状态参数被调整得满足材料的非线性本构关系，则最终得到问题的解答。

几何非线性问题比较复杂，它涉及非线性的几何关系和依赖于变形的平衡方程等问题，因此，表达格式和线性问题相比，有很大的改变，这将在下一章专门讨论。

这两类非线性问题的有限元格式都涉及求解非线性代数方程组，所以在本章开始对非线性代数方程组的求解作—一般性的讨论。

这对下一章也是必要的准备。

正如在前面已指出的，材料非线性问题可以分为两类。

材料非线性有限元分析材料非线性有限元分析是一种重要的计算力学方法，用于研究在载荷作用下，材料会发生非线性行为的情况。

这种分析方法已经被广泛应用于工程领域，例如建筑结构、航空航天以及汽车工业等。

本文将详细介绍材料非线性有限元分析的原理、方法和应用。

首先，我们来介绍一下材料非线性。

在工程领域，材料的非线性行为主要包括弹塑性、损伤、断裂、破坏等。

这些非线性行为往往在高载荷作用下会显著增加结构的应力和应变，从而导致结构的失效。

因此，准确地预测和分析这些非线性行为对于工程设计和结构优化具有重要意义。

材料非线性有限元分析是一种基于有限元方法的计算机模拟技术，用于模拟和分析复杂结构在非线性载荷下的力学行为。

它通过将结构离散为许多小的有限元单元，并以数学模型描述每个单元的材料行为，从而建立了结构的有限元模型。

然后，结构的力学行为可以通过求解相应的离散形式的力学方程得到。

在材料非线性有限元分析中，有两个关键问题需要解决。

首先是材料本构模型的建立。

材料本构模型是描述材料应力和应变关系的数学模型，常用的包括弹性模型、塑性模型、损伤模型等。

选择合适的材料本构模型对准确预测和分析结构的非线性行为至关重要。

其次是数值方法的选择。

对于材料非线性问题，通常需要使用迭代算法，如牛顿-拉夫森法，来求解非线性方程。

此外，还需要选择适当的数值积分方法，以解决离散形式的力学方程。

材料非线性有限元分析在许多领域都有广泛的应用。

在结构工程领域，它可以用于分析钢筋混凝土结构、大跨度桥梁以及高层建筑等的受力性能。

在航空航天领域，材料非线性有限元分析可用于研究飞机机翼、航天器的结构强度和振动特性。

在汽车工业中，它可以用于分析车辆的碰撞、耐久性和振动特性。

总结起来，材料非线性有限元分析是一种重要的计算力学方法，能够准确地模拟和分析结构在非线性载荷下的力学行为。

它在工程领域有着广泛的应用，能够为工程设计和结构优化提供科学依据。

未来随着计算机硬件和数值方法的不断发展，材料非线性有限元分析将在更多领域得到应用，并为解决工程实际问题提供更准确和高效的方法。