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

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

2、。输出结果中包括输入的数据（以便校核），角位移的值（以1/EI为单位）以及每个单元的左右两端弯矩值。四. 程序试算1. 算例1 算力图示：12kN8kN.m8kN10kN.m 6kN6kN/m5.3 L3rrlEl输入数据：结点：1 2 3 4 5 6 0;Fp： 0 12 8 0 6 0;运行程序如下：单元：1 2 3 4 5 6; q： 8 0 0 4 0 6;长度：4 6 6 8 4 6;m: 0 0 -8 0 10 0 0El : 1 1.5 1 2 1 1.5;ca TurboC+* IDEThe data you )put in: joint :ljengtli1,52.01,01

3、,5Fp12-B10rrJie angtle<l/EI> :11.303730-9*7805341.053733-1010.34002?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左端弯矩0.00000-14.92439-7.30243-12.3

4、7565-8.16809-7.95197右端弯矩14.92439-0.6975712.3756518.168097.9519723.02401弯矩图如下:2. 算例2算例图示:8kN/m1fnn12kNI 88.m10kN.m6kN6kN/mEl 1.5EI2m2m6m -输入数据：结点：0 1 2 3 4 5 6;Fp： 0 12 8 0 6 0;单元：1 2 3 4 5 6; q： 8 0 0 4 0 6;长度：4 6 6 8 4 6;m: 0 0 -8 0 10 0 0EI : 1 1.5 1 2 1 1.5;运行程序如下:cT TiirboC+* IDEHXThe ddta you 卫

5、in: Ju ill l :0123456unit:123456Leng'tli :4G6S46El:i*01.520l.e1.5Pp:8120PGQq:0Q&40G00-e010 8The angle<l/EI>：1.6«&721-10.06021814.871010-12,18322117,195206-26.597603unit1245£left M-9.82331-12.35339-7.76314-12-S533-6.585&2-14,10360M12.3F339T一环丽12.553931A结果为：角位移为： 丄(0,1.

6、686721,-10.080218,14.871010,-12.183221,17.195206,-26.597603)单元编号123456左端弯矩-9.82331-12.35339-7.76314-12.55393-6.58562-14.10360右端弯矩12.35339-0.2368612.5539316.5856214.103600.00000弯矩图如下: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.5 1 2 1 1.5;Fp： 0 12 8 0 6 0;q： 8 0 0 4

7、0 6;m: 0 0 -8 0 10 0 0运行程序如下:joiJitM12345unit:1 23456lengtli46Ge4GEI-1.01*51.02.0i01*5Fp:0 128S0qHB0040Gneq-0QleQThe angle<l/EI>-1.653925一丫14.2&呛203-10,24025310.BG2OG3unit number123456left H-3.B397S-12,3285?-7-B7793-12-1?318-B.217Z2-7.93794rigfht M12.32059-0.1220712.1931818,217227.?37?423.

8、03103駅 TurboC+* IDEIT hedata <ou putin - X结果为：角位移为： 右(0,1.653925,-9.949034,14.264283,-10.248253,10.062063,0)单元编号123456左端弯矩-9.83978-12.32059-7.87793-12.19318-8.21722-7.93794右端弯矩12.32059-0.1220712.1931818.217227.9379423.03103弯矩图如下:4. 算例4算例图示：12kN818kN .m4El4kN/10kN.m 6kN6kN/m EI 1.5EI.4m 3m 3m3m. 3

