有限元基础教学课件PPT.ppt

连续钢桁梁桥力学分析，2010年 钢桥最危险工况的有限元分析 拼装最后4个节间钢桁梁 成桥状态 * 1. 结构的离散化 结构离散成有限个单元，相邻单元体仅在结点处连接以替代原结构。 (a) 弹性体原形 悬 臂 深 梁 求： ， ， ， ， 。 (b) 离散模型(三角形单元) 三、 有限单元法分析过程 * 位移法有限元以结构的位移作为基本解函数，其单元上的位移试函数取为低阶多项式函数，其项数等于单元节点上的自由度数。 为单元节点上的位移列阵，基本未知量。 为位移形状函数（设为已知形式） 对平面问题： （a） 2. 选择单元位移模式 (a)代入几何方程 （b） 3. 单元的力学特性分析 再代入物理方程 （c） 最小势能原理： （e） 建立单元节点力 与节点位移 的关系 （d） * 计算等效节点力：单元边界上的面力和单元体力按虚功等效原则转化为等效结点力。 建立结构整体刚度方程 将各单元刚度方程迭加 （f） 式中 为结构的整体结点位移列阵， 为结点载荷列阵， 为总体刚度矩阵。 再将 代入（c），得到各单元的应力值。 引入已知位移边界条件， 为非奇异方程，解此 得到节点位移向量 。 （d） * * * * 一、位移矢量 ， ， 取笛卡尔坐标轴 ，对空间上任一点处的任一 方向用矢量 表示，其单位方向矢量为 。 在外力下， 点，P点的位移矢 P点 §0.1 应变分析 * 引入算子 为梯度矢 （几何线性） 在单连通域中： 一一对应,多连通域中未必一一对应. 几何方程（6个）：位移应变之关系 P点处在 轴的三个微段的变化，得到变状态的6个分量 二、应变 * 取P点处一微平行六面体与xyz平行， 决定P点应力状态的6个分量记为 体力（外力）： 一、平衡方程：（由微六面体平衡所致） （要求 可导） §0.2 应力分析 一个面上应力可分解为一个正应力，二个剪应力分量 物体表面上面力矢 应力与外力在表面 上平衡，则 其中矩阵 此表面处外法线方向 二、应力边界条件 * 对均质各向同性线弹性体，有广义Hooke定律 称为弹性矩阵 （物理线性） 或 §0.3 本构方程(应力－应变关系) * 物理量： ——Lame常数 G——剪切模量 E——弹性模量 ——泊松比 或 （位移型， ） 或 连续性 物理性 方程（15个） 边界条件 未知函数（15个） 类型 (应力型， 自由边， ） 静力平衡 ， 弹性力学平衡问题 微分方程边值问题 （15个方程求解15个未知量，在 和 上） 解法：（1）位移法；（2）应力法；（3）混合法 §0.4 弹性力学平衡问题的微分方程提法 * （位移表示的应力边界条件） 由上解出 显式： 式中 物体表面 , 取未知函数 ，经代换 弹性力学位移法定解问题： * ， ， ， 2. 对平面应变问题有 1. 对平面应力问题有 应力应变关系： 应力应变关系： * 一、弹性体的形变势能 引例 1. 弹性力的功 弹性（簧）力： 方向： —弹簧刚度系数 或 力大小 ,（0-5-1） §0.5 弹性力学变分原理 （原长） （弹性力） * 变形过程是静平衡状态, 弹簧从A1运动到A2位置， 弹性力做功： 弹性（簧）力性质： 1 弹力做功： 仅与弹簧初、末变形量 、 有关， 与路径无关； 2 可正、可负； 3 特别 时， 注：外力做功 （原长） （弹性力） * 以变形量为 处为势能零点，则在 处，弹簧力做功为 特别：以弹簧的自然位置为零势点，则弹性能 （弹簧力做功的负值，因为外力与弹簧力反向，故弹性能相当于外力的功） 弹性能=外力做功 2. 弹性力场中的势能：弹簧力做功的负值 （外力功） (0-5-3) 弹簧力的势能 * 例：A为杆截面积，L—杆长。 材料试件单向拉伸曲线 产生均匀的轴向应力 ，轴向应变 . (0-5-4) ， 弹性杆在轴向力P作用下， 3. 弹性杆的形变势能（应变能） （原长） （弹性力） * 假定：弹性体在受力的过程中始终保持平衡 外力P做功: 以弹性应变能储存于杆中。 杆的应变能 等于外力功 （例：射箭） 定义：单位体积中具有的应变能称为(应变)比能，记为 对杆： 特别： 即 为线性关系(Hooke定律) (与路径无关) （如 轴沿杆