Wear performance prediction method and system for hydrogen-barrier coatings in rubber dynamic seals for hydrogen energy equipment

Abstract

The present invention relates to the prediction of the wear performance of hydrogen barrier coatings on the surface of materials, and aims to provide a method and system for predicting the wear performance of hydrogen barrier coatings on rubber dynamic seals for hydrogen energy equipment. The method comprises: establishing a finite element model of friction and wear of a rubber material substrate, a hydrogen barrier coating on the surface of the substrate, and a wear part above the coating; establishing a bilinear cohesive force unit for the superimposed friction effect between the substrate and the hydrogen barrier coating on the surface of the substrate; establishing coupling constraints for realizing the coordinated rotational motion of the substrate and the hydrogen barrier coating on the surface of the substrate; based on the ABAQUS-UMAT user subroutine module, using the UMESHMOTION subroutine combined with the ALE grid adaptive technology to calculate the wear of the hydrogen barrier coating on the surface of the substrate to obtain the wear depth; obtaining the solved coating wear depth, surface stress and cross-sectional stress to predict the wear performance of the coating. The present invention can accurately obtain the wear depth, bonding surface and surface stress response results of the hydrogen barrier coating on the surface of the rubber sealing material during the rotational wear process, and predict the wear performance of the coating.

Description

Technical Field

[0001] The present invention relates to the field of wear performance prediction of material surface coatings, and in particular to a method and system for predicting wear performance of hydrogen-barrier coatings of rubber dynamic seals for hydrogen energy equipment. Background Art

[0002] As the most promising sustainable clean energy in the 21st century, hydrogen energy can help achieve deep decarbonization in multiple fields such as energy, transportation, and petrochemicals. Accelerating the pace of hydrogen energy development has become a global consensus. However, the hydrogen damage problem faced by hydrogen-prone materials has become a key threat to the safe operation of hydrogen energy equipment. By preparing a hydrogen-resistant coating on the surface of the material, the penetration of hydrogen can be slowed down without destroying the mechanical properties of the material itself, thereby achieving protection against hydrogen damage to the material. However, polymer materials such as rubber used in dynamic sealing components of hydrogen pressure-bearing equipment not only face the problem of hydrogen damage caused by hydrogen intrusion, but also face failure problems caused by friction and wear. In actual working conditions, the hydrogen-resistant coating will have rotating, reciprocating and other frictional contacts with the wear parts, which will not only cause continuous wear of the coating, but may even cause peeling between the coating and the substrate under the action of shear force. Therefore, in order to ensure the long-term and reliable operation of the material, the wear performance of the hydrogen-resistant coating on its surface must be predicted.

[0003] Although some studies have proposed finite element simulation schemes for the wear of materials or their surface coatings [Wang He, Jin Yifan, Wang Yanxiang, et al. Simulation study of wear mechanism of AIN insulation coating based on ABAQUS [J]. Mechanical and Electrical Product Development and Innovation, 2023, 36(4): 34-38.], [Wang Kaimo, Shen Huoming, Wang Yuxing, et al. Research on numerical prediction model of tangential micro-wear of coated zirconium alloy cladding tube [J]. Journal of Chongqing University of Technology, 2024, 38(3): 117-122.], [A wear simulation method based on ABAQUS 202111291454.5], their evaluation and prediction methods for the wear performance of materials and components are all based on the wear simulation method of the UMESHMOTION subroutine in ABAQUS. This method is relatively effective in calculating and evaluating the wear and performance of a single structural component or material. However, when faced with the wear simulation of a multi-layer material combination model, it can only obtain the wear process of the top layer of material, and cannot accurately reflect the stress response and bonding state of the material interlayer interface during the friction and wear process. Therefore, it is impossible to accurately predict the interlayer peeling behavior of the hydrogen barrier coating on the material surface caused by stress concentration and excessive shear force during the wear stage. In addition, since the wear performance test of the material usually adopts the working condition of rotational wear, this working condition has serious nonlinear boundary constraint problems compared to the reciprocating wear motion working condition carried out in the existing literature, especially the coordinated rotation constraint conditions of the substrate component and the coating component, which are more difficult to set.

[0004] Therefore, how to calculate the coating wear under highly nonlinear boundary-constrained rotating wear conditions while introducing the stress response calculation of the interface between the coating and the substrate is an urgent problem to be solved in the process of predicting the wear performance of hydrogen barrier coatings on the surface of rubber sealing materials. Summary of the Invention

