弹性直梁问题变分原理及有限元素法.docVIP

弹性直梁问题的变分原理及有限元素法 讨论的问题：一变剖面的梁，一端固支，另一端简支。承受轴向拉力N，分布横向载荷以及端点弯矩的作用。 控制微分方程及边界条件（以梁的挠度w表示） 称谓：把满足方程及全部边界条件的挠度叫真实挠度，精确解；把满足基本边界条件但不满足微分方程和自然边界条件的挠度叫（变形）可能挠度。 最小势能原理（变分原理） 把载荷看作是不变的已知函数，把挠度看作是可变的自变函数。 整个系统的势能包括三部分： 梁的应变能： (2) 轴向应变能： 横向载荷势能： 系统总势能： * 除w为可变外，其余变量假定为已知的不变量。 最小势能原理：在所有变形可能的挠度中，精确解使系统的总势能取最小值。 由于是w的二次函数，不用变分法而用较初等的方法也能作出数学证明。 证明过程： 设是精确解，它满足微分方程及所有边界条件。 设是某一变形可能的挠度，仅知道它满足： 令： 当时，才是w的变分。由上式关系知满足： 在处， 在处， 与相应的总势能： 其中： Note：满足基本边界条件，可得 (功的互等定理) 功的互等定理：第一组力在第二组位移上所做的功等于第二组力在第一组位移上所做 的功。（真实力状态关于虚拟位移作功；虚拟位移产生的虚拟力关于真实位移作功）此式表明： 由的算式，当时， 或 当（轴受压）未达到临界压力时，。 所以： 式中的等号只有在为刚体位移时才能成立（即弹性能仅为零的情况）。 以上即是证明的最小势能原理。 上述的证明可普遍适用于其他类型的边界条件；也适用于其他复杂弹性力学体系。 精确解既然是使总势能取最小值，那么必有： 即： （正好补足了可能挠度尚未满足的边界平衡条件） 所以，最小势能原理与平衡条件（边界及内部）完全等价。 说明：尽管变分法与原问题的微分方程系统等价，但具体求解时，变分法涉及的导数阶次要低（这是能量法的优点，求得的解可能满足微分方程的连续性要求，也可能不满足）。 能量逼近解（不是真解）不能满足力的自然边界条件（选取可能位移时一般很难取得满足力的边界条件，故一般情况有限元不能获得满足力边界条件的解）。从这点上看，只能是近似解。当然，如果在满足力的边界条件的那些函数集合中选，则解的精度要大大提高（即更为逼近）。 强调三点： 上述证明的是真实挠度使 系统 势能取最小值。 又通过变分法证明了系统能量的极值曲线满足梁的微分控制方程（平衡方程）。需要深入认识的是：系统能量泛函中要求的挠度为可能挠度（即满足连续性与位移边界条件的曲线，但不满足微分方程）；变分的结果恰使极值曲线应满足微分方程及自然边界条件，补足了真实解的所需条件（理论上的等价性）。 微分方程解与能量泛函解对函数连续性的要求是不同的，故能量方法对函数的连续性要求较放松。以后在广义变分原理中会看到对解函数的连续性要求更低。 用Ritz法求解梁的弯曲问题 考虑：端固支，端简支，变剖面梁在轴向拉力和横向载荷联合作用下的平衡问题，求挠度函数（当EJ为变数时很难求解精确解）。 求解：设挠度表达式为： 式中，是变形可能的某一特解，即和是x的连续函数，且满足下列非齐次的位移边界条件： 在处， 在处， 为n个适当选定的变形可能的齐次解，即和都是x的连续函数，并且满足齐次的位移边界条件，即： 在处， 在处， 是n个待定常数。 改成用矩阵表示（便于计算机运算及与有限元方法对比） 记： 未定系数矢量由最小势能原理确定。 代入能量泛函： 代入w的计算举例： ） ） Note：① 代入积分式后成为常数，对变分无意义，故可在变分意义下忽略掉。 ② ） ③ 同理得其他项的计算结果。最终有：  计算 注意求导运算的规则：在行乘列的数中，对行变量求导，列向量不变；对列向量求导得行向量的转置。同时注意到上式中的矩阵对称性，即： 于是得： 上面的计算推导显粗，细做： 令： Ritz法中的关键取决于可能挠度的选取是否恰当。 若问题中的位移边界条件是齐次的，则； 若问题中的位移边界条件是非齐次的，不可少。 若级数的前n项已颇接近精确解，级数的后几项只起“修正”作用，那么少取几项也能解决问题。反之，级数的前n项与精确解相差颇远，则加重了后面各项的“修正”负担，那么级数项只好取得多一些。 iii).有限元素法求解梁的弯曲问题 基本步骤： ① 先将梁分割成若干个（如n个）有限单元。 ② 构造单元的无量纲局部坐标系（用于几何构造） 取第e个元素分析（如右图） 取无量纲局部坐标系: 线性构造几何坐标： (后者为规范式) 局部坐标实际上为坐标基函数（也称形状函数，或参数坐标）。从坐标的观点来理解，参数坐标并不是独立的，也不具有描述一个任一