Parameter optimization of a motorized spindle lubrication system using biogeography-based optimization

Abstract

Using response surface methodology, two-target parameter mathematical models for a 100MD60Y4 motorized spindle are established using biogeography-based optimization. Air supply pressure, fuel supply rate, and oil supply interval are the independent variables, while the spindle temperature increase and oil mist PM_2.5 cumulative volume distribution percentage are the response variables. Biogeography-based optimization is used to find a set of optimal lubrication system parameters with the constraint of a small temperature increase in the motorized spindle at 30,000 r/min while considering the influence of oil mist particle size on environmental quality. After lubrication system parameter optimization, comparing the particle size distribution density and temperature increase between the biogeography-based optimization model’s calculated values and the experimental test results, the errors were 1.8% and 1.17%, respectively.

Keywords

Oil–air lubrication system motorized spindle biogeography-based optimization ambient air quality Box–Behnken

Introduction

Computer numerical control (CNC) machine tools occupy a pivotal position in national economic development. Automation, intelligence, and green development are hallmarks of CNC machine tools. In order to meet the demands of low-carbon manufacturing and sustainable development, the design, manufacture, and use of CNC machine tools are necessary.^1,2 An important development is high-speed spindle technology, which is used in advanced cutting operations. However, these spindles experience elevated temperatures that affect equipment life. Accordingly, most motorized spindles employ oil–air lubrication which simultaneously provides bearing lubrication and rotor cooling. Although oil–air lubrication uses less oil than oil mist and atomize the oil, small oil droplets remain in the air, which we consider “oil mist” for purposes of environmental air quality and safety. Research shows that the health risk created by oil mist is closely related to the size of the oil particles. Oil particles less than 2.5 μm diameter (PM_2.5) enter the the bronchi and alveoli, allowing harmful substances in the oil to infiltrate the bloodstream, causing harm to the human body.^3–5 For this reason, the European machine tool manufacturing industry is relatively focused on the environmental performance of machine tools. Technologies such as dry cutting and hard cutting—processes based on using less cutting fluid and lubricant—are targeted at reducing machine tool operation environmental emissions. Weinert⁶ provided a detailed and comprehensive discussion on dry-cutting technology and cutting technology using less cutting fluid. The Michigan University of Technology’s Environmental Awareness Project has studied the waste generated by various mechanical processes of machine tools, developing various methods for predicting waste and proposals to reduce emissions via process optimization.⁷ Recently, scholars have studied oil–air lubrication systems and motorized spindles.^8–10 Lin¹¹ developed a genetic algorithm optimization approach to search for optimal bearing locations on a motorized spindle shaft. Jeng and colleagues compared different conditions of oil–air lubrication with mixed ceramic ball bearings and steel ball bearings while increasing temperature. The study accounted for the load, spindle speed, fuel supply rate, gas flow rate of oil, and operating parameters, studying how they impacted steel bearing temperature increase.^12,13 Huang and Zhang¹⁴ analyzed the frictional heating of a high-speed motorized spindle and the heat radiation of an oil–air lubrication system. Jiang and Mao¹⁵ investigated the effects of the operational parameters of an oil–air lubrication system on the temperature rise of both the hybrid ceramic ball bearing and steel ball bearing.

The above-mentioned studies separately considered environmental performance and processing efficiency; however, there has been no combined discussion of these topics. A challenge in the machinery manufacturing industry is how to achieve high-quality processing efficiency and effective ambient air quality control. The 100MD60Y4 motorized spindle, manufactured by Luoyang Bearing Research Institute Co., Ltd., is the subject of this study. Based on a single-factor experiment using response surface methodology, two-target parameter mathematical models are established. Air supply pressure, fuel supply rate, and oil supply interval are the independent variables, while the motorized spindle temperature increase and oil mist PM_2.5 cumulative volume distribution are the response variables, which are optimized using biogeography-based optimization (BBO). With the constraint of ensuring a small motorized spindle temperature increase, the influence of oil mist particle size on environmental quality is evaluated to find a set of optimal lubrication system parameters. The results are compared with those from the “fmincon” function, which can be used to solve the minimum value of a nonlinear multi-function in MATLAB.

Experimental method

In this article, a TC-2008 multi-channel thermometer was used to measure the temperature increase in the motorized spindle bearing outer races, and a laser particle size analyzer was used to measure the oil mist released from the spindle to the environment.

