《Abaqus非线性教程》由会员分享，可在线阅读，更多相关《Abaqus非线性教程（21页珍藏版）》请在金锄头文库上搜索。
1、第七章非线性这一章讨论ABAQUS中的非线性结构分析。线性分析与非线性分析的区别概述如下：线性分析到目前为止所讨论的分析实例均为线性情形，也就是施加的载荷和系统响应间存在线性关系。例如：如果一线性弹赞在10N的载荷下伸长lm,那么施加20N的载荷就会伸长2m。这意味着在线性分析中，结构的柔度阵只需计算一次（将刚度阵集成并求逆即可得到）。其它载荷情形下，结构的线性响应可通过将新的载荷向量与刚度阵的逆相乘得到。此外，结构对不同载荷情形的响应，可以用常数來进行比例变换或相互叠加的方式，來得到结构对一种完全新的载荷的响应，这要求那种新载荷是先前各载荷的线性组合。载荷的叠加原则假定所有的载荷的边界条件相
2、同。ABAQUS在线性动力学模拟中使用了载荷的叠加原理，这将在第九章进行讨论。非线性分析结构的非线性问题指结构的刚度随其变形而改变的分析问题。实际上所有的物理结构均为非线性的，而线性分析只是一种方便的近似，这对设计來说通常是足够精确的。显然，线性分析对包括加工过程的许多结构模拟來说是不够的，象锻造或冲压、压溃分析、及橡胶部件、轮胎和发动机垫圈分析等问题。一个简单的例子就是具有非线性刚度响应的弹簧（见图7-1）o线性弹簧.刖度杲常数图7-1线性和非线性弹簧特性由丁刚度依赖丁位移，所以不能再用初始柔度乘以所施加的载荷的方法來计算任意载荷时弹簧的位移了。在非线性分析中结构的刚度阵在分析过程中必须进行
3、许多次的生成、求逆，这使得非线性分析求解比线性分析昂贵得多。由丁非线性系统的响应不是所施加载荷的线性函数，因此不可能通过叠加来获得不同载荷的解。每种载荷都必须作为独立的分析进行定义及求解。7#7.1非线性的来源在结构力学模拟中有三种非线性的來源:材料非线性边界非线性儿何非线性7.1.1材料非线性这种非线性也许是人们最熟悉的，并将在第八章中进行更深入的讨论。大多数金属在小应变时都具有良好的线性应力/应变关系，但在应变较大时材料会发生屈服，此时材料的响应变成了非线性和不可逆的（见图7-2）o图7-2弹一塑性材料轴向拉伸应力一应变曲线橡胶可以近似认为具有非线性的、可逆的（弹性）响应的材料（见图7-3
4、）o材料的非线性也可能与应变以外的其它因素有关。应变率相关材料的材料参数和材料失效都是材料非线性的表现形式。材料性质也可以是温度和其它预先设定的场变量的函数。7.1.2边界非线性若边界条件随分析过程发生变化，就会产生边界非线性问题。考虑图74所示7-2的悬臂梁，它随施加的载荷发生挠山I,氏至碰到障碍。图7-4将碰到障碍物的悬臂梁梁端部的竖向挠度与载荷在它接触到障碍以前是线性关系。在端部碰到障碍时梁端部的边界条件发生突然的变化，阻止竖向挠度继续增大，因此梁的响应将不再是线性的。边界非线性是极度不连续的，在模拟分析中发生接触时，结构的响应特性会在瞬时发生很大的变化。另一个边界非线性的例子是将板材材
5、料冲压入模具的过程。在与模具接触前,板材在压力下的伸展变形是相对容易产生的，在与模具接触后，由丁边界条件的改变，必须增加压才能使板材继续成型。边界条件非线性将在第十一章进行讨论。7.1.3几何非线性第三种非线性的來源是与分析过程中模型的儿何改变相联系的。儿何非线性发生在位移的大小影响到结构响应的情形。这可能由于：大挠度或转动。“突然翻转”。初应力或载荷硬化。例如，考虑端部受竖向载荷的悬臂梁（见图7-5）o若端部挠度较小，分析时可以认为是近似线性的。然而若端部的挠度较大，结构的形状乃至丁其刚度都会发生改变。另外，若载荷不能保持与梁垂直，载荷对结构的作用将发生明显的改变。当悬臂梁由端部挠曲时，载荷
6、可以分解为一个垂苴丁梁的分量和另一个沿梁的长度方向作用的分量。所有这些效应都会对悬臂梁的非线性响应作出贡献（也就是梁的刚度随它所承受载荷的增加而不断变化）。图7-5悬臂梁的人挠度可以预料大挠度和转动对结构承载方式有重耍影响。然而，并非位移相对丁结构尺寸很大时，儿何非线性才显得重要。考虑一块很大的弯板在所受压力下的“突然翻转”现象，如图76所示。在此例子中板的刚度在变形时会产生戏剧性的变化。当平板突然翻转时，刚7#度就变成了负的。这样，尽管位移的量值相对丁板的尺寸來说很小，在模拟分析中仍有严重的儿何非线性效应，这是必须加以考虑的。变形前形狀7#7#图7-6人平板的突然翻转7.2非线性问题的求解结
7、构的非线性载荷一位移曲线见图77。分析的目标是确定其响应。ABAQUS使用NewtomRaphson法來求解非线性问题。在非线性分析中的求解不能象线性问题中那样，只求解一组方程即可，而是逐步施加给定的载荷，以增量形式趋于最终解。因此AEAQUS将计算过程分为许多载荷增量步，并在每个载荷增量步结束时寻求近似的平衡构形。ABAQUS通常耍经过若干次迭代才能找到某一载荷增量步的可接受的解。所有增量响应的和就是非线性分析的近似解。考虑作用在物体上的外部作用力尸和内部(节点)作用力Z,(分别见图7-8(a)与图7-8(b)。作用于一节点上的内部作用力是由包含此节点的各单元中的应力引起的。(a)模拟计算中
8、的外部载荷(b)作用于节点上的内部作用力图78物体上的外部载荷和内部作用力为了使物体处于平衡状态，每个节点上施加的净作用力必须为零。因此平衡的基本判据为内部作用力/和外部作用力P必须互相平衡：7.2.1分析步，增量步和迭代步这一节将引入一些新词汇來描述分析中的不同组成部分。清楚地理解分析步,载荷增量步和迭代步的区别是很重耍的。模拟计算的加载过程包含单个或多个步骤，所以要定义分析步。它一般包含分析过程选择，载荷选择，和输出耍求选择。而且每个分析步都可以采用不同的载荷、边界条件、分析过程和输出耍求。例如：步骤一：将板材夹于刚性夹具上。步骤二：加载使板材变形。步骤三：确定变形板材的H然频率。增量步是
9、分析步的一部分。在非线性分析中，一个分析步中施加的总载荷被分解为许多小的增量，这样就可以按照非线性求解步骤來进行计算。当提出初始增量的大小后，ABAQUS会动选择后继的增量大小。每个增量步结束时，结构处于(近似)平衡状态，结果可以写入输出数据库文件、重启动文件、数据文件或结果文件中。选择某一增量步的计算结果写入输出数据库文件的数据称为帧(frames)o迭代步是在一增量步中找到平衡解的一种尝试。如果模型在迭代结束时不是处丁平衡状态，ABAQUS将进行另一轮迭代。随着每一次迭代，ABAQUS得到的解将更接近平衡状态；有时ABAQUS需耍进行许多次迭代才能得到一平衡解。当平衡解得到以后一个增量步才
10、完成，即结果只能在一个增量步的末尾才能获得。7.2.2平衡迭代和收敛性结构对丁一个小的载荷增量AP的非线性响应示于图7-9中。ABAQUS利用基丁时构形的结构初始刚度Ko,和增量AP來计算结构的位移修正值c。利用c.将结构的构形更新为U.o7#Load图7-9增量步中的首次迭代收敛性基丁结构新的构形u，ABAQUS形成新的刚度心利用K來计算更新后的构形中结构的内部作用力/,o所施加的总载荷P和/的差值可如下计算：&=P-L其中R是迭代的作用力残差值。如果&在模型的每一F1由度上均为零，图7-9中的a点将位于载荷一挠度曲线上，结构将处丁平衡状态。在非线性问题中，儿乎不可能使R等于零，因此ABAQ
11、US将R与容许残差进行比较。如果R比作用力容许残差小，ABAQUS就接受结构的更新构形作为平衡结果。默认的容许残差设置为结构中对时间进行平均的作用力的0.5%oABAQUS在整个模拟过程中鬥动从空间分布和对时间平均的角度计算这个值。若R比目前的容许残差小，就认为P和你处于平衡状态，U就是结构在当前载荷下合理的平衡构形。而ABAQUS在接受此解前，还耍检查位移修正值c.与总的增量位移AUa=Ua-Uo相比是否是一小量。若c大丁增量位移的1%,ABAQUS将重新进行迭代。只有这两个收敛性检査都得到满足，才认为此载荷增量下的解是收敛的。上述收敛判断规则有一个例外，即所谓线性增量情况。线性增量的定义是
12、指增量步内最大的力残差小丁时间平均力乘以10的增量步，凡严格满足这个定义的增量步无需再进行迭代，无需进行任何检査即可认为其解是可接受的。若迭代结果不收敛，ABAQUS将进行另一种迭代以使内部和外部作用力达到平衡。第二种迭代采用前而迭代结束时计算得到的刚度K和R.起来确定另一位移修正值5,这使得系统更加接近平衡状态（见图7-10中的点b）。ABAQUS利用结构新构形山中的内部作用力计算新的作用力残值氐，再次将任意IH由度上的最大作用力残值与作用力残值容许值进行比较，将第二种迭代的位7#移修正值Cb与增量位移从二从-山进行比较。如果需耍的话ABAQUS将进行进一步的迭代。对丁非线性分析中的每次迭代
13、，ABAQUS耍重新形成模型的刚度矩阵并求解方程组。从计算费用的角度來说，这意味着每次迭代等价于进行一次完整的线性分析。现在可以清楚地看到非线性分析的计算费用可能要比线性问题大许多倍。可以在每一收敛的增量步上保存结果，所以非线性模拟计算中得到的输出数据量将是线性分析中可得到数据量的很多倍。因此在规划计算机资源时，就应考虑这些因素及所想进行的非线性模拟计算的类型。7.2.3自动增量控制ABAQUS鬥动调整载荷增量步的大小，因此它能便捷而有效地求解非线性问题。用户只需在每个分析步计算中给出第一个增量的大小，ABAQUS会H动调整后续增量的大小。若用户未提供初始增量大小，ABAQUS会试图将该分析步
14、的全部载荷都作为第一增量步载荷來施加，这样在高度非线性的问题中ABAQUS不得不反复减小增量大小，从而导致CPU时间的浪费。一般來说，提供一个合理的初始增量大小（见7.4.1节例题）将是有利的；只有在很平缓的非线性问题中才可能将一分析步中的所有载荷施加于一个增量步中。在一个载荷增量里得到收敛解所需的迭代步数会随系统的非线性程度而变化。默认情况下，如果在16次迭代中仍不收敛或出现发散，AEAQUS会放弃当前增量步，并将增量大小置为先前值的25%,重新开始计算，即利用比较小的载荷增量来尝试找到收敛的解。若此增量仍不收敛，ABAQUS将再次减小增量大小。ABAQUS允许一增量步中最多有五次增量减小，
15、否则就会中止分析。如果增量步的解在少丁五次迭代时就收敛，这表明找到解答相对很容易。因此如果连续两个增量步只需少丁五次的迭代就可以得到收敛解，ABAQUSH动将增量大小提高50%。H动载荷增量方案的详细内容在信息文件中给出，第7.4.2节中有更详细的描述。7.3用ABAQUS进行非线性分析我们现在讨论怎样用ABAQUS分析非线性问题，以儿何非线性问题为主进行讨论。7.3.1几何非线性只需对模型做些小的修改就可以将儿何非线性效应包含丁分析中。首先要在定义分析步时考虑儿何非线性效应，耍给出分析步中允许的最大增量步的数目。如果ABAQUS需耍比此数目更多的增量步來完成分析，它将中止分析并给出出错信息。分析步中默认的增量步数是100,但如果题目有显著的非线性，可能会需要更多的增量步。用户给出的增量步数目是ABAQUS可以采用的增量步数的上限，而不是它所必须使用的增量步数。在非线性分析中，一个分析步是发生于一段有限的“时间”内的，除非惯性效应或率相关效应作为重要因素进入分析，否则这里的“时间”并没有物理含义。用户是在这个理解背景下指定初始时间增量叫门和此分析步的总时间的。这些数据也指定了第一个增量步中所施加的载荷的比例。初始载荷增量如下给定：仏丝
