课程作业曲柄摇杆优化设计：宋*学号：2012138229班级：20121057三峡大学机械与动力学院目录1.曲柄摇杆机构优化设计题目要求 (1)2.课题描述 (2)3.数学模型的建立 (3)3.1设计变量的确定 (3)3.2约束条件的建立 (3)3.3目标函数的建立 (6)4.数学模型的建立 (8)5.用matlab优化计算程序及分析讨论 (8)5.1讨论及结果分析 (9)5.2.程序代码过程 (11)6.参考文献 (10)小结 (12)1. 0(3πψψ+=式中0ϕ和0ψ得小于45 ≥[min γγl 1=1l 4 =5。

2.在曲柄输入角从0ϕ到20πϕ+的过程中，使摇杆输出角尽量满足一个给定的函数)(0ϕf 即公式（1）。

对此我将0ϕ到20πϕ+等分为m 分，当然输出角也将对应的分为m 分，然后我将输出角对应的数值与期望函数进行拟合，如果误差降到最小，那么得到的结果将会是优化的解，这是将连续型函数转化为离散型的问题，利用matalab 编程计算，从而求解。

运动模型如图（1）所示图（1）曲柄摇杆机构运动模型图3.数学模型的建立3.1设计变量的确定定义：设计变量是除设计常数之外的基本参数，在优化设计过程中不断地进行修改、调整、一直处于变化的状态，这些基本参数都叫做设计变量。

对于本课题，设计常量为21,l l 长度，分别为1和5。

决定机构部分杆长尺寸32,l l ，以及摇杆按照已知运动规律开始运动时曲柄所处的位置角0ϕ应该列为设计变量即为 X=[]Tx x x 321=[]Tl l 032ϕ由于整个机构的杆长都是按比例来设计的，他们都是1的倍数，按照题目要求曲柄的初始位置为极位角，即0ϕ。

则可以根据曲柄摇杆机构各杆长度关系得到0ϕ和相应的摇杆3l 位置角0ψ的函数，关系式为⎥⎦⎤⎢⎣⎡+-++=42123242210)(2)(arccos l l l l l l l ϕ （2）⎥⎦⎤⎢⎣⎡--+=43232422102)(arccos l l l l l l ψ （3）由已知条件可知21,l l 长度分别为1和5，而根据公式（2）（3）可知，0ϕ0ψ 是由32,l l 的长度来决定，所以32,l l 为独立变量，则可以确定本课题的设计变量 X=[][]TTl l x x 3221=,这是一个二维优化问题。

3.2约束条件的建立定义：如果一个设计满足所有对它提出的要求，成为可行设计；一个可行设计必须满足某些设计限制条件，这些限制条件做为约束条件。

对本题分析可知机构要满足两个约束条件即①杆长条件满足曲柄摇杆机构存在条件②传动角满足最小传动角大于45度︒=≥45][min γγ（1）杆长条件满足曲柄摇杆机构存在条件则有a.最短杆与最长杆长度之和应小于或等于其余两杆之和b.连架杆与机架中至少有一杆是最短杆当最短杆为曲柄时即满足曲柄摇杆存在条件，得到以下约束条件00)(121≤-=≥=x l x g （4）00)(232≤-=≥=x l x g （5）060)(2132413≤--=≤--+=x x l l l l x g （6） 040)(2143124≤--=≤--+=x x l l l l x g （7） 040)(1242135≤--=≤--+=x x l l l l x g （8）（1）传动角满足最小传动角大于45度︒=≥45][min γγ(注：以本机构为例，传动角为32,l l 之间所夹的锐角；机械原理，西工大版） ①当曲柄在),0[π时，如图（2）所示图（2）左极限最小传动角示意图相应的传动角约束条件为[]02)(arccos 180)(32241232206≤-+-+-=γl l l l l l x g （9）②当曲柄在]2,(ππ区间上运动时，相应的传动角约束条件为，如图(3)图（3）右极限最小传动角示意图[]02)(arccos )(3221423227≤⎥⎦⎤⎢⎣⎡--+-=l l l l l l x g γ （10）这是一个具有2个设计变量，7个不等式约束条件的优化设计问题，可以选用约束优化方程成语来计算。

3.3目标函数的建立定义：满足所有约束条件的设计方案是可行设计方案，优化设计的任务就是要对各个设计方案进行比较，从而找出那个最佳的设计方案。

而对设计方案进行优劣比较的标准就是目标函数，或称为评价指标、评价函数。

针对本课题，目标函数可根据已知的运动规律和机构实际运动规律之间的偏差最小作为指标来建立，即取机构的期望输出角)(00ϕψf =和实际输出角)(0ϕψi i f =的平方误差积分最小作为目标函数，表达式为ϕψψπϕϕd i Ei 22)(00-⎰+，而这时一个连续型函数，为了方便计算，我们将这个问题转化为离散型的问题。

把输入角度ϕ取m 个点进行数值计算，它可以化约（4)表达式最小来求解。

2121)(),()(i mi Ei x x f x f ψψ-==∑= （11）Ei ψ--------期望输出角，Ei ψ=)(i E ϕψ；m--------输入角的等分数；i ψ-------实际输出角，由公式（1)可知；由曲柄的运动情况，可以分成三种运动模型，一种是在曲柄在机架之上运动，另一种是曲柄在机架下面运动，最后一种是二者都满足。

我将分别对此讨论，写出相应的目标函数并分析前两种结果对最终结果的影响。

（1)当πϕ≤≤i 0时，如图(4)图（4）曲柄在πϕ≤≤i 0区间模型图实际输出角为i i i βαπψ--= )0(πϕ≤≤i （12）⎥⎦⎤⎢⎣⎡-+=⎥⎦⎤⎢⎣⎡-+=221223222322arccos 2arccos x x x l l l i i i i i ρρρρα (13)⎥⎦⎤⎢⎣⎡+=⎥⎦⎤⎢⎣⎡-+=i i i i i l l l ρρρρβ1024arccos 2arccos 2421242 (14)i i i l l l l ϕϕρcos 1026cos 2412421-=-+= (15)由于我们将2~00πϕϕ+等分为m 分，则实际的输入角i ϕ可以用函数表示出来为mi i 20πϕϕ+= 这里我将输出角的等分数设置成30，则 可以表示出实际输入角的函数为600πϕϕi i += （16）（2）当πϕπ2≤≤i 时，如图(5)图（5）曲柄在πϕπ2≤≤i 区间运动模型图实际输出角为i i i βαπψ+-= )2[πϕπ≤≤i (17)i i i i ϕρβα,,,表达式如（13）（14）（15）（16）所示。

（3）当（1）（2）两种情况都综合考虑进去时，则应该表示为⎩⎨⎧≤≤+-=≤≤--=)2[)0[πϕπβαπψπϕβαπψi i i ii i i i （18） i i i i ϕρβα,,,表达式如（13）（14）（15）（16）所示。

4.数学模型的建立通过上面的分析后，将输入角分成 30 等分（m=30），经过转化为标准形式得到曲柄摇杆机构优化设计标准数学模型为目标函数：min )(),()(2121→-==∑=i mi Ei x x f x f ψψ设计变量：[][]TTl l x x 3221=约束条件：见公式（4）（5）（6）（7）（8）（9）（10)机械优化设计中的问题，大多数属于约束优化问题，此为非线性约束优化问题，运用 MATLAB 优化工具箱的命令函数 fmincon 来处理有约束的非线性多元函数最小化优化问题。