Single-factor experiments were done, with each experiment lasting 40 min. The experimental setup is shown in Figure 1.

Figure 1.

Experimental setup.

The spindle bearing temperature rise (ΔT) as a function of air supply pressure, fuel supply rate, and oil supply interval is shown in Figure 2. In each experiment, the spindle bearing temperature rise reached a steady state before the end of the 40 min. The cumulative distribution (Q3) and distribution density (q3*) of oil mist particles with air supply pressure, fuel supply rate, and oil supply interval are shown in Figures 3 and 4, respectively.

Figure 2.

Bearing temperature increase: effect of (a) air pressure when the oil supply interval is 0.2 s and the fuel supply rate is 24 mL/h, (b) oil supply interval when the air supply pressure is 0.29 MPa and the fuel supply rate is 24 mL/h, and (c) fuel supply rate when the oil supply interval is 0.2 s and the air supply pressure is 0.29 MPa.

Figure 3.

Cumulative distribution of oil mist particles: effect of(a) air pressure when the oil supply interval is 0.2 s and the fuel supply rate is 24 mL/h, (b) oil supply interval when the air supply pressure is 0.29 MPa and the fuel supply rate is 24 mL/h, and(c) fuel supply rate when the oil supply interval is 0.2 s and the air supply pressure is 0.29 MPa.

Figure 4.

Distribution density of oil mist particles: effect of (a) air pressure when the oil supply interval is 0.2 s and the fuel supply rate is 24 mL/h, (b) oil supply interval when the air supply pressure is 0.29 MPa and the fuel supply rate is 24 mL/h, and (c) fuel supply rate when the oil supply interval is 0.2 s and the air supply pressure is 0.29 MPa.

The impact of these parameters is shown to have a significant effect on air quality. The oil mist particle size PM_2.5 cumulative volume with time under different conditions is shown in Figure 5. The oil mist particle size PM_2.5 cumulative volume does not change with time under certain conditions.

Figure 5.

Oil mist particle size PM_2.5 cumulative volume with time.

Model design

Referencing the previous single-factor experiment studying lubrication system parameters and using the Box–Behnken experimental design criteria, three levels (high, medium, and low) for each of the three factors were chosen. These values were 0.29, 0.35, and 0.41 MPa for air supply pressure; 6.5, 15.25, and 24 mL/h for fuel supply rate; 0.6, 90.6, and 180 s for oil supply interval; and the spindle speed was fixed at 30,000 r/min. According to the Box–Behnken center combination using three response-point analyses with three factors and three levels, a total of 17 experimental points from groups 1 to 12 were used for the analysis; the last five sets of data for the central experiment were used to estimate the error (Table 1). The advantage of this design is that the number of experiments is relatively small; since all factors are not at a high level at the same time, the experimental design ensures that all design points are in a safe range.

Table 1.

Design of a three-factor, three-level Box–Behnken experiment.

Experiment number	Independent variable			Response variable		Normalization
Experiment number	Air supply pressureG (MPa)	Fuel supplyrate B (mL/h)	Oil supply interval t (s)	Bearing outer ring temperatureincrease $Δ T$ (°C)	PM_2.5 cumulativevolumedistribution $Q 3$	$Δ T^{*}$ (°C)	$Q 3^{*}$
1	0.29	6.5	90.6	20.4	0.287	0.8833	0.2845
2	0.41	6.5	90.6	15.5	0.077	0.0667	0.074
3	0.29	24	90.6	19	0.578	0.65	0.5763
4	0.41	24	90.6	16	1	0.15	1
5	0.29	15.25	0.6	21.1	0.005	1	0.0017
6	0.41	15.25	0.6	16.2	0.766	0.183	0.7648
7	0.29	15.25	180	19	0.652	0.65	0.6516
8	0.41	15.25	180	15.1	0.003	0	0
9	0.35	6.5	0.6	16.2	0.305	0.183	0.3034
10	0.35	24	0.6	15.7	0.054	0.1	0.0517
11	0.35	6.5	180	16.1	0.02	0.167	0.0171
12	0.35	24	180	16.3	0.147	0.2	0.1442
13	0.35	15.25	90.6	16.6	0.097	0.25	0.094
14	0.35	15.25	90.6	16.6	0.141	0.25	0.138
15	0.35	15.25	90.6	16.7	0.158	0.267	0.1559
16	0.35	15.25	90.6	16.4	0.162	0.217	0.1593
17	0.35	15.25	90.6	16.8	0.141	0.283	0.138

