基于MATLAB的高次多项式凸轮型线优化设计8页.docx

机械装备优化设计三级项目基于MATLAB的高次多项式凸轮型线优化设计班级：12级机械装备-2班设计人员（按贡献大小排序）：唐俊杰（12011010093）卫健行（120101010092）王荟博（120101010095） 摘要动力凸轮型线的设计十分重要,以高次多项式凸轮型线为例,在基于丰满系数和磨损设计多目标函数情况下,利用MATLAB及其优化工具箱(optimizationtoolbox)对目标函数数学模型进行优化求解应用MATLAB的优化函数提供的强大计算功能,确定了凸轮型线高次五项式函数中的系数,并能快速找到目标函数的优化值显著提高了型线优化设计的速度和精度,还可根据实际情况灵活地调整权重系数W1和W2的值的大小,计算方便快捷在与传统设计方法比较结果表明,经优化设计,提高了动力凸轮的丰满系数,降低了凸轮型线的磨损关键词:优化设计; MATLAB; 高次多项式; 凸轮型线前言随着工程技术的发展和最优化技术的广泛应用,科学技术领域中需要求解很多优化问题而工程中的最优化问题,绝大多数是有约束的,且多属于非线性规划问题[1]。

内燃机配气系统中凸轮型线的设计就是这一种有约束非线性规划问题由于各种限制条件的复杂性,传统凸轮曲线设计方法难以找到合适的曲线参数,本文讨论此类情况下,利用MATLAB优化工具箱(optimizationtoolbox)解决多目标函数情况下凸轮型线优化设计问题通过MAT-LAB强大的计算功能,确定了凸轮型线高次五项式函数中的系数,利用其优化工具箱寻找了以丰满系数和磨损设计为数学模型的最优值小组分工：唐俊杰 负责MATLAB编程和WORD的制作王荟博 负责WORD的制作和PPT的制作卫健行 负责WORD的制作，查阅资料一、 高次五项式凸轮型线凸轮曲线采用高次多项式型线,多项式项位向右影响渐弱,通常高次多项式取五项至七项第六项以后对动力性能指标的影响力已很小,随着高次多项式项数的增加,凸轮丰满系数减小,加速度的绝对值变大,凸轮机构工作性能下降而对于磨损设计的关键位置桃尖而言,因α=αB,β=(1-α/αB)= 0,高次多项式项数的增加不起作用,磨损与h″α有关,h″α=pCpω 2/α 2B[2]故对一般内燃机配气机构的凸轮,在以丰满系数和磨损为目标函数来说,为设计计算方便,取五项即可满足设计要求。

高次五项式作为气门升程曲线可表达为h(α)= C0+Cpβp+Cqβq+Crβr+Csβs式中:h(α)表示气门升程,简写为hα;α表示凸轮转角,将基本段始点取作α= 0;β=(1-α/αB),αB为基本段半包角;C0,Cp,Cq,Cr,Cs表示待定系数;p,q,r,s表示幂指数;取p= 2,q= 2n,r= 2n+ 2m,s=2n+ 4m,式中一般取n为3至20之间的实数;m为1至20之间的实数二、 凸轮型线优化设计目标函数凸轮型线对系统的工作质量有着重要的影响,所以配气凸轮的优化设计,在很大程度上能提高内燃机的工作性能[4]为使配气机构具备良好的充、排气性能要求凸轮型线具有较大的丰满系数ξ=∫2αB0h(α)dα2αBhmax式中:hmax为凸轮的最大升程;ξ表示凸轮型线的丰满系数,其物理意义为高次曲线和基本工作段包角所围的面积与最大升程和基本工作段包角所围面积之比它反映配气机构的气体通过能力,丰满系数大,进排气效率高,动力性能好第一目标函数是建立以丰满系数为最大目标的函数对高次五项式方程推导,得丰满系数为[5]ξ=1C0C0+Cp1+p+Cq1+q+Cr1+r+Cs1+s(2)第二目标函数的建立要求凸轮型线磨损量为最小,对第二目标函数磨损设计的求解,实际上演化为对高次五项式2阶求导,h″α=pCpω 2/α 2B,其物理含义为:表示与高次曲线形状、初始条件等因素相关的在桃尖处的垂直加速度大小。

第二目标函数是建立以磨损量为最小目标的函数凸轮型线多目标函数合成为minF(x)= W1(-ξ)+W2(h″α) minF(x)=W1-1C0C0+Cp1+p+Cq1+q+Cr1+r+Cs1+s+W2(pCpω 2α 2B) 式中:W1、W2分别为第一、第二目标函数权重;系数C0、Cp、Cq、Cr、Cs为含有n、m变量的系数所以x=[x1,x2]T=[n,m]T,目标函数F(x)为含有n、m变量的函数C0= hmax(最大升程) (5)Cp=-C0qrs+v(qr+qs+rs-q-r-s+ 1)(q-p)(r-p)(s-p) (6)Cq=-C0prs+v(pr+ps+rs-p-r-s+ 1)(p-q)(r-q)(s-q) (7)Cr=-C0pqs+v(pq+ps+qs-p-q-s+ 1)(p-r)(q-r)(s-r) (8)Cs=-C0pqr+v(pq+pr+qr-p-q-r+ 1)(p-s)(q-s)(r-s) (9)式中:v=vgαBω,vg为阀门落座速度;αB为工作段半包角;ω为凸轮角速度。

