有限单元法电子课件(桁架).ppt

有限单元法(Finite Element Method) * 基本内容 第一章 绪论 第二章 杆系有限元分析 第三章 弹性力学基本方程 第三章 弹性力学平面问题有限元 第四章 弹性力学空间问题有限元 第六章 ANSYS简介 绪 论 ● 将连续体（结构）离散为若干子区域 子区域由结点连接为等效的组合体 ● 每个单元内假设场变量为多项式(系数不同） 用分区域连续场函数近似全区域的连续场函数 无穷自由度问题转化为有限自由度问题 ● 利用变分原理得到离散场变量的大代数方程组 将微分方程边值问题转化为代数方程来求解 杆系结构 连续体 ▲ 有限元法的应用 ● 工程分析中应用最广泛的一种数值计算方法 ● 科学与工程计算、计算机辅助设计(CAD)和制造（CAM) ▲ 有限元法的要点 绪 论 ● Courant (1943) 三角形单元分片连续函数结合最小时能远离三角形 单元求成功解 St. Venant 扭转问题。 ● Turner, Clough 等（1956）首次采用 三角形单元求解平面应力问题。● Clough(1960)进一步求解弹性力学平面问题，提出“有限单元法” （Finite Element Method）的名称。 ▲ 有限元的特性 ● 对于复杂几何形状的适应性 ● 对于各种物理问题的可应用性 ▲ 有限元法发展历史 ● 建立于严格理论基础上的可靠性 ● 适合计算机实现的高效性(同步发展） 绪 论 ▲ 有限元的发展现状 ● 单元种类更加丰富，适应于更多更复杂的工程和科学计算 ● 理论基础更加完善(广义变分原理） ● 计算机科学与技术的发展为高效快速计算提供了保证 ● 大型商业软件的开发和推广和升级促进了有限元法的广泛应用 ANASYS, ABAQUS, SAP, MARC, NASTRAN …… ● 工程技术领域应用最广泛、成效最显著的数值分析方法。 计算机辅助设计(CAD)和制造（CAM) 杆系结构 (Structure of line elements) 桁架(truss) 平面桁架(plane truss) 空间桁架(spacial truss) 刚架(frame) 平面刚架(plane frame) 空间刚架(spacial frame) 杆系结构 桁架结构的特点 结构特点 : 由光滑铰链连接的接直杆 载荷特点：载荷只作用在结点处 约束特点：光滑铰链，只约束线位移 变形特点：各杆只有轴向拉伸和压缩 单元特点：轴向拉（压）杆 类型：空间桁架，平面桁架 平面桁架（Trusses）有限元分析（1） 1、单元的位移形函数 、局部坐标系中的单元分析 图2 局部坐标系中的杆单元 图3 单元形函数（线性）示意图 形函数 （shape function） 图1 平面桁架示意图 单元的结点位移 平面桁架（Trusses）有限元分析（2） 2、单元应变 --- 几何矩阵 3、单元应力 其中E为弹性模量, [C]=E[B] --- 应力矩阵 4、单元刚度矩阵 局部坐标系中的单元刚度矩阵 单元结点位移 单元的结点力 平面桁架（Trusses）有限元分析（3） 整体坐标系中的结点载荷向量 1、结点力矢量的坐标转换 或 图4 单元结点力的坐标转换 坐标转矩阵 二、整体坐标系中的单元分析 2、结点位移矢量的坐标转换 整体坐标系中的结点位移 平面桁架（Trusses）有限元分析（4） 3、反变换 4、整体坐标系中的单元刚度矩阵 整体坐标系中的单元刚度矩阵 可记为： 或： 平面桁架（Trusses）有限元分析（5） 1、整体刚度矩阵的形成 三、整体结构的分析 图5 三个杆件组成的简单桁架 整体分析的目的是建立整体刚度方程 整体结点力 整体结位移 整体刚度矩阵 第 j 点的位移对第 i 点力的贡献 平面桁架（Trusses）有限元分析（6） 1、整体刚度矩阵的形成 表1 单元的局部码和整体码的对应关系 2 3 2 j 3 1 1 i (3) (2) (1) 两种编码 对应关系 同阶）分别为 形成规律：整体刚度是所有单元刚度对整 体刚度的贡献。 如何考虑这个贡献？ 图5 三个杆件组成的简单桁架 平面桁架（Trusses）有限元分析（8） 整体刚度矩阵 同阶）分别为 整体刚度方程： 该方程不可解，矩阵奇异！ 图5 三个杆件组成的简单桁架 为什么？？？ 支承条件的必要性。 平面桁架（Trusses）有限元分析（9） 边界条件的引入： 同阶）分别为 整体刚度方程： *