合工大机械优化设计课程实践报告

1、文档来源为：从网络收集整理.word版本可编辑.欢迎下载支持合肥工业大学机械优化设计课程实践班级：研究报告机械设计制造及其自动化12-3班学号:姓名：授课教师：王卫荣日期：2015 年11 月14 日一、一维搜索程序作业 31. 入=0.618的证明32. 编写0.618法程序并计算 4二、单位矩阵程序作业 6三、连杆机构问题和自选工程优化问题 71. 连杆机构问题72. 自选工程优化问题 .14四、课程实践心得体会 18一、一维搜索程序作业1. 入=0.618的证明黄金分割法，又称作0.618法，适用于a, b区间上的任何单谷函数求极小值问题。黄金分割法是建立在区间消去法原理基础上的试探方法

2、，即在搜索区间a,b内适当插入两点 a、a，并计算其函数值。a、a将区间分成三段。应用函数的 单谷性质，通过函数值大小的比较，删去其中一段，使搜索区间得以缩短。然后再在 保留下来的区间上做同样的位置，如此迭代下去，使搜索区间无限缩小，从而得到极小点的数值近似值。黄金分割法要求插入点 a、a的位置相对于区间a，b两端点具有对称性，即图1-1黄金分割法ai = b -入（b -a ）a2 = a + 入（b _a ）（3-1）其中，入为待定常数。下面证明入=0.618。除对称性要求外，黄金分割法还要求保留下来的区间内再插入一点所形成的区间新三段，与原来区间的三段具有相同的比例分布。设原有区间a,

3、b长度为1如图1-1所示，保留下来的区间a，b长度为入，区间缩短率为入。为了保持相同的比 例分布，新插入点 a应在入（1 -入）位置上，a在元区间的1 -入位置应相当于在 保留区间的入2位置。故有1 -入二入2即入2入-1 = 0取方程正数解得Vs1k =超 0.6182若保留下来的区间为a，b，根据插入点的对称性，也能推得同样的入的值。2. 编写0.618法程序并计算 (1)0.618法程序:#i nclude #in clude float m=0.618;float fun (float t) float y;y=cos(t);return y;mai n() float a,b,eps

4、;prin tf(mi n=”)；输入函数下限%输入函数上限%输入精度%scan f(%f,&a);%prin tf(max=);scan f(%f,&b);%float t1,t2,t,f1,f2,mi n; prin tf(eps=);scan f(%f,&eps);%文档来源为 :从网络收集整理 .word 版本可编辑 .欢迎下载支持while(b-a)/b=eps) t1=a+(1-m)*(b-a);t2=a+m*(b-a);f1=fun(t1);f2=fun(t2);if(f1=f2) a=t1;t1=t2;f1=f2;t2=a+m*(b-a);f2=fun(t2); else b=

5、t2;t2=t1;f2=f1;t1=a+(1-m)*(b-a);f1=fun(t1);t=(a+b)/2;min=fun(t);printf( 最优点 t=%fn,t); % 输出最优点 t% printf( 最优值 f= %fn,min); % 输出最优值 f%(2) 计算 :1.a= 0, b = 2 n, Y= cos(t)2. a=0,b=10, y= (t-2)*(t-2)+3二、单位矩阵程序作业编写生成单位矩阵的程序： 程序文本#include void main( void ) int a100100;int N,i,j;printf( 请输入所要输出矩阵的阶数(最多 100 阶

6、) :); scanf(%d,&N);printf( 输出的矩阵阶数为 %dn,N);printf( N );/* 制作表头 */for (i=0;iN;i+)/* 横行序号 */printf(%3d,i+1);printf(n);for (i=0;iN+1;i+)printf(-);printf(n);for (i=0;i100;i+)for (j=0;j100;j+)if (i=j)aij=1;elseaij=0;for (i=0;iN;i+)printf (%2d:,i+1);for (j=0;jN;j+) printf(%3d,aij);printf(n);/*/*分割线 */数组赋值

7、 */* 输出所需数组 */* 纵列序号 */结果显示从键盘输入 9，显示 9 阶单位矩阵，结果如下三、连杆机构问题和自选工程优化问题1. 连杆机构问题 ：问题描述：图1现优化一曲柄连杆机构，如图1所示，已知曲柄长度L1为44mm机架长度L4 为220mm,要求当曲柄的转角在 0, 0+n 12 时，对应的摇杆的输出角为屮i , 且两者满足对应函数关系屮i= W0+ ( 0 - i )2 , 0和屮0分别对应于四连杆在 初始位置时曲柄和摇杆的位置角。要求机构传动角的范围是n /4 , 3n /4 ，优化该问题使得从动件的一系列实际输出角与期望实现函数屮=f ()的值的平方偏差之和最小。模型建立

8、文档来源为：从网络收集整理.word版本可编辑.欢迎下载支持1、设计变量曲柄摇杆机构按照原动件和从动件的对应关系可知其有 5个独立参数，对于图1 分别为曲柄长度L1，连杆长度L2，摇杆长度L3,机架长度L4,曲柄初始角 0和摇 杆的初始角屮0,由于L1和L4已知，且由图1的几何关系知：所以 0和屮0已不再是独立参数，而是杆长的函数。经上分析独立变量只有 L2和L3。因此，选择连杆长度L2和摇杆长度L3作为设计变量。即：X = L2 L3T = X1 X2T2、目标函数由上面图2和图3中机构的几何关系可得如下的运动规律:S为角度区间的分段数；W si为机构的实际输出角，计算式为:si =咒卩“匕

