矩阵位移法程序设计2.pdf

第 五 讲 第二章第二章 平面刚架静力分析的程序设计平面刚架静力分析的程序设计 §2.1 概 述 在上一章 ，我们详细讲述了矩阵位移法 (直接刚 度法) 分析平面杆系结构的基本原理和计算方法 。 本章着重介绍如何把矩阵分析的过程变成计算机 的可执行程序，实现计算机自动化分析与计算结构力 学的问题 。这一过程就是程序设计 。 传统的程序设计通常分两步进行 ： 第一步: 做程序的框图设计 ，把矩阵位移法的计 算过程用流程框图表示 ； 第二步: 是用计算机语言 (例如 FORTRAN ) 编写 程序 。 我们介绍一种新的程序设计方法 —— PAD 软件 设计方法 ，用 PAD 设计代替传统的程序框图设计 。 所谓 PAD ，是 Problem Analysis Diagram 的缩 写 ，它是用一些框和线去描述所计算的过程。 优 点 与传统的程序框图相比 ，它更能简捷 ，明了地表 现程序的逻辑过程 ，与人们的解题思路相一致 ，便于 编写程序 ，因此更易于初学者掌握 。 如图 2.1 列出了几种常用的 PAD 标识符号 ，关 于 PAD 的进一步介绍请参见附录 Ⅰ。 一般处理框 (框中写出处理名或各种语句) 重复框 (后判断循环 ，框中写出循环的条件) S 选择框 (若 Q 条件为真时 ，执行 S ， 1 1 Q 否则执行 S ；若 S 为空时 ，可以省略 2 2 S 2 该框 ) 接口 子程序调用框 (框中写出子程序名) def. 或 定义 (用于添加或分解 PAD ) 定义框 (框中写出定义名) 图 2.1 程序设计的过程分为 四步完成： 1. 把计算过程模块化 ，给出总体程序结构的 PAD 设 计: 所谓计算过程的模块化，就是把矩阵位移法的计算 公式中的每一个相对独立的计算部分当作为一个模块， 每一模块有其独立的功能，各模块之间又有联系。某 些大的模块又可由几个较小的模块组成。 若将每一模块的任务由一子程序来完成，则通过 各子程序之间的接口将各模块连接起来，就得到总体 程序结构的 PAD 设计(如图2.2 所示)。它表现的是程 序总体的结构关系，与理论计算过程相对应。 模块的划分直接影响到程序的简捷程度。一般地， 应将计算过程中重复执行的部分作为小模块，通过对 小模块(子程序)的多次调用完成重复计算。这样不仅 使得程序简捷易读，而且还提高了计算效率。因此， 在进行模块划分时，要对计算理论和基本公式熟练地 掌握。 2. 主程序的 PAD 设计: 主程序的 PAD 与程序结构的 PAD 不同，它不 是直接与所有模块联系，而是仅与某些较大的模块相 联系。通过这种联系，主程序完成对整个结构的分析 计算任务。 3. 子程序的 PAD 设计: 按各个模块所规定的任务，把相应的计算步骤具 体化的表示出来，就是子程序的 PAD 设计(详见以下 各节)。各个子程序完成所规定的具体任务。 4. 根据主程序和各子程序的 PAD 设计 ，用程序语言 编写计算程序 。 §2.2 平面刚架计算的主要标识符和程序结构 1. 数组和变量标识符的说明 为了便于阅读源程序 ，现将程序中所使用的主要 标识符分为整型类和实型类说明如下 ： (1) 整型变量 NN: 结构中的结点总数 (包括所有支座处的结点) ； NF: 固定支座个数； ND: 非固定支座个数 (包括有支座位移的固定支座) ； NDF: 非固定支座中的约束总个数 ； NE: 结构中 的单元总数 ；