求弧齿锥齿轮齿面参数方程

就砍三刀

求弧齿锥齿轮齿面参数方程，我自己算的结果有点点夸张！
f =(-sin(w*t*n+atan(V/H))*w*n*(H^2+V^2)^(1/2)*cos(a)*cos(w*t)-((cos(w*t*n+atan(V/H))*(H^2+V^2)^(1/2)+(rw+ah*tan(i))*cos(k))*cos(a)-ah*sin(a))*sin(w*t)*w-cos(w*t*n+atan(V/H))*w*n*(H^2+V^2)^(1/2)*sin(w*t)-(sin(w*t*n+atan(V/H))*(H^2+V^2)^(1/2)+(rw+ah*tan(i))*sin(k))*cos(w*t)*w)*(((cos(w*t*n+atan(V/H))*(H^2+V^2)^(1/2)+(rw+ah*tan(i))*cos(k))*cos(a)-ah*sin(a))*cos(w*t)-(sin(w*t*n+atan(V/H))*(H^2+V^2)^(1/2)+(rw+ah*tan(i))*sin(k))*sin(w*t)-((cos(w*t*n+atan(V/H))*(H^2+V^2)^(1/2)+(rw+ah*tan(i)+cos(i))*cos(k))*cos(a)-(ah-sin(i))*sin(a))*cos(w*t)+(sin(w*t*n+atan(V/H))*(H^2+V^2)^(1/2)+(rw+ah*tan(i)+cos(i))*sin(k))*sin(w*t))-sin(w*t*n+atan(V/H))*w*n*(H^2+V^2)^(1/2)*sin(a)*((cos(w*t*n+atan(V/H))*(H^2+V^2)^(1/2)+(rw+ah*tan(i))*cos(k))*sin(a)+ah*cos(a)-(cos(w*t*n+atan(V/H))*(H^2+V^2)^(1/2)+(rw+ah*tan(i)+cos(i))*cos(k))*sin(a)-(ah-sin(i))*cos(a))+(-sin(w*t*n+atan(V/H))*w*n*(H^2+V^2)^(1/2)*cos(a)*sin(w*t)+((cos(w*t*n+atan(V/H))*(H^2+V^2)^(1/2)+(rw+ah*tan(i))*cos(k))*cos(a)-ah*sin(a))*cos(w*t)*w+cos(w*t*n+atan(V/H))*w*n*(H^2+V^2)^(1/2)*cos(w*t)-(sin(w*t*n+atan(V/H))*(H^2+V^2)^(1/2)+(rw+ah*tan(i))*sin(k))*sin(w*t)*w)*(((cos(w*t*n+atan(V/H))*(H^2+V^2)^(1/2)+(rw+ah*tan(i))*cos(k))*cos(a)-ah*sin(a))*sin(w*t)+(sin(w*t*n+atan(V/H))*(H^2+V^2)^(1/2)+(rw+ah*tan(i))*sin(k))*cos(w*t)-((cos(w*t*n+atan(V/H))*(H^2+V^2)^(1/2)+(rw+ah*tan(i)+cos(i))*cos(k))*cos(a)-(ah-sin(i))*sin(a))*sin(w*t)-(sin(w*t*n+atan(V/H))*(H^2+V^2)^(1/2)+(rw+ah*tan(i)+cos(i))*sin(k))*cos(w*t))

w角速度，n运动的速度比  H  V  rw 刀盘直径  ah 刀盘上点的高度  a齿根高  b旋转角度 i 刀盘角度 k 刀盘上点的坐标点角度 y0齿坯顶端高度
如果有人有研究的话可不可以给我传份？邮箱chenxin5236987410@163.com