第十四章弹性力学的变分解法.ppt

第十四章 弹性力学的变分解法 　14.1 弹性体的虚功原理 　14.2 贝蒂互换定理 　14.3 位移变分方程 最小势能原理 　14.4 最小势能原理推导以位移表示的平衡微分方程 及边界条件的实例 　14.5 基于最小势能原理的近似计算方法 　14.6 应力变分方程 最小余能原理 　14.7 基于最小余能原理的近似计算方法 　14.8 最小余能原理在平面问题和扭转问题中的应用 ＊14.9 弹性力学的广义变分原理 ＊14.10 哈密顿变分原理 ＊14.11 作为古典变分法革新和发展的有限单元法 第十四章 弹性力学的变分解法 14.1　弹性体的虚功原理 14.2　贝蒂互换定理 14.3　位移变分方程 最小势能原理 14.5　基于最小势能原理的近似计算方法 14.6　应力变分方程 最小余能原理 14.7　基于最小余能原理的近似计算方法 14.8　最小余能原理在平面问题和扭转问题中 的应用 ＊14.9　弹性力学的广义变分原理 ＊14.10　哈密顿变分原理 ＊14.11　作为古典变分法革新和发展的有限单元法 10.应变分量的变换公式： 在求近似解时，可取应力函数 （14-24） 求出　　后，代入（14-24）即得问题的近似解 将（14-24）代入（14-22），余能　　成为　　　　　　　　　　的二次函数，取极值的条件为 例1 边界条件为 设 （d） 取 为具有因子 因此应力函数可取为 （e） 矩形薄板，两端受按抛物线分布的拉力作用，求应力 解： 在（e）中仅取一项，即令 代入方程（14-22）得 若 得 相应的应力分量为 2、最小余能原理在扭转问题上的应用 在§14-4中已经写出柱形杆扭转时的应变能（因考虑的是线弹性体，故即为余能）为 按应力解法，杆的横截面上的切应力可表示为 外力功为 于是 可以证明，这里，变分方程 等价于 在求近似解时，取应力函数为下列形式： 将(14-27)代入(14-26)，使总余能能成为 的二次函数，其取极值的条件为 求出 后代入(14-27)，得到要求的应力函数的近似解。由此可求出应力分量。 例２ 解：　 此函数在边界上为零。由薄膜比拟可知，m和n必为偶数，先取一项 ，即 ，于是 矩形截面杆，边长分别为2a和2b，求最大切应力 取应力函数为 代入（14-26）中积分得到 由　　　　　解得 于是 代入到（9-12） 最大切应力发生在长边中点处，注意到　　　　　，其值为 最小势能原理，以位移作为独立的自变函数 最小余能原理，以应力作为独立的自变函数 一类变量的变分原理 以位移、应变和应力等三类变量作为独立的自变函数 以位移和应力等两类变量作为独立的自变函数 三类变量的广义变分原理 两类变量的广义变分原理 1、胡海昌—鹫津久一郎变分原理： （1)由总势能出发，用拉格朗日乘子法导出广义势能原理 体内满足几何方程——6个条件 现在要引入9个拉格朗日乘子 6个 3个 满足几何边界条件——3个条件 位移变分法中最小势能原理 为有条件的驻值问题 即位移 以及相应 ，为几何上可能： 使条件驻值泛函变成无条件驻值泛函 令广义势能为 变成 和 18个变量都可以任意变分的无条件驻值的泛函 从 ，并利用高斯积分定理，有 由于 ， ， ， 彼此独立且完全任意，故上式成立的条件为 在V内分别有 在 上，有 在 上，有 从(b)、(c)和(e)，若取 则分别表示物理方程、平衡微分方程和应力边界条件 于是 于是 如果取 重复上面的步骤，则可得到弹性力学的全部基本方程和全部边界条件 （2)由总余能出发，用拉格朗日乘子法导出广义余能原理 体内满足平衡微分方程——3个条件 满足应力边界条件——3个条件 位移变分法中最小余能原理 为有条件的驻值问题 这使条件驻值泛函变成无条件驻值泛函 现在要引入6个拉格朗日乘子 3个 3个 即应力分量 总余能 于是，广义余能为 由 ，并利用高斯积分定理，有 在V内分别有 在 上，有 在 上，有 显然，(j), (l),和(m)分别为物理方程、平衡微分方程和应力边界条件 由于 ， ， ， 彼此独立且完全任意，故上式成立的条件为 由式(k)和(o)可