有限元课件 第2讲 矩阵分析及弹性力学基础.ppt

对张量的理解 张量不随坐标系的改变而改变 例如位移矢量ui ：无论从哪一个坐标系观察，它反映的总是A点移动到B点的客观事实，不随观察者所在的坐标系而改变。 再如应力张量σij和应变张量εij ，尽管在不同的坐标系中具有不同的分量，但是它们所描述的却是某点的同一个应力和应变状态。 Voigt标记 含义：在有限元编程中，常常将对称的二阶张量写成列向量，将非常棘手的对称四阶张量（如弹性系数矩阵Dijkl）转换成二阶张量。这种转换过程称为voigt标记。 转换规则： 应力张量（动力学量）的转换 应变张量（运动学量）的转换 应力张量的Voigt标记 应变张量的Voigt标记 剪切应变需要乘以2，这是源于能量表达式的需要。 弹性系数矩阵的Voigt标记 平面问题及其基本方程 弹性体在满足一定条件时，其变形和应力的分布规律可以用在某一平面内的变形和应力的分布规律来代替，这类问题称为平面问题。平面问题分为平面应力问题和平面应变问题。 平面应力问题 很薄的等厚薄板，只在板边上受到平行于板面并且不沿厚度变化的面力，体力也平行于板面且不沿厚度变化。 平面应变问题 很长的柱形体，支承情况不沿长度变化，在柱面上受到平行于横截面而且不沿长度变化的面力，体力也如此分布。 三大类基本方程 在弹性力学中针对微小的单元体建立基本方程，把复杂形状弹性体的受力和变形分析问题归结为偏微分方程组的边值问题。弹性力学的基本方程包括 平衡方程：内力和外力的关系 几何方程：应变和位移的关系 物理方程（本构方程）：应力和应变的关系 平衡方程 ab=dxad=dy 习惯上 张量指标形式： 单位体积力 几何方程 张量指标形式： 物理方程 平面应力问题： 平面应变问题： 张量指标形式： 边界条件 o x y ds 外法线n的方向余弦  l=dy/ds m=dx/ds 位移BC 力 BC 张量指标形式： 三维问题基本方程及边界条件 可以将平面问题的基本方程推广到三维问题。基本变量如下： 位移： 应变： 应力： 平衡方程 几何方程 物理方程 边界条件 三维问题基本方程的张量指标形式 平衡方程： 几何方程： 物理方程： 边界条件： 作业：一维拉杆问题(忽略体力) 用弹性力学的基本方程和边界条件求解拉杆应力、应变及位移分布。 E=210000MPa，A=Pi*502 mm2，l =1000mm , P=600N 第2讲 矩阵算法及弹性力学基础 2.1 矩阵算法 线性方程组的表示 行向量和列向量 矩阵加、减、乘法运算 矩阵的转置、对称矩阵、单位矩阵 矩阵行列式 矩阵求逆 矩阵的微分和积分 正定矩阵（正定二次型） 线性方程组的表示 求解方法：高斯消元法、迭代法 行向量和列向量 矩阵加、减、乘法运算 对称方阵 矩阵转置、对称矩阵、单位矩阵 或 矩阵行列式 奇异矩阵(方阵) 如果方阵A的行列式 则其逆存在，记为 A的伴随矩阵 矩阵的逆 对于： 线性方程组的求解，变为求解系数矩阵的逆矩阵 矩阵的微分和积分 二次型：含有n个变量的二次齐次多项式 若取 正定二次型 则 利用矩阵及其运算，二次型可表示为 A: 对称矩阵 正定二次型：设 为实二次型，如果对于 任意的非零实向量X，都有 A: 正定矩阵 关于正定矩阵 正定矩阵是特殊的对称实矩阵 正定矩阵的对角元aii0 正定矩阵的行列式|A|0 A为正定矩阵的充要条件是A的所有顺序主子式皆大于0 二次型的微商 对向量x各元素的偏导数 2.2 弹性力学基础 关于弹性力学 五个基本假定 外力和内力 应力、应变、位移 指标记法和求和约定 张量及Voigt标记 平面问题基本方程及边界条件 三维问题基本方程及边界条件 关于弹性力学 弹性力学是研究弹性体在约束和外载荷作用下内力和变形分布规律的一门学科。 力学学科 研究对象 特征 中学力学 质点 无变形 理论力学 质点系及刚体 无变形 材料力学 简单变形体(构件) 小变形 结构力学 数量众多的简单变形体 小变形 弹性力学 任意变形体 小变形 弹塑性力学 任意变形体 任意变形 力学学科各分支的关系 五个基本假定 连续性：无空隙，能用连续函数描述 均匀性：各个位置物质特性相同 各向同性：同一位置的物质各个方向上具有相同特性 线弹性：变形和外力的关系是线性的，外力去除后，物体可恢复原状 小变形：变形远小于物体的几何尺寸，建立基本方程时可以忽略高阶小量。 外力和内力 体力—分布在物体体积内的力，例如重力和惯性力。 面力—分布在物体表面上的力，例如接触压力、流体压力。 分布力：连续分布在表面某一范围内 集中力：分布力的作用面积很小时的简化 内力—外力作用下，物体内部相连各部分之间产生的相互作用力