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

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

ID:42164422

大小:107.25 KB

页数:4页

时间:2019-09-09

优化设计作业2—螺栓优化_第1页
优化设计作业2—螺栓优化_第2页
优化设计作业2—螺栓优化_第3页
优化设计作业2—螺栓优化_第4页
资源描述:

《优化设计作业2—螺栓优化》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库

1、基于MATLAB优化工具箱优化设备的螺栓连接1刖a机械优化设计,就是在给定的环境条件下,在对机械产品的形态、几何尺寸关系以及其他因索的约束范围内,以机械系统的功能、强度和经济性等为优化对象,选取设计变量,建立目标函数和约束条件,并使目标函数获得最优值的一种现代设计方法。机械优化设计广泛的应用于齿轮、轴承、连杆、凸轮、机床等产品的设计屮。处理优化问题,主要有两个步骤:首先要针对工程实际问题,建立数学模型;然后根据数学模型的特点选择优化方法及其计算程序,作必要的简化和加工,用计算机求得最优设计方案。目前,已冇很多成熟的优化方法程序可供选择,

2、但它们各冇自己的特点和适用范围;实际应用屮必须注意优化方法或初始参数选择而带来的收敛性问题等。而MATLAB的优化工具箱则选用最佳方法求解、初始参数输入简单、语法符合工程设计语言要求,编程工作量少,优越性明显。2MATLAB优化工具箱寻优的优点MATLAB语言是MathWorks公司开发的软件产品,是一种面向科学与工程的高级语言,运用它所提供的优化工具箱求解机械优化问题与传统的求解机械优化问题的方法相比有着很人的优越性:(1)利用MATLAB优化工具箱来求解机械优化问题,可以避免Ft?丁我们优化方法选择不当而造成无法得到最优解或所求最优

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

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

5、料为45Cr,抗拉强度=lOOOMpa,屈服强度<rs=320Mpa,拉压疲劳极限a,=330Mpa,许用疲劳安全系数[Sa]=1.7,取残余预紧力F=1.6F,采用铜皮石棉密封垫片,螺栓相对刚度Kc=0.8o从安全、可靠、经济的角度来选择螺栓的个数n和螺栓的直径d。3.2目标函数取螺栓组连接经济成本Cn最小为目标。当螺栓的长度、材料和加工条件一定时,螺栓的总成本与n,d值成正比,故木问题优化设计的目标函数为minf(x)=Cn=nd=x{x2由此可见,设计变量为螺栓个数n和直径d为x={nd]1=(X]x2]J3.3确定约束条件(1)

6、强度约束条件:螺栓在脉动载荷下工作,因此螺栓组连接须满足疲劳强度条件_2cr_j-min、「1—半j其中6为应力幅值。罰为最小应力;K”为疲劳极限综合影响系数,取心=4.4;鸭”为应力折算系数,取儿=0.23.气缸最大载荷P=£兀D;p螺栓最大工作载荷F。=F”+£,nSnax螺栓最小工作载荷=FQ-KaF=Fo-O.Sp/n.

7、)=^-8d=^^-8x2<0一nx,(2)扳手工作空间约束条件:考虑扳手工作空间,螺距间距应大于2d,故扳手工作空间约束条件为q(X)»毗2勺-叫0nx}(3)非约束条件g](x)=-坷<0,g5(X)=-%2SO螺栓连接的优化数学模型综上所述,木问题的数学模型可表达如下设计变量:x=[x^x2Y目标函数:min/(x)=xtx2约束条件:Mg“(x)<0,w=1,2,3,4,54.应用MATLAB工具箱求解functionf=stud_obj(x)f=x(1)*x(2);globalpKsigampsaisigam_lp=8.5*l

8、e6;Ksigam=4.4;psai=0.23;sigam_l=330;D2=240;p=l/4*pi*D2A2*le-6*p;function[czceq]=stud_conl(x)globalpKsig

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。