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

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

2、（以便校核），角位移的值（以1/EI為單位）以及每個(gè)單元的左右兩端彎矩值。四程序試算1 .算例1算力圖示：10kN.m輸入數(shù)據(jù):結(jié)點(diǎn)：1 2 3 4 5 6 0;單元：1 2 3 4 5 6 ; 長(zhǎng)度：4 6 6 8 4 6 ; EI: 1 1.5 1 2 1 1.5Fp : 0 12 8 0 6 0;m : 0 0 -8 0 10 0 0運(yùn)行程序如下: Turbo C+* IDEThe data you put in： Ju lilt :unitljciiy L li =2.01.-elafTJie angle：11,3C3738-1.434142-8-98053414.053733-10.

3、19210710.348027Linit nufnhvleft Fl-8.16839-7.95197-fl.fi97K77.9K197PH.結(jié)果為：角位移為： 丄(11.383738,-1.434142,-8.980504,14.053733,-10.192107,10.048027,0)單兀編號(hào)123456左端彎矩0.00000-14.92439-7.30243-12.37565-8.16809-7.95197右端彎矩14.92439-0.6975712.3756518.168097.9519723.02401彎矩圖如下:算例圖示:10kN.m12kN8kN8kN/m8kN.m1 :|24E

4、m 松5E3mf SnEJ 3mf4kN/mFTTTT12EI8m6kN7 EI1-5EI /；：|2m|2m6kN/m6m輸入數(shù)據(jù):結(jié)點(diǎn)：O 1 2 3 4 5 6;單元：1 2 3 4 5 6 ; 長(zhǎng)度：4 6 6 8 4 6 ; EI: 1 1.5 1 2 1 1.5m : 0 0 -8 0 10 0 0Fp : 0 12 8 0 6 0;運(yùn)行程序如下:| TurboC+* IDEQThe data tjou plit in *joint :12345bunit:123456Lengtli :4GGG46EI:1.51.02.01.01.5Fp:G128S60q:B00406rt:eQ-

5、8QiaQrThe angle-1.606721-10,0002101-1.071010-12.10322117.17520G26.577603unit number123456left n-7.B2331-12,3533?-7-76314-12.5533-6.5B562-14.10360r isrht M12.3533?-B.23&8&12.5539316,5856214.103&0結(jié)果為:角位移為： 丄(0,1.686721,-10.080218,14.871010,-12.183221,17.195206,-26.597603)單兀編號(hào)123456左端彎矩-9.82331-12.3533

6、9-7.76314-12.55393-6.58562-14.10360右端彎矩12.35339-0.2368612.5539316.5856214.103600.00000彎矩圖如下:3.算例3算例圖示：10kN.m8k4m 一 3m 3m輸入數(shù)據(jù)：結(jié)點(diǎn)：0 1 2 3 4 5 0;單元：Fp : 0 12 8 0 6 0；q4kN/m2EI8m6kN /6kN/m6昭驚I1 2 3 4 5 6 ; 長(zhǎng)度：4 6 6 8 4 6 ; EI: 1 1.5 1 2 1 1.5:8 0 0 4 0 6 ;m : 0 0 -8 0 10 0 0運(yùn)行程序如下:joint 12J45unit1 23456

7、lengtli466e4GEl-1.01*51.02.0 1.01*5Fp0 128S60q-B0040Gneq-0Q10QThe ongtle-1.653925一丫94903414.2&吩203-10,24025310.HG2OG3unit number12345left ri03979-12.32B5?-7.B77S3-12-1?318-B.21722rigfht M12.32059-0-1220712.1931630,217227.?37?4寅 TurboC+* IDEIT liedata ijdu putin :甘23.S3103結(jié)果為：角位移為：古(0,1.653925,-9.949

8、034,14.264283,-10.248253,10.062063,0)單兀編號(hào)123456左端彎矩-9.83978-12.32059-7.87793-12.19318-8.21722-7.93794右端彎矩12.32059-0.1220712.1931818.217227.9379423.03103彎矩圖如下:4 .算例4算例圖示：10kN.m輸入數(shù)據(jù):N6kN4kN/m Hr 6 i6kN/m2EI 丄El丄8m 2m2miEL7結(jié)點(diǎn)：1 2 3 4 5 6 7;單元：1 2 3 4 5 6 ; 長(zhǎng)度：4 6 6 8 4 6 ; EI: 1 1.5 1 2 1 1.5m : 0 0 -8

