Problèmes du couplage hydromécanique des sols. Cas de la consolidation unidimensionnelle de Terzaghi dans un calcul par la méthode des éléments finis

Annexes

Ecriture du programme de calcul de la consolidation unidimensionnelle de Terzaghi dans Matlab :

global R invM

N =20;

h =1/(N+1) ;

z=(0: h :1)';

M=h/6*(4*diag(ones(N,1))+ diag(ones(N-1,1),1)+diag(ones(N-1,1),-1));

invM=inv(M);

R=1/h*(2*diag(ones(N,1))-diag(ones(N-1,1),1)-diag(ones(N-1,1),-1));

eta0=u0(z(2:N+1));

[t,V]=ode23('surpres',[0,1],eta0);

V=[zeros(length(t),1),V, zeros(length(t),1)];

[ZZ,TT]=meshgrid(z,t);

mesh(ZZ,TT,V)

surf (ZZ,TT,V,'EdgeColor','none')

xlabel('Z','Fontsize',12)

ylabel ('Tv','Fontsize',12)

zlabel ('u(Z,Tv)/u0(Z)','Fontsize',12)