有限元法(杆系)2014版.pptVIP

第2章 杆系结构有限元法 1. 有限元法应用的简例（直杆） 2. 平面桁架系统 材料力学解答 有限元法解答 离散化 将直杆划分成n个有限段 两段轴线之间的连接点称为节点 每个有限段称为一个单元 第i个单元的长度为Li，包含第i和第i+1两个节点 用单元节点位移表示单元内部位移 把外载荷集中到节点上 建立结点的力平衡方程 * 1. 有限元法的简例（直杆） 例：自重作用下的等截面直杆 分析杆中位移、应变和应力分布 ①材料力学解答 ②有限元法解答 杆的长度为L，截面积为A， 弹性模量为E，单位长度的重量为q N(x) x u x L-x q(L-x) 3单元情况 x 节点2： 节点3： 节点4： 平衡方程 约束条件 a a a a a a ② ① ③ 1 2 3 4 有限元方法分析的一般步骤 d 节点力（内、外力）平衡条件； a 离散，建立有限元模型：节点和单元 b 确定单元位移分布（由结点位移进行插值） c 单元特性分析； e 引入位移边界条件后求解； 例： 2 1 3 （1） （2） F 已知 直杆长度为l，两杆夹角为 45度。两杆的横截面面积 为A，材料弹性模量为E。 求：图示桁架各杆端位移。 2 平面桁架系统 步骤1：离散和选择单元类型 解： 将系统离散为2个单元，3个节点（其中1个共用） 2 1 3 （1） （2） 2 u为单元中x方向位移（v是y方向位移，不产生变形），杆中任一点位移，可假设为(两节点，两参数）： (2-1) 待定系数 位置坐标 步骤2：选择单元位移函数 i j l 代入节点的位移, 步骤2：选择单元位移函数 所以： (2-2) 受轴力作用杆横截面上的应力为： 步骤3：单元特性分析（确定应变、应力、轴力） 步骤3：单元特性分析 应力与应变的关系是： 目标：确定节点位移与节点力的关系 杆端力为： 矩阵形式表示 (2-3) 步骤3：单元特性分析 局部坐标下的节点力列向量 局部坐标下的节点位移列向量 单元刚度矩阵 步骤3：单元特性分析 式（2-3）还可写成 （2-4） 式（2-3）或（2-4）称为单元刚度方程 须首先进行坐标转换 在杆系结构有限元法中，每个单元都有自己的局部坐标系， 但对整个结构而言需建立统一的整体坐标系。 把任意的一个单元取出来，放在整体坐标系下，考察一下该单元在两种坐标系下的物理量的转换（变换）关系。 符号约定：局部坐标系下的物理量用加上画线来标记。 首先定义一下整体坐标系x-y与局部坐标系 的夹角符号： x-y坐标系沿逆时针转动到与 坐标系重合，则x-y坐标系转过的角度 为正。 步骤4：节点力平衡方程 步骤4：构造总体刚度矩阵 杆端力的坐标变换 矩阵形式表示 或 （2-5） 同理 （2-6） A.局部坐标系下的物理量： B.结构坐标系下的物理量： 步骤4：构造总体刚度矩阵 单元刚度矩阵的坐标变换: 在局部坐标系下，单元的刚度方程为： 由式（2-5）可得： 由式（2-6）可得： 带入可得： 使： （2-7） （2-8） 结构坐标系下的单元刚度矩阵 结构坐标系下的单元刚度方程 步骤4：构造总体刚度矩阵 所以：对单元1，α＝0 对单元2，α＝45 构造结构刚度矩阵（节点力合成） 步骤4：构造总体刚度矩阵 对于节点1有： 节点1的节点力 单元(1)中节点1的节点力（杆端力） 同理有：对节点2： 对节点3： 任一单元(e)的单元刚度方程为： 或 * 第i个单元中的位移用所包含的结点位移来表示： （1- 1） 其中为第i结点的位移，为第i结点的坐标。 第i个单元的应变为，应力为，内力为： （1- 2） （1- 3） （1- 4） 把第i单元和第i+1单元重量的一半，集中到第i+1结点上。 集中单元重量 对于第i+1结点，由力的平衡方程可得： （1-5） 令，并将（1- 5）代入得： （1-6） 根据约束条件，。对于第n+1个结点， （1-7） 建立所有结点的力平衡方程，可以得到由n+1个方程构成的方程组， 可解出n+1个未知的接点位移。