[0005] In order to solve at least one of the problems existing in the prior art, the present invention provides a high-precision method for predicting the wear performance of hydrogen-resistant coatings on rubber dynamic seals for hydrogen energy equipment, which simultaneously considers the calculation of coating wear depth and stress response.

[0006] To achieve the purpose of the present invention, the present invention provides a method for predicting the wear performance of a hydrogen barrier coating on a rubber dynamic seal for hydrogen energy equipment, comprising:

[0007] Establish a finite element model of friction and wear of a rubber-containing substrate, a hydrogen barrier coating on the substrate surface, and a wear part above the coating;

[0008] Establish a bilinear cohesive force unit with superimposed friction effect between substrate and substrate surface hydrogen barrier coating;

[0009] Establish the coupling constraint conditions to achieve the coordinated rotational motion of the substrate and the hydrogen barrier coating on the substrate surface;

[0010] Based on the ABAQUS-UMAT user subroutine module, the UMESHMOTION subroutine combined with the ALE grid adaptive technology is used to calculate the wear of the hydrogen barrier coating on the substrate surface and obtain the wear depth.

[0011] The solved coating wear depth, surface stress and cross-sectional stress are obtained to predict the coating wear performance.

[0012] Furthermore, the steps of establishing the friction and wear finite element model are as follows: based on ABAQUS, a rubber material substrate, a hydrogen barrier coating on the substrate surface, and a wear member above the coating are established, material properties and mesh division are set respectively, and analysis steps and general contact properties are set after assembly.

[0013] Furthermore, the method of establishing the bilinear cohesive force unit with superimposed friction effect between the substrate and the hydrogen barrier coating on the substrate surface is to insert a bilinear cohesive force zero-thickness viscous interface unit with superimposed friction effect between the substrate and the hydrogen barrier coating on the substrate surface in the friction and wear finite element model to describe the interlayer damage between the substrate and the substrate surface coating. The bonding strength constitutive law of the interface unit is as follows:

[0014] The final effective relative displacement δ of the mixed damage mode under power law distribution is calculated as follows:

[0015]

[0016] Where β is the mixed damage mode parameter, δ0 is the damage relative displacement, K is the interface stiffness, G nc and G sc represent the critical values ​​of damage mode I and damage mode II, respectively, and k is a constant;

[0017] When friction effects are taken into account, the normal and tangential components of the cohesive force are given by:

[0018]

[0019] in is the average normal relative displacement, d is the damage variable, is the average tangential relative displacement vector, T t,f It represents the shear stress acting on the damaged part of the interface due to the friction effect.

[0020] Furthermore, the shear stress T acting on the damaged part of the interface due to the friction effect is t,f The expression is

[0021] T t,f =dt t,f (11)

[0022] in

[0023]

[0024] Where μ is the friction coefficient, ζ is a constant, s represents the sliding rate, x is the maximum effective separation in the shear mode under the effect of friction, g(x) is the damage variable under the effect of friction, and t t,f is the interface friction traction.

[0025] Furthermore, the coupling constraint conditions for achieving coordinated rotational motion between the substrate and the hydrogen barrier coating on the substrate surface are as follows: first, a reference point is set on a preset range area at the bottom of the substrate, and the axis perpendicular to the substrate plane is taken as the z direction. The preset range area is coupled and constrained with the reference point, and the angular velocity of rotation along the z direction is set.

[0026] Furthermore, a reference point is set at the center of the preset range area.

[0027] Furthermore, the preset range area is the area of ​​10% to 15% of the bottom thickness of the substrate.

[0028] Furthermore, the wear calculation of the hydrogen barrier coating is performed based on the ABAQUS-UMAT user subroutine module using the UMESHMOTION subroutine combined with the ALE grid adaptive technology as follows:

[0029] The subroutine uses the Archard wear model:

[0030]

[0031] Where V represents the wear volume of the material, σ is the wear coefficient, a is the relative sliding distance between the contact units, F is the normal load at the wear contact point, and H represents the hardness of the material. Dividing both sides of Equation (5) by the wear contact area A yields:

[0032]

[0033] Where h represents the wear depth, P represents the contact pressure, and λ is the specific wear rate;

