有限元非线性.pptVIP

8.2非线性问题基本解决思路材料非线性：方程形式不变，将材料本构关系线性化，分段求解，将线性问题的方程推广用于非线性问题。几何非线性：通常采用增量分析法，建立变化位的平衡方程。有两种表达格式：（1）在整个分析过程中参考位保持不变，始终取初始位，称为完全Lagrange格式；（2）在整个分析过程中参考位不断被更新，参考前面每一步荷载步开始的位形，称为修正Lagrange格式。求解方法：直接迭代法(割线刚度法)、N-R法或mN-R法(切线刚度法)、初应力法、初应变法和增量法（切线、初应力、初应变增量，主要用于弹塑性分析）。结构整体平衡方程：（1）假定初始近似解：（2）由本构关系求出（3）由平衡方程求得下一步近似解：（4）重复（2）和（3），直到两次结果非常接近。1、直接迭代法（割线刚度）初始线弹性解RP-δ凸时收敛，凹时可能发散。2、N-R法（切线刚度）010203任何具有一阶导数的连续函数Ψ(x)，在xn点的一阶Taylor展开：非线性方程Ψ(x)=0在xn附近的近似方程是线性方程Newton-Raphson迭代公式每次迭代需要修改K。针对结构平衡方程：Ψ(δ)=[K]{δ}-{R}={F(δ)}-{R}=0利用N-R公式，有：每次迭代不改变它的刚度值始终取初始刚度，计算量小，但收敛慢些。荷载增量法：把荷载分成很多小的荷载步，在每一个荷载步上使用一次或都多次N-R方法。实质上是分段线性化。3、初应力法1如果材料的应力-应变关系可以表示成2即由应变确定应力，设想用具有初应力的线弹性物理方程代替上式：3调整初应力值，可以做到上述两式得到的应力相同。式中：[KT]0为结构的起始切线刚度矩阵，{F}为与初应力等价的节点荷载改写平衡方程弹性解8.4.1材料弹塑性行为弹塑性：卸载后存在不可恢复的残余变形。它与非线性弹性材料有显著区别：加载同，卸载不同。硬化：屈服后应力随应变继续增加；卸载后再加载屈服应力提高，一般等于卸载时的应力。各项同性硬化：反向加载的屈服应力与正向卸载点应力数字上相等。随动硬化：卸载点应力与反向加载的屈服应力绝对值之和等于2倍初始屈服应力。混合硬化：介于上二者。各种硬化塑性特征循环硬松弛循环塑性一般表现循环硬化循环蠕变8.4.2塑性力学的基本法则式中等效应力为米赛斯（VonMises）屈服准则材料在复杂应力状态下的等效应力达到单向拉伸的屈服极限时，材料开始屈服。于是，米赛斯屈服条件可写成：几何上以σ1=σ2=σ3为轴线的圆柱面。0102或用一般应力表示01等效应力还可用应力偏量表示为02式中032、应变强化01假定材料进入屈服后，总应变增量可分成弹性的和塑性两部分02对应于等效应力，定义等效应变为各项同性硬化运动硬化Prager，Zeigler修正对应于塑性应增量的等效应变称为塑性等效应增量，因为塑性变形不产生体积变化，泊松比为0.5，故有01材料进入屈服以后，进行卸载或部分卸载后在加载，新的屈服应力仅与卸载前的等效塑性应变总量有关。新的屈服只有当等效应力适合02才发生。上式为等向硬化材料的米赛斯准则，反映等向强化材料屈服和强化之间的关系。H’为强化阶段的曲线斜率。0301普朗特-路斯（Prandtl-Reuss）塑性流动理论02如果将等向强化米赛斯准则式写成03则F可以看成n维应力空间的一个曲面，称为屈服面。04对于金属一类材料，塑性应变增量和屈服面之间存在如下关系