一维Lagrangian和Eulerian有限元.pptVIP

5 求解方法 5 求解方法 为了求解非线性问题，最简单的方法，即时间显式积分。最广泛应用的显式方法是中心差分方法，采用对角或集中质量矩阵。 速度 加速度 在时间间隔中点的导数值由在间隔端点处函数值的差得到，顾名思义为中心差分公式。 从t＝0出发，取时间步长Δt 5 求解方法 5 求解方法 对于位移的更新不需要代数方程的任何解答，因此，在某种意义上，显式积分比静态线性应力分析更加简单,不需要矩阵求逆解刚度方程。 如在流程图中看到，对于控制方程和时间积分公式，大多数的显式程序是直接向前赋值。程序从施加初始条件开始, 第一个时间步多少与其它时间步的不同在于它仅取半步，这使程序能正确地解释关于应力和速度的初始条件。 大部分程序运算时间是在计算单元节点力，尤其是内部节点力。节点力是逐个单元进行计算的。在开始计算前，从总体的列矩阵中集合出单元节点速度和位移。如流程图所示，内部节点力的计算包括应变方程和本构方程的应用。通过应力为内部节点力赋值。当完成了单元节点力的计算，根据它们的节点编号将其离散到总体列矩阵。 5 求解方法 稳定性准则 显式积分的缺陷在于时间步长必须低于一个临界值，否则由于数值不稳定将使解答隆起。对于采用对角质量的2节点单元的临界时间步长为 是单元的初始长度 是波速 收敛性控制（结果稳定性）准则 能量守恒(断裂力学中为能量平衡) 增加时间步长：放大质量，调整单元尺寸。 6 Eulerian格式的控制方程 6 Eulerian格式的控制方程，弱形式，有限元方程 在Eulerian格式中，节点在空间固定，相关变量为Eulerian空间坐标 x 和时间 t 的函数，应力度量为Cauchy（物理的）应力，变形度量为变形率，运动由速度描述。在Eulerian格式中，因为不能建立未变形、初始的构形，所以不能将运动表示为参考坐标的函数。 作业: 2.14 练习: 1, 2, 3, 6 问题：1. TL和UL的异同，基本变量的区别。 2. 为什么固体力学分析应用L格式，而不采用E格式？ 3. TL格式的控制方程。 4. TL格式的虚功原理及每一项的物理意义。 5. 如何从弱形式到强形式？ 6. 从TL格式的虚功原理建立有限元离散方程，方程中每一项的表达式？ 7. 显式积分求解的条件稳定性准则和收敛性控制准则。 TL格式比UL格式需要更多的存储空间，以存储形函数及其导数值；而UL则需要在每一个时间步重复搜索和计算形函数，也会影响计算效率。因此，在实际问题中应有所选择，例如对于大变形的瞬态问题，或者路径无关材料，可采用TL，而与变形历史有关的路径相关材料，如弹塑性和粘弹塑性材料，则可采用UL。 谢谢 Written by 庄茁 ? Copyright 2009 * * * 变形梯度的偏微分是没有物理意义的。 * 体力b与加速度的量纲一致 * [[f]]表示在f(x)中的跳跃。 * 较少的强制连续性要求－跳跃函数给出内部连续条件 * 问：为什么消除关于应力P 的导数？ 消除关于应力P 的导数，否则，力边界条件不得不施加在试函数上。建立弱形式中的关键步骤是分部积分，d(uv)/dx=udv/dx+vdu/dx * 直接写出虚能量方程，不用变分项乘以方程得到弱形式。 * U上杠是指定位移边界。 * 初始条件回到材料坐标上。 * 运动方程－半离散方程 3 有限元离散，单元和总体矩阵 从弱形式的一致性推导出的质量矩阵称为一致质量矩阵。在许多应用中，采用对角质量矩阵(集中质量矩阵)更有优势。质量矩阵对角化的过程是相当特殊的，这些过程没有理论。最常用的一种过程是对矩阵的行求和 式中用到了这样的事实，形函数对行求和必须等于1。这种对角化的过程使物体的总动量守恒，例如对于任意的节点速度，对角质量的系统动量应该等价于一致质量的系统动量： 对角质量矩阵也可以由下式赋值 运动方程在空间是离散的，在时间上是连续的，有时简称离散方程。在有限元离散中，质量矩阵常常为非对角阵(一致质量矩阵)，此时运动方程区别于牛顿第二定律，当MIJ≠0时，节点J处的力可以在节点I处产生加速度。而集中质量矩阵的运动方程等价于牛顿第二定律。 3 有限元离散，单元和总体矩阵 上式为在质点I上的静力。由牛顿第三定律，作用在节点上的力大小相等，而方向相反，因此内部节点力需要一个负号。 运动方程－半离散方程(矩阵形式) 单元和总体矩阵 在有限元程序中，通常以一个单元水平计算节点力和质量矩阵，将单元节点力结合入总体矩阵，称为离散或矢量组合。 组