机械优化设计作业

1、合肥工业大学机械优化设计课程实践研究报告班 级： 机设12-6班 学 号： 2012216281 姓 名： 丁雷鸣 授课老师： 王卫荣 日 期： 2015年 11月 10 日 目 录1、 =0.618的证明12、 一维搜索程序作业1（1）例1程序文本1（2）例1输出结果截图2（1）例2程序文本2（2）例2输出结果截图3三、单位矩阵程序作业4（1）程序文本4（2）输出结果截图4四、连杆机构问题6（1）目标函数6（2）约束条件7（3）选择方法7（4）程序文本7（5）数据输入截图8 （6）输出结果9 5、 自行选择小型机械设计问题或其他工程优化问题10（1）设计变量10（2）目标函数10（3）约束条

2、件10（4）程序文本10（5）数据输入截图11 （6）输出数据11 六、机械优化设计课程实践心得体会 131一、=0.618的证明在实际计算中，最常用的一维搜索方法是黄金分割法。黄金分割法是建立在区间消去法原理基础上的试探方法，即在搜索区间内适当插入两点，。并且计算其函数值。黄金分割法要求插入点，的位置相对于区间两端点具有对称性，即、其中为待定常数。除对称要求外，黄金分割法还要求保留下来的区间内再再插入一点，所形成的区间新三段与原来的区间三段具有相同的比例分布。设原区间长度为1，保留下来的区间长度为，区间缩短率为。为了保持想相同的比例分布，新插入点应该在位置上，在原区间的1-位置应该相当于在保

3、留区间的位置。故有取方程正数解，得二、一维搜索程序作业例1、a=0,b=,f(x)=cosx （1）例1程序文本#includeincludevoid main ()float a,b,c=0.618,aa3,y3,d;scanf(“%f,%f,%f”,&a,&b,&d):aa1=b-c*(b-a);aa2=a+c*(b-a);y1=cos(aa1);y2=cos(aa2);doif（y1y2）a=aa1;aa1=aa2;y1=y2;aa2=a+c*(b-a);elseb=aa2;aa2=aa1;y2=y1;aa1=b-c*(b-a);y1=cos(aa1);while(fabs(b-a)/b

4、d);aa0=(a+b)/2;y0=cos(aa0);printf(“a=%fn”,aa0);printf(“y=%fn”,y0);（2）例1输出结果截图：输入a=0,b=2,精度d=0.000001,输出极小值点和函数极小值如下：例2、a=0,b=10,f(x)=(x-2)+3（3）例2、程序文本#include#includevoid main() float a,b,c=0.618,aa3,y3,d; scanf(%f,%f,%f,&a,&b,&d); aa1=b-c*(b-a); aa2=a+c*(b-a); y1=(aa1-2)*(aa1-2)+3; y2=(aa2-2)*(aa2-

5、2)+3; do if(y1y2) a=aa1;aa1=aa2;y1=y2; aa2=a+c*(b-a); y2=(aa2-2)*(aa2-2)+3; else b=aa2;aa2=aa1;y2=y1; aa1=b-c*(b-a); y1=(aa1-2)*(aa1-2)+3; while(fabs(b-a)/b)d); aa0=(a+b)/2; y0=(aa0-2)*(aa0-2)+3; printf(a*=%fn,aa0); printf(y=%fn,y0); （4）例2输出结果截图：输入a=0,b=10,精度d=0.000001，输入极小值点和函数极小值如下：三、单位矩阵程序作业作业：编写

6、生成单位矩阵的程序。要求：通用、输出美观、语言少为佳。（1）程序文本#includevoid main (void)#define m 500int i,j,n,amm; printf(please input a number);scanf(%d,&n);for(i=1;i=n;i+) for(j=1;j=n;j+) if(i=j) aij=1; else aij=0; printf(%d,aij); printf(n); （2）输出结果截图：当n=4时，输出结果如下：当n=12时，输出结果如下：四、连杆机构问题设计一曲柄摇杆机构，要求曲柄从转到=+时，摇杆的转角最佳再现已知的运动规律：=+

