I recently wrote Liu Jing2012 paper in the bearing fault dynamics, encountered difficulties, do not know how to determine the inner ring and outer ring raceway contact connect

L=0.002;

B=0.002;

H=0.75405;

Eta_bd=75.405;

Eta__d=1;

a_1=0.85;

b_1=0.96;

e_1=0.76;

p_1=3.34;

p_2=1.71;

for i=1:Nb

Theta_d(i)=2*pi(i-1)/z+w_c*t ;%计算球的位置

end

m=mod(Theta_d(i),2*pi);

if m>=Theta_0&& m<=Theta_0+T

H_2=h_2*sin(0.5*pi/T*(m-Theta_0));

else

H_2=0;

end

T=ascn(L/D_0);

H_d=0.5*d-((0.5*d)^2-(0.5*B)^2)^0.5;

if H<H_d

H_sd=H;

else

H_sd=H_d ;

end

Theta_ox=0;

for j=1:Nb

Theta_d(j)=2*pi(j-1)/z+w_c*t+Theta_ox;

end

delta_(j)=(x(1)-x(5))*cos(Theta_d(j))+(x(3)-x(7))*sin(Theta_d(j))-Gamma-H-P;

if delta_(j)>0

xi_i=1;

else

xi_i=0;

end

Gamma=sin(0.5*pi/DeltaT(mod(Theta_d(i),2*pi)-Theta_0));

p=p_1+(p_2-p_1)*sin(Gamma);