平面问题有限元PPT讲稿.ppt
《平面问题有限元PPT讲稿.ppt》由会员分享,可在线阅读,更多相关《平面问题有限元PPT讲稿.ppt(78页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、平面问题有限元第1页,共78页,编辑于2022年,星期日3-13-1有限单元法的计算步骤有限单元法的计算步骤3-23-2平面问题的常应变平面问题的常应变(三角形三角形)单元单元3-33-3单元刚度矩阵单元刚度矩阵3-43-4单元刚度矩阵的物理意义及其性质单元刚度矩阵的物理意义及其性质3-53-5平面问题的矩形单元平面问题的矩形单元3-63-6六节点三角形单元六节点三角形单元3-73-7单元载荷移置单元载荷移置3-83-8整体分析整体分析3-93-9整体刚度矩阵的形成整体刚度矩阵的形成3-103-10整体刚度矩阵的特点整体刚度矩阵的特点3-113-11支承条件的处理支承条件的处理3-123-12
2、应力计算应力计算1/19/2023第2页,共78页,编辑于2022年,星期日弹性力学平面问题的有限单元法包括五个主要步骤:弹性力学平面问题的有限单元法包括五个主要步骤:1 1、所分析问题的数学建模、所分析问题的数学建模 2 2、离散化、离散化 3 3、单元分析、单元分析 4 4、整体分析与求解、整体分析与求解 5 5、结果分析、结果分析1/19/2023图 3-1第3页,共78页,编辑于2022年,星期日有限单元法的基础是用所谓有限个单元的集合体来有限单元法的基础是用所谓有限个单元的集合体来代替原来的连续体,因而必须将连续体简化为由有限个单代替原来的连续体,因而必须将连续体简化为由有限个单元组
3、成的离散体。元组成的离散体。对于平面问题,最简单,因而最常用的单元是对于平面问题,最简单,因而最常用的单元是三角三角形单元形单元。因平面问题的变形主要为平面变形,故平面上所有的节因平面问题的变形主要为平面变形,故平面上所有的节点都可视为平面铰,即每个节点有两个自由度。单元与单元在点都可视为平面铰,即每个节点有两个自由度。单元与单元在节点处用铰相连,作用在连续体荷载也移置到节点上,成为节节点处用铰相连,作用在连续体荷载也移置到节点上,成为节点荷载。如节点位移或其某一分量可以不计之处,就在该节点点荷载。如节点位移或其某一分量可以不计之处,就在该节点上安置一个铰支座或相应的连杆支座。如图上安置一个铰
4、支座或相应的连杆支座。如图3-13-11/19/2023第4页,共78页,编辑于2022年,星期日1/19/20231、位移函数、位移函数 如果弹性体的位移分量是坐标的已知函数,则可用几何方程求应变分如果弹性体的位移分量是坐标的已知函数,则可用几何方程求应变分量,再从物理方程求应力分量。但对一个连续体,内部各点的位移变化情量,再从物理方程求应力分量。但对一个连续体,内部各点的位移变化情况很难用一个简单函数来描绘。况很难用一个简单函数来描绘。有限单元法的基本原理是分块近似,即将弹性体划分成若干细有限单元法的基本原理是分块近似,即将弹性体划分成若干细小网格,在每一个单元范围内,内部各点的位移变化情
5、况可近似地小网格,在每一个单元范围内,内部各点的位移变化情况可近似地用简单函数来描绘。对每个单元,可以假定一个简单函数,用它近用简单函数来描绘。对每个单元,可以假定一个简单函数,用它近似表示该单元的位移。这个函数称为似表示该单元的位移。这个函数称为位移函数位移函数,或称为,或称为位移模式、位位移模式、位移模型、位移场移模型、位移场。对于平面问题,单元位移函数可以用多项式表示,对于平面问题,单元位移函数可以用多项式表示,多项式中包含的项数越多,就越接近实际的位移分布,越精确。多项式中包含的项数越多,就越接近实际的位移分布,越精确。但选取多少项数,要受单元型式的限制。但选取多少项数,要受单元型式的
6、限制。第5页,共78页,编辑于2022年,星期日1/19/2023 三结点三角形单元三结点三角形单元六六个个节节点点位位移移只只能能确确定定六六个个多多项项式式的的系系数数,所所以以平平面面问问题题的的3节节点点三三角角形形单单元元的的位位移移函数如下,函数如下,该该位位移移函函数数,将将单单元元内内部部任任一一点点的的位位移移设设定定为为坐坐标标的的线线性性函函数数,该该位位移移模模式式很很简简单单。其其中中 为为广广义义坐坐标标或或待待定定系系数数,可可据节点据节点i、j、m的位移值和坐标值求出。的位移值和坐标值求出。位移函数写成矩阵形式为:位移函数写成矩阵形式为:第6页,共78页,编辑于
7、2022年,星期日1/19/2023最终确定六个待定系数最终确定六个待定系数其中其中为为2A第第1行行各各个个元元素素的的代代数余子式,数余子式,a(i,j,m),b(i,j,m),c(i,j,m)只是记号,表示此方阵仅与x(i,j,m),y(i,j,m)有关。第7页,共78页,编辑于2022年,星期日1/19/2023令令 (下标(下标i,j,m轮换)轮换)简写为简写为I是单位矩阵,是单位矩阵,N称为形函数矩阵,称为形函数矩阵,Ni只只与与单单元元节节点点坐坐标标有有关关,称称为为单单元元的的形状函数形状函数第8页,共78页,编辑于2022年,星期日据弹性力学几何方程得据弹性力学几何方程得
8、单元的应变分量单元的应变分量由于三节点三角形单元的位移函数为线性函数,则由于三节点三角形单元的位移函数为线性函数,则单元的应变分量均为常量,故这类三角形单元称为单元的应变分量均为常量,故这类三角形单元称为常应变单元(位移在单元内和边界上为线性变化,常应变单元(位移在单元内和边界上为线性变化,应变为常量)应变为常量)1/19/2023第9页,共78页,编辑于2022年,星期日2 2、形函数的特点及性质、形函数的特点及性质1)1)形函数形函数N Ni i为为x x、y y坐标的函数,与位移函数有相同的阶次。坐标的函数,与位移函数有相同的阶次。2)2)形函数形函数N Ni i在在i i节点处的值等于
9、节点处的值等于1 1,而在其他节点上的值为,而在其他节点上的值为0 0。即即1/19/20233)单元内任一点的三个形函数之和恒等于单元内任一点的三个形函数之和恒等于1。4)形函数的值在形函数的值在01间变化。间变化。第10页,共78页,编辑于2022年,星期日1/19/20233、收敛性分析、收敛性分析 选择单元位移函数时,应当保证有限元法解答的收敛性,即当网格选择单元位移函数时,应当保证有限元法解答的收敛性,即当网格逐渐加密时,有限元法的解答应当收敛于问题的正确解答。因此,逐渐加密时,有限元法的解答应当收敛于问题的正确解答。因此,选用的位移模式应当满足下列条件:选用的位移模式应当满足下列条
10、件:(1)位移函数必须含单元常量应变。位移函数必须含单元常量应变。(2)单单元元必必须须能能反反映映单单元元的的刚刚体体位位移移(即即单单元元应应变变为为0时时的的位位移移)。前前面面位移函数改写为(注意:位移函数改写为(注意:为为0)则单元刚体位移为则单元刚体位移为显然,位移函数包含显然,位移函数包含了单元的刚体位移了单元的刚体位移(平动和转动)(平动和转动)第11页,共78页,编辑于2022年,星期日(3)(3)位移函数在单元内部必须连续位移。因为线性位移函数在单元内部必须连续位移。因为线性函数,内部连续函数,内部连续 (4)(4)位移函数必须保证相邻单元在公共边界处的位移函数必须保证相邻
11、单元在公共边界处的位移协调(即在公共边界上位移值相同)。位移协调(即在公共边界上位移值相同)。如右图如右图 设公共边界直线方程为设公共边界直线方程为y=Ax+By=Ax+B,代入位移,代入位移函数可得:边界上位移为函数可得:边界上位移为显然,显然,u,vu,v仍为线性函数,即公共边界上位仍为线性函数,即公共边界上位移连续协调。移连续协调。综上所述,常应变三角形单元的位移函数满综上所述,常应变三角形单元的位移函数满足解的收敛性条件,称此单元为协调单元足解的收敛性条件,称此单元为协调单元1/19/2023y=Ax+B 边界不协调产生裂缝边界不协调产生重迭第12页,共78页,编辑于2022年,星期日
12、1/19/2023例题:图示等腰三角形单元,求其形函数矩阵例题:图示等腰三角形单元,求其形函数矩阵N。第13页,共78页,编辑于2022年,星期日1/19/2023 由三角形的面积由三角形的面积把上步求得的a(i,j,m),b(i,j,m),c(i,j,m)代入第14页,共78页,编辑于2022年,星期日4 4、应力、应变矩阵、应力、应变矩阵将位移函数代入平面问题几何方程,得应变矩阵:将位移函数代入平面问题几何方程,得应变矩阵:1/19/2023第15页,共78页,编辑于2022年,星期日应力矩阵应力矩阵由平面问题物理方程得:由平面问题物理方程得:应变矩阵应变矩阵BB反映了单元内任一点的应变与
13、节点位移间的关系反映了单元内任一点的应变与节点位移间的关系应力矩阵应力矩阵SS反映了单元内任一点的应力与节点位移间的关系反映了单元内任一点的应力与节点位移间的关系显然,常应变三角形单元的应变矩阵显然,常应变三角形单元的应变矩阵BB为常量矩阵,说明在该单为常量矩阵,说明在该单元上的应力和应变为常值。由此可见,在相邻单元的边界处,应元上的应力和应变为常值。由此可见,在相邻单元的边界处,应变及应力不连续,有突变。变及应力不连续,有突变。1/19/2023第16页,共78页,编辑于2022年,星期日1/19/2023 yiFixmF xjF xiF ymFy jFmj*yiFi*xmF*xjF*xiF
14、*ymF*yjFmjys*xy*y*xgeexytxs(a)节点力、内部应力(b)虚位移、虚应变 讨论单元内部的应力与单元的节点力的关系,导出用节点位移表示讨论单元内部的应力与单元的节点力的关系,导出用节点位移表示节点力的表达式。节点力的表达式。由应力推算节点力,需要利用平衡方程。采用虚功方程表示出由应力推算节点力,需要利用平衡方程。采用虚功方程表示出平衡方程,即外力在虚位移上所作的虚功等于应力在虚应变上作的平衡方程,即外力在虚位移上所作的虚功等于应力在虚应变上作的虚应变功。虚应变功。第17页,共78页,编辑于2022年,星期日1/19/2023 考虑上图三角形单元的实际受力,节点力和内部应力
15、为:考虑上图三角形单元的实际受力,节点力和内部应力为:任意虚设位移,节点位移与内部应变为任意虚设位移,节点位移与内部应变为第18页,共78页,编辑于2022年,星期日1/19/2023 令实际受力状态在虚设位移上作虚功,外力虚功为令实际受力状态在虚设位移上作虚功,外力虚功为第19页,共78页,编辑于2022年,星期日1/19/2023 计算内力虚功时,从弹性体中截取微小矩形,边长为计算内力虚功时,从弹性体中截取微小矩形,边长为dx和和dy,厚度为,厚度为t,图示微小矩形的实际应力和虚设变形。图示微小矩形的实际应力和虚设变形。第20页,共78页,编辑于2022年,星期日1/19/2023 微小矩
16、形的内力虚功为微小矩形的内力虚功为整个弹性体的内力虚功为整个弹性体的内力虚功为第21页,共78页,编辑于2022年,星期日1/19/2023 根据虚功原理,得根据虚功原理,得这就是弹性平面问题的虚功方程,实质是外力与应力之间的平衡方这就是弹性平面问题的虚功方程,实质是外力与应力之间的平衡方程。程。虚应变可以由节点虚位移求出:虚应变可以由节点虚位移求出:代入虚功方程代入虚功方程第22页,共78页,编辑于2022年,星期日1/19/2023 接上式,将应力用节点位移表示出接上式,将应力用节点位移表示出 有有 令令实际上,单元刚度阵的一般格式可表示为实际上,单元刚度阵的一般格式可表示为 则则建立了单
17、元的节点力与节点位移之间的关系,建立了单元的节点力与节点位移之间的关系,称为单元刚度矩阵。称为单元刚度矩阵。它是它是6*6矩阵,其元素表示该单元的各节点沿坐标方向发生单位位移时矩阵,其元素表示该单元的各节点沿坐标方向发生单位位移时引起的节点力,它决定于该单元的形状、大小、方位和弹性常数,而与引起的节点力,它决定于该单元的形状、大小、方位和弹性常数,而与单元的位置无关,即不随单元或坐标轴的平行移动而改变。单元的位置无关,即不随单元或坐标轴的平行移动而改变。第23页,共78页,编辑于2022年,星期日1/19/2023 由于由于D中元素是常量,而在线性位移模式下,中元素是常量,而在线性位移模式下,
18、B中的元素也是常量,且中的元素也是常量,且因此因此可以进一步得出平面应力问题和平面应变问题中的单元刚度矩阵。可以进一步得出平面应力问题和平面应变问题中的单元刚度矩阵。第24页,共78页,编辑于2022年,星期日单元刚度矩阵单元刚度矩阵 可记为可记为分块矩阵形式分块矩阵形式将应变矩阵将应变矩阵BB的分块阵代的分块阵代入单元刚度矩阵,可得其子入单元刚度矩阵,可得其子块计算式:块计算式:对于常应变三角形单元,对于常应变三角形单元,考虑平面应力问题弹性考虑平面应力问题弹性矩阵矩阵DD,可得,可得1/19/2023第25页,共78页,编辑于2022年,星期日1/19/2023 上述推导单元刚度矩阵的过程
19、可归纳为上述推导单元刚度矩阵的过程可归纳为单元刚阵单元刚阵K的物理意义是单元受节点力作用后抗变形的能力。其元的物理意义是单元受节点力作用后抗变形的能力。其元素素 的意义为:当第的意义为:当第j个自由度发生单位位移,而其他自由度的位个自由度发生单位位移,而其他自由度的位移为移为0时,在第时,在第i个自由度上所施加的力。若按节点来说明,则刚阵中每个自由度上所施加的力。若按节点来说明,则刚阵中每个子块个子块 表示:当节点表示:当节点j处发生单位位移,而其他节点固定时,在节点处发生单位位移,而其他节点固定时,在节点i上所施加的力。上所施加的力。s tABT tA BD BKTe=ede eF D B
20、BDS=(6)(6)(3)(3)(3)(3)(6(63)3)(3(33)3)(3(36)6)(3(36)6)(6(66)6)第26页,共78页,编辑于2022年,星期日1/19/2023 节点力和节点位移的关系:节点力和节点位移的关系:(以简单平面桁架为例以简单平面桁架为例)平面问题中,离散化的单元组合体极为相似,单元组合体在节点平面问题中,离散化的单元组合体极为相似,单元组合体在节点载荷的作用下,节点对单元、单元对节点都有作用力与反作用力载荷的作用下,节点对单元、单元对节点都有作用力与反作用力存在,大小相等方向相反,统称为节点力。存在,大小相等方向相反,统称为节点力。节点力和节点位移的关系前
21、面已经求出:节点力和节点位移的关系前面已经求出:第27页,共78页,编辑于2022年,星期日1/19/2023 单元刚度矩阵的物理意义:单元刚度矩阵的物理意义:将将 写成分块矩阵写成分块矩阵写成普通方程写成普通方程其中其中 表示节点表示节点S(S=i,j,m)产生单位位移时,在节点产生单位位移时,在节点r(r=i,j,m)上所上所需要施加的节点力的大小。需要施加的节点力的大小。第28页,共78页,编辑于2022年,星期日1/19/2023 单元刚度矩阵的物理意义:单元刚度矩阵的物理意义:将节点力列矩阵将节点力列矩阵 与节点位移列矩阵与节点位移列矩阵 均展开成均展开成(6*1)阶列矩阵,单元阶列
22、矩阵,单元刚度矩阵相应地展开成刚度矩阵相应地展开成(6*6)阶方阵:阶方阵:元素元素K的脚码,标有的脚码,标有“-”的表示水平方向,没有标的表示水平方向,没有标“-”的表示垂直方的表示垂直方向。向。第29页,共78页,编辑于2022年,星期日1/19/2023 单元刚度矩阵的物理意义:单元刚度矩阵的物理意义:单元刚度矩阵的每一个元素都有明显的物理意义。单元刚度矩阵的每一个元素都有明显的物理意义。表示节点表示节点S(S=i,j,m)在水平方向、垂直方向产生单位在水平方向、垂直方向产生单位位移时,在节点位移时,在节点r(r=i,j,m)上分别所要施加的水平节点力和垂直节上分别所要施加的水平节点力和
23、垂直节点力的大小。例如点力的大小。例如 表示节点表示节点j在垂直方向产生单位位移时,在在垂直方向产生单位位移时,在节点节点i所需要施加的水平节点力的大小。所需要施加的水平节点力的大小。第30页,共78页,编辑于2022年,星期日1 1)单元刚度矩阵是对称阵,)单元刚度矩阵是对称阵,(只要证明只要证明 )2 2)单元刚阵主对角线元素恒为正值;因为主对角元素)单元刚阵主对角线元素恒为正值;因为主对角元素 表示表示力的方向和位移方向一致,故功总为正值。力的方向和位移方向一致,故功总为正值。3 3)单元刚阵是奇异阵,即)单元刚阵是奇异阵,即|K|=0|K|=0,这是因为计算单元刚阵时没有,这是因为计算
24、单元刚阵时没有对单元的节点加以约束,虽然,单元处于平衡状态,但容许单元产对单元的节点加以约束,虽然,单元处于平衡状态,但容许单元产生刚体位移,故从单元刚度平衡方程不可能得到唯一位移解生刚体位移,故从单元刚度平衡方程不可能得到唯一位移解 ,只能得到唯一的节点力解。,只能得到唯一的节点力解。4 4)单元刚阵所有奇数行的对应元素之和为零,所有偶数行的)单元刚阵所有奇数行的对应元素之和为零,所有偶数行的对应元素之和也为零。由此可见,单元刚阵各列元素的总和为对应元素之和也为零。由此可见,单元刚阵各列元素的总和为零。由对称性可知,各行元素的总和也为零。零。由对称性可知,各行元素的总和也为零。1/19/20
25、23第31页,共78页,编辑于2022年,星期日例题:求下图所示单元例题:求下图所示单元的刚度矩阵,设的刚度矩阵,设1/19/20231、求、求B2、求、求D3、求、求S4、求、求第32页,共78页,编辑于2022年,星期日几点说明:几点说明:1 1)单元刚度方程是满足节点力平衡条件而建立的,即有限元方程是一组)单元刚度方程是满足节点力平衡条件而建立的,即有限元方程是一组节点力平衡方程组。节点力平衡方程组。2 2)单元内任一点位置的平衡条件往往不满足,即微分平衡方程可能不满足。对)单元内任一点位置的平衡条件往往不满足,即微分平衡方程可能不满足。对于非线性单元,位移函数常不满足以位移为未知量的平
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 平面 问题 有限元 PPT 讲稿
限制150内