3D Finite Element Model Research Articles

Prediction and control of profile for silicon steel strip in the whole tandem cold rolling based on PSO-BP algorithm

In the rolling process of silicon steel strip of the five-stand UCMW tandem cold rolling line, the quality of finished profile is the result of cooperative control of five stands, the exit profiles of the first four stands jointly decide the exit profile of the fifth stand, so it is necessary to accurately control the strip profile in the rolling process. However, the shapemeter is only typically equipped at the exit of fifth stand, making it impossible to obtain the strip profile of the first four stands. A profile prediction model for the whole tandem cold rolling based on PSO-BP algorithm is established, compared with the finite element method, this model can satisfy the demand for fast prediction of strip profile and quickly analyzing the profile change characteristics for five stands according to the adjustment of profile control parameters in the industrial field. Firstly, a 3D elastic-plastic finite element model of rolls and strip for five stands is built, and a large number of profile calculation data are as the training sample sets to build the PSO-BP model. Through the transfer calculation of strip profile between adjacent stands, the PSO-BP model of five stands is established, and the validity is verified by comparing with the actual measured data. Secondly, to improve the quality of product profile, focusing on the unreasonable setting of production parameters, many improvement schemes of WRB and WRS for five stands are designed, and the PSO-BP model is used to quickly calculate the strip profile, then comparing the result of each scheme, the optimum improvement scheme is selected, which is furtherly carried out rolling experiments in the industrial field, and the quality of strip profile is well improved. At present, the new scheme has been successfully applied to industrial production.

Biomechanical behavior of implant retained prostheses in the posterior maxilla using different materials: a finite element study.

The aim of this study is to evaluate the biomechanical behavior of the mesial and distal off-axial extensions of implant-retained prostheses in the posterior maxilla with different prosthetic materials using finite element analysis (FEA). Three dimensional (3D) finite element models with three implant configurations and prosthetic designs (fixed-fixed, mesial cantilever, and distal cantilever) were designed and modelled depending upon cone beam computed tomography (CBCT) images of an intact maxilla of an anonymous patient. Implant prostheses with two materials; Monolithic zirconia (Zr) and polyetherketoneketone (PEKK) were also modeled .The 3D modeling software Mimics Innovation Suite (Mimics 14.0 / 3-matic 7.01; Materialise, Leuven, Belgium) was used. All the models were imported into the FE package Marc/Mentat (ver. 2015; MSC Software, Los Angeles, Calif). Then, individual models were subjected to separate axial loads of 300N. Von mises stress values were computed for the prostheses, implants, and bone under axial loading. The highest von Mises stresses in implant (111.6MPa) and bone (100.0MPa) were recorded in distal cantilever model with PEKK material, while the lowest values in implant (48.9MPa) and bone (19.6MPa) were displayed in fixed fixed model with zirconia material. The distal cantilever model with zirconia material yielded the most elevated levels of von Mises stresses within the prosthesis (105MPa), while the least stresses in prosthesis (35.4MPa) were recorded in fixed fixed models with PEKK material. In the light of this study, the combination of fixed fixed implant prosthesis without cantilever using a rigid zirconia material exhibits better biomechanical behavior and stress distribution around bone and implants. As a prosthetic material, low elastic modulus PEKK transmitted more stress to implants and surrounding bone especially with distal cantilever.

Full-Scale Model Tests of Two Box-Type Soil-Steel Structures with Different Crown and Haunch Radii.

Compared with circular, arched, and pipe-arched soil-steel structures, box-type soil-steel structures (BTSSSs) have the advantages of high cross-section utilization and low cover depth. However, the degree of influence of the crown and haunch radii on the mechanical performance of BTSSSs is still unclear. Therefore, two full-scale BTSSS models with a span of 6.6 m and a rise of 3.7 m but with different crown and haunch radii were established, and the mechanical properties during backfilling and under live load were tested. Afterward, 2D finite element models (FEMs) were established using the ABAQUS 2020 software and verified using the test data. The influence of cross-section geometric parameters on mechanical performance was analyzed by using the FEM, and a more accurate formula for calculating the bending moment during backfilling was proposed. The results show that the BTSSS with a smaller crown radius has a stronger soil-steel interaction, which promotes more uniform stress on the structure and makes the structure have smaller relative deformations, bending moments, and earth pressure. The span and arch height greatly influence the bending moment and deformation of the structure. Based on the CHBDC, the crown and haunch radii were included in the revised calculation formula.

