六跨连续梁内力计算程序程序_1.docx

《六跨连续梁内力计算程序程序_1.docx》由会员分享，可在线阅读，更多相关《六跨连续梁内力计算程序程序_1.docx（7页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、六跨连续梁内力计算程序程序六跨连续梁内力计算程序讲明文档程序适用范围本程序用来解决六跨连续梁在荷载作用下的弯矩计算。荷载能够是集中力Fp作用于跨中、分布荷载q分布全垮、集中力偶m作用于结点的任意组合情况。端部支承可为铰支或固支。程序编辑方法使用TurboC按矩阵位移法的思路进行编辑，用TurboC中的数组来完成矩阵的实现，关键的求解K=P的步骤用高斯消元法。程序使用方法运行程序后，根据提示，依次输入结点编号，单元编号，单元长度，抗弯刚度EI的倍数，集中力，均部荷载，集中力偶，各个数据间用空格隔开，每一项输入完毕后按回车键，所有数据输入完毕后按任意键输出结果。输出结果中包括输入的数据以便校核，角

2、位移的值以1/EI为单位以及每个单元的左右两端弯矩值。程序试算算例1算力图示：输入数据：结点：1234560；单元：123456；长度：466846；EI：1121；Fp：0128060；q：800406；m：00-801000运行程序如下：结果为：角位移为：1(11.383738,-1.434142,-8.980504,14.053733,-10.192107,10.048027,0)EI单元编号123456左端弯矩右端弯矩弯矩图如下：算例2算例图示：6EI8kN/m4m3m2m8mkN/m12365474kN/m3m3m3m2m6m12kN8kN8kN.m6kN10kN.mEIEIEI1.

3、5EI1.52EI输入数据：结点：0123456；单元：123456；长度：466846；EI：1121；Fp：0128060；q：800406；m：00-801000运行程序如下：结果为：角位移为：1(0,1.686721,-10.080218,14.871010,-12.183221,17.195206,-26.597603)EI单元编号123456左端弯矩右端弯矩弯矩图如下：算例3算例图示：6EI8kN/m4m3m2m8mkN/m1236544kN/m3m3m3m2m6m12kN8kN8kN.m6kN10kN.mEI1.5EIEIEI2EI1.5输入数据：结点：0123450；单元：12

4、3456；长度：466846；EI：1121；Fp：0128060；q：800406；m：00-801000运行程序如下：结果为：角位移为：1(0,1.653925,-9.949034,14.264283,-10.248253,10.062063,0)EI单元编号123456左端弯矩右端弯矩弯矩图如下：算例4算例图示：6EI8kN/m4m3m2m8mkN/m12365474kN/m3m3m3m2m6m12kN8kN8kN.m6kN10kN.mEIEIEI1.5EI1.52EI输入数据：结点：1234567；单元：123456；长度：466846；EI：1121；Fp：0128060；q：800

5、406；m：00-801000运行程序如下：结果为：角位移为：1(11.364772,-1.396211,-9.113262,14.660626,-12.126579,17.179023,-26.589511)EI单元编号123456左端弯矩右端弯矩弯矩图如下：五源程序#include#include定义变量intjoint7;结点编号intunit6;单元编号floatlength6;长度floatEI6;抗弯刚度floatP6;集中力floatq6;均部荷载floatm7;集中力偶doubleI6;线刚度intnumber=7,i,j;doubleK77=0;整体刚度矩阵doublek62

6、2;单元刚度矩阵doubleMP62,Mq62,Mlast62,M7=0,tatleM7,mm7,Mqueue7=0,antiM7=0;依次为：集中力、均布荷载引起的固端弯矩，最终杆端弯矩，全部荷载引起的弯矩，总的结点荷载，按结点编号排列的结点集中力偶，按结点编号排列的固端弯矩，等效结点弯矩doubleangle7=0;角位移voidinput();输入函数voidsolve();求解函数voidoutput();输出函数voidGauss();高斯消元法voidmain()input();solve();output();voidinput()clrscr();printf(Pleaseinputdata.nnThejointnumber:);for(i=0;i=1;i+)Mqueuejointi-1=Mi;mmjointi-1=mi;for(i=0;i=1&jointi+1=1&i=j;m-)Klm+=-Kjm*Klj;tatleMnumber-1/=Knumber-1number-1;Knumber-1number-1=1;for(i=0;i=0;i-)for(j=number-1;ji;j-)anglei=anglei-Kij*anglej;