弹性力学及有限元法：第4章 平面问题的有限元法.pptVIP

第四章 平面问题的有限元法 第4章 平面问题的有限元法 4.1 平面三角形单元矩阵推导 4.1 平面三角形单元的单元刚度矩阵推导 4.1 平面三角形单元的单元刚度矩阵推导 4.1 平面三角形单元的单元刚度矩阵推导 4.1 平面三角形单元的单元刚度矩阵推导 4.1 平面三角形单元的单元刚度矩阵推导 4.1 平面三角形单元的单元刚度矩阵推导 4.1 平面三角形单元的单元刚度矩阵推导 4.1 平面三角形单元的单元刚度矩阵推导 4.1 平面三角形单元的单元刚度矩阵推导 4.1 平面三角形单元的单元刚度矩阵推导 4.1 平面三角形单元的单元刚度矩阵推导 4.1 平面三角形单元的单元刚度矩阵推导 4.1 平面三角形单元的单元刚度矩阵推导 4.1 平面三角形单元的单元刚度矩阵推导 4.2 利用平面三角形单元进行结构整体分析 4.2.1 单元的组集 4.2.1 单元的组集 4.2.1 单元的组集 4.2.1 单元的组集 4.2.1 单元的组集 4.2.1 单元的组集 4.2.1 单元的组集 4.2.2 边界条件的引入 4.2.2 边界条件的引入 4.2.2 边界条件的引入 4.2.2 边界条件的引入 4.2.2 边界条件的引入 4.2.2 边界条件的引入 例4-1 如图4-3所示由两个三角形单元所组成的结构系统，试分别采用直接组集法及转换矩阵法进行刚度矩阵的组集，并进行边界条件的引入构件新的有限元方程。 4.2.2 边界条件的引入 解：①采用直接组集法进行整体刚度矩阵的组集 整个系统中共有4个节点，每个节点有两个自由度，则整体刚度矩阵的维数为8。三角形单元刚度矩阵的维数为8*8，因而需要将每个单元扩展成的矩阵再进行叠加。用分块矩阵表示，扩展完成的单元（1）和单元（2）的刚度矩阵可表示为 4.2.2 边界条件的引入 矩阵中各元素的下标表示整体节点编号， 和 则表示单元（1）和单元（2）对整个系统刚度矩阵的贡献。 两者相加则得到总刚矩阵，表示为 4.2.2 边界条件的引入 ②采用转换矩阵法进行整体刚度矩阵的组集 采用转换矩阵法进行整体刚度矩阵组集的关键是获得每个单元的转换矩阵。单元转换矩阵的行数为单元自由度数，转换矩阵的列数为整体刚度矩阵的维数。对于上述结构，每个转换矩阵的维数为，两个单元的转换矩阵分别为 4.2.2 边界条件的引入 获得了每个单元的转换矩阵，则可按 进行整体刚度矩阵的求解。 ③边界条件的引入 图4-3所示结构，引入约束前有限元方程可表示为 4.2 利用平面三角形单元进行整体分析 边界条件分析：由于1, 2节点为固定约束，位移边界条件为u1=0,v1=0,u2=0,v2=0，另外，在节点4的x轴负方向有作用力，对应在整体刚度矩阵、整体位移列阵和整体节点力列阵中对应去掉边界条件中位移为0的行和列，并引入节点作用力，则最终引入边界条件的有限元方程为。 4.3 有限元法的实施步骤 4.3 有限元法的实施步骤 4.3 有限元法的实施步骤 （2）单元位移模式的选择（有限元编程选用步骤） 根据分块近似的思想，选择一个简单的函数来近似地构造每一单元内的近似解。例如，若以节点位移为基本未知量，为了能用节点位移表示单元体的位移、应变和应力，在分析求解时，必须对单元中位移的分布做出一定的假设，即选择一个简单的函数来近似地表示单元位移分量随坐标变化的分布规律，这种函数称为位移模式。位移模式的选择是有限单元法分析中的关键。由于多项式的数学运算比较简单、易于处理，所以通常是选用多项式作为位移模式。多项式的项数和阶数的选择一般要考虑单元的自由度和解答的收敛性要求等。 确定形函数 4.3 有限元法的实施步骤 （3） 单元刚度分析（有限元编程） 通过分析单元的力学特性，建立单元刚度矩阵。首先利用几何方程建立单元应变与节点位移的关系式，然后利用物理方程导出单元应力与节点位移的关系式，最后由虚功原理或最小势能原理推出作用于单元上的节点力与节点位移之间的关系式，即单元刚度矩阵。（实际上有多种方法，直接利用刚度矩阵公式，利用已推导的刚度矩阵元素解析表达式，高斯积分，从ANSYS导入，。。。） （4） 等效节点力计算 分析对象经过离散化以后，单元之间仅通过节点进行力的传递，但实际上力是从单元的公共边界上传递的。为此，必须把作用在单元边界上的表面力，以及作用在单元上的体积力、集中力等，根据静力等效的原则全都移置到节点上，移置后的力称为等效节点力。 4.3 有限元法的实施步骤 (5)整体结构平衡方程建立 建立整体结构的平衡方程也叫做结构的整体分