虚拟土木工程程序设计大赛.docVIP

“虚拟土木工程”程序设计大赛 ―――设计文档及程序说明 项目名称: VB平面珩架静力分析器 项目内容： 该程序是按照矩阵位移法的后处理法的基本原理和分析过程，并借助 Visual Basic 所设计的，适用于平面珩架在节点静力荷载作用下的内力与位移计算问题。 主要功能包括：节点位移，杆件轴力的计算；珩架节点位移的动态演示； work plane上的增添、修改、查询等可视化操作；输入文档与输出文档 项目意义: 一、 自主研发的结构力学计算工具，培养团队意识 二、 在矩阵位移法的基本原理上加入便捷的可视化操作，以及前处理，后处理 三、 针对教学辅助――帮助理解矩阵位移法以及平面珩架受力、变形形态 四、 建立输入文档、输出文档，方便储存以及读入 五、 程序模块具有通用性，子程序可直接用于平面刚架、空间杆系结构分析。 程序设计： 由于先处理法在计算和形成总刚度矩阵时先已将支座约束条件作了考虑，单位刚度矩阵的元素在送入总刚度矩阵之前，某些位移为零的行与列被删除，故总刚度矩阵的阶数较之后处理法低，占用的计算机存储量也低，相应地，线性方程的求解速度相对更快。 一、INPUT： 子程序INPUT共包括五个输入语句，先后输入五组原始数据。 依照用户在对话窗口中输入的数据给这五组数据赋值，并同时生成输入文档。 也支持用户在文档中输入后直接读取该文档的快捷操作 ( 1 ) 基本参数 对应输入语句 依次填入珩架的节点总数，单位总数，受外荷载作用的节点总数，支座节点总数，材料的弹性模量 （2）节点坐标 对应输入语句 以节点为序依次填入节点的编号和节点的x，y坐标值 （3）单元两端节点号和单元横截面面积 对应输入语句 以单元为序依次填写单元的编号，单元两端节点的编号和横截面面积 （4）节点荷载 对应输入语句 对有外荷载作用的节点，以节点为序填入节点的编号和x，y方向的节点荷载（方向定义为与结构坐标系的方向一致为正），无荷载时，数值填入零 （5）支座约束信息 对应语句 对支座节点，以节点为序填入支座节点的编号，节点x，y方向位移状态的指示信息（方法为：1为自由端，0为约束端。例如则IC（1）=1 ; IC（2）=0） 二、ASSEM： 该子程序功能为在TK数组中生成总刚度矩阵。先后调用子程序STIFF和 ELASS 总刚度矩阵采用二维等带宽存放，以节省存储单元，提高计算效率 上图中，平行线之间宽度方向可含的元素个数称为总刚度矩阵的带宽，为M 由总刚度矩阵主对角线元素开始至其中一条平行线之间宽度可含元素的个数为半带宽，为MS 当采用高斯消去法求解线性方程组时，半带宽M以外的零元素可不必计算，又因为是对称矩阵，故实际只要计算半带宽MS范围内。 总刚 单刚 三、STIFF 该子程序为按照下式计算出结构坐标系中珩架单元的刚度矩阵，并存放在二维数组ELST中 [k]= 四、ELASS 该程序用于将单元矩阵送入总矩阵 单元刚度矩阵的元素是以节点子块为单位送入总刚度矩阵 当前子块在单元刚度矩阵中的为位置由变量 的值确定，它们分别为单元刚度矩阵中位于该子块之前的行数和列数 子块送入总刚度矩阵后的位置由变量 的值确定，它们分别为总刚度矩阵中位于该子块之前的行数和列数，如下图所示： 五、BOUND 该子程序功能为在总刚度方程中引入位移边界条件 若IB对应元素的值为1则表示该项位移并未受到约束，不需作对应的边界条件处理 若IB对应元素的值为0则表示该项位移受到约束，需对总刚度方程作对应的边界条件处理 六、SLBSI 该子程序功能为求解线性方程组 方程组的求解按照高斯消去法的基本原理进行。并且由于采用二维等带宽存放，避免了不必要的计算，从而提高了计算的效率 七、FORCE 该子程序是根据已求得的节点位移AL，计算杆件的轴力和节点合力，分别存放于数组FORCE和REAC中。对于支座节点而言，节点合力即为支座反力 八、OUTPUT 该子程序输出计算结果（包括节点位移、杆件轴力、节点合力），分别存放在AL, FORCE, REAC中便于后处理时调用。同时生成输出文档，以后直接读取该文档即可恢复数据以及work plane中的工作状态。 主程序 INPUT输入 ASSEM形成总刚 BOUND边界条件 STIFF 与 ELASS 7 SLBSI解