优化设计作业2—螺栓优化

1、课程编号 12S15C0208 课程名称 化工设备优化设计 学期 2016年春学位层次 硕士 适合专业 化工过程机械 共 4 页第 3 页基于MATLAB优化工具箱优化设备的螺栓连接1 前言机械优化设计，就是在给定的环境条件下，在对机械产品的形态、几何尺寸关系以及其他因素的约束范围内，以机械系统的功能、强度和经济性等为优化对象，选取设计变量，建立目标函数和约束条件，并使目标函数获得最优值的一种现代设计方法。机械优化设计广泛的应用于齿轮、轴承、连杆、凸轮、机床等产品的设计中。处理优化问题，主要有两个步骤：首先要针对工程实际问题，建立数学模型；然后根据数学模型的特点选择优化方法及其计算程序，作必要

2、的简化和加工，用计算机求得最优设计方案。 目前，已有很多成熟的优化方法程序可供选择，但它们各有自己的特点和适用范围；实际应用中必须注意优化方法或初始参数选择而带来的收敛性问题等。而MATLAB的优化工具箱则选用最佳方法求解、初始参数输入简单、语法符合工程设计语言要求，编程工作量少，优越性明显。2 MATLAB 优化工具箱寻优的优点MATLAB 语言是Math Works 公司开发的软件产品,是一种面向科学与工程的高级语言,运用它所提供的优化工具箱求解机械优化问题与传统的求解机械优化问题的方法相比有着很大的优越性:(1) 利用MATLAB 优化工具箱来求解机械优化问题,可以避免由于我们

3、优化方法选择不当而造成无法得到最优解或所求最优解并不理想的情况。在这个工具箱中,对每一种函数每一步的求解都是通过选择一种最佳方法来进行的。例如在求解约束优化问题时,我们一般的方法总是将其变换为较容易的子系统问题,然后求解,这种方法效率不高,在工具箱中此方法已被集中于对KT方程进行求解的方法所取代。在求解KT 方程时,选择的是序列二次规划(SQP) 方法,并通过BFGS算法来更新Hessian 矩阵。(2) 利用MATLAB 优化工具箱来求解最优化问题,可以节省编制优化程序的时间。在用此工具箱解优化问题时,我们只需利用文件编辑器来编写目标函数及约束函数的M文件,然后调用相应的优化函数,系统即可自

4、动运行求出最优解,对于无约束的优化问题只需在命令窗口中输入相应的目标函数及初值,直接调用相应函数即可。(3)MATLAB 工具箱还提供给我们各种形式的输出结果。如我们将options 设置为1 时,就可以以表格的形式输出优化结果,其中包含了迭代次数、各个迭代阶段的函数值等。(4) 利用option 参数还可以实现选择主要算法、选择搜索方向算法、控制有限微分梯度计算中变量x 扰动的水平等功能。3. 算法举例3.1 问题的提出图示为一压气机气缸与缸盖连接的示意图。已知气缸外径D1=400mm,内径D2=240mm，缸内工作压力p=8.5Mpa，螺栓材料为45Cr，抗拉强度，屈服强度，拉压疲劳极限，

5、许用疲劳安全系数，取残余预紧力，采用铜皮石棉密封垫片，螺栓相对刚度。从安全、可靠、经济的角度来选择螺栓的个数n和螺栓的直径d。3.2 目标函数 取螺栓组连接经济成本最小为目标。当螺栓的长度、材料和加工条件一定时，螺栓的总成本与n，d值成正比，故本问题优化设计的目标函数为由此可见，设计变量为螺栓个数n和直径d为 3.3 确定约束条件 （1）强度约束条件：螺栓在脉动载荷下工作，因此螺栓组连接须满足疲劳强度条件其中为应力幅值;为最小应力；为疲劳极限综合影响系数，取=4.4；为应力折算系 数，取=0.23.气缸最大载荷螺栓最大工作载荷， 螺栓最小工作载荷, 螺栓应力幅值 对于普通螺纹，小径，于是疲劳强

6、度约束条件为 （2）密封约束条件：考虑密封安全，螺栓间距应小于8d，故密封约束条件为 （3）扳手工作空间约束条件：考虑扳手工作空间，螺距间距应大于2d，故扳手工作空间约束条件为 （4）非约束条件 螺栓连接的优化数学模型 综上所述，本问题的数学模型可表达如下 设计变量： 目标函数：约束条件：4. 应用MATLAB工具箱求解 function f=stud_obj(x) f=x(1)*x(2); global p Ksigam psai sigam_1 p=8.5*le6; Ksigam=4.4;psai=0.23;sigam_1=330;D2=240; p=1/4*pi*D22*le-6*p;

7、function c,ceq=stud_conl(x) global p Ksigam psai sigam_1 p,x F0=(1.6*p+p)/x(1); d1=0.85*x(2); A=1/4*pi*d12; F1=F0-0.8*p/x(1); sigma_max=F0/A; sigma_min=F1/A; sigma_1=(sigma_max-sigma_min)/2; Sa=(2*sigma_1+(Ksigma-psai)*sigma_min)/(Ksigma+psai)*(2*sigma_a+sigma_min); c=1.7-Sa; 400*pi/x(1)-8*x(2); -400*pi/x(1)+2*x(2); 16-x(1); -x(2); ceq=; x0=7,20; x,feval=fmincon(stud_obj,x0,stud_conl) 运行结果： x = 16.00 28.84 fval = 461.39根据实际问题的意义取整、标准化：n=16，d=30，经验证n和d的取值满足约束条件。5. 结语应用MATLAB优化工具箱进行优化设计问题求解,编程