三、 优化函数约束条件a) 边界约束s.t.　g1(X)= [3-x1]≤ 0g2(X)= [x1- 20]≤ 0g3(X)= [1-x2]≤ 0g4(X)= [x2- 20]≤ 0b) 最小曲率半径约束g5(X)= [rmin]-rmin　　　　[rmin]-r0-C0-amin(ω/57.3)2≤ 0式中:amin= pCpω 2/α 2B;rmin为凸轮外形最小曲率半径,mm;[rmin]许用最小曲率半径c) 最大加速度约束本应对最大加速度进行约束,但由于在以磨损设计作为第二目标函数时,其数学推导结果,实际上是对加速度进行间接限制,故不再在约束条件中对加速度重复约束经分析以上式(4)的数学模型,可知是一个多目标约束非线性的最小化问题,应用MATLABOpti-mizationToolbox中fmincon函数[1],编写m代码,实现上述优化数学模型的求解,通过给定条件参数,由MATLAB计算出幂指数p、q、r、s,及相应的系数Cp、Cq、Cr、Cs,从而实现凸轮高次五项式型线幂指数和系数的确定d) MATLAB程序计算实例凸轮机构的设计参数如下:凸轮基圆半径r0=21mm,最大升程hmax= 13mm,转速n= 1800r/min,工作段半包角αB= 62℃,阀门落座速度vg= 0.2m/s,权重系数W1、W2取等值0.5,权重系数也可根据具体情况进行权重调整。

为了避免某一单个目标函数值过大而影响到整个函数的取值,在编程时引入K值概念,计算时将各目标函数分别除以K1、K2值K值一般取该函数约束区域内的最大值e) 约束条件程序f) 程序运行及计算结果在MATLAB程序的命令状态下输入两个目标函数的权重W1,W2的值,并可根据实际情况变化很方便地进行调整,此处以W1、W2都为0.5为例,计算结果如下:从结果可看出,优化设计的凸轮型线丰满系数ξ高于常规设计;反映桃尖磨损设计的h″α的值也优于常规设计设计时可以在MATLAB中调整权重系数W1和W2值的大小,定义两目标函数的重要程度如设置W1等于1,W2等于0,则转化为只考虑丰满系数为目标函数的单目标设计的特殊情况四、 结束语函数凸轮中的高次多项式动力凸轮被普遍采用,但用一般设计方法设计函数凸轮得到的丰满系数较小,而采用优化设计,且将丰满系数作为目标函数,则可得到理想的设计结果[8]本文基于丰满系数和磨损设计双目标函数情况下,应用MATLAB的优化函数提供的强大计算功能,能快速找到目标函数的优化值,最优解为n = 20,m = 11.3707,丰满系数ξ= 0.6329,最小曲率半径rmin= 9.2820mm,磨损值h″α=-267.4890m/s2,相应的幂指数p= 2,q= 40,r= 62.7414,s= 85.4829,各项系数为Cp=- 14.4696,Cq= 3.5386,Cr=- 2.8226,Cs=0.7537。

满足设计要求,显著提高了型线优化设计的速度和精度,还可根据实际情况灵活地调整权重系数W1和W2的值的大小,计算方便快捷五、 感悟在这次《优化设计》三级项目中，我们小组三个人各自分工协作，互相帮助，共同完成了这次项目并且算出来了结果，对优化设计的数学建模和在机械设计中的重要作用有了更新的认识书山有路勤为径，学海无涯苦作舟只有我们不断的学习新的知识，我们才能更好解决我们之后所遇到的问题参考文献[ 1]　褚洪生,杜增吉,阎金华等编著.MATLAB 7.2优化设计实例指导教程[M].北京:机械工业出版社,2006[ 2]　何振俊,周俊冬.配气凸轮轮廓曲线设计的多目标函数建模[J].机械设计,2008,(5):34～ 36[ 3]　孙明,付军.柴油机高次方配气凸轮型线设计[J].青岛远洋船员学院学报,2005,(4):1～ 3[ 4]　李兵,顾德裕.基于MATLAB的配气凸轮优化设计[J].山东内燃机,2006,(1):31～ 33[ 5]　李金林,谭建荣.改进模拟退火算法在配气凸轮机构优化设计中的应用[J].工程图学学报,2006,(1):19～ 23[ 6]　游国忠,赵晓丹,王忠,梅德清.遗传算法在内燃机凸轮优化设计中的应用[J].江苏理工大学学报(自然科学版), 2000,(3):31～ 34[ 7]　管荣法.凸轮与凸轮机械[M].北京:国防工业出版社,1993[ 8]　詹樟松,杨正军,刘兴春.高次多项式动力凸轮优化设计及MATLAB算法实现[J].内燃机,2004,(1):4～ 6。