有限元分析平面问题有限单元法

有限元分析平面问题有限单元法第一页，共二十七页，2022年，8月28日内容Chp.4平面问题的有限元法1

4.1连续体的离散化4.2常应变三角形单元4.3形函数的性质4.4刚度矩阵要求理解：连续体有限元分片插值的含义

有限元分析中的几个常用概念，包括：

形函数，应变矩阵，应力矩阵，刚度矩阵等

掌握：常应变三角形单元分析流程课后作业推导常应变三角形单元各矩阵有限元分析——4－1,2,3,4第二页，共二十七页，2022年，8月28日回顾刚架有限元分析的基本流程

整体离散梁单元分析单元组装整体解算刚架结构结构整体刚度方程PP1432123456eijTjmjqjTimiqiO’单元内力eij△jfj△ifiy’x’O’单元变形单元刚度方程材料力学第三页，共二十七页，2022年，8月28日

离散化方法

节点：i，j单元编号ePP1432123456单元(element):节点(node):1,2,3,44123456eij其中：e=1,2,3,4,5,6;i,j=1,2,3,4回顾天然划分出各梁单元第四页，共二十七页，2022年，8月28日梁单元分析由材料力学得：单元节点力与单元节点位移之间的关系局部坐标系单元受力eijTjmjqjTimiqiO’局部坐标系单元变形eij△jfj△ifiy’x’O’回顾第五页，共二十七页，2022年，8月28日总体刚度矩阵的组装PP1432123456组装原理：位移协调条件节点平衡条件位移协调条件：各单元共享节点的位移相等节点平衡条件：各节点单元内力与节点外力构成平衡力系回顾最终数学模型：第六页，共二十七页，2022年，8月28日基本概念

回顾单元（element）节点(node)单元节点位移(nodedisplacement)单元节点内力(nodeforce)单元刚度矩阵(elementstiffnessmatrix)整体刚度矩阵组装原理总体刚度矩阵，载荷向量，刚度方程第七页，共二十七页，2022年，8月28日新问题：连续体的有限元分析如何离散化？

如何单元分析？

如何组装单元？

第八页，共二十七页，2022年，8月28日一、连续体的离散化人工节点逼近性离散节点位置，单元形状、大小都会影响对原始结构的逼近描述平面问题:三角形、矩形、任意四边形…空间问题:四面体、正方体、任意六面体…第九页，共二十七页，2022年，8月28日单元分析的理论基础：二、平面问题三角形单元分析建立三角形单元节点力与单元节点位移之间的关系单元分析的目的：弹性力学基本理论vmumvjviuii(xi,yi)j(xj,yj)m(xm,ym)eujyxoFymFxmFyjFyiFxii(xi,yi)j(xj,yj)m(xm,ym)eFxjyxo单元节点位移单元等效节点力(基本变量，基本方程，边界条件)～？第十页，共二十七页，2022年，8月28日二、平面问题三角形单元分析平面问题弹性力学基本理论回顾：基本变量几何方程平衡方程物理方程平面应力问题平面应力问题第十一页，共二十七页，2022年，8月28日单元变形分析（位移）问题：平面三角形问题中仍包含无穷多内部连续点，如何用节点位移表征内部连续位移？二、平面问题三角形单元分析解决思路：插值，即用节点位移插值近似其它点的位移。vmumvjviuii(xi,yi)j(xj,yj)m(xm,ym)eujyxo如何插值？多项式中包含的项数越多，就越接近实际的位移分布，越精确。但选取多少项数，要受单元型式的限制。第十二页，共二十七页，2022年，8月28日二、平面问题三角形单元分析单元变形分析（位移）vmumvjviuii(xi,yi)j(xj,yj)m(xm,ym)eujyxo插值系数的确定：待定系数法6个方程6个系数第十三页，共二十七页，2022年，8月28日二、平面问题三角形单元分析矩阵表达和运算令A为三角形单元的面积

[T]*为[T]的伴随矩阵：第十四页，共二十七页，2022年，8月28日二、平面问题三角形单元分析矩阵表达和运算（续）令同理，由回代，得第十五页，共二十七页，2022年，8月28日二、平面问题三角形单元分析三角形单元形函数令得单元内各点位移：单元节点位移：形函数第十六页，共二十七页，2022年，8月28日

例题：图示等腰三角形单元，求其形函数矩阵[N]。ijm第十七页，共二十七页，2022年，8月28日

由三角形的面积第十八页，共二十七页，2022年，8月28日二、平面问题三角形单元分析单元变形分析（应变）位移应变代入记为[B]矩阵称为应变矩阵。

常数第十九页，共二十七页，2022年，8月28日二、平面问题三角形单元分析单元应力分析应变应力代入记为[S]矩阵称为应力矩阵。

[D]为弹性矩阵第二十页，共二十七页，2022年，8月28日单元平衡分析二、平面问题三角形单元分析FymFxmFyjFyiFxii(xi,yi)j(xj,yj)m(xm,ym)eFxjyxo问题：假设隔离出的各单元，仅在节点上有内力的作用（若单元边界上和体内有力的作用，可事先等效到节点上，如何等效？），则单元节点内力与变形产生的应力构成平衡力系。单元满足平衡条件平衡方程用虚功方程代替。第二十一页，共二十七页，2022年，8月28日单元平衡分析二、平面问题三角形单元分析平面问题虚功原理代入得内力虚功＝第二十二页，共二十七页，2022年，8月28日二、平面问题三角形单元分析单元平衡分析外力虚功＝载荷等效公式第二十三页，共二十七页，2022年，8月28日二、平面问题三角形单元分析单元平衡分析外力虚功＝内力虚功令[K]e单元刚度矩阵单元刚度方程第二十四页，共二十七页，2022年，8月28日连续体三角形单元分析小结单元位移：节点位移插值单元刚度方程：由虚功方程建立离散化：单元分析：单元节点人工单元、节点单元应变：由几何方程建立单元应力：由物理方程建立下次课内容：单元组装和单元等效节点外载第二十五页，共二十七页，2022年