clc
clear all
clear Figure1
nd=1.3;
A=[2211,2005,1974,2911,932];
m=[0.145,0.168,0.108,0.478,0.064];
E=[203.4,203.4,203.4,193,103.4];
G=[82.7,77.2,77.2,69,41.4];
G=G*10^9;
E=E*10^9;
for i=1:5
d=6.5;
F=351.5/2;
Sut=A(i)/d^m(i);
Sr=0.45*Sut;
Se=0.5*Sr*(1-(0.45/2)^2)^(-1);
c=linspace(-30,30,1000);
O=F*((4*c.^2-c-1)*16./(4*(c-1)*3.14*d^2)+4/(3.14*d^2));
Q=nd*O/Se+(nd*O/Sut).^2-1;
%Q=(nd*F*4/(3.14*Sut*d^2))^2*((4*c.^2-1)./(c-1)).^2+nd*F*4/(3.14*Se*d^2)*((4*c.^2-1)./(c-1))-1;

Qmin=min(abs(transpose(Q)));
for j=1:1000
    if abs(Q(j)-Qmin)<0.005
        C=c(j);
    end
end

hold on;
plot(c,Q);
hold on;
z=linspace(0,0,1000);
plot(c,z,'r');
hold on;
axis([0 30 -6 5]);

d=d/1000;
Lmin=0.377;
Lmax=0.89;
Ymax=Lmax-Lmin;
D=C*d;
k=500/Ymax;
Na=G(i)*d^4/(8*k*D^3);
Nb=Na-G(i)/E(i);
L0=(2*C-1+Nb)*d;
Li=Lmin-L0;
Fi=Li*k;
fprintf('\nQmin= %g C=%g\n',Qmin,C);
fprintf('d=%g 4<C=%g<12 D=%g 3<Na=%g<15 L0=%g \nSut=%g Sr=%g Se=%g k=%g Li=%g Fi=%g\n\n',d*1000,C,D*1000,Na,L0,Sut,Sr,Se,k,Li,Fi);
end