结构有限元.docVIP

 --确定列高的方法：ID表 --主对角线元素地址 数组[A] --计算主元地址的公式可写为：MAXA j+1 MAXA j + hj式中， j 1，2，…，N；hj——刚度矩阵[K]第j列的列高。一维数组A的总长度 S ，即刚度矩阵K按变带宽存贮的总存贮量：S MAXA N+1 - MAXA 1 ----某一元素在一维数组中的地址由其相应主元地址决定。 --1. 结构刚度矩阵之所以能化为带状矩阵，主要 利用了其 性和 性。 2. 结点约束信息表ID数组的作 用： 。 3. 结构刚度矩阵采用变带宽存储时，一维数组 A存储 ，一维数组MAXA存 储 。 4.在等带宽存贮中，最小带宽UBW由 确定。 作业题：4. 图示结构为一平面桁架结构，按图示方式对结点 编码，若总体刚度矩阵采用等带宽存贮，总自由度 数NDF为 ，最大半带宽UBM为 ；若总 体刚度矩阵采用变带宽存贮，其总体刚度矩阵的第 5个自由度，对应于第 结点的第 个自 由度，第5列的列高为 。 -- 结点荷载总列阵的组集 输入文件中要包含 Pc 的数据信息：读入结点力向量N, F1, F2, F3, F4, F5, F6 --涉及到两个子程序： （1）计算整体坐标系下由跨间荷载引起的等效结点荷载。 （2）计算整体坐标系下由跨间温度引起的等效结点荷载。 ---1. 刚度矩阵的性质 对称矩阵、正定矩阵、 稀疏矩阵。 2. 三种存储方式 ? 方阵存储 ? 等带宽存储 ? 变带宽存储 压缩存储?利用对称性，等带宽存贮和变带宽存贮只存 贮刚度矩阵的上三角部分。即主对角线元素 以及上半带区（或顶线以下）内的元素。 3. 组集结构刚度矩阵工作流程： 单刚 -----?判断约束情况-----?对号入座组集总刚 不去存贮与刚性支承约束相应的行和列 4. 结点编号应得当 无论采用等带宽存贮或变带宽存贮，都要求结 点号码编得恰当，使得刚度矩阵的非零元素都聚集 在主对角线附近狭窄的带形区内。对于大型结构应 该特别注意选取恰当的结点编号方式。这是因总存 贮量大体上与带宽成正比的原故（对等带宽存贮， 则准确地与最大半带宽成正比）。 -----线性方程组的若干直接解法 有限元法分析到最后总是归结为解线性代数方程组：[K] δ P 有限元法的求解效率很大程度上取决于线性代数方程组的解法（包括矩阵的存储方式），因为在有限元法的整个求解时间中解线性代数方程组的时间占了很大比重，当单元增多、网格加密、未知数成倍增加时，尤为如此。 --线性方程组的数值解法一般有两种。直接解法：经过有限步算术运算，可求得方程组的精确解的方法（若在计算过程中没有舍入误差）。迭代解法：用某种极限过程去逐步逼近线性方程组精确解的一种近似计算方法。 迭代法具有占存储单元少，程序设计简单，原始系数矩阵在迭代过程中不变等优点，但存在收敛性及收敛速度等问题。 --本章介绍最基本的两类直接解法： （1）高斯消元法 针对[K]的等带宽存储 （2）LDLT分解法 针对[K]的变带宽存储 ---高斯消元法的基本思想 基本思想：通过逐步消元（行的初等变换），把方程组化为系数矩阵为三角形矩阵的同解方程组，然后用回代法解此三角形方程组（简单形式）得原方程组的解。 两大步： 1. 消元 2. 回代 --用高斯消元法求解方程组 （1）如果[K]矩阵写成方阵的形式，则高斯消元采用方阵存储算法。 （2）如果[K]矩阵写成带状矩阵的形式，则高斯消元法采用等带宽存储算法。 --高斯消元等带宽存储 1 每一次消元运算，只需在带宽范围内进行。 （2）元素aij （i≤ j）在矩形数组中的行号和列号为：ROW i COL j-i+1即行号不变，列号左移i-1 （3）消元公式 （4）待修改的方程号码k的变化范围 （5）方程中需修改的元素的下标变化范围 （6）回代过程 -- 实对称正定矩阵的LDLT分解 结构刚度矩阵为实对称正定矩阵 [A]对称：[A]T [A] ? [A]正定：[A]的各阶顺序主子式均大于零。即 --对称矩阵的乔来斯基（Cholesky）分解 定理： 设[A]为对称正定矩阵，则存在唯一分解[A] [L][D][L]T，其中[L]为单位下三角阵，[D] diag d11,d22,…,dNN 且di 0 i 1,…,N . --变带宽存储下LDLT分解算法 [A]为N阶正定矩阵，假设已按变带宽存储方式进行了存储。即：一维数组A存贮矩阵[A]上三角变带宽内的元素； 一维数组MAXA存储矩阵主对角元在数组A中的地址。 一维数组A的总长度 S ，即刚度矩阵[K]按变带宽存贮的总存贮量：S MAXA N+1 - MAXA 1 ---编程思路 1、程序的功能 ?把一维数组存储的N