几何非线性与屈曲分析
非线性屈曲分析
ansys workbench非线性屈曲分析(2013-08-26 21:26:29)转载▼标签：ansys很多旋转受压结构必须进行屈曲分析，常规结构屈曲分析软件有nastran、abaqus和ansys，nastran对线性大型模型分析效率较高；abaqus屈曲分析使用较少；ansys使用比较频繁，其快速建模，与CAD软件的良好借口及有限元模型前处理的便捷性（WB界面）很有吸引力，屈曲分析功能较为完善，可以进行线性、非线性和后屈曲分析。
ansys学习资料中介绍较多的是线性屈曲分析。
线性屈曲分析在工业实际中预测的值偏高，有的甚至超过实际实验测试值的几十倍，线性分析唯一优势是其分析速度较快。
但在实际中其预测值参考价值不大，仅给定结构屈曲失效的上限值。
非线性屈曲分析考虑其他因素，包括结构加工缺陷（几何），材料非线性等，因此较为接近实际情况，但计算耗时较长。
针对最艰难学习情况归纳总结非线性屈曲分析时技术要点及应注意事项。
对于规则旋转壳，承受外压载荷作用，进行非线性屈曲分析时，必须加上几何缺陷，关键步是添加APDL语句/prep7upgeom,0.1,1,1,file,rstcdwrite,db,file,cdb/solu该步引入屈曲模态情况下的几何缺陷，缺陷为屈曲模态变形相对值的0.1倍，该值可以根据实际加工水平等其他条件确定，上述语句保存在txt文档中，在workbench流程APDL模块调用。
分析详细流程为，static structure模块导入几何，施加载荷和边界条件，分析求解，将linear buckling拖入流程中，共享static structure模块数据，进行线性屈曲模块分析，Mechanial APDL模块调用屈曲分析结果，并调入（addinput）上面内含几何缺陷命令语句命令的txt文件，更新，将Mechanical结果导入Finite Element modeler模块，更新，此时在缺陷附近的单元节点位置发生改变。
屈曲分析
问题概述：
一个可靠的产品设计，不仅强度要满足设计要求，而且结构要有足够 的刚度来保证产品性能。现代电子产品（其他产品也一样）已经越来越小， 电子元件之间的空隙非常狭小，在刚度不够的情况下往往导致零部件之间 的干涉。例如，手机从高处坠落，有可能会出现摔坏的情况，这可能是因 为外壳变形过大破坏了内部结构。通常，运行一个静态分析就可以得到结 构在载荷作用下的变形。在某些结构，如承受压应力的部件，在压力载荷 到达一定程度以后会发生于静态分析相比大的多的不可思议的变形，这就 是由于结构已经在这一载荷作用下发生了失稳，这时就需要稳定性分析即 屈曲分析。 实际上结构发生失稳也是由于应力刚度矩阵在影响，应力刚度矩阵可 以加强或减弱结构刚度，这与应力是拉应力还是压应力有关。正如前面计 算出的结果一样，拉应力会使结构的横向刚度增强；结构受压时，会导致 结构的刚度减弱，当压力越来越大时，刚度弱化超出了结构固有的刚度， 结构就表现的很脆弱，位移急剧增大，发生屈曲。
半径i的数值就能使 减小。可见，如果不增加截面面积，尽可能的把材料放在离截 面形心较远处以取得较大的I和i值，就能提高临界应力。
改进措施：在相同截面积下，将杆的结构改为空心杆，截面见图1.其余 各项设置与实心杆相同，计算出的BLF值和实心杆的BLF对比，见图2，失稳 临界载荷因子有明显提高，说明上述分析是正确的。
图1 空心杆截面图
图2 实心杆与空心杆BLF值对比
屈曲分析示例
一端固定一端自由的薄壁圆筒屈曲模态振型
屈曲分析示例
细长圆杆失稳分析及改进
细长圆杆如下图，直径15mm，长200mm，一端固定，一端自由，且受 到100N的压力作用，进行失稳分析并改进。
圆杆三维模型
各阶失稳临界载荷因子（BLF） 失稳的屈曲模态振型（10阶）
线性屈曲与非线性屈曲
1．线性屈曲特征值屈曲分析属于线性分析，它对结构临界失稳力的预测往往要高于结构实际的临界失稳力，特载值分析得到的是第一类稳定问题的解，只能得到屈曲荷载和相应的失稳模态，它的优点就是分析简单，计算速度快。
事实上在实际工程中应用还是比较多的，比如分析大型结构的温度荷载，而且钢结构设计手册中的很多结果都是基于特征值分析的结果，例如钢梁稳定计算的稳定系数，框架柱的计算长度等。
它的缺点主要是：不能得到屈曲后路径，不能考虑初始缺陷如初始的变形和应力状态，不能考虑材料的非线性。
但特征值屈曲分析作为非线性屈曲分析的初步评估作用是非常有用的。
特征值屈曲分析所预测的结果我们只取最小的第一阶，所以你所得出的特征值临界失稳力的大小应为F＝实际施加力*第一价频率。
2. 非线性屈曲非线性屈曲分析的第一步最好进行特征值屈曲分析，特征值屈曲分析能够预测临界失稳力的大致所在，因此在做非线性屈曲分析时所加力的大小便有了依据。
由于非线性屈曲分析要求结构是不“完善”的，比如一个细长杆，一端固定，一端施加轴向压力。
若次细长杆在初始时没有发生轻微的侧向弯曲，或者侧向施加一微小力使其发生轻微的侧向挠动。
那么非线性屈曲分析是没有办法完成的，为了使结构变得不完善，你可以在侧向施加一微小力。
这里由于前面做了特征值屈曲分析，所以你可以取第一阶振型的变形结果，并作一下变形缩放，不使初始变形过于严重，上步完成后，加载计算所得的临界失稳力，打开大变形选项开关，采用弧长法计算，设置好子步数，计算。
后处理，主要是看节点位移和节点反作用力（力矩）的变化关系，找出节点位移突变时反作用力的大小，然后进行必要的分析处理。
非线性分析比较好的是能够得到结构和构件的屈曲后特性，可以考虑初始缺陷还有材料的非线性包括边界的非线性性能。
但是在分析的时候最好是在线性特征值的基础上，因为这种方法的结果依赖所加的初始缺陷，如果所加的几何缺陷不是最低阶，可能得到高阶的失稳模态。
midas几何非线性理论知识 (1)
当结构的变形相对杆件长度已不能忽略时，为了在结构变形后的形状上建立平衡，并考虑初始缺陷对结构屈曲承载力的影响，必须对结构进行基于大挠度理论的非线性屈曲分析。
在midas中可以这样处理：对于索结构或张悬梁结构中，定义的只受拉索单元并不能进行特征值分析，因为其只能定义在几何非线性分析中。
如要进行特征值分析，那么要将只受拉索单元转换为只受拉桁架单元。
先对该结构进行几何非线性，得出自重作用下的初始索力，然后将索单元定义为只受拉桁架单元，将计算所得的索力按初始荷载加到单元中：荷载－>初始荷载－>小位移－>初始单元内力加入张力。
1、问：在MIDAS 中如何计算自重作用下活荷载的稳定系数(屈曲分析安全系数)?答：稳定分析又叫屈曲分析，所谓的荷载安全系数(临界荷载系数)均是对应于某种荷载工况或荷载组合的。
例如：当有自重W 和集中活荷载P 作用时，屈曲分析结果临界荷载系数为10 的话，表示在10*(W+P)大小的荷载作用下结构可能发生屈曲。
但这也许并不是我们想要的结果。
我们想知道的是在自重(或自重+二期恒载)存在的情况下，多大的活荷载作用下会发生失稳，即想知道W+Scale*P 中的Scale 值。
我们推荐下列反复计算的方法。
步骤一：先按W+P 计算屈曲分析，如果得到临街荷载系数S1。
步骤二：按W+S1*P 计算屈曲，得临界荷载系数S2。
步骤二：按W+S1*S2*P 计算屈曲，得临界荷载系数S3。
重复上述步骤，直到临街荷载系数接近于1.0，此时的S1*S2*S3*Sn 即为活荷载的最终临界荷载系数。
(参见下图)midas官方网站的说话，供大家参考：考虑几何非线性同时进行稳定分析可以实现。
方法如下：1、将进行稳定分析所用荷载定义在一个荷载工况下；2、定义非线性分析控制，选择几何非线性，在非线性分析荷载工况中添加此荷载工况，并对其定义加载步骤；3、分析；4、查看结果中的阶段步骤时程图表，查找变形发生突变的位置点，及加载系数，即可推知发生失稳的极限荷载。
屈曲分析
屈曲（失稳）征值屈曲分析与非线性屈曲分析：很多现有的ANSYS资料都对特征值屈曲分析进行了较为详细的解释，特征值屈曲分析属于线性分析，它对结构临界失稳力的预测往往要高于结构实际的临界失稳力，因此在实际的工程结构分析时一般不用特征值屈曲分析。
但特征值屈曲分析作为非线性屈曲分析的初步评估作用是非常有用的。
以下是我经过多次计算得出的一些分析经验，欢迎批评。
1. 非线性屈曲分析的第一步最好进行特征值屈曲分析，特征值屈曲分析能够预测临界失稳力的大致所在，因此在做非线性屈曲分析时所加力的大小便有了依据。
特征值屈曲分析想必大家都熟练的不行了，所以小弟不再罗嗦。
小弟只说明一点，特征值屈曲分析所预测的结果我们只取最小的第一阶，所以你所得出的特征值临界失稳力的大小应为F＝实际施加力*第一价频率。
2. 由于非线性屈曲分析要求结构是不“完善”的，比如一个细长杆，一端固定，一端施加轴向压力。
若次细长杆在初始时没有发生轻微的侧向弯曲，或者侧向施加一微小力使其发生轻微的侧向挠动。
那么非线性屈曲分析是没有办法完成的，为了使结构变得不完善，你可以在侧向施加一微小力。
这里由于前面做了特征值屈曲分析，所以你可以取第一阶振型的变形结果，并作一下变形缩放，不使初始变形过于严重，这步可以在Main Menu> Preprocessor> Modeling> Update Geom中完成。
3. 上步完成后，加载计算所得的临界失稳力，打开大变形选项开关，采用弧长法计算，设置好子步数，计算。
4. 后处理，主要是看节点位移和节点反作用力（力矩）的变化关系，找出节点位移突变时反作用力的大小，然后进行必要的分析处理。
屈曲的特征理解：当结构轴向（梁，板，壳）承受压缩载荷作用时，若压缩载荷在临界载荷以内，给结构一个横向干扰，结构就会发生挠曲，但当这个横向载荷消除时，结构还会恢复到原有的平衡状态，此时杆的直的形式的弹性平衡是稳定的。
屈曲分析分析原理
屈曲分析分析原理屈曲分析原理字数 765预计阅读时间 5min1、小位移和大位移小位移：在利用欧拉公式计算时，属于线弹性计算，忽略了结构的变形对结构的影响，结构的刚度矩阵是不变的。
而实际上，结构的变形是可以影响荷载的作用效应的。
如下图所示。
对杆件施加一定的荷载后，杆件会产生相应的变形，在这个变形的基础上，荷载会继续作用在这个（刚度矩阵）已经改变的杆件上从而导致二阶变形。
为了更好理解，我用银行利息的例子比喻一下这个现象。
比如我拿一万元钱作为荷载，施加到银行这个杆件上，那么它会产生相应的利息。
之后我这个本金加利息的基础上再次对银行施加荷载以获取进一步的利息。
这就是大位移：几何非线性的，考虑了结构变形的影响。
小位移和大位移的计算公式：2、几何刚度在大位移计算中，考虑了结构变形对荷载作用效应的影响，也就是结构刚度的改变，于是引入几何刚度的概念。
同样用一个比喻来帮助大家理解几何刚度的概念，就是拔河。
在大家的感性认识中，绳子在张紧（受拉）状态下的刚度是不是要比松弛（不受力）状态下的刚度大呢？而实际上，绳子的弹性刚度是没有改变的，所以随着外力的改变，我们引入几何刚度来描述这一现象。
3、计算原理Midas的线性屈曲分析可计算包含桁架单元、梁单元、板单元、实体单元的结构的临界荷载系数和相应的屈曲模态。
结构的静力平衡方程如下：结构的几何刚度矩阵由各单元的几何刚度矩阵构成，各单元的几何刚度矩阵与构件的内力相关。
将几何刚度矩阵用临界荷载系数与使用初始荷载计算的几何刚度矩阵的乘积表示如下：上述平衡方程失稳的条件是存在奇异解，即等效刚度矩阵的行列式的值为零。
即线性屈曲分析就是解下式的特征值，屈曲分析中的特征值就是临界荷载系数。
所谓临界荷载就是初始荷载乘以临界荷载系数的荷载值，表示结构作用临界荷载时结构会发生屈曲(失稳)。
结构失稳时常伴随大位移变形和材料屈服，所以屈曲分析常要求考虑几何非线性线或材料非线性。
ANSYS结构非线性分析指南(一至三章)
ANSYS结构非线性分析指南(一到三章)屈服准则概念：1.理想弹性材料物体发生弹性变形时，应力与应变完全成线性关系，并可假定它从弹性变形过渡到塑性变形是突然的。
2.理想塑性材料（又称全塑性材料）材料发生塑性变形时不产生硬化的材料，这种材料在进入塑性状态之后，应力不再增加，也即在中性载荷时即可连续产生塑性变形。
3.弹塑性材料在研究材料塑性变形时，需要考虑塑性变形之前的弹性变形的材料这里可分两种情况：Ⅰ.理想弹塑性材料在塑性变形时，需要考虑塑性变形之前的弹性变形，而不考虑硬化的材料，也即材料进入塑性状态后，应力不再增加可连续产生塑性变形。
Ⅱ.弹塑性硬化材料在塑性变形时，既要考虑塑性变形之前的弹性变形，又要考虑加工硬化的材料，这种材料在进入塑性状态后，如应力保持不变，则不能进一步变形。
只有在应力不断增加，也即在加载条件下才能连续产生塑性变形。
4.刚塑性材料在研究塑性变形时不考虑塑性变形之前的弹性变形。
这又可分两种情况：Ⅰ.理想刚塑性材料在研究塑性变形时，既不考虑弹性变形，又不考虑变形过程中的加工硬化的材料。
Ⅱ.刚塑性硬化材料在研究塑性变形时，不考虑塑性变形之前的弹性变形，但需要考虑变形过程中的加工硬化材料。
屈服准则的条件：1.受力物体内质点处于单向应力状态时，只要单向应力大到材料的屈服点时，则该质点开始由弹性状态进入塑性状态，即处于屈服。
2.受力物体内质点处于多向应力状态时，必须同时考虑所有的应力分量。
在一定的变形条件（变形温度、变形速度等）下，只有当各应力分量之间符合一定关系时，质点才开始进入塑性状态，这种关系称为屈服准则，也称塑性条件。
它是描述受力物体中不同应力状态下的质点进入塑性状态并使塑性变形继续进行所必须遵守的力学条件，这种力学条件一般可表示为）＝Cf（σij又称为屈服函数，式中C是与材料性质有关而与应力状态无关的常数，可通过试验求得。
屈服准则是求解塑性成形问题必要的补充方程。
1.1 什么是结构非线性在日常生活中，经常会遇到结构非线性。
abaqus屈曲载荷因子
Abaqus屈曲载荷因子是用于衡量结构在受力过程中发生屈曲的临界载荷大小的参数。
屈曲分析主要是为了确定支撑柱、承重梁等受压结构的屈曲模态和临界载荷。
在Abaqus中，可以通过线性屈曲分析和非线性屈曲分析来计算屈曲载荷因子。
线性屈曲分析主要用于预估临界失稳载荷和失稳模态。
在Abaqus中，进行线性屈曲分析的方法是通过Buckle进行的。
一般线性屈曲分析只需要关注第一阶屈曲模态，并根据计算所得的第一阶屈曲载荷因子预估使结构发生屈曲所需要的临界载荷是多大。
但通常而言，线性屈曲分析得到的临界失稳载荷可能会大于结构实际的承载能力。
非线性屈曲分析则可以更准确地模拟结构在屈曲过程中的实际响应。
在Abaqus中，对于较复杂的模型，可以考虑使用Riks法进行后屈曲计算，从而获取屈曲后的结构响应情况。
此外，对于涉及接触脱开等特别复杂的问题，可能需要借助Explicit来实现；而对于局部褶皱问题，需要借助Static、Stabilize来实现。
总的来说，Abaqus屈曲载荷因子是用于衡量结构在受力过程中发生屈曲的临界载荷大小的参数。
在实际工程应用中，需要根据结构的具体情况选择合适的屈曲分析方法，并关注屈曲载荷因子的计算结果，以确保结构的安全性能。
更深入地探讨其在工程应用中的价值以及屈曲分析在结构设计中的重要性。
1. 工程应用价值：屈曲载荷因子是衡量结构屈曲性能的重要参数，可以用于评估结构在受力过程中是否容易发生屈曲。
在工程实践中，通过计算屈曲载荷因子，可以帮助工程师发现潜在的屈曲问题，从而优化结构设计，提高结构的安全性能和稳定性。
2. 屈曲分析的重要性：屈曲分析对于结构设计至关重要，尤其是在承受高压、重载或复杂应力状态的结构中。
屈曲分析可以揭示结构在受力过程中可能出现的问题，如局部或整体失稳、褶皱等。
这些问题可能会导致结构性能下降、甚至失效，因此，在结构设计过程中进行屈曲分析是确保结构安全、可靠和高效的重要手段。
3. 非线性屈曲分析：在实际工程应用中，许多结构的受力特性呈现非线性，因此，进行非线性屈曲分析是十分必要的。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
33
主讲：练章华 教授
4.应力刚化
结构的面外刚度可能极大地受到结构中面内应力状态的影响。面内 应力和横向刚度之间的耦合，称为应力刚化
P
尽管应力刚化理论 假定单元的转动和 应变是小的。
刚化应力仅可以通过 进行大挠度分析得到 P
刚化应力可采用小挠 度或线性理论得到
它在薄的、高应力的结构中， 如缆索或薄膜中，是最明显的。 一个鼓面，当它绷紧时会产生 垂向刚度，这是应力刚化结构 的一个普通的例子。
Lzh_CAE
39
主讲：练章华 教授
/PREP7 /TITLE, BUCKLING …. ET , 1 , BEAM3 R , 1 , 0.25 , 52083E-7 , 0.5 MP , EX , 1 , 30E6 N,1 N , 11 , , 100 FILL E,1,2 EGEN , 10 , 1 , 1 FINISH Save , buckle1 , db Resume , buckle1 , db /SOLU ANTYPE , STATIC PSTRES ,ON D, 1 , ALL
Lzh_CAE
32
主讲：练章华 教授
3. 具有小应变的大挠度（转角）
在一个大挠度分析中，单元的转动可 以任意地大，但是应变假定是小的。
大挠度效应
用NLGEOM，ON命令来激活那些支 持这一特性的单元中的大挠度效应。
GUI：Main Menu＞Solution＞Analysis Options
Lzh_CAE
应力硬化梁
34
Lzh_CAE
5．旋转软化
在小位移分析中这种调整近 似于由于大的环形运动而导 致几何形状改变的效应
主讲：练章华 教授
旋转软化是指动态质量效应调整（软化）旋转物体的刚度矩阵 通常它和预应力[PSTRES] 一起使用，这种预应力由旋 转物体中的离心力所产生。 它不应和其它变形非线性、 大挠度和大应变一起使用 旋转软化用OMEGA命 令中的KPSIN来激活 GUI：Main Menu>Preprocessor>Load>-Load-Apply >-Structura1-Other>Angular Velocity
分叉临界点的载荷－位移平衡路径
极值临界点的载荷－位移平衡路径 Lzh_CAE
10
主讲：练章华 教授
1．线性稳定分析
如果失稳前结构处于小变形状态，可以不考虑几何非线性对平衡方程和几 何方程的影响。如同时假定材料仍处于弹性状态，则失稳前可采用线弹性 分析来求解结构内的位移和应力。
t
ut pu
t
t p
1．大应变效应 2．关于大应变的特殊建模讨论 3. 具有小应变的大挠度（转角） 4．应力刚化 5．旋转软化
Lzh_CAE
29
主讲：练章华 教授
1．大应变效应
一个结构的 总刚度 依赖于 它的组成部件 单元的方向 和单元的刚
单元的形 状改变 大变形 单元刚度 将改变
单元的取 向改变
大转角
当一个单元的节点经历 位移后，那个单元对总 体结构刚度的贡献可以 以两种方式改变。
约束
Lzh_CAE
allsel time , 8e-5 autots , 1 deltim , 4.4e-7 kbc , 1 outres , all , 4 solve fini /postl set , last top_node = 26 * get , deform , node , top_node , u ,y /dscale , 1 , 1 pldisp , 2 plnsol , epto , eqv fini /post26 nsol , 2 , top_node , u , y , displace deriv , 3 , 2 , 1 , , velocity , , , 1 /axlab , y , velocity plvar , 3 fini 21
来自非线性静态分析的结果主要由位移、 应力、应变以及反作用力组成。可以用 POST1 或者用POST26来检查这些结果。
Lzh_CAE
19
主讲：练章华 教授
非线性瞬态分析实例
一个铜弹以给定的速度射向壁面，壁面假定是刚性的 和无摩擦的。目的是研究弹和壁面接触后80s内的现 象，即弹的整个变形，速度历程以及最大等效Von Mises应变。求解使用SI单位系统。
t
t t
p p)
t
12
Lzh_CAE
主讲：练章华 教授
三、非线性结构分析简介 1．非线性结构
订书机
木架 ①几何非线性 ②材料非线性 ③状态变化(包括接触)
轮胎
非 线 性 Lzh_CAE
13
主讲：练章华 教授
2．非线性分析的特殊性 （1）逐步递增载荷和平衡迭代
非线性问题需要 一系列带校正的 线性近似来求解
(1).非线性 屈曲分析 (2).特征值 屈曲分析
非线性屈曲分析比线性屈曲分析更精确，故 建议用于对实际结构进行的设计或估计中。
特征值屈曲分析用于预测一个理想弹性结构的理论 屈曲强度（歧点）。特征值屈曲分析经常产生非保 守结果，通常不能用于实际现实生活中的工程分析。
特征值（或线性）屈曲分析
Lzh_CAE
36
主讲：练章华 教授
Lzh_CAE
37
主讲：练章华 教授
Lzh_CAE
38
主讲：练章华 教授
特征值屈曲分析实例
如图5一15所示，一根长为L，两端铰支的细 长杆受到轴向加载。此杆截面的高度为h，面 积为A。由于对称性，只给杆的上端建模，则 上半部分的边界条件变为自由一固支。为了 描述屈曲模态，在x方向取10个主自由度。杆 的惯性矩为 I=Ah2／12＝0.0052083 in4详细参数如下： 材料性质： 模型尺寸： 载荷。 EX=30E6 psi L=200 in F=1 lbf A=0.25 in2 h＝0.5 in
主讲：练章华 教授
合位移m
等效应变
Lzh_CAE
22
主讲：练章华 教授
X方向位移
等效应力Pa
Lzh_CAE
23
主讲：练章华 教授
(s)
0 -50
0.E+00 2.E-05 4.E-05 6.E-05 8.E-05 1.E-04
(m/s)
-100 -150 -200 -250
Lzh_CAE
24
主讲：练章华 教授
35
GUI：Main Menu＞Solution ＞Analysis Options
Lzh_CAE
主讲：练章华 教授
六、 屈曲分析 1.概念
屈曲分析是一种用于确定结构开始变得不 稳定时的临界载荷和屈曲模态形状（结构 发生屈曲响应的特征形状）的技术。
ANSYS提供两种分析结构屈曲载荷和屈曲模态的技术
11
10
10 9 9 8 8 7 7 6 6 5 5 4 4 3 3非线性问题与屈曲分析
一、几何非线性问题基本理论 二、几何非线性有限元方程的建立 三、非线性结构分析简介 四、非线性分析的基本过程 五、几何非线性 六、屈曲分析 七、材料非线性
Lzh_CAE
1
主讲：练章华 教授
Lzh_CAE
2
主讲：练章华 教授
1 有限应变与应力
(1) 格林(Green)应变，也称格林－拉格朗日 (Green－Largrange)应变
大转动 Lzh_CAE
它的局部刚度转化到全 局部件的变换也将改变
30
主讲：练章华 教授
2．关于大应变的特殊建模讨论 （1）应力－应变
在大应变求解中，所有应力－应变输入和结 果将采用真实应力和真实（或对数）应变。
对真实应变超过50％的塑性分析，应使用 大应变单元（VISCO106，107及108）
Lzh_CAE
u Ke1Q
DBL0u
Lzh_CAE
11
主讲：练章华 教授
2. 非线性稳定分析
此类分析是指失稳前结构处于大变形状态，这时结构的刚度矩阵是 载荷幅值p和位移向量u的非线性函数。
tcr t t t t K K ( K K) t
t
SeeP128-130
pcr p 1 (
rectng , 0 , 0.0032 , 0 , 0.0324 lesize , 2 , , , 20 , 1 lesize , 1 , , , 4 , 1 建模
mshape , 0 , 2d mshkey , 1 amesh , 1 网格 fini /solu antype , 4 trnopt , full lumpm , on nlgeom , 1 ic , all , uy , ,-227 d , all ,, 0 , , , , uz nsel , s , loc , x , 0 d,all , ,0, , , ,ux nsel , s , loc , y , 0 d ,all , , 0 , , , ,uy
（1）子步数或时间步长 （2）自动时间步长
以载荷增量加 载，程序在每 一步中进行平 衡迭代
Lzh_CAE
18
主讲：练章华 教授
四、非线性分析的基本过程
非线性分析的步骤
非线性分析建模过程中可能包含 特殊的单元或非线性材料性质。 除此之外，建模过程与线性分析 类似。
1．建模 2．加载并求解 3．检查结果
31
主讲：练章华 教授
(2) 单元的形状
变形前的单元
变形后的大应变单元 (不好的边长比)
变形后的大应变单元 (较好的边长比) 修改后变形前的单元
在大变形分析的任何迭代中粗劣的单元形状（即大的横纵比， 过度的顶角及具有负面积的已扭曲单元）是不好的，所以必须 象注意单元的原始形状一样注意单元已扭曲后的形状。