预览加载中,请您耐心等待几秒...
1/7
2/7
3/7
4/7
5/7
6/7
7/7

在线预览结束,喜欢就下载吧,查找使用更方便

如果您无法下载资料,请参考说明:

1、部分资料下载需要金币,请确保您的账户上有足够的金币

2、已购买过的文档,再次下载不重复扣费

3、资料包下载后请先用软件解压,在使用对应软件打开

机床主轴结构优化设计 指导老师: 姓名: 学号: 机床主轴结构优化设计 一.机械优化设计的一般过程 ①建立优化设计的数学模型 ②选择适当的优化方法 ③编写计算机程序 ④准备必要的初始数据并上机计算 ⑤对计算机求得的结果进行必要的分析 其中,建立优化设计的数学模型是首要的和关键的一步,其基本原则仃: 1、设计变量的选择 在充分了解设计要求的基础上,应根据各设计参数对目标函数的彤响程度认真分析其 主次,尽最减少设计变最的数目,以简化优化设计问题。另外,还应注意设计变量应 当相互独立,否则会使目标函数出现“山脊”或“沟谷”,给优化带来困难。 2、目标函数的确定 常取其中最主要的指标作为目标函数,而其余的指标列为约束条件。 3、约束条件的确定 在选取约束条件时应当避免出现相互矛盾的约朿。因为相互矛盾的约束必然导 致可行域为一空集,使问题的解不存在。另外应当尽量减少不必要的约束。 不必要的约束不仅增加优化设计的计算臺,而且可能使可行域缩小,影响优化 结果。 二、优化实例 机床主轴是机床中重耍零件之般为多支承空心阶梯轴。为了便丁使用材料力学公式进 行结构分析,常将阶梯轴简化成以当量直径表示的等截面轴。下面以两支承主轴为 例,说明优化设计的全过程。 右图所示的是一个己经简化的机床主轴。已知主轴内径 d=30mm,外力215000N,许用挠度yO=O.05mm«主轴材料是 铸钢。密度6Kg/mm3,弹性模*E=210GPa设计变 量数n=3,约束函数个数m=5,设计变量的初值、上下限列于 表8-1初始数据 设计变量XIX2X3-----―一殳 初始值480100120 下限值3006090 上限值650140150 表8-1中。 三.实例分析-数学模型的建立 在设计这根主轴时,有两个重要因素需要考虑。一是主轴的自重:一是主轴伸出 端的挠度。对于普通机床,并不追求过高的加工精度,对机床主轴的优化设计,以选 取主轴的口重最轻为目标,外伸端的挠度是约束条件。 (1) 设计变量的确定 当主轴的材料选定时,其设计方案由四个设计变量决定。即孔径d、外径 D、跨距1及外伸端长度a。由丁•机床主轴内孔常用于通过待加工的棒 料,其大小由机床型号决定,不能作为设计变量。 ]T 故设计变饋取为x=[xxx2X3T=0Da] 目标函数的确定 机床主轴优化设计的目标函数则为 22 f(x)=|irp(xi+x3)(x2一d) 式中,p-材料的密度。 约束条件的确定 主轴刚度是一个重要性能指标,其外伸端的挠度y不得超过规定值y。,据 此建立性能约束 g(x)=y-y0<0 在外力F给定的情况下,y是设计变星x的函数,其值按下式计算 _Fa2(l+a) y3EI 式中,I^^-d4) 则血)曙誅勢—y。兰0 此外,通常还应考虑主轴内最人应力不得超过许用应力。由丁•机床主轴对 刚哎要求比较高,当满足耍求时,强皮尚仃相当富裕,因此应力约束条件可 不考虑。边界约束条件为设计变量的取值范圉,即 】min—amin—a-amax 综上所述,将所有约束函数规格化,主轴优化设计的数学模型可表示为 22) minf(x)丄+x3)(x2一d 2Fbn C/X_64FX3(X1+X3)/ 3—y° 阳乃皿(X产沖 )<0 g2(X <0 g3(x)=l-x2/Dniin g4(X)^/Dniax-1。 <0 g5(x)=l-x3/aniin ,这是因为无论从减小伸出端挠度上 这里未考虑两个边界约束,X]<和X3<amax 看,都要求主轴跨距X丄、伸出端长度X3往小处变化,所以对其上限可以不作限 制。这样可以减少一些不必要的约束,有利于优化计算。 四、优化方法的选择 该实例中,设计变量数n二3,约束函数个数m二5,口有一个约束函数是非线性的, 故属丁•非线性多变量约束优化问题,可以选择MATLAB中的fmincon()函数。因为 fmincon函数是实现多变量约束优化,解决非线性多变量约束问题的一种优化方法。当 然也可以选择遗传算法ga()0 五、程序及编译结果 1.编写一个M文件返回目标函数在点x处的值: functionf=myl(x) p=7.8*10^(-6); d=30; f=0・25*pi*p*(x(l)+x(3))*(x(2)/s2-dzv2); 2.将非线性约束编进一个M文件内 function[c,ceq]=nonlconl(x) F=15000; d=30; E=2•丄*10A5; c=(64*F*x(3)A2*(x(l)+x(3)))/(3*E*pi*(x(2)^4-dA4)*0・05)-1; ceq=[]; 3.(1)调用fmincon()函数: »lb=[3006090]'; %设计变量的下限 ub=[6501