平面刚架程序的框图设计和源程序.pdf

附录 A-2 平面刚架程序的框图设计和源程序 附录 A-2 平面刚架程序的 框图设计和源程序 本附录共有三部分内容。第一部分为平面刚架程序的框图设计，此项 框图设计主要是针对 FORTRAN77 语言而编写的，对于 Fortran90 语言亦 可以作为参考。第二部分为平面刚架源程序（FORTRAN77 语言）。第三部 分为平面刚架源程序（Fortran90 语言）。 A-2-1 平面刚架程序的框图设计 一、总框图与程序标识符 (一)总框图( 框图 1) 框图 1 总框图 平面刚架的整个计算过程参见上面的总框图。 总框图由两级子框图组成。平面刚架的整个计算过程由一级子框图 1、 2 、3 、4 和 5 组成。根据这 5 个子框图可以编写出平面刚架的主程序。二级 - 1 - 附录 A-2 平面刚架程序的框图设计和源程序 子框图 01、02、03、04 和 05 分别对应一个子程序，总框图中亦标明了子程 序的名称，它们在运算过程中被主程序多次调用。主程序与子程序之间的 数据传递均通过形式参数来实现。 按照上述总框图编制的平面刚架静力计算程序的适用范围是： 1．结构形式 由等截面直杆组成的具有任意几何形状的平面杆系结构：刚架、组合 结构、桁架、排架和连续梁。 平面刚架单元之间的连接结点可以是刚结点、铰结点或刚铰混合结点。 2 ．支座形式。 结构的支座可以是固定支座、铰支座、滚轴支座和滑动支座。 3 ．荷载类型。 作用在结构上的荷载包括结点荷载与非结点荷载，各种非结点荷载类 型见《结构力学 I — 基本教程》的表 9 －1。 4 ．材料性质。 结构的各个杆件可以用不同的弹性材料组成。 在平面刚架的矩阵分析中，考虑了杆件的弯曲变形和轴向变形，而忽 略了剪切变形的影响。 (二)程序标识符 现将子框图和源程序中主要标识符的意义说明如下： 1．整型变量 NE——单元数； NJ——结点数； N ——结点位移未知量总数； NW——最大半带宽； NPJ——结点荷载数； NPF——非结点荷载数； IND——非结点荷载类型码； M——单元序号。 2 ．双精度型变量 BL——单元长度； SI——单元的 sin 值； CO——单元的 cos 值。 3 ．整型数组 JE (2,NE)——单元杆端结点编号数组； JN (3,NJ)——结点位移分量编号数组； - 2 - 附录 A-2 平面刚架程序的框图设计和源程序 JC (6)——存放单元定位向量的数组。 4 ．实型数组 EA (NE)、EI (NE)——单元的 EA 、EI 数值； X (NJ)、Y (NJ)——结点坐标数组； PJ(2,NPJ)——结点荷载数组； PF (4,NPF)——非结点荷载数组。 5 ．双精度型数组 KD (6,6)——存放局部坐标系中的单刚k e 的数组； KE (6,6)——存放整体坐标系中的单刚k e 的数组； T (6,6)——存放单元坐标转换矩阵 T 的数组； KB (N,NW)——存放整体刚度矩阵 K 的数组； P (N)——结点总荷载数组，后存结点位移； F0 (6)——局部坐标系中的单元固端力数组； F (6)——先存放整体坐标