Sophia asked . 2020-11-17

help with fzero function

3*(1-f1/alphax)*A1+(f1/3*alphax)*(A2+2*A3/A4)=0 with 

A1 = (sigma_e - sigma_m)/(2*sigma_e - sigma_m);
A2 = (sigma_e - sigma33_com)/(sigma_e + 0.14*(d/L)*(sigma33_com-  sigma_e));
A3 = (sigma_e - sigma11_com);
A4 = sigma_e + ((1-0.14*(d/L))/2)*(sigma11_com - sigma_e);

x3 = fzero(@(sigma_e) trial0(sigma_e,L,d,sigma_m,f1),0.015)

where trial0 is the function given below. Is this correct? I am new to MATLAB and any help will be most appreciated.

function sigma_e1 = trial0(sigma_e,L,d,sigma_m,f1)

global sigma_m
global L
global d
global sigma_e
global f1

A1 = (sigma_e - sigma_m)/(2*sigma_e - sigma_m);
A2 = (sigma_e - sigma33_com)/(sigma_e + 0.14*(d/L)*(sigma33_com-sigma_e));
A3 = (sigma_e - sigma11_com);
A4 = sigma_e + ((1-0.14*(d/L))/2)*(sigma11_com - sigma_e);

sigma_e1 = 3*(1-f1/alphax)*A1+(f1/3*alphax)*(A2+2*A3/A4);
end

I am using the following command lines:

 >> sig3=sigma33_com(0.0025,1.5e-6,1.85,1850,2.5e-7)
sig3 =
         936.23
 >> sig1=sigma11_com(2.5e-3,1.5e-6,1.85,2.5e-7,7.5e-7,1850)
sig1 =
         6.5845
>> x3=fzero(@(sigma_e) trial0(2.5e-3,1.5e-6,1e-9,8e-4,sig3,sig1,sigma_e),0.02)

and the functions are:

function sigma_e1 = trial0(L,d,sigma_m,f1,sigma33_com,sigma11_com,sigma_e)
A1 = (sigma_e - sigma_m)/(2*sigma_e - sigma_m);
A2 = (sigma_e - sigma33_com)/(sigma_e + 0.14*(d/L)*(sigma33_com-sigma_e));
A3 = (sigma_e - sigma11_com);
A4 = sigma_e + ((1-0.14*(d/L))/2)*(sigma11_com - sigma_e);
sigma_e1 = 3*(1-f1/alphax)*A1+(f1/3*alphax)*(A2+2*A3/A4);
end

function sigma11 = sigma11_com(L,d,sigma_s,t,a,sigma_c)
B1 = d^2*sigma_c + 2*(sigma_c + sigma_s)*(t^2 + 2*a*t);
B2 = d^2*sigma_s + 2*(sigma_c + sigma_s)*(t^2 + 2*a*t);
sigma11 = sigma_s/(L+2*t)*(L*B1/B2+2*t);
end

function sigma33 = sigma33_com(L,d,sigma_s,sigma_c,t)
B3 = (L+2*t)*sigma_s*(sigma_c*d^2 + sigma_s*(4*d*t+4*t^2));
B4 = 2*t*sigma_c*d^2+2*t*sigma_s*(4*d*t+4*t^2)+sigma_s*L*(d+2*t)^2;
sigma33 = B3/B4;
end

and alphax is obtained from the function:

function alpha1 = alphax(L,d,t)
alpha1 = d^2*L/((d+2*t)^2*(L+2*t));
end

matlab , matlabsolutions