9、 0 10 0 0Fp : 0 12 8 0 6 0;運(yùn)行程序如下:結(jié)果為:角位移為：(11.364772,-1.396211,-9.113262,14.660626,-12.126579,17.179023,-26.589511)單兀編號(hào)123456左端彎矩-0.00000-14.25984-7.18863-12.73600-6.53707-14.11573右端彎矩14.95284-0.8113712.7360016.5370714.115730.00000彎矩圖如下:五源程序#in clude#in clude定義變量int joi nt7;纟口點(diǎn)編號(hào)int uni t6;單兀編號(hào)floa

10、t le ngth 6;長(zhǎng)度f(wàn)loat El 6;抗彎剛度f(wàn)loat P6;集中力float q6;均部何載集中力偶線剛度整體剛度矩陣單元?jiǎng)偠萬(wàn)loat m7;double I6;int number=7,i,j;double K77=0;double k622; 矩陣doubleMP62,Mq62,Mlast62,M7=0,tatleM7,mm7,Mqueue7=0,antiM7=0;依次為：集中力、均布荷載引起的固端彎矩，最終桿端彎矩，全部荷載引起的彎矩，總的結(jié)點(diǎn)荷載，按結(jié)點(diǎn)編號(hào)排列的結(jié)點(diǎn)集中力偶，按結(jié)點(diǎn)編號(hào)排列的固端彎矩，等效結(jié)點(diǎn)彎矩double angle7=0;角位移void inp

11、ut();輸入函數(shù)void solve();求解函數(shù)void output();輸出函數(shù)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;i6;i+)scanf(%d,&uniti);printf(nThe length:);for(i=0;i6;i+)scanf(%f,&l

12、engthi);printf(nThe EI:);for(i=0;i6;i+)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=0;i6;i+)MPi0=-Pi*lengthi/8;MPi1=Pi*lengthi/8;for(i=0;i6;i

13、+)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+)Mqueue jointi-1=Mi;mmjointi-1=mi;for(i=0;i7;i+)antiMi=-Mqueuei;for(i=0;i7;i+)tatleMi=antiMi+mmi;for(i=0;i6;i+)Ii=(EIi/lengthi);if( join

14、t0=0) Kjoint1-1joint1-1+=4*I0;if( joint6=0) Kjoint5-1joint5-1+=4*I5;for(i=0;i=1&jointi+1=1&i6;i+)K jointi-1jointi-1+=ki00;K jointi-1jointi+1-1+=ki01;K jointi+1-1jointi-1+=ki10;K jointi+1-1jointi+1-1+=ki11;getch();void output()clrscr();printf(The data you put in:ntjoint:);for(i=0;i7;i+)printf(%d,join

15、ti);printf(nntunit :);for(i=0;i6;i+)printf(t);printf(%d,uniti);printf(nntlength :);for(i=0;i6;i+)printf(t);printf(%1.0f,lengthi);printf(nntEI :);for(i=0;i6;i+)printf(%1.1f,EIi);printf(nntFp :);for(i=0;i6;i+)printf(t);printf(%1.0f,Pi);printf(nntq :);for(i=0;i6;i+)printf(t);printf(%1.0f,qi);printf(nnt

16、m :);for(i=0;i6;i+)printf(t);printf(%1.0f,mi);Gauss();for(i=0;i6;i+)Mlasti0=ki00*anglejointi-1+ki01*angle0;Mlasti1=ki10*anglejointi-1+ki11*angle1;printf(n);printf(nThe angle(1/EI):nn);for(i=0;inumber;i+)printf(%12.6f,anglei);jointi+1-1+MPi0+Mqijointi+1-1+MPi1+MqiII);printf(nunit number);for(i=0;i6;i+)printf(%6d ,uniti);printf(nnleft M t);for(i=0;i6;i+) printf(%12.5f,Mlasti0);printf(nnright Mt);for(i=0;i6;i+) printf(%12.5f,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(Ki j!=0)BOXm=Ki