您的位置：360文档中心› 机械原理大作业(用复数矢量法求连杆速度、加速度+c语言程序实现)

机械原理大作业(用复数矢量法求连杆速度、加速度+c语言程序实现)

合集下载

C语言机械原理编程,连杆运动分析图线

C语言机械原理编程,连杆运动分析图线
t=-k*40/20-120;
setcolor(GREEN);
line(m11,t11,m1,t);
m11=m1,t11=t;
k=l3;
m1=n*600/360;
t=-k*40/(0.05)-80;
setcolor(MAGENTA);
line(m12,t12,m1,t);
m12=m1,t12=t;
n+=0.1;
i+=10;
}
outtextxy(-85,-420, "θ2、θ3/度");
}
int main()
{
//创建大小为800 * 600的绘图窗口
initgraph(800, 600);
//设置原点(0, 0)为屏幕中央（Y轴默认向下为正）
setorigin(100, 500);
//使用白色填充背景
setbkcolor(WHITE);
//输出屏幕提示
settextcolor(BLACK);
setfont(18,0,"宋体");
danwei();
danwei1();
settextcolor(RED);
danwei2();
settextcolor(RGB(201, 62, 0));
outtextxy(240,-450, "角度关系图");
void matrixmultiplication(float a_2[2][2],float b_2[2][1],float y[2][1])
{
int i,j,k,a_column=2,a_row=2,b_row=1;
float temp;
for(i=0;i<a_column;i++)

机械原理大作业一连杆机构参考例子

机械原理大作业一连杆机构参考例子

机械原理大作业课程名称：机械原理设计题目：机械原理大作业院系：汽车工程学院车辆工程班级：1101201姓名：。

学号：。

指导教师：游斌弟大作业1 连杆机构运动分析1、运动分析题目如图所示机构，已知机构各构件的尺寸为280mm AB =,350mm BC =,320mm CD =,160mm AD =，175mm BE = 220mm EF =，25mm G x =，80mm G y =，构件1的角速度为110rad/s ω=，试求构件2上点F 的轨迹及构件5的角位移、角速度和角加速度，并对计算结果进行分析。

图 12、对机构进行结构分析该机构由I 级杆组RR （原动件1）、II 级杆组RRR （杆2、杆3）和II 级杆组RPR （滑块4及杆5）组成。

I 级杆组RR ，如图2所示；II 级杆组RRR ，如图3所示；II 级杆组RPR ，如图4所示。

图2 图 2图4 3、建立坐标系建立以点A为原点的固定平面直角坐标系4、各基本杆组运动分析的数学模型（1）同一构件上点的运动分析：如图5所示的构件AB,,已知杆AB 的角速度=10/rad s ω，AB 杆长i l =280mm,可求得B 点的位置B x 、B y ，速度xB v 、yB v ，加速度xB a 、yB a 。

=cos =280cos B i x l ϕϕ; =sin =280sin B i y l ϕϕ;图 3==-sin =-BxB i B dx v l y dt ωϕω； ==cos =;B yB i B dyv l x dt ωϕω222B 2==-cos =-BxB i d x a l x dt ωϕω；2222==-sin =-ByB i Bd y a l y dtωϕω。

图 4（2）RRRII 级杆组的运动分析如图6所示是由三个回转副和两个构件组成的II 级组。

已知两杆的杆长2l 、3l 和两个外运动副B 、D 的位置（B x 、B y 、D x 、D y ）、速度（ xB yB xD yD v v v v 、、、 ） 图6和加速度（xB yB xD yD a a a a 、、、）。

哈工大机械原理大作业1 连杆机构运动分析 牛头刨床 学号5号

哈工大机械原理大作业1 连杆机构运动分析 牛头刨床 学号5号

大作业1 连杆机构运动分析运动分析题目图1-13 牛头刨床机构简图一、按比例画出机构运动简图见A3图纸二、机构的结构分析及基本杆组的划分（1）机构的结构分析该机构为平面连杆机构，活动构件数n=5；有A、B、C、E、F五个转动副和B、F两个移动副，没有高副，故低副数PL=7，高副数PH=0。

机构的自由度F=3n-2PL-PH=3*5-2*7=1。

机构有确定运动轨迹。

（2）基本杆组的划分图1-13中，原动件构件1为Ⅰ级杆组，移除之后按照拆杆组法对机构进行拆分，可得由滑块2和杆3组成的RPRⅡ级杆组，以及由杆4和滑块5组成的RRPⅡ级杆组。

拆分图如下：图1-13-1 Ⅰ级杆组图1-13-2 RPRⅡ级杆组图1-13-3 RRPⅡ级杆组三、各杆组的运动分析数学模型图1-13-1 Ⅰ级杆组 x B =x A +l AB cos ϕy B =y A +l AB sin ϕ图1-13-2 RPR Ⅱ级杆组 x E =x B +(l CE -s)cos ϕj y E =y B +(l CE -s)sin ϕjϕj =arctan(B 0/A 0)A 0=xB -x CB 0=y B -y CS =2020B A + 图1-13-3 RRP Ⅱ级杆组F 的位置 x F =x E -l EF cos ϕiy F =y E +l EF sin ϕi ϕi =arcsin(A 1/l EF )A 1=H-H 1-y EF 的速度F 的加速度四、建立坐标系以A 为原点建立坐标系，则A （0,0）五、编程使用MATLAB 编写程序如下lAB=108;lCD=620;lEF=300;H1=350;H=635;xC=0;yC=-350;syms t;f=(255*pi/30)*t;xB=lAB*cos(f);yB=lAB*sin(f);A0=xB-xC;B0=yB-yC;S=sqrt(A0^2+B0^2);f1=atan(B0/A0);xE=xB+(lAB-S)*cos(f1);yE=yB+(lAB-S)*sin(f1);A1=H-H1-yE;f2=asin(A1/lEF);xF=xE-lEF*cos(f2);yF=yE+lEF*sin(f2);a=0:0.001:0.5;xF=subs(xF,t,a);plot(a,xF)title('位移x随时间t变化图像') xlabel('t(s)'),ylabel(X')lAB=108;lCD=620;lEF=300;H1=350;H=635;xC=0;yC=-350;syms t;f=(255*pi/30)*t;xB=lAB*cos(f);yB=lAB*sin(f);A0=xB-xC;B0=yB-yC;S=sqrt(A0^2+B0^2);f1=atan(B0/A0);xE=xB+(lAB-S)*cos(f1);yE=yB+(lAB-S)*sin(f1);A1=H-H1-yE;f2=asin(A1/lEF);xF=xE-lEF*cos(f2);yF=yE+lEF*sin(f2);vF=diff(xF,t);a=0:0.001:0.5;vF=subs(vF,t,a);plot(a,vF)title('速度v随时间t变化图像') xlabel('t(s)'),ylabel(V')lAB=108;lCD=620;lEF=300;H1=350;H=635;xC=0;yC=-350;syms t;f=(255*pi/30)*t;xB=lAB*cos(f);yB=lAB*sin(f);A0=xB-xC;B0=yB-yC;S=sqrt(A0^2+B0^2);f1=atan(B0/A0);xE=xB+(lAB-S)*cos(f1);yE=yB+(lAB-S)*sin(f1);A1=H-H1-yE;f2=asin(A1/lEF);xF=xE-lEF*cos(f2);yF=yE+lEF*sin(f2);aF=diff(xF,t,2);a=0:0.001:0.5;aF=subs(aF,t,a);plot(a,aF)title('加速度a随时间t变化图像') xlabel('t(s)'),ylabel(A')六、位移、速度、加速度图像位移xE随时间t变化图像速度vE随时间t变化图像加速度aE随时间t变化图像七、运动分析由MATLAB编程后绘制所得图线分析可知，牛头刨床的刨头F行程约为500mm，在0.06s~0.17s之间速度和加速度波动较大，在0.17s~0.29s 之间速度和加速度波动较小，这就保证了刨头在空行程时有急回运动，而在工作行程是由较均匀的切削速度，运动连续，且运动性能良好稳定。

机械原理大作业平面连杆机构报告

机械原理大作业平面连杆机构报告

平面连杆机构的运动分析（题号：平面六杆机构）一、题目说明图示为一平面六杆机构。

设已知各构件的尺寸如表1所示，又知原动件1以等角速度ω1=1 rad/s沿逆时针方向回转，试求各从动件的角位移、角速度及角加速度以及E点的位移、速度及加速度的变化情况。

已知其尺寸参数如下表所示：二、题目分析1、建立封闭图形：L 1 + L 2= L 3+ L 4组号1ll2l 3 l 4 l 5 l 6 α A B C2-A2-B 2-C26.5 67.5 87.5 52.4 43.0 60l 2=116.6 l 2=111.6l 2=126.6L 1 + L 2= L 5+ L 6+AG2、机构运动分析：a 、角位移分析由图形封闭性得：⎪⎪⎩⎪⎪⎨⎧⋅-⋅+=+-⋅-⋅+⋅⋅-⋅+=+-⋅-⋅+⋅⋅=⋅+⋅+⋅=⋅+⋅55662'2221155662'222113322114332211sin sin )sin(sin sin cos cos )cos(cos cos sin sin sin cos cos cos θθθαπθθθθθαπθθθθθθθθL L y L L L L L x L L L L L L L L L L G G 将上式化简可得：⎪⎪⎩⎪⎪⎨⎧=⋅-⋅+-⋅+⋅-=⋅-⋅+-⋅+⋅⋅-=⋅-⋅⋅-=⋅-⋅G G y L L L L L x L L L L L L L L L L L 66552'233466552'2331133221143322sin sin )sin(sin cos cos )cos(cos sin sin sin cos cos cos θθαθθθθαθθθθθθθθb 、角速度分析上式对时间求一阶导数，可得速度方程：化为矩阵形式为：⎪⎪⎩⎪⎪⎨⎧=⋅⋅-⋅⋅+⋅-⋅+⋅⋅=⋅⋅+⋅⋅-⋅-⋅-⋅⋅-⋅⋅-=⋅⋅-⋅⋅⋅⋅=⋅⋅+⋅⋅-0cos cos )cos(cos 0sin sin )sin(sin cos cos cos sin sin sin 66655522'233366655522'2333111333222111333222ωθωθωαθωθωθωθωαθωθωθωθωθωθωθωθL L L L L L L L L L L L L L⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⋅-⋅⋅=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⋅⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⋅-⋅⋅-⋅⋅⋅-⋅--⋅-⋅-⋅⋅⋅-00cos sin cos cos cos )cos(sin sin sin )sin(00cos cos 00sin sin 1111165326655332'26655332'233223322θθωωωωωθθθαθθθθαθθθθθL L L L L L L L L L L L L Lc 、角加速度分析：矩阵对时间求一阶导数，可得加速度矩阵2233222333'223355665'22335566622332233'22sin sin 0cos cos 00sin()sin sin sin cos()cos cos cos cos cos 00sin sin 00cos(L L L L L L L L L L L L L L L L L θθεθθεθαθθθεθαθθθεθθθθθα-⋅⋅⎡⎤⎡⎤⎢⎥⎢⎥⋅-⋅⎢⎥⎢⎥⋅=⎢⎥⎢⎥-⋅--⋅-⋅⋅⎢⎥⎢⎥⋅-⋅⋅-⋅⎣⎦⎣⎦-⋅⋅-⋅⋅⋅-211221123123355665'2223355666cos sin )cos cos cos 0sin()sin sin sin 0L L L L L L L L L θωθωωθθθωθαθθθω⎡⎤⋅⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⋅⎢⎥⎢⎥⎢⎥⋅+⋅⎢⎥⎢⎥⎢⎥⋅⋅-⋅⎢⎥⎢⎥⎢⎥⋅-⋅⋅-⋅⎢⎥⎣⎦⎣⎦⎣⎦d 、E 点的运动状态位移：⎩⎨⎧⋅-⋅+=⋅-⋅+=55665566sin sin cos cos θθθθL L y y L L x x G EG E速度：⎪⎩⎪⎨⎧⋅⋅-⋅⋅=⋅⋅+⋅⋅-=555666555666cos cos sin sin ωθωθωθωθL L v L L v yx E E 加速度：⎪⎩⎪⎨⎧⋅⋅-⋅⋅+⋅⋅+⋅⋅-=⋅⋅+⋅⋅+⋅⋅-⋅⋅-=5552555666266655525556662666cos sin cos sin sin cos sin cos εθωθεθωθεθωθεθωθL L L L a L L L L a yx E E三、流程图四、源程序#include<stdio.h>#include<stdlib.h>#include<math.h>#define PI 3.1415926#define N 4void Solutionangle(double [18],double ); /*矢量法求角位移*/void Solutionspeed(double [N][N],double [N],double [18],double ); /*角速度求解*/void Solutionacceleration(double [N][N],double [N][N],double [N],double [18]);/*角加速度求解*/void GaussianE(double [N][N],double [N],double [N]);/*高斯消去*/void FoundmatrixA(double [18],double [N][N]); //创建系数矩阵Avoid FoundmatrixB(double [18],double ,double [N]);//创建系数矩阵Bvoid FoundmatrixDA(double [18],double [N][N]);//创建矩阵DAvoid FoundmatrixDB(double [18],double ,double [N]);//创建矩阵DB//定义全局变量double l1=26.5,l2=111.6,l3=67.5,l4=87.5,l5=52.4,l6=43.0;double l2g=65.0,xg=153.5,yg=41.7,inang=60*PI/180,as1=1.0;//主函数void main(){int i,j;FILE *fp;double shuju[36][18];double psvalue[18],a[N][N],da[N][N],b[N],db[N],ang1;//建立文件，并制表头if((fp=fopen("filel","w"))==NULL){printf("Cann't open this file.\n");exit(0);}fprintf(fp,"\n The Kinematic Parameters of Point 5\n"); fprintf(fp," ang2 ang3 ang5 ang6"); fprintf(fp," as2 as3 as5 as6"); fprintf(fp," aas2 aas3 aas5 aas6");fprintf(fp," xe ye vex vey aex aey\n");//计算数据并写入文件for(i=0;i<36;i++){ang1=i*PI/18;Solutionangle(psvalue,ang1);FoundmatrixB(psvalue,ang1,b);FoundmatrixA(psvalue,a);Solutionspeed(a,b,psvalue,ang1);FoundmatrixDA(psvalue,da);FoundmatrixDB(psvalue,ang1,db);Solutionacceleration(a,da,db,psvalue);for(j=0;j<4;j++){shuju[i][j]=psvalue[j]*180/PI;}for(j=4;j<18;j++){shuju[i][j]=psvalue[j];}fprintf(fp,"\n");for(j=0;j<18;j++)fprintf(fp,"%12.3f",shuju[i][j]);}fclose(fp);//输出数据for(i=0;i<36;i++){ang1=i*PI/18;printf("\n输出ang1=%d时的求解\n",i*10);printf("angle angspeed angacceleration E:\n"); for(j=0;j<4;j++)printf("%lf\t",shuju[i][j]);printf("\n");for(j=4;j<8;j++)printf("%lf\t",shuju[i][j]);printf("\n");for(j=8;j<12;j++)printf("%lf\t",shuju[i][j]);printf("\n");for(j=12;j<18;j++)printf("%lf\t",shuju[i][j]);printf("\n");}}/*矢量法求角位移*/void Solutionangle(double value[18],double ang1){double xe,ye,A,B,C,phi,alpha,csn,ang5g,d2,d,ang2,ang3,ang5,ang6; A=2*l1*l3*sin(ang1);B=2*l3*(l1*cos(ang1)-l4);C=l2*l2-l1*l1-l3*l3-l4*l4+2*l1*l4*cos(ang1);ang3=2*atan((A+sqrt(A*A+B*B-C*C))/(B-C));if(ang3<0)//限定ang3大小{ang3=2*atan((A-sqrt(A*A+B*B-C*C))/(B-C));}ang2=asin((l3*sin(ang3)-l1*sin(ang1))/l2);xe=l4+l3*cos(ang3)+l2g*cos(ang2-inang);ye=l3*sin(ang3)+l2g*sin(ang2-inang);phi=atan2((yg-ye),(xg-xe));d2=(yg-ye)*(yg-ye)+(xg-xe)*(xg-xe);d=sqrt(d2);csn=(l5*l5+d2-l6*l6)/(2.0*l5*d);alpha=atan2(sqrt(1.0-csn*csn),csn);ang5g=phi-alpha;ang5=ang5g-PI;ang6=atan2(ye+l5*sin(ang5g)-yg,xe+l5*cos(ang5g)-xg);value[0]=ang2;value[1]=ang3;value[2]=ang5;value[3]=ang6;value[12]=xe;value[13]=ye;//限定角度大小for(int i=0;i<4;i++){while(value[i]>2*PI)value[i]-=2*PI;while(value[i]<0)value[i]+=2*PI;}}/*角速度求解*/void Solutionspeed(double a2[N][N],double b2[N],double value[18],double ang1) {double ang2,ang3;ang2=value[0];ang3=value[1];double p2[N];GaussianE(a2,b2,p2);value[4]=p2[0];value[5]=p2[1];value[6]=p2[2];value[7]=p2[3];value[14]=-l3*value[5]*sin(ang3)-l2g*value[4]*sin(ang2-inang);value[15]=l3*value[5]*cos(ang3)+l2g*value[4]*cos(ang2-inang);}/*角加速度求解*/void Solutionacceleration(double a3[N][N],double da3[N][N],double db3[N],double value[18]){int i,j;double ang2,ang3;ang2=value[0];ang3=value[1];double bk[N]={0};double p3[N];for(i=0;i<N;i++){for(j=0;j<N;j++){bk[i]+=-da3[i][j]*value[4+j];}bk[i]+=db3[i]*as1;}GaussianE(a3,bk,p3);value[8]=p3[0];value[9]=p3[1];value[10]=p3[2];value[11]=p3[3];value[16]=-l3*value[9]*sin(ang3)-l3*value[5]*value[5]*cos(ang3)-l2g*value[8]*si n(ang2-inang)-l2g*value[4]*value[4]*cos(ang2-inang);value[17]=l3*value[9]*cos(ang3)-l3*value[5]*value[5]*sin(ang3)+l2g*value[8]*cos (ang2-inang)-l2g*value[4]*value[4]*sin(ang2-inang);}/*高斯消去法解矩阵方程*/void GaussianE(double a4[N][N],double b4[N],double p4[N]){int i,j,k;double a4g[N][N],b4g[N],t;for(i=0;i<N;i++)for(j=0;j<N;j++)a4g[i][j]=a4[i][j];for(i=0;i<N;i++)b4g[i]=b4[i];//使主对角线上的值尽可能大if(a4g[0][0]<a4g[1][0] && a4g[0][1]>a4g[1][1]){for(j=0;j<N;j++){t=a4g[0][j];a4g[0][j]=a4g[1][j];a4g[1][j]=t;}t=b4g[0];b4g[0]=b4g[1];b4g[1]=t;}if(a4g[2][2]<a4g[3][2] && a4g[2][3]>a4g[3][3]){for(j=0;j<N;j++){t=a4g[2][j];a4g[2][j]=a4g[3][j];a4g[3][j]=t;}t=b4g[2];b4g[2]=b4g[1];b4g[3]=t;}//初等行变换for(j=0;j<N;j++)for(i=0;i<N;i++){if(i!=j){for(k=0;k<N;k++)if(k!=j){a4g[i][k]-=a4g[i][j]/a4g[j][j]*a4g[j][k];} b4g[i]-=b4g[j]*a4g[i][j]/a4g[j][j];a4g[i][j]=0;}}for(i=0;i<N;i++)b4g[i]/=a4g[i][i];p4[0]=b4g[0];p4[1]=b4g[1];p4[2]=b4g[2];p4[3]=b4g[3];}//创建系数矩阵Avoid FoundmatrixA(double value5[18],double a5[N][N]){double ang2,ang3,ang5,ang6;ang2=value5[0];ang3=value5[1];ang5=value5[2];ang6=value5[3]; a5[0][0]=-l2*sin(ang2);a5[0][1]=l3*sin(ang3);a5[1][0]=l2*cos(ang2);a5[1][1]=-l3*cos(ang3);a5[2][0]=-l2*sin(ang2)-l2g*sin(ang2-inang);a5[2][2]=l5*sin(ang5);a5[2][3]=l6*sin(ang6);a5[3][0]=l2*cos(ang2)+l2g*cos(ang2-inang);a5[3][2]=-l5*cos(ang5);a5[3][3]=-l6*cos(ang6);a5[0][2]=a5[0][3]=a5[1][2]=a5[1][3]=a5[2][1]=a5[3][1]=0;}//创建系数矩阵Bvoid FoundmatrixB(double value6[18],double ang1,double b6[N]){b6[0]=b6[2]=l1*sin(ang1)*as1;b6[1]=b6[3]=-l1*cos(ang1)*as1;}//创建矩阵DAvoid FoundmatrixDA(double value7[18],double da7[N][N]){double ang2,ang3,ang5,ang6,as2,as3,as5,as6;ang2=value7[0];ang3=value7[1];ang5=value7[2];ang6=value7[3]; as2=value7[4];as3=value7[5];as5=value7[6];as6=value7[7];da7[0][0]=-l2*as2*cos(ang2);da7[0][1]=l3*as3*cos(ang3);da7[1][0]=-l2*as2*sin(ang2);da7[1][1]=l3*as3*sin(ang3);da7[2][0]=as2*(-l2*cos(ang2)-l2g*cos(ang2-inang));da7[2][2]=as5*l5*cos(ang5);da7[2][3]=as6*l6*cos(ang6);da7[3][0]=as2*(-l2*sin(ang2)-l2g*sin(ang2-inang));da7[3][2]=as5*l5*sin(ang5);da7[3][3]=as6*l6*sin(ang6);da7[0][2]=da7[0][3]=da7[1][2]=da7[1][3]=da7[2][1]=da7[3][1]=0; }//创建矩阵DBvoid FoundmatrixDB(double value8[18],double ang1,double db8[N]){db8[0]=db8[2]=l1*as1*cos(ang1);db8[1]=db8[3]=l1*as1*sin(ang1);}五、计算结果及相关曲线图：A组：数据图像角位移角速度角加速度E点的位移、速度及加速度六、体会及建议：这次的作业给了我们一次锻炼机会，让我们基本掌握了平面六杆机构解析算法的基本原理，从理论到程序实践，每一步都环环相扣，每一个结果都息息相关，最终在我们共同的努力下完成了本次作业。

机械原理课程大作业

机械原理课程大作业

机械原理课程大作业基于MATLAB平面连杆机构运动学和动力学分析指导老师：王玉丹目录作业一：平面连杆机构运动学分析第2页作业二：平面连杆机构动力学分析第15页作业一L(AE)=70mm,L(AB)=40mm,L(EF)=60mm,L(DE)=35mm,L(CD)=75m m,L(BC)=50mm,原动件以等角速度W1=10rad/s回转。

试以图解法求在θ1=50°时C点的速度和加速度.对机构进行运动分析，写出C点的位置、速度及加速度方程。

解题过程：令AB=r1, BC=r2, CD=r3, DE=r4，AE=r6，EF=r8, AF=r7，角EAF=θ1。

分析：对机构进行位置分析由封闭形ABCDEA可得：r1+r2=r6+r3+r4 （1）由封闭图形AEFA可得：r7=r6+r8 （2）将（1）（2）两式整理可得：r2-r3-r4=-r1+r6-r8+r7=r6【一】（1）位置方程：【二】速度方程：【三】加速度方程：【四】根据位置方程式编制如下函数：【五】进行数据输入，运行程序进行运算。

根据上面分析的θ1 的极限位置取θ1 的范围为40°-55°并均分成15个元素：输出的P、矩阵的第二列到第四列分别是θ2 、θ3 、4θ4 的值，第一列是AF杆的长度r1’。

【六】第二步根据速度方程式编写如下函数：根据第一步得到的数据进行数据输入，运行程序计算各速度值。

程序如下：程序运行得到q矩阵，第一行到第三行分别是a2、a3、a4 的值，第四行是杆AF上滑块运动的速度，即F点的速度。

【七】第三步编写加速度计算函数：【八】根据第一步和第二步输入数据，运行程序得到各加速度的值：【1】计算C点在θ1 =55°，w1 =10rad/s时的速度，加速度：总结数据绘出各构件的位置、速度和加速度的表格如下：【2】输出图像1）角位置程序及输出的图像：2）F点速度程序及输出的图像：3）角加速度程序及输出的图像：4）F点的加速度程序及输出图像：作业二在图示的正弦机构中，已知：L(AB)=100mm,h1=120mm,h2=80mm, W1=10rad/s(常数)，滑块2和构件3的重量分别为，G2 =40 N 和G3 =100 N,质心S2 和S3 的位置如图所示，加于构件3上的生产阻力Fr=400 N，构件1的重力和惯性力略去不计。

平面连杆机构运动分析c语言版 - 机械原理课程设计

平面连杆机构运动分析c语言版 - 机械原理课程设计

机械原理课程设计平面连杆机构运动分析（VC++语言版）第一章机构运动分析的目的和方法一、 运动分析的目的当设计新的机械时，或是分析现有机械的工作性能时，都必须先计算其机构的运动参数，也就是进行机构的运动分析。

通过位移或轨迹分析，可以确定机构的外廓、某些构件运动所需的空间或判断它们运动时是否相互干涉；可以考察某构件能否实现预定位置变化的要求，或某构件上某点能否实现预定轨迹变化的要求。

通过速度分析，可以确定机构中从动件的速度变化规律能否满足工作要求。

进行生产能力和功率的计算；还能对进一步做机构的加速度分析和受力分析提供必要的数据，也是机构运动方程求解的基础。

通过加速度分析，可以确定机构各构件及构件上某些点的加速度，了解机构加速度的变化规律；是计算惯性力进行动力分析、考察机构动力平衡、防止振动和噪音的基础。

二、 基于直角坐标法的运动分析直角坐标法的要点是以点的运动为基础贯穿整个分析过程。

作运动分析时先求机构中某些特殊点的运动，再根据点的运动计算机构的角运动。

第二章 连杆机构运动分析的直角坐标法原理一、平面机构的组成原理我们知道，机构具有确定运动的条件是其原动件的数目应等于其所具有的自由度的数目。

因此，如将机构及和机架相联的原动件与其余构件拆分开，则有其余构件构成的构件组必然是一个自由度为零的构件组。

而这个自由度为零的构件组，有时还可以再拆成更简单的自由度为零的构件组。

我们把最后不能再拆的最简单的自由度为零的构件组称为组成机构的基本杆组。

含有高副的平面机构，可以按一定条件用低副替代高副而化成平面低副机构。

这样，任何机构都可以看成是有若干个基本杆组依次联结于原动件和机架上而构成的。

最简单的基本杆组是有两个构件和三个低副构成，称为Ⅱ级组，是应用最多的基本杆组，我们把同一构件上两点的运动分析及常见的Ⅱ级组的运动分析分别编写成独立的子程序，对机构进行运动分析时，仅需调用相应的子程序即可。

二、连杆机构运动分析的直角坐标法原理将机构放在直角坐标系中，选择其上的某些特殊点如铰链点、几何关节、机构质心等作为考察对象。

哈工大机械原理大作业连杆机构运动分析29

哈工大机械原理大作业连杆机构运动分析29

机械原理大作业一课程名称：机械原理设计题目：连杆机构运动分析-29院系：机电工程学院班级：设计者：学号：指导教师：设计时间：哈尔滨工业大学一、设计题目(29)二、VB编程得到动画、轨迹、位移、速度、加速度线图1.用VB编写程序得到的机构动画2.构件2上的E点轨迹4.构件7的速度线图三、计算说明书(一)建立机构运动分析数学模型对机构进行结构分析，找出基本杆组1.基本杆组划分a)杆件1为原动件b)构件2、3为RRR型II级杆组c)构件4、5为RRR型II级杆组d)构件6、7为RRP型II级杆组2.结构分析与自由度计算机构为平面连杆机构，活动构件数n=7,L P =10,转动副为9个，移动副为1个，高副为0个，则机构自由度为F=3n-2P L -P H =3x7-2x10=1。

