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

图示结构受到三个约束：节点1在铅直方向，节点2在水平和铅直方向都受到约束，即v1=u2=v2=0 ，这样就可以消除其刚体位移。 R4y 5 1 2 3 4 (1) (2) (3) (4) R3x R5x R5y 支撑条件 前例总刚是10×10阶，总刚方程是10阶线性方程组，由于v1=u2=v2=0，故可把v1、u2、v2所在的第2、3、4行和列统统划去。 （3）处理方法 如用“×”表示非零元素，并写出对称矩阵的下三角阵，则处理前的总刚方程为 10×10阶 处理后的总刚方程为 边界条件处理后，总刚方程变为7阶线性方程组，总刚变成非奇异矩阵。 7个线性方程中含7个未知位移量，故可求得唯一解。 7×7 阶 总刚处理后有二个变化： 一是总刚降阶，如整个结构原有n个节点和r个约束，则总刚在处理前为2n×2n阶，而处理后变为(2n-r) ×(2n-r)阶。 二是总刚变为非奇异矩阵。 边界约束处理后的总刚，其对称性和稀疏性并无影响。 求解(2n-r)阶的总刚方程即可得(2n-r)个节点位移分量，代入式{σ}=[S]{δ}e就可求得各单元的应力分量。 7.平面应力问题有限元法计算步骤 根据具体的结构及给定条件，给出计算简图，标明尺寸，外载及支承条件。 计算单刚。可按公式计算单元面积，应变矩阵[B]，弹性矩阵[D]，最后形成单刚。 选定坐标系，划分单元，编出单元号与节点号(注意节点编号对总刚带宽的影响)，准备各节点坐标值，材料的弹性模量、泊松比，载荷值及其作用的节点号与作用方向，约束的节点及约束方向或已知位移值。 将单刚组集成总刚［K］。 将载荷移置到节点上，并形成载荷向量｛R｝。 边界条件处理。 求解线性方程组以得到节点的位移分量，并可从而求得单元应力 8.3 轴对称问题 若弹性体的几何形状、物理性质、约束条件及外载荷都是轴对称的，则物体的应力、应变、位移等就是轴对称的。 在轴对称问题中，过对称轴任意面的应力、应变和位移是相同的，因而这类空间问题可以简化为类似平面问题处理。 轴对称问题有限元法与平面问题有限元法基本类似，但在数学上要繁琐一些。 1. 应力应变关系 轴对称问题采用柱坐标(r,θ, z)表示一点的位移。 取z轴为对称轴，这时，应变和位移只与z和r有关，与θ无关，即 轴对称体的三角形单元 u=u(r,z) w=w(r,z) 轴对称问题需考虑四个应力分量 对应的应变分量为 几何方程 物理方程 弹性矩阵 2. 三角形截面环单元 轴对称问题中，采用三角形截面的圆环单元，它是由roz面上的三角形ijm环绕对称轴z旋转一周得到的。 相邻单元在其棱边互相连接，单元的棱边都是圆，称为节圆。每个节圆与roz平面的交点就是节点。 各单元在roz平面上形成三角形网格。基本未知量仍取节点位移。 单元节点位移可表示为 仿照平面问题，取线性位移函数 求得单元位移 将单元位移代入几何方程，得单元体内的应变 单元的εr、εz、γrz应变分量都是常量；但环向正应变εθ 不是常量，它与 fi、fj、fm中的r有关。 Δ为三角形单元面积 把单元三节点i, j, m的坐标及位移分别代入式(8-28) ，得 i j m a,b,c与节点坐标有关的常数 把系数αi代入式(8-28) ，并按ui、uj、um和vi、vj、vm集项，得 写成矩阵形式 式中 Ni—三角形单元形函数，反映单元的位移形态 i j m 单元形函数矩阵 (2) 单元刚度矩阵 根据设定的单元位移函数，依次利用几何方程、物理方程与虚功方程，即可导出节点位移与节点力间的关系式，亦即单元刚度方程，从而建立单元刚度矩阵。 几何方程 将位移函数代入并化简 节点位移→应变 i j m [B]—应变矩阵，各元素只是与单元节点坐标有关的常量，因此单元内任意点的应变分量为常量，故称为常应变单元。 简记为 应变矩阵 式中 弹性矩阵[D]和应变矩阵[B]中各元素都是常量，故单元中的应力分量也是常量。 不同的单元往往具有不同的应力。因此在相邻单元的公共边界上，将存在应力突变。 物理方程 节点位移→应力 或写成 式中 应力矩阵 i j m 外力(节点力)所做虚功等于应力所做虚功 节点位移→节点力（推导单刚） 设厚度为t的三角形单元发生某种虚位移，节点虚位移为{δ*}e，引起的虚应变为{ε*}e，则虚功方程为 将 {ε*}=[B]{δ*}e 和{ε*}T = ({δ*}e)T[B]T代入上式，得 σx τxy σy i j m Ui Uj Um Vi Vj Vm σy σx 三角形单元受力状态 （外力与内应力） 单元的刚度方程 可写成 单元刚度矩阵(简称单刚) 式中 或 [D]和[B]都是常量 [K]e表示单元的刚度，与弹簧刚度的物理本质相同，即抵抗变形的能力。 [K]e为6×6阶方阵，共