Help with Delaunay Triangulations plot

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

[KSSV]

What is question?

[Hoang Van Hieu]

Thanks for respone but there is no question. I just want to write a program to demomstrate stresses neer the well bore, that's all.

[KSSV]

What data you have?

[Hoang Van Hieu]

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

[Hoang Van Hieu]

wad i want to plot sigma-e, sigma-d,sigma-z,t-ed,t-ez,t-dz

[darova]

Do you have original equations?

[Hoang Van Hieu]

sadly not, That's why i have to start from the beginning.

[Hoang Van Hieu]

Do you know which type of plot or model is this.

[darova]

I'm afraid i don't

Answer 1

[KSSV]

I suggest to use this package for Delaunay Triangulation.

https://in.mathworks.com/matlabcentral/fileexchange/25555-mesh2d-delaunay-based-unstructured-mesh-generation