有限元基础讲稿第2讲之3(免费阅读).ppt

? Copyright 2000 PTC. Visit Visit ? Copyright 2000 Parametric Technology Corporation Thanks for your listening * ? Copyright 2000 PTC. Visit 第三节 2.1 有限元基本方法 2.2 平面问题 2.3 轴对称问题 2.4 空间问题 2.5 等参元和板壳元 三、轴对称问题 有 限 元 分 析 轴对称问题的有限元分析，实际上是利用对称性把空间问题转化为平面问题来分析。因而轴对称问题的处理方法与平面问题类似，但比一般平面问题稍显复杂。 一般情况下，采用圆柱坐标系，以对称轴为z轴，任一对称面取为rz平面。轴对称问题中只有径向位移u和轴向位移w，它仅与坐标r、z有关，而与θ无关，因而只需考虑坐标平面rz上的位移。 三、轴对称问题 有 限 元 分 析 右图是采用三角形截面的整圆环单元，它由rz面上的三角形ijm绕对称轴z回转一周而得。 单元的基本未知量仍取为结点位移，用列阵表示为： 仿照平面问题，仍取线性位移模式，即： 三、轴对称问题 有 限 元 分 析 利用结点位移插值，仍可得到与平面问题类似的关系式： 其中形函数Ni为： 式中： Nj，Nm依此类推 三、轴对称问题 有 限 元 分 析 矩阵形式 代 入 单元内的应变 dj、dm依次类推 三、轴对称问题 有 限 元 分 析 简写 ?应变分量中的周向应变 不是常量，它与di、dj、dm中的r和z有关。这也是轴对称问题与平面问题的一个重要的不同点。 三、轴对称问题 有 限 元 分 析 单元的应力分量仍可表示为： 弹性应力应变矩阵 令[S]=[D]×[B] 三、轴对称问题 有 限 元 分 析 其中： 三、轴对称问题 有 限 元 分 析 实际上，应力列阵{σ}中只有 在单元体中是常量，其余都不是常量。为了简化后面的积分计算和消除对称轴上的r=0（作分母）所引起的麻烦，实际应用中把各单元中的r和z近似地当作常量，即用各单元的形心坐标，即： 三、轴对称问题 有 限 元 分 析 代 入 di dj dm 经上述近似处理后，各式中di、dj、dm也变为常量。三角形环单元也可以近似地当作常应变单元来处理，求出的应力应变是单元形心处的应力应变近似值。上述处理方法称为常应变处理。 三、轴对称问题 有 限 元 分 析 轴对称问题中，单元刚度矩阵的构成与平面问题类似，与虚功方程相似的轴对称问题的虚功方程为： 单元等效结点力{R}e的虚功 整个三角形环单元中应力的虚功 ?与平面问题不同，这里的结点力是指整个结圆上的力。 三、轴对称问题 有 限 元 分 析 用与平面问题相同的方法推得： 单 元 刚度矩阵 可见，单元刚度方程也得到与平面问题相同的形式： 式中 ： 三、轴对称问题 有 限 元 分 析 单元刚度矩阵 矩阵子块为： （r=i,j,m; s=i,j,m） 由于进行了近似常应变处理，积分号内的变量r和z都是用形心坐标 、 代替，于是： 三、轴对称问题 有 限 元 分 析 实际计算时有： r,s=i,j,m 在等效结点力的移置计算中，一点重要的区别就是，轴对称问题中的所谓结点力实际上是作用在整个结圆上的力。 三、轴对称问题 有 限 元 分 析 值得说明的是，在轴对称问题中进行的近似常应变处理得到的刚度矩阵也是近似的。在要求计算结果非常精确的情况下，常采用增加修正项的办法把近似的刚度矩阵补足为精确的刚度矩阵，需要增加较多的积分计算。 根据实际应用的情况看，在工程误差要求（5％）的范围内，采用近似的刚度矩阵是可行的。当然，误差的大小还与单元的大小有关，因而在使用近似的刚度矩阵时可划分小些的单元来缩小计算结果的最终误差，使之更逼近真解。 ? Copyright 2000 PTC. Visit 下次内容 2.1 有限元基本方法 2.2 平面问题 2.3 轴对称问题 2.4 空间问题 2.5 等参元和板壳元