Research on the meshing load-bearing characteristics of the micro-textured thermoelastic hydrodynamic lubrication interface for a warship PRTS

Gears are subjected to external excitation from alternating loads, and the micro-textured meshing interface (MTMI) will share larger contact stresses during the warship power rear transmission system (PRTS) torsional process. The current elastoplastic interface load-bearing contact model ignores the time-dependent changes in the textured element micro-convex peak (MCP) base diameter, which is usually regarded as a certain constant value, and which is extremely inconsistent with the time-varying characteristics of the MCP matrix diameter of actual MTMI, leading to deviation of the load-bearing analytical values determined by the current contact model from actual data. A generalised thermoelastic hydrodynamic lubrication (TEHL) contact load-bearing model with interface micro-texture (IMT) is proposed, where the contact area between all MCPs across the MTMI is represented by the equivalent scale factor parameter and the shape distribution density function is modified to ensure that the MCP is solved integrally. A mathematical model of meshing anti-scuffing load-bearing capacity (ASLBC) in a TEHL steady state is derived to reveal the correlation between contact stiffness and damping of the meshing MTMI under the influence of alternating loads, which provides a theoretical basis and data reference for homogeneous interface enriched lubrication (IEL) effect improvement and meshing ASLBC enhancement of the contact IMT for the PRTS.

1 Introduction

A micro-textured gear unit is the warship power rear transmission system (PRTS), a high-speed and heavy-duty gear transmission system for large warships, as shown in Fig. 1 (Source: National Key Research and Development Program Project, Grant No. 2023YFB3406301). The excellent lubrication and load-bearing regulation of the micro-textured meshing interface (MTMI) have been well recognised, and the teeth surface micro-texture seems to be a penalty area for gears, which is attributed to the traditional perception of similarity between meshing interface micro-pitting and texturing units (Özdemir et al., 2017; Liu et al., 2020; Guilbault and Lalonde, 2019). The gear surface micro-texturing is applied to the interface to modulate time-varying meshing stiffness and enhance the lubrication enrichment effect, thereby improving meshing load-bearing performance (Croccolo et al., 2020). The current elastoplastic interface load-bearing contact model ignores the time-dependent changes in the base diameter of the texturing micro-convex peaks (MCP), which is usually recognised as a constant value, and which is extremely inconsistent with the time-varying characteristics of the base diameter of the MCP elements in the actual MTMI, thus leading to the deviation of the load-bearing resolved values determined by the current contact model from the actual data.

Figure 1 A high-speed and heavy-duty gear transmission system for large warships (Source: National Key Research and Development Program Project, Grant No. 2023YFB3406301).

As the core component for torque transmission, gears possess diverse sizes, high load-bearing capacity, and high efficiency, which are of vital importance in determining the reliability and accuracy of the PRTS (Tu et al., 2021; Marques et al., 2021). In the case of a PRTS operating under special environments, it is confronted with extreme usage conditions such as high and low temperatures, deep-sea pressure, strong acids and alkalis, and other extreme conditions of use. The maximum pressure in the deep sea reaches over 200 to 300 MPa, which is 2000 to 3000 times that of atmospheric pressure. The marine environment brings challenges such as acid and alkali corrosion to the PRTS.

Figure 2 Contact fatigue and scuffing affecting load-bearing capacity arising from the use of deep-sea gears (Source: MPRD, Grant No. MG20220203).

The deep-sea gearbox is exposed to alternating environments of high temperature, high humidity, high salt fog and seawater spray for a long time in the alternating environment, which causes serious electrochemical corrosion problems. Deep-sea gears are subjected to contact fatigue and scuffing affecting load-bearing capacity in service, as shown in Fig. 2 (Source: MPRD, Grant No. MG20220203). High-end equipment requires gears with greater load-bearing capacity and serviceability, and controlling the load-bearing capacity of deep-sea gears has become critical. Although advanced processing technologies such as surface hardening and precision machining improve the gear load-bearing capacity, they still cannot meet the ever-increasing service life and load-bearing capacity requirements, and frequent equipment accidents caused by insufficient gear strength are common, which makes the gear load-bearing capacity become a key bottleneck restricting the performance and reliability of the PRTS.

The gear directly affects the overall performance of the PRTS. The meshing domain, where the frictional heat causes a sharp increase in teeth surface temperature, leads to the rupture of the interfacial oil film layer and consequently lubrication failure. This results in rigid contact between the metal MTMI, an increase in the friction coefficient between the teeth surfaces and a rapid rise in the instantaneous contact temperature of the meshing interfaces. The gear profile and the lubrication state between the teeth surfaces change, and the metal teeth surfaces are torn off along the high-speed sliding direction, resulting in scuffing load-bearing failure (Chen et al., 2021, 2023; Zhao et al., 2022a). Even short-term overloading, improper lubrication or excessive temperature rise causes serious damage to the teeth profile and further deteriorates the lubrication state between the meshing interfaces. Considering the harsh real-time working conditions of warship PRTSs, gear damage and failure often occur or originate from the microscopic interface of the teeth surfaces. Scuffing load-bearing fatigue of the MTMI frequently leads to gear lubrication failure. Traditional design methods, which are based on Hertzian contact theory and assumptions such as smooth surfaces, no lubrication, isothermal elastohydrodynamic lubrication and homogeneous elastic materials, are no longer capable of accurately predicting the lubrication characteristics and load-bearing performance of gears in warship PRTSs. Therefore, it is necessary to deeply explore the evolution of lubrication characteristics in heavily loaded line contact interfaces and the correlation mechanism with the load-bearing performance of MTMI. This will provide a scientific basis and theoretical reference for enhancing the thermoelastic hydrodynamic lubrication (TEHL) characteristics of the meshing interfaces under real-time full working conditions of the PRTS and for accurately assessing the anti-scuffing load-bearing capacity (ASLBC) of the MTMI.

