1、第二章 薄板小挠度弯曲的变分方程及近似解法,薄板小挠度弯曲的变分方程 Ritz法 Ritz法的应用举例 Galerkin法 Galerkin法的应用举例,第二章 薄板小挠度弯曲的变分方程及近似解法,建立薄板小挠度弯曲问题的变分方程，用变分法推导弹性力学问题的基本方程和边界条件，并在此基础上发展一系列的近似解法，是解决弹性力学问题的一个重要途径。 薄板小挠度弯曲问题的变形能与余变形能 弹性体每单位体积变形能增量为：对于薄板小挠度弯曲问题，上式变为：沿厚度积分，用曲率、扭率和弯矩、扭矩表示板单位面积上的应变能密度增量把板的应变能密度看成是 、 和 的函数，应满足,2.1 薄板小挠度弯曲的变分方程,
2、由弯矩、扭矩以及曲率、扭率的表达式：得到单位面积上的应变能密度为：把内力分量 、 和 看成是自变量，板的余变形能密度满足：板的余变形能密度增量为：,2.1 薄板小挠度弯曲的变分方程,板的余应变能密度应满足：所以有：对于线弹性问题，变形能密度W与余变形能密度We相等。 在线弹性情况下，余变形能密度可以表示为：, 2.1 薄板小挠度弯曲的变分方程,几何可能解 ：在面内关于x，y连续，且二阶可导，满足所有的位移边界条件，可以由几何方程得到对应的可能变形。 静力可能解：在面内满足与外力q平衡，满足所有的力边界条件，由物理方程的到的变形不一定是几何可能的。 真实解: 满足所有方程和边界条件。 运用最小位
3、能原理推导薄板小挠度弯曲问题的微分方程和边界条件 最小位能原理：在一切几何可能解中，真实解使位能取得最小值。 对应几何可能解的位能表达式为：U是每单位体积的变形能，Xi是体积力， 为给定的边界力, 2.1 薄板小挠度弯曲的变分方程,对于薄板，板上作用与z方向相同的分布载荷q，在给定位移边界上，有 或 。在力的边界上，给定弯矩 ， 给定等效剪力 ，在k个未给定位移的角点上，有给定的集中力 。 用几何可能解表示的薄板的位能为：, 2.1 薄板小挠度弯曲的变分方程,真实解使位能的变分为零，对上式变分得薄板小挠度弯曲的变分方程。其中单位面积变形能（变形能密度）的变分为：, 2.1 薄板小挠度弯曲的变分
4、方程,由格林公式得：带入位能的变分得：, 2.1 薄板小挠度弯曲的变分方程,在位移边界上，几何可能解满足： ，上一个式子可以进一步简化为：真实解应使： ，由于变分 的任意性，真实解w应满足：变分方程与薄板弯曲微分方程及相应边界条件是等价的。, 2.1 薄板小挠度弯曲的变分方程,运用最小余能原理 ，我们同样可以得到薄板弯曲基本微分方程：最小位能原理和最小余能原理得到的变分方程，与薄板小挠度弯曲的微分方程及边界条件是等价的，对于微分方程，我们可以采用双三角级数或单三角级数法来求解，对于变分方程，可以采用Ritz法、Galerkin法等近似解法求解。, 2.1 薄板小挠度弯曲的变分方程,最小位能原理
5、指出：在一切几何可能解中，真实解使得薄板的位能取得最小值。假设真实解的表达式为：其中 为满足所有位移边界条件的设定函数，而ai是待定常数，它应该满足：因为 是任意的，故满足上式的条件是：如果 是完备函数族，且自由度 ，则所得到的解是精确解。一般情况下，若给出合理的设定函数，则只需要几项上式就可以得到足够精度的解。,2.2 Ritz法,方程变为代数方程把求w的问题转化为求ai的问题，,位能表达式写成关于挠度w的形式：如果每个边都满足至少下面边界条件之一：则有：得,2.2 Ritz法,例1：四边固支矩形薄板，利用Ritz法求解，承受均布载荷q，挠度假设为：满足边界条件将挠度表达式带入到位能表达式中
6、，得到,2.3 Ritz法的应用,2.3 Ritz法的应用,aij应满足：即：这是关于aij代数方程组，当m=1, n=1, 时有,2.3 Ritz法的应用,对于方板a=b：最大挠度发生在板的中点 当m=3, n=3时有只要位移函数取得合适，Ritz法具有足够的精度。增加m，n精度可以得到进一步改善。,2.3 Ritz法的应用,例2：两对边简支，一边固定，一边自由，承受均布载荷q，挠度假设为：满足所有位移边界条件 代入到位能的表达式,2.3 Ritz法的应用,真实解使得 ，若m=1，得到关于a1的代数方程 ：求解得： 所以挠度为：,2.3 Ritz法的应用,在一切几何可能解中，真实解应是使位能
7、的变分为零的解。薄板小挠度弯曲的变分方程为： 如果几何可能解在力的边界上满足所有力的边界条件，则位能的变分可进一步表达为 ：,2.4 Galerkin法,假设位移边界条件和力的边界条件均为齐次，Galerkin法假设板的挠度为：满足所有的边界条件， ai是待定的系数，满足变分方程，即：因为 是任意的所以有：这是一个关于ai的代数方程组，求解可得ai。,2.4 Galerkin法,物理含义： 在域内，近似的挠度不可能使得平衡方程 处处精确满足，但应该使残差值在各个位移位移函数所做的功在整个域上的积分为零。,2.4 Galerkin法,例1：仍以四边固支矩形薄板为例，利用Galerkin法求解，承受均布载荷q，挠度假设为：满足全部边界条件,2.4 Galerkin法的应用,若只取：m=1, n=1代入得到：,2.4 Galerkin法的应用,作业14-5,
