资源描述：
《机械原理大作业一(平面连杆机构的运动分析).doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、大作业（一）平面连杆机构的运动分析（题号：_10B_）学校：西北农林科技大学学院：机械与电子工程学院指导老师：郭红利一．题目及原始数据；二、牛头刨床机构的运动分析方程三．计算程序框图；四．计算源程序；五．计算结果；六．运动线图及运动分析七．参考书；一、题目及原始数据；图b所示的为一牛头刨床（Ⅲ级机构）。假设已知各构件的尺寸如表2所示，原动件1以等角速度ω1=1rad/s沿着逆时针方向回转，试求各从动件的角位移、角速度和角加速度以及刨头C点的位移、速度和加速度的变化情况。F'Fy5C6Ｇ2BAh1h2h13x4EDb）表2牛头刨床机构的尺寸参数（单位：mm）题号lAB
2、lCDlDEhh1h2ABC7—A180960160900460110h2=120h2=135h2=140要求：每三人一组，每人一个题目，每组中至少打印出一份源程序，每人计算出原动件从0゜～360゜时（N=36）各运动变量的大小，并绘出各组对应的运动线图以及E点的轨迹曲线。二、牛头刨床机构的运动分析方程1）位置分析建立封闭矢量多边形由图可知=，故未知量有、、、。利用两个封闭图形ABDEA和EDCGE，建立两个封闭矢量方程，由此可得：把（式Ⅰ）写成投影方程得：（式Ⅱ）由以上各式用型转化法可求得，解：高斯消去法求解2.速度分析对（式Ⅱ）求一次导数得：（式Ⅲ）矩阵式：=（
3、Ⅳ）采用高斯消去法可求解（式Ⅳ）可解得角速度ω3,ω4;3.加速度分析把式Ⅳ对时间求导数得矩阵式：=+（式Ⅴ）采用高斯消去法可求解（式Ⅴ）可得角加速度三．程序流程图J=1，N打印结果结束调用高斯消去法子程序求解加速度方程（3）求得α2，α3，α4，α5及α6（或，αC）再求出aEx及aEyB（K）=-DA（K，Ⅱ）ω1（Ⅱ）+DB（K）Ⅱ=1，NDB（K）=DB（K）ω1K=1，Nω(1)=ω2，ω(2)=ω3ω(3)=ω4，ω(4)=ω5调用系数矩阵A子程序，并计算其矩阵DA调用系数矩阵B子程序，并计算其矩阵DB调用高斯消去法子程序求解速度方程（2）求得ω2，ω3
4、，ω4，ω5及ω6（或，vC）再求出vEx及vEyB（J）=B（J）ω1J=1，N调用系数矩阵A子程序，并计算A调用原动件位置参数列阵B子程序，并计算BθI=（I-1）×10゜调用牛顿迭代法子程序求解位置方程（1）求得θ2，θ3，θ4，θ5及θ6（或lBC，s5）并计算xG及yG读入：l1，l2，l2′，l3，l4，l5，l6，xG及yG和θ2，θ3，θ5及θ6（或lAB，lCD，lDE，h，h1，h2及lBD，θ2，θ3，θ4及θ5）的初值，N，ω1，E开始位置分析速度分析加速度分析迭代次数IT=0调用位置方程（1）子程序代入θi的初值，并计算fi停止求得θi值
5、
6、fi
7、≤E调用系数矩阵AIT＞ITmax?调用高斯消去法子程序求解A△θi=fi，求出△θiθi=θi+△θiIT=IT+1YNYN四、计算源程序#include#include#include#definePI3.#defineN4#defineE0.0001#defineT1000voidSolutionangle(double[12],double);/*迭代法求角位移*/voidSolutionspeed(double[N][N],double[N],double[12],double);/*角速度求解*
8、/voidSolutionacceleration(double[N][N],double[N][N],double[N],double[12]);/*角加速度求解*/voidGaussianE(double[N][N],double[N],double[N]);/*高斯消去*/voidFoundmatrixA(double[12],double[N][N]);//创建系数矩阵AvoidFoundmatrixB(double[12],double,double[N]);//创建系数矩阵BvoidFoundmatrixDA(double[12],double[N][N
9、]);//创建矩阵DAvoidFoundmatrixDB(double[12],double,double[N]);//创建矩阵DB//定义全局变量doublel1=180,l3=960,l4=160,h=900,h1=460,h2=110,as1=1.0;//主函数voidmain(){inti,j;FILE*fp;doubleshuju[36][12];doublepsvalue[12],a[N][N],da[N][N],b[N],db[N],ang1;//建立文件，并制表头if((fp=fopen("shuju","w"))==NULL){printf("