9、乞卩乞2小riJ耳十g 巒根据图中的角度关系求得：a 耳）*洛所以根据本机构设计问题，以机构实际输出角屮 si与理论输出角屮i的平方偏差最小原则来建立目标函数。优化目标函数表达式：3、约束条件根据已知条件，该机构的约束条件有两方面：一是传动运动过程中的传动角应大于 45度且小于 135 度；二是保证四杆机构满足曲柄存在条件。( 1)保证传动角约束图4图5根据图 4 和图 5 中机构处于最大传动角和最小传动角时的连杆几何关系， 由余弦 定理知将L2=X1, L3=X2 L1=44, L4=220代入上面两式得 (2)、曲柄存在条件由机械原理的知识可知，曲柄存在的条件为： 将已知的杆长和设计变量代

10、入上述条件得：经分析上述杆长条件不起约束作用，实际起作用的约束只g1(X)和g2(X),所以最终的数学模型为：优化设计1 、程序运行结果Command Windo：wWorkspace:根据程序运行结果可知， exitflag 的值为 1，说明目标函数收敛到局部最优解， 优化效果较为理想， 此时目标函数 f 的值为 0.0121rad2 ，连杆的长度为 181.5602mm， 摇杆的长度为 102.4337mm。2 、结果分析当曲柄在 0, 0+ 范围内转动时，摇杆输出角与期望实现函数屮 =f ()的平方偏差值之和最小为 0.0121rad2 ，最优点位于约束条件 g2(X)=0 上。2. 自

11、选工程优化问题(1)问题描述：设计一个压缩圆柱螺旋弹簧，要求其质量最小。弹簧材料为65Mn最大工作载荷为Fmax 40N，最小工作载荷为0,载荷变化频率fr 25Hz，弹簧寿命为10 不发生共振现象，要求g4(x) 0.356 106x21x32 375 0 弹簧旋绕比的限制g5(x)怡帚4.0 0h，弹 簧钢丝直径d的取值范围为1-4mm中径D2的取值范围为10-30mm 工作圈数n不应 小于4.5圈，弹簧旋绕比C不应小于4，弹簧一端固定，一端自由，工作温度为50C， 弹簧变形量不小于10mm(2)数学模型的建立设计变量本题优化目标是使弹簧质量最小，圆柱螺旋弹簧的质量可以表示为：式中，-弹簧

12、材料的密度，对于钢材=7.8 10对d,n,D2的限制且应取标准值，即 1.0，1.2，1.6，2.0，2.5，3.0，3.5，4.0mm 等。kg/mm3;n-工作圈数；n2-死圈数，常取n2=1.5-2.5，现取压=2;D2-弹簧中径(mm ;d -弹簧钢丝直径(mrh；将d，n，D2作为设计变量，即泣-打-J.=_ ,V目标函数将已知参数代入公式，进行整理后得到问题的目标函数为约束条件根据弹簧性能和结构上的要求，可写出问题的约束条件：(1)强度条件gO350 163.0x12.86x0.860(2)刚度条件g2(x)0.4 10 2x14x2x310.0 0(3)稳定性条件g3(x)24

13、33.7x3 （X2 1.5）为 0.44 10 捲 X2X3 0建立优化模型由上可知， 该压缩圆柱螺旋弹簧的优化设计是一个三维的约束优化问题， 其数学 模型为：min f (x) M 0.192457 10 4(x2 2) x12 x3 (计算时系数无影响可舍去) (3)程序编制并运行结果从上面的分析，以重量最轻为目标的汽门弹簧的优化设计问题共有 3个设计变量， 11个约束条件。按优化方法程序的规定，编写数学模型的程序如下：FX=0.0000192457*(X(2)+2)*X(1)*2*X(3)RETURNENDSUBROUTINE GGX(N,KG,X,GX)DIMENSION X(N),

14、GX(KG)GX(1)=350*X(1)*(2.86)-163*X(3)*(0.86)GX(2)=0.004*X(2)*X(3)*3-10*X(1)*4GX(3)=3.7*X(3)*X(1)*4-(X(2)+1.5)*X(1)*5-0.0044*X(2)*X(3)*3GX(4)=356000*X(1)-375*X(2)*X(3)*2GX(5)=X(3)-4*X(1)GX(6)=X(1)-1GX(7)=4-X(1)GX(8)=X(2)-4.5GX(9)=50-X(2)GX(10)=X(3)-10GX(11)=30-X(3)RETURNEND利用惩罚函数法(SUM法)计算，得到的计算结果如下：优化

15、结果一：优化结果二：两种不同结果比较：优化结果比较d/mmn/圈数D2/mm第一次优化.176910E+01.574960E+01.162097E+02第二次优化.176910E+01.574960E+01.162097E+02不同初值所得到的结果相同。优化的结果为：x=1.7691,5.7496,16.2097，圆整为标准值为：x=1.6,6,16， 经检验合格，相应弹簧的最轻质量 6.3g。四、课程实践心得体会尊敬的王老师：由于我最近忙于考研，只剩下四十多天的时间，在时间紧张的情况下，作业完成 得不是很好，请您原谅。但是您上的课我有去听，虽然有时上课写写作业，但是您上 课的激情还是感染到了我，正如您所说，我也觉得您是一位把教学当作自己的一种爱 好的好老师。除了大一高数老师经常板书以外，您也是我遇到的板书最多的一位老师， 哈哈。谈到机械优化设计，您让我对优化设计问题有了新的深刻的认识，以前觉得高数深奥难懂，没什么实际用途，但是上了您的优化设计课，发现用处还是非常大的，而 且通过您结合工程实际问题的讲解，更加让我对数