六跨连续梁内力计算程序程序

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

2、据（以便校核），角位移的值（以1/EI为单位）以及每个单元的左右两端弯矩值。四程序试算1.算例1算力图示：10kN.m结点：1 2 3 4 5 6 0; 单元：1 2 3 4 5 6; 长度：4 6 6 8 4 6; EI : 1 1 2 1Fp： 0 12 8 0 6 0;q： 8 0 0 4 0 6;m 0 0 -8 0 10 0 0运行程序如下:ca TurboC+* IDEThe data you )put in: joint :ljengtli1,52.01,01,5Fp12-B10rrJie angtle :11.303730-9*7805341.053733-1010.34002

3、?kin it numberl&Et ri6.0B008-14.9243?-8.16B09(rigtit M14.92439-0.6975712.375&510.168B97.951J723.02401结果为:角位移为：(11.383738,-1.434142,-8.980504,14.053733,-10.192107,10.048027,0)单兀编号123456左端弯矩右端弯矩弯矩图如下:2.算例2算例图示:10kN.m输入数据:结点：0 1 2 3 4 5 6; 单元：1 2 3 4 5 6; 长度：4 6 6 8 4 6; EI : 1 1 2 1Fp： 0 12 8 0 6 0q：

4、8 0 0 4 0 6m 0 0 -8 0 10 0 0运行程序如下: TurboC+* IDEThedata you putin :joint m12345bim it:123456le rigtli -46ge46EI:1.01,51.0 2.01.01*5Pp-S128Q60X-B004Q6eB-8Q18QTJicangle :L.606721-10.G82iali,071010-12,10322117-ZG.5?C03unit number12345left H-7.B2331-12.3533?-7.7I&314-12.5533-G.5B562-14 .1S36Qr igfht 111

5、2.35339-O.236BI12.5539316.&B5(2L4.103&I3结果为:角位移为：(0,1.686721,-10.080218,14.871010,-12.183221,17.195206,-26.597603)单兀编号123456左端弯矩右端弯矩弯矩图如下:3.算例3算例图示:10kN.m输入数据:结点：0 1 2 3 4 5 0; 单元：1 2 3 4 5 6; 长度：4 6 6 8 4 6; EI : 1 1 2 1Fp： 0 12 8 0 6 0;q： 8 0 0 4 0 6;m 0 0 -8 0 10 0 0运行程序如下：g TurboC+* IDEThe data

6、you putin ijoint =1im it:12lengtli -4EI:1.a1.5Fp:G12XsB0R:eBTJie ansle!I.653925-?i-4283234534566e4G1.02.01.01,58S6040G-0Q10Sunit numberJ23456left H-9-03770-12.3205?-7.B7793-12-19318-B.217Z2-? .93794riht E12.32059-0.1220712.1931B18.217927.?37?423.03103-10.24U25310.862063结果为:角位移为：E| (0,1.653925,-9.949

7、034,14.264283,-10.248253,10.062063,0)单元编号123456左端弯矩右端弯矩弯矩图如下:4.算例4算例图示：结点：1 2 3 4 5 6 7单元：1 2 3 4 5 6; 长度：4 6 6 8 4 6; El : 1 1 2 1;Fp： 0 12 8 0 6 0q： 8 0 0 4 0 6;m 0 0 -8 0 10 0 0运行程序如下:| TurboC+* IDEThe dataflit in ：joint :1 234bb7unit:1 2345ELengtli :46b4bEl:l.U1*51上2Wl.U1*5Fp:U12H甘bM(X:0k*4Mbn:U

8、0-U1U0The angle：11.364772-1.396211-,1132&214.660626-12,12657?17.179023-26-581Vbllunit number123456left M-14.95284-7.188&3-12_7360-6.53707-14.11573riqht M1勺.?5284-0.8113712.7360016-5370714.11573Q结果为:角位移为：1 (11.364772,-1.396211,-9.113262,14.660626,-12.126579,17.179023,-26.589511)单元编号123456左端弯矩右端弯矩弯矩图如

9、下:五源程序#in clude #in clude定义变量int joint7;结点编号int unit6;编号float length6; 长度float EI6;刚度float P6;float q6;float m7;double I6;刚度int number=7,i,j; double K77=0; double k622; 矩阵 double单元抗弯集中力均部荷载集中力偶线整体刚度矩阵单元刚度MP62,Mq62,Mlast62,M7=0,tatleM7,mm7,Mqueue7=0,antiM7=0; 依次为：集中力、均布荷载引起的固端弯矩，最终杆端弯矩，全部荷载引起的弯矩，总的结

10、点荷载，按结点编号排列的结点集中力偶，按结点编号排列的固端弯矩，等效结点弯矩 double angle7=0; 角 位移void input();输入函数void solve();求解函数输出高斯消void output();函数void Gauss();元法void main()input();solve();output();void input()clrscr();printf(Please input data.nnThe joint number:); for(i=0;i7;i+)scanf(%d,&jointi);printf(nThe unit number:); for(i=0

11、;i6;i+)scanf(%d,&uniti);printf(nThe length:); for(i=0;i6;i+)scanf(%f,&lengthi); printf(nThe EI:);scanf(%f,&EIi);printf(nThe Fp:);for(i=0;i6;i+)scanf(%f,&Pi);printf(nThe q:);for(i=0;i6;i+)scanf(%f,&qi);printf(nThe m:);for(i=0;i6;i+)scanf(%f,&mi);void solve()for(i=0;i7;i+)if(jointi=0) number-=1;for(i=

12、0;i6;i+)MPi0=-Pi*lengthi/8;MPi1=Pi*lengthi/8;for(i=0;i6;i+)Mqi0=-qi*lengthi*lengthi/12;Mqi1=qi*lengthi*lengthi/12;for(i=1;i6;i+)Mi=MPi-11+MPi0+Mqi-11+Mqi0;M0=MP00+Mq00;M6=MP51+Mq51;if(joint0=0)j=1;else j=0;for(i=j;i=1;i+)Mqueuejointi-1=Mi;mmjointi-1=mi;for(i=0;i7;i+)antiMi=-Mqueuei;for(i=0;i7;i+)tat

13、leMi=antiMi+mmi;for(i=0;i6;i+)Ii=(EIi/lengthi);if(joint0=0) Kjoint1-1joint1-1+=4*I0;if(joint6=0) Kjoint5-1joint5-1+=4*I5; for(i=0;i=1&jointi+1=1&i6;i+)Kjointi-1jointi-1+=ki00;Kjointi-1jointi+1-1+=ki01;Kjointi+1-1jointi-1+=ki10;Kjointi+1-1jointi+1-1+=ki11; getch();void output()clrscr();printf(The dat

14、a you put in:ntjoint :);for(i=0;i7;i+)printf(t);printf(%d,jointi);printf(nntunit :);for(i=0;i6;i+)printf(t);printf(%d,uniti);printf(nntlength :);for(i=0;i6;i+)printf(%,lengthi);printf(nntEI :);for(i=0;i6;i+)printf(t); printf(%,EIi);printf(nntFp :);for(i=0;i6;i+)printf(t); printf(%,Pi);printf(nntq :)

15、;for(i=0;i6;i+)printf(t); printf(%,qi);printf(nntm :);for(i=0;i6;i+)printf(t);printf(%,mi);Gauss();for(i=0;i6;i+)Mlasti0=ki00*anglejointi-1+ki01*anglejointi+1-1+MPi0 +Mqi0;Mlasti1=ki10*anglejointi-1+ki11*anglejointi+1-1+MPi1 +Mqi1;printf(n);printf(nThe angle(1/EI):nn);for(i=0;inumber;i+)printf(%,ang

16、lei);printf(n);printf(nunit number);for(i=0;i6;i+)printf(%6d ,uniti);printf(nnleft M t);for(i=0;i6;i+) printf(%,Mlasti0);printf(nnright Mt);for(i=0;i6;i+) printf(%,Mlasti1); getch();void Gauss()int l,m;double box;double BOX7=0;for(j=0;j(number-1);j+)for(i=j;inumber;i+)if(Kij!=0)for(m=0;mnumber;m+)BOXm=Kim;Kim=Kj