9、m丄 2EIEl亠8m pm|2m6m输入数据:结点：1 2 3 4 5 6 7;Fp： 0 12 8 0 6 0;运行程序如下：单元：1 2 3 4 5 6;长度：4 6 6 8 4 6;El : 1 1.5 1 2 1 1.5;q： 8 0 0 4 0 6;m: 0 0 -8 0 10 0 0riPiSdata, t/au Jputin - joint :12unit:12lengtli :4bEI:1.5Fp甘12(1 :Vm:0ca TurboC+* IDEThe anglB<l/EI>：11.364772-1.396211-¥,11326214.663626 -

10、12,12657917.179023 -265Wybllun it nuinber123456left M-0.-1-4.95284-7.188&312.-6.53707-14-11573r itjlit M14.95284-0.8113712.7360016,5376714.11573Q结果为： 角位移为:(11.364772,-1.396211,-9.113262,14.660626,-12.126579,17.179023,-26.589511)单元编号123456左端弯矩-0.00000-14.25984-7.18863-12.73600-6.53707-14.11573右端弯

11、矩14.95284-0.8113712.7360016.5370714.115730.00000弯矩图如下:五源程序#in clude<stdio.h>定义变量 结点编号 单兀编号 长度抗弯刚度#in clude<c oni o.h> int joi nt7; int uni t6; float le ngth6;float EI 6;float P6;float q6;float m7;double I6;int number=7,i,j; double K77=0;集中力均部荷载 集中力偶 线刚度整体刚度矩阵double k622; 单元刚度矩阵double MP6

12、2,Mq62,Mlast62,M7=0,tatleM7,mm7,Mqueue7=0,antiM7=0;依次为：集中力、均布荷载引起的固端弯矩，最终杆端弯矩，全部荷载引起的弯矩，总的结点荷载，按结点编号排列的结点集中力偶，按结点编号排列的固端弯矩，等效结点弯矩double angle7=0; void input(); void solve();void output(); void Gauss();角位移输入函数求解函数输出函数 高斯消元法void main()input();solve();output();void input()clrscr();printf("Please i

13、nput data.nnThe joint number:");for(i=0;i<7;i+)scanf("%d",&jointi);printf("nThe unit number:");for(i=0;i<6;i+)scanf("%d",&uniti);printf("nThe length:"); for(i=0;i<6;i+)scanf("%f",&lengthi);printf("nThe EI:");for(i=

14、0;i<6;i+)scanf("%f",&EIi);printf("nThe Fp:");for(i=0;i<6;i+)scanf("%f",&Pi);printf("nThe q:");for(i=0;i<6;i+)scanf("%f",&qi);printf("nThe m:"); for(i=0;i<6;i+)scanf("%f",&mi);void solve() for(i=0;i<7

15、;i+)if(jointi=0) number-=1;for(i=0;i<6;i+)MPi0=-Pi*lengthi/8;MPi1=Pi*lengthi/8; for(i=0;i<6;i+)Mqi0=-qi*lengthi*lengthi/12;Mqi1=qi*lengthi*lengthi/12; for(i=1;i<6;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<7&&jointi>=1;i+)Mqueu

16、ejointi-1=Mi; mmjointi-1=mi; for(i=0;i<7;i+) antiMi=-Mqueuei; for(i=0;i<7;i+)tatleMi=antiMi+mmi;for(i=0;i<6;i+)Ii=(EIi/lengthi);if(joint0=0) Kjoint1-1joint1-1+=4*I0;if(joint6=0) Kjoint5-1joint5-1+=4*I5; for(i=0;i<6;i+) ki00=4*Ii,ki01=2*Ii,ki10=2*Ii,ki11=4*Ii;for(i=j;jointi>=1&&

17、;jointi+1>=1&&i<6;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 data you put in:ntjoint :");for(i=0;i<7;i+)printf("t");printf("%d",jointi);printf(&quo

18、t;nntunit :");for(i=0;i<6;i+)printf("t");printf("%d",uniti);printf("nntlength :");for(i=0;i<6;i+)printf("t");printf("%1.0f",lengthi);printf("nntEI :");for(i=0;i<6;i+)printf("t");printf("%1.1f",EIi);printf(&

19、quot;nntFp :");for(i=0;i<6;i+)(+9vo.!i)04 ossneo 宀 兰E-J0L%=)tu_dM=F)tu_d(+9vo.!i)04EAUWrfuK 宀 兰一 b-J0L%=)t=u_dM=F)tu_d(+9vo.!i)04 bnuwrfuK 宀 esrJOL%=)tu_dM=F)tu_d011006二 L=_=se_Nq-gcxlL%=c_d (+9v_o")04JylAI lllguwrfuK-(o=_=se_Nq-gcxlL%=c_d (+9v_o")04 =9 IAI UQUWrfuKewun-= P9%=)上 u_d(+9v_o.!i)04-(= qiunu 七 unurfuK =urfuK宀 二一o_6u5J9cxlL%=£u_d(+ qEnuv-o")04M=u 一三(LI/Lo_6ue laurfuK =urfuK宀 A Id三bl/l+- L