Hydrodynamic Bearing Analysis Back
Tribology

Hydrodynamic Bearing Analysis

Visualize journal bearing oil film pressure distribution and shaft locus in real time. Compute Sommerfeld number, minimum film thickness, and friction coefficient to optimize bearing design.

Bearing Parameters
Shaft Radius R (mm)50
Radial Clearance c (μm)100
Eccentricity Ratio ε0.50
Viscosity μ (mPa·s)46
Speed N (rpm)3000
Load W (kN)10.0
L/D Ratio1.00
Results
Sommerfeld No. S
Min. Film h_min (μm)
Attitude Angle φ (°)
Friction Coeff. f

Theory (Short-Bearing Approx.)

$$S = \frac{\mu N}{P}\left(\frac{R}{c}\right)^2$$ $$h(\theta) = c(1 + \varepsilon\cos\theta)$$ $$h_{\min}= c(1-\varepsilon)$$ $$\tan\phi = \frac{\pi\sqrt{1-\varepsilon^2}}{4\varepsilon}$$

What is Hydrodynamic Bearing Analysis?

🧑‍🎓
What exactly is the "Sommerfeld number" that this simulator calculates? It sounds complicated.
🎓
Basically, it's a single number that tells you if your bearing will work smoothly or fail. It combines all the key parameters you see in the simulator—viscosity (μ), speed (N), load (W), radius (R), and clearance (c). A high Sommerfeld number means a thick, protective oil film. Try lowering the load (W) slider and watch the Sommerfeld number increase; you'll see the shaft center move closer to the bearing center, indicating more stable operation.
🧑‍🎓
Wait, really? So the eccentricity ratio (ε) shows how off-center the shaft is. If ε is 0, does that mean it's perfectly centered?
🎓
Exactly! When ε = 0, the shaft is perfectly centered in the bearing, and the oil film thickness is uniform all around. In practice, that only happens with no load or very high speed. For instance, in a real car engine bearing at startup (low speed, high load), ε is high—the shaft is pushed almost against the bearing wall. Drag the eccentricity slider to 0.9 in the simulator and see how the minimum film thickness (h_min) becomes tiny, risking contact.
🧑‍🎓
So the oil film pressure is what actually lifts the shaft? How does changing the L/D ratio affect that?
🎓
Yes, the high pressure in the converging wedge of oil generates the lifting force. The Length-to-Diameter (L/D) ratio is crucial. A "shorter" bearing (low L/D, like 0.5) lets oil escape from the sides more easily, which typically results in lower peak pressure but a different stability profile. A common case is a crankshaft main bearing, often with L/D around 0.4. Change the L/D control here and observe how the pressure distribution plot changes shape—it becomes less peaked with a lower ratio.

Physical Model & Key Equations

The core of hydrodynamic lubrication is the Reynolds equation. For a "short bearing" approximation (length much less than diameter), it simplifies to describe how pressure (p) builds up along the bearing circumference (θ) and length (z). The solution gives the pressure distribution that supports the load.

$$\frac{\partial}{\partial z}\left(h^3 \frac{\partial p}{\partial z}\right) = 6\mu U \frac{\partial h}{\partial \theta}$$

Here, $h$ is the local oil film thickness, $\mu$ is dynamic viscosity, and $U$ is the surface speed of the journal. This equation shows that pressure generation ($\partial p/\partial z$) depends on the changing film thickness ($\partial h/\partial \theta$)—that's the converging "wedge" that pumps oil.

The performance is summarized by the dimensionless Sommerfeld Number (S). It's derived from scaling the Reynolds equation and is the primary design parameter.

$$S = \frac{\mu N}{P}\left(\frac{R}{c}\right)^2$$

Where $P = W / (L \cdot D)$ is the average bearing pressure, $N$ is rotational speed (rps), $R$ is journal radius, and $c$ is radial clearance. A high S (>1) indicates thick-film, stable lubrication. A very low S (<0.01) suggests boundary lubrication with potential for wear.

Real-World Applications

Internal Combustion Engines: Crankshaft main and connecting rod bearings are classic hydrodynamic bearings. At high RPM, they operate with a high Sommerfeld number and a stable oil film. At startup (low N, high W), the simulator shows S plummets, explaining why most engine wear occurs during cold starts.

