机械优化设计-作业

合肥工业大学《机械优化设计》课程实践研究报告班级：机械设计制造及其自动化2009-5班学号：20090495姓名：王**授课老师：王卫荣日期：2012年5月29日目录机械优化设计研究报告概述21作业一31.1题目〔1〕31.2求解步骤〔1〕31.3结果〔1〕31.4题目〔2〕31.4求解步骤〔2〕31.5结果〔2〕52作业二52.1题目52.2求解步骤52.3结果63作业三63.1题目63.2求解步骤73.3结果84心得体会9机械优化设计研究报告概述优化设计是20世纪60年代初开展起来的一门新学科，它是将最优化原理和计算技术应用于设计领域，为工程设计提供一种重要的科学设计方法。利用这种新的设计方法，人们就可以从众多的设计方案中寻找出最正确设计方案，从而大大提高设计效率和质量。因此优化设计是现代设计理论和方法的一个重要领域，它已广泛应用于各个工业部门。优化方法的应用领域很多，开展也很迅速。今年来开展起来的计算机辅助设计〔CAD〕，在引入优化设计方法后，使得在设计过程中既能够不断选择设计参数并评选出最优化设计方案，又可以加快设计速度，缩短设计周期。在科学技术开展要求机械产品更新周期日益缩短的今天，把优化设计方法与计算机辅助设计结合起来，使设计过程完全自动化，已成为设计方法的一个重要开展趋势。通过本学期的课程学习，我们已经掌握了一些常用的优化方法的原理和计算过程的理论知识。本次实践是稳固学到的理论知识的绝佳方法，通过实践，学生可以将理论知识运用到具体问题当中，培养分析问题和解决问题的能力，同时也能学到优化方法的操作步骤。1作业一1.1题目〔1〕以与课本51页相反的方向证明0.618法。ａｂ1.2求解步骤〔1〕1-ｂ在区间[a,b]内ｂ区间[a,b]两端点具有对称性，即1〔1-1〕其中，是待定常数。由几何关系知，图1-1即〔1-2〕取方程正数解，得1.3结果〔1〕以上证明步骤与课本51页证明方法中所删去的区间不同，但最终结果相同，均可证明出。1.4题目〔2〕编制用0.618法求解函数在区间[0,5]上极小值的程序。1.4求解步骤〔2〕这属于一维搜索的问题，0.618法的程序框图如图1-2所示。根据程序框图的求解思想，利用C语言编制0.618法的程序如下：图1-2#include<stdio.h>voidmain(void){ floata,b,c,d,a1,a2,y1,y2,x,y; c=0.01; d=0.618; a=0.0; b=10.0; printf("%d",a); printf("%d",b); printf("%d",c); a1=b-d*(b-a); a2=a+d*(b-a); /*................赋初值..................*/ y1=(a1-2)*(a1-2)+3; y2=(a2-2)*(a2-2)+3; /*................求初值..................*/ while(((b-a)/b)>=c&&((y2-y1)/y2>=c||-((y2-y1)/y2)>=c)) { if(y1>=y2) { a=a1; a1=a2; y1=y2; a2=a+d(b-a); y2=(a2-2)*(a2-2)+3; } else { b=a2; a2=a1; y2=y1; a1=b-d(b-a); y1=(a1-2)*(a1-2)+3; } } /*................循环体.................*/ x=(a+b)/2; y=(x-2)*(x-2)+3; printf("a*=%f\n",x); printf("y=%f\n",y); /*................最优解.................*/}1.5结果〔2〕将以上程序编译，连接，运行后的结果如下所示。2作业二2.1题目编制能生成n阶单位矩阵的程序，要求n由人工输入。2.1求解步骤思路：利用C语言程序编程，首先是输入n。n阶单位矩阵在C语言程序中就相当于是一个二维数组a[i][j]，当i=j时，a[i][j]=1；当i!=j时，a[i][j]=0。按此思路编程时要定义一个二维数组，并且i=j=n，需要在程序运行时定义i，j，也就是说二维数组[]中的数是个变量，但是在C中要求定义数组时，[]中必须为定值常量或表达式，因此，这种思路缺乏可操作性。现采用以下思路，即利用两个for循环语句，构成嵌套循环格式，并且两个循环次数相同。每行输出n个数后换行，一共输出n列，即能形成n阶矩阵。当行数i等于列数j时，输出1，否那么输出0，即能输出n阶单位矩阵。程序如下：#include<stdio.h>voidmain(void){ inti,j,n; printf("%d",n); scanf("%d",&n); for(i=0;i<=n;i++) { for(j=0;j<=n;j++) { if(i==j) printf("1"); printf("0"); } printf("\n");}2.3结果例如输入n=5，那么运行结果如下：3作业三3.1题目试对图3-1所示主轴进行优化设计，主轴内劲，外力，许用挠度。设计变量数n=3，约束函数个数m=5，收敛精度，初始惩罚因子，惩罚因子缩减系数c=0.2。图3-13.2求解步骤当主轴的材料选定时，其设计方案由四个设计变量决定。即孔径、外径、跨距以及外伸端长度。由于机床主轴内孔常用于通过待加工的棒料，其大小由机床型号决定，不能作为设计变量。故设计变量取为机床主轴优化设计的目标函数那么为式中——材料密度。再确定约束条件。主轴的刚度是一个重要性能指标，其外伸端的挠度不得超过规定值据此建立性能约束。在外力F给定的情况下，y是设计变量x的函数，其值按下式计算式中那么此外，通常还应考虑主轴内最大应力不得超过许用应力。由于机床主轴对刚度要求比拟高，当满足刚度要求时，强度尚有相当富裕，因此应力约束条件可不考虑。边界约束条件为设计变量的取值范围，即综上所述，将所有约束函数规格化，主轴优化设计的数学模型可表示为这里未考虑两个边界约束：和，这是因为无论从小减小伸出端挠度上看，还是从降低主轴重量上看，都要求主轴跨距、伸出端长度往小变化，所以对其上限可以不作限制。这样可以减少一些不必要的约束，有利于优化计算。该问题利用惩罚函数法求解，编制程序段如下。C ====================== SUBROUTINEFFX(N,X,FX)C ====================== DIMENSIONX(N) COMMON/ONE/ITE,KTE,ILI,NPE,NFX,NGR NFX=NFX+1 FX=PE*p*(X(1)+X(3))*(X(2)*X(2)-d*d)/4 RETURN ENDC ========================= SUBROUTINEGGX(N,KG,X,GX)C ========================= DIMENSIONX(N),GX(KG) GX(1)=64X(3)*X(3)*(X(1)+X(3))/3*E*PE*(X(2)^4-d^4)*Y(0)-1 GX(2)=1-X(1)/L(MIN) GX(3)=1-X(2)/D(MIN) GX(4)=X(2)/D(MAX)-1 GX(5)=1-X(3)/A(MIN) RETURN ENDC ========================= SUBROUTINEHHX(N,KH,X,HX)C ========================= DIMENSIONX(N),HX(KH) X(1)=X(1) RETURN END3.3结果最终求得最优解4心得体会通过三次作业时机的实践，增强了我分析问题，解决问题的能力。在此之前，我们遇到的问题大多都是已经数学化了的问题，题目给出式子，然后求解，这样的题目，或许只能考考我们演算技巧而已，并不能考出我们的解决问题的能力。而只有独立的分析实际中的具体问题，然后将其数学化，建立数学模型，再选用适宜的解决方法，求解答案。这样的过程或许才能更好的培养或者考察学生的理论知识水平和实践能力。解决具体的实际问题，思路要灵活。面对不同的问题，可能存在着多种解决方案，选其最优的方案可能效率和效果要更为显著。此外，解决问题的方法尽量简洁，例如在编制n阶矩阵的程序时，我最初的思路是利用二维数组，但由于定义数组元素个数时出现了问题，即假设定义为固定数值，那么n就会有取值范围，适用性就会下降，如果