Effect of matrix material property on the composite tibia fracture plate: a biomechanical study.

For the purpose of fixing tibia fractures, composite bone plates are suggested. Metal plates cause stress shielding, lessen the compression force at the fracture site, and have an impact on the healing process because they are significantly more rigid than bone. To prevent excessive shear strain and consequent instability at the fracture site, it is imperative to reduce stiffness in the axial direction without lowering stiffness in the transverse direction. Only a carefully crafted fiber reinforced composite with anisotropic properties will suffice to accomplish this. The purpose of the current study is to examine the impact of axial and shear movements at the fracture site on the fixing of metal and composite bone plates. After modeling the tibia with a 1 mm fracture gap, titanium plates, carbon/epoxy, carbon/PEEK, and carbon/UHMWPE composite bone plates were used to fix it. There are 6 holes on each of the 103 mm long plates. To determine the stresses and axial movement in the fracture site, anatomical 3D Finite Element (FE) models of the tibia with composite bone plates are built. The simulations that were run for various composite plate layouts and types give suggestions for selecting the best composite bone plate. Although the matrix material causes some variations in behaviors, most of the plates perform as well as or even better than metal plates. Thus, the appropriate composite combinations are recommended for a given fracture structure.

Effect of Case Depth and Hardness Distribution on the Rolling Contact Fatigue Performance of G20CrNi2MoA Carburized Steel

This study aims to investigate the effects of the case depth and hardness distribution of carburizing on the rolling contact fatigue (RCF) properties of G20CrNi2MoA steel under deeper carburizing conditions. Four groups of specimens of different case depths are prepared and subjected to RCF tests. A 2D finite element model based on continuum damage mechanics is developed to further analyze these two factors. As the depth of the carburized case increases, the RCF life first increases and then decreases. The optimum case depth is approximately one‐quarter of the thickness of the specimen. A gentler hardness distribution due to a change in the diffusion process parameters ensures a minimal increase (8%) in the RCF performance. The proper case depth has a more notable effect on the RCF life compared with those of the gentle or steep hardness distribution. It is expected that the results can help to promote a better understanding of the heat treatment process of carburized bearings and ensure more reliable carburized bearings are produced.

Investigation of the Seismic Performance of Historical Masonry Aqueducts with Eulerian-Lagrangian Coupled Analysis

ABSTRACT Throughout history, aqueducts have been built to exceed obstacles on the waterway while preserving the energy of water. Our priceless cultural heritage, aqueducts must be protected to be passed on to future generations. Filling the channel with water can increase earthquake resistance due to the damping effect of water. Investigating of fluid and structure interaction of aqueducts through numerical analysis is crucial to understand their behavior under seismic forces. The aqueducts, for which 3D finite element models were prepared, were analyzed using the combined Eulerian-Lagrangian approach, taking into account two main scenarios: with and without fluid-structure interaction. With this approach, it’s aimed to observe the damping effect of the fluid, the sloshing on its surface and possible water overflows. Using the CDP material model for masonry units, the analyzes were diversified with various ground motion records at different Peak Ground Acceleration values. The presence of water in the channel significantly reduced displacements, stresses and damage rates. In addition, it was seen that the spectral acceleration value of the earthquake record was more effective in causing damage rather than the PGA value. Graphics and contour diagrams are given to understand and compare the analysis results of aqueduct models under different conditions.

Experimental and numerical analysis of tunnel primary support using recycled, and hybrid fiber reinforced shotcrete

Shotcrete plays a pivotal role in the construction of tunnels and underground structures; however, its inherent brittleness necessitates reinforcement to enhance ductility. This research explores the use of fiber-reinforced shotcrete as primary tunnel support to enhance ductility and reduce brittleness. Traditional steel mesh reinforcement complexities have led to the investigation of alternative materials. The research evaluates different fiber mix designs, including industrial steel, recycled steel fibers from tires, and Forta fibers, examining their strength parameters and deformation performance. A 3D finite element model is used to simulate a horseshoe-shaped tunnel with optimal mix design and plain shotcrete in a soil environment. The study finds that hybrid industrial and recycled fibers are more effective than single fibers, enhancing compressive, tensile, and flexural strength and reducing ground surface settlement and tensile damage. The optimal mix design of this study has increased compressive, tensile, and flexural strength, as well as flexural toughness, compared to plain shotcrete. Numerical modeling reveals that utilizing fiber reinforced shotcrete made out of optimal fiber mix design as primary support results in a significant reduction in ground surface settlement and tensile damage value. Furthermore, the study shows a significant reduction in the damaged zone area under tensile stresses. The results of the study highlight the potential of fiber reinforced shotcrete as a primary support for tunnels, leading to improved performance and sustainability in tunnel construction.

