六跨連續(xù)梁內(nèi)力計(jì)算程序程序

1、六跨連續(xù)梁內(nèi)力計(jì)算程序說明文檔程序適用范圍本程序用來解決六跨連續(xù)梁在荷載作用下的彎矩計(jì)算。荷載可以是集中力Fp （作用于跨中）、分布荷載q （分布全垮）、集中力偶m （作用于結(jié)點(diǎn)）的任意組合情況。端部支承可為狡支或固支。程序編輯方法使用Turbo C按矩陣位移法的思路進(jìn)行編輯，用Turbo C中的數(shù)組來完成矩陣的實(shí)現(xiàn)，關(guān)鍵的求解K/二P的步驟用高斯消元法。程序使用方法運(yùn)行程序后，按照提示，依次輸入結(jié)點(diǎn)編號，單元編號，單元長度，抗彎剛度（EI的倍數(shù)），集中力，均部荷載，集中力偶，各個(gè)數(shù)據(jù)間用空格隔開，每一項(xiàng)輸入完畢后按回車鍵，所有數(shù)據(jù)輸入完畢后按任意鍵 輸出結(jié)果。輸岀結(jié)果中包括輸入的數(shù)據(jù)（以便校

2、核），角位移的值（以1/EI為單位）以及每個(gè)單元的左右兩端彎矩值。程序試算算例1算力圖示：12kN8kN. m8kN10kN. m6kN/m輸入數(shù)據(jù)：結(jié)點(diǎn)：1 2 3 4 5 6 0;單元:Fp： 0 12 8 0 6q： 8 0 0 4 0 6 ;0 ;運(yùn)行程序如1 2 3 4 5 6 ;長度：4 6 6 8 4 6 ; EI: 1 1 2 1 ;m： 0 0 -8 0 10 0 0TurboC+* IDEThe data vou. DUt in *joint -1 23456甘im it:1 23456lengtli -466e46El:1.01*51.02.0l.e1*5Fp:0 128

3、S6 0q:B0040 6eQ-0Q10 0The angtle :11.303730434142-10-15210710. 34002?unit number123456left ric. oeoea-14. 5243?-7. 302-13-12. 37565-8.16839-?.?51?7M14,92439-0.6975712. 375G530,168097. 9519723-02401結(jié)果為: 角位移為:古(11. 383738,-1. 434142, -8. 980504, 14. 053733,-10. 192107, 10. 048027, 0)單兀編號123456左端彎矩右端彎矩

4、彎矩圖如下:算例2算例圖示:12kNEl丄E1kN8kN/miTT75ET4m3ni 3m3m 3mlOkN.m一 6kNT腳1. 5EI8m 2m2m 6m 輸入數(shù)據(jù)：結(jié)點(diǎn)：0 1 2 3 4 5 6 ;單元：1 2 3 4 5 6 ;長度：4 6 6 8 4 6 ; El: 112 1;Fp： 0 12 8 0 6 0 ；q： 800406；m： 0 0 -8 0 10 0 0TheddtH you u in: Ju01unit:12Leng tli :4GEl:2*01.5Pp:812q:0Q00Theangle：1. 6&721-10. 06021814. 871010運(yùn)行程序如下:c

5、T TiirboC+* IDE2left M-9.82331-12.35339M12.3F339T 環(huán)麗unit1hEE2345634566S46201. e1.50PGa&40G-eQ10 8-12,18322117,195206-26. 59760345-7.76314-12- S533-6. 585&2-14,1036012 553931A結(jié)果為:角位移為：古(0, 1. 686721,-10. 080218, 14. 871010,-12. 183221, 17. 195206, -26. 597603)單兀編號123456左端彎矩右端彎矩彎矩圖如下:算例3算例圖示:kN/m212kN

6、I 8 dm3m3m|輸入數(shù)據(jù)：lOkN.m4KN- KNftElI I I I Ik7|1.5EI結(jié)點(diǎn)：0 1 2 3 4 5 0；單元：1 2 3 4 5 6；長度：4 6 6 8 4 6； El: 112 1；q： 800406;m： 0 0 -8 0 10 0 05ENm6k6kN/m1.5EI駅 TurboC+* IDEIT hedata ou putin - joijitMl2345unitl 23456lengtli46Ge4GEll. 0l*5l.O2.0i0l*5Fp0 128S0qHB0040GneQ-0QleQThe anglel. 653925Y 14. 2&嗆 203

7、-10, 24025310.BG20G3unit numberl23456left H-3.B397S-12, 3285?-7-B7793-12-I?3l8-B.217Z2-7.93794rigfht M12.32059-0.1220712.1931818,217227. ?37?423. 03103Fp： 0 12 8 0 6 0 ; 運(yùn)行程序如下：結(jié)果為:角位移為：(0, 1. 653925, -9. 949034, 14. 264283,-10. 248253, 10. 062063, 0)算例4算例圖示:輸入數(shù)據(jù): 單兀編號左端彎矩右端彎矩彎矩圖如下:結(jié)點(diǎn)：1 2 3 4 5 6 7;

