荷载移置等技术总结.ppt

人有了知识，就会具备各种分析能力， 明辨是非的能力。 所以我们要勤恳读书，广泛阅读， 古人说“书中自有黄金屋。 ”通过阅读科技书籍，我们能丰富知识， 培养逻辑思维能力； 通过阅读文学作品，我们能提高文学鉴赏水平， 培养文学情趣； 通过阅读报刊，我们能增长见识，扩大自己的知识面。 有许多书籍还能培养我们的道德情操， 给我们巨大的精神力量， 鼓舞我们前进。 * Company Logo Company Logo 弹性体的应变能和余能 单位体积的应变能(应变能密度)为 应变能是个正定函数，只有当弹性体内所有的点都没有应变时( )，应变能才为零。 单位体积的余能(余能密度)为 余能也是个正定函数。在线性弹性力学中 弹性体的应变能等于余能。 有限单元位移法的理论基础是变分原理。因此，研究有限变形问题的变分原理和弹塑性力学中的变分原理以及它们的广义变分原理以及有限元离散都是必要的。 线弹性力学的变分原理 弹性力学变分原理包括基于自然变分原理的最小位能原理和最小余能原理，以及基于约束变分原理的胡海昌-鹫津久广义变分原理和Hellinger-Reissner混合变分原理等。 最小位能原理的建立可以从虚位移原理出发。虚位移原理的表达式是 (3.2) 其中是 应力张量，如将弹性力学的物理方程式代入，则可得到 (3.1) 最小位能原理 Company Logo 因为 是对称张量，且 则有 (3.3) 由此可见式(3.1)中体积分的第一项就是单位体积应变能的变分。 在线弹性力学中，假定体积力 和边界上 面力 的大小和方向都是不变的，即可从 位势函数 和 导出，则有 (3.4) 其中 (3.6) 将式(3.3)和式(3.4)代入式(3.2)，就得到 (3.5) Company Logo 是系统的总位能，它是弹性体变形位能和外力位能之和。式(3.5)表明：在所有区域内连续可导的并在边界上满足给定位移条件式的可能位移中，真实位移使系统的总位能取驻值。还可以进一步证明在所有可能位移中，真实位移使系统总位能取最小值，因此式(3.5)所表述的称为最小位能原理。 最小余能原理 最小余能原理的推导步骤和最小位能原 理的推导类似，只是从虚应力原 理出发。 (3.6) 其中 它是弹性余能和外力余能的总和，即系统的总余能。式(3.6)表明，在所有在弹性体内满足平衡方程，在边界上满足力的边界条件的可能应力中，真实的应力使系统的总余能取驻值。还可以用与证明真实位移使系统总位能取最小值类同的步骤，证明在所有可能的应力中，真实应力使系统总余能取最小值，因此式(3.6)表述的是最小余能原理。 线性方程组解法 线性方程组的解法： 直接解法 包括高斯消去法、等带宽高斯消去法、三角分解法，以及适用于大型方程组求解的分块算法和波前法等。 迭代解法 雅可比迭代法、高斯-赛德尔迭代法、超松弛迭代法和共轭梯度法等。 利用矩阵的对称、稀疏、带状分布等特点提高方程求解效率是关键。 在方程组的阶数不是特别高时，通常采用直接解法。 当方程组的阶数过高时，为避免舍入误差和消元时有效数损失等对计算精度的影响，可以选择迭代方法。 高斯消去法的一般公式 第m次消元时，以第m-1次消元后的第m行元素为主元行，对第i行元素(im)的消元公式为， 若原系数矩阵是对称矩阵，则在消元过程中的待消矩阵仍然保持对称。 如果 雅可比迭代的原理 A为非奇异矩阵，主元非零 D为对角矩阵 通过迭代逐步逼近方程组的解， 波前法的特点 刚度矩阵K和载荷列阵P不按照自然编号进入内存，而按照参加运算的顺序排列；以集成完毕的自由度作为主元对其它列的元素进行消元修正。不形成整体刚度矩阵。 共轭梯度法 求解最优化问题的搜索算法。 迭代公式为， 搜索方向， 线性方程组， 单元载荷移置 有限元法的求解对象是单元的组合体，待求解的未知变量都定义在单元的结点上，因此作用在弹性体上的外力，需要移置到相应的单元结点上成为结点载荷。 载荷移置要满足静力等效原则。静力等效是指原载荷与结点载荷在任意虚位移上做的虚功(virtual work)相等。 结点的虚位移 单元的虚位移 令结点载荷为 处于平衡状态的弹性体，真实位移分量发生了位移边界条件所允许的微小改变，这个微小的改变就是虚位移(virtual displacement)。 1）集中力的移置 在单元内任意一点p作用集中力 由于虚位移是任意的， 例1、在均质，等厚的三角形单元ijm的一点p（0.25a, 0.5a）上作用有集中载荷Px, Py。 2）体力的移置 令单元所受的均匀分布体力为， 虚功