Hybrid self-centering braces with NiTi-SMA U-shaped and frequency-dependent viscoelastic dampers for structural and nonstructural damage control

This study proposes a novel hybrid self-centering brace (HSB) with structural and nonstructural damage-control capabilities for the seismic resilience of building structures. NiTi-based shape memory alloy (SMA) U-shaped dampers and frequency-dependent viscoelastic dampers are included in HSBs to provide self-centering capacity for avoiding structural damage and velocity-proportional damping for nonstructural damage control, respectively. This study includes two major parts, i.e. component-level mechanical behaviour and system-level dynamic response evaluation of HSBs. The working principle and desired nonlinear behaviour of the HSB were first presented. Then, the preliminary design procedures were developed. After validating the efficiency of the modelling of SMA U-shaped and viscoelastic dampers using prior test results, comprehensive 3D finite element model simulations were performed to study the hysteretic responses of HSBs under cyclic loadings, where the influences of loading frequency, loading amplitude and viscoelastic damper’s contribution were considered. Analysis results demonstrated that the designed HSB can attain the desired deformation mechanism and nonlinear hysteretic behaviour during dynamic cyclic loadings. Higher loading frequency and amplitude and more contributions of the viscoelastic damper can introduce higher loading and energy-dissipation capacity. Subsequent dynamic analyses of the single-degree-of-freedom systems were conducted to assess the dynamic performance and advantages of HSBs. Results showed that compared with traditional self-centering systems without velocity-proportional damping, the systems with HSBs can achieve zero residual deformations, lower base shear demands and smaller absolute acceleration responses, confirming the promising potential of HSBs in developing seismic resilient building structures by reducing structural and nonstructural damage.

Accurate Time-Efficient Thermal Modeling and Discretization Methodology Applied to Electric Machines

The increasing power density of electric machines for electric vehicles necessitates accurate thermal analysis in design and real-time condition monitoring through simulations with computation times not exceeding one minute. Existing thermal models are either accurate and complex while requiring large computation times or they are fast but lack accuracy. To address these challenges, a novel hybrid modeling methodology is proposed, combining 3D lumped parameter thermal networks with thermal resistances extracted from 2D finite element models (FEM) to reduce computation demands while preserving accuracy. Further, a methodology is introduced to define the thermal network discretization for a desired threshold accuracy while minimizing computation time. Analysis of this discretization methodology demonstrates that the desired accuracy can be reliably achieved with limited computational effort. Verification of the modeling methodology against 3D FEM in OpenFOAM shows key component temperature deviations of less than 0.67 °C in steady-state simulations and 0.75 °C in transient simulations, comparable to the 3D FEM uncertainty of 0.39 °C. The computation times for steady-state and transient simulations are 18.3 s and 43.5 s respectively compared to 44 h and 175 h for 3D FEM. This highlights the time efficiency of the proposed methodology while maintaining accuracy for key component temperatures.

3D nonlinear finite element simulation of RC beams corroding under service loads

ABSTRACT In the present study, a 3D nonlinear finite element (FE) model has been developed to simulate reinforcement corrosion in RC beams under sustained load. The corrosion effects were modelled as bond deterioration at the rebar-concrete interface, damage in cover concrete and mechanical degradation of corroded rebars. A simplified temperature field-based methodology was proposed in the commercial program ABAQUS to simulate the evolution of corrosion degradation while a sustained load is maintained in the beam. The FE model was validated using available experimental results in terms of cracking behaviour and load-deflection response and found to be in good agreement. Parametric investigations were carried out to analyse the performance of beams under simultaneous corrosion and sustained loads. The findings demonstrated that pitting corrosion under simultaneous loading conditions can significantly compromise a beam’s structural integrity. Specifically, a scenario with 20% mass loss and 60% sustained load resulted in a 33.3% increase in concrete strain and a significant increase of 738% in the tensile steel strain leading to yielding. Furthermore, the influence of load and corrosion sequence on the behaviour of the beam was assessed by comparing its response when corroded under sustained load with a beam corroded and subsequently loaded.