[0034] In order to obtain the wear depth from the contact pressure at the wear contact node, it is necessary to redefine the Archard formula at the local scale. In addition, since wear is a dynamic process, the time-varying differential Archard formula can be obtained:

[0035]

[0036] Where dt is the infinitesimal time for the wear part to move, v is the sliding velocity, and Equation (7) is integrated over the sliding time using Euler integration to obtain the wear depth at each wear contact node for a given wear cycle:

[0037] h i+1 =h i +λvP i Δt (16)

[0038] where h i and h i+1 is the wear depth of a given contact node at the i-th and i+1-th wear cycle increments, P i is the contact pressure at the contact node, Δt is the calculation time step;

[0039] After calculating the wear depth of a single wear node at each cycle increment, the ALE mesh adaptation technology is used to reconstruct the wear node and the node positions below it to obtain a smooth wear profile.

[0040] Furthermore, the coating wear depth, surface stress and cross-sectional stress obtained by the solution are obtained to predict the coating wear performance. Specifically, in the process in which the rubber material substrate and the hydrogen barrier coating on the substrate surface come into contact with the upper wear member and produce relative displacement, the coating wear depth, surface stress and cross-sectional stress data are obtained by combining the constitutive law of the cohesive force unit and the constitutive law of wear to simultaneously calculate and obtain the coating wear performance prediction.

[0041] The present invention also provides a wear performance prediction system for hydrogen-resistant coatings of rubber dynamic seals for hydrogen energy equipment, which is used to implement the above method. The system includes the following modules:

[0042] Friction and wear model building module, used to establish the friction and wear finite element model of the rubber material substrate, the hydrogen barrier coating on the substrate surface, and the wear parts above the coating;

[0043] Cohesion unit establishment module, used to establish a bilinear cohesion unit with superimposed friction effect between substrate and substrate surface hydrogen barrier coating;

[0044] A constraint condition establishment module is used to establish coupling constraint conditions for realizing the coordinated rotational motion of the substrate and the hydrogen barrier coating on the substrate surface;

[0045] Wear calculation module, which is used to calculate the wear of hydrogen barrier coating on substrate surface based on ABAQUS-UMAT user subroutine module, using UMESHMOTION subroutine combined with ALE grid adaptive technology to obtain wear depth;

[0046] The prediction module is used to obtain the solved coating wear depth, surface stress and cross-sectional stress to predict the coating wear performance.

[0047] Compared with the prior art, the present invention can at least achieve the following beneficial effects:

[0048] First, the present invention establishes a bilinear cohesive force unit with superimposed friction effect between substrate and hydrogen barrier coating. By defining the constitutive law of the bonding strength of the interface unit, it can accurately predict the interlayer damage between substrate and coating and the change of bonding strength of the interface unit.

[0049] Second, by coupling a preset area below the substrate with a regional reference point and defining rotational constraints, the present invention achieves a coordinated rotational motion state between the substrate and the coating without affecting the stress response accuracy of the coating-substrate interface above, accurately simulating rotational wear conditions.

[0050] 3. The present invention is based on the use of the UMESHMOTION subroutine combined with the ALE grid adaptation technology to accurately predict the wear depth of the coating surface during the wear stage. By combining the surface stress and cross-sectional stress response results of the coating, the wear behavior of the coating can be accurately reflected and predicted. BRIEF DESCRIPTION OF THE DRAWINGS

[0051] Figure 1 A finite element model diagram of the friction and wear of a rubber material substrate, a hydrogen barrier coating on the substrate surface, and a wear part above the coating provided in an embodiment of the present invention.

[0052] Figure 2 This is a flow chart of numerical simulation of the proposed bilinear cohesive force unit with superimposed friction effect and the use of the UMESHMOTION subroutine combined with the ALE grid adaptation technology for calculation in an embodiment of the present invention.

[0053] Figure 3 1 is a comparison chart of the simulation results and experimental results of the coating wear depth in an embodiment of the present invention.

[0054] Figure 4 This is a diagram showing the Mises stress results on the coating surface in an embodiment of the present invention.

[0055] Figure 5 This is a diagram showing the Mises stress results of the coating cross section in an embodiment of the present invention. DETAILED DESCRIPTION

[0056] The present invention will be described in further detail below with reference to the embodiments and drawings, but the embodiments of the present invention are not limited thereto.

