两铰拱的内力计算

1、两铰拱的内力计算姓名：王训 班级：1133002学号：1113300225两饺拱的内力计算，应用visio画出简图，matlab计算的方式计算出拱弯矩并画出弯矩。 每个整数点的弯矩在matlab中算出，为易于观察，用excel制表表示出来。一、计算原理。M。F = cHP f1511-1 , c、M = 1 -3 - 8)264=j F2d +j M2, dS EA s S EI sA M M A =p1-d EIX = MA1p511M = M1X1 + Mp拱的尺寸受力情况计算简图4kN/m8m4mm4kN/m基本体系12kNST 1/4X=1X=1单独作用1/4Ml图M1=1-(x-8)

2、A2/648m4m-|- m荷载单独作用4kN/m12kN35KN代梁受力示意图Mc0=128KN.mKTTnB41KN厂Mo 128F =c =32 K二 4 N二、计算方法利用力法解此超静定拱由mat lab的notebook功能导入word生成m-book文件。三、运算过程和结果在计算简图中已经算出支座反力等基本量用matlab计算过程如下clearsyms xFn=cos(atan(1-x/8)/4;S1=sqrt(1+(1-x/8)人2);M1=1-(x-8)人2/64;Mp1=3*x;Mp2=144-9*x;d111=1/5680000*int(Fn人2*S1,x,0,16);d1

3、12=1/28000*int(M1人2*S1,x,0,16);d11=d111+d112;d1p1=1/28000*int(Mp1*M1*S1,x,0,12);d1p2=1/28000*int(Mp2*M1*S1,x,12,16);d1p=d1p1+d1p2;X1=-d1p/d11;Mc=X1;Mc1=vpa(Mc,6)Mn1=M1*X1+Mp1t1=0:1:12;f1=subs(Mn1,x,t1)Mn2=M1*X1+Mp2t2=12:1:16;f2=subs(Mn2,x,t2)plot(t1,f1)axis(0,16,-10,20)hold onplot(t2,f2)title(弯矩图)ho

4、ld offMc1 =-26.8202Mn1 =3*x - (x-8)八2/64 -1)*(36*log(2八(1/2) + 1)/4375 - (51*log(8 0八(1/2) + 4)/8750 - (9*log(24*2八(1/2) + 34)/(3*5八(1/2) + 7)/1750 - (27*log(5八(1/2)/2 + 1/2)/4375 + (93*log(8*2八(1/2) - 8)/17500 + (9*log(8*2八(1/2)+8)/17500+(7*2八(1/2)/625-(17*5八(1/2)/2 8 00)/(236700305 92834295 63973*

5、log(2八(1/2) 1)/3305 656538 0087516495 872 000 - (11*log(8*2八(1/2) - 8)/56000 + (11*log(8*2八(1/2) + 8)/56000 - 2八(1/2)/16800) f1 =0-3.2860-5.7338-7.3436-8.1151-8.0486-7.1439-5.4011-2.82020.59894.85619.951415.8849Mn2 = 144 - (x-8)八2/64 -1)*(36*log(2八(1/2) + 1)/4375 - (51*log(8 0八(1/2) + 4)/8750 - (9*l

6、og(24*2八(1/2) + 34)/(3*5八(1/2) + 7)/1750 - (27*log(5八(1/2)/2 + 1/2)/4375 + (93*log(8*2八(1/2) - 8)/17500 + (9*log(8*2八(1/2)+8)/17500+(7*2八(1/2)/625-(17*5八(1/2)/2 8 00)/(236700305 92834295 63973*log(2八(1/2)+1)/3305 656538 0087516495 872 000 - (11*log(8*2八(1/2) - 8)/56000 + (11*log(8*2八(1/2) + 8)/56000 - 2八(1/2)/16800) - 9*x f2 =15.884910.65646.26622.71400四、计算结果讨论1、弯矩图已在matlab中画出，由弯矩图可以得到每个截面的弯矩。2、为了粗略简洁地观察0-16整数的弯矩，又matlab中fl (x=0-12m)和f2(x=12-16m)计 算出来，用excel制表如下。x/m01234567M/KN.m0-3.286-5.7338-7.