平面问题的三角形单元课件.pptx

《平面问题的三角形单元课件.pptx》由会员分享，可在线阅读，更多相关《平面问题的三角形单元课件.pptx（45页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、有限元分析有限元分析的基本的基本步骤步骤（1 1）结构离散化结构离散化：将结构分割成有限个单元体：将结构分割成有限个单元体（2 2）选择位移模式选择位移模式：假定内任一点位移可以用单元节：假定内任一点位移可以用单元节点位移来表达，它们之间存在某种简单函数关系。点位移来表达，它们之间存在某种简单函数关系。第1页/共45页由虚位移原理由虚位移原理 可以得到单元的刚度矩阵可以得到单元的刚度矩阵（3 3）计算计算单元刚度矩阵单元刚度矩阵求单元节点位移求单元节点位移与与 节点节点内力的关系内力的关系求出节点内力：求出节点内力：有限元分析的基本步骤第2页/共45页（5 5）由于上述总刚度矩阵常常是）由于上

2、述总刚度矩阵常常是奇异矩阵奇异矩阵，无，无法求解。法求解。引入边界条件引入边界条件，求解整个结构的所有，求解整个结构的所有单元单元节点的位移节点的位移 节点应变节点应变 节点节点应力应力。（4 4）集合所有节点的平衡方程集合所有节点的平衡方程，形成整个，形成整个结构的平衡方程组，结构的平衡方程组，有限元分析的基本步骤有限元分析的基本步骤第3页/共45页有限元的单元分析单元分析第4页/共45页有限元分析实例求解有限元分析实例求解通过材料力学，弹性力学和有限元法分别求解对比：例：等截面直杆在自重作用下的拉伸 图(a)(a)单位杆长重量为q q，杆长为L L，截面面积为A A，弹性模数为E E 第5

3、页/共45页材料力学求解方法材料力学求解方法第6页/共45页材料力学求解方法材料力学求解方法根据力平衡条件有：内力 N(x)=q(L-x)取微元 dx，则其伸长为 x截面上的位移：根据几何方程求应变，物理方程求应力。这里 应变 根据本构方程求应力第7页/共45页材料力学求解方法材料力学求解方法所以三个节点处的位移函数如下：所以三个节点处的应变函数如下：所以三个节点处的应力函数如下：第8页/共45页弹性力学求解方法弹性力学求解方法+dN+dNqdx第9页/共45页弹性力学求解方法弹性力学求解方法微元微元力平衡方程力平衡方程：微元微元几何方程几何方程：材料材料本构方程本构方程：给出给出边界方程边界

4、方程：第10页/共45页有限元法求解有限元法求解有限单元法求解直杆拉伸：1 1、离散化（节点和单元、离散化（节点和单元）2 2、外载荷集中到节点上，即把阴、外载荷集中到节点上，即把阴影部分的重量作用在节点影部分的重量作用在节点i i上上 第11页/共45页有限元法求解有限元法求解3 3、假设线单元上的位移为线性函数、假设线单元上的位移为线性函数 由几何方程：由本构方程：节点内力：第12页/共45页有限元法求解有限元法求解有限单元法求解直杆拉伸：4 4、以、以i i节点节点为对象，列力的平衡方程为对象，列力的平衡方程令令 将位移和内力的关系代入得将位移和内力的关系代入得 用节点位移表示的平衡方程

5、，其中用节点位移表示的平衡方程，其中i=1i=1，2 2，n n有有n n个方程个方程未知数也有未知数也有n n个，解方程组，得出个，解方程组，得出节点节点位移，进而计算应力位移，进而计算应力 第13页/共45页2-1 有限单元法的概念有限单元法的概念有限单元法求解直杆拉伸：假设线单元数为假设线单元数为3 3个的情况，个的情况，平衡方程有平衡方程有3 3个：个：i=1i=1时，时，i=2i=2时时,i=3i=3时，时，联立解得联立解得 与材料力学的精确解答在节点处与材料力学的精确解答在节点处位移位移完全相同，回代可以得到完全相同，回代可以得到各点各点应变值应变值，继续回代可以得到各点，继续回代

6、可以得到各点应力值。应力值。第14页/共45页有限元法求解有限元法求解设取n=3，求解含节点位移的线性方程组，得各点位移如下 第15页/共45页第16页/共45页有限元的单元分析单元分析第17页/共45页1 1 三角形单元位移插值函数三角形单元位移插值函数假设已知假设已知如何求单元内(x,y)点位移？第18页/共45页1 1 三角形单元位移插值函数三角形单元位移插值函数将i，j，m节点坐标（已知）代入上式得含待定系数的方程组选择位移插值函数如下：选择位移插值函数如下：第19页/共45页代入上述位移函数可得：求解6个待定系数第20页/共45页其中A为三角形面积将将待定系数待定系数代入单元内部位移

