第二章有限元法模拟第一节有限元模拟数学原理主要内容有限元发展历史沿革塑性力学知识简单回顾刚塑性有限元基本原理刚粘塑性有限元基本原理有限元简介有限元发展历史40年代萌芽;50年代中期创始;(M.J.Turner和R.W.Clough);60年代诞生;(R.W.Clough,O.C.Zienkiewicz);有限元简介塑性成型有限元发展历史1967年,P.V.Marcal和I.P.King[3]1968年Y.Yamada[4]弹塑性小变形问题;1970年H.D.Hibbit,1975年R.M.McMeeking大变形弹塑性问题;1973年C.H.Lee和S.Kobayashi1979年,O.C.Zienkiewicz刚塑性有限元法;80年代O.C.Zienkiewicz;S.Kobayashi和S.Oh推导了刚粘塑性有限元方法……..有限元简介刚塑性有限元特点优点没有应力的误差累积可以取较大的增量步历史变量少有限元列式简单缺点不适合计算弹性问题不能处理回弹不能处理残余应力不适合计算小变形问题塑性力学知识回顾塑性力学的基本方程-----应力平衡微分方程000zyxzyxzyxzyzxzzyyxyzxyxx塑性力学知识回顾塑性力学的基本方程----小变形几何方程)(21)(21)(21zuxwzwywzvyvxvyuxuxzzxzzyyzyyxxyx)(21)(21)(21zuxwzwywzvyvxvyuxuxzzxzzyyzyyxxyx塑性力学知识回顾塑性力学的基本方程----屈服准则Tresca、Mises准则Tresca、Mises准则s31s塑性力学知识回顾塑性力学的基本方程----本构方程zxzxyzyzxyxyyxzzxzyyzyxx232323)](21[)](21[)](21[zxzxyzyzxyxyyxzzxzyyzyxxdddddddddddd232323)](21[)](21[)](21[塑性力学知识回顾塑性力学的基本方程----上限原理TisiDSijvUisidSuTdSuKdVdSuTTDijU****][TisiDSijvUisidSduTdSuKdVddSduTTDijU****][塑性力学知识回顾基本问题描述1、有一块刚塑性材料,其边界S围成的体积为V内满足平衡方程、协调方程、体积不可压缩条件和本构关系。{基本方程}2、在边界S上,一部分表面上给定位移速度,只一部分表面,给定表面力。{力学边界条件}3、在体积V内可能生热;通过表面S,变形体与工具或与周围环境之间可能发生热交换q,但维持能量平衡。{热力学边界条件}刚塑性有限元基本原理基本问题求解15个方程组成的方程组上限原理;应力平衡微分方程(3)小变形几何方程(6)本构方程(6)(屈服准则)TisiDSijvUisidSuTdSuKdVdSuTTDijU****][刚塑性有限元基本原理数学基础讨论极值的求解方法!TisiDSijvUisidSuTdSuKdVdSuTTDijU****][优化穷举法求导数方法刚塑性有限元基本原理数学基础函数:对于某一值域中的任意数X,F都有唯一的值与之对应,我们称F是一个函数。泛函:对于某一值域中的任意一个函数F(x),II都有唯一的值与之对应,我们称II是一个泛函。实质:泛函是自变量为函数的函数。刚塑性有限元基本原理数学基础函数:F(X)自变量:X微分:df(x)极值:一阶微分为0泛函:II宗量:F(X)变分:&II(F(x))极值:一阶变分为0刚塑性有限元基本原理Markov变分原理设一刚塑性物体.体积为V,边界表面为S,整个物体处于塑性变形状态。在表面的SP部分给定表面力,在表面的SV部分给定速度。在满足几何方程、体积不可压缩条件和速度边界条件的一切动可容速度场,使泛函psiivdsvpdv取得极少值的速度场,即是真实速度场。刚塑性有限元基本原理自然变分原理0Vsiivpdsvpdv刚塑性有限元基本原理广义变分原理----拉格朗日乘子法),(yxfz0),(yxg),(),(yxgyxfz0),(00yxgzygyfyzxgxfxz刚塑性有限元基本原理广义变分原理----拉格朗日乘子法0VsiivpdsvpdvdsvpdVdVpsiivvvm刚塑性有限元基本原理广义变分原理----拉格朗日乘子法m速度场应变速率应力场刚塑性有限元基本原理广义变分原理----罚函数法0VsiivpdsvpdvdsvpdVdVpsiivvv22)(vm刚塑性有限元基本原理广义变分原理----可压缩法认为材料是可压缩的g刚塑性有限元基本原理广义变分原理----各种方法的比较1)拉格朗日乘子法引进了附加的未知数,所以总自由度增加约50%,计算消耗就增加很多。2)虽然三种方法的总刚度矩阵都是对称、稀疏的,但拉格朗日乘子法的刚度阵不呈带状,所以拉格朗日法需要更多的存储空间和计算时间。3)如果罚函数法的罚数选得过大,或可压缩性法的g值选得过小,有可能使刚度方程病态。4)理论上,拉格朗日乘子法更精确刚粘塑性有限元基本思想刚粘塑性的概念刚塑性:忽略变形速度的影响;(低温)粘塑性:考虑变形速度的影响。(高温)),,(T刚粘塑性有限元基本思想刚粘塑性模型1、过应力模型),,(T])(1)[(nrY刚粘塑性有限元基本思想刚粘塑性模型2、Backofen模型),,(Tmc刚粘塑性有限元基本思想刚粘塑性模型3、Rosserd模型),,(Tnmk刚粘塑性有限元基本思想刚粘塑性模型4、Norton-Hoff模型),,(TijmijK1')3(2mmK1)3(刚粘塑性有限元基本思想刚粘塑性模型5、Sellars-Tegart模型),,(TmA)(sinh11刚粘塑性有限元基本思想刚粘塑性有限元列式VijijSiiVdVdSuFdVFVijijSiiVijdVdSuFdVEF)(ddEijijijij00')(刚粘塑性有限元基本思想刚粘塑性有限元列式dVdSuFdVvvSiiVF22)(dVdSuFdVEvvSiiVijF22)()(ddEijijijijij00')(刚粘塑性有限元基本思想刚粘塑性有限元列式])(1)[(nrY])(1)[(])(1)[()(100nnijrnrYdrYdE小结Markov变分原理广义变分原理刚粘塑性变分原理