《机械优化设计》大作业Word文件下载.doc

1、机械优化工程实例 连杆机构11 自选机构16课程实践心得 20附列程序文本 21进步，努力，建议25一、0.618的证明、一维搜索程序作业关于的证明黄金分割法要求插入点，的位置相对于区间两端具有对称性，即其中为待定常数。此外，黄金分割法还要求在保留下来的区间内再插入一点所形成的区间新三段，与原来的区间三段具有相同的比例分布。黄金分割法还要求在保留下来的区间内再插一点所形成的区间新三段，与原来的区间三段有相同的比例分布。b1211-图一设原区间的长度为1，如图一所示，保留下来的区间长度为，区间缩短率为。为了保持相同的分别比例。插入新点应在位置上，在原区间的1位置应相当于在保留区间的位置。故有：

2、a23（1-）2图二解得编写0.618的程序，并计算下列问题程序框图（1）采用MATLAB进行编程%fun.m:黄金分割法求极值点%输入数据% a 搜索区间下限% b 搜索区间上限% e 精度%输出数据：% x 极小值点%其他常量：% c1,c2,c3 区间% r 黄金分割比例0.618%定义函数function x=fun（a,b,e） r=0.618; c1=b-r*（b-a）;c2=a+r*（b-a）; y1=f（c1）;y2=f（c2）; while （abs（b-a）/b）e）&（abs（y2-y1）/y2）e） if y1=y2 a=c1;c1=c2;y1=y2; c2=a+r*（

3、b-a）; else b=c2;c2=c1;y2=y1; c1=b-r*（b-a）;y1=f（c1）; end end x=0.5*（a+b）;end 对f函数的确立10 function y=f（x） y=（x-2）2+3;11 function y=f（x） y=cosx如果要计算y=（x-2）2+3;的黄金分割法，则需要将图所示的f脚本中的函数写成如所示如果要计算y=cosx；黄金分割法，则需要将图所示的f脚本中的函数写成如所示fun函数表示对matlab的主程序语言。26函数NO.2运算结果：函数NO.1运算结果：（2）采用C语言进行编程#include #define K 0.618

4、double f（double）; /*函数值计算函数声明*/void main（void）double a,b,size;double a1,a2;int I;printf（“请输入区间两端点（端点值应大于0）:”）;scanf（“%lf,%lf”,&a,&b）; /*输入端点值*/printf（“请输入精度:scanf（“%lf”,&size）; /*输入精度*/printf（“区间为（%lf,%lf）,精度为%lfn”,a,b,size）;printf（“序号t a1tt a2tt f（a1）tt f（a2）n”）;for（i=0;isize） /*用精度控制循环次数*/i+; a1=b

5、-K*（b-a）; /*按0.618法插入两点*/ a2=a+K*（b-a）; printf（“%2d:t%fta%ft%ft%fn”,I,a1,a2,f（a1）,f（a2）; /*输出每次计算后a1,a2,f（a1）,f（a2）的值*/if（f（a1）=f（a2）a=a1;elseb=a2;printf（“所求极小值点为:x=%lft极小值f（x）=%fn”,a,f（b）;double f（double x） /*函数值计算函数*/double f;f=（x-2）*（x-2）+3;return（f）;对于y=cosx，须在程序中加一个#include（math.h）头程序，以示我要调用函数语

6、句。再将倒数第三行函数改为y=cos（x），再次编译运算即可。运行如图所示我们可以看出，大约在x=2处取到极值。与理论相符。我们可以看出，大约在处取到极值。二、用简单的语句写一个单位矩阵采用matlab构造：Matlab构造的矩阵最是简便，因为本身就有一个单位矩阵的函数。由于matlab是基于C语言而设立的一个数学运用软件，所以他的集成度非常高，而且具有很好的开放性，于是我们得知eye函数，将其构造。Eg：若要构造一个3维的单位矩阵，则输入eye（3）若要构造一个n为的单位矩阵，这输入eye（n）用C语言构造程序如图所示： int main（） int I,j,n; while（scanf（“

7、%d”,&n） = 1） for（I = 0; I n; +i） for（j = 0; j +j） printf（“%d “,I = j）; printf（“n”）; return 0;三、机械优化设计工程实例 连杆机构问题（1）连杆机构问题描述图 1 机构简图设计一曲柄连杆摇杆机构，要求曲柄从时，摇杆的转角最佳再现已知的运动规律：且=1，=5，为极位角，其传动角允许在范围内变化。（2）数学模型的建立设计变量：这里有两个独立参数和。因此设计变量为目标函数：将输入角分成30等分，并用近似公式计算，可得目标函数的表达式约束条件：GX（1）=-X（1）0GX（2）=-X（2） 0GX（3）=-（X（