7、模式得到代入单元内部位移模式得到任意点位移：任意点位移：第21页/共45页式中：式中：进一步进一步简化，简化，令令位移形函数位移形函数单元内部位移模式可以简写为：单元内部位移模式可以简写为：第22页/共45页单元内部位移模式的矩阵表达式：单元内部位移模式的矩阵表达式：单元内部位移模式的矩阵表达式可以简记为：单元内部位移模式的矩阵表达式可以简记为：位移位移转换转换矩阵函数矩阵函数或或位移形函数矩阵位移形函数矩阵第23页/共45页单元内部位移模式的矩阵表达式：单元内部位移模式的矩阵表达式：位移形位移形函数函数Ni物理物理含义含义故故,Ni称为称为位移位移形形函数。函数。第24页/共45页单元内部位

8、移模式单元内部位移模式必须满足三个条件才能必须满足三个条件才能保证收敛：保证收敛：注意注意：如如第25页/共45页2 2 由由节点位移节点位移求求应变应变 几何方程几何方程第26页/共45页2 2 由由节点位移节点位移求求应变应变 几何方程几何方程第27页/共45页式中：式中：第28页/共45页3 3 由由应变应变求求应力应力 本构方程本构方程将应变矩阵代入上式将应变矩阵代入上式第29页/共45页已已知知：单单元元节节点点力力和和节节点点虚虚位位移移，节点力所做的虚功节点力所做的虚功W为：为：4 4 由由应力应力求求节点力节点力 虚功方程虚功方程第30页/共45页已已知知：单单元元内内部部应应

9、力力和和虚虚应应变变，则则整个弹性体内的变形虚功整个弹性体内的变形虚功U为为第31页/共45页第32页/共45页单元刚度矩阵为单元刚度矩阵为 将虚应变矩阵将虚应变矩阵代入上式并整理代入上式并整理再将应力矩阵代入得再将应力矩阵代入得第33页/共45页第34页/共45页式中：式中：第35页/共45页位移函数位移函数 几何方程几何方程 物理方程物理方程 虚功方程虚功方程单元刚度矩阵单元刚度矩阵k单元分析小结假设假设节点位移已知已知可以表达可以表达节点内力第36页/共45页单元分析小结也就是说也就是说每个单元每个单元对对节点节点贡献贡献的的力力知道了知道了第37页/共45页取取节节点点i分分析析，单单

10、元元，和和的的共共用用节节点点i，所所以以节节点点i的的内内力力为为三三个个单单元元分力在该点之和：分力在该点之和：5 5 节点平衡方程组节点平衡方程组 整体刚度矩阵整体刚度矩阵第38页/共45页由由于于每每个个节节点点有有两两个个未未知知位位移移分分量量，所所以以根根据据下下列列两两个个平平衡衡方程方程理论上，可以求解。理论上，可以求解。5 5 节点平衡方程组节点平衡方程组 整体刚度矩阵整体刚度矩阵第39页/共45页 从从网网格格节节点点的的整整体体来来看看，本本问问题题共共计计6个个节节点点，每每个个节节点点有有两两个个位位移移分分量量，共共计计12个个（未未知知）位移位移分量分量。每每个

11、个单单元元分分析析都都可可以以用用12个个节节点点位位移移中中的的6个个，描描述述该单元的节点内力：节点内力：也也就就是是说说所所有有单单元元的的节节点点内内力力都都能用能用12个位移未知量来表达。个位移未知量来表达。5 5 节点平衡方程组节点平衡方程组 整体刚度矩阵整体刚度矩阵第40页/共45页 列列出出所所有有节节点点的的内内、外外力力平平衡衡方方程程：准准确确的的说说是是12个个方方程程可可以以求求解解12个个未未知知量量（可可能能是是位移也可能是外力位移也可能是外力）。）。注注意意：边边界界上上的的节节点点，有有些些位位移移是是已已知知的的，有有些些是是外外力力已已知知的的。如如果果没

12、没有有边边界界条条件件，方方程程会有无穷多个解。会有无穷多个解。5 5 节点平衡方程组节点平衡方程组 整体刚度矩阵整体刚度矩阵第41页/共45页平面平面问题三角形单元问题三角形单元有限元分析小结有限元分析小结有限元分析的基本步骤：有限元分析的基本步骤：（1 1）结构离散化结构离散化：将结构分割成有限个单元体：将结构分割成有限个单元体（2 2）选择位移模式选择位移模式：假定内任一点位移可以用单元节：假定内任一点位移可以用单元节点位移来表达，它们之间存在某种简单函数关系。点位移来表达，它们之间存在某种简单函数关系。第42页/共45页由虚位移原理由虚位移原理 可以得到单元的刚度矩阵可以得到单元的刚度

13、矩阵（3 3）计算计算单元刚度矩阵单元刚度矩阵求单元节点位移求单元节点位移与与 节点节点内力的关系内力的关系求出节点内力：求出节点内力：平面问题三角形单元有限元分析小结第43页/共45页（5 5）由于上述总刚度矩阵常常是）由于上述总刚度矩阵常常是奇异矩阵奇异矩阵，无，无法求解。法求解。引入边界条件引入边界条件，求解整个结构的所有，求解整个结构的所有单元单元节点的位移节点的位移 节点应变节点应变 节点节点应力应力。（4 4）集合所有节点的平衡方程集合所有节点的平衡方程，形成整个，形成整个结构的平衡方程组，结构的平衡方程组，平面平面问题三角形单元问题三角形单元有限元分析小结有限元分析小结第44页/共45页感谢您的观看。第45页/共45页