有限元分析(基础).pptVIP

矩阵代数复习 1、矩阵定义 一组元素按行、列次序排列成的矩形阵列称为矩阵。若矩阵的元素排列为m 行和n列，称为mn 阶矩阵。 2、方阵 3、行矩阵和列矩阵 一个单独的行组成的矩阵称为行矩阵，如： 由单列组成的矩阵称为列矩阵，如： 1 4、纯量(标量) 仅由一个单独的元素所组成的11阶矩阵称为纯量。 5、矩阵乘法 两个规则： （1）两个矩阵仅当他们是共形时才能相乘，即 （2）不具有交换律，即 2 6、转置矩阵 将一个阶矩阵的行和列依次互换，所得的阶矩阵称之为原矩阵的转置矩阵，如： 其转置矩阵为 当连乘矩阵的乘积被转置时，等于倒转了顺序的各矩阵的转置矩阵之乘积。若 A=B C D 7、零矩阵 元素全部为零的矩阵称为零矩阵，用0表示。 3 任意矩阵与单位矩阵相乘仍等于原矩阵，即 4 10、逆矩阵 在矩阵运算中，没有矩阵的直接除法， 除法运算由矩阵求逆来完成。例如，若 一个矩阵的逆矩阵由以下关系式定义： 矩阵求逆时必须满足两个条件： （1）矩阵是一个方阵。 （2）矩阵的行列式不为零，即矩阵是非奇异矩阵（行列式为零的矩阵称为奇异矩阵）。 正交矩阵的逆矩阵等于其转置矩阵，即 5 1 结构离散与向量表示 工程上许多由金属构件所组成的结构，如塔式桁构支承架、起重机起重臂架、钢结构桥梁、钢结构建筑等可以归结为杆系结构。杆系结构按各杆轴线及外力作用线在空间的位置分为平面杆系和空间杆系结构。 杆系结构可以由杆单元、梁单元组成。 (a) 塔式起重机 (b) 履带式起重机 (c) 钢结构桥梁 (d) 埃菲尔铁塔 图1 杆系结构 6 1.1 结构离散化 由于杆系结构本身是由真实杆件联接而成，故离散化比较简单，一般将杆件或者杆件的一段( 一根杆又分为几个单元 )作为一个单元，杆件与杆件相连接的交点称为结点。 杆系结构的离散化的要点可参考如下： a. 杆件的转折点、汇交点、自由端、集中载荷作用点、支承点以及沿杆长截面突变处等均可设置成结点。这些结点都是根据结构本身特点来确定的。 b. 结构中两个结点间的每一个等截面直杆可以设置为一个单元。 c. 变截面杆件可分段处理成多个单元，取各段中点处的截面近似作为该单元的截面，各单元仍按等截面杆进行计算。 7 d. 对曲杆组成的结构，可用多段折线代替，每端折线为一个单元。如若提高计算精度，也可以在杆件中间增加结点。 e. 在有限元法计算中，载荷作用到结点上。当结构有非结点载荷作用时，应该按照静力等效的原则将其变换为作用在结点上的等效结点载荷。 (a) 结点载荷处理方式 (b) 等效结点载荷处理方式 图2杆系结构离散化示意图 8 1.2 坐标系 图3 坐标系示意图 为了建立结构的平衡条件，对结构进行整体分析，尚需要建立一个对每个单元都适用的统一坐标系，即结构坐标系或称之为整体坐标系、总体坐标系。 9 1.3 向量表示 在有限单元法中力学向量的规定为：当线位移及相应力与坐标轴方向一致时为正，反之为负；转角位移和力矩，按右手法则定出的矢量方向若与坐标轴正向相一致时为正。对于任意方向的力学向量，应分解为沿坐标轴方向的分量。 刚架结构示意图 (b) 结点位移和结点力分向量 图4 平面刚架分析示意图 10 结点位移列向量为 单元e结点位移列向量为 结点力向量为 单元e结点力列向量为 11 2 位移函数及要求 有限单元法分析中，虽然对不同结构可能会采取不同的单元类型，采用的单元的位移模式不同，但是构建的位移函数的数学模型的性能、能否真实反映真实结构的位移分布规律等，直接影响计算结果的真实性、计算精度及解的收敛性。 为了保证解的收敛性，选用的位移函数应当满足下列要求: a. 单元位移函数的项数，至少应等于单元的自由度数。它的阶数至少包含常数项和一次项。至于高次项要选取多少项，则应视单元的类型而定。 b. 单元的刚体位移状态和应变状态应当全部包含在位移函数中。 c. 单元的位移函数应保证在单元内连续，以及相邻单元之间的位移协调性。 12 13 a. 把连续体离