插床机构导杆机构(

机械原理课程设计

编程说明书

设计题目：插床机构-导杆机构（1）

目 录

1. 设计任务及要求-----------------------------------------------------------------------1

2. 数学模型-------------------------------------------------------------------------------2

3. 程序框图--------------------------------------------------------------------------------3

4. 程序清单--------------------------------------------------------------------------------4

5. 运行结果--------------------------------------------------------------------------------11

6. 心得与体会-----------------------------------------------------------------------------12

7.

参考文献--------------------------------------------------------------------------------13

一．设计任务及要求

已知：程速比系数K=2，滑块5的冲程H=100mm ，中心距23

o o l

=150mm,比值

3O B

BC

l L

=1，各构件重心S 的位置，曲柄每分钟转数1n =60r/min 。

要求：

1. 设计导杆机构；

2. 作机构两个位置的动画显示；

3. 作滑块的运动线图（编程设计）；

4. 编写说明书；

二．数学模型

1. 极位夹角 =60˚

2. 杆长

杆1的长 1l =32o o l )cos αββ-⋅ 杆2的长 2l =αsin 232123221⋅⋅++o o o o l l l l 杆3 4的长 γγγsin 2/)cos (sin 22243-==x l l 3．运动分析

杆1的角速度 60/211n πω= 滑块2的速度 )sin(112βαω-⋅⋅-=l v

滑块2的加速度 22221112)cos()sin(ωβαωβαε⋅+-⋅--⋅⋅-=l l l a 杆3与y 轴夹角 )cos sin arctan(1132α

α

γ⋅⋅+=l l l o o

杆3的角速度 )cos(12

1

3βαωω-⋅=

l l 杆3的角加速度 ]2)sin()cos([1

3511112

3ωβαωβαεεv l l l ----=

杆4与y 轴夹角 )sin arcsin(

4

3l x

l -=γψ 杆4的角速度 ψγωωcos /cos 4334l l =

杆4的角加速度 ψ

ψ

ωγωγεεcos sin cos cos 42

44233334l l l l +-=

滑块5的速度 ψ

ψγωcos )

sin(3

35-=l v 滑块5的加速度 γωγεψωψεcos sin cos sin 23333244445l l l l a --+=

三.程序框图

四、程序清单

#include

#include

#include

#include

#include

#define pi 3.1415926

#define N 600

void init_graph(void);

void initview();

void draw();

void cur();

double weit1[N],weit2[N],weit3[N];

double sita1[N],sita2[N],sita3[N];

double omigar1[N],omigar2[N],omigar3[N]; double

a=50.0,d=160.0,e=130.0,f=115.0,g=115.0,w1=6.2831 852;

main()

{

int i;

double alf=0, detat=0,theta1=0;

double weit=0,w3=0,ekq3=0,b=0,va=0;

double theta2=0,w4=0,ekq4=0;

double s=0,vc=0,ac=0;

double q=0,j=0,u=0;

detat=2*pi/(N*w1);

for(i=0;i

{

alf=w1*detat*i;

weit=atan((a*sin(alf)-d)/(a*cos(alf)));

if(weit<0)

weit=weit+2*pi;

else

if(0

weit=weit+pi;

theta1=weit-pi;

b=sqrt(a*a+d*d-2*a*d*sin(alf));

w3=a*w1*cos(alf-weit)/b;

va=-a*w1*sin(alf-weit);

ekq3=(-a*w1*w1*sin(alf-weit)-2*va*w3)/b; theta2=asin((f*sin(theta1)-e)/g);