各基本杆组的运动分析数学模型1.原动件（I 级构件）图1（1）位置分析B 点位置方程cos sin B A i iBA i i x x l y y l ϕϕ=+⎧⎨=+⎩ （1-1）（2）速度分析把（1-1）对时间t 求导得B 速度方程sin cos B B A i i i BB A i i i dx x x l dt dy y y l dt ϕϕϕϕ⎫==-⎪⎪⎬⎪==+⎪⎭g g g g g g （1-2）（3）加速度分析把（2-2）对时间他求导得B 加速度方程222222cos sin sin cos B B A i i i i i i BB A i i i i i i d x x x l l dt d y y y l l dt ϕϕϕϕϕϕϕϕ⎫==--⎪⎪⎬⎪==-+⎪⎭g g g g g g g g g g g g g g （1-3） 上式中 i i d dt ϕϕω==g22i i d dtϕϕα==g g 2. RRR 杆组的运动分析3.RRP杆组的运动分析图2（1） 位置分析内运动副C 的位置方程cos cos sin sin sin cos C B i i K j j jCB i i K j j j x x l x s l y y l y s l ϕϕϕϕϕϕ=+=+-⎧⎪⎨=+=++⎪⎩ （2-1）i l 的转角 0arcsin()ji j iA l l ϕϕ+=+ （2-2）式中0()sin ()cos B k j B K j A x x y y ϕϕ=---满足装配条件0j i A l l +≤ 滑块D 相对参考点K 的位移为sin cos cos sin C k j jC K j jjjx x l y y l s ϕϕϕϕ-+--==（2-3）滑块D 的位置方程 cos sin D K iD K ix x s y y s ϕϕ=+⎧⎨=+⎩ （2-4）（2） 速度分析杆i l 的角速度123sin cos j ji i Q Q Q ϕϕωϕ-+==g（2-5）滑块D 沿导路的移动速度 123cos sin i i i iQ l Q l s Q ϕϕ-+=g（2-6）式中1(sin cos )K B j j j j Q x x s l ϕϕϕ=--+ggg2(cos sin )K B j j j j Q y y s l ϕϕϕ=-+-g g g3sin sin cos cos i i j i i j Q l l ϕϕϕϕ=+内运动副C 的速度为sin cos C B i i iC B i i ix x l y y l ϕϕϕϕ⎧=-⎪⎨⎪=+⎩gg gg g g （2-7）外移动副D 的速度为cos sin sin cos D K i j jD K i j j x x s s y y s s ϕϕϕϕϕϕ⎧=+-⎪⎨⎪=++⎩gg g g g g g g（2-8） （3） 加速度分析 杆i l 的角加速度453sin s j ji i Q Q co Q ϕϕεϕ-+==g g（2-9）滑块D 沿导路移动的加速度为453cos sin i i i i Q l Q l s Q ϕϕ--=g g （2-10） 式中224cos (sin cos )(cos sin )2sin K B i i i j j j j j j j j j j Q x x l s l s l s ϕϕϕϕϕϕϕϕϕϕ=-+-+---g g g g g g g g g g225sin (cos sin )(sin cos )2cos K B i i i j j j j j j j j j j Q y y l s l s l s ϕϕϕϕϕϕϕϕϕϕ=-++--++g g g g g g g g g g内运动副C 的加速度22cos sin sin cos C B i i i j i j C B i i i j i jx x l l y y l l ϕϕϕϕϕϕϕϕ⎧=--⎪⎨⎪=-+⎩g gg g g g g g gg g g g g（2-11） 滑块D 的加速度22cos sin cos 2sin sin cos sin 2cos D K j j j j j j jD K j j j j j j jx x s s s s y y s s s s ϕϕϕϕϕϕϕϕϕϕϕϕϕϕ⎧=+---⎪⎨⎪=++-+⎩g g g g g g g g g g gg g g g g g g gg g g（2-12） 以上为基本杆组的通用分析，根据已知条件进行赋值均在程序中体现(二) 计算流程框图（见下页）(三)VB编写程序清单'定义变量Dim xA As Double 'A点的坐标、速度、加速度Dim yA As DoubleDim vxA As DoubleDim vyA As DoubleDim axA As DoubleDim ayA As DoubleDim xB As Double 'B点Dim yB As DoubleDim vxB As DoubleDim vyB As DoubleDim axB As DoubleDim ayB As DoubleDim xD As Double 'D点Dim yD As DoubleDim vxD As DoubleDim vyD As DoubleDim axD As DoubleDim ayD As DoubleDim xE As Double 'E点Dim yE As DoubleDim vxE As DoubleDim vyE As DoubleDim axE As DoubleDim ayE As DoubleDim xG As Double 'G点Dim yG As DoubleDim vxG As DoubleDim vyG As DoubleDim axG As DoubleDim ayG As DoubleDim xH As Double 'H点Dim yH As DoubleDim vxH As DoubleDim vyH As DoubleDim axH As DoubleDim ayH As DoubleDim xM As Double 'M点Dim yM As DoubleDim vxM As DoubleDim vyM As DoubleDim axM As DoubleDim xN As Double '参考点NDim yN As DoubleDim vxN As DoubleDim vyN As DoubleDim axN As DoubleDim ayN As DoubleDim xK As Double 'K点Dim yK As DoubleDim vxK As DoubleDim vyK As DoubleDim axK As DoubleDim ayK As DoubleDim delt1 As Double 'AB杆初始转角Dim L1 As Double 'AB杆长Dim f1 As Double 'AB杆转角Dim w1 As Double 'AB杆角速度Dim e1 As Double 'AB杆角加速度Dim L2 As Double 'BC杆Dim f2 As DoubleDim w2 As DoubleDim e2 As DoubleDim L3 As Double 'CD杆Dim f3 As DoubleDim w3 As DoubleDim e3 As DoubleDim LBE As Double 'BE杆Dim delt2 As DoubleDim L4 As Double 'EF杆Dim f4 As DoubleDim w4 As DoubleDim e4 As DoubleDim LFG As Double 'FG杆Dim f5 As DoubleDim w5 As DoubleDim e5 As DoubleDim deltGH As Double 'GH杆Dim LGH As DoubleDim L6 As Double 'HM（HK）杆Dim f6 As DoubleDim w6 As DoubleDim e6 As DoubleDim Lj As Double 'MK杆Dim wj As DoubleDim ej As DoubleDim ss As Double '移动副K的位移Dim vss As Double '移动副K的速度Dim ass As Double '移动副K的加速度Dim pi As DoubleDim pa As DoubleDim i As DoubleDim fj1 As DoublePrivate Sub Command1_Click() '点E的轨迹Picture1.Scale (-50, 250)-(600, -100)Picture1.Line (-50, 0)-(600, 0) 'XPicture1.Line (0, 250)-(0, -100) 'YFor i = -50 To 600 Step 50 'X轴坐标Picture1.DrawStyle = 2Picture1.Line (i, 250)-(i, -100)Picture1.CurrentX = i - 0: Picture1.CurrentY = 0 Picture1.Print iNext iFor i = -100 To 250 Step 50 'Y轴坐标Picture1.DrawStyle = 2Picture1.Line (-50, i)-(600, i)Picture1.CurrentX = -20: Picture1.CurrentY = i + 7 Picture1.Print iNext iFor fj1 = 0 To 360 Step 0.01f1 = fj1 * paCall RR1Call RRR1Call RR2Picture1.PSet (xE, yE)Next fj1End SubPrivate Sub Command2_Click() '求点G的位移Picture2.Scale (-20, 300)-(380, 200)Picture2.Line (-20, 0)-(380, 0) 'XPicture2.Line (0, 300)-(0, 200) 'YFor i = 0 To 360 Step 30 'X轴坐标Picture2.DrawStyle = 2Picture2.Line (i, 300)-(i, 0)Picture2.CurrentX = i - 10: Picture2.CurrentY = 0Picture2.Print iNext iFor i = 200 To 300 Step 10 'Y轴坐标Picture2.Line (0, i)-(380, i)Picture2.CurrentX = -25: Picture2.CurrentY = i Picture2.Print iNext iFor fj1 = 0 To 360 Step 0.01f1 = fj1 * paCall RR1Call RRR1Call RR2Call RRR2Call RR3Call RRP1Picture2.PSet (fj1, ss)Next fj1End SubPrivate Sub Command3_Click() '求点G的速度Picture3.Scale (-30, 600)-(380, -600)Picture3.Line (-30, 0)-(380, 0) 'XPicture3.Line (0, 600)-(0, -600) 'YFor i = 0 To 360 Step 30 'X轴坐标Picture3.DrawStyle = 2Picture3.Line (i, 600)-(i, -600)Picture3.CurrentX = i - 10: Picture3.CurrentY = 0 Picture3.Print iNext iFor i = -600 To 600 Step 50 'Y轴坐标Picture3.Line (0, i)-(380, i)Picture3.CurrentX = -28: Picture3.CurrentY = i Picture3.Print iNext iFor fj1 = 0 To 360 Step 0.01f1 = fj1 * paCall RR1Call RRR1Call RR2Call RRR2Call RR3Call RRP1Picture3.PSet (fj1, vss)Next fj1End SubPrivate Sub Command4_Click() '求点G的加速度Picture4.Scale (-20, 20000)-(380, -10000)Picture4.Line (-20, 0)-(380, 0) 'XPicture4.Line (0, 20000)-(0, -10000) 'YFor i = 0 To 360 Step 30 'X轴坐标Picture4.DrawStyle = 2Picture4.Line (i, 20000)-(i, -10000)Picture4.CurrentX = i - 10: Picture4.CurrentY = 0 Picture4.Print iNext iFor i = -10000 To 20000 Step 1000 'Y轴坐标Picture4.Line (0, i)-(380, i)Picture4.CurrentX = -25: Picture4.CurrentY = i + 5 Picture4.Print iNext iFor fj1 = 0 To 360 Step 0.01f1 = fj1 * paCall RR1Call RRR1Call RR2Call RRR2Call RR3Call RRP1Picture4.PSet (fj1, ass)Next fj1End SubPrivate Sub Form_Load() '赋初值xA = 0yA = 0vxA = 0vyA = 0axA = 0ayA = 0xD = 730yD = -138vxD = 0vyD = 0axD = 0ayD = 0xG = 465yG = -248vxG = 0vyG = 0axG = 0ayG = 0xN = 545yN = 0vxN = 0vyN = 0axN = 0ayN = 0delt1 = 0L1 = 200w1 = 10e1 = 0L2 = 670L3 = 350delt2 = 0LBE = 335L4 = 380LFG = 130deltGH = 56 * paLGH = 100L6 = 486Lj = 0fj = pi / 2wj = 0ej = 0pi = 3.14pa = pi / 180End SubPublic Sub RR1() '杆ABxB = xA + L1 * Cos(f1 + delt1)yB = yA + L1 * Sin(f1 + delt1)vxB = vxA - w1 * L1 * Sin(f1 + delt1)vyB = vyA + w1 * L1 * Cos(f1 + delt1)axB = axA - w1 ^ 2 * L1 * Cos(f1 + delt1) - e1 * L1 * Sin(f1 + delt1) ayB = ayA - w1 ^ 2 * L1 * Sin(f1 + delt1) + e1 * L1 * Cos(f1 + delt1) End SubPublic Sub RRR1() '杆BCDDim xC As DoubleDim yC As DoubleDim vxC As DoubleDim vyC As DoubleDim axC As DoubleDim ayC As DoubleDim LBD As DoubleDim JCBD As DoubleDim val1 As DoubleDim fDB As DoubleDim C2 As DoubleDim C3 As DoubleDim S2 As DoubleDim S3 As DoubleDim G11 As DoubleDim G12 As DoubleDim G13 As DoubleLBD = Sqr((xD - xB) ^ 2 + (yD - yB) ^ 2)If LBD > L2 + L3 And LBD < Abs(L2 - L3) ThenIf MsgBox("RRR杆组杆长不符合要求", vbOKOnly, "提示") = 1 Then EndElseEnd IfElseEnd IfIf LBD < L2 + L3 And LBD > Abs(L2 - L3) Thenval1 = (L2 ^ 2 + LBD ^ 2 - L3 ^ 2) / (2 * L2 * LBD)JCBD = Atn(-val1 / Sqr(-val1 * val1 + 1)) + 2 * Atn(1)ElseEnd IfIf LBD = L2 + L3 ThenJCBD = 0ElseEnd IfIf LBD = Abs(L2 - L3) ThenIf L2 > L3 ThenJCBD = 0ElseEnd IfIf L2 < L3 ThenJCBD = piElseEnd IfElseEnd IfIf xD > xB And yD >= yB Then '第一象限fDB = Atn((yD - yB) / (xD - xB))ElseEnd IfIf xD = xB And yD > yB ThenfDB = pi * 2ElseEnd IfIf xD < xB And yD >= yB Then '第二象限fDB = pi + Atn((yD - yB) / (xD - xB)) ElseEnd IfIf xD < xB And yD < yB Then '第三象限fDB = pi + Atn((yD - yB) / (xD - xB)) ElseEnd IfIf xB = xD And yD < yB ThenfDB = 3 * pi / 2ElseEnd IfIf xD > xB And yD <= yB Then '第四象限fDB = 2 * pi + Atn((yD - yB) / (xD - xB)) ElseEnd Iff2 = fDB + JCBDxC = xB + L2 * Cos(f2)yC = yB + L2 * Sin(f2)If xC > xD And yC >= yD Then '第一象限f3 = Atn((yC - yD) / (xC - xD))ElseEnd IfIf xC = xD And yC > yD Thenf3 = pi / 2ElseEnd IfIf xC < xD And xC >= xD Then '第二象限f3 = pi + Atn((yC - yD) / (xC - xD))ElseEnd IfIf xC < xD And xC < xD Then '第三象限f3 = pi + Atn((yC - yD) / (xC - xD))ElseEnd IfIf xC = xD And xC < xD Thenf2 = 3 * pi / 2ElseEnd IfIf xC > xD And xC <= xD Then '第四象限f2 = 2 * pi + Atn((yC - yD) / (xC - xD))ElseEnd IfC2 = L2 * Cos(f2)S2 = L2 * Sin(f2)C3 = L3 * Cos(f3)S3 = L3 * Sin(f3)G11 = C2 * S3 - C3 * S2w2 = (C3 * (vxD - vxB) + S3 * (vyD - vyB)) / G11w3 = (C2 * (vxD - vxB) + S2 * (vyD - vyB)) / G11vxC = vxB - w2 * L2 * Sin(f2)vyC = vyB + w2 * L2 * Cos(f2)G12 = axD - axB + w2 ^ 2 * C2 - w3 ^ 2 * C3G13 = ayD - ayB + w2 ^ 2 * S2 - w3 ^ 2 * S3e2 = (G12 * C3 + G13 * S3) / G11e3 = (G12 * C2 + G13 * S2) / G11axC = axB - e2 * L2 * Sin(f2) - w2 ^ 2 * L2 * Cos(f2)ayC = ayB - e2 * L2 * Cos(f2) - w2 ^ 2 * L2 * Sin(f2)End SubPublic Sub RR2() '杆BExE = xB + LBE * Cos(f2 + delt2)yE = yB + LBE * Sin(f2 + delt2)vxE = vxB - w2 * LBE * Sin(f2 + delt2)vyE = vyB + w2 * LBE * Cos(f2 + delt2)axE = axB - w2 ^ 2 * LBE * Cos(f2 + delt2) - e2 * LBE * Sin(f2 + delt2) ayE = ayB - w2 ^ 2 * LBE * Sin(f2 + delt2) + e2 * LBE * Cos(f2 + delt2) End SubPublic Sub RRR2() '杆EFGDim xF As DoubleDim yF As DoubleDim vxF As DoubleDim vyF As DoubleDim axF As DoubleDim ayF As DoubleDim LEG As DoubleDim JFEG As DoubleDim val2 As DoubleDim fGE As DoubleDim C4 As DoubleDim C5 As DoubleDim S4 As DoubleDim S5 As DoubleDim G21 As DoubleDim G22 As DoubleDim G23 As DoubleLEG = Sqr((xG - xE) ^ 2 + (yG - yE) ^ 2)If LEG > L4 + LFG And LEG < Abs(L4 - LFG) ThenIf MsgBox("RRR杆组杆长不符合要求", vbOKOnly, "提示") = 1 Then EndElseEnd IfElseEnd IfIf LEG < L4 + LFG And LEG > Abs(L4 - LFG) Thenval2 = (L4 ^ 2 + LEG ^ 2 - LFG ^ 2) / (2 * L4 * LEG)JFEG = Atn(-val2 / Sqr(-val2 * val2 + 1)) + 2 * Atn(1)ElseEnd IfIf LEG = L4 + LFG ThenJFEG = 0ElseEnd IfIf LEG = Abs(L4 - LFG) ThenIf L4 > LFG ThenJFEG = 0ElseEnd IfIf L4 < LFG ThenJFEG = piElseEnd IfElseEnd IfIf xG > xE And yG >= yE Then '第一象限fGE = Atn((yG - yE) / (xG - xE))ElseEnd IfIf xG = xE And yG > yE ThenfGE = pi * 2ElseEnd IfIf xG < xE And yG >= yE Then '第二象限fGE = pi + Atn((yG - yE) / (xG - xE)) ElseEnd IfIf xG < xE And yG < yE Then '第三象限fGE = pi + Atn((yG - yE) / (xG - xE)) ElseEnd IfIf xE = xG And yG < yE ThenfGE = 3 * pi / 2ElseEnd IfIf xG > xE And yG <= yE Then '第四象限fGE = 2 * pi + Atn((yG - yE) / (xG - xE)) ElseEnd Iff4 = fGE - JFEGxF = xE + L4 * Cos(f4)yF = yE + L4 * Sin(f4)If xF > xG And yF >= yG Then '第一象限f5 = Atn((yF - yG) / (xF - xG))ElseEnd IfIf xF = xG And yF > yG Thenf5 = pi / 2ElseEnd IfIf xF < xG And xF >= xG Then '第二象限f5 = pi + Atn((yF - yG) / (xF - xG))ElseEnd IfIf xF < xG And xF < xG Then '第三象限f5 = pi + Atn((yF - yG) / (xF - xG))ElseEnd IfIf xF = xG And xF < xG Thenf4 = 3 * pi / 2ElseEnd IfIf xF > xG And xF <= xG Then '第四象限f4 = 2 * pi + Atn((yF - yG) / (xF - xG)) ElseEnd IfC4 = L4 * Cos(f4)S4 = L4 * Sin(f4)C5 = LFG * Cos(f5)S5 = LFG * Sin(f5)G21 = C4 * S5 - C5 * S4w4 = (C5 * (vxG - vxE) + S5 * (vyG - vyE)) / G21w5 = (C4 * (vxG - vxE) + S4 * (vyG - vyE)) / G21vxF = vxE - w4 * L4 * Sin(f4)vyF = vyE + w4 * L4 * Cos(f4)G22 = axG - axE + w4 ^ 2 * C4 - w5 ^ 2 * C5G23 = ayG - ayE + w4 ^ 2 * S4 - w5 ^ 2 * S5e4 = (G22 * C5 + G23 * S5) / G21e5 = (G22 * C4 + G23 * S4) / G21axF = axE - e4 * L4 * Sin(f4) - w4 ^ 2 * L4 * Cos(f4)ayF = ayE - e4 * L4 * Cos(f4) - w4 ^ 2 * L4 * Sin(f4)End SubPublic Sub RR3() '杆GHxH = xG + LGH * Cos(f5 + deltGH)yH = yG + LGH * Sin(f5 + deltGH)vxH = vxG - w5 * LGH * Sin(f5 + deltGH)vyH = vyG + w5 * LGH * Cos(f5 + deltGH)axH = axG - w5 ^ 2 * LGH * Cos(f5 + deltGH) - e5 * LGH * Sin(f5 + deltGH) ayH = ayG - w5 ^ 2 * LGH * Sin(f5 + deltGH) + e5 * LGH * Cos(f5 + deltGH) End SubPublic Sub RRP1() '杆组H(M)KDim A0 As DoubleDim Q1 As DoubleDim Q2 As DoubleDim Q3 As DoubleDim Q4 As DoubleDim Q5 As DoubleDim val3 As DoubleA0 = Lj + (yN - yH) * Cos(fj) - (xN - xH) * Sin(fj)val3 = A0 / L6fi = Atn(val3 / Sqr(-val3 * val3 + 1)) + fjxM = xH + L6 * Cos(fi)yM = yH + L6 * Sin(fi)ss = (xM - xN) * Cos(fj) + (yM - yN) * Sin(fj)xK = xN + ss * Cos(fj)yK = yN + ss * Sin(fj)Q1 = vxN - vxH - wj * (ss * Sin(fj) + Lj * Cos(fj))Q2 = vyN - vyH + wj * (ss * Cos(fj) - Lj * Sin(fj))Q3 = L6 * Sin(fi) * Sin(fj) + L6 * Cos(fi) * Cos(fj)w6 = (-Q1 * Sin(fj) + Q2 * Cos(fj)) / Q3vss = -(Q1 * L6 * Cos(fi) + Q2 * L6 * Sin(fi)) / Q3vxM = vxH - w6 * L6 * Sin(fi)vyM = vyH + w6 * L6 * Cos(fi)vxK = vxN + vss * Cos(fj) - ss * wj * Sin(fj)vyK = vyN + vss * Sin(fj) + ss * wj * Cos(fj)Q4 = axN - axH + w6 ^ 2 * L6 * Cos(fi) - ej * (ss * Sin(fj) + Lj * Cos(fj)) - wj ^ 2 * (ss * Cos(fj) - Lj * Sin(fj)) - 2 * vss * wj * Sin(fj)Q5 = ayN - ayH + w6 ^ 2 * L6 * Sin(fi) + ej * (ss * Cos(fj) - Lj * Sin(fj)) - wj ^ 2 * (ss * Sin(fj) + Lj * Cos(fj)) + 2 * vss * wj * Cos(fj)e6 = (-Q4 * Sin(fj) + Q5 * Cos(fj)) / Q3ass = (-Q4 * L6 * Cos(fi) - Q5 * L6 * Sin(fi)) / Q3axM = axH - e6 * L6 * Sin(fi) - w6 ^ 2 * L6 * Cos(fi)ayM = ayH + e6 * L6 * Cos(fi) - w6 ^ 2 * L6 * Sin(fi)axK = axN + ass * Cos(fj) - ss * ej * Sin(fj) - ss * wj ^ 2 * Cos(fj) - 2 * vss * wj * Sin(fj)ayK = ayN + ass * Sin(fj) + ss * ej * Cos(fj) - ss * wj ^ 2 * Sin(fj) + 2 * vss * wj * Cos(fj)End Sub(四)计算结果分析用VB编程可得构件6的位移、速度、加速度图像，由图像可得机构运动一个周期内每个时刻的运动参数值，较为直观。

