8有限元法基础(白).pptxVIP

第八章 有限元法基础;弹性力学问题 平面问题 轴对称问题 有限元法基本概念 有限元法解题基本方法、步骤 形函数、应力矩阵、应变矩阵、单元刚度矩阵、整体刚度矩阵等;有限元法发展历史;;8.1 有限元法的简单引例(一维有限元法) ;杆单元;对于每一个单元，可以用线性函数去近似地描述实际的单元位移。设单元(e)ij位移函数为：;;记{f}=ue，则;根据虚功原理，移置前单元载荷所作虚功应等于移置后节点载荷所作虚功，所以有;固定端的节点1还受约束反力R=-qL，于是单元载荷移置后，各节点的节点载荷R1、R2、R3、R4分别为;把 {f}=[N] {δ}e 代入几何方程ε= d u / dx，得 ;把{ε}=[B]{δ}e 代入物理方程σ=Eε，得;利用虚功方程分析单元的节点力与节点位移的关系。;;(r, s = i, j。r＝s 时，取“＋”；r≠s 时，取“－”);;对每个节点列出力平衡方程，得;把各单元刚度矩阵补零升阶到4×4阶，则;;把[K]和{R}具体表达式代入 [K]{δ}={R} ，得;7.求解方程组;有限元法解题步骤;总体结构合成 通过节点平衡方程建立以节点位移为未知量、以总体刚度矩阵为系数的线性代数方程组Kδ=R。;8.2 平面问题有限元法 ;(2) 平面应变问题;应力列阵 {σ}={σx σy τxy}T;2. 连续体的离散化 ;平面问题常用单元的类型 ;3. 单元分析;位移法中，以六个节点位移分量作为基本未知量。对应的物理量是六个节点力分量.;单元分析的任务，就是推导基本未知量{δ}e与其对应量{F}e之间的关系，即;(1) 单元位移函数 ;Δ为三角形单元面积;;形函数的几何解释;;i j m;;外力(节点力)所做虚功等于应力所做虚功;单元的刚度方程;[K]e可写成分块形式;对称性。单刚是对称矩阵，即[K]e=([K]e)T，或 Kij=Kji.;4. 载荷移置 ;设单元ijm中任意一点M(x, y)受集中力P，其分量为Px和Py，即{P}={Px Py}T。;体力的移置 ;面力的移置 ;单元上同时受到集中力、体力和面力时，可按上述方法分别移置，然后迭加。;整体分析;5. 总体刚度矩阵 ;设离散后的单元组合体有五个节点四个单元，节点与单元的编号及总体节点载荷如图所示。;建立节点平衡方程时要解除节点约束，以约束反力代之。节点约束反力与节点载荷具有同样性质，可广义地称为节点载荷。;节点的平衡方程 ;用分块矩阵表示 ;(C);将式(B)~(E)代入式(F) ，得;(2) 总体刚度矩阵的形成方法;1) 主对角线子块Krr （r=s）由环绕节点r的诸单元刚度矩阵相应对角子矩阵[Krr]e迭加而成。;2) r≠s时，若r、s是相邻单元公共边上的两个节点，则总刚子矩阵[Krs]就是公用该边的两个相邻单元的单刚子矩阵[Krs]e的迭加。;3) 当r≠s时，若r、s仅是某单元一个边上的两节点，则它们等于该单元刚度矩阵相应的一个子块。;4) 当r≠s时，若r、s不同属于一个单元的两个节点，即互不相关的两个节点，则总刚的这些子矩阵[Krs]=0。;总刚子矩阵[Krs]物理意义(与单刚子矩阵一样) r—节点力作用的节点编号， s—产生位移的节点编号。 当r=s时，该节点的位移与所有公用该节点的单元在该节点的节点力有关。;1) 对称性。总刚为对称矩阵，所有对称于主对角线的元素都相等。利用其对称性可在计算机中只存贮矩阵的上三角或下三角部分。;;;6. 边界条件的处理 ;;;处理后的总刚方程为 ;;7.平面应力问题有限元法计算步骤 ;将单刚组集成总刚［K］。 将载荷移置到节点上，并形成载荷向量｛R｝。 边界条件处理。 求解线性方程组以得到节点的位移分量，并可从而求得单元应力;8.3 轴对称问题 ;1. 应力应变关系 ;;2. 三角形截面环单元 ;单元节???位移可表示为;;;;3. 单元刚度矩阵 ;矩阵［B］中含有坐标r、z，使积分运算比较复杂。为简化计算取单元形心的坐标 、 代替矩阵[B]中的坐标r、z，得到一个近似单元刚度矩阵;4. 载荷向节点移置 ;;设rz平面上单元ijm的ij边有线性分布的径向表面力，在节点i, j的集度分别为 ， ，ij边的边长l。此时有， ， ，于是节点i的等效节点载荷为 ;;采用位移法有限元分析弹性问题的步骤是什么？ 什么是形函数、应力矩阵、应变矩阵、单元刚度矩阵，形函数有什么性质，说明单元刚度矩阵的特点和物理意义。 单元刚度矩阵如何叠加成整体刚度矩阵？为什么应该这样叠加? 叠加成的整体刚度矩阵又有什么特点? 如图(a)所示为一厚度t的等厚矩形薄板，一端固定，一端承受均匀拉力q。该