7、且已知=1，=5，为极位角,其传动角允许在40范围内变化。设计变量该机构的运动简图如上图所示。在这个问题中，已知=1，=5且和不是独立参数，它们可由下式求出：=arccos=arccos所以还问题只有两个独立参数和，因此设计变量为x=（1）目标函数将输入角分成30等分，并用近似公式计算，可得目标函数的表达式f(x)=式中当时的机构实际输出角，其计算公式为=式中=arccos=arccos= arccos=arccos=为当时的理想输出角，其值由下式计算=+（2）约束条件 平面铰链四杆机构常用的约束条件有曲柄存在条件和传动角条件。由此得约束条件为、（3）选择方法采用惩罚函数法进行计算。（4）程序

8、文本procedure ffx; /目标函数var p0,q0,p,ri,a,b,q,qi,k :real; i:integer ;begin with form1.sumt do begin p0:=arccos(sqr(1.0+x1)-x2*x2+25.0)/(10.0*(1.0+x1); q0:=arccos(sqr(1.0+x1)-x2*x2-25.0)/(10.0*x2); k:=90.0/30.0*(3.1415926/180.0); fx:=0;for i:=0 to 30 do;beginp:=i*k+p0;qi:=q0+2.0*sqr(p-p0)/(3.0*3.1415926

9、);ri:=sqrt(26.0-10.0*cos(p);a:=arccos(ri*ri+x2*x2-x1*x1)/(2.0*ri*x2);b:=arccos(ri*ri+24.0)/(10.0*ri);q:=3.14159-a-b;fx:=fx+sqr(q-qi)*k;end;end;end;procedure ggx; /约束函数begin with form1.sumt do begin gx1:=-x1; gx2:=-x2; gx3:=6.0-x1-x2; gx4:=x1-x2-4.0; gx5:=x2-x1-4.0; gx6:=x1*x1+x2*x2-1.414*x1*x2-16; g

10、x7:=36-x1*x1-x2*x2-1.414*x1*x2; end;end;procedure hhx;begin with form1.sumt do beging hx1:=hx1; end;end;end.（5）数据输入截图输入初始点为（5,5），精度为0.001 （6）输出结果五、自行选择小型机械设计问题或其他工程优化问题某车间生产甲、乙两种产品。生产甲种产品每件需要材料9kg、3个工时、4kw电，可获利60元。生产乙种产品每件需材料4kg、10个工时、5kw电，可获利120元。若每天能供应材料360kg，有300个工时，能供200kw电，问每天生产甲、乙两种产品各多少件，才能够获

11、得最大的利润。（1）设计变量设每天生产的甲、乙两种产品分别为、件。因此设计变量为x=（2）目标函数此问题的数学模型为f(，)=60+120所以目标函数的表达式为minf(x)= -60-120（3）约束条件依题意得约束条件为：0、9+4-360（材料约束）3+10-300（工时约束）4+5-200（电力约束）（4)程序文本procedure ffx; /目标函数begin with form1.hfgd do begin nfx := nfx+1; fx:=-60*x1-120*x2; end;end;procedure ggx; /约束函数begin with form1.hfgd do b

12、egin gx1:=9*x1+4*x2-360; gx2:=3*x1+10*x2-300; gx3:=4*x1+5*x2-200; gx4:=-x1; gx5:=-x2; end;end;（5）输入数据截图(6)输出数据 常用优化方法 约束随机法 一、初始数据= 设计变量个数 n = 2 不等式约束个数 kg = 5 - 随机方向个数 nsr = 6 - 初始步长 t0 = 0.001 收敛精度 eps = 0.0001 - 设计变量初始点 x0:x1=10x2=10 - 设计变量下界 bl:bl1=0bl2=0 - 设计变量上界 bu:bu1=100bu2=100 - 初始点目标函数值 f(

13、x0)= -1800 - 初始点处的不等约束函数值 g(x0):gx1= -2.300000e+02gx2= -1.700000e+02gx3= -1.100000e+02gx4= -1.000000e+01gx5= -1.000000e+01- 二、计算过程_数据= 设计变量迭代点 x: 迭代次数 ite = 1x1= 1.000018e+01x2= 1.000088e+01 - 目标函数值 f(x)= -3811.30860498082 - 三、优化结果_数据= 迭代次数 ite = 4 目标函数计算次数 ifx = 125 - 设计变量最优点 x*:x1= 2.113596e+01x2= 2.309058e+01 - 最优值 f(x*)= -4039.05576008917 - 最优点处的不等约束函数值 g(x*):gx1= -7.740408e+01gx2= -5.686514e+00gx3= -1.395986e-04gx4= -2.113725e+01gx5= -2.309017e+01- st