The warship PRTS is a high-speed and heavy-load gear transmission system. The MTMI are typically in a mixed lubrication state characterised by the coexistence of boundary and thin-film fluid lubrication (Cheng et al., 2024; Sánchez et al., 2019; Černe et al., 2024). The state cannot generate a sufficient lubricating film to completely separate the MTMI. This is especially true under high dynamic operating conditions or during start–stop operations, where the characteristics of the MTMI have a significant impact on lubrication behaviour and load-bearing capacity. The meshing pair interface consists of two parts: direct contact of the MCP and the load-bearing capacity of the TEHL film (Wang et al., 2022; Chang et al., 2023; Yang et al., 2023). The contact morphology interface micro-texture (IMT) is a new idea and method for analysing the scientific issues related to gear lubrication and load-bearing capacity in warship PRTSs. The lubrication and load-bearing effects in the meshing time domain process are not negligible (Cao et al., 2024; Zhang et al., 2025). The key technology is that the appropriate IMT can generate a coupling effect of “secondary enrichment lubrication” and “micro-hydrodynamic load-bearing”. However, the heat flux generated by alternating friction of the MTMI causes a significant increase in teeth temperature, which in turn alters the stable state mechanism of the warship PRTS meshing interface TEHL (Ammosov et al., 2020; Zhao et al., 2022b). This includes multi-scale parameters such as IMT lubrication enrichment characteristics, time-varying thermoelastic deformation, transient thermal expansion, thermomechanical coupling of the stress structure, interface morphology IMT configuration characteristics and distribution density (Li et al., 2025).

In the field of warship PRTS research, an urgent and critical scientific issue that needs to be explored is how to reasonably design the IMT of the MTMI under real-time full operating conditions to improve the interface enriched lubrication (IEL) characteristics and enhance the key performance of the ASLBC. It is also essential to construct a universal mathematical model for mixed lubrication of the MTMI, simulate and solve the core variables, and verify and calibrate them through full operating condition experiments. This research aims to reveal the mechanism of the IEL effect of mixed TEHL meshing IMT and optimise the multi-scale parameter design of the microelement configuration. Additionally, it seeks to explore the influence of mechanisms and patterns of the geometric characteristics and directional parameters of the MTMI on IEL properties in order to improve the ASLBC to adapt to the high contact line loads during the alternating meshing time domain. The research on the IEL mechanism and ASLBC of the mixed TEHL meshing IMT is of great theoretical significance and practical application value for the development of the tribodynamics of the MTMI in high-speed and heavy-load gear transmission systems.

2 Steady-state TEHL textured interface load-bearing contact modelling

2.1 Constructing an interface load-bearing contact model based on fractal theory

The MTMI of the warship PRTS is in continuous contact, and the actual micro-contact area is much smaller than the macroscopic theoretical meshing domain area. This research is based on the meshing IMT fractal theory and investigates the temporal variation of single micro-convex peaks (MCPs) within the load-bearing contact domain, unveiling the evolving trends of actual contact changes in MCPs on non-macroscopically smooth interfaces. The mathematical relationship between the actual contact interface domain and the true MTMI loads is determined by solving the distribution density function of the differential MCP elements using integration. Figure 3 illustrates the mathematical model of the correlation between height deformation and individual differential MCP element protrusions on the MTMI. The MCP elements deform transiently under the cyclic excitation loads on the textured interface, considering the different heights of the MCP elements and the different loads, and the true deformation values of the MCP elements are limited within the range of 0−δ. The teeth profile of the meshing gear pair of a single MCP element within the contact interface domain is usually approximated by a cosine function (with the single base length set as l). Based on the description of the three-dimensional profile characteristics of the MTMI using a bivariate function (Codrignani et al., 2020; Wang et al., 2023; Meng et al., 2019), it is expressed as

$$\begin{array}{}\text{(1)}& \begin{array}{rl}z\left(x,y\right)& =h{\left({\displaystyle \frac{G}{h}}\right)}^{D-\mathrm{2}}\cdot {\left({\displaystyle \frac{\ln \mathit{\phi }}{M}}\right)}^{\mathrm{1}/\mathrm{2}}\sum _{m=\mathrm{1}}^M\sum _{n=\mathrm{0}}^{n_{max}}{\mathit{\phi }}^{n(D-\mathrm{3})}\\ & \left\{\cos {\mathit{\varphi }}_{m,n}-\cos \left[{\displaystyle \frac{\mathrm{2}\mathit{\pi }{\mathit{\phi }}^n\sqrt{x^{\mathrm{2}}+y^{\mathrm{2}}}}{h}}\right]\right.\\ & \left.\cdot \cos \left(\arctan {\displaystyle \frac{y}{x}}-{\displaystyle \frac{\mathit{\pi }m}{M}}\right)+{\mathit{\varphi }}_{m,n}\right\},\end{array}\end{array}$$

where z(x,y) represents the height function of the three-dimensional textured interface profile (the profile height of a single MCP element); G is the amplitude of the height of the MCP elements of the MTMI; D represents the fractal dimension; φ denotes a constant value that is not less than 1.0; and ϕ_m,n denotes a uniformly distributed random phase within the interval [0, 2π], generated by a random number generator. Its role is to ensure that the frequency components of the meshing interface profile do not spatially overlap. h is the sampling profile length of the MTMI (the diameter of the base of a single MCP element); M represents the number of superimposed ridgelets used to construct the surface profile morphology, the specific value of which can be determined based on the power spectrum of the imaged surface morphology; n is the distribution index of a random frequency of the interface profile; and n_max is the random phase when the interface profile is uniformly distributed.

Figure 3 Mathematical model of the correlation between height deformation and individual MCP element protrusions on the MTMI.

$$\begin{array}{}\text{(2)}& z\left(x\right)={\displaystyle \frac{G^{(D-\mathrm{1})}}{h^{(\mathrm{2}-D)}}}\cos \left({\displaystyle \frac{\mathit{\pi }x}{h}}\right),-h/\mathrm{2}<x<h/\mathrm{2}\end{array}$$

As illustrated in Fig. 3, when the MCPs of the textured contour are subjected to the interface normal contact pressure, the actual deformation height of the MCPs is designated as ω, the distance between the apex of the MCP and its base is defined as δ, the base diameter of the MCP is represented as h, the transverse contour diameter is denoted as h_t, the actual contour diameter is h_r, and the transverse contact area a of a single MCP is described as

$$\begin{array}{}\text{(3)}& a=\mathit{\pi }\mathit{\rho }\mathit{\omega }={\displaystyle \frac{h^{\mathrm{2}}}{\mathit{\pi }}}\left\{\mathrm{1}-\cos \left({\displaystyle \frac{\mathit{\pi }{h}_{\mathrm{t}}}{\mathrm{2}h}}\right)\right\},\end{array}$$

where ρ represents the radius of curvature of a single MCP, thereby allowing the calculation of the radius of curvature at position x=0 of the MCP (see Expression 2).

