Abstract
We address the problem of optimizing the performance of lubricated devices by means of artificial texturing. We consider a slider (or equivalently a thrust bearing) and minimize the friction using optimization tools such as sensitivity analysis and genetic algorithms. We show that textures that perform significantly better than the smooth (untextured) one can be found, and that the optimized texture depends on the working conditions (load, velocity). The GENESIS code that we used, with no fine tuning of algorithmic variables, proved a valuable tool in the identification of improved shapes.
Nomenclature
h | = | =Local film thickness |
hp | = | =Depth of dimples |
hmin | = | =Minimum film thickness |
Ng | = | =Number of groups of dimples |
ac | = | =Square cell side |
ap | = | =Edgelength of the dimples |
as | = | =Half distance between dimples |
k | = | =Steepness parameter |
x | = | =coordinate |
y | = | =coordinate along film thickness |
z | = | =coordinate |
p | = | =Local pressure |
U | = | =Velocity |
W | = | =Load capacity |
Wa | = | =Applied load |
F | = | =Friction force |
1. Introduction
Static and dynamic characteristics of lubricated devices are very sensitive to geometrical changes in the surfaces. The question then arises as to what improvements can be brought to, for example, a slider (which is also a model of a thrust bearing) by texturing its surface by some of the newly available methods such as indentation, chemical etching, laser ablation, etc. A few recent works have suggested that significant improvement can be brought to lubricated devices by artificial texturing Citation1–3.
This report addresses a preliminary study in this direction that makes use of optimization tools such as sensitivity analysis and genetic algorithms. We consider the simple case of a slider loaded with a static, vertical force. Only a texture on the static surface is considered, so that the texture distribution does not move with respect to the load and the problem has a time-independent steady state.
This simple case is important to understand the basic mechanisms that could bring improvement to the slider's performance. We study square texture cells with square microdimples, and optimize the dimples' size and depth.
The methodology is the solution of the Reynolds' equation by finite element method coupled to a Newton-like method to find the equilibrium position. At this position the relevant quantities, such as friction force and minimal film thickness, can be evaluated so as to drive the optimization algorithm.
The conclusions of this preliminary assessment are, basically, that optimization tools (either deterministic or stochastic ones) can indeed find improved designs of lubricated devices, even when the design variables are in fact at the micro-scale of the surface's texture. We considered it important to assess evolutionary algorithms because in some areas of tribology the current trend is towards very sophisticated shapes that are not parameterized variations of previous ones. This article reports on a first step towards a systematic treatment of “optimum tribological surfaces”, or “engineered tribological surfaces”, by means of evolutionary algorithms.
2. Theory and definitions
2.1. Governing equation
The basic relation governing the slider lubrication problem is obtained by introducing the Reynolds assumptions into the Navier--Stokes equation. This relation can be written in a non-dimensional form for the stationary case, in the (x, z) coordinate system (see ) as: (1) where Ω is the integration domain and its boundary is defined as Ω = Γ D\cup Γ N}, being ΓD the Dirichlet part of the boundary and ΓN the Neumann one.
We consider two geometries. The first one is periodic in the width direction, and thus models a slider of infinite width, while the second one corresponds to a finite-width slider with atmospheric pressure along the sides. The corresponding boundary definitions are:
• | Infinite width | ||||
• | Ω = [qb;0,1]qb; ; × (- (ac/2(ac/2)) | ||||
• | Γ _N = \{(x,z) \,|\, x ∈ [qb;0,1]qb; ; , z = \pm ({a_{\rm c}}/{2})\} | ||||
• | Finite width | ||||
• | Ω = [qb;0,1]qb; ; × [qb;Vh ;0, B]qb;Vh ; | ||||
• | Γ N = {(x,z) ,|, x ∈ [qb;0,1]qb; ;, z = (B/2) |
The thickness of the fluid film can be expressed as: (2) where hU is the height of the upper surface, in which the load is applied, hL is the height of the lower (textured) surface, but omitting the texture contribution, and hp is the contribution to the fluid film thickness coming from the texture (see ). In the case considered here hU is the following quadratic function: (3)
The equilibrium position hU of the upper surface is found by Newton--Raphson iterations on the equilibrium condition (4) where Wa is the applied load, assumed along the y direction, and (5)
is the load capacity. The friction force F in the sliding direction (x) is calculated by integrating the shear stress along the fluid film: (6)
2.2. Dimples
The dimples we consider are square and are arranged in a square pattern. They can be characterized by few parameters, as shown in . The edgelength of the square cell that contains the dimple is denoted by ac. The edgelength of the dimple is ap, so that the area fraction covered by dimples is simply . The third parameter that characterizes a dimple is hpo. In cell coordinates (ξ , η) the shape of the dimples is given by (7) where Ω o = [Vh ;as,as+ap}]Vh ;× [as,as+ap], a s = (1/2)\, (ac-ap) and k is a “steepness parameter”.
In the optimization process, the domain Ω is divided into regions. The dimples in each region are considered a “group” (Ng). The shape parameters (ap, hpo and k) are varied keeping their values uniform for each group.
2.3. Numerical formulation
The numerical treatment of equation (1) is a quite standard Galerkin finite element formulation. Considering a discrete space Vh, which satisfies the Dirichlet conditions, one looks for a pressure field ph∈Vh satisfying (8) where ph- = min {ph,0 }, and ε is the penalization parameter, that we take as 10-10. The term involving ε models cavitation phenomena. In the limit ε → 0 it tends to satisfy the Elrod conditions p = ∂ p/∂ n = 0 at the boundary of the cavitation zone. The penalization is quadratic so as to render the left-hand side of (8) continuously differentiable with respect to ph, which is convenient for the Newton--Raphson treatment of the nonlinearity. Though cavitation does not occur in the optimized shapes found in the cases reported here, it is important to include it because along the optimization process the evolutionary algorithms may propose shapes that produce cavitation.
2.4. Optimization -- genetic algorithm
Code GENESIS Version 5.0 of J. Grefenstette (http://www.aic.nrl.navy.mil/galist/src/#C) was used as the evolutionary algorithm in this study. A brief overview of the code is as follows:
Let f:Rn → R be the objective function, and (X1,X2, ··· ,Xn) its arguments (an n-vector). The algorithm works on a population (initially random) (9) making it evolve as shown in the diagram (see ). Notice the classical steps in genetic algorithms Citation4:
• | Evaluation: Compute the objective function f. This is the most time-consuming step in our case since it involves a non-linear finite-element analysis. | ||||
• | Selection: Randomly choose those individuals that survive to the next generation, assigning more probability to those with smaller f. | ||||
• | Crossover: Randomly exchange genes (bits) between two individuals. The crossover rate can be controlled. | ||||
• | Mutation: Randomly invert bits of an individual. The mutation rate can be controlled. |
The GENESIS code allows for the mentioned rates and the population size to be adjusted so as to have a reasonable convergence rate without saturation or instability.
3. Results and discussion
3.1. Rayleigh step -- test case
As first application a case with known optimum was chosen: The so-called Rayleigh step Citation5. In this problem the objective function is the load capacity. The problem reads:
Find h*(x), h*(x) ≥ 1∀x ∈ [qb;0,1]qb; such that (10) where W[qb;h]qb; denotes the integral in (5) when p is the solution of (1). The exact optimum is given by (see ): (11)
The code GENESIS was used in three different ways. Notice that in this “validation” problem there are no dimples; we are simply optimizing the shape h(x) in the set of piecewise constants.
In Case 1 the domain [qb;0,1]qb; ; was divided into 100 subregions, the value of h in each subregion being a state variable. The variation range for each variable was set to be the interval [qb;1,2]qb;.
In Case 2 the subdivision was the same, but the first 65 variables were assigned the range [qb;1.5,2]qb;, the following 10 variables the range [qb;1,2]qb; and the remaining 25 variables the range [qb;1,1.5]qb;.
Finally, in Case 3 it was assumed that the profile had stepped shape. The state variables were, thus, the film thickness on the left of the step, hleft, the film thickness on the right, hright, and the discontinuity location, D.
The results from the first two cases after 30,000 evaluations of the objective function can be seen in , compared to the exact optimum. Results from Case 3 can be observed in , as a function of the number of evaluations. The load value obtained for Case 1, corresponding to the shape shown in is 0.9789 of the exact optimum, while that of Case 2 is 0.9921 of the exact optimum. In these cases a population size of 50 individuals was taken.
Table 1. Results of the GENESIS code compared to the exact ones, indicating the number of evaluations (Case 3).
The mutation and crossover rates used in the calculations were 0.001 and 0.6 respectively.
The above results show that in general the genetic algorithm leads to “rough” shapes, especially in regions where h is much greater than one. The lubrication equation is not very sensitive to changes in the film thickness at regions where this thickness is large. Since the “rough” shapes have loads very close to the optimal one, no driving force exists for smoothing out the shape. Of course, if the algorithm is allowed to proceed indefinitely, the exact shape is eventually attained. For Case 2 above this happens after 200,000 function evaluations, approximately, and only then the obtained shape is “smooth”. It could be interesting to design a smoothing procedure specific for partially optimized lubrication shapes. It has also been shown that adding information, as is expected, narrows the state space and allows the algorithm to get close to the optimum in fewer generations.
Another interesting point is the algorithm's behavior as a function of the number of state variables. In we show the best load value (relative to the optimal one) as attained after 30,000 evaluations of the cost function. Compare the results with different numbers of state variables Nvar = 50, 100 and 200. It is clear that the fewer the state variables, the fewer the number of iterations to get to within some given tolerance from the optimum. In all three cases the ranges of the variables were set according to Case 2 described earlier.
Table 2. Relative load capacity obtained with the GENESIS code using 50, 100 and 200 variables to define the film thickness h.
Finally, let us briefly comment about algorithmic parameters. The GENESIS code has default values for the crossover and mutation rate constants. Varying any of these constants does not have an impact on the convergence behavior of the algorithm. An example is shown in , in which the mutation rate constant is varied keeping the crossover rate constant fixed. It is observed that increasing the mutation rate constant by an order of magnitude does not affect practically the convergence behavior of the algorithm. Only after setting the constant to 100 times the default value () does the mutation inhibit convergence.
4. Textured sliders
Going now to the technological problem of texture design in the slider configuration (), in this case the objective function is f = F, the friction force, while the state variables {Xj} are texture-shape coefficients.
4.1. Infinite width case
We consider a square array of dimples, periodic in the width (z) direction. We first conducted a sensitivity analysis in order to narrow the search space, studying the variations of the friction and of the minimal spacing around the untextured state (hpo=0 everywhere) with respect to dimple depth.
Let denote the value of hpo for the Ith group of dimples. Since all other input data are fixed, the solution of the problem (i.e., the pressure field, the eccentricity, the minimal film thickness, the friction force, etc.) depends just on . We can thus calculate derivatives of any quantity, such as for example the friction force F, by numerical differentiation, that is (12)
Whether it is convenient or not to carve dimples in some region of the bearing results from the derivative of F and hmin with respect to . These derivatives are shown in . Based on them, we consider the following dimple-depth distribution: (13) so as to allow dimples only where F is expected to decrease while hmin is expected to increase.
The dimple width is set to 1/100, and the domain is divided into Ng = 20 subregions, so that each subregion contains a group of 5 dimples. From criterion (15) the search space is narrowed to groups 12–17. In these groups we set as variables the dimple depth hpo, width as, and the steepness parameter k, within the ranges The total number of variables is thus 18.
The evolution of the friction and of hmin along the evolution process is shown in for a non-dimensional applied load of 0.06. Notice that the results are shown relative to the smooth case. A significant improvement of 4% is obtained. Moreover, this improvement comes with a 2% increase in the film thickness, which is important in terms of wear. The texture corresponding to the best individual found can be seen in . Also shown in the figure is the optimized shape corresponding to a smaller applied load, W = 0.006. It is clear that the optimized textures strongly depend on the operating conditions of the device.
4.2. Finite width case
In this case we simulate half the width of the slider because of the symmetry. The cell width adopted is (1/24), so that the maximum number of dimples is 288 = 24× 12. Dimples are only allowed in the central portion of the slider (see ), since from the previous paragraph we know that this portion is the one in which dimples can be beneficial.
In we illustrate the pressure field corresponding to a uniform dimple depth of hp = 0.5. The dimples were divided into 35 groups of 2× 2 dimples each. Two free parameters were assigned to each group of dimples, the depth and the size, with the ranges: hpo ∈ (0,1.3) - a s ∈ [qb;0.0025,0.007]qb;.
The total number of parameters is thus 70. Again the genetic algorithm found a design with a 4.7% reduction in friction force. The best texture found after 3000 evaluations of the objective function is shown in .
The result after 18,000 evaluations is given in and . The reduction in friction force corresponding to this pseudo-optimal texture is 5.2/% .
5. Conclusions
The recent development of advanced surface texturing techniques, such as laser ablation, chemical etching, etc., opens the possibility of “engineered surfaces” for lubrication devices Citation2,Citation3. There exists significant activity in this direction, since important gain can be achieved in terms of friction, wear and lifetime extension.
It should not, however, be believed that some magical texture exists that is optimal no matter what the load and velocity of the device is. For each specific working condition, there exists a specific texture that improves the performance, but it depends on the load and velocity. Optimal design is thus necessary to identify the surface shape that is best for each device. It is thus appropriate to assess different optimization tools in this problem.
In this article we have focused on evolutionary algorithms, which have the potential of identifying global optima, and the results have shown that quite standard genetic procedures work rather well for the problem under study. Further investigation is however required since the cases considered do not involve crucial phenomena such as massive cavitation and viscous heating.
Notice, to conclude, that we have assumed throughout this analysis that the Reynolds approximation holds. This requires that the length scales in the sliding direction are much greater than the dimple depth. If the dimple length is comparable to its depth, the Navier--Stokes equations need to be solved to calculate the effect of the dimples and eventually optimize the design. This is the subject of ongoing research.
Acknowledgement
The authors acknowledge support from Agencia Nacional de Promoción Científica y Tecnológica (Argentina) through grants PICT'99 12-6337 and PICT'00 12-9848.
References
- Ronen Etsion, A, Etsion, I., and Kligerman, Y, 2001. Friction-reducing surface-texturing in reciprocating automotive components., Tribol. Trans 44 (2001), pp. 359–366.
- Etsion, I, and Halperin, G, 2002. A laser surface textured hydrostatic mechanical seal, Tribol. Trans 45 (2002), pp. 430–434.
- Ryk, G, Kligermand, Y., and Etsion, I, 2002. Experimental investigation of laser surface texturing for reciprocating automotive components, Tribology Transactions 45 (2002), pp. 444–449.
- Pham, D.T, and Karaboga, D, 2000. Intelligent Optimization Techniques -- Genetic Algorithms Tabu Search Simulated Annealing and Neural Networks. Springer: London; 2000.
- Rayleigh, L, 1918. Notes on the theory of lubrication, Philosophical Magazine 35 (1918), pp. 1–12.