哈工大机械原理大作业连杆机构9

哈工大机械原理大作业连杆机构9

机械原理大作业大作业一：连杆机构运动分析学生姓名：学号：指导教师：丁刚完成时间：机电工程学院机械设计系制二〇一八年四月连杆机构运动分析1题目（9）图1 设计题目在图1所示的机构中，已知l AB=60mm，l BC=180mm，l DE=200mm，l CD=120mm，l EF=300mm，h=80mm，h1=85mm，h2=225mm，构件1以等角速度ω1=100rad/s 转动。

求在一个运动循环中，滑块5的位移、速度和加速度曲线。

2分析结构1、杆1为RR主动件，绕A以ω1 转动，自由度1.2、4杆和滑块5为RRP II级杆组.，自由度0.3、2，3杆组成II级杆组RRR，自由度0.总共自由度为F=5*3-2*7=1 .由上述的杆组类型，确认出所需运动分析数学模型：同一构件上的点、RRP、RRR。

3.杆组法对平面连杆机构进行运动分析3.1对主动件杆1 RR I级构件的分析主动杆1转角：φ= [0°,360°) δ=0°,则φ’=ω1=100 rad/s角加速度φ’’=0.已知h2=225mm, h=80mm, l AB=60mm 所以A（225mm,80mm）A点速度(0,0),加速度(0,0)B点位置(x A+l AB*cos(φ), Y A+l AB*sin(φ))B点速度(-l AB*sin(φ), l AB*cos(φ)),加速度(-l AB*cos(φ), -l AB*sin(φ))3.2RRRII 级杆组分析（模型参考教材P37-38）图3 如图所示两个构件组成II 级杆组。