$$\begin{array}{}\text{(4)}& \mathit{\rho }={\left|{\displaystyle \frac{\sqrt[\mathrm{3}]{\mathrm{1}+{z^{\prime }}^{\mathrm{2}}}}{z^{\prime \prime }}}\right|}_{x=\mathrm{0}}={\displaystyle \frac{h^D}{G^{(D-\mathrm{1})}{\mathit{\pi }}^{\mathrm{2}}}}\end{array}$$

Based on the actual textured interface contact area and combined with related thermo-elastoplastic deformation studies (Rajput et al., 2021; Galda et al., 2019; Liu et al., 2020), it is demonstrated that when the textured interface is subjected to external alternating excitation loads, it inevitably coexists with a single MCP at the critical stage of elastic deformation. The critical deformation amount ω_l of a single MCP is expressed as (Wang et al., 2021)

$$\begin{array}{}\text{(5)}& {\mathit{\omega }}_l=\mathit{\rho }\left({\displaystyle \frac{\mathit{\pi }K{\mathit{\sigma }}_y}{\mathrm{2}E}}\right),\end{array}$$

where E represents the base material comprehensive elastoplastic modulus and $\frac{\mathrm{1}}{E}=\frac{(\mathrm{1}-{\mathit{\nu }}_{\mathrm{1}}^{\mathrm{2}})}{E_{\mathrm{1}}}+\frac{(\mathrm{1}-{\mathit{\nu }}_{\mathrm{2}}^{\mathrm{2}})}{E_{\mathrm{2}}}$ holds true, with ν_1,2 denoting Poisson's ratio. K describes the correlation coefficient between the hardness value and the yield strength σ_y.

By combining Eqs. (2) and (4), ω_l is described as

$$\begin{array}{}\text{(6)}& {\mathit{\omega }}_l={\displaystyle \frac{h^D}{G^{(D-\mathrm{1})}}}{\left({\displaystyle \frac{K{\mathit{\sigma }}_y}{\mathrm{2}E}}\right)}^{\mathrm{2}}.\end{array}$$

The critical contact area a_l of the cross-section of a single MCP is expressed as

$$\begin{array}{}\text{(7)}& a_l=\mathit{\pi }\mathit{\rho }{\mathit{\omega }}_l={\displaystyle \frac{\mathrm{1}}{\mathit{\pi }}}\left\{{\left({\displaystyle \frac{h^D}{G^{(D-\mathrm{1})}}}\right)}^{\mathrm{2}}{\left({\displaystyle \frac{K{\mathit{\sigma }}_y}{\mathrm{2}E}}\right)}^{\mathrm{2}}\right\}.\end{array}$$

According to the classical Hertz contact theory, the elastic contact load P_E(ω) of a single MCP subjected to alternating excitation loads leading to elastoplastic deformation is expressed as (Petare et al., 2019)

$$\begin{array}{}\text{(8)}& P_{\mathrm{E}}\left(\mathit{\omega }\right)={\displaystyle \frac{\mathrm{4}E\sqrt[\mathrm{3}]{\mathit{\omega }}\sqrt{\mathit{\rho }}}{\mathrm{3}}}.\end{array}$$

By combining Eqs. (2) and (4) and substituting Eq. (8), we obtain

$$\begin{array}{}\text{(9)}& P_{\mathrm{E}}\left(\mathit{\omega }\right)={\displaystyle \frac{\mathrm{4}E\sqrt[\mathrm{3}]{\mathit{\omega }}\sqrt{\mathit{\rho }}}{\mathrm{3}}}={\displaystyle \frac{\mathrm{4}E{G}^{(D-\mathrm{1})}{h}^{(\mathrm{3}-D)}}{\mathrm{3}\mathit{\pi }}}\sqrt[\mathrm{3}]{\mathrm{1}-\cos \left({\displaystyle \frac{\mathit{\pi }{h}_{\mathrm{t}}}{\mathrm{2}h}}\right)}.\end{array}$$

2.2 Single MCP load-bearing contact model

The aim is to establish the correlation between each MCP and its maximum peak value in order to achieve an integrable solution for the total elastic load of all MCPs across the entire contact area of the textured interface. Set the maximum MCP contact base diameter of the textured interface as H_max, and let H_max t represent the maximum truncated diameter of the MCP profile. By considering the cosine expression of the height profile of the MCP, the formula for the maximum value of the MCP profile is obtained as follows:

$$\begin{array}{}\text{(10)}& z\left({\displaystyle \frac{H_{max\mathrm{t}}}{\mathrm{2}}}\right)=G^{(D-\mathrm{1})}{H}_{max}^{(\mathrm{2}-D)}\cos \left({\displaystyle \frac{\mathit{\pi }{H}_{max\mathrm{t}}}{\mathrm{2}{H}_{max}}}\right).\end{array}$$

Considering the maximum micro-unit peak contact base diameter H_max of the textured interface and combining it with the calculation formula for the contact area of the micro-unit peaks, the actual bearing contact area of the textured interface is expressed as

$$\begin{array}{}\text{(11)}& a_H={\displaystyle \frac{H_{max}^{\mathrm{2}}}{\mathit{\pi }}}\left\{\mathrm{1}-\cos \left({\displaystyle \frac{\mathit{\pi }{H}_{max\mathrm{t}}^{\mathrm{2}}}{\mathrm{2}{H}_{max}}}\right)\right\},\end{array}$$

where a_H represents the contact area of the textured interface with the maximum MCP contact base diameter H_max and the contact area of the remaining micro-unit peaks is indirectly calculated using H_max t.

By introducing the contact equivalent proportionality factor ${\mathit{\epsilon }}_b=\frac{h_{\mathrm{t}}}{h}$, Expression (3) is transformed into

$$\begin{array}{}\text{(12)}& a={\displaystyle \frac{h^{\mathrm{2}}}{\mathit{\pi }}}\left\{\mathrm{1}-\cos \left({\displaystyle \frac{\mathit{\pi }{\mathit{\epsilon }}_b}{\mathrm{2}}}\right)\right\}.\end{array}$$

From the above equation, it is concluded that the contact area of micro-unit peaks with different morphological scales can necessarily be obtained by the corresponding equivalent proportionality factor ε_b, that is, $a=\frac{h^{\mathrm{2}}}{\mathit{\pi }}\left\{\mathrm{1}-\cos \left(\frac{\mathit{\pi }{h}_{\mathrm{t}}}{\mathrm{2}h}\right)\right\}=\frac{h^{\mathrm{2}}}{\mathit{\pi }}\left\{\mathrm{1}-\cos \left(\frac{\mathit{\pi }{\mathit{\epsilon }}_b}{\mathrm{2}}\right)\right\}=\frac{H_{max}^{\mathrm{2}}}{\mathit{\pi }}\left\{\mathrm{1}-\cos \left(\frac{\mathit{\pi }{H}_{max\mathrm{t}}^{\mathrm{2}}}{\mathrm{2}{H}_{max}}\right)\right\}$ holds true.

