塑性理论及有限元要点.ppt

Newton – Raphson迭代法 u o P u1 u2 u3 un 刚度方程 N-R迭代方法每一次迭代都修正其切线刚度，一般情况下收敛性良好，较适合于高度非线性问题的求解。 Taylor转开 切线刚度 修正的N－R迭代法：采用第一次的切线刚度代替每一次的切线刚度进行计算，计算量减少，但收敛性不于N－R法。 u o P u1 u2 u3 un 拟牛顿法：每次采用迭代后用 一个方式修正矩阵[K]，既不像N－R法每次迭代重新计算一个新的矩阵，也不像 mN－R不改变矩阵；计算效率高，而收敛性也好。 增量法：即变刚度法。该方法中，载荷是逐步施加的，每一步计算切线刚度；一般与N-R方法等迭代方法结合适用； 弧长方法：该方法中，增加了作为变量的载荷大小因子λ，附加约束方程进行求解；该方法可用于非线性屈服的极值点附近的分析和加工软化的情况。 u o P u1 u2 u3 un r1 r2 r3 2、线性代数方程组求解方法 高斯消去法－简单，但稳定性差 三角分解法 (改进)平方根法 波前法－适合大型有限元刚度方程求解 可用于不 对称矩阵 用于对称矩阵 直接解法 稀疏、对称、带宽小 迭代解法 线性方程组的迭代解法是建立一种从已知近似解计算得到新的近似解方法，新近似解更接近精确解。该方法能保持矩阵的稀疏性，计算简单，程序编制容易，在很多情况下收敛较快，能得到较高精度得结果。主要有一下几种方法： 雅可比迭代法 高斯－赛德尔迭代法 松弛法－弹性力学分析中松弛因子取1.8左右收敛较快 最速下降法－有时收敛非常缓慢 共轭梯度法－存贮量小、计算简单、充分利用矩 阵得稀疏性，且具有超线性收敛。 十、材料属性 1、屈服准则 Von Mises屈服准则(各向同性) 平面应力问题 平面应变问题 轴对称问题 三维问题主应力状态 三维问题一般应力状态 Hill各向异性屈服准则(主要针对板料问题) 三维问题一般应力状态 为单轴拉伸屈服应力 为单向剪切屈服应力 板料成形中，如认为只存在厚度方向的各向异性，而板平面内性质各方向不变时，取3向为厚度方向，1为长度方向，2为宽度方向，则该准则可作简化。 塑性应变比(或厚向异性系数) 简化的屈服条件 Barlat各向异性屈服准则(主要针对板料问题) S1 , S2 , S3 是对称矩阵Sij 的主值 a, b, c, d, f, g, h 是材料参数, 当都为1时, 材料为各向同性 板料成形中，平面应力问题Barlat各向异性屈服准则 p 根据45 °方向的塑性应γ45变比求得。 2、加工硬化与材料模型 材料加工硬化规律的两种假设 单一曲线假设 ：流动应力是等效应变的函数，与应力状态无关，决定于材料的性质；该假设方便应用，有待进一步证实。 能量条件假设：材料的硬化程度取决于变形过程中的塑性变形功，与应力状态和加载路径无关；该假设具有一般性，较复杂，不便应用。 材料加工硬化类型： 等向强化 随动强化 混合强化 ?2 ?1 ?2 ?1 ?2 ?1 加工硬化模型 Y ? o ?1 ?2 ?3 Y1 Y2 Y3 Y 1 0 0 2 0.01 60 3 0.1 100 4 0.5 400 5 1.0 500 6 5.0 600 一般形式(不含卸载) 采用公式形式表达 采用线图或数表形式表达 多线性随动强化 Y ? o ?1 ?2 2?s ?2 ?1 (?S) ?1 (?S) Y ? o ?1 ?2 2?1 ?2 ?3 ?3 2?2 双线性随动强化 非等向强化 多线性等向强化 双线性等向强化 理想弹(刚)塑性 Y ? o ?xc ?yt ?yc ?xt Y ? o ?1 ?2 ?max ?1 (?S) 2?max ?1 (?S) Y ? o ?1 ?2 ?2 ?3 ?max 2?max Y ? o ?s DYNAFORM软件材料属性的说明： 材料18采用与应变相关的各向同性塑性模型，其硬化采用能量条件假设。 材料24采用各向同性弹塑性模型，应力与应变和应变速率关系可任意定义。 材料36采用平面应力状态下的Barlat三参数平面应力各向异性屈服准则。 材料37只考虑板料厚度方向各向异性Hill屈服条件的弹塑性模型，应力与应变曲线可任意定义。 材料39与材料37类似，但增加了成形极限图的定义。 十一、成形极限图的分析 –40 –30 –20 –10 0 +10 +20 +30 +40 +20 +40 +60 +80 +100 安全区 破裂区 临界区 可能起皱区 起皱区 横向应变ε2(%) 纵向应变ε1(%) 材料硬化指数n 对成形极