《基于Matlab实现曲柄摇杆机构的运动设计.docx》由会员分享，可在线阅读，更多相关《基于Matlab实现曲柄摇杆机构的运动设计.docx（7页珍藏版）》请在三一文库上搜索。
1、颜幽以上交上学实用数值方法(Matlab)小论文题目：基于Matlab实现曲柄摇杆机构的运动设计小组成员姓名：毛晓雯 学号:201202070607班级： 机自6班2014-2015(1学期提交日期：2014年12月29日基于Matlab实现曲柄摇杆机构的运动设计1问题提出与数学模型的建立曲柄摇杆机构是钱链四杆机构中的一种，在实际工程应用中，该机构应用广泛，如缝纫机踏板机构、搅拌机机构等。现要求设计一曲柄摇杆机构，能同时实现以下几个要素：1）为提高机构的急回特性，极位夹角应尽可能大（017 ）;2）为改善机构的传力性能，当该机构曲柄与连杆重叠共线时，最大压力角尽可能小（0 max 55 ）；3
2、）该摇杆摆角=60 。1.1 设计变量的确定设11,12,13,14分别为该四杆机构的杆长，考虑计算的方便性，令1 1, 1 X1 , ； X 2 ,l1l1l11411Tx 3 ,于是设计变量为 X= x 1 x 2 x 3图1-1曲柄摇杆机构简图1.2 目标函数的建立当曲柄摇杆机构的各杆长度确定后，该机构的摇杆摆角、 最大压力角及极位夹角都会确定下来，即该机构的各项性能也能确定下来。这里，将摇杆摆角=60这个目标处理为无限接近 60这个目标值，定义为一目标函数，求之同 60之差的绝对值的最小值。由图1-1及设计要求，可列出该设计的分目标函数。因为=C2ADQAD,由余弦定理可得极位夹角的目
3、标函数：fi(x) I I 222222x 11 x 3 x 2x 11 x 3 x 2arccos arccos2x 3 x 112x 3 x 11对于曲柄摇杆机构来说，当主动曲柄与机架共线的两位置之一处，压力角达到最大。在这里,1-2所示)。故最大压力角我们很容易知道当 &点位于A和D之间的这种情况，压力角最大(如图f2 x的目标函数:max222x 1 x 2 x 3 1一 arccos22x 1 x 2摇杆的实际最大摆角为=ADC3 - ADC2 ,故222222x 2 x 3 x 11x 2 x 3 x 1 1=arccos arccos2x 2 x 32x 2 x 3所以实际最大摆
4、角与理论最大摆角60间偏差的目标函数:图1-2曲柄摇杆压力角达到最大位置这里，根据设计要求，极位夹角 尽可能大为好。为了能实现用 Matlab编程求解，做一下处理, 定义f1 x 。这样就实现了 f1 x , f2 x , f3 x min ,满足Matlab优化工具箱对目标函数的统一 要求。下面对各目标函数的范围做进一步明确：1)为保证极位夹角尽可能大，确定极位夹角不小于10 ,即1017 ,于是转换成弧度，有0.2967 f1(x)0.1745;2)由设计要素的条件 0 max 55得:0 f2(x) 0.95993;编辑版word3 )为保证摇杆的摆角 尽可能接近60 ，使 1 f 3
5、x 1 ，转换成弧度，有：0 f2(x) 0.01745。1.3 约束条件的确定曲柄摇杆机构要满足曲柄存在的条件。其中最长杆与最短杆的长度之和应小于等于其它两杆长度之和；最短杆为连架杆(此处为li)o为了便于计算，这里取li 1,因此有:g1x1x 10g2x1x 20g3x1x(3)0g4 x1x 3x 1x 20g5 x1x 1x 2x 302使用Matlab进行优化求解这里曲柄摇杆机构的数学模型属于一般非线性多目标规划问题。其基本模型1min F (x) weight.y goalx,ys.t. c(x) 0Ceq(x) 0Ax bAeq.x beq调用 Matlab 优化工具箱中非线性
6、规划求解函数fgoalattain 来求解。其命令的基本格式：X,fval,attainfactor,exitflag=fgoalattain(ObjFun,x0,goal,weight,A,b,Aeq,beq,lb,ub); 其中各个参数的含义如下： ObjFun 为定义求解的多目标函数； x0 为变量初值； A ， b 为约束函数中线性不等式约束的参数；Aeq,beq为约束函数中线性等式约束的参数；c(x),ceq(x)i匀为非线性函数组成的向量；lb,ub分别是变量 X 的下限和上限值； goal 是目标函数边界组成的向量； Weight 是目标函数的权重值。对于判断求解是否有效可通过a
7、ttainfactor 和 exitflag 两个参数来确定： 当 exitflag0 时， 表示目标值没有溢出goal。本例求解最优解，分两个步骤进行：1)先建立M文件ObjFun.m,定义各个目标函数输入： function f=Objfun(x)f1=-abs(acos(x(1) -1)A2+x(3)A2-x(2F2)/(2*x(3)*(x(1)- 1)-acos(x(1)+1)A2+x(3F2-x(2)A2)/(2*x(3)*(x(1)-1);% 极位夹角的目标函数f2=pi/2-acos(x(1)A2+x(2)A2-(x(3)-1)A2)/(2*x(1)*x(2); % 最大压力角的
8、目标函数f3=abs(acos(x(2)A2+x(3)A2-(x(1)+1)A2)/(2*x(2)*x(3)-acos(x(2)A2+ x(3)A2-(x(1)-1)A2)/(2*x(2)*x(3)-pi/3); %摇杆摆角的目标函数f=f1;f2;f3;2)根据非线性多目标规划求解函数fgoalattain的形式，建立主程序输入：goal=-0.1745 0.95993 0.01745; % 给定目标函数的上边界向量weight=0.3333 0.3333 0.3333; % 个目标的权重值x0=3 3 3;A=-1 0 0;0 -1 0;0 0 -1;1 -1 -1;-1 -1 1;b=-
9、1; -1; -1; -1; -1;Aeq=;beq=;lb=1 1 1;ub=;x,fval=fgoalattain(Objfun,x0,goal,weight,A,b,Aeq,beq,lb,ub)运算结果：x=2.95132.14332.8873fval= 0.27840.87890.0000对于多目标优化设计，各个目标之间往往存在互相制约和矛盾的情况，因此其很难准确给出最优结果。在使用Matlab优化工具箱求解时，只有在明确了各个目标函数的权重值前提下，才有可能。曲柄摇杆机构的优化设计含有三个目标函数，在满足不同目标函数权重值的前提下，理论上满足条件的解应该是无穷多个。综合考虑，上述程序
10、是令各目标函数权重值Weight=0.3333,即在各目标函数同等重要的情况下，并且在杆长初始值为3 3 3的前提下，得出的优化解。为了得到较为可靠的优化解，带入不同的杆长初始值，将计算结果统计成表2-1如下：表2-1计算结果记录A组值第二组值第三组值第四组值杆长初始值2345l11111l21.85392.95133.97034.1310l32.02332.14332.23952.1749l42.37912.88733.55283.812011.252915.951115.068812.2326max48.615550.357253.319450.609460.000060.000060.0
11、00060.0000从上表数据可以看出，满足该多目标优化解的结果很多。对比各组数据，综合考虑，第二组解是相对最优解，其中：li 1, 12 2.9513, I3 2.1433, I4 2.8873,这组数据使得机构有良好的急 回特性和传力特性。3结论在进行曲柄摇杆机构及其他类型机构的多目标优化设计时，应对各个目标函数做统一边界处理，明确各目标函数在整个数学模型中的权重值和所设计的机械性能特性，使用Matlab优化工具箱，进行程序编写，找出满足条件的各组解，在进行性能相对比较分析，从而找到最优解。这种方法方便、 快捷、高效并且准确，为较好地解决连杆机构多目标优化设计问题提供了一个途径。编辑版word参考文献1阳明盛，熊西文.Matlab基础及数学软件M.大连：大连理工大学出版社，20122 何正风 . Matlab 在数学方面的应用 M. 北京：清华大学出版社， 2012（此文档部分内容来源于网络，如有侵权请告知删除，文档可自行编辑修改内容，供参 考，感谢您的配合和支持）编辑版 word