已知了B 的位置（x B ，y B ）= (x A +l AB *cos(φ), Y A +l AB *sin(φ)),速度（x ’B ，y ’B ） 和加速度（x ’’B ,y ’’B ）, 已知运动副D （0,0）， 还可知，x ’D =y ’D =0, x ’’D =y ’’D =0. l BC =180 mm, l CD = 120mm所以，x c =x D +l CD *cos(φi)= x B +l BC *cos(φj) y c =x D +l CD *sin(φi)= x B +l BC *sin(φi) 对于φ的求解： A 0=2*l CD （x B -x D ） B 0=2*l BC （y B -y D ） C 0=l CD 2+ l BD 2- l BC 2为了保证机构的装配正常：l BD ≤l CD + l BC AND l BD ≥Abs （l CD - l BC ）可求3杆的转角φi=2*arctan(（B 0±sqrt （A 02 + B 02- C 02））/（A 02+ C 02)），角速度w3=φi ’和角加速度α3= φi ’’3.3 同一构件上的点(模型参考书P35-36)Φiφjφi已知D(0,0),速度(0,0),加速度(0,0),3杆转角φi 角速度φi’角加速度φi’’，Φi和它的导数在3.2都有体现LDE= 200mm可求出E的坐标,速度,加速度.x E =x c+lCE*cos(φi)y E =x C+lCE*sin(φi)同样地，速度、加速度通过求导即可得出算式，可以编出程序。

  1. 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

