通稿-第11章.docVIP

11.10 平面杆件结构的计算程序与算例 本节介绍按矩阵位移法基本理论编制的平面杆件结构计算程序。目前，已出版的结构力学教材给出的平面杆件结构计算程序几乎都是采用BASIC语言或FORTRAN语言编制的，本教材提供C语言编制的计算程序，主要是考虑到学习C语言读者的需要。 11.10.1 平面杆件结构计算程序的适用范围 该程序适用于等截面直杆组成的平面杆件结构，包括：连续梁、刚架、桁架及组合结构。拱形结构，可以沿拱轴线划分若干曲杆单元，应用该程序以直杆单元代替曲杆单元进行近似计算。该程序考虑了弯曲变形和轴向变形，不考虑剪切变形。结构内部的结点可以是刚结点、铰结点或组合结点。结构外部的支座可以是滚轴支座、铰支座、固定端支座或定向滑动支座。作用在结构上的荷载可以是一般荷载（包括结点荷载和非结点荷载）、支座移动或两种荷载的组合。对于温度改变的计算，读者可以在熟悉该程序的基础上编制相应的程序段完成计算，教学过程中可结合上机实验环节进行。 11.10.2 结点位移分量的编码 该程序要求按结点顺序对每个结点按方向、方向、转动方向的顺序编三个位移总码。桁架结点的角位移编零码，桁架单元只有轴向变形，抗弯刚度赋值为零。遇到铰结点时，要按照多结点处理、多位移编码的原则处理。为减少基本未知量的数目，节约计算机内存，可以按下列规则进行结点处理、结点位移编码： （1）刚结点：单结点处理、编三个位移码。 （2）铰结点：若连接单元均为一般单元，需要多结点处理，多位移编码。线位移编同码、角位移编异码；若连接单元是一般单元和桁架单元，可以单结点处理、编三个位移码。角位移总码表示的是一般单元的角位移，并不是桁架单元的角位移，因为桁架单元的抗弯刚度赋值为零；若连接单元均为桁架单元，单结点处理、编三个位移码，角位移编零码。 （3）组合结点：单结点处理，编三个位移码。角位移总码表示的是一般单元的角位移，并不是桁架单元的角位移，因为桁架单元。 例如，图11-31(a)所示刚架，结点编号以及结点位移的编码如图11-31(b)所示。图11-32(a)所示结构，结点编号以及结点位移的编码如图11-32(b)所示。 11.10.3 结构整体刚度矩阵的等带宽存贮 前面介绍过，结构整体刚度矩阵为对称矩阵且表现为稀疏性与带状性。结构中结点未知位移愈多，稀疏性与带状性就愈明显。为节约计算机的内存，程序中存贮的元素采用等带宽存贮，也叫半带宽存贮，只存贮中上三角部分半带宽范围内的元素，带状区域以外的零元素不再存贮。具体为：将中上半带宽范围内的元素，存贮在中。为的对称矩阵、为的矩阵。为结点未知位移的总个数，称为整体刚度矩阵的最大半带宽。的特点是：中主对角线上的元素存贮在中的第一列上；中第行、第列的元素放到中第行、第列的位置。 程序设计中，的最大半带宽作为原始数据输入给计算机，因此要根据结点位移的总码计算。计算方法是：设用、分别表示单元定位向量中的最大分量、最小非零分量，则当向整体刚度矩阵中累加单元刚度矩阵中的元素时所产生的半带宽可按下式计算： (11-66) 对每一个单元按式(11-66)计算，其中最大的即为整体刚度矩阵的最大半带宽。 11.10.4 计算程序的总框图 计算程序的总设计如框图1所示。总框图由两级子框图组成。根据一级子框图1、2、3、4、5编制计算的主程序。二级子框图01、02、03、04、05分别对应一个子程序。主程序与子程序之间的数据传递通过形式参数来实现。 11.10.5 计算程序的应用 计算程序应用时，需要对计算的题目建立一个输入数据文件。在输入数据文件中按照输入的顺序包括以下信息： （1）结构信息：，，，，， 为单元个数；为结点个数；为结点未知位移个数；为整体刚度矩阵的最大半带宽；为结点荷载的个数；为非结点荷载以及支座移动产生的单元在局部坐标系下已知杆端位移的总 个数。 （2）结点信息：，， 、为结点坐标信息，为结点位移分量编码数组。按照结点号的顺序依次输入结点的五个信息，分别为：坐标，坐标，方向位移编码，方向位移编码，角位移编码。 （3）单元信息：，， 为单元杆端结点编号数组，为单元抗拉压刚度数组，为单元抗弯刚度数组。按照单元号的顺序依次输入单元的四个信息，分别为：单元始端结点号，单元终端结点号，单元抗拉压刚度，单元抗弯刚度。遇到、情况，、应取一般单元的倍，例如，可取、。 （4）结点荷载信息： 为结点荷载数组。每个结点荷载输入