平面刚架计算(APF).docVIP

平面刚架内力计算程序 APF 一、程序功能及编制方法 平面刚架内力计算程序（APF），能够计算任意形状的平面刚架和连续梁的内力和结点位 移。能够适应任意荷载、温度变化及支座移动等复杂情况。但不能处理刚架有内部铰结点的情况。 程序用直接刚度法集成总刚度矩阵，等带宽二维数组存储总刚度矩阵，后处理法引进支座条件，高斯消去法解刚度方程。 本程序图形显示功能只显示刚架的计算简图, 以作为判断输入数据是否正确的参考。 二、程序使用方法 1．填写数据文件 将要计算的题目按规定的格式用Word或Windows中的” 记事本”或”写字板” 编辑写成数据文件（数据文件的填写格式见后），并给数据文件一个恰当的名字后存入程序包所在的文件夹中。 2．运行 按程序包提示的步骤运行。计结果全部数据部将写在这个文件中。计算结束后可用Windows中的“记事本“打开结果文件的内容，或用打印机打印输出文件以供使用。 三、数据文件填写格式 可用Word或Windows中的“记事本”编辑器建立数据文件。 数据文件的内容为： 总信息：NJ，NE，NZ，NB，NP，NF，E0，DS 其中：NJ——结点数。通常把杆件的交点作为结点，也可将单元上任意点当成结点； NE——单元数； NZ——支座约束数; NB——杆的类型数。不同截面尺寸的杆就属于不同的类型； NP——结点荷载数； NF——非结点荷载数； EO——弹性模量； DS——是否显示热打印结点位移的信息。如要求显示和打印结点位移，则DS 1否则 DS＝0。 2．结点坐标数组XY NJ，2） 每个结点填一行，每行两个数分别是结点的x坐标和y坐标。 3．单元信息数组JM NE，3） 每个单元填一行，每行三个数分别是单元左右两端的结点号和单元的类型号。 4．单元截面信息数组AI NB，2） 每类单元填一行，每行两个数分别是单元横截面面积和惯性矩。 5．约束信息数组ZC NZ，2） 每个支座约束填一行，每行2个数分别为约束所在的结点号和约束方向。约束方向规定沿X方向的约束为“1”，y方向的约束为“2”，力矩约束为“3”。 6．结点荷载信息PJ NP，3） 每个结点荷载填一行，每行3个数分别为荷载作用结点号，荷载的方向和荷载的数值。荷段作用方向以沿x方向为“1”，y方向为“2”力矩为“3”。荷载数值以与整体坐标方向一致的力为正值，反之为负值。以 绕结点顺时针转的力矩为正值，反之为负值。 如果没有结点荷载，该数组不填。 7．非结点荷载信息PF NF，5） 每个非结点荷载填一行，每行5个数分别为荷载作用的单元号、荷载类型号和常数c、a、G。程序中提供了9种不同类型的荷载，任何复杂的荷载都可由这9种荷载组合而成。非结点荷载的正负号与单元坐标系有关，与单元坐标系方向一致的荷载为正值，反之为负值。各类型荷载如下图所示。 如果结构上未作用非结点荷载则此数组不填。 四、计算例题 例APFEX1 刚架如图所示，已知其横梁尺寸为： 40cm ×80cm（A1 0.32 m2，I1 0.01707m4）。柱的截面尺寸为40cm0cm（A2 0.24m2，I2 0.0072m4），弹性模量 E0 1E7kN/m2。试求图示荷载作用下的各柱端内力。 b）所示，然后建立名为APFEX1.DAT的数据文件，内容为： 5，4，7，2，1，3，1E7，1 0，0，0，4，6，0，6，4，12，4 1，2，2，2，4，1，3，4，2，4，5，1 0.32，0.01707，0.24，0.0072 1，1，1，2，1，3，3，1，3，2，3，3，5，2 2，1，150 2，1，6，3，25，4，4，0，2，400，4，4，0，4，250 说明: （1 结点荷载的方向与整体坐标方向一致者为正。 （2 非结点荷载的方向与单元坐标方向一致者为正。当单元坐标方向由左向右时,正向应向上,所以,这时q向下为负。 输出文件命名为APFexI.OUT。计算完成后，此文件中记录了全部输入的数据以供核对，并记录了算得的各结点的位移和各杆端的内力。 部分计算结果: ELE. FN I FN J FS I FS J M I M J 1 -29.484 29.484 101.463 -101.463 -221.875 -183.979 2 48.537 -48.537 -29.484 179.484 183.979 442.928 3 596.029 -596.029 48.537 -48.537 -150.487 -43.660 4 0.000 0.000 416.545 233.455 -399.268 0.000 例 APFEX2 试求图示刚架由于支座位移所引起的内力。已知各杆EI 37500k