Thus, the base diameter h of each micro-unit peak is expressed as (Xing et al., 2021)

$$\begin{array}{}\text{(13)}& h=H_{max}\sqrt{\mathrm{1}-\cos \left({\displaystyle \frac{\mathit{\pi }{H}_{max\mathrm{t}}^{\mathrm{2}}}{\mathrm{2}{H}_{max}}}\right)/\mathrm{1}-\cos \left({\displaystyle \frac{\mathit{\pi }{\mathit{\epsilon }}_b}{\mathrm{2}}}\right)}.\end{array}$$

By substituting the above equation into Eq. (9), the elastic contact load P_E(ω) caused by the elastoplastic deformation of a single micro-unit peak under alternating excitation loads can be solved (Rosenkranz et al., 2019):

$$\begin{array}{}\text{(14)}& \begin{array}{rl}{P}_{\mathrm{E}}\left(\mathit{\omega }\right)& ={\displaystyle \frac{\mathrm{4}E{G}^{(D-\mathrm{1})}{h}^{(\mathrm{3}-D)}}{\mathrm{3}\mathit{\pi }}}{\left\{\mathrm{1}-\cos \left({\displaystyle \frac{\mathit{\pi }{\mathit{\epsilon }}_b}{\mathrm{2}}}\right)\right\}}^{D/\mathrm{2}}\\ & {\left\{\mathrm{1}-\cos \left({\displaystyle \frac{\mathit{\pi }{H}_{max\mathrm{t}}^{\mathrm{2}}}{\mathrm{2}{H}_{max}}}\right)\right\}}^{(\mathrm{3}-D)/\mathrm{2}}.\end{array}\end{array}$$

Let $\mathit{\chi }={\left\{\mathrm{1}-\cos \left(\frac{\mathit{\pi }{\mathit{\epsilon }}_b}{\mathrm{2}}\right)\right\}}^{D/\mathrm{2}}$, where χ represents the correlation function with the contact equivalent proportionality factor ε_b. By simplifying the above Eq. (14), we obtain

$$\begin{array}{}\text{(15)}& P_{\mathrm{E}}\left(\mathit{\omega }\right)={\displaystyle \frac{\mathrm{4}E{G}^{(D-\mathrm{1})}\mathit{\chi }}{\mathrm{3}{\mathit{\pi }}^{(D-\mathrm{1})/\mathrm{2}}}}{a}^{(\mathrm{3}-D)/\mathrm{2}}.\end{array}$$

2.3 Micro-unit peak morphology distribution density function correction model

In the alternating meshing region of the warship PRTS, the teeth surface contact area with the micro-unit peak morphology distribution cannot accurately reflect the actual IMT contact area. Based on the fractal model, the contact distribution characteristics of the MCP are described in the form of a micro-unit density function, which is expressed as

$$\begin{array}{}\text{(16)}& J\left(\mathit{\alpha }\right)={\displaystyle \frac{D{\mathit{\alpha }}_{H_{max}}^{D/\mathrm{2}}}{\mathrm{2}{\mathit{\alpha }}^{\left(D/\mathrm{2}+\mathrm{1}\right)}}},\end{array}$$

where ${\mathit{\alpha }}_{H_{max}}$ represents the contact area of the maximum micro-unit peak base diameter H_max and J(α) is the distribution density function of the non-curved surface contact area.

Considering the curvature variation characteristics of non-flat interfaces, their contact area will also change accordingly, as shown in Fig. 4. Unlike the distribution density model of non-curved surface contact areas, the curved surface contact area is not less than α, and the number of micro-unit peaks decreases, meaning the total number does not exceed ${\mathit{\alpha }}^{-D/\mathrm{2}}$. By introducing the curved micro-unit peak contact coefficient ε_Jand deriving the spherical interface contact area distribution formula J(α), the initial area distribution density function J^′(α) is further corrected (Wos et al., 2020).

Figure 4 Schematic contact morphology of microelementary peaked bodies at a nonplanar texture interface.

The curved surface microelementary peaked body contact coefficient ${\mathit{\epsilon }}_J={\left(\frac{S_L}{\sum S_n}\right)}^{K_{\mathrm{0}}}$, where S_L represents the theoretical non-flat interface contact area; here, $S_L=B_{\mathrm{1}}{\left(\frac{\mathrm{\Omega }(F,E)C}{K_{\mathrm{0}}}\right)}^{B_{\mathrm{2}}}$ holds, with B₁ and B₂ denoting parameters related to the geometric characteristics of the non-flat interface morphology. F represents the contact width of the microelement convex peak body. Ω(F,E) represents a function related to the interface material (parameter E) and its contact form (parameter F), and C represents the normal contact load between the interfaces. ∑S_n represents the total contact area of any number of non-flat interfaces, and its expression is $\sum S_n=S_{\mathrm{1}}\pm S_{\mathrm{2}}\pm \mathrm{\cdots }\pm S_n$ (where the plus sign indicates external meshing and the minus sign indicates internal meshing). K₀ represents the comprehensive curvature parameter, and $K_{\mathrm{0}}=X_{\mathrm{1}}+X_{\mathrm{2}}+X_{\mathrm{3}}+X_{\mathrm{4}}$ exists, with X₁ and X₂ denoting the principal curvatures of the normal cross-section of the upper non-flat interface and X₃ and X₄ denoting the principal curvatures of the normal cross-section of the lower non-flat interface.

3 Analysis of lubrication load-bearing characteristics considering the textured interface contact model

3.1 Analysis of contact teeth surface stiffness and damping

Based on the transient process of gear meshing in the warship PRTS, and considering the textured interface contact model, the multi-scale parameter correlation of the TEHL textured interface is determined to analyse the time-varying meshing load-bearing characteristics (i.e. meshing stiffness and viscous damping characteristics) (Sivayogan et al., 2022; Gupta et al., 2022). Figure 5 shows the geometric parameter analysis model of an involute gear with textured interface meshing contact, aiming to determine the characteristic parameters of the textured interface and the analytical parameters of lubrication properties at key meshing positions during the alternating contact process of the teeth surface. This model is used to solve the contact stiffness and viscous damping of the meshing interface.