[0057] See also Figure 2 The present invention provides a method for predicting the wear performance of a hydrogen barrier coating of a rubber dynamic seal for hydrogen energy equipment, comprising the following steps:

[0058] Step 1: Establish a finite element model of friction and wear of a rubber-containing substrate, a hydrogen barrier coating on the substrate surface, and a wear member located above the coating.

[0059] In this step, when establishing the friction and wear finite element model, the rubber material substrate, the hydrogen barrier coating on the substrate surface, and the wear part above the coating are established based on ABAQUS, and the material properties and mesh division are set respectively. Then, after assembling them (the rubber material substrate, the hydrogen barrier coating on the substrate surface, and the wear part), the analysis step and general contact properties are set.

[0060] Step 2: Establish a bilinear cohesive force unit with superimposed friction effect between the substrate and the hydrogen barrier coating on the substrate surface.

[0061] The bilinear cohesive force unit for establishing the superimposed friction effect between the substrate and the hydrogen barrier coating is specifically:

[0062] A bilinear cohesive zero-thickness viscous interface element with superimposed friction effect is inserted between the substrate and the hydrogen barrier coating on the substrate surface in the friction and wear finite element model to describe the interlayer damage between the substrate and the substrate surface coating. The constitutive law of the bonding strength of the interface element is as follows:

[0063] The final effective relative displacement δ of the mixed damage mode under power law distribution is calculated as follows:

[0064]

[0065] Where β is the mixed damage mode parameter, δ0 is the damage relative displacement, K is the interface stiffness, G nc and G sc represent the critical values ​​of damage mode I and damage mode II respectively, and κ is a constant.

[0066] When the friction effect is considered, the normal component of the cohesive force T coh,n and the tangential component T coh,t It is given by:

[0067]

[0068] in is the average normal relative displacement, d is the damage variable, is the average tangential relative displacement vector, T t,f It represents the shear stress acting on the damaged part of the interface due to the friction effect.

[0069] T t,f =dt t,f (19)

[0070] in

[0071]

[0072] Where μ is the friction coefficient, ζ is a constant. In some embodiments of the present invention, ζ is 0.03, s represents the slip rate, x is the maximum effective separation in the shear mode under the effect of friction, g(x) is the damage variable under the effect of friction, and t t,f is the interface friction traction.

[0073] Step 3: Establish coupling constraint conditions to achieve coordinated rotational motion between the substrate and the hydrogen barrier coating on the substrate surface.

[0074] In this step, the coupling constraint conditions for achieving coordinated rotational motion between the substrate and the hydrogen barrier coating are specifically as follows: first, a reference point is set at the center position of a preset range area at the bottom of the substrate, and then the area is coupled with the reference point, and the angular velocity of rotation along the z direction (with the axis perpendicular to the substrate plane as the z direction) is set.

[0075] In some embodiments of the present invention, the preset range area is the area within the lower 10% to 15% of the substrate thickness. Preferably, the preset range area is the area within the lower 10% of the substrate thickness.

[0076] Step 4: Based on the ABAQUS-UMAT user subroutine module, the UMESHMOTION subroutine combined with the ALE grid adaptive technology is used to calculate the wear of the hydrogen barrier coating.

[0077] The wear calculation of the hydrogen barrier coating on the substrate surface is performed based on the ABAQUS-UMAT user subroutine module using the UMESHMOTION subroutine combined with the ALE grid adaptive technology as follows:

[0078] The subroutine uses the Archard wear model:

[0079]

[0080] Where V represents the wear volume of the material, σ is the wear coefficient, which reflects the wear resistance of the contact unit, a is the relative sliding distance between the contact units, F is the normal load at the wear contact point, and H represents the hardness of the material. Dividing both sides of Equation (5) by the wear contact area A yields:

[0081]

[0082] in Indicates the wear depth, represents the contact pressure, and λ = K / H is defined as the specific wear rate.

[0083] In order to obtain the wear depth from the contact pressure at the wear contact node, it is necessary to define the differential Archard formula at the local scale. In addition, since wear is a dynamic process, the differential Archard formula with time t can be obtained:

[0084]

[0085] Where dt is the infinitesimal time for the wear part to move, v is the sliding velocity, and Equation (7) is integrated over the sliding time using Euler integration to obtain the wear depth at each wear contact node for a given wear cycle:

[0086] h i+1 =h i +λvP i Δt (24)

