有限元分析及应用 第二讲.ppt

变分法 1、泛函与变分 函数 y=f(x) 求y 的极值，即求微分，由dy=0 可得。 泛函J=J [y(x)] 函数y(x)为自变量，J为函数y的函数，称J为y的泛函，求泛函的极值，即求变分， 由 可得。 例：平面上AB两点，连接AB的曲线很多，要求一条曲线使重物靠自重由A沿此曲线滑到B所需的时间最短，即求最速下降曲线。 显然，AB间直线路径最短，但重物运动的速度增长并不是最大，即下滑的时间并非最短。 设AB间有n条曲线 ，每条曲线对应一个时间 ，即T是y(x)函数，即泛函，求变分的极值则可得最速下降曲线 有关泛函的具体构造可参考相关教材 x y v p B A 变分法 2、变分原理 即微分方程边值问题的解等价于相应泛函极值问题的解。 变分法就是利用变分原理求解边值（即微分方程+边界条件）的一种方法。 微分方程 边值问题 泛函极值 问题 平面应力问题 变分原理 值得注意：构造微分问题的相应泛函并非容易。 里兹法 实例：求满足以下微分方程及边界条件的函数 构造相应的泛函： 假设满足边界条件的近似函数： 泛函极值转化为多元函数的极值： N个代数方程组： 里兹法 求解过程： 1、据微分方程及边界条件构造相应的泛函； 2、假设满足边界条件的近似函数，使泛函转化为以待定系数为变量的多元函数； 3、利用多元函数极值条件，求解线性方程组得到待定系数； 4、得到泛函的极值函数。 虚功原理、最小势能原理及变分法之关系 微分方程 边值问题 泛函极值 问题 变分原理 虚功方程 最小势能极值 能量守恒 纯数学描述 物理描述 回顾  1 场问题的一般描述 2 有限元法是一种数值方法 3 有限元法的基本思想（二次近似） 4 有限元法的基本步骤（5步） 5 有限元法的基本定义（节点、单元、节点力、节点载荷） 6 位移函数的构造方法及基本条件 7 位移函数的收敛性条件（协调元、非协调元） 8 材料力学与弹性力学区别 9 弹性力学的几个基本概念（位移、应力、应变等） 10 弹性力学的基本方程（平衡方程、几何方程、物理方程） 11虚功原理、最小势能原理及变分法（里兹法） 12 两种平面问题及轴对称问题 物理方程----应力应变关系 将应变分量表为应力分量的函数，可称为物理方程的第一种形式。若将式(2-10)改写成应力分量表为应变分量的函数的形式，并将式(2-9)代入，可得物理方程的第二种形式： 式(2-11)可用矩阵的形式表示如下： 式(2-12)可简写为： [D]称为弹性矩阵，它完全决定于弹性常数E和 2-4 两种平面问题 弹性力学可分为空间问题和平面问题，严格地说，任何一个弹性体都是空间物体，一般的外力都是空间力系，因而任何实际问题都是空间问题，都必须考虑所有的位移分量、应变分量和应力分量。但是，如果所考虑的弹性体具有特殊的形状，并且承受的是特殊外力，就有可能把空间问题简化为近似的平面问题，只考虑部分的位移分量、应变分量和应力分量即可。 平面应力问题 平面应变问题 特点;独立分量都只有三个。 平面应力问题 厚度为t的很薄的均匀木板。只在边缘上受到平行于板面且不沿厚度变化的面力，同时，体力也平行于板面且不沿厚度变化。 以薄板的中面为xy面，以垂直于中面的任一直线为Z轴。由于薄板两表面上没有垂直和平行于板面的外力，所以板面上各点均有： 另外由于平板很薄，外力又不沿厚度变化，可认为在整个薄板内各点均有： 于是，在六个应力分量中，只需要研究剩下的平行于XOY平面的三个应力分量，即 ，所以称为平面应力问题。 x y 0 t/2 t/2 z y 图 2-10 q q 平面应力问题 应力矩阵(2-2) 可以简化为： 平面应力问题 物理方程(2-10)中后两式可见，这时的剪应变： 由物理方程(1-10)中的第三式可见： 一般 ， 并不一定等于零，但可由 及 求得，在分析问题时不必考虑。于是只需要考虑 三个独立应变分量即可，于是应变矩阵(2-3-2)简化为： 平面应力问题 物理方程(2-10)简化为： 转化成应力分量用应变分量表示的形式： 平面应力问题 将(2-21)式用矩阵方程表示： 它仍然可以简写为： 弹性矩阵[D]则简化为： 平面应力问题 只有 三个应变分量需要考虑，所以几何方程(2-3) 简