数值模拟与仿真基础——3 平面问题的有限元法.pptVIP

数值模拟与仿真技术基础 麻向军 华南理工大学 聚合物新型成型装备国家工程研究中心 聚合物成型加工工程教育部重点实验室 3、平面问题的有限元法 3.0 引言 3.1 位移函数 3.2 单元刚度方程 3.3 载荷移置与等效节点载荷 3.4 结构刚度方程 3.5 位移边界条件的处理 3.6 应力计算 3.7 解题示例与公式推广 3.8 斜边界问题的处理 3.9 六节点三角形单元 3.0 引　言 1、为什么先进行平面问题的有限元法 2、选用的单元类型及特点 3、三角形单元的网格剖分原则 4、节点编号的约定 5、三角形单元划分的示例 1、为什么先进行平面问题的有限元法 (1)、平面问题的有限元分析较简单，具有典型性 (2)、在工程应用中有其实际意义，主要表现在在满足工程精度的要求下，降低问题的复杂性，提高分析问题的效率。 (3)、平面问题的有限元分析是今后进一步分析轴对称问题、三维问题及板壳问题的基础。 从平面问题的有限元法分析入手，可有利于有限元基本概念、方法、理论的理解与掌握。 2、选用的单元类型及特点 进行平面问题研究时，选用三角形单元较简单。三节点的三角形单元又是最简单而又被广泛采用的一种单元类型。 　　由于在平面问题分析中，结构发生的是平面变形，三角形的三个节点可以看作是平面铰，每个节点具有两个自由度，这样共有３个节点６个自由度，如果节点位移或其中某一个分量为零时，可在该节点处设置一个平面铰支座或连杆支座，以限制其位移。由三角单元离散的结构是由三角形单元的节点铰接而成的。 3、三角形单元的网格剖分原则 (1)、各节点必须相连。 (3)、单元的大小、数目取决于计算精度的要求和计算容量的限制 划分网格时，首先要满足计算精度的要求，同时可利用结构的对称性，循环对称性的特点，从厚结构中取出一部分进行分析，或者对有应力集中的构件，采用疏密不同的网格剖分。也可以采用子结构法。 4、节点编号的约定 (1)、节点编号分为局部节点编号和总体节点编号两种。如下图中的矩形，分为 5 个节点，4 个单元，其中1、2、3、4、5 为总体节点编号。 对于任一单元如 ②中，3、5、4 为局部节点编号，在公式推导中用 i、j、m 编号，约定其为逆时针顺序。主要是要保证用 i、j、m 节点计算的单元面积为正值。 5、三角形单元划分的示例 3.1 位移函数 1、位移函数的一般形式 2、三节点三角形单元的位移函数 3、位移函数及其性质 4、位移函数与解的收敛性 1、位移函数的一般介绍 (1)、位移函数的一般介绍 结构离散化后，要对单元进行力学特性分析，即确定单元节点力与节点位移之间的关系，这时就需要把单元内的任一点的位移分量表示成坐标的某种函数。这种函数就叫位移函数。 (5)、位移函数一般形式为多项式： ① 进行数学运算（如微分，积分）较简单； ② 任意阶次的多项式可以近似地表示精确解。 2、三节点三角形单元的位移函数 (1)、位移函数形式 六个节点位移只能确定六个多项式的系数，所以平面问题的 3 结点三角形单元的位移函数如下： (2)、待定系数 β1, β2, … β6 的求解 克来姆法则 　 设有一线性方程组： 由克来姆法则可知：当 2A?0，上述方程有唯一解： 为了描述方便，引入系数 (4)、位移函数的插值函数形式 写成矩阵形式： (5) 例题：图示等腰三角形单元的形态矩阵[N]。 3、位移函数及其性质 当节点位移一定时，单元形态完全决定于Ｎi ,Ｎj ,Ｎm。因而形函数就具有如下的性质： 　　性质1、形函数 Ｎi 在节点 i 处的值为１，而在其它两个节点 (j , m) 处的值为零。 性质2、在单元任一节点处，三个形函数之和等于1。 4、位移函数与解的收敛性 (1)、位移函数应满足的条件： 选择位移函数时，为保证有限元法的收敛性，必须满足以下4个条件： ① 位移函数必须包含单元的常量应变 ② 位移函数必须包含单元的刚体位移 ③ 位移函数在单元内部必须是连续函数(连续性要求) ④ 位移函数应使相邻单元间的位移协调(保续性要求) 　　上述四个条件中，若