杆梁问题的有限元法.pptVIP

有限元法基础Finite Element Method 第二章 杆件静力学问题有限元法 要 点 了解杆件系统有限元法的解是精确解； 杆件系统中的单元都在单元局部座标中讨 论，然后再投影到整体座标中去； 杆件系统中的单元是梁单元，其变形以弯曲 和拉伸变形为主，常忽略剪切变形的影响。 第二章 杆件静力学问题有限元法 铰连接 拉压力 与截面大小有关 简支、固支等连接 拉压、弯曲、扭转 与截面大小、形状有关 梁： 桁杆： 第一节 引 言 [K]称为对应于施加在系统上各节点力的刚度矩阵 引 言 引 言 一、单个弹簧的刚度矩阵 引 言 二、组合弹簧的刚度矩阵 引 言 所研究的对象是细长的杆件，即轴线方向的尺寸远比其他二个方向的尺寸大的多，如轴支柱、加强筋等； 杆件系统可分为平面杆件系统和空间杆件系统，若组成结构物的杆件都在同一个平面内，外力也在这一平面内，则称为平面杆件系统； 若组成结构物的杆件不在同一个平面内，则称为空间杆件系统； 杆件系统中所用的单元主要为梁单元，如果各杆件只受拉压作用则采用桁架单元。 第二节 杆件问题有限元法 2.1 桁架结构有限元法的研究对象 杆系结构单元 2.2 概述 杆系结构：梁、拱、框架、桁架等，它们常可离散成杆单元和梁单元。 拱 框架 桁架 结构离散 梁单元离散，各杆有局部坐标系，方向为杆的轴线方向。结构的刚度方程是在统一的坐标系（总体坐标系）中建立并求解，因此需将单元局部坐标各量转换到总体坐标系中。 取杆件与杆件交点、集中力作用点、杆件与支承的交点为节点。相邻两节点间的杆件段是单元。节点编号时力求单元两端点号差最小。 坐标系 有限元中的坐标系有整体坐标系和局部坐标系。对于一个结构，整体坐标系一般只有一个；而局部坐标系有很多个，一个单元就有一个局部坐标。对不同局部坐标系每一个单元的规定都是相同的，这样，同类型单元刚度矩阵相同。 杆系结构单元主要有铰接杆单元和梁单元两 种类型。它们都只有2个节点i、j。 约定：单元坐标系的原点置于节点i；节点i到j的杆轴（形心轴）方向为单元坐标系中x轴的正向。 y轴、z轴都与x轴垂直，并符合右手螺旋法则。 对于梁单元， y轴和z轴分别为横截面上的两个惯性主轴。 单元分析 下图示出了一维铰接杆单元，横截面积为A，长度为l，弹性模量为E，轴向分布载荷为px。单元有2个结点i，j，单元坐标为一维坐标轴x。 1、一维杆单元 （1）位移模式和形函数 ① 位移模式 因为只有2个结点，每个结点位移只有1个自由度，因此单元的位移模式可设为： （2-3） 式中a1、a2为待定常数，可由结点位移条件 x=xi 时， u=ui x=xj 时， u=uj xj-xi=l 确定。再将由此确定的a1、a2 其代入式（2-3），得 ② 形函数 将式（2-4）改写为下列形式 式中形函数[N]为 （2）应变矩阵 一维铰接杆单元仅有轴向应变 将式（2-5）、（2-6）代入上式，得 上式也可写为 式中[B]为应变矩阵 由应力应变关系 （3）应力矩阵 将式（2-7）代入上式，得 式中[S]为应力矩阵 (4) 单元刚度矩阵 单元刚度矩阵由虚功原理推出 对于等截面铰接杆单元(截面积为A ) ，dv=Adx，故有： (5) 等效节点力 单元上作用分布力px，则等效节点力计算公式仍为以下形式 当分布力集度px为常数时，有 （2-13） 将式（2-8）代入上式，得 例2.1 一维拉杆 图示阶梯形直杆，各段长度均为l，横截面积分别为3A，2A，A，材料重度为γ，弹性模量E。求结点位移和各段杆中内力。 离散化：将单元划分为3个单元，4个结点。 单元刚度矩阵： 等效结点荷载：按静力等效原则，有: 对号入座，组成总刚，形成整体结构平衡方程: 设结点1的约束反力为F1，则有: 整体结构平衡方程 划去节点1所对应的第1行、行1列 。 解得结点位移 2、平面桁架杆单元 （1）单元坐标单元位移向量 看成局部坐标下的拉压杆 （2-15） 应力矩阵 （2-16） 应变矩阵 单元刚度矩阵 （2-14） 形函数 等效节点力 静力等效 3、空间杆单元 （1）单元坐标单元位移向量 (2) 形函数 (3) 应变矩阵 (4) 应力矩阵 (5) 等价节点力 (6) 单元坐标单元刚度矩阵 对于等