8、1）+X（2）+6.00GX（4）=-（X（2）+4.0）+X（1） 0GX（5）=-（4.0+X（1）+X（2） 0GX（6）=-（1.4142*X（1）*X（2）-X（1）*2-X（2）*2）-16.00GX（7）=-（X（1）*2+X（2）*2+1.4142*X（1）*X（2）+36.00（3）程序编制C=SUBROUTINE FFX（N,X,FX）DIMENSION X（N）COMMON /ONE/I1,I2,I3,I4,NFX,I6NFX=NFX+1P0=ACOS（1.0+X（1）*2-X（2）*2+25.0）/（10.0*（1.0+X（1）Q0=ACOS（1.0+X（1）*2-X（

9、2）*2-25.0）/（10.0*X（2）T=90.0*3.1415926/（180.0*30.0）FX=0.0DO 10 K=0,30 PI=P0+K*TQE=Q0+2.0*（PI-P0）*2/（3.0*3.1415926）D=SQRT（26.0-10.0*COS（PI）AL=ACOS（D*D+X（2）*X（2）-X（1）*X（1）/（2.0*D*X（2）BT=ACOS（D*D+24.0）/（10.0*D） IF（PI.GE.0.0 .AND. PI.LT.3.1415926） THEN QI=3.1415926-AL-BT ELSE QI=3.1415926-AL+BT END IF IF

10、（K.NE.0 .OR. K.NE.30）THEN FX=FX+（QI-QE）*2*T FX=FX+（QI-QE）*2*T/2.012 CONTINUERETURNENDC=SUBROUTINE GGX（N,KG,X,GX）DIMENSION X（N）,GX（KG）GX（1）=-X（1）GX（2）=-X（2）GX（3）=-（X（1）+X（2）+6.0GX（4）=-（X（2）+4.0）+X（1）GX（5）=-（4.0+X（1）+X（2）GX（6）=-（1.4142*X（1）*X（2）-X（1）*2-X（2）*2）-16.0GX（7）=-（X（1）*2+X（2）*2+1.4142*X（1）*X（2

11、）+36.0SUBROUTINE HHX（N,KH,X,HX）DIMENSION X（N）,HX（KH）X（1）=X（1）输入数据如下：2,7,04.3,3.22.,0.2,0.01,0.00001,0.000010,1,03.5,2.8,5.5,4.5（4） 可执行程序的生成13 将程序文本保存成文件t001.txt，存放在OPT_EXAMFORTRAN0-fortran目录下。14 使用DOS操作界面，OPT_EXAMFORTRAN0-fortran目录下，输入命令 asumt.for+t001.txt t001.for 回车。生成t001.for文件。15 for1 t001; 回车；f

12、or2 回车；link t001; 回车，即生成可执行程序t001.exe。从t001.dat文件中加载数据，运行t001.exe。t001t001.rtf 回车得到结果。操作过程如图：操作步骤（5）结果与分析 结果显示如下：= PRIMARY DATA = N= 2 KG= 7 KH= 0 X : .4300000E+01 .3200000E+01 FX: .7847605E-02 GX: -.4300000E+01 -.3200000E+01 -.1500000E+01 -.2900000E+01 -.5100000E+01 -.6729390E+01 -.1218939E+02 HX:

13、PEN = .3974390E+01 R = .2000000E+01 C = .2000000E+00 T0= .1000000E-01 EPS1= .1000000E-04 EPS2= .1000000E-04= ITERATION COMPUTE =* IRC = 1 R = .2000000E+01 PEN= .3974390E+01* IRC = 2 R = .4000000E+00 PEN= .3119361E+01* IRC = 3 R = .8000001E-01 PEN= .6407908E+00* IRC = 4 R = .1600000E-01 PEN= .1439420

14、E+00* IRC = 5 R = .3200000E-02 PEN= .3733451E-01* IRC = 6 R = .6400000E-03 PEN= .8986510E-02* IRC = 7 R = .1280000E-03 PEN= .2346609E-02* IRC = 8 R = .2560000E-04 PEN= .9042178E-03* IRC = 9 R = .5120001E-05 PEN= .5396702E-03* IRC =10 R = .1024000E-05 PEN= .4410887E-03* IRC =11 R = .2048000E-06 PEN=

