-《机械优化设计》课程实践报告.doc

-《机械优化设计》课程实践报告.doc合肥工业大学《机械优化设计》课程实践研究报告班级: 机设10-5班学号: 姓名: 授课老师: 王卫荣日期:2013年5月8日目录—・一维搜索程序 +自行选择小型机械设计问题四•课程心得体会一、(a,b)区间上的极小值。用C语言编写程序如下:1•当a=o,b=2兀,y=cosx时。#include<>#include<>#(){floata,b,m[3],y[3],c;scanf(n%f,%f,%f\&a,&b,&c);m[l]=b-N*(b-a);m[2]=a+N*(b・a);y[l]=cos(m[l]);y[2]=cos(m[2]);do{if(y[l]>y[2]){a=m[ll;m[ll=m[2];y[l]=y⑵;m[2]=a+N*(b・a);y[2]=cos(m[l]);}elsem[l]=b-N*(b-a);y[l]=cos(m[l]);j}while(fabs((b-a)/b)>c);m[0]=(a+b)/2;y[0]=cos(m[0]);pFintfCa*=%f\n",m[O]);printf(ny=%f\nM,y[O]);输入a,b,及精度£的值,输出极小点,函数极小值。如下2•当a=0,b=10,y=(x-2)2+3.#incIude<>#incIude<>#(){fIoata,b,m[3],y[3],c;scanf("%f,%f,%f",&a,&b,&c);m[1]=b-N*(b~a);m[2]二a+N*(b-a);y[1]=pow((m[1]-2),2)+3;y[2]二pow((m[2]-2),2)+3;do{if(y[1]>y[2])ta=m[1] ;y[1]=y[2];m[2]二a+N*(b-a);y[2]二pow((m⑴-2),2)+3;}eIse{b二m⑵;m⑵=m[1];y⑵二y[1];m[1]二b-N*(b-a);y[1]二pow((m⑴-2),2)+3;}}whiIe(fabs((b~a)/b)>c);m[0]=(a+b)/2;y[0]二pow((m[0]-2),2)+3;printf(,,a*=%f\n,,,m[0]);printfCy=%f\rT,y[0]);输入a,b,及精度£的值,输出极小点,函数极小值。如下:二、用最简捷的语言生成一个n维单位矩阵,it为输入量。用c语言编写程序如下:#incIude<>voidmain(){inti,j,k=1,m=0,n;scanf("%d",&n);for(i=1;i〈二n;i++){for(jh;j〈=n;j++)if(i二二j)printf("%3d",k);eIseprintf("%3d”,m);printf(H\nn);}}运行结果:当11=7时运行结果:当n=10时一皿•10厶1000000000二0100000000-001000000010001000000100001000000000010000|+自行选择小型机械设计问题。要求:(1)分析优化对象,根据设计问题的要求,选择设计变量,确立约束条件,建立日标函数,建立优化设计的数学模型并编制问题程序;(2)选择适当的优化方法,进行优化计算;(3)分析计算结果,并加以说明。1、连杆机构问题题目:设计一曲柄连杆摇杆机构,要求曲柄11从h从久=%+90°时,摇杆A的转角最佳再现已知的运动规律:妙严必+二(0-©J且h二1,/4=5,久为极位角,其传动角允许在45</<135°范围内变化。1、数学模型的建立该机构的运动简图如上所示,在这个问题中,一直L1二1丄5二5,其他非独立参数可由公式求出,:设计变量为:x=lXixJ=[l2l.],将输入角分成30等分,并用近似公式计算,可得30目标函数的表达式:/&)=工[@-0訂@-加)]/=!约束函数按曲柄存在条件及对传动角的限制来建立,得gl(x)=-xl<0g2(x)=-x2<0g3(x)=6-xl~x2<0g4(x)=xl-x2~4<0g5(x)=x2-xl~4<0g6(x)二xl*xl+x2*x2--16<0g67x)二36-xl*xl-x2*x2T・414xlx2<02、编写子程序如下:[Delphi程序语言-惩罚函数法】unitsumt_fgh;interfaceusesWindows,Messages,SysUtils,Classes