A comparative evaluation of masticatory load distribution in different types of prosthesis with varying number of implants: A FEM analysis

AimTo identify the optimal number and position of implants to reduce stress concentration on the implant, denture, and attachment system for sustaining an overdenture prosthesis. Materials and methodsBy incorporating one to eight indigenous implants with bar-type attachments, eight 3D finite element models of mandibular overdentures were created. All models received a 200 N vertical load, and the biomechanical characteristics of peri-implant bone were assessed. ResultThe study observed that with a vertical load of 200 N, the maximum equivalent stress around peri-implant tissue in all models was within the physiological tolerance threshold of bone. The von Mises stress values ranged from 116.18 MPa to 536.7 MPa. ConclusionThe three-implant-supported overdenture model revealed superior peri-implant stress, stability, cost-effectiveness, and hygiene maintenance outcomes. Placing a third implant in the mid-symphysis region may offer a practical solution to reduce rotations in two-implant-supported overdentures.

Stress distribution on teeth restored with veneers using various incisal preparation designs: A 3D finite element analysis study

Statement of problemThe type of veneer preparation is often chosen according to the patient’s tooth structure and occlusion. Taking biomechanics into account in this decision-making process provides the clinician with more technical information on how to improve the clinical longevity of restorations. However, biomechanical analyses of veneer preparation designs are sparse. PurposeThe purpose of this 3-dimensional (3D) finite element analysis with microcomputed tomography (µCT) data obtained from realistic models was to assess the influence of different preparations for ceramic and composite resin veneers on restoration and resin layer stress distribution. Material and methodsFour replicas of a central incisor were printed and prepared for laminate veneers with 4 different incisal edge designs: shoulder (SH), palatal chamfer (PC), palatal chamfer and oblique fracture involving the distal angle (OF-PC), and palatal chamfer involving horizontal incisal fracture (IF-PC). After fabrication and cementation of the veneers, the restored replicas were assessed with µCT, and 3D finite element models were built. A 100-N load was applied on the palatal surface at 60 and 125 degrees relative to the longitudinal axis. Maximum principal stress and stress distribution on the veneers, cement layer, and tooth structure were calculated and analyzed. ResultsThe SH preparation exhibited better stress distribution than the PC preparation, and the cement layer and the veneer were subjected to lower stress. The IF-PC preparation had better stress distribution than the OF-PC. The shoulder and IF-PC showed higher stress on laminate veneers, but lower stress on the cement layer. Ceramic veneers exhibited lower stress than composite resin veneers. ConclusionsThe different incisal preparations for laminate veneers influenced stress distribution on restorations and on the resin cement layer. The shoulder type preparation showed better stress distribution and the composite resin veneers showed unfavorable results compared with the ceramic veneers.

Seismic response of a ground-bridge structure system in horizontal and inclined liquefiable site to near-fault pulse-like ground motions

Previous earthquake events have indicated that bridge structures are more vulnerable to severe damage when subjected to near-fault pulse-like (NF-P) ground motions. Structural seismic damage will be exacerbated when bridges are located in liquefiable site. Meanwhile, the seismic response of bridges located in inclined liquefiable site is different from that of bridges in horizontal liquefiable site. This study focuses on exploring the seismic response of a ground-bridge structure system located in horizontal and inclined liquefiable site to the NF-P ground motions. The seismic response of the whole system to the near-fault non-pulse-like (NF-NP) and far-field (FF) ground motions is studied for comparison. Two three-dimensional (3D) finite element (FE) models of the ground-structure system considering soil-pile interaction are developed. They are associated with the horizontal and inclined liquefiable site, respectively. The seismic responses of the ground-bridge structure system to the three different types of ground motions (i.e., NF-P, NF-NP and FF) are comprehensively evaluated from two perspectives. On the one hand, the seismic time history responses of the ground-bridge structure system under single representative ground motion are comprehensively assessed. On the other hand, the average responses of the ground-bridge structure system under multiple types of ground motions are explored. The results reveal that in contrast to the NF-NP and FF ground motions, the NF-P ground motions have more significant effect on various responses of the ground-bridge structure system. In addition, compared to the bridge in horizontal liquefiable site, the bridge in inclined liquefiable site exhibits larger seismic responses. Therefore, the velocity pulse effect of the NF-P ground motions and the liquefaction effect of the inclined site need to be emphasized in structural seismic design.

