安徽建筑工程学计算结构力学6.pptVIP

第六章 线性代数方程组的求解 　　因为消元分解法是按行进行消元的，故一般采用一维变带宽上三角按行存贮(LDLT法采用)。 　　这时只需141+30=171个存贮单元，对于大型复杂结构来说，节约的存贮空间则非常可观。 　　KD数组的程序设计见讲义所附结构静动力综合分析程序PSTDY中的SUB.QKD(P106)，或参考教材P196。 3. 带状矩阵的解法 (1)形成一维变带宽上三角按行存贮的结构刚度矩阵ZK(NZY)，NZY:存贮的刚度系数的总数。其框图设计同SUB. KJX1，只不过增加: 　　①调用KD数组； 　　②由KD数组寻址。 　　程序见后附结构静动力综合分析程序PSTDY中子程序SUB.KJX(P110)，或参考教材P199。 * 首页 上页 返回 下页 * 计算结构力学 第六章 线性代数方程组的求解 结构刚度方程 　　　　　　[K]{Δ}={P} 　　求解方法有直接法和迭代法两大类。直接法从Gauss消元法衍生而来，在程序设计中又称为消元分解法。 解:由5—3中例可知刚度方程为 Gauss消元法解线性方程组的步骤： 第一步： (令[K](1)=[K]) 　　除K11外，使第1列的元素为零，即(1)式中第二式+第一式 6-1 Gauss消元法 例:消元法求解 结构刚度方程 　　由于 ，故第三行类似的运算就不必进行，式中上标(1)、(2