[0087] where h i and h i+1 is the wear depth of a given contact node at the i-th and i+1-th wear cycle increments, P i is the contact pressure at the contact node, Δt is the calculation time step;

[0088] After calculating the wear depth of a single wear node at each wear cycle increment, the ALE mesh adaptation technology is used to reconstruct the positions of the wear node and the nodes below it to obtain a smooth wear profile.

[0089] Step 5: Obtain the solved coating wear depth, surface stress, and cross-sectional stress to predict the coating wear performance.

[0090] In this step, the coating wear depth, surface stress and cross-sectional stress obtained by the solution are obtained to predict the coating wear performance as follows:

[0091] During the process of contact and relative displacement between the rubber material substrate and the hydrogen-barrier coating on the substrate surface and the upper wear member, the constitutive equation of the bilinear cohesive force unit and the constitutive equation of wear are combined to synchronously calculate the coating wear depth, surface stress and cross-sectional stress data, thereby predicting the coating wear performance.

[0092] The wear performance of a material mainly refers to its wear resistance, which is mainly judged by the wear characteristics (such as wear depth, wear volume, etc.) obtained through tribological tests or numerical simulations under specific working conditions. This parameter determines the degree of material wear during the wear process; by extracting and analyzing the stress response of the material during the wear process, it is believed that obvious stress concentration will lead to local failure of the material, corresponding to the deterioration of wear performance.

[0093] In some embodiments of the present invention, the effectiveness of the present invention is verified by using specific examples:

[0094] A friction and wear finite element model including the grinding ball 1 (i.e., the grinding member), the hydrogen barrier coating 2 on the substrate surface, and the rubber substrate 3 was established in ABAQUS 6.14, as shown in the following figure: Figure 1 As shown. The mechanical properties of the substrate and the hydrogen barrier coating on the substrate surface, the friction coefficient and wear rate of the coating surface, and the bonding performance parameters between the coating and the substrate are all obtained through experiments. The hydrogen barrier coating 2 on the substrate surface and the substrate 3 are both meshed with full hexahedral meshes, and the mesh properties are defined using the C3D8R element type. The grinding ball 1 is defined as an analytical rigid body. The Surface to Surface contact type is selected, and the Coulomb friction model and the "Penalty" tangential algorithm are applied, where the grinding ball 1 is regarded as the Master face and the coating specimen is regarded as the Slave face. The normal behavior is set to "hard contact" with pressure interference, and standard contact control and automatic stabilization functions are used to ensure the stability of the contact.

[0095] ABAQUS / Standard is used to calculate the surface wear depth and stress response of the coating during the rotational wear process, such as Figure 2As shown in the figure, after completing the three-dimensional wear model, material property definition and load boundary constraint conditions, the simulation calculation is immediately carried out. When the upper grinding steel ball comes into contact with the coating, the contact stress response begins, and then the incremental step is applied. At the beginning of each incremental step iteration, the user material subroutine UMAT is called to simultaneously perform interface wear calculation and bonding strength calculation. In the wear calculation, the main program first inputs the current contact point information (contact stress, relative displacement, node position), and then calculates the wear depth of the node based on the contact point information, and further updates the node geometric features, and then returns the updated state to the subroutine. In the interface bonding strength calculation, the normal and tangential node displacement jumps are first calculated based on the current node position change information, and then the normal and tangential stresses at the node are calculated. Then, based on the bonding strength constitutive law of the interface unit, it is judged whether the cohesive element unit has failed, so as to determine whether to further calculate the damage amount of superimposed friction and slip. Finally, it is determined whether the iteration has converged. If the iteration does not converge, it returns to iterate again. If the iteration converges, it solves the next incremental step until the incremental step ends, thereby obtaining the wear results and stress response results of the hydrogen barrier coating on the given substrate surface.

[0096] Depend on Figure 3 It can be seen that the numerical simulation results of the surface wear depth of the hydrogen barrier coating of this embodiment under the rotating wear condition are consistent with the experimental results.

[0097] The surface and cross-sectional Mises stress results of the coating specimens at the initial stage of wear, during wear, and after wear are shown as follows: Figure 4 and Figure 5 As shown in the figure, by setting up a bilinear cohesive force unit with superimposed friction effect, the stress response of the hydrogen barrier coating on the substrate surface and the coating-substrate interface can be accurately obtained.