Figure 5 Geometric parameter analysis model of textured interface meshing contact involute gear.

The time-varying key parameters during the meshing process are the meshing dynamic force, sliding velocity and radius of curvature. Point O represents any contact position on the teeth profile of the meshing surface. V₁ and V₂ are the velocities of the pinion and gear at the position of point O, respectively. Points O₁ and O₂ are the geometric centres of the pinion and gear, respectively. Line K₁K₂ represents the common tangent of the base circles. r_b1 and r_b2 are the base circle radii of the two gears, respectively. α and β are the meshing rotation angles of the two gears, respectively. The radius of curvature at the contact point O is expressed as

$$\begin{array}{}\text{(17)}& \left\{\begin{array}{l}{\mathit{\rho }}_{\mathrm{1}}=K_{\mathrm{1}}O=\mathit{\alpha }{r}_{b\mathrm{1}}=\mathit{\alpha }{K}_{\mathrm{1}}{O}_{\mathrm{1}}\\ {\mathit{\rho }}_{\mathrm{2}}=K_{\mathrm{2}}O=\mathit{\beta }{r}_{b\mathrm{2}}=\mathit{\beta }{K}_{\mathrm{2}}{O}_{\mathrm{2}}=K_{\mathrm{1}}{K}_{\mathrm{2}}-K_{\mathrm{1}}O,\end{array}\right.\end{array}$$

where the equivalent radius of curvature at the meshing point is ${\mathit{\rho }}_{\mathrm{e}}={\mathit{\rho }}_{\mathrm{1}}{\mathit{\rho }}_{\mathrm{2}}/({\mathit{\rho }}_{\mathrm{1}}+{\mathit{\rho }}_{\mathrm{2}})$. Given the rotational speeds of the two gear shafts as n₁ and n₂, the tangential velocity at the meshing point on the teeth surface is expressed as

$$\begin{array}{}\text{(18)}& \left\{\begin{array}{l}{v}_{\mathrm{1}}={\mathit{\omega }}_{\mathrm{1}}{r}_{\mathrm{1}}=\mathit{\pi }{n}_{\mathrm{1}}{r}_{\mathrm{1}}/\mathrm{30}\\ v_{\mathrm{2}}={\mathit{\omega }}_{\mathrm{2}}{r}_{\mathrm{2}}={\mathit{\omega }}_{\mathrm{2}}{r}_{\mathrm{2}}(z_{\mathrm{1}}/z_{\mathrm{2}})\mathit{\pi }{n}_{\mathrm{1}}/\mathrm{30}.\end{array}\right.\end{array}$$

Figure 6 Variation of the load distribution coefficient of the MTMI. (a) Time-varying load of the meshing interface. (b) Relation of the load distribution coefficient with the rotation angle.

Under alternating excitation loads, the meshing teeth surfaces are in a time-varying contact state, and the load borne during the transmission process originates from multiple meshing teeth. This leads to the exploration of the load distribution law among the meshing teeth surfaces (Kaneta et al., 2022). The variation curve of the load distribution coefficient between the meshing teeth surfaces of the warship PRTS is shown in Fig. 6. The load distribution coefficient at the initial point of the actual contact profile of the meshing teeth surface is the smallest, with two teeth participating in the meshing state. In the double-teeth contact (DTC) region, the load between the teeth surfaces increases and exhibits a linear distribution growth. The double-teeth meshing ends when a single tooth enters the initial meshing state. From the lowest point to the highest point of single-tooth meshing, only one pair of teeth bears the steady-state load, and the meshing node is located in the middle of the single-tooth contact (STC) region. The end of the single-tooth meshing state marks the beginning of the meshing-in of a new tooth pair (i.e. the double-teeth meshing state), and the meshing process terminates at the tooth tip. During this transient process, the load between the meshing teeth surfaces decreases and exhibits a linear trend.

Based on the proposed TEHL analytical model, the contact stiffness and damping of the meshing tooth surface interface are solved. Using the gear parameters provided in Table 1 as research data, the influence of the slip velocity and alternating load factors on the contact stiffness and damping of the meshing interface is explored. The variations in tooth surface contact stiffness and damping under different meshing interface slip velocities are analytically determined, as shown in Fig. 7.

Table 1 Basic parameters of the non-smooth contact interface.

As shown in Fig. 7, as the slip velocity at the meshing interface increases, the contact stiffness and damping of the tooth surface exhibit a decreasing trend. When the meshing region transitions from the single-tooth contact (STC) domain to the double-teeth contact (DTC) domain, the contact stiffness and damping of the tooth surface undergo a sudden step-like change. Specifically, the contact damping initially increases and then decreases, while in the meshing-in STC region, the contact stiffness shows a decreasing rather than an increasing trend. When the DTC transitions back to the STC, the contact stiffness and damping of the tooth surface sharply decrease, reaching a minimum point, after which they gradually increase and then decrease again, forming a local minimum within a specific range.

Figure 7 Time-varying contact stiffness and damping of the MTMI under different slip velocities are illustrated as follows: (a) time-varying stiffness characteristics of the MTMI and (b) time-varying damping behaviour of the MTMI.

This phenomenon can be explained as follows: when the relative slip velocity between the transient interfaces of the meshing tooth surfaces approaches 0, the hydrodynamic fluid film becomes extremely thin, and the contact area of the micro-peaks increases. As a result, the influence of the micro-peaks dominates the contact stiffness and damping. Generally, the direct contact stiffness of the tooth surface is much smaller than the stiffness of the hydrodynamic fluid film. Moreover, a decrease in the slip velocity of the tooth surface leads to an increase in the amplitude of time-varying fluctuations in the contact stiffness and damping of the meshing interface.

Figure 8 Time-varying contact stiffness and damping of the MTMI under different interface loads. (a) Time-varying stiffness of the MTMI. (b) Time-varying damping of the MTMI.

The variations in tooth surface contact stiffness and damping under different meshing interface loads are analytically determined, as shown in Fig. 8. From the figure, it can be observed that as the meshing interface load increases, the tooth surface contact stiffness and damping also increase. Similar to the influence of slip velocity on tooth surface contact stiffness and damping, the time-varying trends of both stiffness and damping are nearly identical under varying loads. Compared to the effects of slip velocity, the influence of different meshing interface loads on tooth surface contact stiffness and damping is more consistent, with similar amplitude fluctuations. This indicates that the distribution characteristics of the hydrodynamic fluid film and its thickness dominate the values of tooth surface contact stiffness and damping. Furthermore, within the meshing domain, the distribution of hydrodynamic film pressure and film thickness has a more significant impact on slip velocity.

3.2 Analytical solution for tooth surface meshing stiffness and interface contact damping

Based on the aforementioned numerical analysis of tooth surface contact stiffness and damping, combined with the initial contact tooth surface meshing stiffness and damping, the tooth surface meshing stiffness and interface contact damping under TEHL conditions are further solved. The influence of slip velocity and alternating load on the variation trends of tooth surface meshing stiffness and interface contact damping is analysed. The time-varying patterns are illustrated in Fig. 9. As shown in the figure, when the slip velocity increases, the interface meshing stiffness and damping decrease. The variation amplitude of the interface meshing stiffness is less pronounced compared to the significant fluctuations in damping. Additionally, abnormal fluctuations in amplitude are observed during the transition between STC and DTC, or at the initial meshing-in and meshing-out stages. The transition from DTC to STC intensifies the time-varying damping, and during the alternating stages of STC and DTC, the slip velocity decreases, accompanied by a reduction in damping and large amplitude fluctuations.

Figure 10 illustrates the time-varying relationship between the interface meshing stiffness and damping under different tooth surface loads.

Figure 9 Meshing time-varying stiffness and damping of different slip velocities. (a) Time-varying interface meshing stiffness. (b) Time-varying interface meshing damping.

Figure 10 Meshing time-varying stiffness and damping of different loads. (a) Time-varying interface meshing stiffness. (b) Time-varying interface meshing damping.

As shown in the figure, as the load on the tooth surface increases, the time-varying meshing stiffness and damping of the interface also increase. Compared to the results of the slip velocity analysis, the variation amplitude of the meshing stiffness and damping exhibits larger and more pronounced fluctuations. Particularly in Fig. 10a, when the meshing region transitions from the DTC to the STC, a sudden change in meshing stiffness occurs.

This change becomes more significant as the tooth surface load decreases, which is attributed to the reduction in meshing stiffness under TEHL conditions. Additionally, meshing damping is significantly influenced by the tooth surface load. During the alternating transition between STC and DTC (or at the initial meshing-in and meshing-out stages), the amplitude of damping shows abnormal fluctuations, as depicted in Fig. 10b. When transitioning from DTC to STC, the interface meshing damping undergoes a sudden change, and the amplitude variation trends under different teeth surface loads are almost similar. Overall, the time-varying patterns of interface meshing damping under different tooth surface loads are generally consistent, with only slight differences observed in the time-varying trends within the DTC region.

Under TEHL conditions, modifications such as profile fitting optimisation or edge trimming are often performed along the meshing length direction of the tooth surface interface to reduce and suppress the stress concentration. Additionally, considering the finite length of the actual contact line on the meshing tooth surface, these meshing characteristics will influence the analysis of the contact interface ASLBC.

4 Simulation analysis of the influence of texture parameters on the load-bearing performance of the MTMI

Considering the characterisation parameters of surface textures, the equations for contact load and contact area for the MTMI are derived. The influence of material parameters, friction coefficients and peak amplitudes of MCP bodies on the normal stiffness, viscous damping and ASLBC of the meshing interface is analysed. The intrinsic relationship between the microscopic morphology of textures and the load-bearing characteristics of meshing tooth surfaces is established. Based on the derived mathematical model of textured interface contact and the simulation using MATLAB, the time-varying trends of the ASLBC under different texture parameters are revealed. This approach aims to determine the optimal texture parameters, thereby guiding the design of micro-textured patterns on meshing teeth surfaces to enhance the operational stability and service life of warship PRTSs under alternating and complex working conditions.

4.1 Simulation analysis of the influence of fractal dimension on meshing load-bearing performance

Considering different fractal dimensions D and using MATLAB simulations, the relationship between the textured tooth surface load and the actual interface contact area within the meshing region is analysed.

Figure 11 Correlation law between textured surface loads and actual interfacial contact area in meshing domains with different fractals. (a), (b) Effect of actual interface contact area on load characteristics variation of textured surface. (c), (d) Time-varying trend of elastomer share of contact area in meshing domains.

As shown in Fig. 11a, under the same meshing load, as the contact interface area increases, the unit micro-area decreases, and the load-bearing capacity of the tooth surface improves. The results indicate a proportional linear relationship between the textured tooth surface load and the actual interface contact area. The analysis also reveals that as the fractal dimension D increases, the slope of the variation curve gradually rises and reaches its maximum at fractal dimension D=1.50, thereby demonstrating that the textured tooth surface exhibits optimal ASLBC at this fractal dimension. From Fig. 11b, it can be observed that at fractal dimension D>1.50, if the slope of the variation curve shows a declining trend, it indicates the existence of an optimal constant value for this fractal dimension. Specifically, when the fractal dimension of the interface texture is D=1.50, the meshing tooth surface achieves the best ASLBC.

From Fig. 11c and d, it can be concluded that as the contact area of the textured interface increases, the proportion of elastic microelements in the contact area of the meshing region gradually rises. However, as the fractal dimension increases, the proportion of elastic elements in the contact area shows a declining trend. This decline is particularly significant at fractal dimension D<1.50, further illustrating that an increase in fractal dimension leads to a sharp rise in the number of plastic deformations caused by the peaks of the microelements.

4.2 Simulation analysis of the influence of the interface friction coefficient on meshing load-bearing performance

The time-varying influence of the interface friction coefficient on the load-bearing characteristics of the contact area is not significant, as shown in Fig. 12. Based on traditional contact mechanics theory, the tangential contact load of the meshing tooth surface is significantly constrained by the interface friction coefficient, whereas the normal contact load is less affected by it. Therefore, optimising the interface friction coefficient does not significantly improve the meshing load-bearing performance.

Figure 12 Effect of time-varying interfacial friction coefficients on contact area load-bearing characteristic curves law. (a) Effect of actual interface contact area on load characteristics variation of textured surface. (b) Time-varying trend of elastomer share of contact area in meshing domains.

5 Simulation analysis of the influence of micro-texture-scale parameters on meshing load-bearing performance

Based on the numerical simulation of the tooth surface meshing model, the TEHL performance of the textured interface under sliding contact in warship PRTS gears is analysed for different micro-texture-scale parameters. The characteristics of the medium layer film, friction coefficient and medium film pressure at the meshing interface are examined in relation to the geometric parameters of the texture morphology, such as pit radius, depth, peak height of micro-convex bodies and width-to-diameter ratio of concave–convex microelements. The time-varying trends of these parameters are analysed to determine the mechanism by which micro-texture-scale parameters influence the lubrication properties and load-bearing performance of the MTMI. This analysis aims to guide the micro-texture topology design and functional optimisation of the contact interface in warship PRTS meshing teeth surfaces.

5.1 Simulation analysis of the influence of pit microelement radius on load-bearing performance

Figure 13a illustrates the time-varying relationship of the minimum medium layer film thickness ratio at the interface under different pit microelement radii during the sliding contact transient of textured meshing teeth surfaces. As shown in the figure, the time-varying curve of the minimum medium layer film thickness ratio at the MTMI exhibits a gradually decreasing trend. The rate of change slows down at the initial stage of TEHL meshing and then accelerates until the end of the meshing process. Comparing the time-varying curves of the minimum medium layer film thickness ratio under different pit microelement radii, it is evident that when the depth of the pit microelements and the height of the convex peak microelements are held constant, an increase in the pit microelement radius (PMER) leads to a gradual increase in the minimum medium layer film thickness ratio at the MTMI. This is particularly evident under TEHL conditions, where an increase in the PMER promotes the accumulation of hydrodynamic fluid within the textured microelements, resulting in an increase in the medium layer film thickness. However, the rate of increase diminishes significantly as the PMER further increases.

Figure 13 Bearing characteristics of the pit microelement radius. (a) Time-varying trend of the interface minimum dielectric layer film thickness ratio for different pit microelement radii. (b) Time-varying trend of the dielectric layer film pressure and pressure of microelement peaks for different pit microelement radii. (c) Cloud view of the interfacial dielectric layer film pressure for different pit microelement radii.

Figure 13b depicts the time-varying relationships of the hydrodynamic medium layer film pressure and the MCP pressure within the meshing time domain. The results indicate that the hydrodynamic medium layer film pressure initially increases gradually and then decreases sharply during the meshing process. In contrast, the MCP pressure increases rapidly within the steady-state meshing region. This behaviour is attributed to the gradual reduction of the meshing gap as the steady-state meshing transient progresses, with the medium layer film pressure playing a dominant role. As the medium layer film pressure increases, the MCP pressure initially emerges and then rises sharply, while the medium layer film pressure gradually diminishes. Comparing the two pressures, it is observed that as the PMER increases, the initial emergence of both the medium layer film pressure and the MCP pressure is progressively delayed. This can be explained by the fact that larger pit microelements store more lubricant, which provides sufficient medium layer film pressure to meet the normal load requirements of the TEHL interface during the meshing transient. Further increasing the PMER does not significantly improve the lubrication and friction performance but may instead inhibit the torque transmission capacity of the meshing pair. Specifically, when the PMER exceeds a certain threshold, the stored lubricant within the pit microelements increases the medium layer film pressure, leading to a gradual reduction in the peak pressure of the micro-convex elements.

In Fig. 13c, an appropriate PMER can increase the hydrodynamic medium layer film thickness to some extent and improve lubrication and friction performance. However, an excessively large PMER results in the stored lubricant within the texture still bearing a portion of the normal load at the end of the meshing transient, thereby reducing the frictional torque at the interface and affecting the torque transmission capacity.

In summary, the results indicate that an appropriate pit microelement configuration can effectively reduce the interface friction coefficient under TEHL conditions during the meshing transient, thereby minimising sliding friction and contact losses. Based on the scale factor of the PMER, if the PMER is too small, the moment at which the medium layer film pressure decreases occurs earlier, weakening the normal load-bearing capacity provided by the medium layer film and exacerbating interface friction and wear. Conversely, if the PMER is too large, it increases the shear friction coefficient of the medium layer film during the initial meshing stage, leading to higher frictional energy losses. Additionally, at the end of the meshing transient, the medium layer film provides a certain load-bearing capacity, which reduces the peak contact pressure of the MCP and negatively impacts torque transmission.

5.2 Simulation analysis of the influence of pit microelement depth on load-bearing performance

Assuming the PMER and the height of the MCP texture are constant, the time-varying relationships of the minimum medium layer film thickness ratio, medium layer film pressure, and MCP texture pressure during the alternating contact transient of meshing teeth surfaces under different pit microelement depths are illustrated in Fig. 14.

Figure 14 Bearing characteristics of the pit microelement depth. (a) Time-varying trend of the interface minimum dielectric layer film thickness ratio for different pit microelement depths. (b) Time-varying trend of the dielectric layer film pressure and pressure of microelement peaks for different pit microelement depths. (c) Cloud view of the interfacial dielectric layer film pressure for different pit microelement depths.

As shown in Fig. 14a, when the pit microelement depth increases, the minimum medium layer film thickness ratio at the interface exhibits a gradually decreasing trend, which is opposite to the effect of the PMER. An increase in the PMER prolongs the contact transient during which the micro-convex peak texture bears pressure. This trend is similar to the influence of the PMER on the medium layer film thickness ratio. Specifically, a larger PMER allows for the accumulation of more hydrodynamic fluid, enhancing the load-bearing capacity of the medium layer film. Consequently, this delays the time domain during which the MCP texture bears the load.

As shown in Fig. 14b and c, when the meshing tooth surface operates under TEHL conditions and the PMER is relatively small, the medium layer film pressure is not significantly affected by changes in the pit microelement depth. This can be attributed to the fact that, under TEHL conditions, the medium layer film thickness at the meshing interface provides sufficient load-bearing capacity, thereby mitigating the influence of the PMER on the minimum medium layer film pressure. Consequently, the effect of the pit microelement depth becomes more pronounced. Additionally, as the PMER increases, the influence of the pit microelement depth on the pressure of the MCP texture gradually diminishes.

5.3 Simulation analysis of the influence of micro-convex peak texture height on load-bearing performance

Assuming the PMER and depth are constant, the time-varying relationships among the minimum medium layer film thickness ratio, medium layer film pressure, and micro-convex peak (MCP) texture pressure during the alternating contact transient of meshing tooth surfaces under different MCP texture heights are illustrated in Fig. 15. As shown in Fig. 15a, as the height of the MCP texture increases, the minimum medium layer film thickness ratio at the interface exhibits a gradually decreasing trend. This is attributed to the fact that the MCP texture increases the surface roughness of the tooth, leading to a reduction in the minimum medium layer film thickness ratio. Excessively high MCP textures can cause the tooth surface to prematurely enter TEHL contact at higher sliding speeds during the meshing transient, thereby exacerbating friction and wear. When the PMER is very small, excessively high MCP textures can also result in an excessively low minimum medium layer film thickness ratio at the interface during the meshing transient, making it difficult to meet the required medium layer film load-bearing capacity under TEHL conditions.

Figure 15 Bearing characteristics of pits with microelement convex peak height. (a) Time-varying relationship of interface minimum dielectric layer film thickness ratio. (b) Time-varying trend of dielectric layer film pressure and pressure of microelement peaks. (c) Cloud view of interfacial dielectric layer film pressure.

As shown in Fig. 15b and c, when the height of the MCP texture increases, the pressure of the MCP emerges prematurely, thereby activating the anchoring effect of the MCP to enhance frictional contact and increase the friction coefficient at the MTMI. This is intended to achieve greater torque transmission capacity. When the PMER increases, the pressure of the MCP is delayed, indicating that optimising the matching values of the MCP texture height and radius through coordinated control can improve the torque transmission capacity of the meshing pair. Additionally, appropriately increasing the height of the MCP texture can mitigate the influence of the PMER scale on the peak pressure of the micro-convex texture. As the height of the MCP texture continues to increase, the rate of decline in the maximum peak pressure difference gradually stabilises.

Future research could employ microbial machining to fabricate linear textured surfaces on meshing interfaces, thereby validating the performance of microelement IEL and ASLBC in slip-line contact textured gear pairs. This would provide experimental verification for the theoretical model (Ruan et al., 2022; Wang et al., 2023, 2025).

6 Discussion of the results

This study is the first to introduce microelement texture (MET) into the homogeneous interface enriched lubrication (IEL) regulation of the micro-textured meshing interface (MTMI) of warship power rear transmission systems (PRTSs) to control meshing anti-scuffing load-bearing capacity (ASLBC). The proposed novel lubrication and load-bearing evaluation method is original and is expected to break through the limitations of heavy-load gear fatigue and scuffing failure in the warship PRTS, making it immune to the exacerbation of friction and wear caused by changes in lubrication characteristics at the meshing interface. Guided by the actual needs and problems of the warship PRTS gears, this research proposes a new method for evaluating ASLBC with interface micro-texture (IMT) based on the regulation of homogeneous IEL characteristics on the meshing tooth surface. This method has the potential to form a systematic and independent intellectual property system, providing a theoretical foundation and analytical tools for the study of structural design, processing methods, transmission efficiency, lubrication characteristics, fatigue life and reliability of existing warship PRTS gears.

7 Conclusions

This topic investigates the load-bearing characteristics of textured interfaces under TEHL conditions. A generalised contact load-bearing model for MTMI under TEHL is proposed, incorporating an equivalent contact scale factor. This factor represents the contact area of all micro-convex peaks (MCPs) on the meshing tooth surface using the peak with the highest micro-convex body and modifies the morphological distribution density function to ensure integrability for micro-convex bodies. A distribution density function for the MET in the contact area of the meshing interface is established, enabling more accurate analysis of the actual contact area of the MTMI. A mathematical model for the load-bearing performance of textured meshing interfaces under steady-state TEHL conditions is derived, revealing the relationship between load and the contact stiffness of textured teeth surfaces as well as the damping of the meshing interface. This provides a theoretical foundation and data reference for optimising the IEL effect of MTMI and enhancing the ASLBC of the contact IMT. Based on the research findings, the following conclusions are drawn:

1. The meshing interface load and sliding velocity govern the tooth surface contact stiffness and interface contact damping with similar time-varying patterns. The distribution characteristics of the hydrodynamic fluid medium layer and its film thickness dominate the values of tooth surface contact stiffness and interface contact damping. Meanwhile, the distribution of hydrodynamic fluid film pressure and the layer thickness have a more pronounced influence on the sliding velocity within the meshing region.

2. An increase in the fractal dimension leads to a sharp rise in the number of plastic deformations caused by the peaks of the MET. As the fractal dimension increases, the proportion of elastic elements in the contact area shows a declining trend.

3. Under TEHL conditions with a small pit microelement radius (PMER), the medium layer film pressure is minimally affected by changes in the pit microelement depth. The medium layer film thickness at the MTMI provides sufficient load-bearing capacity, thereby reducing the influence of the PMER on the minimum medium layer film pressure. Consequently, the effect of the pit microelement depth becomes more pronounced. Additionally, as the PMER increases, the influence of the pit microelement depth on the pressure of the micro-convex peak texture gradually diminishes.

This study has the following limitations:

1. The experimental validation remains incomplete, lacking supporting data from actual working conditions of the warship power rear transmission system (PRTS).

2. The research methodology primarily relies on theoretical derivation, and empirical studies need further reinforcement.

For these issues, this study proposes a generalised thermoelastic hydrodynamic lubrication (TEHL) contact load-bearing model with interface micro-texture (IMT). The model employs an equivalent scale factor parameter to characterise the global micro-asperity contact area of the micro-textured meshing interface (MTMI) and ensures the accuracy of micro-asperity calculations by modifying the shape distribution density function. Building upon this foundation, a mathematical model for the anti-scuffing load-bearing capacity (ASLBC) of meshing interfaces under steady-state TEHL conditions is derived. This systematically elucidates the coupling mechanism between contact stiffness and damping of meshing MTMI under alternating loads.

The research findings will help address the limitations of previous studies on the time-varying friction TEHL characteristics of warship PRTSs, which did not fully consider the micro-topography of the rolling/sliding transition points/lines in the contact interface. The study explores how to regulate the multi-scale characterisation parameters of interface micro-texture (IMT) to achieve optimal interface lubrication characteristics and load-bearing performance. This breakthrough overcomes the poor effectiveness of traditional lubrication under harsh operating conditions. Gear surfaces with MTMI maintain a low-friction state at the contact interface for an extended period. This provides the possibility for enriched lubrication and anti-scuffing load-bearing capacity of meshing gear surfaces under high-speed and heavy-load extreme operating conditions in warship PRTSs. Investigating the evolution mechanism of IEL characteristics under the influence of multi-scale characterisation parameters, and coupling the stress field, film thickness field and temperature field with macroscopic parameters such as configuration and load distribution, as well as microscale parameters such as micro-asperity peaks and texture direction, will provide theoretical and technical support for further improving the ASLBC of meshing gear surfaces in warship PRTSs. This is a research work of practical significance.