凸轮理论轮廓曲线设计c语言程序

/************

习题

4-1

凸轮轮廓曲线绘制程序

**********/

#include

#include

#define PI 3.1415926

void main()

{

double e=10.0,ro=40.0,rt=10 .0,h=20.0,phi1=150.0,phis=30.0,phi_1=120.0,phi_s=6 0;

/*

各字母含义：

e

ro,

滚子半径

< br>rt,

行程

h,

推程运动角

phi1,

远休止角

phis,

运动角

phi_1,

近休止角

phi_s*/

double s,alp[73],x[73],y[73],x1,y1,xr[73],yr[73];

/*< /p>

分别代表从动件位移，压力角，理论轮廓曲线点的位置

x,y

，

x

，

y

对运动角的微分，

实际轮廓曲线的位置

x,y* /

double

ic,ic1, ic2,ic3,ic4,so,s1,cop,sip,phi,gam,bel,del=5.0,q,t;

/*

定义

变量

*/

int i;

gam=phi1+phis;

bel=phi1+phis+phi_1;

ic=(int)(360.0/del);

ic1=(int)(phi1/del);

ic2=(int)(gam/del);

ic3=(int)((phi1+phis+phi_1/2.0)/del);

ic4=(int)(bel/del);

so=sqrt(ro*ro-e*e);

printf(

THETA

Theoretical contour curve. Actual contour curve of cam.

pressure

angle

n

/*

在屏幕上输出文件头

*/

printf(

deg

x/mm

y/mm

X/mm

Y/mm

radn

for(i=0;i<=ic;i++)

/*< /p>

每隔五度建立循

环

*/

{

phi=i*del*PI/180.0;

/*

将角 度化为幅

度

*/

cop=cos(phi);

sip=sin(phi);

if(i<=ic1)

/*

推程运动判断

*/

{

s=h/2.0*(1.0-cos(180.0*phi/phi1));

s1=h*PI*sin(180.0*phi/phi1)/2.0/phi1;

x1=-(s+so)*sip+s1*cop-e*cop;

y1=(s+so)*cop+s1*sip-e*sip;

}

else if(i<=ic2)

/*

远休止角

*/

{

s=h;

s1=0;

x1=-(s+so)*sip-e*cop;

y1=(so+s)*cop-e*sip;

}

else

if(i<=ic3)

/*

回程等加速运

动判断

< br>*/

{

s=h-2.0*h*(phi-gam)*(phi- gam)/phi_1/phi_1;

s1=-4.0*h*(phi- gam)/phi_1/phi_1;

x1=s1*cop-(s+so)*sip-e*cop;

y1=s1*sip+(so+s)*cop-e*sip;

}

else

if(i<=ic4)

/*

回程等减速运

动判断

< br>*/

{

s=2.0*h*(bel-phi)*(bel- phi)/phi_1/phi_1;

s1=-4.0*h*(bel- phi)/phi_1/phi_1;

x1=s1*cop-(so+s)*sip-e*cop;

y1=s1*sip+(so+s)*cop-e*sip;

}

else

/*

近休止角

*/

{

s=0;

s1=0;

x1=-(s+so)*sip-e*cop;

y1=so*cop-e*sip;

}

t=fabs(s1-e);

alp[i]=atan2(t,(s+so));

/*

计算压力角

*/

q=rt/sqrt(x1*x1+y1*y1);

x[i]=(s+so)*cop-e*sip;

/*

理论 轮廓曲线位

置

*/

y[i]=(s+so)*sip+e*cop;

xr[i]=x[i]-q*y1;

/*< /p>

实际轮廓曲线位

置

*/

yr[i]=y[i]+q*x1;

phi=i*del;

printf(

}

}

凸轮理论轮廓曲线和实际轮廓 曲线计算数据（以教材东南大学机械学科组郑文伟，吴克坚

主编《机械原理》

p521

页题

4-1

为例）