matlab concrete Stake soil mechanics force calculation program

This week a fellow brother idle nothing to graduate the way to help him write a program by the force of a pile,

According to his reference papers to write formulas always felt made out graduation atmosphere,

Thesis was Nonlinear analysis of pile load-settlement behavior in layered soil are interested go and see only a section

Saying so nice inside curve in the end is how come. . .

Soil 3 into three main functions, plus two iterative procedure

Currently calculated strange to complement posted on a program Cato reference see it.

Because according to the paper in an iterative one pile up too tired to count 65W times, matlab operational efficiency we are all aware of. . . So we saved a lot of

function [ U,P ] = diedaihanshu( u,Z,P1)
%UNTITLED3 此处显示有关此函数的摘要
%   此处显示详细说明
%第三阶段基本参数
format long g
a=1.018e-5;
b=5.683e-8;
fsu=55.525e3;
emtsu=0;
Ai=1.0465e7;
if u<=4.626e-3
    fu=u/(b+a*u);
else
    fu=emtsu*u+fsu;
end
fu;
X=[];
fudao= 1/(b + a*u) - (a*u)/(b + a*u)^2;
fsi=fudao*u+(fu-fudao*u);
X=[X;fu,fudao,fsi];
B=-X(:,2);
emt1=sqrt(-B/Ai);
emt2=-sqrt(-B/Ai);
P1=2*1.196*207e3*u/(1-(0.3)^2);
Ci=-fu+fudao*u;
c1i=exp(-emt1*(Z))*(emt2*u +P1/Ai+emt2*Ci/B)/(emt2-emt1);
c2i=exp(-emt2*(Z))*(emt1*u +P1/Ai+emt1*Ci/B)/(emt1-emt2);
U=c1i*exp(emt1*(Z-0.1))+c2i*exp(emt2*(Z-0.1))-Ci/B;
syms x
y=c1i*exp(emt1*x)+c2i*exp(emt2*x)-Ci/B;
diff(y);
p=subs(diff(y),x,Z-0.1);
P=-35e6*(pi/4)*(1.196)^2*eval(p);
end

The following iterative procedure

p=[];
u=[];
X=[];
Y=[];
for i=2e-4:2e-4:4e-3
    x=i;
    P1=2*1.196*207e3*x/(1-(0.3)^2);
    for j=46.7:-0.1:30.3
        Z=j;
        [U,P]=diedaihanshu1(x,Z,P1);
        u=[u;U];
        p=[p;P];
        x=U;
        P1=P;
    end
    X=[X;u(length(u))]
    Y=[Y;p(length(p))]
end
u,p,X,Y;
plot(X,Y,'o')

The abscissa is an iterative pattern ordinate displacement U of each segment before the force reaches the growth curve can be seen to a certain value is setting up. 

After my analysis is the reason 

if u<=4.626e-3
    fu=u/(b+a*u);
else
    fu=emtsu*u+fsu;
end

Here circulation discrimination caused by a particular theory I did not care ////

However, this curve further glance like thing, no mutations substantially smooth.

Iteration 2

X=[];
Y=[];
u1=1e-6
p1=2*1.196*207e3*u1/(1-(0.3)^2);
for Z=46.7:-0.1:0
    if Z>30.3
        [U,P]=diedaihanshu( u1,Z,p1);
        X=[X,Z];
        Y=[P,Y];
        u1=U;
        p1=P;
    elseif Z<=10.0
        [U,P]=diedaihanshu3(u1,Z,p1);
        X=[X,Z];
        Y=[P,Y];
        u1=U;
        p1=P;
    else 
        [U,P]=diedaihanshu2(u1,Z,p1);
        X=[X,Z];
        Y=[P,Y];
        u1=U;
        p1=P;
    end
end

plot(X,Y,'o')

Outputs of the three sections a strange line

This I can only think it gives initial data is more toxic, poor soil is relatively large?

In fact, I have always felt very messed soil mechanics, a bunch of empirical formula really okay.

Hey tomorrow to my brother who took the break to see the results of mentors, let us look forward to the next bar.