5.用matlab 优化计算程序及分析讨论5.1讨论及结果分析(1)当曲柄在πϕ≤≤i 0运动时猜想 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;-1 1];b = [-6;4;4];% 使用多维约束优化命令fmincon[x,fn] = fmincon(optimfun,x0,a,b,[],[],lb,ub,confun);disp '连杆机构实现函数优化设计最优解'fprintf('连杆相对长度a=%3.4f \n',x(1))fprintf('摇杆相对长度b=%3.4f \n',x(2))fprintf('输出角平方误差之和f*=%3.4f \n',fn)% 计算最优点x*的性能约束函数值g = confun(x);disp '最优点的性能约束函数值'fprintf('最小BCD夹角约束函数值g1*=%3.4f\n',g(1))%fprintf('最大BCD夹角约束函数值g2*=%3.4f\n',g(2))（2）调用目标函数及画图function f=optimfun(x)s=30;qb=1;jj=5;fx=0;fa0=acos(((qb+x(1))^2-x(2)^2+jj^2)/(2*(qb+x(1))*jj));%曲柄初始角pu0=acos(((qb+x(1))^2-x(2)^2-jj^2)/(2*x(2)*jj));%摇杆初始角for i=1:sfai=fa0+0.5*pi*i/s;pu(i)=pu0+2*(fai-fa0)^2/(3*pi);%摇杆期望角ri=sqrt(qb^2+jj^2-2*qb*jj*cos(fai));alfi=acos((ri^2+x(2)^2-x(1)^2)/(2*ri*x(2)));bati=acos((ri^2+jj^2-qb^2)/(2*ri*jj));if fai>0 & fai<=pips(i)=pi-alfi-bati;%elseif fai>pi & fai<=2pi%ps(i)=pi-alfi+bati;endfx=fx+(pu(i)-ps(i))^2;endi=1:30;f=fx;%输出角平分误差之和plot(i,ps(i),'r-.',i,pu(i),'b-*');legend('期望曲线','实际曲线');grid on（3）调用约束条件function [c,ceq]=confun(x)qb=1;jj=5;m=45*pi/180;n=135*pi/180;%c(1)=x(1)^2+x(2)^2-(jj-qb)^2-2*x(1)*x(2)*cos(m);%最小BCD夹角传动角约束c(2)=-x(1)^2-x(2)^2+(jj+qb)^2+2*x(1)*x(2)*cos(n);%最大BCD夹角约束ceq=[];（4）传动角与输入角关系代码%曲柄摇杆机构运动过程中传动角变化曲线图源代码l2=4.0483 ;l3=2.6550;l4=5;a1=acos(((1+l2)^2+l4^2-l3^2)/(2*l4*(1+l2)));a2=pi-acos((l3^2+l4^2-(1+l2)^2)/(2*l4*l3));for i=1:360;fai(i)=a1+(pi/2)*(i/360);rou(i)=sqrt(1+l4^2-2*l4*cos(fai(i)));if rad2deg(acos((l2^2+l3^2-rou(i)^2)/(2*l2*l3)))<=90CDJ(i)=acos((l2^2+l3^2-rou(i)^2)/(2*l2*l3));elseCDJ(i)=pi-acos((l2^2+l3^2-rou(i)^2)/(2*l2*l3));endendx=fai;y=CDJ;plot(x,y)xlabel('曲柄输入角/rad','fontsize',12,'fontname','宋体');ylabel('传动角/rad','fontsize',12,'fontname','宋体');title('给定区间的传动角变化曲线图','fontsize',12,'fontname','宋体');grid on6.参考文献【1】机械原理第七版；西北工业大学机械原理及机械零件教研室编；主编桓做模文杰【2】机械优化设计；工业大学靖民主编7.小结通过对《工程优化与matlab实现》的学习，我初步了解了matlab软件的的使用方法，接触了相关的规准则及设计方法，最主要的，我学习了一种新的思维方式，对我产生了不小的影响。