找回密码
 注册
Simdroid-非首页
查看: 49|回复: 1

计算桩的非线性沉降

[复制链接]
发表于 2011-9-14 09:48:13 | 显示全部楼层 |阅读模式 来自 江苏南京
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=[p r;q t];
    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=');    %容许误差
上面是我写的一个小程序,用来计算桩的非线性沉降,程序运行不了,还请达人帮忙纠错?谢谢
 楼主| 发表于 2011-10-9 16:11:14 | 显示全部楼层 来自 江苏南京
Simdroid开发平台
没有人帮忙啊!
回复 不支持

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-3-28 22:54 , Processed in 0.030783 second(s), 11 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表