c弹性力学hapter6.pptVIP

第六章 弹性力学平面问题的有限元分析 三角形单元中的节点位移如下： 根据克莱姆法则，可求出?1， ?2 ， ?3 同理可得 vi = ?4+?5xi+?6yi vj = ?4+?5xj+?6yj vm = ?4+?5x m+?6ym 解出?4 , ?5 , ?6 ? 写成矩阵形式 单元的应变 由于 单元的应力 根据弹性方程 对于平面应力情况 单元刚度矩阵 单元内有虚位移 对于平面应力问题 整体刚度矩阵 每一个单元的子矩阵如下 形函数的性质 Ni 在其余二节点上的值等于零 把面积的行列式以第一行展开乘第二行的代数余子式 把面积的行列式以第一行展开乘第三行的代数余子式 在单元上任一点的三个形函数之和等于 1 面积坐标 面积坐标与形函数的关系 面积坐标对直角坐标的积分 用虚功原理确定等效节点力 若三角形单元上作用有集中力g、分布力q(力/面积)和体积力p(力/体积)，则根据静力等效原理，节点力所做的虚功等于三种力所做的虚功。 集中力的等效节点力计算 体积力的等效节点力 形成载荷列阵{F} §6-4 边界条件的处理和整体刚度矩阵 的修正，计算实例 根据约束情况若在第一点的水平位移为: u1= β1，在第二节点的水平位移为: u2 = β3，把节点所对应刚度矩阵的行和列第一行和第一列及第三行和第三列, 除主对角元改成1，其余的元素都改成零，同时把左端的{F}载荷列阵中对应的行改为己知位移值β1,β3 ，其余的行都减去节点位移值与原来刚度矩阵该行的相应列元素的乘积。 乘大数的方法 同理可得 对于平面应力问题，三角形单元的刚度矩阵 面积坐标对直角坐标的积分 矩形单元的局部座标 矩形单元的位移函数 把 α5 , α6 , α7 , α8 代入 矩形单元的应变 由于u和v都是的ξ, η 函数，微分后代入{ε}得 同样可计算平面应变下的应力[Si]，只需把上式中的 由于矩形单元有四个节点，单元刚度矩阵有四行四列的 子矩阵 由于[Bi], [Bj]都是 ξ , η 的函数，它不能对 x, y积分， 所以需要换元 由集中力引起的等效节点力 由表面力引起的等效节点力 由体积力引起的等效节点力 最后形成载荷列阵{F} 收敛准则 位移函数的选择 它可以是三角函数，也可以是有理分式，在一般情况下取多项式，多项式中各项系数的个数取决於单元的自由度数，在多项式中由於材料的各向同性，多项式的各项必须使 x和 y 对位移有相同的影响，所以就以杨辉三角形选择多项式。 1 x y x2 xy y2 x3 x2y xy2 y3 1.节点应力的计算 把与节点相关各单元上同点的应力取平均值 2. 边界上应力的处理 边界上节点的应力是由内点的应力推算出来的， 一般用抛物线公式插值而得 x1 , x2 , x3 为节点1, 2, 3的座标，xo表示节点o的座标. 第六章 弹性力学平面问题的有限元分析 J 称为雅可比行列式 换元后积分区间为[-1, 1] 第六章 弹性力学平面问题的有限元分析 矩形单元的等效节点力计算 {F0} 表示作用在各节点上的集中力 第六章 弹性力学平面问题的有限元分析 收敛准则，多项式位移函数幂次的选择 1.位移函数必须包含单元的刚体位移 由于刚体位移不产生应变，则{ε}=0，可把下式改写 若 {ε}=0, 则 α2=α6=α3+α5=0 2. 位移函数必须包含单元的常应变值 εx=α2 ,εy= α6 ,, γxy=α3+ α5 由于 第六章 弹性力学平面问题的有限元分析 3. 位移函数在单元内必须连续，并使相邻单元的位移协调. 第六章 弹性力学平面问题的有限元分析 整体刚度矩阵[K]是奇异阵，必须考虑边界约束条件，排除弹性体的刚体位移。消除了整体刚度矩阵的奇异性之后，才能从方程组 中求解节点位移。 一般情况下，所考虑问题的边界往往已有一定的位移约束条件，排除了刚体运动的可能性。否则，应当适当指定某些节点的位移值，以避免出现刚体运动。在引用这些边界条件以后，待求节点未