2.一维Lagrangian和Eulerian有限元讲解.ppt

例2.5 3节点二次位移单元 - 检查网格畸变 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和UL格式的控制方程。 4. TL和UL格式虚功原理及每一项物理意义。 5. 如何从弱形式到强形式？ 6. 从TL和UL格式的虚功原理建立有限元离散方程，方程中每一项的表达式？ 7. 显式积分求解的条件稳定性准则和收敛性控制准则。 TL格式比UL格式需要更多的存储空间，以存储形函数及其导数值；而UL则需要在每一个时间步重复搜索和计算形函数，也会影响计算效率。因此，在实际问题中应有所选择，例如对于大变形的瞬态问题，或者路径无关材料，可采用TL，而与变形历史有关的路径相关材料，如弹塑性和粘弹塑性材料，则可采用UL。 * * * 式中的第二个定义有时不明确，因为它出现了涉及一个独立变量x对于另一个独立变量X 的偏微分，这是没有物理意义的。因此，务必理解为无论何时x在文中出现，它意味着是一个函数 x=(X,t)。 * 体力b与加速度的量纲一致 Rou W(int)每单位体积的内能，其中W(int)是每单位质量的内能 * * [[f]]表示在f(x)中的跳跃，在材料界面处，应变、即位移的导数不连续，需要补充该条件。跳跃函数给出内部连续条件 * 较少的强制连续性要求－跳跃函数给出内部连续条件 * 问：为什么消除关于应力P 的导数？ 由于在这种弱形式中会出现应力的导数，位移和速度就不得不是C1函数；在高于一维的情况下C1函数是不容易构造的。而且，不得不随之构造试函数以便于满足面力边界条件，这也是困难的。通过分部积分消去应力的导数，在线性化方程中也导致了某些对称性，这将在第6章中见到。因此，分部积分是建立弱形式的关键步骤。 消除关于应力P 的导数，否则，力边界条件不得不施加在试函数上。另外，对应力函数降低了平滑性，可以C0函数。 建立弱形式中的关键步骤是分部积分，d(uv)/dx=udv/dx+vdu/dx * 如果我们允许较低平滑的变分项和试函数，在强形式中将附加一个方程，即内部连续条件(5.3.13c)；如果选取的变分项和试函数满足经典的平滑条件，在强形式中则没有内部连续条件。 * 二阶导数唯一定义的意义是一阶导数连续，二阶导数不为零（存在） * 建立弱形式的等价性 * * 直接写出虚能量方程，不用变分项乘以方程得到弱形式。 * U上杠是指定位移边界。 注意与书中的定义区别 * 单元形函数相对于母单元形函数，总体形函数是域内所有与该节点有关的形函数求和。 * * 初始条件回到材料坐标上。 (1-x)**pai—无理函数，(1-x)**n/(1-y)—有