有限元5-矩阵位移法复习.pptVIP

* §2-1 概 述 矩阵位移法的理论基础是传统的位移法，只是它的表达形式采用矩阵代数，而这种数学算法便于编制计算机程序，实现计算过程的程序化。 一、矩阵位移法的基本思路 矩阵位移法又可以称为杆件结构的有限元法； 矩阵位移法的两个基本步骤是 （1）结构的离散化；（2）单元分析；（3）整体分析， 任务 意义 单元分析 建立杆端力与杆端位移间的刚度方程，形成单元刚度矩阵 用矩阵形式表示杆件的转角位移方程 整体分析 由变形条件和平衡条件建立结点力与结点位移间的刚度方程，形成整体刚度矩阵 用矩阵形式表示位移法基本方程 指杆件除有弯曲变形外，还有轴向变形和剪切变形的单元，杆件两端各有三个位移分量， 这是平面结构杆件单元的一般情况。 二、杆端位移、杆端力的正负号规定 一般单元： 凡是符号上面带了一横杠的就表示是基于局部座标系而言的。 现在讨论单元刚度方程。单元刚度方程是指由单元杆端位移求单元杆端力时的一组方程，可以用“ ”表示，由位移求力称为正问题。 在单元两端加上人为控制的附加约束，使基本杆单元的两端产生任意指定的六个位移，然后根据这六个杆端位移来推导相应的六个杆端力。 我们忽略轴向受力状态和弯曲受力状态之间的相互影响，分别推导轴向变形和弯曲变形的刚度方程。 §2-2 单元刚度矩阵(局部座标系) 进行单元分析，推导单元刚度方程和单元刚度矩阵。 一、一般单元 首先，由两个杆端轴向位移 可推算出相应的杆端轴向力 e e e e e 1 2 其次，由杆端横向位移 可以用角变位移方程推导出相应的杆端 横向力 e e e e e e e 将上面六个方程合并，写成矩阵形式： EA l 6EI l2 6EI l2 EA l 12EI l3 12EI l3 4EI l 2EI l 上面的式子可以用矩阵符号记为 e e e e 这就是局部座标系中的单元刚度方程。 e 可求单元杆端力 e e = (1) (2) (3) (4) (5) (6) (1) (2) (3) (4) (5) (6) 0 0 0 0 0 0 6EI l2 0 6EI l2 0 -EA l -6EI l2 -6EI l2 EA l -12EI l3 12EI l3 2EI l 4EI l 0 0 0 0 0 0 -6EI l2 0 6EI l2 0 只与杆件本身性质有关而与外荷载无关 通过这个式子由单元杆端位移 局部座标系的单元刚度矩阵 二、单元刚度矩阵的性质 （1）单元刚度系数的意义 e —代表单元杆端第j个位移分量等于1时所引起的第i个杆端力分量。 例如 代表单元杆端第2个位移分量 时所引起的第5个杆端力分量 的数值。 （2）单元刚度矩阵 是对称矩阵， e 即 。 （3）一般单元的刚度矩阵 是奇异矩阵； e 从数学上可以证明一般单元的刚度矩阵 e 的行列式 e =0 因此它的逆矩阵不存在 从力学上的理解是，根据单元刚度方程 e e e e e e e 由 有一组力的解答(唯一的)，即正问题。 由 如果 e 不是一组平衡力系则无解；若是一组平衡力系，则解答不是唯一的，即反问题。 三、特殊单元 若单元六个杆端位移中有某一个或几个已知为零，则该单元称为特殊单元，其刚度方程是一般单元刚度方程的特例。 e 以连续梁为例： 1 2 e e e e 1 2 e e e e e e e e e 为了程序的标准化和通用性，不采用特殊单元，只用一般单元，如果结构有特殊单元，可以通过程序由一般单元来形成。 *