15、.4128061E-03* IRC =12 R = .4096001E-07 PEN= .4071001E-03* IRC =13 R = .8192003E-08 PEN= .4014154E-03* IRC =14 R = .1638401E-08 PEN= .3992769E-03* IRC =15 R = .3276801E-09 PEN= .3987240E-03* IRC =16 R = .6553602E-10 PEN= .3984505E-03* IRC =17 R = .1310720E-10 PEN= .3983411E-03* IRC =18 R = .2621441E-

16、11 PEN= .3982889E-03= OPTIMUM SOLUTION = IRC= 18 ITE= 35 ILI= 94 NPE= 447 NFX= 482 NGR= 0 R= .2621441E-11 PEN= .3982696E-03 .4135127E+01 .2315381E+01 .3982623E-03 -.4135127E+01 -.2315381E+01 -.4505072E+00 -.2180254E+01 -.5819746E+01 -.7079845E+01 -.3633157E-03Stop - Program terminated.从结果中得知当X1=4.13

17、5127，X2=2.315381时，目标函数取得最小值, 摇杆的转角最佳再现已知的运动规律：。 自选结构模型设计某带式输送机减速器的高速级齿轮传动。已知高速级输入功率P1 = 10kW,小齿轮转速n1 =960 r /min,传动比i = 3. 2。齿轮材料和热处理:大齿轮45号钢（调质）硬度为217255HBS,小齿轮40Cr（调质）硬度为241286HBS,工作寿命15 年,假设每年工作300天,两班制,带式输送机工作平稳,转向不变。常规设计方案采用直齿圆柱齿轮: m=2.5, z1=30, d=1。解：（1）设计变量，单级直齿圆柱齿轮传动的中心距 :齿宽：将m,作为设计变量，即： =2

18、）目标函数根据多目标优化的线性加权法建立体积最小的目标函数: f （ x） =1f1 （ x） +2f2 （ x）=1+2其中:1 ,2 是加权系数,且1 +2 = 1,分别根据设计时径向和轴向安装位置的要求设定;取1 = 1表示要求中心距最小,取2 = 1则表示要求齿宽最小。3 ）约束条件模数的限制:对于传递动力的齿轮,通常要求模数不少于1. 5-2,得约束条件: （2）小齿轮齿数的限制:小齿轮齿数应不大于产生根切的最小齿数17 ,得约束条件: （3）齿宽系数的限制:由于min max ,约束条件为: （4）齿面接触强度的限制，根据公式并查表得约束条件：（5）齿根弯曲强度的限制，根据公式查表

19、得约束条件：4 建立数学优化模型高速级齿轮传动多目标优化设计的数学模型为:（1 取0.6，2取0.4）Fun（x）=min1 +2 =5编写程序并运行结果目标函数M文件：function f=zhwm（x）f=0.6*2.1*x（1）*x（2）+0.4*x（1）*x（2）*x（3）;约束函数M文件：function c ceq=zhwy（x）c（1）=1.04*107-2.916*105*（x（1）*x（2）3*x（3）;c（2）=1.04*107-8.95*106*（x（1）*x（2）3*x（3）;c（3）=1.51*106-303.57*x（1）3*x（2）2*x（3）;c（4）=1.42

20、*106-2445.92*x（1）3*x（2）2*x（3）;ceq=;优化函数M文件：x0=2 32 1;lb=1.5 17 0.7;ub=2 inf 1.15;u=;x,fval=fmincon（zhwm,x0,lb,ub,zhwy）约束函数目标函数经过Matlab优化并圆整后的齿轮参数如下：设计参数模数m/mm齿数齿宽系数中心距/mm普通设计2.530157.5MATLAB优化1.791127.43771.149984.5254圆整后参数1.75117.6经过计算，最小体积为87.15。四、课程学习心得从懒懒散散的寒假到匆匆忙忙的考试周，仿佛只是在弹指一挥间。然而我们就在这段时间内，又学了几门课程，而这些课程中，最有趣，最好玩的当然就是机械优化设计了机械优化设计是一门理论性非常强的一门课，刚开始的时候不得要领，但是在王卫荣老师的带领下，逐渐登堂入室，渐渐有所了解。首先，这门课所依托的最强武器便是数学。数学即是方法，思维，和逻辑。其次所依托的便是算法程序语言。归根结底的还是数学，因为数学赋予了算法以逻辑的力量，分析的过程，以及所能解决的一切问题。还好我的数学基础还算可以，逻辑思维还算通顺，再加上老师所给予的方法