模拟齿轮的时变啮合刚度

喜阳阳

我想模拟一下齿轮啮合的时变刚度，可是为什么出来图像是一条直线呢。大神帮忙看一下，这是哪里出了问题

clc
clear
k1=10^9;%单对齿轮啮合刚度
k2=1.6*10^9;%两对齿轮啮合刚度
e=1.3;%重合度
k0=k1.*(e-1)+k2.*(2-e);%平均啮合刚度                                                                                                                                                   
w=2*pi*16*2000/60;%齿轮的啮合频率
n=sym('n');
t=sym('t');
fai=(2-e)*n*pi;%齿轮刚度谐波项相位
fai
kn=2.*(k1-k2)/(n*pi)*sin((2-e)*n*pi);%齿轮刚度谐波项
kn
s1=symsum(kn*cos(n*w*t+fai),n,1,10);
s1
s=k0+k1;%齿轮的时变啮合刚度
t=0:1/10000:1/100;
figure
plot(t,s)