基于MATLAB算法的深沟球轴承的优化设计

基于MATLAB算法的深沟球轴承的优化设计

摘要:基于MATLAB提出了深沟球轴承优化设计方法。以摩擦力矩最小作为优化设计的目标函数，利用MATLAB优化工具箱算法对深沟球轴承进行优化设计，得到最优设计结果。计算实例表明，与轴承的设计参数相比，优化后，轴承的摩擦力矩减小了25%，降低了摩损和发热量。从而说明了利用MATLAB优化设计方法是一种有效的优化设计方法。

关键词:MATLAB;摩擦力矩;深沟球轴承;优化设计

MATLAB是矩阵实验室(Matrix Laboratory) 的简称，和Mathematical、Maple 并称为三大数学软件。它在数学类科技应用软件中，在科学计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、图像处理、信号检测、金融建模设计与分析等领域。在机械领域，深沟球轴承运用在各个方面。轴承的磨损和发热是轴承失效的重要原因。因此，我们寻求一种方法来优化轴承的结构，降低轴承失效率，提高商用轴承的质量和使用寿命。

1、MATLAB优化工具箱描述

MATLAB优化工具箱可以求解线性、非线性最小化、最大和最小化、二次规划、半无限问题、线性、非线必方程(组)以及线性和非线性最小二乘问题。另外，MATLAB优化工具箱还提供线性和非线性最小化、方程求解、曲线拟合以及二次规划问题的求解方法。

1.1 线性规划

线性规划标准形式为:

求X=[x1 x2 L x n]T

使f(X)=c1x1+ c2x2+L+ c n x n最小

且满足a11x1+ a12x2+L+ a1n x n=b1

a21x1+ a22x2+L+ a2n x n=b2

说明：

1）当m=n时有唯一解a m1x1+ a m2x2+L+ a mn x n=b m

2) 当m>n时无解x i≥0 （i=1，2，l，n）

3) 当m

1.2 非线性规划

非线性规划问题包括无约束多元函数极小值和有约束函数极小值问题。1.2.1 无约束多元函数极小值

无约束多元函数极小值标准形式为:

minf (x)

x [ x1，x2，...，x n] 求解无约束多元函数极小值的fminsearch函数标准形式为:[x ，fval]=fminsearch (fun，x o，options，P1，P2，⋯⋯) 其中，P1，P2为函数fun的附加参数。

1.2.2 有约束多元函数极小值

有约束多元函数极小值标准形式为:

minf (x)

sub.to C(x)≤0

Ceq(x) 0

A·x≤b

Aeq ·x b eq

lb ≤x ≤ub

其中，Aeq—等式约束函数的系数矩阵;beq—等式约束函数的常数向量;lb—设计变量下限;ub—设计变量上界；C ( x ) —等式约束函数的返回向量; Ceq(x)—不等式约束函数的返回向量。

求解有约束多元函数极小值的fmincon函数标准形式为:

[x，fval]=fmincon(fun，x0，A，b，Aeq，beq，lb，ub，nonlcon，options，P1，P2，⋯⋯)

其中，nonlcon—非线性约束函数的调用函数。

2、深沟球轴承的优化设计

深沟球轴承是最具代表性的滚动轴承，用途广泛。适用于高转速甚至极高转速的运行，而且非常耐用，无需经常维护。该类轴承摩擦系数小、极限转速高、结构简单、制造成本低，易达到较高制造精度。尺寸范围与形式变化多样，应用在精密仪表、低噪音电机、汽车、摩托车及一般机械等行业，是机械工业中使用最为广泛的一类轴承。在电器行业领域，深沟球轴承也可以用在商用空调风机、电机等部件中。

2.1 优化设计变量和目标函数

深沟球轴承的摩擦力矩是其磨损和发热量的体现，把摩擦力矩作为目标函数，影响这个目标函数值的主要因素有滚珠直径、滚珠所在圆的直径、滚珠数量以及滚道位置角等几个参数，把这些参数作为设计变量，寻找摩擦力矩最小时的设计变量的最优解。

2.1.1 设计变量

X [d，D，Z，α]T [x1，x2，x3，x4]T

其中:d —滚珠直径; D —滚珠所在圆的直径; Z—滚珠数量; α—滚道位置角。

2.1.2 目标函数

目标函数是设计问题寻优的准则和基础。本文以摩擦力矩最小为目标函数，则有:

式中，M —摩擦力矩;F a —轴承的轴向载荷;Fr—轴承的径向载荷。于是，目标函数为：

2.2 约束条件的建立

(1)滚珠的接触强度条件：

式中，HB为布氏硬度，于是，约束函数为:

(2)设计结构尺寸条件:

6≤d ≤ 14.5

1.8(D) ≤Z ≤

2.38(D)

30 ≤ ≤ 60

40 ≤D ≤ 250

约束函数为:

g2(x) 6 x1≤0

g3(x) x1 14.5≤0

g4(x) 1.8(x2/x1) x3 ≤0

g5(x) x3 2.38((x2/x1 )≤0

g6(x) 30 x4≤0

g7(x) x4 60 ≤0

g8(x) 40 x2 ≤0

g (x) x 250≤0

3 优化实例和结果分析[4]

已知，轴承轴向载荷Fa =1000N，轴承径向载荷Fr=2500N。应用MATLAB优化工具箱函数fmincon来求解最优解。

3.1 编写MATLAB优化设计程序

3.1.1 主程序

x0=[13 220 40 43];%设计变量的初使值;

lb=[6 40 5 30];%设计变量的下限值;

ub=[14.5 250 200 60];%设计变量的上限值;

[x，fval]=fmincon(@xiul8_4，x0，[]，[]，[]，[]，lb，ub，@

fcon). %求解程序

3.1.2 目标函数子程序

function f=xiul8_4(x)

f=0.04*((1.5*1000+1.5*2500*tan(pi/180

*x(4)))/x(3)/sin(pi/180*x(4)))*x(1)/x(2);

3.2 优化结果

x(1)= 12.9697;%滚珠直径%

x(2)= 250.0000;%滚珠所在圆的直径%

x(3)= 45.8762;%滚珠数量%x(4)= 35.3921;%滚道位置角%

f= 0.3252;%摩擦力矩%