平面杆系计算(FASB).docVIP

平面任意杆件结构内力计算程序（FASB） 一、程序的功能及编制方法 能计算由不同单元构成的平面任意杆件组合结构的内力及结点位移。 程序用直接刚度法集成总刚度矩阵，变带宽一维数组存储总刚度矩阵，先处理法引进支座条件。三角分解法解刚度方程。 二、程序使用方法 使用方法与“APF”程序和“TRUSS”程序相同。 三、数据文件填写格式 数据文件填写格式与“APF”程序基本上相同。输入的变量和数组为： 总信息：NJ，NE，NZ,NB,NP，NF，E0，DS 结点坐标数组：XY(NJ，2） 单元信息数组：JM(NE，5） 单元截面数组：AI(NB，2） 约束信息数组：ZC(NZ，2） 集中荷载信息数组：PJ(NP，3） 非结点荷载信息数组：PF(NF，4） 以上输入数据中只有JM数组和PF数组与“APF”程序不同，其他数组均与“APF”程序 的相应数组的格式和意义完全相同。 单元信息数组JM(NE,5)共有NE行，每行5个数分别填写该单元的左端结点号，右端节点号，左端约束情况和右端约束情况以及杆的类型号。杆端的约束情况是刚性约束填“1”，铰约束填“0”。 非结点荷载信息数组PF（NF，4）为每个非结点荷载填一行，每行四个数分别为荷载作用单元号，荷载类型号和参数C,G值。 非结点荷载的类型共有二类，如图所示。(a)图为类型1,(b)图为类型2。 这里没有设立半跨荷载或局部分布荷载的情况，若遇到此类情况，将荷载变化点设为结点，便可解决问题。若遇到斜杆的情况，可将荷载分解为垂直杆轴和沿杆轴两种情况的叠加来处理。 四、计算例题 例FASBEX1 试求图示刚架各杆内力。 解：将各结点和单元编号如图示。 由于结构为静定，截面尺寸的弹性模量均不 影响内力计算结果。故设E0=1，I=1。考虑 到计算结果与手算结果相比较，应忽略杆轴 向变形影响，故设A=10000。输入数据文件 名为FASBEX1.DAT。计算结果输出文件名 为FASBEX1.OUT。 输入数据为： 8,7,6,1,2,2,1,0 0,0,0,4,2,4,4,0,4,4,6,4,8,0,8,4 1,2,0,1,1,2,3,1,0,1,3,5,0,1,1 4,5,0,1,1,5,6,0,0,1,6,8,0,1,1,7,8,0,1,1 10000,1 1,1,1,2,4,1,4,2,7,1,7,2 6,2,-40,8,1,40 2,1,-20,-20,3,1,-20,-20 计算的部分结果: 单元 N(i) N(j) Q(i) Q(j) M(i) M(j) 1 20.000 -20.000 0.000 0.000 0.000 0.000 2 0.044 -0.044 20.000 20.000 0.000 0.000 3 0.044 -0.044 -20.000 60.000 0.000 80.000 4 60.000 -60.000 20.000 -20.000 0.000 -80.000 5 -19.945 19.945 0.000 0.000 0.000 0.000 6 -19.934 19.934 -40.000 40.000 0.000 80.000 7 40.000 -40.000 20.000 -20.000 0.000 -80.000 （单位：） 例FASBEX2 平面刚架如图所示。已知EI=常数。 试求各杆内力。 解: 建立名为FASBEX2.DAT的数据文件如下： 5,4,4,1,0,4,1,0 0,0,0,6,4,8,8,6,8,0 1,2,0,1,1,2,3,1,0,1,3,4,0,1,1,4,5,1,0,1 10000,1 1,1,1,2,5,1,5,2 2,1,-16,-16,2,2,-8,-8 3,1,-16,-16,3,2,8,8 说明：由于分布荷载不垂直于单元轴线,计算 时应将其分解为垂直于轴线的分量和平行于轴线的 分量 计算的部分结果: 单元① M(i)