The results of these experiments are listed under “response variable” in Table 1. In the experimental results, the variation range of the motorized spindle temperature was 15.1°C–21.1°C, and the oil mist particle size PM_2.5 cumulative volume distribution was 0.003–1. The data range was relatively wide, and the order of magnitude was not uniform, leading to inaccurate predictions. Therefore, the min–max normalization method was used to normalize the data to reduce the data range and improve prediction and model accuracy. The original data were linearly transformed with equation (1) to map the resulting value between 0 and 1, where max and min are the maximum and minimum values of the response variable data, respectively. The results are listed in Table 1 under the “normalization” heading; $x^{*}$ is denoted by $Δ T^{*} Q 3^{*}$

x^{*} = \frac{x - min}{max - min}

(1)

The response surface method is a regression method for function estimation that provides a functional relationship between the factors in a multi-factor experiment and the experimental response values. Surface analysis of the function and quantitative analysis of each factor was conducted and their influence on the response values was determined. For the multifactor experiments, ΔT and Q3 were measured at t = 40 min (the end of the experiment).

When the regression model of $Q 3^{*}$ was analyzed, the quadratic term was significant. Therefore, the quadratic polynomial was chosen to fit oil mist particle sizes less than 2.5 µm in the cumulative volume distribution, $Q 3^{*}$ , for the spindle lubrication system outlet. The resulting response surface formula is shown in equation (2)

\begin{matrix} Q 3^{*} (G, B, t) = 9.25 - 53.61 \times G - 0.13 \times B \\ + 1.34 \times t + 0.30 \times GB - 3.94 \times Gt + 0.007 \times Bt \\ + 79.46 \times G^{2} + 7.92 \times 10^{- 4} \times B^{2} - 0.03 \times t^{2} \end{matrix},

(2)

where G is the air supply pressure, B is the fuel supply rate, t is the oil supply interval, and $Q 3^{*}$ is the oil mist PM_2.5 cumulative volume distribution.

To evaluate the reliability of the experimental results and the fitting accuracy of the model, the results were tested for significance. First, a variance experiment was conducted, determining the significance of each factor. Second, the significance of the overall experiment was determined. Table 2 shows the partial regression analysis of the oil mist PM_2.5 cumulative volume distribution percentage for the spindle lubrication system outlet.

Table 2.

Variance of the PM_2.5 cumulative volume distribution regression equation.

Project	Sum of squaresof deviations	Freedom W	Mean squaredeviation	F experimentvalue	P experimentvalue
Model	1.19	9	0.13	4.02	0.040
B	0.013	1	0.013	0.40	0.5467
G	0.15	1	0.15	4.55	0.0705
t	0.012	1	0.012	0.36	0.5661
B × G	0.10	1	0.1	3.06	0.1237
B × t	0.50	1	0.5	15.23	0.0059
G × t	0.036	1	0.036	1.09	0.3308
B ²	0.34	1	0.34	10.48	0.0143
G ²	0.015	1	0.015	0.47	0.5147
t ²	0.020	1	0.020	0.60	0.4631
Residual	2.705E−003	4	6.764E−004	–	–
Correct square sums	1.42	16	–	–	–

Multiple correlation coefficient R² = 0.9298; correction coefficient $R_{adj}^{2} = 0.93805$ ; model resolution = 67.77.

From the analysis of variance in Table 2, the experimental value P = 0.04 < 0.05 indicates that the model in equation (2) was significant. The correction coefficient of the model was 0.93805, which indicates that the model can explain a change of 93.80% in the response value; that is, only 6.20% of the total variation cannot be explained by it. The multiple correlation coefficient was 0.9298, indicating a strong correlation between the predicted value and the measured value; the model fits well. The model resolution was 67.66 > 4, indicating sufficient resolution. Thus, equation (2) is suitable and can be used to analyze and predict the oil mist PM_2.5 cumulative volume distribution percentage.

When the regression model of $Δ T^{*}$ was analyzed, the quadratic term was significant, so the quadratic poly-nomial was chosen to fit the bearing outer ring temperature increase in the spindle $Δ T^{*}$ . The response surface formula is shown in equation (3)