3D non-linear finite element modeling of one-way RC slab strengthened with concrete overlay

Reinforced concrete (RC) slabs often require strengthening due to deterioration, increased loads, or design changes. Concrete overlay is a commonly used technique for strengthening, but predicting the behavior of overlays can be challenging. To address this challenge, this study presents a three-dimensional non-linear finite element model using Abaqus/CAE to analyze the structural performance of one-way reinforced concrete (RC) slabs strengthened with concrete overlays. The model considers three different bonding conditions at the interface: friction, epoxy, and shear connectors. The predictions of the model were compared to experimental results, and it was found that the model accurately captures the load–deflection response, cracking behavior, and interfacial slip in the strengthened slabs. Additionally, the model accurately predicts the new structural capacity of the strengthened slabs, showing good agreement with experimental data. These findings offer valuable insights into the effectiveness of concrete overlays as a traditional strengthening technique for RC slabs. Moreover, the developed finite element model provides a useful tool for optimizing the design of these overlays.

BEHAVIOR OF A STEEL STRUCTURE RAILWAY BRIDGE UNDER DYNAMIC LOADINGS

Most old steel structure railway bridges in Indonesia have deteriorated throughout their service life since they were constructed almost a century ago. However, those bridges' performance must be maintained to have essential safety issues and live extension of the railway bridge structure. Therefore, inspecting and evaluating those steel railway bridges is necessary to maintain the service requirement. Vertical deformation of the steel railway bridge caused by dynamic loadings needs to be observed. The objective of the study is to assess the old steel railway bridge by evaluating the strength characteristics of the structures against the working forces, particularly the moving, wind, and seismic loads. In order to understand the phenomena impacted by the dynamic loadings, the steel structure railway bridge was instrumented using deformation sensors, strain gages, accelerometers, and passive infrared. The steel railway bridge was analyzed using a 3D finite element model. This study discussed the influence of dynamic loadings on the steel structure railway bridge. This paper elaborates and provides suggestions to solve problems and recommended action in practice for future study. This paper may be useful for researchers and practicing engineers.

Modelling and numerical analysis for rotatory friction welding of U75V steel rails

This article presents a rotatory friction welding (RFW) method for U75V steel rail, aiming to mitigate challenges related to property discrepancies between the as-welded joints and the rail base metal (BM) and to narrow the heat-affected zone (HAZ) in conventional flash-butt welding (FBW) joints. A rotational intermediate plate is designed for rails with non-axisymmetric cross-sections, necessitating stationary during RFW. Advantages include achieving a relatively uniform welding heat input and maintaining the peak temperature of the contact interface near A1. To implement these concepts, a 2D finite element (FE) model for the RFW process of U75V rail steel rods was established and validated through experiments with identical process parameters. Microstructure predictions derived from continuous cooling transformation diagram confirm that ferrite microstructure is formed near A1 through rail steel RFW. Subsequently, a 3D FE model for intermediate plate RFW steel rails is developed to explore appropriate process parameter combinations. A suitable process parameters combination was identified, ensuring the peak temperature of the majority model contact interface does not exceed A1, resulting in a 76.7% reduction in HAZ (from ∼50 to 11.66 mm), and axial shortening of 8.10 mm, a significant decrease compared to the usual burn-off (30–40 mm) during FBW. These findings underscore the efficacy of this innovative welding solution and emphasize the significance of simulation technology in process optimization.

Calculation of Equivalent Thermal Conductivity of Motor Winding Based on Conductor Distribution

Abstract To improve the calculation accuracy of motor thermal field, this paper proposes a method for calculating the equivalent thermal conductivity of winding based on conductor distribution, and a 24-slot 26-pole outer-rotor permanent magnet synchronous motor is taken as an example to verify the method. First, the equivalent thermal conductivity of winding is calculated by traditional equivalent methods, namely the equivalent insulation layer method and the area weighting method. Second, the equivalent thermal conductivity of the slot winding and end winding in the three-dimensional (3D) direction is derived using the conductor distribution method. Third, the 3D finite element models corresponding to three equivalent methods are established, and the temperature of the motor under different operating conditions is simulated and compared. Finally, an experimental platform is built to verify the accuracy of the proposed conductor distribution method.

