能量原理与与变分法.pptVIP

固体力学非线性数值方法 西安交通大学航天航空学院 宋亚勤 yqsong@mail.xjtu.edu.cn 2013年9-10月 第一章弹性力学简介 第二节：能量原理与变分法 谢谢！ 请提问、进行讨论！ 将上式整理、移项、合并，可得： 完全任意，且互相独立， 要使上式成立，则须有： —— Ritz 法方程 或称 Rayleigh- Ritz 法方程 说明： （1） 由 U 的位移表达式可知， U 是系数 的二次函数， 因而，上式为各系数的线性方程 组。 互不相关，因而，总可以求出全部的系数。 （2） 求出了系数 就可求得其它量，如位移、应力等 （3） 在假定位移函数时，须保证其满足全部位移边界条件。 2. 伽辽金（Galerkin）法 设取位移的表达式如下： 同时满足： （1）位移边界条件； （2）应力边界条件； 位移的变分： 将其代入伽辽金变分方程： 得到： 完全任意，且互相独立， 要使上式成立，则须有： 将物理方程和几何方程代入，有 —— 伽辽金（Galerkin）法方程 说明： （1） 与 Ritz 法类似，得 3m 阶的线性方程组，可求出3m个系数。 （2） 伽辽金（Galerkin）法与 Ritz 法的区别：在于设位移函数时，前者要求同时满足应力、位移边界条件，而后者只要求满足位移边界条件。 （1）位移变分方程 （2）虚位移方程 位移变分方程小结： ——也称 Lagrange 变分方程： （3）最小势能原理 说明： （1）只要求：虚位移满足位移边界条件； （2）对虚位移方程，也适用各种材料的物理方程。 如：塑性材料、非线性弹性材料等。 （4）伽辽金（Galerkin）变分方程 要求：可能（虚）位移满足： （1）位移边界条件； （2）应力边界条件。 x §4 应力变分方程 余能密度 P l0 ?l ? ? O （1）单向应力状态 设： —— 一般的应力应变关系 形变势能： ? d? ? ?0 ?0 —— 单位体积的形变势能 余能密度： —— 单位体积的形变余能 对线弹性体，显然有： ——应变能密度等于余能密度 表明：余能密度在数值上等于图中矩形面积减去 U1 后余下的面积。 一般情形： —— 单位体积的形变势能 —— 单位体积的形变余能 （2）三向应力状态 对线弹性体，有： 弹性体余能： 对线弹性体： 物体余能常用应力表示： （3）余能的变分 对照余能密度的表达式，有： 若将上式中应变分量利用几何方程表示成位移形式，有： 代入余能的变分表达式，有： 应力变分方程 设有任一弹性体，在外力的作用下处于平衡状态。其应力和位移分别为： —— 实际的应力和位移 建立：物体余能的变分与应力变分之间的关系。 （1）应力的变分 假设：作用于物体的体力不变，而应力分量发生如下变分： —— 常称为虚应力 变化后应力状态： （2）应力变分方程（假定可能应力是问题的解） 都满足平衡方程 并作用于同样的体力， 将其分别代入平衡微分方程，并进行比较，应有： 此应力状态满足平衡方程以及应力边界条件（基本边界条件）. （a） 张量表示 在位移给定的边界上， 由于应力的变分（增量）将引起一个允许表面力： 由边界上应力与边界面法向余弦关系，在位移给定边界上，应有： （b） 张量表示 在应力边界上，满足无外力的边界条件： （c） 张量表示 由余能的变分： 利用奥-高公式，将上式每一项作变换，如： 将其代入余能的变分，并整理有： 0 0 0 得到： 上式表明： 由于应力的变分，余能的变分等于允许表面力的变分在实际位移上所做的功（虚功）。 —— 应力变分方程， 也称Castigliano变分方程。 * 1、弹性体形变势能 2、泛函与变分 —— 最小势能原理、里兹（Ritz）法、伽辽金（Galerkin）法 3、位移变分方程 4、应力变分方程 —— 最小余能原理、卡氏（Castigliano）定理 5、自然变分原理和广义变分原理 6、弹性力学修正变分原理 1. 弹性力学问题的微分提法及其解法： （1）平衡微分方程 （2）几何方程 （3）物理方程 （4）边界条件 应力边界条件； 位移边界条件； 定解问题 求解方法： （1）按位移求解 基本方程： （a）以位移为基本未知量的平衡微分方程; （2）按应力求解 基本方程： （a）平衡微分方程； （b）边界条件。 （b） 相容方程； （c） 边界条件。 （a） 归结为求解联立的微分方程组； 求解特点： （b） 难以求得解析解。 从研究微小单元体入手，考察其平衡、变形、材料性质，建立基本方程： （3）混合解法 2. 弹性力学问题的变分提法及其解法： 基本思想： 在所有可能的解中，求出最接近于精确解的解； 将定解问题转变为求解线性方程组。 弹性力学中的变分原理 —— 能量原理