[0098] In some embodiments of the present invention, a system for predicting the wear performance of hydrogen-barrier coatings on rubber dynamic seals for hydrogen energy equipment is provided, which is used to implement the methods provided in the aforementioned embodiments. The system includes the following modules:

[0099] Friction and wear model building module, used to establish the friction and wear finite element model of the rubber material substrate, the hydrogen barrier coating on the substrate surface, and the wear parts above the coating;

[0100] Cohesion unit establishment module, used to establish a bilinear cohesion unit with superimposed friction effect between substrate and substrate surface hydrogen barrier coating;

[0101] A constraint condition establishment module is used to establish coupling constraint conditions for realizing the coordinated rotational motion of the substrate and the hydrogen barrier coating on the substrate surface;

[0102] Wear calculation module, which is used to calculate the wear of hydrogen barrier coating on substrate surface based on ABAQUS-UMAT user subroutine module, using UMESHMOTION subroutine combined with ALE grid adaptive technology to obtain wear depth;

[0103] The prediction module is used to obtain the solved coating wear depth, surface stress and cross-sectional stress to predict the coating wear performance.

[0104] As for the wear performance prediction system of the hydrogen-barrier coating of the rubber dynamic seal for hydrogen energy equipment disclosed in the embodiment, since it corresponds to the method disclosed in the embodiment, the description is relatively simple, and the relevant parts can be referred to the description of the method part.

[0105] The above examples are merely specific implementation cases of the present invention and are not intended to limit the present invention in any form. Although the present invention has been disclosed above with preferred implementation cases, they are not intended to limit the present invention. Any person skilled in the art can, without departing from the scope of the technical solution of the present invention, make certain changes or modifications to the above-disclosed structures and technical contents to create equivalent implementation cases with equivalent changes. However, any simple modifications, equivalent changes, and modifications made to the above implementation cases based on the technical essence of the present invention without departing from the content of the technical solution of the present invention are still within the scope of the technical solution of the present invention.

Claims

1. A method for predicting the wear performance of hydrogen barrier coatings on rubber dynamic seals for hydrogen energy equipment, characterized in that: The following steps are involved: Establish a finite element model of friction and wear of a rubber-containing substrate, a hydrogen barrier coating on the substrate surface, and a wear part above the coating; Establish a bilinear cohesive force unit with superimposed friction effect between substrate and substrate surface hydrogen barrier coating; Establish the coupling constraint conditions to achieve the coordinated rotational motion of the substrate and the hydrogen barrier coating on the substrate surface; Based on the ABAQUS-UMAT user subroutine module, the UMESHMOTION subroutine combined with the ALE grid adaptive technology is used to calculate the wear of the hydrogen barrier coating on the substrate surface and obtain the wear depth. Obtain the solved coating wear depth, surface stress and cross-sectional stress to predict the coating wear performance; The method for establishing the bilinear cohesive force unit with superimposed friction effect between the substrate and the hydrogen barrier coating on the substrate surface is to insert a bilinear cohesive force zero-thickness viscous interface unit with superimposed friction effect between the substrate and the hydrogen barrier coating on the substrate surface in the friction and wear finite element model to describe the interlayer damage between the substrate and the substrate surface coating. The bonding strength constitutive law of the interface unit is as follows: The final effective relative displacement δ of the mixed damage mode under power law distribution is calculated as follows: Where β is the mixed damage mode parameter, δ0 is the damage relative displacement, K is the interface stiffness, G nc and G sc represent the critical values ​​of damage mode I and damage mode II, respectively, and k is a constant; When friction effects are taken into account, the normal and tangential components of the cohesive force are given by: in is the average normal relative displacement, d is the damage variable, is the average tangential relative displacement vector, T t,f It represents the shear stress acting on the damaged part of the interface due to the friction effect.

2. The method for predicting the wear performance of hydrogen barrier coating of rubber dynamic seal for hydrogen energy equipment according to claim 1, characterized in that: The steps of establishing the friction and wear finite element model are as follows: based on ABAQUS, a rubber material substrate, a hydrogen barrier coating on the substrate surface, and a wear member above the coating are established, material properties and mesh division are set respectively, and analysis steps and general contact properties are set after assembly.

