第三章弹性力学有限元法

.第三章弹性力学有限元法.第三章弹性力学有限元法3.1有限元法求解问题的基本步骤3.2连续体离散化3.3单元分析3.4整体分析3.5边界约束条件处理3.6求解、计算结果的整理和有限元后处理.3.1有限元法求解问题的基本步骤A.工程结构的几何简化B.载荷简化C.边界条件的简化1.建立力学模型2.连续体离散化用合适的单元将连续体划分为有限个具有规则形状的的单元集合，单元的选取应视所分析问题的性质、规模和精度要求而定。.3.1有限元法求解问题的基本步骤3.单元分析A.位移模式(位移函数、插值函数）选取B.单元分析单元的位移模式一般采用多项式，因为多项式计算简便，并且随着项数的增加，可以逼近任何一段光滑的函数曲线。a.建立单元刚度矩阵b.建立单元节点力列阵弹性力学几何方程和物理方程静力等效原则.3.1有限元法求解问题的基本步骤4.整体分析和有限元方程求解A.建立整体刚度矩阵B.建立整体节点力列阵C.代入边界条件D.选择适当的代数方程求解a.高斯消元法b.三角分解法c.波前法d.雅克比迭代法.3.1有限元法求解问题的基本步骤5.结果后处理和分析A.应力误差的减小B.结果输出方式C.结果分析.3.2连续体离散化1.杆状单元uivjviujji平面杆单元只能承受轴向的拉压载荷，平面杆单元每个节点有两个自由度。空间杆单元只能承受轴向的拉压载荷，平面杆单元每个节点有三个自由度。uivjviujjiwiwjA.轴力单元.3.2连续体离散化B.梁单元平面梁单元平面梁单元每个节点有三个自由度，两个线位移一个角位移，可承受平面内的体力，集中力、分布力和垂直平面的弯矩的作用。空间梁单元空间梁单元每个节点有六个自由度，三个线位移三个角位移，可承受各个方向的体力，集中力、分布力和弯矩的作用。juiujwjwiθiθjiuivjviujjiwiwjθziθxiθyiθxiθzi1.杆状单元.3.2连续体离散化2.平面单元ijluivjviujulvl平面单元属于二维单元，只能承受单元平面内的分布力和集中力，不能承受面外载荷。三节点三角形单元每个节点有两个自由度，因此只能采用线性模式，应变值为常量，也称为常应变单元或常应力单元。.3.2连续体离散化四节点四边形单元四边形单元有四个节点，每个节点也有两个自由度，采用双线性位移模式，计算精度较高。iuivijvjujlulvlmumvm2.平面单元.3.2连续体离散化3.薄板弯曲单元和薄板单元A.薄板弯曲单元wiθyiθxiijml四边形弯曲单元四边形单元有四个节点，每个节点有三个自由度，主要承受横向载荷和绕水平轴的弯矩。.3.2连续体离散化A.薄板弯曲单元三角形弯曲单元三角形单元有三个节点，每个节点有三个自由度，主要承受横向载荷和绕水平轴的弯矩。jwiθyiθxiim3.薄板弯曲单元和薄板单元.3.2连续体离散化3.薄板弯曲单元和薄板单元B.薄板单元四边形薄板单元四边形薄板单元有四个节点，每个节点有五个自由度，可承受各个方向的载荷和绕水平轴的弯矩。wiθyiθxiijmluivi.3.2连续体离散化3.薄板弯曲单元和薄板单元A.薄板单元三角形薄板单元三角形单元有三个节点，每个节点有五个自由度，可承受各个方向的载荷和绕水平轴的弯矩。jwiθyiθxiimviui.3.2连续体离散化4.多面体单元四面体单元四面体单元有四个节点，每个节点有三个自由度。wiiuivijml.3.2连续体离散化六面体单元六面体单元有八个节点，每个节点有三个自由度。wiiuivimlrpqsj4.多面体单元.3.2连续体离散化5.等参单元求解实际问题时人们总希望用最少的单元实现比较高的计算精度，而且所选用的单元对复杂结构也有比较好的适应性。计算期望：手段：单元形状的变化和单元内位移函数的变化用相同数目的结点参数和相同的插值函数进行变换。四边形四节点等参元ηη=1lξξ=1jη=-1iξ=-1mxyo1ijlmξηξ=1ξ=-1η=1η=-1o.4141),(),(iiiiiiyNyxNx3.2连续体离散化四边形四节点等参元插值函数注意:总体坐标系适用于整体结构，局部坐标系只适用于具体某个单元。常用的对于平面问题还有八节点等参元，空间问题有八节点空间等参元，二十节点等参元等。5.等参单元.3.2连续体离散化5.轴对称单元对于回转结构，如果约束条件和载荷都对称于回转轴，其应力、应变和位移也都对称于回转轴线，这类应力应变问题称为轴对称问题，通常用柱坐标来描述应力、应变和位移，单元为实心圆环体，仅截面不同Oyzx三角形环单元.3.2连续体离散化四边形环单元OyzxOyzx回转圆锥薄壳单元5.轴对称单元.3.3单元分析1.单元的插值函数（各种多项式）yaxaavyaxaauee6543212121121098726524321yaxyaxayaxaavyaxyaxayaxaaueejmilrsjmi.3.3单元分析1.单元的插值函数（各种多项式）xyayaxaavxyayaxaauee119874321四节点矩形单元的插值多项式mmlljjiiemmlljjiievNvNvNvNvuNuNuNuNu)by)(ax(N)by)(ax(N)by)(ax(N)by)(ax(Nmlji1141114111411141.axby3.3单元分析1.单元的插值函数（各种多项式）令：)1)(1(41),(iiiNmlji,,,矩阵形式：eemmlljjiimljimljieevuvuvuvuNNNNNNNNvuqN00000000.3.3单元分析1.单元的插值函数（各种多项式）形函数特点：A.本点为1，他点为0),,,(),(mljiijijyxNiji当当01B.在单元内任一点各形函数之和等于11mljiNNNNC.单元任意一条边上的形函数，仅与该边两端节点的坐标有关，而与其他节点无关.3.3单元分析1.单元的插值函数（各种多项式）单元形函数必须满足的条件：A.位移模式在单元内连续，在单元的公共边界处协调.B.位移模式必须包括单元的刚体位移.C.位移模式还必须包括单元的常应变状态..3.3单元分析2.单元分析弹性力学平面问题的几何方程：eemljixyyxBqqBBBBεm,l,j,ixNyNyNxNNNxyyxiiiiiii000000BB称为应变矩阵，其分块子矩阵为：.3.3单元分析2.单元分析代入无量刚插值函数，应变矩阵为：i0i0m,l,j,i)(b)(a)(a)(babiiiii000011100141B令：由平面问题物理方程单元内任意一点的应力可表示为eeexyyxeSqDBqDεσS为应力矩阵.3.3单元分析2.单元分析D为弹性矩阵其表达式为：2110010112EDE是杨氏模量，是泊松比应力矩阵S的分块子矩阵为xNyNyNxNyNxNEiiiiiii212112S)(m,l,j,i.3.3单元分析2.单元分析)1(21)1(21)1()1()1()1()1(40000002iiiiiiiaaabababES代入无量刚插值函数，应力矩阵为：)(m,l,j,i对于平面应变问题：E21E1.3.3单元分析2.单元分析由虚功原理得：eATedxdytδDBBFATedxdytDBBKeeeδKF单元刚度矩阵可分块表示为：emmemlemjemielmelleljeliejmejlejjejieimeileijeiieKKKKKKKKKKKKKKKKK对于平面应力问题每一个子快为：222112112111114rprprprppTrerpKKKK)(EtddtabSBK.3.3单元分析2.单元分析)(ba)(abKprprprprrp3112131111其中：prprrpK2112)(ba)(abKprprprprrp3112131122prprrpK2121r和p遍历i、j、l和m得到单元刚度矩阵对于平面应变问题：E21E1.3.3单元分析2.单元分析单元刚度矩阵的性质：(1)单元刚度矩阵与所选单元的位移模式，几何形状、大小及单元的材料性质有关(2)单元刚度矩阵具有对称性(3)单元刚度矩阵是主元恒为正的奇异矩阵，即单元矩阵没有逆矩阵且0rrK.3.3单元分析3.载荷移置(1)集中力TememelelejejeieieQYXYXYXYXF由虚功原理得到QNFTeQ或yiixiiQNYQNX)(m,l,j,i设在矩形单元体上：体力T),(YXVF分布面力集中力T),(YXAFT),(yxQQQ.3.3单元分析3.载荷移置(3)分布面力由虚功原理得到或)(m,l,j,iTememelelejejeieieAYXYXYXYXFsAeAtdsFNFTsiisiitdsYNYtdsXNXs为单元上作用有外载荷的边(2)体力TememelelejejeieieVYXYXYXYXF由虚功原理得到或)(m,l,j,idtdVAeVFNFTAiiAiidYtdNYdXtdNX.3.3单元分析4.整体分析A.位移节点位移按总体编码由小到大排列起来得到TTT2T121nnδδδδiiivuqn,,,i21其中：B.节点力节点位移按总体编码由小到大排列起来得到n,,,i21其中：TTT2T121nnFFFFNeeiNeeiiiiYXYX11F.1.3.4整体分析C.刚度矩阵nmljikkkkkkkkkkkkkkkknmljimmmlmjmilmllljlijmjljjjiimilijiinn1122k单个单元的扩充.N个单元进行叠加求和dtdNeNee1T1DBBkKnnnmnlnjninmnmmmlmjmimlnlmllljliljnjmjljjjijinimilijiiinmljiKkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk1111111111111.3.4整体分析.单元分析4.整体分析22rpk是的子矩阵，按下式计算Neerprp1kkn,,,r21()21n,,,pFKδ整体刚度K矩阵具有以下性质：(1)整体刚度矩阵是对称的稀疏矩阵(2)整体刚度矩阵的主对角线元素必然大于零(3)未经约束条件处理的刚度矩阵是奇异矩阵.5.