Help with Delaunay Triangulations plot

Question

0 votes

i'm started learning Matlab 1 week ago and the last 5 days i tried to write this but fail. Does any one know has any document of this.

P/s: Sorry for may bad English

9 Comments
Show 7 older comments Hide 7 older comments

Hoang Van Hieu on 26 May 2020

Open in MATLAB Online

here

%nhap thong so dau vao
al = 0.8
gamma = 90
phi = 0
sigma_H = 48
sigma_h = 36
sigma_v = 53
Pwf = 36
Pp = 28
a = 0.4
v = 0.25
%Delta---------------------------------------------------------
d= linspace(0,360,100);
%--------------------------------------------------------------
%khoang cach bat ki tu tam gieng den diem bat ki quanh gieng khoan
e = linspace(a,a+1,50);
e = e';
sigma_x = (sigma_H*cosd(phi)^2 + sigma_h* sind(phi)^2)*cosd(gamma)^2 + sigma_v*sind(gamma)^2;
sigma_y = sigma_H*sind(phi)^2 * sigma_h*cosd(phi)^2;
sigma_zz = (sigma_H*cosd(phi)^2 + sigma_h*sind(phi)^2)*sind(gamma)^2 + sigma_v*cosd(gamma)^2;
t_xy = 0.5*(sigma_h - sigma_H)*sind(2*phi)*cosd(gamma);
t_xz = 0.5*(sigma_H*cosd(phi)^2 + sigma_h*sind(phi)^2 - sigma_v)*sind(2*gamma);
t_yz = 0.5*(sigma_h - sigma_H)*sind(2*phi)*sind(gamma);
sigma_e = ((sigma_H + sigma_h).*0.5).*(1-a.^2./e.^2) + ((sigma_H - sigma_h).*0.5).*(1 + 3.*(a.^4./e.^4) - 4.*(a.^2./e.^2)).*cosd(2.*d) + sigma_x.*(1 + 3.*(a.^4./e.^4) - 4.*(a.^2./e.^2)).*sind(2.*d) + Pwf.*a.^2./e.^2
sigma_d = 0.5.*(sigma_H - sigma_h).*(1+a.^2./e.^2) - 0.5.*(sigma_H -sigma_h).*(1+3.*a.^4./e.^4).*cosd(2.*d) - sigma_x.*(1+3.*a.^4./e.^4).*sind(2.*d) - Pwf.*a.^2./e.^2
sigma_z = sigma_v - v.*(2.*(sigma_H - sigma_h).*(a.^2./e.^2).*cosd(2.*d) + 4.*sigma_x .*(a.^2./e.^2).*sind(2.*d))
t_ed = (-0.5.*(sigma_H - sigma_h).*(1 - 3.*(a.^4./e.^4) + 2.*(a.^2./e.^2)).*sind(2.*d))
t_ez = 0
t_dz = 0