| 查看: 321 | 回复: 1 | ||
13097320684新虫 (小有名气)
|
[求助]
matlab 已有1人参与
|
|
M=input('请输入弯矩:'); H=input('请输入剪力:'); h_1=input('请输入h_1'); h_2=input('请输入h_2'); q_1=input('请输入q_1'); q_2=input('请输入q_2'); q_3=input('请输入q_3'); q_4=input('请输入q_4'); z=input('请输入桩的入土深度z:'); g=input('请输入桩的变形系数g:'); E=input('请输入下段弹性模量:') I=input('请输入下段惯性距:') E_1=input('请输入上段弹性模量:') I_1=input('请输入上段惯性距:') A=[2.60000000000000,0.0331460000000000,1.75473000000000,2.90670000000000,2.72365000000000,-1.81479000000000,-0.926020000000000,1.33485000000000,2.74972000000000,-2.62126000000000,-3.59987000000000,-1.87734000000000,0.916790000000000,-2.43695000000000,-5.14023000000000,-5.35541000000000,-2.82106000000000;2.80000000000000,-0.385490000000000,1.49037000000000,3.12843000000000,3.28769000000000,-2.38756000000000,-1.17548300000000,0.841170000000000,2.86653000000000,-3.10341000000000,-4.71748000000000,-3.10791000000000,0.197290000000000,-2.34558000000000,-6.02299000000000,-6.99007000000000,-4.44490000000000;3,-0.928090000000000,1.03679000000000,3.22471000000000,3.85838000000000,-3.05319000000000,-2.82410000000000,0.0683700000000000,2.80406000000000,-3.54058000000000,-5.99979000000000,-4.68788000000000,-0.891260000000000,-1.96928000000000,-6.76460000000000,-8.84029000000000,-6.51972000000000;3.50000000000000,-2.92799000000000,-1.27172000000000,2.46304000000000,4.97982000000000,-4.98062000000000,-6.70806000000000,-3.58647000000000,1.27018000000000,-3.91921000000000,-9.54367000000000,-10.3404000000000,-5.85402000000000,1.07408100000000,-6.78895000000000,-13.6924000000000,-13.8261000000000;4,-5.85333000000000,-5.94097000000000,-0.926770000000000,4.54780000000000,-6.53316000000000,-12.1581000000000,-10.6084000000000,-3.76647000000000,-1.61428000000000,-11.7306600000000,-17.9186000000000,-15.0755000000000,9.24368000000000,-0.356720000000000,-15.6105000000000,-23.1404000000000]; t=abs(A(i,1)-z*g); indexi=find(t==min(t)); i=i(indexi); A1=A(i,2); B1=A(i,3); C1=A(i,4); D1=A(i,5); A2=A(i,6); B2=A(i,7); C2=A(i,8); D2=A(i,9); A3=A(i,10); B3=A(i,11); C3=A(i,12); D3=A(i,13); A4=A(i,14); B4=A(i,15); C4=A(i,16); D4=A(i,17); M0=M+H*(h_1+h_2); HO=H; HH=(B2*D1-B1*D2)/((g^3*E*I)*(A2*B1-A1*B2)); MH=(A2*D1-A1*D2)/((g^2*E*I)*(A2*B1-A1*B2)); HM=(B2*C1-B1*C2)/((g^2*E*I)*(A2*B1-A1*B2)); MM=(A2*C1-A1*C2)/((g*E*I)*(A2*B1-A1*B2)); X0=HO*HH+MO*HM; J0=-(H0*MH+MO*MM); M_z=g^2*E*I*(X0*A3+(J0*B3)/g+(MO*C3)/(g^2*E*I)+(HO*D3)/(g^3*E*I)); Q_z=g^3*E*I*(X0*A4+(J0*B4)/g+(MO*C4)/(g^2*E*I)+(HO*D4)/(g^3*E*I)); WYO=(H*(1/3*(n*h_1^3+h_2^3)+n*h_1*h_2*(h_1+h_2)))/(E_1*I_1)+(M*(h_2^2+n*h_1*(2*h_2+h_1)))/(2*E_1*I_1); WY=X0-J0*(h_1+h_2)+WY0; fprintf('桩柱顶的水平位移=%6.3f',WY); 为什么提示超出索引范围啊 |
回复此楼» 猜你喜欢
一志愿重庆大学085700资源与环境总分308求调剂
已经有6人回复
-
268求调剂
已经有7人回复
-
302求调剂
已经有11人回复
-
306求调剂
已经有4人回复
-
材料学学硕080502 337求调剂-一志愿华中科技大学
已经有4人回复
-
求调剂
已经有6人回复
-
313求调剂
已经有3人回复
-
南昌大学材料专硕311分求调剂
已经有6人回复
-
316求调剂
已经有6人回复
-
346求调剂[0856]
已经有7人回复
霜小妹二
木虫 (正式写手)
哈哈
- 应助: 154 (高中生)
- 金币: 1806.7
- 红花: 25
- 帖子: 531
- 在线: 66小时
- 虫号: 3237838
- 注册: 2014-05-27
- 专业: 机械动力学
2楼2016-05-16 22:43:48
