matlab弹性地基梁法程序

%《基础工程》第三章“弹性地基梁法”例题3-1，P86 EI=input('条形基础抗弯刚度EI=\n');

l=input('底面长l=\n');

b=input('底面宽b=\n');

sm=input('预估平均沉降sm=\n');

F1=input('左起1m处集中力F1=\n');

F2=input('左起5.5m处集中力F2=\n');

F3=input('左起11.5m处集中力F3=\n');

F4=input('左起16m处集中力F4=\n');

M1=input('左起1m处弯矩M1=\n');

M4=input('左起16m处弯矩M4=\n');

p0=(F1+F2+F3+F4)/(l*b);

k=p0/(sm/1000);

t=(k*b/(4*EI))^0.25;

r=t*l;s=pi/4;

if r<=s

fprintf('该梁属于短梁\n');

elseif r>s&r

fprintf('该梁属于有限长梁\n');

else

fprintf('该梁属于长梁\n');

end