弹塑性有限元方法.docVIP

第三章 弹塑性有限元方法的实施 §3.1 增量平衡方程和切线刚度矩阵 1、 分段线性化的求解思想 塑性变形的特点决定了塑性本构关系的非线性和多值性，上面由塑性增量理论给出了塑性应力—应变关系 其中 说明当前应力状态不仅与当前应变有关，而且和达到这一变形状态的路径（加载历史）有关。这里包含了屈服准则、强化条件和加卸载准则。 由此，对物理非线性问题，通常采用分段线性化的纯增量法和逐次迭代的方法求解。即将加载过程分成若干个增量步，选择其中任意一个增量步建立它的增量平衡方程并求解，对整个过程的求解有普遍意义。 2、 增量平衡方程和切线刚度矩阵 设t时刻（加载至i-1步终），结构（单元）在当前载荷（广义体力和表面力）的作用下处于平衡状态，此时物体内一点的应力、应变状态为。在此基础上，施加一个载荷增量，即从时刻，则在体内必然引起一个位移增量和相应的、，只要足够小，就有。 倘若初始状态已知，加载过程已知，则可以确定（即可以确定，然后可在硬化曲线上得到所对应的硬化系数）于是上面的方程成为线性的。在这一增量过程中，应用于虚功原理可得到如下虚功方程： （1） 根据小变形几何关系，再由虚位移的任意性，并设，展开后，其中单元在t时刻载荷等效节点力：；内增量载荷的等效力。 这样，由方程（1）可得平衡方程： （2） 即： 因为t时刻（第i步终）结构处于平衡状态 （3） 这样（2）式变为: 即： （4） 将和代入上式得增量平衡方程： （5） 对增量位移求导： （6） 于是（5）式成为 （7） 为单元切向刚度矩阵。集合所有单刚后得到结构总的增量平衡方程 （8） 方程（8）是线性的，可以直接求解。 硬化系数的数值表示 根据单一曲线定理，对于一般稳定性硬化材料，在其简单加载过程中，和之间存在着一一对应的确定的函数关系，这一关系可用单向拉伸实验来确定。 例如，对于Mises各向同性硬化材料 （8） 在有限元分析中，作为初始参数应把这一曲线输入（用函数或数字的形式），在加载过程中弹塑性矩阵不断地修改，根据当前的应力或应变来确定。目前，硬化曲线的输入格式有两种： 解析表达式 根据单一曲线定理，由单向拉伸试验曲线直接得出硬化曲线的解析式。例如： （a）Mises各向同性线性硬化材料 单向拉伸曲线有： 当 当 （9） 则有 （10） -----------------------------------------------------------------------------------------------------------------附：对于一般材料的硬化曲线的求法（求） 如单拉曲线 则硬化曲线 根据 ===》 ===》 其中单拉时等效应变为 因为 ，，平均应变为 所以 ，当 时 （b）Mises各向同性幂硬化材料 单向拉伸曲线有： 当 当 （11） 由屈服点条件：得 据（8）式得 （12） 其中： 根据离散的单拉实验数据，采用样条插值计算 （参看清华大学孟凡中教材：弹塑性有限变形理论和有限元方法） 3.3 过渡单元弹塑性矩阵的确定 三种变形状态 弹塑性变形体中，在一个载荷增量步内可能有三种变形状态： 弹性区：加载前后均处于弹性状态，故采用弹性阵不变。 塑性区：加载前后均处于塑性状态，其弹塑性矩阵由塑性增量理论确定（与当前应力水平和塑性变形增量的总量有关） 过渡区：加载前处于弹性状态，加载后进入塑性状态，所以，在这一过程中采用弹性矩阵或最终的都不合适，必须寻找一个合适的弹塑性矩阵。 加权平均的弹塑性矩阵 过渡单元在加载后的应力计算（以单拉状态为例） 在时间步内施加一个增量载荷后，讨论某单元的应力应变状态。 设某单元加（卸）载前的应力状态，相应的应变（A点）处于弹性状态（弹性区间O’C）。加载后，按弹性计算得到应变增量，到达B点。显然B点不是实际的应力状态，因为已经超过了C点，进入了