\begin{matrix} Δ T^{*} (G, B, t) = 11.44 - 55.45 \times G - 0.04 \\ \times B - 0.19 \times t + 0.15 \times G \times B + 0.47 \\ \times G \times t + 0.002 \times B \times t \\ + 66.65 \times G^{2} - 7.29 \times 10^{- 4} \times B^{2} - 0.02 \times t^{2} \end{matrix} .

(3)

From the analysis of variance in Table 3, the experimental value $P = 0.0004 < 0.01$ indicates that the model in equation (3) was very significant. The correction coefficient of the model was 0.9078, indicating that it can explain a change of 90.78% in the response value; 9.22% of the total variation cannot be explained by this model. The multiple correlation coefficient was 0.9597, which indicates strong correlation between the predicted value and the measured value; the model fits well. The model resolution was 27.17 > 4, indicating that the model had sufficient resolution. Thus, equation (3) is suitable and can be used to analyze and predict the spindle bearing temperature increase.

Table 3.

Variance of the motorized spindle temperature increase regression equation.

Project	Sum of squares of deviations	FreedomW	Mean squaredeviation	F experimentvalue	P experimentvalue
Model	1.29	9	0.14	18.51	0.0004
B	0.97	1	0.97	125.38	<0.0001
G	0.005	1	0.005	0.65	0.4475
t	0.025	1	0.025	3.26	0.1139
B × G	0.025	1	0.025	3.24	0.1147
B × t	0.007	1	0.007	0.90	0.3737
G×t	0.003	1	0.003	0.44	0.5304
B ²	0.24	1	0.24	31.38	0.0008
G ²	0.013	1	0.013	1.70	0.2337
t ²	0.005	1	0.005	0.67	0.4398
Residual	2.409E−003	4	6.023E−004	–	–
Correct square sums	1.34	16	–	–	–

Multiple correlation coefficient R² = 0.9298; correction coefficient $R_{adj}^{2} = 0.93805$ ; model resolution = 67.77.

The optimal design of the response surface model requires further statistical feature diagnosis to determine the accuracy. One method is the normal probability graph of the student’s residual, which is used to show whether the residuals conform to the normal distribution. If the points are concentrated on the plane and distributed along the line, then they conform to a normal distribution. The spindle temperature increase and oil mist PM_2.5 cumulative volume distribution represent the two target parameters of the response surface model; the relevant student’s residual normal probability diagram is shown in Figure 6, which indicates that the response surface model established was statistically accurate.

Figure 6.

Normalized plot of residuals for the target parameter response surface models: (a) M* and (b) T*.

Two-objective optimization

The objective functions $Δ T^{*} (G, B, t)$ and $Q 3^{*} (G, B, t)$ were established using the response surface method to evaluate the influence of the operation parameters of the motorized spindle lubrication system on the temperature increase and the oil mist PM_2.5 cumulative volume distribution percentage. These two parameters were arbitrarily given weights of 0.6 and 0.4 for $Δ T^{*} (G, B, t)$ and $Q 3^{*} (G, B, t)$ , respectively. The minimum value of the nonlinear function was given as $min y = F (x)$ under the constraint condition, and the resulting mathematical model is described as follows

\begin{matrix} min y = F (x) = 0.6 \times Δ T^{*} (x) + 0.4 \times Q 3^{*} (x) \\ 6.5 ⩽ B ⩽ 24 \\ 0.29 ⩽ G ⩽ 0.41 \\ 0.01 ⩽ t ⩽ 3 \\ x = (G, B, t) \end{matrix}

(4)

Design procedure

First stage

There are many methods for solving constrained optimization.^16,17 Traditional global optimizations such as genetic algorithms and particle swarm algorithms are the best tools to solve such problems,^18–20 but their search and convergence accuracy have yet to be improved. BBO²¹ is one of the best intelligent optimization algorithms for global optimization. It is superior to the genetic algorithm and other optimization algorithms in terms of convergence accuracy and velocity and is especially suitable for optimization in practical engineering applications. In function optimization, production scheduling, heart disease diagnosis, and other fields, BBO has been applied successfully.²²

The BBO algorithm seeks an optimal global solution by simulating the mechanism of species migration between habitats. In the BBO algorithm, the biological population lives in different habitats H, and each habitat determines the habitat adaptation characteristics according to the number of species. A habitat suitable for the species has a higher habitat suitability index (HSI). HSI represents the evaluation of fitness according to the objective function. When optimizing, good solutions have a high-HSI habitat, whereas poor solutions have a low-HSI habitat. The BBO algorithm needs to construct multiple groups randomly, and the corresponding suitability index vector (SIV) is the initial solution of the optimization problem. Through evolution of the population, the HSI improves and the optimal solution is obtained.^23–25 The evolution process involves population migration and population mutation. In this article, the relationships between the BBO algorithm and lubrication system parameters are shown in Table 4.