8、單元：1 2 3 4 5 6 ;長度：4 6 6 8 4 6 ; El: 112 1;Fp： 0 12 8 0 6 0 ; 運(yùn)行程序如下：q： 800406;m： 0 0 -8 0 10 0 0ca TurboC+* IDEThe dat H t/au )put： in * joint :lengtli :2. obhi.o121Hhe anglB：11. 364772-1. 396211-, 11326214. 663626 -12,12657917. 179023 -26 5Wblln it nutober結(jié)果為：角位移為：-14. 95284-7.18863-12.7360011|-0.

9、8113712. 736S0 (11. 364772,-1. 396211,-9. 113262, 14. 660626, -12. 126579, 17. 179023, -26. 589511)單兀編號123456左端彎矩右端彎矩彎矩圖如下:五源程序#in cludeO#in cludeO結(jié)點(diǎn)編號 單兀編號 長度抗彎剛度定義變量int joi nt7; int uni t6j;float len gth6_;float EI 6;集中力均部荷載集中力偶線剛度float P6;float q6;float m7;double I6;int number=7, i, j;double整體剛度矩

10、陣K77 = 0;單元?jiǎng)偠染仃嘾ouble MF機(jī)62, Mlast62, M7 = 0, tatleM7, mm7, Mqueue7 = 0, antiM7 = 0;依次為：集中力、均布荷載引起的固端彎矩，最終桿端彎矩，全部荷載引起的彎矩，總的結(jié)點(diǎn)荷載，按結(jié)點(diǎn)編號排列的結(jié)點(diǎn)集中力偶，按結(jié)點(diǎn)編號排列的固端彎矩，等效結(jié)點(diǎn)彎矩double angle7 = 0:角位移void input (); void solve () ; void輸入函數(shù)求解函數(shù)output () ; voidGauss ();輸出函數(shù)高斯消元法void main()(input (); solve ();output ()

11、;void input ()(clrscr ();printf (Please input data nnThe jointnumber:): for(i二0;i7;i+)scanf&jointi);printfCnThe unit number:);for(i=0;i6;i+)scanf(%d，&uniti);printf(nThe length:); for(i=0;i6;i+) scanf (%f，&length i);printf(nThe El:); for(i=0;i6;i+) scanf (%f, &El i);printf(nThe Fp:); for(i=0;i6;i+) s

12、canf&Pi);printf(z/nThe q:);for(i=0;i6;i+)1 scanf&qi）；printfCnThe m:)； for (i=0; i6; i+)scanf (%f, &mi);)void solve ()-for (i=0; i7: i+) if (jointlilnO) number- =4 ;for(i=0;i6;i+)(MPi 0二-Pi*lengthi/8;MPi 1二Pi*lengthi/8;Jfor(i=0;i6：i+)(Mqi 0二-qi*lengthi*lengthi/12;Mqi l=qi*lengthi*lengthi/12;for(i=l;i

13、6;i+)(Mi二MPi-l l+MPi O+Mqi-1 1+Mqi 0;M0=MP0 0+Mq0 0;M6=MP51+Mq5l;if (joint0=0)j=l;else j=0; for (i=j:i=1;i+)(MqueueEjoint; mmjointi-1二mi;for (i=0;i7;i+) antiMi=-Mqueuei; for(i=0;i7;i+)tatleMi=antiMi+mmi; for(i=0;i6;i+)I i = (EI i/lengthi);if (joint 0=0) Kfjoint 11-1 joint l-l+=4*I0;if(joint6=0) Kjoi

14、ntjoint5-1+二4*15; for(i=0;i=1&joint i+1 =l&i6; i+)(Ktjointjoint0 0;KEjointjointi+l-l+=ki 0 1;Kjointi+lT jointi-l+=ki 1 0;Ktjoint i+11-1 joint i+1-1+二ki 1 Ll;getchO ;void output ()clrscr ();printf (z，The data you put in：ntjoint ; for(i=0;i7;i+)printf(t);printfjointi);)printf(nntunit;for(i=0;i6;i+)pri

15、ntf(t);printfuniti);)printf(nntlength;for(i=0;i6;i+)printf(t);printf lengthi);)printf(nntEI ;for(i=0;i6;i+)printf(t);printf (婷:El i);printf (nntFp :“);for(i=0;i6;i+)printf Ct); printfCV,Pi);)printf(nntq;for(i=0;i6;i 卄)printf(t);printf qi)；printf (z，nntm :)； for(i=0;i6;i 卄)printf(t);printf mi);Gauss

16、();for(i=0;i6;i+)MlastEiO=ki00*anglejointi-l+ki01*anglejointi+1-l+MPi0+Mqi0;MlastEi l=ki 1 0*anglejoint1 1*anglejointi+11+Mqi 1;)printf (n);printf (nThe angle(1/EI):nn);for (i二0; inumber; i+)printfangle Li);）；printf(n printf (，znunit numberzz);for(i=0;i6;i+)printf (，z%6d unit Li);printf (/znnleft M t);for(i=0;i6;i+) printfMlasti0);printf(nnright;for(i=0;i6;i+) printfMlasti 1)； getch();void Gauss ()int 1, m;double box; double BOX7二0;for (j=0;j(number-l);j+)for(i二j;inumber;i+)if(Kij!=O)for (m=0; mnumber; m