曲柄摇杆机构优化设计说明

得分课程作业曲柄摇杆优化设计：宋*学号：2012138229班级：20121057三峡大学机械与动力学院目录1.曲柄摇杆机构优化设计题目要求.............................................................................12.课题描述.....................................................................................................................23.数学模型的建立.........................................................................................................33.1设计变量的确定...............................................................................................33.2约束条件的建立...............................................................................................33.3目标函数的建立...............................................................................................64.数学模型的建立.........................................................................................................85.用matlab优化计算程序及分析讨论........................................................................85.1讨论及结果分析...............................................................................................95.2.程序代码过程.................................................................................................116.参考文献..................................................................................................................10小结..............................................................................................................................121.曲柄摇杆机构优化设计题目要求要求设计一曲柄摇杆机构，当曲柄由0转到0+90°时，摇杆的输出角实现如下给定的函数关系：200)(32（1）式中0和0分别为对应于摇杆在右极限位置时曲柄和摇杆的位置角，它们是机架杆l4为原线逆时针度量的角度，见图1。要求在该区间的运动过程中的最小传动角不得小于45°，即：45][min通常把曲柄的长度当成单位长度，即l1=1。另外，根据机构在机器中的许可空间，可以适当预选机架杆的长度，现取l4=5。2.课题描述在曲柄输入角从0到20的过程中，使摇杆输出角尽量满足一个给定的函数)(0f即公式（1）。对此我将0到20等分为m分，当然输出角也将对应的分为m分，然后我将输出角对应的数值与期望函数进行拟合，如果误差降到最小，那么得到的结果将会是优化的解，这是将连续型函数转化为离散型的问题，利用matalab编程计算，从而求解。运动模型如图（1）所示图（1）曲柄摇杆机构运动模型图3.数学模型的建立3.1设计变量的确定定义：设计变量是除设计常数之外的基本参数，在优化设计过程中不断地进行修改、调整、一直处于变化的状态，这些基本参数都叫做设计变量。对于本课题，设计常量为21,ll长度，分别为1和5。决定机构部分杆长尺寸32,ll，以及摇杆按照已知运动规律开始运动时曲柄所处的位置角0应该列为设计变量即为X=Txxx321=Tll032由于整个机构的杆长都是按比例来设计的，他们都是1的倍数，按照题目要求曲柄的初始位置为极位角，即0。则可以根据曲柄摇杆机构各杆长度关系得到0和相应的摇杆3l位置角0的函数，关系式为42123242210)(2)(arccoslllllll（2）43232422102)(arccosllllll（3）由已知条件可知21,ll长度分别为1和5，而根据公式（2）（3）可知，00是由32,ll的长度来决定，所以32,ll为独立变量，则可以确定本课题的设计变量X=TTllxx3221,这是一个二维优化问题。3.2约束条件的建立定义：如果一个设计满足所有对它提出的要求，成为可行设计；一个可行设计必须满足某些设计限制条件，这些限制条件做为约束条件。对本题分析可知机构要满足两个约束条件即杆长条件满足曲柄摇杆机构存在条件传动角满足最小传动角大于45度45][min（1）杆长条件满足曲柄摇杆机构存在条件则有a.最短杆与最长杆长度之和应小于或等于其余两杆之和b.连架杆与机架中至少有一杆是最短杆当最短杆为曲柄时即满足曲柄摇杆存在条件，得到以下约束条件00)(121xlxg（4）00)(232xlxg（5）060)(2132413xxllllxg（6）040)(2143124xxllllxg（7）040)(1242135xxllllxg（8）（1）传动角满足最小传动角大于45度45][min(注：以本机构为例，传动角为32,ll之间所夹的锐角；机械原理，西工大版）当曲柄在),0[时，如图（2）所示图（2）左极限最小传动角示意图相应的传动角约束条件为02)(arccos180)(32241232206llllllxg（9）当曲柄在]2,(区间上运动时，相应的传动角约束条件为，如图(3)图（3）右极限最小传动角示意图02)(arccos)(3221423227llllllxg（10）这是一个具有2个设计变量，7个不等式约束条件的优化设计问题，可以选用约束优化方程成语来计算。3.3目标函数的建立定义：满足所有约束条件的设计方案是可行设计方案，优化设计的任务就是要对各个设计方案进行比较，从而找出那个最佳的设计方案。而对设计方案进行优劣比较的标准就是目标函数，或称为评价指标、评价函数。针对本课题，目标函数可根据已知的运动规律和机构实际运动规律之间的偏差最小作为指标来建立，即取机构的期望输出角)(00f和实际输出角)(0iif的平方误差积分最小作为目标函数，表达式为diEi22)(00，而这时一个连续型函数，为了方便计算，我们将这个问题转化为离散型的问题。把输入角度取m个点进行数值计算，它可以化约（4)表达式最小来求解。2121)(),()(imiEixxfxf（11）Ei--------期望输出角，Ei=)(iE；m--------输入角的等分数；i-------实际输出角，由公式（1)可知；由曲柄的运动情况，可以分成三种运动模型，一种是在曲柄在机架之上运动，另一种是曲柄在机架下面运动，最后一种是二者都满足。我将分别对此讨论，写出相应的目标函数并分析前两种结果对最终结果的影响。（1)当i0时，如图(4)图（4）曲柄在i0区间模型图实际输出角为iii)0(i（12）221223222322arccos2arccosxxxllliiiii(13)iiiiilll1024arccos2arccos2421242(14)iiillllcos1026cos2412421(15)由于我们将2~00等分为m分，则实际的输入角i可以用函数表示出来为mii20这里我将输出角的等分数设置成30，则可以表示出实际输入角的函数为600ii（16）（2）当2i时，如图(5)图（5）曲柄在2i区间运动模型图实际输出角为iii)2[i(17)iiii,,,表达式如（13）（14）（15）（16）所示。（3）当（1）（2）两种情况都综合考虑进去时，则应该表示为)2[)0[iiiiiiii（18）iiii,,,表达式如（13）（14）（15）（16）所示。4.数学模型的建立通过上面的分析后，将输入角分成30等分（m=30），经过转化为标准形式得到曲柄摇杆机构优化设计标准数学模型为目标函数：min)(),()(2121imiEixxfxf设计变量：TTllxx3221约束条件：见公式（4）（5）（6）（7）（8）（9）（10)机械优化设计中的问题，大多数属于约束优化问题，此为非线性约束优化问题，运用MATLAB优化工具箱的命令函数fmincon来处理有约束的非线性多元函数最小化优化问题。5.用matlab优化计算程序及分析讨论5.1讨论及结果分析(1)当曲柄在i0运动时猜想1.由于曲柄的实际输出角i的围不完整，会使)0(i结果拉长到整个2区间，从而产生较大偏差。连杆机构实现函数优化设计最优解连杆相对长度a=4.1286摇杆相对长度b=2.3226输出角平方误差之和f*=0.0076最优点的性能约束函数值最小BCD夹角约束函数值g1*=-7.1214最大BCD夹角约束函数值g2*=-0.0000图（6）当)0(i时迭代拟合图结论：通过对比第三种情况，发现二者结果相同，猜想不成立。（2）当曲柄2i运动时猜想2.由于第一种情况下得到的结论对第三种没有影响，可以猜测在区间)2[i不存在，或者和第三种结果一样。结论：Matlab显示结果运行错误。则可以说，第二种情况是不存在的，对结果不产生影响。（3）当（1）（2）两种情况都综合考虑进去时连杆机构实现函数优化设计最优解连杆相对长度a=4.1286摇杆相对长度b=2.3226输出角平方误差之和f*=0.0076最优点的性能约束函数值最小BCD夹角约束函数值g1*=-7.1214最大BCD夹角约束函数值g2*=-0.0000图（7）20i时迭代你拟合图将连杆长度带入计算，则可以得到传动角的变化规律图（8）传动角随输入角变化的规律图结论：经过matlab优化的曲线跟期望曲线存在细微的差别，输出角平方误差之和f*=0.0076，传动角波动围符合要求，所以此优化方程的解符合要求。5.2.程序代码过程（1）优化设计主程序M文件clc;clear;%铰链四杆机构实现函数的优化设计的主程序%调用目标函数optimfun和非线性约束函数confunx0=[6;4];%设计变量的初始值qb=1;jj=5;%设计变量的下界与上界lb=[1;1];ub=[];a=[-1-1;1-1;-11];b=[-6;4;4];%使用