《弹性力学教学课件》1-有限元位移法的基本概念.ppt

有限元位移法的基本概念有限元位移法的基本概念 一、一个连续梁的简例一、一个连续梁的简例 单元分析单元分析 单元单元：单元单元：整体分析整体分析 即：即：基本思路：基本思路：有限单元法：有限单元法：“先分后合先分后合”先分：首先将连续梁分成单元，进行单元分析，用结点位移先分：首先将连续梁分成单元，进行单元分析，用结点位移B 表示单元内力表示单元内力MBA、MBC。后合：然后将单元再合成结构，进行整体分析，建立结构整体后合：然后将单元再合成结构，进行整体分析，建立结构整体 平衡条件平衡条件MB0，由此可以求出结点位移，由此可以求出结点位移B。（位移法基本方程）（位移法基本方程）单元刚度系数的概念：单元刚度系数的概念：当结点当结点B产生单位转角产生单位转角B 1时，在时，在单元单元所需施加的杆端力矩为所需施加的杆端力矩为k(1)、在单元在单元所需施加的杆端力矩为所需施加的杆端力矩为k(2)。单元单元：单元单元：其中：其中：称为单元称为单元和单元和单元的的 刚度系数。刚度系数。整体刚度系数的概念整体刚度系数的概念 位移法基本方程可写成：位移法基本方程可写成：其中：其中：K称作结构的整体刚度系数，即整个结构在结点称作结构的整体刚度系数，即整个结构在结点B产生单位位产生单位位 移移B 1时，在结点时，在结点B所需施加的结点力偶。所需施加的结点力偶。刚度集成规律刚度集成规律 可由相关单元的单元刚度系数集成。可由相关单元的单元刚度系数集成。二、单元刚度矩阵二、单元刚度矩阵 单元分析：单元分析：单元两个端点重新编码为单元两个端点重新编码为1、2。由转角位移方程得：由转角位移方程得：其中：其中：系数系数k(e)的物理意义：的物理意义：这些系数称作单元这些系数称作单元(e)的刚度系数。的刚度系数。将转角位移方程写成矩阵形式：将转角位移方程写成矩阵形式：杆端力矩列矩阵杆端力矩列矩阵 结点转结点转角列矩角列矩阵阵单元刚度矩阵单元刚度矩阵 三、整体刚度矩阵三、整体刚度矩阵 结结点点转转角角 结结点点力力偶偶 它们之间的转换关系可表示为：它们之间的转换关系可表示为：或：或：K称为整体刚度矩阵称为整体刚度矩阵 K中的元素中的元素Kij称为整体刚度系数，它表示当结点称为整体刚度系数，它表示当结点 j 产生单位转产生单位转角角j 1(其它结点没有转角其它结点没有转角)时，在结点时，在结点 i 所需施加的力偶所需施加的力偶Mi。先分别考虑每个结点单独转动，先分别考虑每个结点单独转动，然后叠加。然后叠加。结点结点1单独转动单独转动1：结点结点2单独转动单独转动2：结点结点3单独转动单独转动3：先分别考虑每个结点单独转动，先分别考虑每个结点单独转动，然后叠加。然后叠加。即：即：或：或：先分别考虑每个结点单独转动，先分别考虑每个结点单独转动，然后叠加。然后叠加。即：即：或：或：先分别考虑每个单元单独变形，先分别考虑每个单元单独变形，然后叠加。然后叠加。首先写出单元刚度矩阵，得出单元刚度方程如下：首先写出单元刚度矩阵，得出单元刚度方程如下：其次进行叠加，叠加前对单元刚度矩阵进行改造：其次进行叠加，叠加前对单元刚度矩阵进行改造：单元的贡献矩阵：单元的贡献矩阵：先分别考虑每个单元单独变形，先分别考虑每个单元单独变形，然后叠加。然后叠加。单元的贡献矩阵：单元的贡献矩阵：最后，将各单元的贡献矩阵最后，将各单元的贡献矩阵叠加，即得出整体刚度矩阵叠加，即得出整体刚度矩阵，刚度方程为：，刚度方程为：四、支承条件的引入四、支承条件的引入 第一步：暂不引入支承条件，第一步：暂不引入支承条件，和荷载情况，先建立整体刚和荷载情况，先建立整体刚度矩阵度矩阵K，写出，写出M与与之之间的转换关系。间的转换关系。第二步：在结点第二步：在结点1和和2引入荷引入荷 载值载值P1和和P2，在固定端引入，在固定端引入支承条支承条 件件30，将整体刚，将整体刚度方程度方程 改写为矩阵方程：改写为矩阵方程：取出前面两个纯量方程：取出前面两个纯量方程：即：即：五、非结点荷载的处理五、非结点荷载的处理 求等效结点荷载求等效结点荷载 第一步：在各结点加约束，第一步：在各结点加约束，阻止结点转动。阻止结点转动。各单元产生的固端力矩：各单元产生的固端力矩：单元单元单元单元各结点的约束力矩分别为该各结点的约束力矩分别为该结点相关单元固端力矩之和：结点相关单元固端力矩之和：第二步：去掉各结点的约束。第二步：去掉各结点的约束。五、非结点荷载的处理五、非结点荷载的处理 求等效结点荷载求等效结点荷载 第一步：在各结点加约束，第一步：在各结点加约束，阻止结点转动。阻止结点转动。各单元产生的固端力矩：各单元产生的固端力矩：单元单元单元单元 求各杆端力矩求各杆端力矩 由于：由于：所以：所以：