一维lagrangian和eulerian有限元.pptxVIP

第2章 一维Lagrangian和Eulerian有限元 ;1 引言 ;1 引言 ;两种格式的主要区别在于： TL，在初始构形上描述变量， UL，在当前构形上描述变量。 不同的应力和变形度量分别应用在这两种格式中。 TL，习惯于采用一个应变的完全度量， UL，常常采用应变的率度量。 这些并不是格式的固有特点，在UL中采用应变的完全度量是可能的，并且在TL中可以采用应变的率度量。; 对于每一种公式，将建立动量方程的弱形式，即虚功原理（或虚功）。这种弱形式是通过对变分项与动量方程的乘积进行积分来建立。在TL格式中，积分在所有材料坐标上进行；在Eulerian和UL格式中，积分在空间坐标上进行。将说明如何处理力边界条件，因此近似（试）解不需要满足力边界条件。这个过程与在线性有限元分析中的过程是一致的，在非线性公式中的主要区别是需要定义积分赋值的坐标系和确定选择应力和应变的度量。 ;2 完全的Lagrangian格式;2.2 TL的控制方程;物体的运动由Lagrangian坐标和时间的函数描述;2 完全的Lagrangian格式;2 完全的Lagrangian格式;质量守恒 ;平衡方程 ;2 完全的Lagrangian格式;2 完全的Lagrangian格式;2 完全的Lagrangian格式;; 自然变分原理是对物理问题的微分方程和边界条件建立对应的泛函，使泛函取驻值得到问题的解答，但是其未知场函数需要满足一定的附加条件。 广义变分原理（或称约束变分方程）不需要事先满足附加条件，采用Lagrange乘子法和罚函数法将附加条件引入泛函，重新构造一个修正泛函，将问题转化为求修正泛函的驻值。称为无附加条件的变分原理。 对于罚函数方法，将罚参数取正值，对修正泛函得到的近似解只是近似地满足附加条件，罚参数值越大，附加条件的满足程度就越好。而在实际计算中，罚函数只能取有限值，所以利用罚函数求解只能得到近似解。; 有限元方法不能直接离散动量方程。为了离散这个方程，需要一种弱形式，称为变分形式，即虚功原理或者虚功率，通过对变分项与动量方程的乘积进行积分来建立的。 虚功原理或者弱形式是等价于动量方程和力边界条件的。后者称为经典强形式。 有限元的理论基础是泛函、变分和虚功原理。; 对于??量方程和力边界条件, 现在建立弱形式，要求：;取动量方程与变分项的乘积并在全域内积分得到弱形式，给出 ;弱形式到强形式; 可以看出，如果允许较低平滑的变分项和试函数，在强形式中将附加一个方程－－内部连续条件。如果选取的变分项和试函数满足经典的平滑条件，在强形式中则没有内部连续条件。对于平滑的变分项和试函数，弱形式仅采用动量方程和力边界条件。 较低平滑性要求的变分项和试函数仅是 连续，需要处;2 完全的Lagrangian格式; 弱形式中包含强形式，强形式中包含弱形式，所以弱和强形式是等价的。对于动量方程，强和弱形式的这种等价称为虚功原理。; 以弱形式作为虚功表达式的观点提供了统一性，对于在不同坐标系上和不同类型问题中建立弱形式是很有用途的：为了获得弱形式，只需要写出虚能量方程。因此，可以避免前面所做的由变分项与方程相乘并进行各种处理的过程。 从数学观点来看，没有必要考虑变分函数作为虚位移：它们是简单的变分函数，满足连续条件和在位移边界上为零。对于有限元方程的离散，方程与变分函数的乘积没有物理意义。  建立弱形式中的关键步骤是分部积分，从而消除了关于应力P 的导数。如果没有这一步，力边界条件就不得不强加在试函数上。作为弱形式，由分部积分和降低对应力和试位移函数平滑性的要求（如C0）是更方便的。;3 有限元离散，单元和总体矩阵;3.1 TL的有限元离散;节点力 ;节点力 ;运动方程－－半离散方程 ;运动方程－半离散方程 ; 运动方程在空间是离散的，在时间上是连续的，有时简称离散方程。在有限元离散中，质量矩阵常常为非对角阵(一致质量矩阵)，此时运动方程区别于牛顿第二定律，当MIJ≠0时，节点J处的力可以在节点I处产生加速度。而集中质量矩阵的运动方程等价于牛顿第二定律。;单元和总体矩阵; 2节点单元一维网格的集合和离散运算的描述，两组单元节点位移的集合：位移根据单元节点编号集合；计算节点力的离散：节点力根据节点编号返回总体力矩阵。; 2节点单元一维网格的单元形函数Ne(X)和总体形函数N(X);3 有限元离散，单元和总体矩阵;例题2.1 ;3 有限元离散，单元和总体矩阵;3 有限元离散，单元和总体矩阵;3 有限元离散，单元和总体矩阵;4 更新的Lagrangian格式;4.1 UL的控制方程;4 更新的Lagrangian