第二章 弹性力学平面问题有限元法精品.ppt

（2）递推公式： 上面的消元反代作法同样适用于系数矩阵为n阶的线性代数方程组，由上过程可知，消第一元时，系数矩阵从f1变到f2，也就是说，每消元一次，要作四件事（建立四个递推公式）。 （1）修改非主导方程的系数项（2）修改主导方程的系数项（3）修改非主导方程的右端项（4）修改主导方程的右端项 其中,相应于(2)(3)(4)项的递推公式容易建立,唯有非主导方程系数项的处理涉及到有关元素如何从刚度矩阵[K]中过渡到[SK]中的问题。 递推公式： 值得注意的是： 1.每消元一次，需修改的方程只涉及到NBD个； 主导方程——1 非主导方程——2，3，4，...NBD 2.随着被修改方程的行数增加，方程中需修正的系数逐个减少（呈三角形，见下图）。3.在总刚[K]中的对角元，在[SK]中变为第一列。 注意: （1）最后一个右端项就是最后一个未知量的解。 （2）反代从第NEQ-1个方程开始，倒推。反代结束。{R}阵是节点位移。 （3）每反代一个方程涉及该方程的NBD-1列元素和NBD个右端项。 2.5 应力计算 全部位移求得后，各单元的节点位移确定后，单元内的应力可由公式（2-14）求得： 其中应力矩阵[S]=[D][B] 在计算程序中,计算某一单元的应力时，首先要确定该单元的三个节点的六个位移分量是已求得的全部节点位移中的那几个。相对应的编号关系建立后，就可以把这些位移分量代入应力分量关系式（2-14）求得σx，σy，τxy。 如仍以图（2-16）为例，当消元的反代过程结束后，节点位移分量全部求得： 如欲确定单元（4）的应力，按（2-14）式: 需从全部位移分量中找出单元（4）相应的位移分量。 应力分量求得后，主应力和主方向按下式计算： 2.6 平面问题有限单元法小结 现在，我们来回顾一下用有限单元法按位移求解弹性平面问题的步骤和方法： （1）根据工程实际问题对计算精度的要求和计算机的容量，将平面弹性体剖分为一定数量的三角形单元，并对单元和节点进行编号，选定右手坐标系XOY（通常取水平X轴向右为正，垂直y轴向上为正），定出所有节点的坐标值： 然后对每个单元按逆时针顺序，确定i,j,m的具体节点编号，称为单元节点信息。 （2）根据节点坐标值，用式（2-2），（2-3）标出各个单元的几何参数bi, ci , bj , cj , bm , cm 和三角形面积Δ值。 （3）再由这些数值，连同各单元的材料常数E和μ，按式（2-18b）或（2-20a）（平面应力问题）或式（2-20b）（平面应变问题）计算每个单元的刚度矩阵[K]e中的各元素。 如为平面应变问题，将式中的E换成 E/(1-μ2)和μ换成μ/(1-μ)。于是得到： (4)按 形成弹性体的总刚度矩阵[K]，即将每个单元刚度矩阵中的元素按节点编号的次序，填到总刚度矩阵相应的行、列位置上去。如果在某个位置上填入若干个数，那就把它们迭加起来。 （5）进行载荷移置，算出弹性体总的节点载荷列阵{R}，从而得到以节点位移{δ}为未知数的基本方程式（2-23a）。 （6）根据弹性体的几何边界条件，对上述基本方程进行处理，得到式（2-25）（2-26）。 （7）解上述线性代数方程组，求得全部节点位移。 （8）用式（2-14）求出各单元的应力分量 σx，σy，τxy ，这些应力在一个单元内是常量，通常把它们当作是作用在单元的形心上。另一种办法是求出节点处的应力；这时我们可认为节点处的应力等于环绕这个节点的若干个单元应力的平均值。 （9）根据单元形心或节点处的应力分量，即可求得其主应力及主方向。 在弹性体的不受外载荷作用的自由边界上，一个主应力的方向应与边界曲线相切，而另一主应力的值应为零。 当弹性体在边界上受到法向分布载荷作用时，一个主应力方向同样应与边界曲线相切，而沿法向的主应力应等于法向载荷的集度。 利用上述方法，也可以对用有限单元法求得的解答的精确度进行估计。 有几个约束条件，就应重复几次上面的修改过程。 处理边界节点零位移条件，还可以采用其它方法，这里不一一列举。 （2）位移不为零的约束处理 如果某些节点的位移值是已知的，如： 我们可将基本方程作如下修改：将第i行的对角元改为1 第i行的其它元素改为零 将右端项改为 式（2-27a） i行 i列 这样修改后，第i个方程变为： 为了保持总刚度矩阵的对称性，以便于以后求解，我们把[K]的第i列也改为零，这时，节点的载荷列阵要作如下修改： i行 i列 式（2-27b） 当将（2-27a）式和（2-27b）式乘开时，两