有限元第一章.pptVIP

（3）物理方程： —— 弹性矩阵 若令： —— 拉梅（Lam’e）系数 其中： —— 柔度矩阵 （4）应力边界（ S? ）条件： —— 边界外法线方向余弦矩阵 （5）位移边界（ Su ）条件： （6）弹性体的应变能与余能 单位体积的应变能（应变能密度）： —— 正定函数 单位体积的余能（余能密度）： —— 正定函数 在线性弹性力学中，有 1.4.2 弹性力学基本方程的张量形式 1. 基本量的张量表示 应力： 应变： —— 二阶对称张量 位移： 体积力： 面力： —— 一阶张量 —— 二阶对称张量 2. 基本方程的张量表示 平衡微分方程： （在域V内） O x y z (x1) (x2) (x3) —— 表示对坐标 xj 的导数 上述方程展开，有 几何方程： （在域V内） 物理方程： （在域V内） —— 四阶张量 （共有34=81个分量，即有81个弹性常数） 由于 的对称性， 因而有 对于等温、绝热的变形过程，有 此时独立的弹性常 数为21个。 进一步对各向同性的线弹性材料，独立的弹性常数只有 2个。 如拉梅常数 G、? 或弹性模量E、泊松比? ， 此时弹性张量为： —— 线弹性材料的物理方程 其中： （2）选取适当的加权余量法，求解； （a）配点法 一项近似解： 取一个配点：x =1/2 ; ∴所求一项近似解为： 两项近似解： 取两个配点： 求解得： ∴所求两项近似解为： 近似解的项数得越多，则解的精度越高。 解得： 即取： 代入： （b）子域法 一项近似解： 取，子域 = 全域 取权函数： 代入： 积分得： 一项近似解为： 两项近似解： 取权函数： 代入： 求解得： 两项近似解为： （c）最小二乘法 将余量的二次方 R2 在域 ? 内的积分： 通过适当选取待定系数 ai ，使 I 达到极小。 即： （9） 将式（9）代入，有 （10） 由此得到 n 个方程，可求解 n 个待定系数 ai 。 将式（10）与式（4）比较： （4） 显然有最小二乘法的权函数： 一项近似解： 代入式（10）： （10） 积分后，解得： 一项近似解为： 两项近似解： 代入式（10）积分，有： 积分后，解得： 两项近似解为： （d）力矩法 一项近似解： 代入式（4）： （4） 求得： 一项近似解为： 与子域法一项近似解相同。 两项近似解： 权函数： 代入式（4）： 求得： 两项近似解为： （e）伽辽金（Galerkin） 法 一项近似解： 权函数： 代入： （4） 解得： 一项近似解为： 两项近似解： 权函数： 代入： （4） 积分后解得： 两项近似解为： 讨论： —— 边界条件 （1） （2） 精确解： 各种方法的近似解与精确解的比较结果： （1）两项近似解均得到了很好的精度，误差 3%; （2）Galerkin 法的结果精度较高; （3）可以推断，随着选取近似解项数的增加，精度会更高; 1.3 变分原理和里兹（Ritz）方法 1.3.1 变分原理的定义和意义 1. 变分原理与变分法 若一连续介质问题存在一标量泛函 ? ： （1.3.1） 则连续介质问题的解 u 一定使泛函 ? 对微小变化? u 取驻值，即，使泛函 ? 的“变分”等于零： （1.3.2） 称为变分原理。 由变分原理求解连续介质问题的方法称为变分法。 说明： （1）要求存在某一标量泛函? 连续介质力学问题； 热传导问题； 流场问题； 电磁场问题等。 （2）是等效积分形式的一种特殊情形。 对式（1.3.1）求变分，有 （2）是等效积分形式的一种特殊情形。 对式（1.3.1）求变分，有 （3）弹性力学中基本变分原理： 最小势（位）能原理 最小余能原理 平衡微分方程+力的边界条件 几何方程+位移边界条件 2. 变分法的求解过程 （1）选取未知函数 u 的近似解； （1.3.3） 注意：使 u 满足强制边界条件。 （2）将函数 u 的近似解代入泛函 ? ( u ) : ~ ~ （3）对泛函 ? ( ai ) 求变分，并令等于零； ~ 2. 变分法的求解过程 （1）选取未知函数 u 的近似解； （1.3.3） 注意：使 u 满足强制边界条件。 （2）将函数 u 的近似解代入泛函 ? ( u ) : ~ ~ （3）对泛函 ? ( ai ) 求变分，并令等于零； ~ （1.3.4） 由于 是任意的， 故上式成立时，必有： 将上式表示成矩阵形式，有： 其中： 得到与待定参数 a 的个数相等的方程组，由此可求得待定参数a 。 —— 里兹（Ritz）法 （1.3.5） 特殊情形： （1.3.6） 式（1.3.6）为一线性方程组。 式中，K 为一对称的常系数矩阵。 若泛