2025年平面連桿機構(gòu)運動學分析與仿真優(yōu)化研究

大作業(yè)(一)(題號：_10B_)二、牛頭刨床機構(gòu)的運動分析方程六.運動線圖及運動分析七.參照書；圖b所示的為一牛頭刨床(Ⅲ級機構(gòu))。假設已知各構(gòu)件的尺寸如表2所示，原動件1以等角速度w?=1rad/s沿著逆時針方向回轉(zhuǎn)，試求各從動件的角位移、角速度和角加速度以及刨頭C點的位移、速度和加速度的變化狀況。b)表2牛頭刨床機構(gòu)的尺寸參數(shù)(單位：mm)hhABC規(guī)定：每三人一組，每人一種題目，每組中至少打印出一份源程序，每人計算出原動件從0"~360*時(N=36)1)位置分析把(式I)寫成投影方程得：由以上各式用型轉(zhuǎn)化法可求得θ?θ?S?S?,θ?=θ?對(式Ⅱ)求一次導數(shù)得：(式Ⅲ)采用高斯消去法可求解(式IV)可解得角速度w3,w4;采用高斯消去法可求解(式V)可得角加速度α?,α4三.程序流程圖讀入：L?,L?,L?,I?,L?,Is,I6,xG及yG和θ2,θ3,θ5及θ6(或LAB,lcp,lDE,h,h,h?及l(fā)BD,θ2,θ3,θ4及θ5)的初值，N,w?,E位置分析位置分析調(diào)用牛頓迭代法子程序求解位置方程(1)求得θ2,θ3,04,θ?及θ6(或lBc,ss)并計算xG及yg調(diào)用系數(shù)矩陣A子程序，并計算A調(diào)用原動件位置參數(shù)列陣B子程序，并計算BJ=1,NB(J)=B(J)w1J=1,N調(diào)用高斯消去法子程序求解速度方程(2)求得w2,W3,J=1,N調(diào)用系數(shù)矩陣A子程序，并計算其矩陣DA調(diào)用系數(shù)矩陣B子程序，并計算其矩陣DB加速度分析加速度分析K=1,NDB(K)=DB(K)w1I=1,NB(K)=-DA(K,IⅡ)w?(Ⅱ)+DB(K)I=1,Nas及α6(或lBC,ac)再求出aEx及aEy0的初值，并計算fNIT>ITmax?NYY求得θ值停止調(diào)用高斯消去法子程序求解A△θ?=fi,求出△θ;解*/voidSolutionaccelerati[12]);/*角加速度求解*/voidFoundmatrixA(double[12]voidFoundmatrixB(double[12],dovoidFoundmatrixDB(double[1//定義全局變量doublel1=180,13=960,14=160,h=900,h1=460,h2=110,asl=1.0;//主函數(shù){doublepsvalue[12],a[N][N],da[N][N],b[N],db[N]//建立文獻，并制表頭if((fp=fopen("shuju","w"{printf("Cann'topenthisfile.\n");}fprintf(fp,"\n//計算數(shù)據(jù)并寫入文獻psvalue[0]=480;psvalue[1]=65*PI/180;psvalue[2]=10*PI/180;psvalue[3]{FoundmatrixB(psvalue,aFoundmatrixDA(psvalFoundmatrixDB(psvalue,anSolutionaccelerationpsvalue[j]=psvalue[j]*180/fprintf(fp,"%12.3f",spsvalue[j]=psvalue[j]*PI/1psvalue[j]+=psvalue[j{printf("%lf\t",shuju[i][jprintf("%1f\t",shuju[i][jprintf("%lf\t",shuju[i][j}}voidSolutionangle(doublevalue[1{s3=value[0];ang3=value[1];ang4=valu{xb=h2+11*cos(ang1);yb=h1+11*sin(ang1);xd=14*cos(ang4);yd=14*sins3=sqrt((xd-xb)*(xd-xb)+(ydang4=atan((yd-yc+h)/(14*value[0]=s3;value[1]=ang3;value[2]=ang4;valu}}voidSolutionspeed(doublea2[N][N],doubleb2[N],doublevalue[12],doubleang1){}voidSolutionaccelerat{{{bk[i]+=-da3[i][j]*value[4}}}voidGaussianE(doublea4[N][N],doubleb4[N],doublep4[N]){if(a4g[0][0]<a4g[1][0]&&a4g[0][1]>{t=b4g[0];b4g[0]=b4g[1];b4g[}if(a4g[2][2]<a4g[3][2]&&a4g[2][3]>{t=b4g[2];b4g[2]=b4g[1];b4g[}{{b4g[i]-=b4g[j]*a4g[i][j]/a4g[j][}p4[0]=b4g[0];p4[1]=b4g[1];p4[2]=b4g[2];p4[3]=b4g[3];}voidFoundmatrixA(doublevalue5[12],doublea5[N][N]){doubles3,ang3,s3=value5[0];ang3=value5[1];ang4=value5a5[0][0]=cos(ang3);a5[0][1]=-s3*sin(ang3);a5[0][2]=-a5[1][0]=sin(ang3);a5[1][1]=s3*cos(ang3);a5[1][2]=a5[2][1]=-13*sin(ang3);a5[2][2]=-14*sin(ang4);a5[3][1]=13*cos(ang3);a5[3][2]=1a5[0][3]=a5[1][3]=a5[2][0]=a5[3][0]}voidFoundmatrixB(doublevalue6[12],doubleangl,doubleb6[N]){b6[0]=-11*sin(angl)*ab6[1]=11*cos(ang1)*a}voidFoundmatrixDA(doublevalue7[12],doubleda7[N][N]){doubles3,ang3,ang4,s5s3=value7[0];ang3=value7[1];ang4=value7s3g=value7[4];as3=value7[5];as4=value7[da7[0][0]=-as3*sin(ang3);da7[0][1]=-s3g*sin(ang3)-s3*cos(ang3)*as3;da7[1][0]=as3*cos(ang3);da7[1][1]=s3g*cos(ang3)-s3*sin(ang3)*as3;da7[2][1]=-13*cos(ang3)*as3;da7[2][2]=-14*cos(angda7[3][1]=-13*sin(ang3)*as3;da7[3][2]=-14*sin(ang4)*as4;da7[0][3]=da7[1][3]=da7[2][0]=da7[2][3]=da7[3][0]=}voidFoundmatrixDB(doublevalue8[12],doubleangl,doubled{}10—B:1AB=180,1CD=960,1DE=160,h=輸出ang1=0時的求解angacceleration:輸出angl=10時的求解輸出angl=20時的求解輸出angl=30時的求解angacceleration:輸出angl=40時的求解輸出angl=50時的求解angleangspeed655.89528783.435335-19.612773260輸出ang1=60時的求解-194.400792輸出ang1=90時的求解-3.858944-199.694701輸出angl=110時的求解angleangspeedangacceleration輸出angl=150時的求解angacceleration:0.6232530.393549輸出angl=170時的求解輸出angl=180時的求解-223.057053輸出ang1=200時的求解輸出ang1=210時的求解297.241960-0.491104輸出angl=220時的求解-0.296427輸出angl=230時的求解27.951831-349.551698輸出angl=240時的求解輸出angl=250時的求解輸出ang1=260時的求解輸出angl=270時的求解輸出angl=280時的求解-0.532397angacceleration:輸出ang1=290時的求解輸出angl=300時的求解輸出ang1=320時的求解輸出ang1=330時的求解輸出angl=340時的求解angacceleration:angacceleration:00五.運動線圖及分析機構(gòu)運動分析：(1)、從θ3-θ1,θ4-θ1曲線圖可知，3桿為擺動導桿，4桿為搖桿，較快，以提高效率；