3. The method for predicting the wear performance of hydrogen barrier coating of rubber dynamic seal for hydrogen energy equipment according to claim 1, characterized in that: The shear stress T acting on the damaged part of the interface due to friction effect t,f The expression is T t,f =dt t,f (3) in Where μ is the friction coefficient, ζ is a constant, s represents the sliding rate, x is the maximum effective separation in the shear mode under the effect of friction, g(x) is the damage variable under the effect of friction, and t t,f is the interface friction traction.

4. The method for predicting the wear performance of hydrogen barrier coatings on rubber dynamic seals for hydrogen energy equipment according to any one of claims 1 to 3, characterized in that: The coupling constraint conditions for achieving the coordinated rotational motion of the substrate and the hydrogen barrier coating on the substrate surface are as follows: first, a reference point is set on a preset range area below the substrate, and the axis perpendicular to the substrate plane is taken as the z direction. The preset range area is coupled and constrained with the reference point, and the angular velocity of rotation along the z direction is set.

5. The method for predicting the wear performance of hydrogen barrier coating of rubber dynamic seal for hydrogen energy equipment according to claim 4, characterized in that: Set the reference point at the center of the preset range area.

6. The method for predicting the wear performance of hydrogen barrier coating of rubber dynamic seal for hydrogen energy equipment according to claim 4, characterized in that: The preset range area is the area at the bottom 10% to 15% of the substrate thickness.

7. The method for predicting the wear performance of hydrogen barrier coating of rubber dynamic seal for hydrogen energy equipment according to claim 1, characterized in that: The wear calculation of the hydrogen barrier coating on the substrate surface is performed based on the ABAQUS-UMAT user subroutine module using the UMESHMOTION subroutine combined with the ALE grid adaptive technology as follows: The subroutine uses the Archard wear model: Where V represents the wear volume of the material, σ is the wear coefficient, a is the relative sliding distance between the contact units, F is the normal load at the wear contact point, and H represents the hardness of the material. Dividing both sides of Equation (5) by the wear contact area A yields: Where h represents the wear depth, P represents the contact pressure, and λ is the specific wear rate; In order to obtain the wear depth from the contact pressure at the wear contact node, it is necessary to redefine the Archard formula at the local scale. In addition, since wear is a dynamic process, the time-varying differential Archard formula can be obtained: Where dt is the infinitesimal time for the wear part to move, v is the sliding velocity, and Equation (7) is integrated over the sliding time using Euler integration to obtain the wear depth at each wear contact node for a given wear cycle: h i+1 =h i +λvP i Δt (8) where h i and h i+1 is the wear depth of a given contact node at the i-th and i+1-th wear cycle increments, P i is the contact pressure at the contact node, Δt is the calculation time step; After calculating the wear depth of a single wear node at each wear cycle increment, the ALE mesh adaptation technology is used to reconstruct the wear node and the node positions below it to obtain a smooth wear profile.

8. The method for predicting the wear performance of hydrogen barrier coating of rubber dynamic seal for hydrogen energy equipment according to claim 7, characterized in that: The coating wear depth, surface stress and cross-sectional stress obtained by the solution are obtained to predict the coating wear performance. Specifically, in the process in which the rubber material substrate and the hydrogen barrier coating on the substrate surface come into contact with the upper wear member and produce relative displacement, the coating wear depth, surface stress and cross-sectional stress data are obtained by combining the constitutive law of the cohesive force unit and the constitutive law of wear to simultaneously calculate and obtain the coating wear performance prediction.

9. A wear performance prediction system for hydrogen energy equipment rubber dynamic seal hydrogen barrier coating, characterized in that: For implementing the method according to any one of claims 1 to 8, the system comprises the following modules: Friction and wear model building module, used to establish the friction and wear finite element model of the rubber material substrate, the hydrogen barrier coating on the substrate surface, and the wear parts above the coating; Cohesion unit establishment module, used to establish a bilinear cohesion unit with superimposed friction effect between substrate and substrate surface hydrogen barrier coating; A constraint condition establishment module is used to establish coupling constraint conditions for realizing the coordinated rotational motion of the substrate and the hydrogen barrier coating on the substrate surface; Wear calculation module, which is used to calculate the wear of hydrogen barrier coating on substrate surface based on ABAQUS-UMAT user subroutine module, using UMESHMOTION subroutine combined with ALE grid adaptive technology to obtain wear depth; The prediction module is used to obtain the solved coating wear depth, surface stress and cross-sectional stress to predict the coating wear performance.

Images