机原大作业
题3-23
班级：
学号：
姓名：
指导老师：
题目
3-23 在图示机构中，已知原动件1以等角速度w1=10rad/s 逆时针转动，l AB =100mm ， l BC =300mm ，e=30mm 。

试用复数矢量法求构件2的转角θ2、角速度w2和角加速度α2，构件3的速度v 3和加速度a 3。

（1）建立封闭矢量多边形
L1+L2=L4+e (a) L1e iθ1+L2e iθ2=L4 (b)
=>{L1cosθ1+L2cosθ2=L4
L1sinθ1+L2sinθ2=e
(c)
联立可得
L4=[l 22
−(e −l 1sinθ1)2]1/2+l 1cosθ1 （d ）
θ2=arctan [e−l 1cosθ1
l
4−l 1cosθ1
] (e)
(2)速度分析
将试(b)对时间t 求导，可得
il 1w 1e iθ1+il 2w 2e iθ2=l 4
(f) 将上式虚部和实部分开，并联解可得 l 4
=v 3=−[w 1l 1sin (θ1−θ2)]/cosθ2 w 2=−
l 1cosθ1
l 2sinθ2w 1
(3)加速度分析
将式（f）对时间进行求导，可得
-l1w12e iθ1+il2w2e iθ1−l2w22e iθ2=l4
将上式虚部和实部分开，并联解可得
w2=α2=−(w12l1sinθ1+w22l2sinθ2)/(l2cosθ2)
l4=a3=[−w12l1cos(θ1−θ2)+w22l2]/ cosθ2
(4)程序代码（c语言）
#include<stdio.h>
#include<math.h>
#define pi 0.1745329
#define pai 3.1415926
int main()
{
double w1=10,l1=0.1,l2=0.3,e=0.03,j1=0,j2,w2,a2,v3,a3,l4;
double j2hudu;
FILE *fp;
fp=fopen("d:\\jiyuandazuoye.txt","w");
int i;
for(i=0;i<=180;i++)
{
l4=sqrt(fabs(l2*l2-((e-l1*sin(j1*pi))*(e-l1*sin(j1*pi)))))+l1*cos(j1*pi);
j2hudu=atan((e-l1*cos(j1*pi))/(l4-l1*cos(j1*pi)));
j2=pai*j2hudu;
v3=-(w1*l1*sin((j1*pi-j2hudu))/cos(j2hudu));
w2=-(l1*w1*cos(j1*pi))/(l2*sin(j2hudu));
a2=-(w1*w1*l1*sin(j1*pi)+w2*w2*l2*sin(j2hudu))/(l2*cos(j2hudu));
a3=(-w1*w1*l1*cos(pi*j1-j2hudu)+w2*w2*l2)/cos(j2hudu);
fprintf(fp,"j1=%.3f l4=%.3f j2=%.3f w2=%.3f a2=%.3f v3=%.3f a3=%.3f\n",j1,l4,j2,w2,a2,v3,a3);
j1+=2;
}
fclose(fp);
return 0;
}
（5）图像。

相关文档
最新文档