清华大学计算固体力学第二次课件一维Lagrangian和Eulerian有限元.ppt

运动方程－半离散方程(矩阵形式) 运动方程在空间是离散的，在时间上是连续的，有时简称离散方程。在有限元离散中，质量矩阵常常为非对角阵(一致质量矩阵)，运动方程区别于牛顿第二定律，当MIJ≠0时，节点I处的力可以在节点J处产生加速度。而集中质量矩阵的运动方程等价于牛顿第二定律。 3 有限元离散，单元和总体矩阵 为在质点I上的净力。由牛顿第三定律，作用在节点上的力大小相等，而方向相反，因此内部节点力需要一个负号。 单元和总体矩阵 在有限元程序中，通常以一个单元水平计算节点力和质量矩阵，将单元节点力结合入总体矩阵，称为离散或矢量组合。 组合单元的质量矩阵和其它方阵到总体矩阵，称为矩阵装配。 通过计算可以从总体矩阵中提取单元节点位移，称为集合。 3 有限元离散，单元和总体矩阵 2节点单元一维网格的集合和离散运算的描述，两组单元节点位移的集合：位移根据单元节点编号集合；计算节点力的离散：节点力根据节点编号返回总体力矩阵。 3 有限元离散，单元和总体矩阵 2节点单元一维网格的单元形函数Ne(X)和总体形函数N(X) 3 有限元离散，单元和总体矩阵 单元节点位移与总体节点位移的关系为 Le为连接矩阵。类似的获得单元节点力。 应用连接矩阵还可以建立单元形函数和总体形函数之间的关系，总体位移场可以由所有单元的位移求和得到： 对单元形函数求和得到总体形函数 3 有限元离散，单元和总体矩阵 例题 3 有限元离散，单元和总体矩阵 3 有限元离散，单元和总体矩阵 3 有限元离散，单元和总体矩阵 3 有限元离散，单元和总体矩阵 4 更新的Lagrangian格式 4.1 UL的控制方程 初始构形 参考构形 当前构形 变形构形 4 更新的Lagrangian格式的控制方程, 弱形式，单元方程 4 更新的Lagrangian格式的控制方程, 弱形式，单元方程 UL格式是TL格式的一个简单转换。在数值上，离散方程是相同的，而实际在同一程序中，对某些节点力我们可以应用TL格式，而对其它的节点力应用UL格式。 为什么采用两种方法，而它们基本上是一致的。 4.1 UL的控制方程 主要原因是它们都在被广泛地应用，因此，为了理解程序和文献，有必要熟悉两种格式。 应变的度量由变形率给出 应力的度量 Cauchy 应力 4 UL格式的控制方程，弱形式，单元方程 4.1 UL的控制方程 以Eulerian坐标表述相关变量，空间坐标 速度应变 UL格式的两个相关变量－速度和Cauchy应力 质量守恒 对于杆 动量守恒 4 UL格式的控制方程，弱形式，单元方程 能量守恒 －热流量 －热源 本构方程 变形度量 4.1 UL的控制方程 边界条件 速度边界等价位移边界 力边界 n 单位法线（＋,－） 一端固定一端自由杆 边界条件满足 初始条件 4 UL格式的控制方程，弱形式，单元方程 4.1 UL的控制方程 4 UL格式的控制方程，弱形式，单元方程 4.2 UL的弱形式 由动量方程乘以变分函数 弱形式－虚功率原理 强形式－虚功率原理的逆过程：动量方程， 力边界条件 内部连续条件 积分在当前域上完成 4 UL格式的控制方程，弱形式，单元方程 4.2 UL的弱形式 内部虚功率 外力虚功率 惯性力虚功率 弱形式 4 UL格式的控制方程，弱形式，单元方程 4.3 UL的单元方程 在一个单元的水平上建立方程，通过装配获得总体方程。相关变量为速度和应力。 建立本构方程、质量守恒方程，动量方程。由于质量守恒是一个代数方程，可以容易地计算任意一点的密度。建立半离散方程。 单元的速度场为 单元的加速度场为 将形函数表示成为材料坐标的函数是非常关键的，它与时间无关。如果将形函数由Eulerian坐标表示为 形函数的材料时间导数不为零（注意与TL区别），并且不能将加速度表示为同样形函数与节点加速度乘积的形式。 4 UL格式的控制方程，弱形式，单元方程 4.3 UL的单元方程 Eulerian坐标与单元坐标 之间的映射为 位移可以由相同的形函数进行插值 4 UL格式的控制方程，弱形式，单元方程 4.3 UL的单元方程 形函数与时间无关，通过位移的导数得到速度和加速度，变分函数由同一形函数给出 变形率可以表示为形函数的形式为 4 UL格式的控制方程，弱形式，单元方程 4.3 UL的单元方程 通过一个B矩阵，将变形率表示为节点速度的形式 变形率可以表示为形函数的形式为 形函数的空间导数由链规则得到 4 UL格式的控制方程，弱形式，单元方程 4.3