滚动轴承内圈模型

function y=gneiguzhangzi（tx）
global mi mb mo c1 c2 c3 c4 dita1 dita2 dita3 dita4;
epson=0.096;
n=16;
l=10;
rb=0.018;
ri=0.060;
Fd=4215;
Fmax=8430;
Fr=10;
wi=2.4.*pi;
wm=wi.*ri./（ri+rb）;
wb=（wi.*ri）./（2.*rb）;
dd=wi.*t+2*pi/5;
for k=1:1:n;
    cc=wm.*t+2.*（k-1）.*pi./16;
    cc=mod（cc2*pi）;
    if 0        Acita=1;
    elseif  2*pi/5        Acita=0;     
    else
        Acita=1;
    end
end
F=Fmax.*（1-（1-cos（cc））./2/epson）.^1.5;

for k=1:1:n;
if 2.*k.*pi./wi     R=1;
else
    （2.*k.*pi+2.*2*pi/5）./wi    R=0;
end
end
T=2.*pi.*（ri+rb）./（n.*wi.*rb）;
for k=0:1:n-1;
if k.*T    S=1;
else
    k.*T+n.*T.*asin（l./2.*ri）./pi