Table 4.

Relationships between BBO and the optimization parameters of the lubrication system.

BBO parameters	Lubrication system
Population migration mutation process	Optimization of lubrication system parameters
Habitat H	F(x) possible solutions
SIV	Lubrication system parameters
Habitat index variable	G, B, t
HSI	Optimize the objective function $F (x)$
Highest HSI	Optimal solution of $F (x)$

BBO: biogeography-based optimization; SIV: suitability index vector; HSI: habitat suitability index.

The BBO algorithm uses the migration mechanism to complete the exchange of information with other groups; a wide scope search is then performed in the solution space of $F (x)$ . The immigration rate $λ_{k}$ and emigration rate $μ_{k}$ , which are functions of the habitat species k, determine the degree of group migration. A high-HSI habitat contains more species, and the number of species is close to saturation; therefore, the corresponding species emigration rate $μ_{k}$ is high and the immigration rate $λ_{k}$ is low. In contrast, a low-HSI habitat corresponds to fewer species, so the immigration rate is high and the emigration rate is low. This article chooses cosine function model with optimal performance.²⁶

In the cosine models, I is the maximum immigration rate, E is the maximum emigration rate, n is the number of species at saturation, and k₀ is the number of species when the group moves to a new equilibrium. The expressions of the immigration rate $λ_{k}$ and the emigration rate $μ_{k}$ in cosine-type functions are as follows

λ_{k} = \frac{I}{2} [\cos (\frac{k π}{n}) + 1]

(5)

μ_{k} = \frac{E}{2} [- \cos (\frac{k π}{n}) + 1]

(6)

Each migration operation is determined by the immigration rate $λ_{k}$ and emigration rate μ_k. The immigration rate $λ_{k}$ of each habitat H is calculated from equation (5), and the habitat $H_{i}$ to be moved into is determined according to $λ_{k}$ . At this point, the SIV components of $H_{i}$ need to be modified. Using equation (6), the emigration rate $μ_{k}$ of each habitat H is calculated, and which habitat $H_{j}$ needs to be moved out of is determined according to $μ_{k}$ . Finally, the SIV of habitat $H_{j}$ replaces the SIV of habitat $H_{i}$ . The fitness of the whole group is evaluated by calculating HSI, that is, $F (x)$ .

Another evolutionary process is the mutation operation. The BBO algorithm simulates the change of the suitability of the population due to a sudden event. The SIV is mutated according to the probability $P_{k}$ , where the number of species is k. $P_{k}$ satisfies the following expression

