I want the value of 'gama1' to be obtained from the following equation

1 view (last 30 days)
following is the code that i have written for obtaining the value of the 'gama1' variable. The issue i am facing is that when i run the code, i am gettin an error message that 'gama1' is an unknown parameter while i am actually trying to find out 'gama1' using the given FZERO function.
zeta=linspace(0,2);
alpha=10;
k=54;
c=500;
rho=7800;
w=0.5;
m=0.019;
n=0.015;
v=-0.38;
V=15;
t1=0.5;
%Tint(1)=500;
beta=1;
Tw=30;
mu0=3750e6;
h=0.025e-3;
gama_p=0.01;
y=500;
Q=zeros(10,1);
i=[1:1:10];
Tint=zeros(10,1);
for i=1:1:10,
Tint(i)=y;
u=i;
y=y+100;
Q(i)=i;
end
%display(Q)
for i=1:10,
if Tint(i)<=955,
mu(i)=1-(3.44e-4)*Tint(i);
else
A=1-((Tint(i)-955)./(1500-955));
mu(i)=0.68*power(A,1.7);
end;
end;
lamda=atan(mu)*180/pi;
phi = 35 + (0.5*(alpha-(atan(mu)*180/pi)));
vc=V*sind(phi)/cosd(phi-alpha);
%lc=t1*sind(phi+lamda-alpha)./(sind(phi)*cosd(lamda));
lc=t1*(sind(phi-alpha)+(tand(lamda).*cosd(phi-alpha)))/sind(phi);
%strain=tand(phi-alpha)+(1./tand(phi));
ref_strain=(V*cosd(alpha))./(h*cosd(phi-alpha));
ref_stress=mu0*(Tw^v).*power(ref_strain,m);
B=(beta*ref_stress)/(rho*c*Tw);
C=(sind(phi).*cosd(phi-alpha))./cosd(alpha);
D=((rho*(V^2))*sind(phi)*cosd(alpha))./(ref_stress.*cosd(phi-alpha));
K=k*cosd(phi-alpha)./(rho*c*h*V*cosd(alpha));
gmma2=(1./tand(phi))+tand(phi-alpha);
T0= (4.5e-4-(8*Tw*K))./((4.5e-4)-20*K);
T1=(4*T0-Tw)/3;
%gama1=0;
der1=(gama1-gama_p)/(2*h);
% AA=m*C.*power(T0,v)*(gama_p^n)*power(der1,(m-1))* (gama1+gama_p)/(h^2);
% BB= v*C.*power(T0,v-1)*(gama_p^(n))*power(der1,m).*((T1-Tw)/(2*h));
% ZZ=n*C*power(der1,m+1).*power(T0,v)*(gama_p^(n-1))-C.*D*der1;
% function F=hola(gama1,gama_p,h,m,C,T0,v,n,T1,Tw,D)
%F=@(gama1) AA+BB+ZZ;
%syms gama1
%F=inline('m*C.*power(T0,v)*(gama_p^n)*power(der1,(m-1))* (gama1+gama_p)/(h^2))+(v*C.*power(T0,v-1)*(gama_p^(n))*power(der1,m).*((T1-Tw)/(2*h)))+(n*C*power(der1,m+1).*power(T0,v)*(gama_p^(n-1))-C.*D*der1')
F=@(gama1) (m*C.*power(T0,v)*(gama_p^n)*power(der1,(m-1))* (gama1+gama_p)/(h^2))+(v*C.*power(T0,v-1)*(gama_p^(n))*power(der1,m).*((T1-Tw)/(2*h)))+(n*C*power(der1,m+1).*power(T0,v)*(gama_p^(n-1))-C.*D*der1);
gama1=fzero(F,0);
M=(power(T0,v));
N=M';
tou0=C*(((gama1-gama_p)/(2*h))^(m))*N*(gama_p^(n));
tou=rho*(power(V*sin(phi),2)*gama1)+tou0;

Accepted Answer

Purushottama Rao
Purushottama Rao on 4 May 2015
gamma1 variable has been used in the equations prior to its declaration. MATLAB does not understand what is gamma1 until it faces gamma1=fzero(F,0)

More Answers (0)

Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!