材料非线性有限元分析.ppt

2000.4 哈尔滨建筑大学 王焕定教授制作 材料非线性有限元分析 什么是塑性 塑性的特性 塑性的特性 弹塑性问题的有限单元法－控制方程 弹塑性问题的有限单元法－增量法 载荷的保守性 平衡迭代 弹性极限 屈服下限 屈服上限 强度极限 强化段 软化段 弹性变形 残余变形 卸载 加、卸载准则 在单向拉伸状态下，J2=σ2/3。在纯剪状态下，J2=τ2。 等向强化的米塞斯屈服准则 在主应力状态下，第二不变量为 由薄壁圆筒的实验研究可得，这种材料的屈服面方程为 式中J2是应力偏张量的第二不变量 ，由于偏张量第一不变量等于零 ，因此 在单向拉伸状态下，J2=σ2/3。在纯剪状态下，J2=τ2。 屈服面式中χ(k)，是由单向应力状态的数据确定的屈服参数。在单向拉伸时为χ2=σB2/3。在纯剪状态下χ=τB。任何情况下χ都是硬化参数塑性功wp的函数。 随动强化的米塞斯屈服准则 　这种材料的屈服面方程为 在塑性力学中，认为材料进入塑性后存在一个势函数（简称塑性势） 。塑性应变增量可由势函数给出： 流动准则又可分为正交（相关）流动准则和非正交（非相关）流动准则两种。前者认为塑性势就是屈服面，因此 。而后者则认为塑性势和屈服面不同。对正交准则，塑性流动方向垂直于屈服面，加、卸载准则取决于非负的尺度因子dλ，它大于零，表示加载，等于零，表示其他情况。 塑性状态的加载和卸载准则 在外部作用下应变点仍在屈服面上，并有新的塑性变形发生，此时称这个过程为塑性加载。 如果应变点离开屈服面退回弹性区，反应是纯弹性的，此过程称塑性卸载。 应变点不离开屈服面，又无新的塑性变形发生，此时称中性变载。 卸载，弹性 加载，塑性 中性变载，塑性 其几何解释为：弹性应力增量指向屈服面内侧或相切时，反应是弹性的。否则是塑性加载，反应是弹塑性的。 塑性加载 塑性卸载 中性变载 对弹塑性问题，本构关系为 弹性矩阵 塑性矩阵 弹塑性矩阵 在应力增量dσij作用下，应变增量dεij 可分成弹性和塑性两部分。 总应变为 正交（相关）流动准则 弹塑性问题本构关系－应力应变关系 非线性的基本特征是结构刚度的改变 以强化的加载状态为例，因为屈服面为 由df =0 确定弹塑性矩阵－等向强化的米塞斯屈服准则 在主应力状态下，第二不变量为 由薄壁圆筒的实验研究可得，这种材料的屈服面方程为 式中J2是应力偏张量的第二不变量， 由于偏张量第一不变量＝0 ，因此 在单向拉伸状态下，J2=σ2/3。在纯剪状态下，J2=τ2。一般情况下，sij=σij -σkkδij /3,所以 . 根据屈服面表达式，可求得 因为 为了求A，需先由屈服面对wp的偏导数求M 式中Gp是 曲线的斜率。同理，对单向拉伸情况， -M=Ep/3。 纯剪 Gp是塑性剪切模量 单向拉伸 Ep是塑性拉伸模量 J2=τ2 由此可得A=G+Gp（或A=G+Ep/3），又因 由此可得弹塑性矩阵为 象线性问题一样，设位移和应变分别为 则全量形式的应力为 增量形式的应力为 涉及路径相关性的材料应力应变非线性（加载及卸载、残余塑性变形等）必须用增量法来求解。 在增量荷载ΔRm作用下，位移、应力、应变和内变量等的增量分别为 下一步迭代时的荷载水平为 设m迭代步的结果已知，位移、应力、应变和内变量等分别记作 2000.4 哈尔滨建筑大学 王焕定教授制作 虚位移原理 1. 虚位移原理和最小势能原理 1） 虚位移原理的虚功方程——矩阵表达 δWe=∫V[Fb]Tδ[u]dV+ ∫Sσ[Fs]Tδ[u]dS =δWi=∫V[σ]Tδ[ε]dV 体积力虚功 表面力虚功 虚变形功 δWe=∫VFbiδuidV+ ∫SσFsiδuidS =δWi=∫VσijδεijdV 虚功方程——张量表达 2） 势能原理的数学表达 Ve=Vε+VP =1/2∫VσijεijdV -∫VFbiuidV- ∫SσFsiuidS = min 总势能 应变能 外力势能 象线性问题一样，设位移和应变分别为 也即单元增量应变为 。象弹性问题一样，建立第m+1步单元刚度方程为 δWe=∫VFbiδuidV+ ∫SσFsiδuidS =δWi=∫VσijδεijdV 精确解时 弹塑性问题非线性方程组求解问题。 将其代入单