Power Generation Turbines: Large steam or gas turbine rotors are supported by journal bearings. The simulator's L/D parameter is critical here; turbines often use L/D ratios near 1.0 to maximize load capacity and control vibration, requiring precise analysis of the pressure field and shaft locus.

Marine Propulsion Shafts: The propeller shaft on a ship is supported by massive stern tube bearings. They handle immense loads and must tolerate slow speeds during maneuvering. Engineers use analysis like this to select the right viscosity (μ) and clearance (c) to ensure the film doesn't collapse under low-S conditions.

High-Speed Machine Spindles: In precision machining, spindle bearings must run with extreme accuracy and minimal vibration. By designing for a specific Sommerfeld number and eccentricity ratio (ε), engineers can position the shaft in a region of high stiffness and damping, as visualized by the shaft locus plot in the simulator.

Common Misunderstandings and Points to Note

When starting to use this tool, there are several points where beginners to CAE, in particular, often stumble. The first is "realistic ranges for input parameters". For example, if you set the clearance C to an extremely small value like "1μm" based purely on theory, the simulation might show a thick h_min, but in reality, considering machining accuracy and thermal expansion, contact would occur immediately. For general industrial machinery, a realistic starting point is around 0.1% of the journal diameter (e.g., about 50μm for a 50mm shaft).

Next is the "interpretation of the Sommerfeld number S". It's easy to think that a larger S is safer, but an excessively large S (e.g., S>10) is also problematic. If the oil film is too thick, shaft vibration (oil whip) becomes more likely, leading to rotational instability. Stable fluid lubrication often occurs with S in the range of 1 to 3, which represents a "safe but not excessive" condition.

Finally, there's the pitfall of "overlooking the temperature dependence of viscosity η". The tool calculates using a constant viscosity, but in actual machinery, friction heat increases the oil temperature, significantly reducing viscosity. For instance, oil with a viscosity of 0.03 Pa·s at 40°C can drop to less than half its value at 80°C inside the bearing. Even if the calculation result is OK, if the oil film breaks down during actual operation, you need to recalculate using the viscosity at the anticipated operating temperature.

Related Engineering Fields

The simulation technology for this journal bearing is actually applied as a foundation in various engineering fields. The most strongly related field is "Tribology". This is the science of friction, wear, and lubrication. Data on oil film pressure and thickness obtained here are utilized in evaluating bearing surface coating technologies and the effects of lubricants with additives.

Next is "Vibration Analysis of Rotating Machinery". The oil film in a journal bearing has spring and damper-like properties (oil film stiffness and damping), which directly affect the critical speed and vibration modes of the rotating shaft. The eccentricity state calculated by this tool becomes an important input parameter for evaluating the shaft's dynamic stiffness.

Furthermore, there is expansion into "Thermal Fluid Analysis (CFD)". While this tool uses one-dimensional approximate calculations, if you want to know the detailed flow of oil inside the bearing, heat generation, and temperature distribution, a three-dimensional CFD simulation is necessary. In such cases, the calculation results from this tool are very useful as boundary conditions or validation data for the CFD model. Also, in the field of Mechatronics, the design philosophy for hydrostatic gas bearings in ultra-precision positioning stages is also based on this fluid lubrication theory.

For Further Learning

If you feel confident with this tool's results and want to take the next step, we recommend following these stages. First, understand the "formulation behind the Raimondi-Boyd chart". The tool internally uses approximation formulas like $$ \varepsilon = \frac{0.21394}{ (S + 0.3)^{0.5} } $$ (the actual ones are more complex) to find ε from S. By tracing this relational formula back to the Reynolds equation $$ \frac{\partial}{\partial x}\left(h^3 \frac{\partial p}{\partial x}\right) + \frac{\partial}{\partial z}\left(h^3 \frac{\partial p}{\partial z}\right) = 6\mu U \frac{\partial h}{\partial x} $$, you can grasp the physical essence of "why pressure is generated by shaft rotation".

Next, consider unsteady states. Actual machinery inevitably passes through the dangerous boundary lubrication region during startup and shutdown. To evaluate wear during this "transient state", you need to learn about the Stribeck curve and the concept of cumulative damage evaluation using it. Furthermore, in actual design, the ultimate goal is often to determine the bearing's stiffness and damping coefficients, integrate them with rotor dynamics analysis, and optimize the design including vibration characteristics. After getting a feel for the parameters by experimenting with the tool, progressing to these more systematic theories will significantly broaden your design perspective.