机械优化设计方案三个案例.docx
- 文档编号:13766641
- 上传时间:2023-06-17
- 格式:DOCX
- 页数:17
- 大小:31.25KB
机械优化设计方案三个案例.docx
《机械优化设计方案三个案例.docx》由会员分享,可在线阅读,更多相关《机械优化设计方案三个案例.docx(17页珍藏版)》请在冰点文库上搜索。
机械优化设计方案三个案例
机械优化设计案例1
1.题目
对一对单级圆柱齿轮减速器,以体积最小为目标进行优化设计。
2.已知条件
已知数输入功p=58kw,输入转速n=1000r/min,齿数比1?
]=550Mpa,许用弯用应力[曲应力u=5,齿轮的许H?
]=400Mpa。
[F3.建立优化模型
3.1问题分析及设计变量的确定
由已知条件得求在满足零件刚度和强度条件下,使减速器体积最小的各项设计参数。
由于齿轮和轴的尺寸(即壳体内的零件)是决定减速器体积的依据,故可按它们的体积之和最小的原则建立目标函数。
单机圆柱齿轮减速器的齿轮和轴的体积可近似的表示为:
222222?
?
)?
0.25(b?
c)(.25Db(d?
d?
dv?
0.25)b(d?
?
d)?
02gzz1g122222222?
?
?
?
d?
)?
0.257l(d8?
dddc?
2112zzzz022222222?
?
)10m(mzu?
d?
b.25?
[m0zb?
d.b?
m8zbub0?
1112zz12222]3228dd6d)?
d?
l?
05bd.?
005b(mzu?
10m?
1..22zz2zz2z121式中符号意义由结构图给出,其计算公式为d?
mz,d?
mz2112D?
umz?
10m12gd?
1.6d,d?
0.25(umz?
10m?
1.6d)2z2g210zc?
0.2b
由上式知,齿数比给定之后,体积取决于b、z、m、l、d和z11
d六个参数,则设计变量可取为z2TT]ddbzmxxxxx]l?
[xx?
[23145z61z213.2目标函数为
222222f(x)?
0.785398(4.75xxx?
85xxx?
85xx?
0.92xx?
xx?
5231116233112222220.8xxxx?
1.6xxx?
xx?
xx?
28x?
32x)?
min6646213316545约束条件的建立3.3.
zz?
17?
,得1)为避免发生根切,应有min0?
?
17?
xg(x)21b?
?
?
?
?
?
?
maxmind的最大值为齿宽系数2)齿宽应满足和,dmaxmin?
?
,,得和最小值,一般取=1.4=0.9maxming(x)?
0.9?
x(xx)?
03212g(x)?
x(xx)?
1.4?
031323)动力传递的齿轮模数应大于2mm,得
g(x)?
2?
x?
0344)为了限制大齿轮的直径不至过大,小齿轮的直径不能大于d,得max1g(x)?
xx?
300?
0352d?
d?
d5)齿轮轴直径的范围:
得maxzminzz0?
100?
xxg()?
560?
x150?
g(x)?
570?
x?
g(x)?
130680200?
x)?
x?
g(69l按结构关系,应距离满足条件:
撑6)轴的支?
b?
2?
?
0.5d?
l=20),得(可取2zminming(x)?
x?
0.5x?
x?
40?
041610)齿轮的接触应力和弯曲应力应不大于许用值,得7.
0550?
xxx)?
1468250g(x)?
(1231170980?
x)?
?
400g(42?
?
2212)x?
0.854?
10xxxx(0.169?
0.6666?
102223170980g(x?
?
400?
)4?
22213)x?
10?
0.?
xxx(0.2824?
0.17710394x23221?
?
][8)齿轮轴的最大挠度,得不大于许用值max440?
.003xxx(xx)?
0g(x)?
117.044521443?
?
][9)齿轮轴的弯曲应力,得不大于许用值ww6x?
102.8512124?
5.5?
2.4?
100g(x)?
()?
153xxx3526x1085?
12.2124?
5.5?
10?
0?
g(x)()?
6163xxx3624.优化方法的选择
由于该问题有6个设计变量,16个约束条件的优化设计问题,采用传统的优化设计方法比较繁琐,比较复杂,所以选用Matlab优化工具箱中的fmincon函数来求解此非线性优化问题,避免了较为繁重的计算过程。
5.数学模型的求解
5.1.1将已知及数据代入上式,该优化设计的数学优化模型表示为:
2222minf(x)?
0.785398(4.75xxx?
85xxx?
85xx?
311132322220.92xx?
xx?
0.8xxxx?
1.6xxx?
xx5631616154312222?
xx?
28x?
32x)6465Subjectto:
0?
17?
x(gx)?
210xx)?
0.9?
x(g(x)?
32210.4?
(xx)?
1g(x)?
x33120?
x?
g(x)?
2340?
300?
x)?
xxg(3520?
x?
(x)?
100g560150?
x)?
x?
g(570x?
x)?
130?
g(680?
?
200(x)g?
x690?
?
405x?
xg(x)?
x?
0.410160550?
xx)?
?
g(x)1468250(x131127098?
400)?
?
0g(x4?
?
22212)10x.854(0.169?
0.6666?
10?
x?
0xxx231227098?
400?
0xg()?
4?
22213)x394?
.2824?
0177?
1010x?
0.xxx(0.2122344g(x)?
117.04x(xxx)?
0.003x?
044514326x?
1012.852124?
5.5?
0)?
2.4?
10(gx)?
(153xxx3526x10.85?
122124?
5.10?
5?
0)g(x?
()?
6163xxx326
5.1.2运用Matlab优化工具箱对数学模型进行程序求解
首先在Matlab优化工具箱中编写目标函数的M文件
myfun.m,返回x处的函数值f:
functionf=myfun(x)
f=0.785398*(4.75*x
(1)*x
(2)^2*x(3)^2+85*x
(1)*x
(2)*x(3)^2-85*x
(1)*x(3)^2+0.92*x
(1)*x(6)^2-x
(1)*x(5)^2+0.8*x
(1)*x
(2)*x(3)*x(6)-1.6*x
(1)*x(3)*x(6)+x(4)*x(5)^2+x(4)*x(6)^2+28*x(5)^2+32*x(6)^2)
由于约束条件中有非线性约束,故需要编写一个描述非线性约束条件的M文件mycon.m:
function[c,ceq]=myobj(x)
c=[17-x
(2)。
0.9-x
(1)/(x
(2)*x(3))。
x
(1)/(x
(2)*x(3))-1.4。
2-。
x(6)-200。
130-x(6)。
x(5)-150。
100-x(5)。
x
(2)*x(3)-300。
x(3).
x
(1)+0.5*x(6)-x(4)-40。
1486250/(x
(2)*x(3)*sqrt(x
(1)))-550。
7098/(x
(1)*x
(2)*x(3)^2*(0.169+0.006666*x
(2)-0.0000854*x
(2)^2))-400。
7098/(x
(1)*x
(2)*x(3)^2*(0.2824+0.00177*x
(2)-0.0000394*x
(2)^2))-400。
117.04*x(4)^4/(x
(2)*x(3)*x(5)^4)-0.003*x(4)。
(1/(x(5)^3))*sqrt((2850000*x(4)/(x
(2)*x(3)))^2+2.4*10^12)-5.5。
(1/(x(6)^3))*sqrt((2850000*x(4)/(x
(2)*x(3)))^2+6*10^13)-5.5]。
ceq=[]。
最后在commandwindow里输入:
x0=[230。
21。
8。
420。
120。
160]。
%给定初始值
[x,fval,exitflag,output]=fmincon(@myfun,x0,[],[],[],[],[],[],@myobj,output)%调用优化过程
5.1.3最优解以及结果分析
运行结果如下图所示:
由图可知,优化后的最终结果为
x=[123.356599.85171.7561147.3157150.4904129.5096]
7f(x)=2.36e*10由于齿轮模数应为标准值,齿数必须为整数,其它参数也要进行圆整,所以最优解不能直接采用,按设计规范,经标准化和圆整后:
x=[1241002148150130]
7f(x)=6.16*106.结果对比分析
73,而优化后mm大约为6.32×10V若按初始值减速器的体积73,优化结果比初始值体积减少为:
10mm的体积V则为6.16×77)×100%10/6.32×=101Δν=-(6.16×2.5%
所以优化后的体积比未优化前减少了2.5%,说明优化结果相对比较成功。
7.学习心得体会
学习机械优化设计课程的心得体会
通过将近一学期的学习,对这门课有了初步的了解和认识,学期伊始,浏览全书,发现全是纯理论知识,觉得这门课会很枯燥,但是又回过头来想想,作为21世纪的大学生,要使自己适应社会需求,首先在做任何事之前都应该有正确的态度看待问题,把这些想法作为促使自己进步的动力,再去学习课本知识,效果应该很不一样,有了想法就付诸行动,随着对课本内容的学习跟老师的讲解,发现并不是像自己在学期初想的那样困难,特别是在老师介绍了一些与机械优化设计相关的计算机语言和计算机软件后,真正体会到科学优化设计的强大跟简洁明了,与传统优化设计方法相比较,大大提高了设计效率和质量。
传统设计方法常在调查分析的基础上,参照同类产品通过估算,经验类比或实验来确定初始设计方案,如不能满足指标要.
求,则进行反复分析计算—性能检验—参数修改,到满足设计指标要求为止。
整个传统设计过程就是人工凑试和定性分析比较的过程,是被动地重复分析产品性能,不是主动设计产品参数。
按照传统设计方法做出的设计方案,有改进余地,但不是最佳设计方案。
而现代化设计工作是借助电子计算机,,应用一些精确度较高的力学数值分析方法,优化软件进行分析计算,找最优设计方案,实现理论设计代替经验设计,用精确计算代替近似计算,用优化设计代替一般的安全寿命可行性设计。
在进行程序求解的过程中,因为是初学Matlab软件,对很多问题的关键点不能够掌握,非线性约束如何书写,上、下限如何选择,函数格式如何书写,变量未定义等等或大或小的问题,但是在一步步排除错误、重新编写程序的过程中,渐渐的对Mtalab熟悉起来,懂得了一些优化方法的简单计算过程和原理,省去了繁琐复杂的优化计算过程
在学完课程之后,反思自己在学习过程中的得失,深深体会到,不论在人生的哪个阶段,都要对自己负责,做任何事都要耐心,细致,“千里之行,始于足下”,学会在物欲横流的社会大潮中,坚持踏踏实实走好人生的每一步。
8.参考文献
[1]孙靖民,梁迎春.机械优化设计.北京:
机械工业出版社,2006.
[2]濮良贵,纪名刚.机械设计.8版.北京:
高等教育出版社,2006.
[3]孙桓,陈作模,葛文杰.机械原理.7版.北京:
高等教育出版社,2006.
[4]李涛,贺勇军,刘志俭.MATLAB工具箱应用指南—应用数学篇[M].北京:
电子工业出版社,2000.
机械优化设计案例2
复杂刀具优化设计数学模型的建立及算法改进
摘要:
目的建立复杂刀具优化的数学模型,提高优化算法速度.方法采用优化设计与CAD相结合的方法.结果与结论解决了传统刀具设计的缺点,改进后的算法速度大幅度提高.
关键词:
数学模型;优化;算法
在传统的刀具设计中,通过查表和经验公式来确定各种结构参数和几何参数,然后,反复计算来得到相对较优的刀具参数.这种方法使设计过程复杂费时,且得不到最优化的参数,设计出的刀具成本高,加工效率低.因而刀具的计算机辅助设计应采用优化设计与CAD相结合的方法,欲进行优化设计,必需首先建立刀具优化设计的数学模型,由于复杂刀具的种类繁多,结构变[1],此篇仅以轮切式拉因而需分门别类地建立模型化多样,优化目标不同,.
刀为例1拉刀优化设计的数学模型
在拉刀参数设计过程中需要选择的主要参数有拉削余量A,齿升量a,f齿距t,容屑槽形状和深度h,容屑系数k,同时工作齿数等,这些参数可分为两类,一类是独立参数,如拉削余量和容屑槽形状等,这些参数基本不受其他如齿升量、齿距、容屑槽深度、容,另一类参数是非独立参数.参数的影响
屑系数等,这些参数既相互限制又相互依赖,第一类参数的选择比较容易.可以用经验公式和数据库来解决.第二类参数比较复杂,只有通过优化的方法才能得到较好的结果.
粗切齿升量的选择是一个比较复杂的问题.增大a可使齿数减少,拉刀f长度变短,但同时又要求容屑槽深度增加.另外齿升量的增加又会引起拉削力的增大,受到拉床和拉刀拉应力的限制.
齿距是决定拉刀长度的一个重要因素,t越大,拉刀越长,同时工作齿数越少.这样会在拉削过程中引起振动,生产效率低,降低刀具的使用寿命。
t过小,又会使容屑空间变小,从而限制了齿升量的增大.
其他参数如同时工作齿数z,容屑槽深度h,容屑系数k都是a和t的fi函数,只有当a和t选择后才能确定.从上述参数分析可知,a和t是拉刀ff设计的关键,在a和t之间应有一最佳组合值,使得a在拉床的额定应力ff和拉刀的许用应力范围内达到最高,即使拉刀的长度最小.
1.1目标函数的建立
确定以af和t为优化的自变量,A为切削余量.拉刀长度是与拉削生产率、成本及其工艺性能有关的参数,拉刀越短对使用和制造越有利,因而取粗[4]L作为优化目标切齿部分长度F=minL(a,t)=tA/(2a).
(1)1.2
ff约束条件的建立
h-1.13(kaL)≥0.
(2)
2/1容屑槽空间的限制1)wf式中h是与t有关的参数。
k为容屑系数,是与t和af有关的参数。
Lw为拉削长度.
F-pπD≥0.(3)
zc/wzi拉床额定拉力的限制2)e式中Fe为拉床额定拉力。
Dw为拉削后孔直径。
p为单位切削力。
zi为同时工作齿数,zi=INT(Lw/t)+1。
zc为组齿数.
[σ]–2pD0.(4)
≥拉刀许用拉应力的限制3)wzi/zcdmin
式中[σ]为拉刀许用拉应力。
dmin为拉刀最小直径.
11-z≥0.(5)最大同时工作齿数的限制4)iz-3≥0.(6)最小同时工作齿数的限制5)i25-t≥0.(7)最大齿距的限制6)t-4≥0.(8)
最小齿距的限制7)8)弧形槽能保证稳定的分屑要求的最大齿升量
h-a≥0。
f(D,n,z)-a≥0.(9)
fzfct-Int(2t)/2=0.(10)
的整数倍齿距应为0.59)1.3优化模型
:
a,t。
自变量f:
F=minL(a,t)=tA/2a。
目标函数ff:
g
(1)=h-1.13(kafLw)≥0。
g
(2)=F-pπD1/2约束方程wzi/zce≥0。
g(3)=[σ]-2pD≥0。
g(4)=11-z≥0。
wzi/zcdminig(5)=z-3≥0。
g(6)=25-t≥0。
g(7)=t-4≥0。
g(8)=h-a≥fi0。
g(9)=f(D,n,z)-a≥0。
g(10)=t-Int(2t)/2=0.fzc2优化算法
2.1标准算法
复合形法是一种采用直接搜索方式求解非线性规划问题的数值计算方法,这个方法可以在N维非线性约束的空间中自动选择并改进设计点,该方法的:
]3.[2一般步骤为{x}(i=1,2,…,n,n+1,…,m)1)m>n+1在可行域内生成个点构i
①初始顶点的形成,可以人工选定这里需要注意两个问题:
成初始复合形,②需要检验初始顶点是否满足约束条件,即检验其可行也可随机产生。
性。
2),f(x)≤将其由小到大的顺序重新编号计算各顶点的目标函数值,1f(x)≤…≤f(x)。
m23)xm-1,即点中心其余个点坏定除去最点的后复合形中确mi?
11?
。
xc=xi
1m?
1?
m4)xxx=x+α(x-x),α为映的映射点对中心点确定最坏点mmccca1.3。
0.9~一般取射系数5)x按则,束条件违果背了某个约检验映射点的可行性:
如a(x+x)/2x,xa把映射点向中心移动一半距离,反复直至映射aacx。
是可行点点a6)f(x),x,完成一次计算新的可行点的函数值用它代替最坏点ma迭代,回到第二步。
7)f(x)-f(x)<ε直到满足,,则终止.以上称为复合重复以上过程m1形法的“标准算法”,由于该算法的概念简单、容易实现,且能有效灵活地处理不等式约束问题,所以在结构化设计中得到广泛的应用.
2.2存在的问题
把上述标准算法应用于工程实际时,就会发现它还存在以下几个问题:
1)过多的可行性检验限制了其在优化设计中的有效应用.初始顶点生成和映射点的确定,都要进行可行性检验,在结构优化设计中,可行性检验其实质上就是结构分析过程,其计算量通常要占总工作量的80%以上,因此结.
必然会导致因计算时间过长而降低算法的效率,构分析次数过多
2)迭代过程中向极值点逼近的速度问题.开始若干次迭代(一次迭代是对于选取一个既满足约束条件又使目标函数值有所改善的新点所需的计算),目标函数值下降得很快,各顶点迅速接近极值点,一般来说,最初的(5~10)次迭代函数值下降得最快.随着迭代次数的增加,函数值的变化却越来越缓慢,也就是说,这时要使目标函数值有微小的改善,都要付出宝贵的计算时间.
3)局部最优点问题.用上述算法得到的最优点有可能是局部最优点,虽然可通过多取几个初始点,经计算后得到几个最优点,然后比较得到全局最优点,但这样必然会导致计算工作量的成倍增加.
2.3分层复合形法
针对标准算法中存在的问题,采用“分层复合形法”,它是对标准复合形法的改进,其基本思想是:
充分利用复合形法开始时目标函数值急剧下降的特点,以迭代次数为控制参数,进行两层优化计算,为避免产生局部最优点,在第一层迭代中,选取多组复合形分别地进行计算,经过若干次有效地迭代,各顶点迅速地逼近最优点,分布在最优点附近.分层复合形法的基本步骤如下:
(n)g,(i=1,…,…,{x},m)n{x1)n}
(1)个初组初始顶点构成选择iiggn=Int[n/2]+1,n为设计变量数.这里只要设计变量所取的值始复合形,g不太小且相互间离得远些,就可不对初始顶点作可行性检验.
2)2)~6)步的计算是第一层迭代,对各初始复合形标准算法第取映射率为αn。
迭代次数为,1t13)n次以后,第一层迭代结束,取两个最好设计点组成当各复合形都迭代t1α(α<α),迭代次数为新的复合形进入第二层迭代,取映射率为221n。
.
第二层迭代得到的最优点可被认为全局最优点t2分层复合形法有以下几个优点:
①迭代次数大大减少。
②以迭代次数为停止准则,可根据需要人工控制计算工作量。
③第二层迭代能有效地产生全局最优点.
3结论
依据本文所述方法,已开发出具有高效率优化CAD系统,证明对传统算法的.
改进是有效的.
参考文献:
[1]唐锡荣.CAD/CAM技术[M].北京:
北京航空航天大学出版社,1994.
18-36.
[2]蔡锁章.计算方法[M].北京:
中国科学技术出版社,1993.54-60.
[3]徐灏.机械设计手册.第二卷[M].北京:
机械工业出版社,1991.40-41.
[4]吴伏家,刘兆华.圆孔拉刀CAD系统研制[J].华北工学院院报,1996,
(增刊):
74-78.
机械优化设计案例3
直齿圆柱齿轮传动的优化设计
摘要:
一、问题描述:
现有一单级渐开线直齿圆柱齿轮减速器,其输入功率N=280kW,输入转速n=980r/min,传动比i=5。
小齿轮为实体结构,大齿轮为腹板式结构1(带有四个减轻孔),两齿轮各部分尺寸的符号如图一所示:
原用常规设计方法的设计结果为:
齿宽B=B=13cm,小齿轮齿数2z=21,模数m=0.8cm,l=42cm,d=12cm,d=16cm。
现要求在保证承载能s21s11力的条件下,通过优选上述有关参数,使减速器的体积达到最小。
二、建立优化设计目标函数:
齿轮传动优化设计中,设计变量一般选为齿轮传动的基本几何参数或性能参数,例如齿数、模数、齿宽系数、传动比、螺旋角、变位系数和中心距分离系数等。
齿轮传动的优化目标,较常见的是体积或质量最小,传动功率最大,工作寿命最长,振动最小,启动功率最小等。
现在选体积最小为优化目标,而减速器的体积主要是取决于内部零件(齿轮和轴)的尺寸大小,在齿轮和轴的结构尺寸确定之后,箱体的尺寸将随之确定,因此将齿轮和轴的总体积达到最小作为优化目标。
减速器内部有两个齿轮和两根轴,为了简化计算,将轴视为光轴,则有
V?
V?
V?
V?
Vg2s1s2g1?
?
22(l?
l?
l)?
)dd?
(l2ss1321144?
?
?
?
2'222222')dC?
?
D?
)(BC)4(?
B)?
?
(?
)dd(?
Bdd(D012s222211s14444
VV3——两轴体积,式中:
cm,;1s2sVV3cm,——两齿轮体积,2g1g
dd——两轴的直径,cm;,1s2slll——轴的长度,cm;,,321dd?
mzd?
mzd;,,——两齿轮的分度圆直径,cm,221121m——两齿轮的模数,cm;
BB?
B?
BB,cm。
——两齿轮的宽度,近似取,2121根据结构设计经验公式,齿轮各部分尺寸关系为:
?
?
5m
'?
1.6dD21sC?
0.2B
'?
2dD?
?
22'')?
0.25(DDd?
102l?
32cml?
28cm并取:
32优化设计中的设计变量取为:
TT]dd,lz,m,x,x,,x],?
[B,,X?
[x,xx24s53s116211将目标函数整理后得到:
2222x85xxx?
0.78539815(4.75xxxx?
?
85(fX)33122131222222)x?
?
28xxxxxx?
1.6xx?
x32?
xxxx?
0.92x?
x0.8?
x63645125163615641三、确定约束条件
zz?
z?
17,于是1()为了避免发生根切,不小于最小齿数,即min11得约束条件
g(X)?
x?
17?
021
(2)为了保证齿轮的承载能力,同时避免载荷沿齿宽分布严重不均,B?
135?
?
16?
16?
x?
X(g)x0要求,由此得:
321m
?
1?
x0?
35?
xg(X)313(3)传递动力的齿轮,模数一般应该大于2mm,并且去标准系列值,所以得:
g(X)?
x?
0.2?
034(4)根据工艺装备条件,要求大齿轮直径不得超过1500cm,于是小mz?
33cm,故得:
330cm,即齿轮直径相应的不能超过1g(X)?
33?
xx?
035210?
d?
15,为围按照经验取5()主、从动轴直径范1s13?
d?
20,所以有2sg(X)?
x?
10?
056g(X)?
x?
13?
067g(X)?
15?
x?
058g(X)?
20?
x?
069l?
B?
2?
?
0.5d?
为箱体)轴的支撑跨距按照结构关系(6,其中2s1?
?
2cm,则有:
内壁到轴承中心线的距离,现取
g(X)?
x?
x?
0.5x?
4?
061410(7)按齿轮的接触疲劳强度条件,有:
3KM1)i?
(1070?
?
1]?
[?
HHaBiKK?
1.3;——载荷系数,取式中:
Mi——小齿轮传递的扭矩,由功率和转速计算可得
Mi?
955000?
280/
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 机械 优化 设计方案 三个 案例