遗传算法 gui matlab,用matlab遗传算法工具gui优化轴承

用matlab遗传算法工具gui优化轴承

设计变量滚珠数目 Z=23 滚珠直径 Db 为变量 1，即 x(1) 节圆直径 dm 为变量 2，即 x(2)fi， fo 分别为轴承内外滚道曲率半径系数 fi=fo=0.515额定静载荷目标函数：32*3.814bioiZaCf/bmd， 分别为钢球与内滚道接触椭圆的无量纲化长半轴和短半轴*iai额定动载荷目标函数: 0.72/31.81.3(cos) (25.4bbCfZDm)./.44c )/bmd 0.3.1.391.720.41/ 0.41/2()39.1[.04()()]()io ic if ff约束条件 1,min ,min()()/0ball bbalgXCdZDC2i0DK3max()()2b40.5gXed5())(0m6.)bD7()03.51igXf8.o运用 matlab 的遗传算法工具 GUI 对适应度函数进行优化。打开 matlab，在 command window 中输入 gatool，打开遗传算法工具。初始界面如下图：在 fitness function 处输入所要求的适应度函数，由于函数过于复杂，先编写它的 M 文件:function y=dynamic(x) %函数名为dynamic ，输入变量为x(1)=D b滚珠直径,x (2 )=d m节圆直径r=x(1)/x(2); % =/(1)/bmDdxa=39.9*(1+(1.04*((1-r)/(1+r))^1.72)^(10/3))^(-0.3);%系数f c的前半部分b=r^0.3*(1-r)^1.39/(1+r)^(1/3);%系数f c的后半部分f=a*b*4.2623;%系数f c，由于=f o=0.515，即(2f i/(2f i-1))^0.41=4.2623y=-1.3*f*23^(2/3)*x(1)^1.8;%式中z=23 ，求出 y的值，即额定动载荷。由于遗传算法工具默认求函数最小值，所及在式子前加个负号，将求最大值的问题转化为求最小值得问题。将其保存为dynamic.m文件。在输入处输入@dynamic。在Number of variables处输入变量数目为2整理一下数学模型中的约束条件，有两个线性不等式，两个边界。以减速器型号为CSF_14_80,轴承内径为 25.2mm，轴承外径为34mm作优化设计。轴承约束条件中的系数取为如下：C ball,min=0.1Db，K Dmin=0.5，K Dmax=0.65，e=0.05， =0.3。约束条件 可转化为1 ,min())/0bbalgxdZDC由于遗传算法工具默认求最小值，所以约束条件中的不等/0.bdmZ式的格式要转化为左边式子小于等于右边不等式。最后可化为：即1.()(2)3x1.()0.36(2)0 xx同理，约束条件 可化为6/mbbgDd即.6()x.()4x所以线性不等式的矩阵为A:[1.1,-0.1366;1.6,1] b:[0;34]将其对应输入到Linear inequalities中。滚珠直径D b的取值范围 可化为：minmax()/2()/2DbDKdKd0.5(342.)/0.6534.b18x同理节圆直径d m的取值范围 即为：()/2(0.5+)mded)29.6()32.56x边界Bounds的矩阵即为：Lower:[2.2;29.6] Upper:[2.86;32.56]输入参数填完后，Options处的参数取默认值，默认的交叉概率为 0.8，变异概率为0.2，二进制编码长度为20。在 Plot functions处勾选best fitness 。点击“start ”按钮运行，可得如下图：通过图可看出极小值为-4159.246，对应的 x(1)=2.751,x(2)=29.6。从 best fitness 图中可看出最小值是在第 51 代时取到，平均适应度 mean fitness 逐渐减小，在第 10-20 代之间稳定在(-4100)----(-4200)之间，最后在 51 代取到最小值。如果将交叉概率设置为 0.5，可得如果将变异概率设置为 0.02，可得用相同方法对减速器型号 CSF_32_80，CSF_40_80，CSF_25_80 进行优化，可得：减速器型号 CSF_14_80 CSF_32_80 CSF_40_80 CSF_25_80轴承内径(mm)25.2 59 71.2 45.1轴承外径(mm)34 79.8 98.2 61.3滚珠直径(mm)2.751 6.501 8.438 5.063节圆直径(mm)29.6 69.4 84.7 53.2额定动载荷(N)4159.246 19593.680 31750.773 12538.740额定静载荷(N)1869.632 10436.588 17508.011 6324.072