资源描述：
《平面杆系结构有限元编制程序.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、平面杆系结构有限元编制程序——斜腿刚构计算一、程序编制说明本程序全部采用VisualBasic语言编写，可以进行简单杆系结构的受力分析，能够精确计算出结构体系的结点位移和内力，并采用文本方式输入和输出数据，方便实用。例子是针对斜腿刚构的计算，全桥共划分为6个单元，计算简图如下。二、子程序划分DATA：数据的准备与输入ZG：总体刚体矩阵DG：单元刚度矩阵TRANSLATE：坐标转换矩阵RESULT：杆端位移与杆端力计算MULT：矩阵相乘MULT1：矩阵与向量相乘GAS：高斯消元法BOUNDARY：边界条件LOAD：外荷载三、标识符说明：1、整数型变量NN：结点
2、总个数ELE：单元总个数GD：固定支座个数JZ：铰支座个数ZC：支承个数MS:系统总自由度数目（总刚阶数）PN：荷载作用的结点数（结点荷载个数）2、一维数组变量EI：杆件抗弯刚度EA：杆件抗拉刚度Cn：cos值Sn：sin值LN：单元左端结点标号数组RN：单元右端结点标号数组FX：水平方向结点外力FY：垂直方向结点外力MOMENT：结点弯矩Z：总体坐标系下结点位移存储P：总体坐标系下结点力存储3、二维数组变量BKE：局部坐标单元刚度矩阵存储TKE：体系总体刚度矩阵存储T：坐标转换矩阵存储四、源程序代码DimI,JAsDoubleDimNNAsInteger,
3、ELEAsInteger,MSAsIntegerDimGDAsInteger,JZAsInteger,ZCAsInteger,PNAsIntegerDimLRN(1To12,1To12)AsDouble,L(1To7)AsDouble,G(1To21,1To21)AsDoubleDimEA(1To7)AsDouble,EI(1To7)AsDouble,CnAsDouble,SnAsDoubleDimX(1To7)AsDouble,Y(1To7)AsDouble,P(1To21)AsDoubleDimTKE(1To21,1To21)AsDouble,TKE1(
4、1To21,1To21)AsDoubleDimBKE(1To6,1To6)AsDouble,T(1To6,1To6)AsDoubleDimTZ(1To6,1To6)AsDouble,Z(1To21)AsDoubleDimLN(1To6)AsInteger,RN(1To6)AsIntegerDimFX(1To7)AsDouble,FY(1To7)AsDouble,MOMENT(1To7)AsDoublePrivateSubForm_Click()Open"e:homeworkgginput.txt"ForInputAs#1Open"e:homework
5、ggoutput.txt"ForOutputAs#2CallData'输入数据CallZG'总刚度矩阵CallBoundary'输入边界条件Callload'输入荷载CallResult'计算结果（杆端力）Close#1Close#2EndSubPublicSubData()'数据的准备与输入Print#2,"平面杆系有限元程序——斜腿刚构计算"Print#2,"========================================="Print#2,"原始数据的输入"Print#2,Print#2,"结点数";Spc(3);"单元数";Spc(3
6、);"固定支座数";Spc(3);"铰支座数";Spc(3);"支承数";Spc(3);"结点荷载数";Spc(3);"系统总自由度数"Input#1,NN,GD,JZ,ZC,PNMS=3*NN:ELE=NN-1Print#2,NN;Spc(6);ELE;Spc(6);GD;Spc(6);JZ;Spc(6);ZC;Spc(6);PN;Spc(6);MSPrint#2,Print#2,"========================================="Print#2,"结点坐标的输入"Print#2,Print#2,"结点";Spc(3);"
7、X";Spc(6);"Y";Spc(3);"单位：米"ForI=1ToNNInput#1,I,X(I),Y(I)Print#2,I;Spc(3);X(I);Spc(6);Y(I)NextIPrint#2,Print#2,"========================================="Print#2,"杆件相关信息的输入"Print#2,Print#2,"杆件号";Spc(3);"左结点号";Spc(3);"右结点号";Spc(3);"抗弯刚度";Spc(3);"抗拉刚度";Spc(3);"杆件长度";Spc(3);"E=2.06E+11
8、N/m2"ForI=1ToELEInput#1,I,
