资源描述：
《两铰拱的内力计算.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、两铰拱的内力计算姓名：王训班级：1133002学号：1113300225两铰拱的内力计算，应用visio画出简图，matlab计算的方式计算出拱弯矩并画出弯矩。每个整数点的弯矩在matlab中算出，为易于观察，用excel制表表示出来。一、计算原理拱的尺寸受力情况计算简图一、计算方法利用力法解此超静定拱由matlab的notebook功能导入word生成m-book文件。二、运算过程和结果在计算简图中已经算出支座反力等基本量用matlab计算过程如下clearsymsxFn=cos(atan(1-x/8))/4;S1=sqrt(1+(1-x/
2、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);d112=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
3、,x,t1)Mn2=M1*X1+Mp2t2=12:1:16;f2=subs(Mn2,x,t2)plot(t1,f1)axis([0,16,-10,20])holdonplot(t2,f2)title('弯矩图')holdoffMc1=-26.8202Mn1=3*x-(((x-8)^2/64-1)*((36*log(2^(1/2)+1))/4375-(51*log(80^(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+
4、(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))/2800))/((2367003059283429563973*log(2^(1/2)+1))/33056565380087516495872000-(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.143
5、9-5.4011-2.82020.59894.85619.951415.8849Mn2=144-(((x-8)^2/64-1)*((36*log(2^(1/2)+1))/4375-(51*log(80^(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
6、800))/((2367003059283429563973*log(2^(1/2)+1))/33056565380087516495872000-(11*log(8*2^(1/2)-8))/56000+(11*log(8*2^(1/2)+8))/56000-2^(1/2)/16800)-9*xf2=15.884910.65646.26622.71400四、计算结果讨论1、弯矩图已在matlab中画出，由弯矩图可以得到每个截面的弯矩。2、为了粗略简洁地观察0-16整数的弯矩，又matlab中f1（x=0-12m）和f2(x=12-16m)计算
7、出来，用excel制表如下。x/m01234567M/KN.m0-3.286-5.7338-7.3436-8.1151-8.0486-7.1439-5.4011x/m8910111213141516M/KN.m-2.82025.59894.85619.951415.884910.65646.26622.7140