P_{k} = {\begin{matrix} - (λ_{k} + μ_{k}) P_{k} + μ_{k + 1} P_{k + 1}, k = 0 \\ - (λ_{k} + μ_{k}) P_{k} + λ_{k - 1} P_{k - 1} + μ_{k + 1} P_{k + 1} \\ 1 ⩽ k ⩽ n - 1 \\ - (λ_{k} + μ_{k}) P_{k} + λ_{k - 1} P_{k - 1}, k = n \end{matrix}

(7)

The mutation rate of $m_{i}$ is

m_{i} = m_{max} \frac{1 - P_{k}}{P_{max}}

(8)

where $m_{max}$ is the maximum mutation rate defined arbitrarily and $P_{max}$ is the maximum population probability of the habitat.

In the mutation operation, according to m_i, group $H_{i}$ , which needs to be modified, is selected; then, the SIV is randomly generated to replace the corresponding SIV in $H_{i}$ to obtain a new HSI value.

Second stage

The BBO algorithm steps are summarized as follows:

Step 1. Set the number of habitats S and the maximum number of species n. Set the maximum value of the immigration rate function I, the maximum value of the emigration rate E, and the maximum mutation rate $m_{max}$ . Set the range of SIV components, randomly select multiple groups, and initialize the suitability vector SIV of each group.

Step 2. Calculate the HSI value of the population and arrange them in ascending order. Based on this, determine the number of species k contained in each group. Calculate the groups’ $λ_{k}$ and $μ_{k}$ according to equations (5) and (6).

Step 3. Perform the species migration operation according to the immigration rate $λ_{k}$ and the emigration rate $μ_{k}$ .

Step 4. Calculate the species probability $P_{k}$ of each group using equation (7) and determine the mutation rate m_i according to equation (8).

Step 5. According to the variation rate of each habitat, the species mutation operation is performed and the HSI is recalculated. The SIV corresponding to the HSI is selected.

Step 6. Determine whether the termination conditions are met. If the conditions are satisfied, end the cycle and output the optimal solution. Otherwise, return to step 2 to continue the search. Using BBO, the objective functions of $Q 3^{*}$ and $Δ T^{*}$ could be used to solve the double-objective problem shown in equation (4). The number of habitats S was 20, the number of generations N was 200, the maximum number of species was 3, the maximum value of immigration rate I was 1, the maximum value of emigration rate E was 1, and the maximum mutation rate $m_{max}$ was 0.05.

After optimization by BBO, the final optimal parameters were G = 0.38 MPa, B = 6.5 mL/h, and t = 2.98 min. The $F (x)$ fitness curve is shown in Figure 7. The curve tended to be smooth when the number of iterations exceeded 100 and tended toward −0.348, where $Δ T^{*}$ was 13.1°C, and $Q 3^{*}$ was 0.0628.

Figure 7.

$F (x)$ fitness curve.

Verification of optimization results

BBO and fmincon function method comparison

The conventional solution to this nonlinear constraint problem is to use the fmincon function provided in MATLAB software. According to the model $F (x)$ established in this article, the optimal objective function value $F (x) = - 0.1655$ was found by the fmincon function. Comparative results are shown in Table 5.

Table 5.

Optimization results using the fmincon function compared with BBO.

Item	fmincon optimum	BBO optimum
G (MPa)	0.3988	0.38
B (mL/h)	6.5	6.5
t (s)	3	2.98
$F (x)$	−0.1655	−0.348

BBO: biogeography-based optimization.

Clearly, the value obtained by the fmincon function was greater than the optimal value $F (x) = - 0.348$ obtained by the BBO algorithm. The fmincon function can solve the nonlinear constraint problem but requires the objective function, the constraint function, and the gradient to be continuous. The function is based on the initial value set by the artificial gradient, and dependence on the initial value is very strong.

Experimental verification

According to the optimal BBO parameters, that is, G = 0.38 MPa, B = 6.5 mL/h, and t = 2.98 min, the optimized values are $Δ T^{*} = 13 . 1^{°} C$ and $Q 3^{*} = 0.0628$ . Experiments were conducted to verify the accuracy of these values. A laser particle size analyzer was used to measure the oil mist under this condition as shown in Figure 8.

Figure 8.

Particle size cumulative distribution and distribution density.

The abscissa in Figure 8 denotes the particle size. The left and right ordinates represent the cumulative distribution Q3 and the density distribution q3*, respectively. As seen in Figure 8, the particle size was 1–40 μm, and the particle size distribution density was maximized at 12.5 μm. The initial temperature of the bearing outer rings of the motorized spindle was 17.8°C then increased to, and remained at, 31.1°C for the duration of the 40-min test. The temperature increase is shown in Figure 9.

Figure 9.

Bearing outer rings of the motorized spindle temperature increase.

The particle size distribution density and temperature increase errors were 1.8% and 1.17%, respectively.

Summary and conclusion

In this article, two-target parameter mathematical models were established based on single-factor experiments, using response surface methodology. The air supply pressure, fuel supply rate, and oil supply interval were the independent variables, while the motorized spindle temperature increase and oil mist PM_2.5 cumulative volume distribution were the response variables. The accuracy of the model was verified by variance and correlation tests.

The proposed model utilizes the BBO algorithm to optimize the parameters of the motorized spindle lubrication system. The results show a target parameter temperature increase of 13.1°C and oil mist PM_2.5 cumulative volume distribution of 0.0628 for lubrication system parameters of G = 0.38 MPa, B = 6.5 mL/h, and t = 2.98 min. The accuracy of the results was experimentally verified.

The proposed model can guide the design of operating parameters during operation of the motorized spindle, ensuring that the spindle system’s dynamic characteristics achieve green environmental goals.

Footnotes

Handling Editor: Xichun Luo

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by the National Science Foundation (51375317 and 51675353), Liaoning Provincial Department of Education Project (LJZ2017005 and LJZ2017035), Science and Technology fund of Liaoning Province (20180550002), and the (local) National Joint Engineering Laboratory for open funds (sjsc-2015-6).

ORCID iD

Dong An

References

Lin

Zhang

Zhong

Application of neural networks to predict the elevated temperature flow behavior of a low alloy steel. Comput Mater Sci 2008; 43: 752–758.

Yan

Liu

et al . Development status and technical system framework of intelligent machine tool. Chinese J Mech Eng 2013; 49: 1–10.

Ali

Umemura

Capillary phenomena and disintegration processes of a square cylindrical liquid column. Int J Heat Mass Tran 2010; 53: 5570–5580.

Park

Olortegui-Yume

Yoon

et al . A study on droplets and their distribution for minimum quantity lubrication (MQL). Int J Mach Tool Manuf 2010; 50: 824–833.

Wlaschitz

Hoflinger

A new measuring method to detect the emissions of metal working fluid mist. J Hazard Mater 2007; 144: 736–741.

Weinert

Inasaki

Sutherland

et al . Dry machining and minimum quantity lubrication. CIRP Ann Manuf Techn 2004; 53: 511–537.

Crittenden

JC.

National center for clean industrial and treatment technologies (Cen CITT). Clean Tech Environ Policy 2000; 2: 80–81.

Lin

Model-based design of motorized spindle systems to improve dynamic performance at high speeds. J Manuf Process 2007; 9: 94–108.

Shin

Analysis of bearing configuration effects on high speed spindles using an integrated dynamic thermo-mechanical spindle model. Int J Mach Tool Manuf 2004; 44: 347–364.

10.

Maeda

Cao

Altintas

Expert spindle design system. Int J Mach Tool Manuf 2005; 45: 537–548.

11.

Lin

CW.

Optimization of bearing locations for maximizing first mode natural frequency of motorized spindle-bearing systems using a genetic algorithm. Appl Math 2014; 5: 2137–2152.

12.

Jeng

Huan

PY.

Temperature rise of hybrid ceramic and steel ball bearing with oil-mist rolling management. Lubr Eng 2000; 56: 18–23.

13.

Jeng

Gao

CC.

Investigation of the ball-bearing temperature rise an oil-air lubrication system. Proc IMechE, Part J: J Engineering Tribology 2001; 215: 139–148.

14.

Huang

Zhang

BL.

FEA of thermal properties for high speed motorized spindle. Aeronaut Manuf Technol 2003; 10: 20–26.

15.

Jiang

Mao

HB.

Investigation of the high speed rolling bearing temperature rise with oil-air lubrication. J Tribol 2011; 133: 021101.

16.

Yang

ZP.

A QP-free algorithm without a penalty function or a filter for nonlinear general-constrained optimization. Appl Math Comput 2018; 316: 52–72.

17.

Alexander

Applied statistics and probability for engineers. Technometrics 2003; 37: 455–457.

18.

Fiandaca

Fraga

Brandani

A multi-objective genetic algorithm for the design of pressure swing adsorption. Eng Optim 2009; 41: 833–854.

19.

Gomes

Silva

Some comparisons for damage detection on structures using genetic algorithms and modal sensitivity method. Appl Math Model 2008; 32: 2216–2232.

20.

Chong

Zak

An introduction to optimization. IEEE Antenna Propag Mag 2014; 38: 60.

21.

Cai

Long

Jiao

JJ.

Hybridizing artificial bee colony with biogeography-based optimization for constrained mechanical design problems. J Cent South Univ 2015; 22: 2250–2259.

22.

Panchal

Singh

Kaur

et al . Biogeography based satellite image classification. Int J Comput Sci Inform Secur 2009; 6: 269–274.

23.

Simon

Ergezer

Population distributions in biogeography-based optimization algorithms with elitism. In: IEEE international conference on systems, San Antonio, TX, 11–14 October 2009, pp.991–996. New York: IEEE.

24.

Simon

Ergezer

Markov models for biogeography-based optimization and genetic algorithms with global uniform recombination. IEEE Trans Syst Man Cybernet Part B 2009; 41: 299–306.

25.

Simon

A probabilistic analysis of a simplified biogeography-based optimization algorithm. Evol Comput 2014; 19: 167–188.

26.

Liu

Wang

Zhang

An improved biogeography-based optimization algorithm for blocking flow shop scheduling problem. Chinese J Electron 2018; 27: 351–358.