弹性力学有限元第五章 变分法解平面问题.pptVIP

河海大学 机电工程学院 力学教研室 第五章 变分法解平面问题 第五章变分法解平面问题 第五章 变分法解平面问题 §5-1 变分法简介 函数的变分 函数 y的微分： 假想函数 y(x)的形式发生改变成为新函数Y(x), 如果对应于x的一个定值，y具有微小的增量： 则称d y为函数 y(x)的变分,如图所示。 假定AB 表示某梁的一段挠曲线，而y是梁截面的真实位移，则CD可表示该梁发生虚位移以后的挠曲线，而虚位移d y就是真实位移y(x)的变分。 第五章 变分法解平面问题 §5-1 变分法简介 函数的变分 函数 y的导数 y’的变分： 函数 y的变分的导数： 导数的变分等于变分的导数 因此总有关系式 第五章 变分法解平面问题 §5-1 变分法简介 泛函及其变分 如果对于某一类函数 y(x)中的每一个函数y(x)，变量I 有一个值和它对应，则变量I称为依赖于函数y(x)的泛函，简单的说，泛函就是函数的函数。 例如：曲线AB的长度可写为： 长度L是依赖于函数 y(x)的泛函。 一般情况下，泛函具有如下的形式： 其中被积函数： 也是 y(x)的泛函。 第五章 变分法解平面问题 §5-1 变分法简介 泛函及其变分 泛函 f 的变分表示为： 高阶项 泛函 f 的变分 当函数y(x)具有变分d y，导数具有变分d y’。对泛函 进行泰勒级数展开，并求其增量： 第五章 变分法解平面问题 §5-1 变分法简介 泛函及其变分 考察泛函 增量的主要部分定义为泛函的变分，则 代入d f，则 显然，存在关系式： 只要积分的上下限不变，变分的运算和定积分运算可以交换次序 第五章 变分法解平面问题 §5-1 变分法简介 泛函的极值问题—变分问题 函数 y(x)在 x =x0取极值的条件是 对于泛函I [ y(x)]，如果在 y =y0(x)的邻近任意一根曲线上的值都不大于或都不小于I [ y0(x)]，也就是： 则称泛函I [ y(x)] 在曲线y =y0(x) 上达到极大值或极小值： 泛函达到极值的必要条件是： 曲线y =y0(x) 称为泛函I [ y(x)] 的极值曲线。 第五章 变分法解平面问题 §5-1 变分法简介 泛函的极值问题—变分问题 设图中 y = y(x)所示曲线通过A,B两点， y(x)具有边界条件： 由泛函 I 的极值条件求出函数y(x)，其中 [解]：首先导出d I的具体形式 考虑其第二部分， 对其进行分部积分 第五章 变分法解平面问题 §5-1 变分法简介 泛函的极值问题—变分问题 代入式 按照边界条件，x =a,x=b, y不变，因此: 则有： 得到： 第五章 变分法解平面问题 §5-1 变分法简介 泛函的极值问题—变分问题 根据d y的任意性,由d I =0 得出极值条件： 由此得出函数 y(x)的微分方程，该微分方程的解答给出函数 y(x) 的表达式。 第五章 变分法解平面问题 §5-1 变分法简介 泛函的极值问题—变分问题 例：求图中AB曲线为最短时的函数 y(x) 解： 由极值条件 ： 解此方程， 因此函数： 直线。 第五章 变分法解平面问题 §5-2 弹性体的形变势能 按照材料力学知识，弹性体只在某一方向上具有均匀的正应力sx和相应的正应变ex，其变形能密度(比能)为: 弹性体只在某两个正交方向受有均匀的切应力txy和相应的切应变gxy，则其比能为： 如果弹性体受到全部六个应力分量, sx, sy, sz, tyz, tzx, txy 弹性体的全部比能为： 平面问题考虑到tyz, tzx为0，平面应力问题有sz =0 ，平面应变问题有ez =0，则必能简化为： 第五章 变分法解平面问题 §5-2 弹性体的形变势能 取z方向为单位长度，考虑到应力应变是坐标的函数，则形变势能为比能在弹性体内的积分： 引入物理方程(平面应力)： 代入比能U1的表达式，得： 对ex, ey, gxy求导： 第五章 变分法解平面问题 §5-2 弹性体的形变势能 弹性体每单位体积中的形变势能对于任一形变分量的改变率等于相应的应力分量： 考虑几何方程，则比能可用位移分量来表示。 积分可得形变势能。 平面应变问题作弹性常数的替换。 第五章 变分法解平面问题 §5-3 位移变分方程 设有平面问题中的任一单位厚度的弹性体，在外力作用下平衡。 u,v为其实际位移分量，假设这些位移分量发生了位移变分(虚位移)d u, d v，成为： 考察其能量方面的变化。 假设弹性体在虚位移过程中没有热能或动能的改变，则形变势能的增加等于外力势能的减少，即外力所做的功，就是所谓的虚功。 外力包括体力和面力。 形变势能的变分写为： 位移变分方程(拉格朗日变分方程) 第五章 变分法解平面问题 §5-3 位移变分方程 由于虚位移微小，因此在虚位