Interfacial Shear Behavior of Composite Concrete Substrate to High-Performance Concrete Overly After Exposure to Elevated Temperature

Basically, the interface shear strength between two concrete layers of varying ages must be sufficient to withstand the applied actions on the structure, specifically fire attack, which may cause the complete collapse of the composite structure. Thus, interfacial shear behavior was investigated and analyzed in this paper under the influence of a set of parameters, including temperature (25, 200, 400, and 600 °C), time exposure (30, 60, 90, 120, and 180 min), concrete type, and fibers type (polypropylene fiber (PPF), steel fiber (SF), and hybrid fiber) by employing a Z-shape push-off test. The test consists of two parts with different ages: normal strength concrete (NCS) and high-performance concrete (HPC). HPC includes high-strength concrete (HSC) and fly ash concrete (FAC). Initially, twenty-five Z-shaped push-off tests were made, four of which were cast as one unit (NSC/or concrete with hybrid (FSP)), and the rest were composite specimens. Furthermore, a 3D finite element model of a composite push-off specimen was developed to simulate and analyze the impact of various time and temperature exposures on the interfacial shear strength of composite specimen N-FSP. The results indicated that temperature degree and exposure time adversely affected the interfacial shear strength. Also, interfacial shear strength is significantly influenced by fiber types. Including combined fiber (SF + PPF) improved the interfacial shear strength by 114% compared to the composite specimen NSC-NSC after exposure to a temperature of 600 °C. In contrast, using PPF negatively affected the interfacial shear strength, recording only 84% of the composite specimen NSC-NSC. In addition, the inclusion of supplementary cementitious material enhanced the interfacial shear strength by 60.5% in the NSC-FAC composite specimen with 30% FA, compared to the NSC-NSC specimen. Finally, a finite element (FE) model was proposed with a satisfactory level of accuracy (0.95 to 1.03) in predicting the maximum shear strength. Additionally, the difference between the FE and experimental stiffness was between 0.92 and 1.07.

Free vibration analysis of composite beams and laminated reinforced panels by refined dynamic stiffness method and CUF-based component-wise theory

This paper presents a comprehensive analysis of the free vibration behavior of composite beams and laminated reinforced panels. Employing high-order theories with displacement variables only, the investigation combines dynamic stiffness method (DSM) and the Carrera unified formulation (CUF). The 3D displacement field can be expanded in the framework of CUF as any order of generic unknown variables over the cross-section, in the case of beam theories. Specifically, Lagrange expansions (LE) of cross-sectional displacement variables are considered, enabling the refinement modeling of complex cross-sections with different layers (layer-wise) and components (component-wise). The governing differential equations and natural boundary conditions are derived using the principle of virtual displacements. Subsequently, an exact dynamic stiffness matrix is developed by establishing a relationship between the amplitudes of harmonically varying loads and the corresponding responses. The Wittrick–Williams algorithm is employed to solve the transcendental eigenvalue problem resulting from this approach. The Lagrange-based CUF(LE)-DSM outperforms other polynomials (such as Taylor)-based CUF in analyzing composite structures, enabling detailed analysis with various geometries, lamination schemes and boundary conditions. It delivers accurate results in reinforced panel analysis, requiring only 13% of the DOFs compared to 3D FE models, thereby reducing computational costs significantly, as confirmed by comparative studies and validations.

Effect of interference magnitude on fretting wear and fatigue strength of scaled press-fitted railway axles

Railway axles represent a typical example of interference fit structures that are vulnerable to fretting fatigue induced by rotational bending loads. This study investigated the influence of interference magnitude on the fretting wear and fatigue of scaled press-fitted railway axles. In this study, the fretting fatigue strength was tested at three interference magnitudes within the standard range, and the fretting wear and fretting cracks were analyzed. The results indicate a trend where the fretting fatigue strength initially rises and then declines with increasing interference magnitude. 3D finite element models, with and without considering the fretting wear, were developed to investigate the variation in contact parameters in the press-fitted region, and the Smith-Watson-Topper multiaxial fatigue model was employed to evaluate the fretting fatigue behavior. The results demonstrate that the competition between contact pressure and axial slip produces relatively slight fretting wear when the interference magnitude takes the middle value, thereby causing a minor stress concentration and ultimately leading to higher fretting fatigue strength.