计算桩的非线性沉降
function uu=disf(u)ks=A1./(A1+B1*u).^2;
kn=A2./(A2+B2*u).^2;
kb=pi*D^2*kn/4;
g1=sqrt(4*D*ks/(E*(D^2-d^2)));
x=linspace(6,-70,40);
xx=x(2)-x(1);
S=ones(2,2);
for n=1:40
p=(exp(g1*xx)+exp(-g1*xx))/2
r=-g2*(exp(g1*xx)-exp(-g1*xx))/2*g1
q=-g1*(exp(g1*xx)-exp(-g1*xx))/2*g2
t=(exp(g1*xx)+exp(-g1*xx))/2
A=;
S=S*A;
end
uu==-N0*(kb*S(3)-S(4))/(kb*S(1)-S(2));
m=uu(1)-u(1);
if m>e %判断容许误差；
disf(uu);
end
其中参数为：D=input('D=');
d=input('d=');
E=input('E=');
N0=input('N0=');
g2=input('g2=');
A1=input('A1=');%A1，B1对应侧向刚度；
B1=input('B1=');
A2=input('A2=');%A2，B2对应底部土体法向刚度；
B2=input('B2=');
e=input('e='); %容许误差
上面是我写的一个小程序，用来计算桩的非线性沉降，程序运行不了，还请达人帮忙纠错？谢谢
没有人帮忙啊！
页:
[1]