弹性力学弹性力学05变分法.pptVIP

例： 图示薄板，宽为 a，高度为 b，左边和下边受连杆支承，右边和上边分别受有均布压力 q1和 q2 作用，不计体力。试求薄板的位移。 解： （1）假设位移函数 （a） 满足边界条件： 试在式（a）中只取两个系数：A1、B1 ，即 （b） （2）计算形变势能 U 将式（b）代入（5-16），有 （平面应力情形下形变势能公式） 积分得： （c） （c） （3）代入Ritz 法方程求解 ∵体力 ∴有 在右边界： 在上边界： 于是有： 将式（c）代入，得 （11-15） 联立求解，得： （f） 代入位移表达式（b），得： （g） 讨论： （1） 如果在位移式（a）中再多取一些系数如：A2、B2等，但是经计算，这些系数全为零。 （2） 位移解（g）满足几何方程、平衡方程和边界条件。 表明：位移解（g）为问题的精确解。 Ritz 法解题步骤： （1）假设位移函数，使其满足边界条件； （2） 计算形变势能 U ； （3）代入Ritz 法方程求解待定系数； （4）回代求解位移、应力等。 * 1. 弹性力学问题的微分提法及其解法： （1）平衡微分方程 （2）几何方程 （3）物理方程 （4）边界条件 应力边界条件； 位移边界条件。 定解问题 求解方法： （1）按位移求解 基本方程： （a）以位移为基本未知量的平衡微分方程; （2）按应力求解 基本方程： （a）平衡微分方程； （b）边界条件。 （b） 相容方程； （c） 边界条件。 （a） 归结为求解联立的微分方程组； 求解特点： （b） 难以求得解析解。 从研究微小单元体入手，考察其平衡、变形、材料性质，建立基本方程： §5-4 弹性体的形变势能和外力势能 2. 弹性力学问题的变分提法及其解法： 基本思想： 在所有可能的解中，求出最接近于精确解的解； 将定解问题转变为求解线性方程组。 弹性力学中的变分原理 —— 能量原理 直接处理整个弹性系统，考虑系统的能量关系，建立一些泛函的变分方程，将弹性力学问题归结为在给定约束条件下求泛函极（驻）值的变分问题。 （变分解法也称能量法） （a）以位移为基本未知量， 得到最小势（位）能原理等。 （b）以应力为基本未知量， 得到最小余能原理等。 （c）同时以位移、应力、应变为未知量， 得到 广义（约束）变分原理。 —— 位移法 —— 力 法 —— 混合法 —— 有限单元法、边界元法、离散元法 等数值解法的理论基础。 求解方法： 里兹（Ritz）法， 伽辽金（Galerkin ）法， 加权残值（ 余量）法等。 3. 弹性力学问题的数值解法： （a）直接求解联立的微分方程组（弹性力学的基本方程） —— 有限差分法； 基本思想： 将导数运算近似地用差分运算代替； 将定解问题转变为求解线性方程组。 典型软件：FLAC 实质： 将变量离散。 （b）对变分方程进行数值求解 —— 有限单元法、边界单元法、离散单元法 等 典型软件： ANSYS，MARC，ADINA，SAP，NASTRAN，ABAQUS 等； —— 基于有限元法的分析软件； UDEC —— 基于离散元法的分析软件； 基本思想： 将求解区域离散， 离散成有限个小区域（单元）， 在小区域（单元）上假设可能解， 最后由能量原理 （变分原理）确定其最优解。 —— 将问题转变为求解大型的线性方程组。 1. 形变势能的一般表达式 P x l0 ?l 单向拉伸： P ?l O P ?l 外力所做的功： 由于在静载（缓慢加载）条件下，其它能量损失很小，外力功全部转化杆件的形变势能（变形能）U： 杆件的体积 令： —— 单位体积的变形能， 称为比能。 三向应力状态： 一点的应力状态： ? x ?y ?z 三向应力状态： 一点的应力状态： ? x ?y ?z 由能量守恒原理，形变势能的值与弹性体受力的次序无关，只取决于最终的状态。 假定所有应力分量与应变分量全部按同样的比例增加，此时，单元体的形变比能： （a） 对于平面问题， 在平面应力问题中， 在平面应变问题中， 因此， （a） （b） 整个弹性体的形变势能： （c） 2. 形变势能的应变分量表示 在线弹性的情况下，由物理方程（2-16） ： 代入式（b），整理得： （d） （e） 将式（e）分别对3个应变分量求导，并将其结果与物理方程(d)比较，得： 从而， （5-15） 表明： 弹性体的比能对于任一应变分量的改变率，就等于相应的应力分量。 3. 形变势能的位移分量表示 只需将几何方程代入式(e)，得： 在上式中，只要将弹模、泊松比代换，即可得到平面应变中的相应公式。 （f） 由式(e)和(f)可知，形变势能是应变分量或位移分量的二次泛函。因此，叠加原理不再适用。