A rubber seal simulation modeling method based on considering permanent deformation compensation

Abstract

The application provides a rubber seal simulation modeling method based on considering permanent deformation compensation, and relates to the technical field of rubber seal performance evaluation. In view of the problems that the prior art cannot accurately quantify the rubber permanent deformation, the parameter calibration is complex, and the leakage prediction depends on experience, the application establishes the quantitative relationship between the short axis diameter ratio b / D and the material yield stress, and reversely corrects the material parameters through the directly measurable long axis and short axis geometric quantities; combined with the b / D temperature-time fitting function of five medium working conditions such as gas-gas and oil-oil, a sealing aging process parameter library and a statics simulation model are constructed; the contact area ratio is extracted to calculate the leakage probability, and the P leakage leakage judgment standard is determined as greater than or equal to 0.95. The application does not need complex material mechanics test, realizes accurate quantitative prediction of the leakage risk of the rubber seal after long-term service, and is suitable for aging analysis and life evaluation of key seals of industrial equipment such as hydropower stations.

Description

Technical Field

[0001] This invention relates to the field of rubber seal performance evaluation and life prediction technology, specifically to a rubber seal simulation modeling method based on permanent deformation compensation, which is particularly suitable for the aging law analysis, sealing performance simulation and leakage analysis of key rubber seals in industrial equipment such as hydropower stations under long-term service conditions. Background Technology

[0002] Rubber seals are indispensable key components in industrial equipment, and their sealing performance directly affects the safe and stable operation of the equipment. During long-term service, rubber materials will age due to factors such as temperature, medium, and stress, leading to molecular chain breakage, changes in cross-linking density, and irreversible permanent deformation. This, in turn, causes a decrease in sealing contact stress and a reduction in contact area, ultimately resulting in leakage failure.

[0003] Existing simulation modeling techniques for rubber seals that consider permanent deformation compensation mainly characterize permanent deformation through hyperelastic-plastic, viscoelastic-plastic coupled constitutive models or phenomenological constitutive models containing aging factors, combined with experimental-simulation iteration, parameter extraction, and measurable geometric quantities. b / D Parameters are calibrated using methods such as mapping, and simulation compensation is achieved through strategies such as geometry-material dual correction, with modeling completed via simulation software. However, existing technologies have the following drawbacks and limitations: (1) The mainstream simulation model generally uses linear elastic constitutive relations to describe the mechanical behavior of rubber materials. It is only applicable to ideal working conditions with small deformation and low stress, and ignores the permanent deformation characteristics of rubber materials during long-term aging. The mechanical response of the seals deviates significantly from that of the actual service process, resulting in large prediction errors of key sealing performance indicators such as leakage pressure.

[0004] (2) Some studies use viscoelastic-plastic constitutive models to characterize the viscoelastic and plastic deformation coupling characteristics of aged rubber. Although this improves the fit of the model to a certain extent, the parameter system of this type of model is complex and requires calibration of multiple key mechanical parameters that change with temperature. Not only is the calibration cycle long and the test cost high, but also the field test conditions are limited, making it difficult to quickly obtain accurate parameters that are suitable for actual working conditions, which restricts its promotion and application in engineering fields.

[0005] (3) Existing commercial simulation software has functional limitations in simulating the aging of rubber seals. It can only achieve simple aging simulations such as "modulus hardening" and cannot achieve synergistic coupling analysis of material performance degradation and geometric deformation of the sealing structure. Meanwhile, the geometric quantities commonly used in the industry for evaluating sealing performance... b / DIt is a key indicator that can be directly measured in engineering, but existing software cannot directly map this measurable geometric quantity into a dual correction coefficient of material parameters and geometric structure. As a result, the judgment of the end of the seal's life can only rely on engineering experience coefficients for estimation, resulting in poor consistency of judgment results. It is easy to have problems with life prediction being too conservative or too aggressive, which not only leads to over-design and cost waste of seals, but may also cause equipment leakage failure due to underestimating the aging risk. Summary of the Invention

[0006] The purpose of this invention is to overcome the shortcomings of the prior art and provide a simulation modeling method for rubber seals based on permanent deformation compensation. This method solves the technical problems of existing viscoelastic-plastic constitutive models, such as complex parameters, long calibration cycles, inability to achieve coordinated correction of materials and geometry, and reliance on empirical coefficients and poor consistency in determining the end of the seal life. This invention corrects material parameters in reverse using directly measurable geometric parameters, eliminating the need for complex material mechanics tests, and achieves accurate quantitative prediction of sealing performance degradation and leakage risk after long-term service of rubber seals.

[0007] To achieve the above-mentioned objectives, this invention provides a simulation modeling method for rubber seals that considers permanent deformation compensation, comprising the following steps: S1, Establish the minor axis diameter ratio b / D The quantitative relationship between yield stress and material yield stress, where b This refers to the short shaft length of the compressed rubber seal. D The original diameter of the rubber seal is... b / D Quantitative characterization of irreversible deformation during rubber aging; S2, Establish the minor shaft diameter ratio under different media conditions. b / D With temperature T Service life t The fitting function is used to calculate the corresponding service time point based on specific operating condition parameters. b / D The yield stress of the material is obtained by reverse calculation using the value and the quantitative relationship established in step S1. s s A parameter library for the aging process of sealing products was constructed, and a corresponding static simulation model was established. S3 extracts the sealing contact area ratio parameter from the static simulation model, calculates the leakage probability using the leakage probability analytical formula, and judges whether the seal has failed according to the preset judgment criteria.

[0008] Preferably, the minor axis diameter ratio is established in step S1. b / D The quantitative relationship with the yield stress of a material specifically includes: S11, the major axis length of the rubber seal after compression is obtained through simulation analysis under different compression ratios and yield stresses. a and minor axis length b ; S12, using formula (1) for the major axis length a and minor axis length b The data was fitted to obtain the fitting parameters corresponding to different compression ratios. k : (1); In formula (1) ; S13, Measure the major axis of the actual rubber seal sample. a and short axis b The original diameter is calculated according to formula (1). D Then the minor axis diameter ratio was calculated. b / D ; S14, Establish the minor axis diameter ratio b / D Compared with yield stress elastic modulus s s The quantitative relationship of / E, where s s The yield stress of the rubber material. E This refers to the elastic modulus of rubber materials.

[0009] Preferably, in step S12, when the compression rate is 25%, the fitting parameters are... k =0.2698, goodness of fit R²=99.41%; when the compression rate is 30%, the fitting parameters are... k =0.2987, goodness of fit R²=97.6%.

[0010] Preferably, in step S14, for a compression ratio of 25%, the minor axis diameter ratio is... b / D Compared with yield stress elastic modulus s s / E The quantitative relationship is as follows: (3); Preferably, in step S2, the minor axis diameter ratio is... b / D With temperature T Service life t The fitting function takes the form of formula (2): (2); In the formula: A , B 1,B 2 represents the fitting parameters. T Absolute temperature t This refers to the length of service.

[0011] Preferably, for different media conditions, the fitting parameters in formula (2) are as follows: Gas operating conditions: A =0.2113, B 1 = 12.5315 B 2=7716, goodness of fit R²=0.91; Oil-oil operating condition: A =0.2136, B 1 = 11.0775 B 2=7137, goodness of fit R²=0.73; Water-oil working condition: A =0.2218, B 1 = 17.3339 B 2=8374, goodness of fit R²=0.80; Water and air conditions: A =0.2210, B 1 = 19.5198 B 2=9126, goodness of fit R²=0.90; Oil and gas operating conditions: A =0.2137, B 1 = 9.9902, B =2=6554, goodness of fit R²=0.70.

[0012] Preferably, in step S3, the extracted sealing contact area ratio parameter is the ratio of the actual contact area to the nominal contact area.

[0013] Preferably, in step S3, the analytical expression for the leakage probability is: ; in, k 1 = 21.97 ， This represents the actual contact area. This refers to the nominal contact area.

[0014] Preferably, in step S3, the preset judgment criterion is: when the leakage probability... P leakage If the value is ≥0.95, the seal is determined to have a probability of leakage; otherwise, the specific leakage probability value is output.

[0015] Preferably, the rubber seal is an O-ring, suitable for sealing systems in hydropower stations, petrochemical plants, and shipbuilding equipment.

[0016] The present invention has the following beneficial effects: 1. This invention constructs a minor axis diameter ratio b / D The quantitative relationship between the yield stress and the material is obtained by inversely correcting the material parameters through directly measurable major and minor axis geometric quantities. This eliminates the need for complex material mechanics tests, greatly simplifies the parameter acquisition process, and improves the convenience and accuracy of parameter correction.

[0017] 2. The present invention is based on b / D As a quantitative evaluation index for irreversible deformation, combined with fitting functions under different media conditions, it can accurately capture the deformation patterns under different temperatures and service times, solving the technical pain point of difficulty in quantifying deformation in traditional technologies.

[0018] 3. This invention, by constructing an aging process parameter library and a static simulation model, and combining contact area ratio calculation with an analytical expression for leakage probability, can quantitatively output the leakage probability and clarify... P leakage A value of ≥0.95 is used as the criterion for determining leakage in terms of probability, which improves the scientific rigor and practicality of the sealing performance evaluation.

[0019] 4. This invention establishes dedicated fitting functions for five common media conditions: gas-gas, oil-oil, water-oil, water-gas, and oil-gas. These functions can meet the sealing simulation requirements under different service environments and have a wide range of applications. Attached Figure Description

[0020] The present invention will be further described below with reference to the accompanying drawings and embodiments.

[0021] Figure 1 This is a flowchart of the simulation modeling method for rubber seals based on permanent deformation compensation, which is the basis of this invention. Figure 2 The minor axis diameter ratio of the material of this invention b / D Test data points for temperature and time (a total of five operating conditions); Figure 3 This is a simulation cloud map of the irreversible deformation under the corresponding compression state in the embodiments of the present invention; Figure 4 This is a cross-sectional view of a rectangular trench in an embodiment of the present invention; Figure 5 This is a graph showing the relationship between leakage probability and contact area ratio in an embodiment of the present invention. Detailed Implementation

[0022] The embodiments of the present invention will be further described below with reference to the accompanying drawings.

[0023] Example 1: Simulation of rectangular groove O-ring rubber seal under water-oil medium conditions Step 1: Establish the minor axis diameter ratio b / D Quantitative relationship with material yield stress; 1. Through finite element simulation analysis, the major axis length of the O-ring after compression is obtained under different yield stresses at two commonly used engineering compression ratios of 25% and 30%. a and minor axis length b data; 2. The above simulation data are fitted using formula (1): (1); Formula (1) can be transformed into a practical engineering formula for calculating the original diameter: ; The fitting results show that when the compression rate is 25%, the fitting parameters are... k =0.2698, goodness of fit R²=99.41%; when the compression rate is 30%, the fitting parameters are... k =0.2987, goodness of fit R²=97.6%; 3. For any actual rubber seal sample, only the length of its major axis after compression needs to be measured. a and minor axis length b The original diameter can then be calculated using the formula described above. D Then the minor axis diameter ratio was calculated. b / D This allows for the quantitative characterization of irreversible deformation during the aging process of rubber. 4. Establish the minor axis diameter ratio b / D Compared with yield stress elastic modulus s s / E The quantitative relationship, for a compression rate of 25% (this compression rate is used in this embodiment), is as follows: (3); in E The elastic modulus of rubber materials can be measured in one step through conventional material mechanics tests.

[0024] Step 2: Construct an aging process parameter library and a static simulation model; 1. Using formula (2) and combined with accelerated aging test data, establish the minor axis diameter ratio under different media conditions. b / D With temperature T Service life t Fitting function: (2); in, A , B 1, B 2 represents the fitting parameters. T Absolute temperature (unit: K). t Service time (unit: days); 2. In this embodiment, the medium pair is water-oil, and its fitting parameters are: A =0.2218, B 1 = 17.3339 B 2=8374, goodness of fit R²=0.80; Table 1 Fitting functions for five media conditions

[0025] 3. Given that the average annual service temperature of the seal is 15.21℃, the elastic modulus of the rubber material is... E =7.8MPa, calculate service life t Material parameters at 1.58 years (578.49 days): Temperature conversion: T =15.21+273.15=288.36K; Substitute into formula (2) to calculate the minor axis diameter ratio: ; Substitute into formula (3) to calculate the ratio of yield stress to elastic modulus: ; Calculate the yield stress: s s =0.234 × 7.8 = 1.8252 MPa; 4. Following the method described above, calculate the service time for 0 years, 1 year, 2 years, 3 years, 4 years, and 5 years in sequence. b / D and s s Values ​​were used to construct a parameter library for the aging process of the sealing product. 5. Establish a two-dimensional axisymmetric static simulation model of the rectangular groove O-ring rubber seal, and incorporate material parameters (elastic modulus) at different service time points. E and yield stress s s Substitute each into the model and perform finite element simulation calculations.

[0026] Step 3: Leakage probability calculation and failure determination; 1. Extract the actual contact area between the seal and the groove from the above static simulation model. S contact and nominal contact areaS nominal Calculate the contact area ratio S contact / S nominal ; 2. Calculate the leakage probability using the aforementioned analytical formula: (4); Among them, parameters k 1 = 21.97; 3. In this embodiment, the simulation yielded a contact area ratio of 0.694401 after 1.58 years of service. Substituting this into formula (4), we can calculate: ; 4. Based on preset judgment criteria: when the leakage probability... P leakage If the value is ≥0.95, the seal is determined to have a probability of leakage; otherwise, the specific leakage probability value is output.

[0027] In this embodiment, the leakage probability is 0.013775 < 0.95, therefore it is determined that the seal has not leaked in terms of probability and can continue to be used.

[0028] Example 2: Simulation of rectangular groove O-ring rubber seal under gas-to-gas medium working conditions This embodiment illustrates the specific implementation of the present invention under gas-to-gas medium conditions.

[0029] Step 1: Establish the minor axis diameter ratio b / D Quantitative relationship with material yield stress; This step is exactly the same as in Example 1, and is as follows: 1. Through finite element simulation analysis, the major axis length of the O-ring after compression is obtained under different yield stresses at two commonly used engineering compression ratios of 25% and 30%. a and minor axis length b data; 2. The above simulation data are fitted using formula (1): (1); Formula (1) can be transformed into a practical engineering formula for calculating the original diameter: ; The fitting results show that when the compression rate is 25%, the fitting parameters are... k =0.2698, goodness of fit R²=99.41%; when the compression rate is 30%, the fitting parameters are... k =0.2987, goodness of fit R²=97.6%; 3. For any actual rubber seal sample, only the length of its major axis after compression needs to be measured. a and minor axis length b The original diameter can then be calculated using the formula described above. D Then the minor axis diameter ratio was calculated. b / D This allows for the quantitative characterization of irreversible deformation during the aging process of rubber. 4. Establish the minor axis diameter ratio b / D Compared with yield stress elastic modulus s s / E The quantitative relationship, for a compression rate of 25% (this compression rate is used in this embodiment), is as follows: (3); in E The elastic modulus of rubber materials can be measured in one step through conventional material mechanics tests.

[0030] Step 2: Construct an aging process parameter library and a static simulation model; 1. Using formula (2) and combined with accelerated aging test data, establish the minor axis diameter ratio under different media conditions. b / D With temperature T Service life t Fitting function: (2); in, A , B 1. B 2 represents the fitting parameters. T Absolute temperature (unit: K). t Service time (unit: days); 2. In this embodiment, the medium pair is a gas-gas operating condition fitting function, and its fitting parameters are: A =0.2113, B 1 = 12.5315 B 2=7716, goodness of fit R²=0.91; 3. Given that the average annual service temperature of this seal is 25℃, the elastic modulus of the rubber material is... E =7.5MPa, calculate service life t Material parameters at 3 years (1095 days): Temperature conversion: T =25 + 273.15 = 298.15K; Substitute into formula (2) to calculate the minor axis diameter ratio: ; Substitute into formula (3) to calculate the ratio of yield stress to elastic modulus: ; Calculate the yield stress: s s =0.1857 × 7.5 = 1.393 MPa; 4. Following the method described above, calculate the b / D ratio for service periods of 0 years, 1 year, 2 years, 4 years, 5 years, and 6 years in sequence. s s Values ​​were used to construct a parameter library for the aging process of the sealing product. 5. Establish a two-dimensional axisymmetric static simulation model of the rectangular groove O-ring rubber seal, and incorporate material parameters (elastic modulus E and yield stress) at different service time points. s s Substitute each into the model and perform finite element simulation calculations.

[0031] Step 3: Leakage probability calculation and failure determination; 1. Extract the actual contact area between the seal and the groove from the above static simulation model. S contact and nominal contact area S nominal Calculate the contact area ratio S contact / S nominal =0.691593; 2. Calculate the leakage probability using the analytical formula for leakage probability: (4); Among them, parameters k 1 = 21.97; 3. In this embodiment, substituting into formula (4) yields the following result: =0.014639; 4. Based on the aforementioned preset judgment criteria: when the leakage probability... P leakage If the value is ≥0.95, the seal is determined to have a probability of leakage; otherwise, the specific leakage probability value is output.

[0032] In this embodiment, the leakage probability is 0.014639 < 0.95, therefore it is determined that the seal has not leaked in terms of probability and can continue to be used.

[0033] Example 3 (Simulation of rectangular groove O-ring rubber seal under water and air medium conditions) This embodiment illustrates the specific implementation of the present invention under water-air medium conditions.

[0034] Step 1: Establish the minor axis diameter ratio b / D Quantitative relationship with material yield stress; This step is exactly the same as in Example 1.

[0035] Step 2: Construct an aging process parameter library and a static simulation model; 1. Using the aforementioned formula (2) and combined with accelerated aging test data, establish the minor axis diameter ratio under different media conditions. b / D Fitting function with temperature T and service time t: ; in, A , B 1, B 2 represents the fitting parameters. T t represents absolute temperature (unit: K), and t represents service time (unit: days). 2. In this embodiment, the medium pair is water-air condition, and its fitting parameters are: A =0.2210, B 1 = 19.5198 B 2=9126, goodness of fit R²=0.90; 3. Given that the average annual service temperature of the seal is 30℃ and the elastic modulus of the rubber material is E=8.0MPa, calculate the service life. t Material parameters at 5 years (1825 days): Temperature conversion: T =30 + 273.15 = 303.15K; Substitute into formula (2) to calculate the minor axis diameter ratio: ; Substitute into formula (3) to calculate the ratio of yield stress to elastic modulus: ; Calculate the yield stress: s s =0.1403 × 8.0 = 1.122 MPa 4. Following the method described above, calculate the service time for 0 years, 1 year, 2 years, 3 years, 4 years, and 6 years in sequence. b / D and s s Values ​​were used to construct a parameter library for the aging process of the sealing product. 5. Establish a two-dimensional axisymmetric static simulation model of the rectangular groove O-ring rubber seal, and incorporate material parameters (elastic modulus) at different service time points. E and yield stress ss Substitute each into the model and perform finite element simulation calculations.

[0036] Step 3: Leakage probability calculation and failure determination; 1. Extract the actual contact area between the seal and the groove from the above static simulation model. S contact and nominal contact area S nominal Calculate the contact area ratio S contact / S nominal ; 2. Calculate the leakage probability using the aforementioned analytical formula: (4); Among them, parameters k 1 = 21.97; 3. In this embodiment, substituting into formula (4) yields the following result: ; 4. Based on the aforementioned preset judgment criteria: when the leakage probability... P leakage If the value is ≥0.95, the seal is determined to have a probability of leakage; otherwise, the specific leakage probability value is output.

[0037] In this embodiment, the leakage probability is 0.328 < 0.95, therefore it is determined that the seal has not leaked probabilistically and can continue to be used.

[0038] Example 4: Simulation of rectangular groove O-ring rubber seal under oil and gas medium conditions This embodiment illustrates the specific implementation of the present invention under oil and gas medium conditions, as well as the process for determining when the seal reaches the end of its service life.

[0039] Step 1: Establish the minor axis diameter ratio b / D Quantitative relationship with material yield stress; This step is exactly the same as in Example 1.

[0040] Step 2: Construct an aging process parameter library and a static simulation model; 1. Using the aforementioned formula (2) and combined with accelerated aging test data, establish the minor axis diameter ratio under different media conditions. b / D With temperature T Service life t Fitting function: ; in, A ,B 1, B 2 represents the fitting parameters. T Absolute temperature (unit: K). t Service time (unit: days); 2. In this embodiment, the medium pair is an oil and gas working condition, and its fitting parameters are: A =0.2137, B 1 = 9.9902, B 2=6554, goodness of fit R²=0.70; 3. Given that the average annual service temperature of this seal is 40℃, the elastic modulus of the rubber material is... E =7.2MPa, calculate service time t Material parameters at 8 years (2920 days): Temperature conversion: T =40 + 273.15 = 313.15K Substitute into formula (2) to calculate the minor axis diameter ratio: ; Substitute into formula (3) to calculate the ratio of yield stress to elastic modulus: ; Calculate the yield stress: s s =0.0167 × 7.2 = 0.120 MPa 4. Following the method described above, calculate the b / D ratio for service years of 0, 2, 4, 6, 7, and 8 years, respectively. s s Values ​​were used to construct a parameter library for the aging process of the sealing product. 5. Establish a two-dimensional axisymmetric static simulation model of the rectangular groove O-ring rubber seal, and incorporate material parameters (elastic modulus) at different service time points. E and yield stress s s Substitute each into the model and perform finite element simulation calculations.

[0041] Step 3: Leakage probability calculation and failure determination; 1. Extract the actual contact area between the seal and the groove from the above static simulation model. S contact and nominal contact area S nominal Calculate the contact area ratio S contact / S nominal ; 2. Calculate the leakage probability using the aforementioned analytical formula: (4); Among them, parameters k 1 = 21.97; 3. In this embodiment, substituting into formula (4) yields the following result: ; 4. Based on the aforementioned preset judgment criteria: when the leakage probability... P leakage If the value is ≥0.95, the seal is determined to have a probability of leakage; otherwise, the specific leakage probability value is output.

[0042] In this embodiment, the leakage probability is 0.962 ≥ 0.95, therefore it is determined that the seal is leaking probabilistically and it is recommended to replace it immediately.

[0043] The above are all preferred embodiments of the present invention and are not intended to limit the scope of protection of the present invention. Therefore, all equivalent changes made in accordance with the structure, shape and principle of the present invention should be covered within the scope of protection of the present invention.

Claims

1. A simulation modeling method for rubber seals considering permanent deformation compensation, characterized in that, Includes the following steps: S1, Establish the minor axis diameter ratio b / D The quantitative relationship between yield stress and material yield stress, where b This refers to the short shaft length of the compressed rubber seal. D The original diameter of the rubber seal is... b / D Quantitative characterization of irreversible deformation during rubber aging; S2, Establish the minor shaft diameter ratio under different media conditions. b / D With temperature T Service life t The fitting function is used to calculate the corresponding service time point based on specific operating condition parameters. b / D The yield stress of the material is obtained by reverse calculation using the value and the quantitative relationship established in step S1. σ s A parameter library for the aging process of sealing products was constructed, and a corresponding static simulation model was established. S3 extracts the sealing contact area ratio parameter from the static simulation model, calculates the leakage probability using the leakage probability analytical formula, and judges whether the seal has failed according to the preset judgment criteria.

2. The method for simulating and modeling rubber seals based on permanent deformation compensation according to claim 1, characterized in that, In step S1, the minor axis diameter ratio is established. b / D The quantitative relationship with the yield stress of a material specifically includes: S11, the major axis length of the rubber seal after compression is obtained through simulation analysis under different compression ratios and yield stresses. a and minor axis length b ; S12, using formula (1) for the major axis length a and minor axis length b The data was fitted to obtain the fitting parameters corresponding to different compression ratios. k : （1）； Official (1) ; S13, Measure the major axis of the actual rubber seal sample. a and short axis b The original diameter is calculated according to formula (1). D Then the minor axis diameter ratio was calculated. b / D ; S14, Establish the minor axis diameter ratio b / D Compared with yield stress elastic modulus σ s The quantitative relationship of / E, where σ s The yield stress of the rubber material. E This refers to the elastic modulus of rubber materials.

3. The method for simulating and modeling rubber seals based on permanent deformation compensation according to claim 2, characterized in that, In step S12, when the compression rate is 25%, the fitting parameters are... k =0.2698, goodness of fit R²=99.41%; when the compression rate is 30%, the fitting parameters are... k =0.2987, goodness of fit R²=97.6%.

4. The method for simulating and modeling rubber seals based on permanent deformation compensation according to claim 2, characterized in that, In step S14, for a compression ratio of 25%, the minor axis diameter ratio is... b / D Compared with yield stress elastic modulus σ s / E The quantitative relationship is as follows: （3）；。 5. The method for simulation modeling of rubber seals considering permanent deformation compensation according to claim 1, characterized in that, In step S2, the minor axis diameter ratio b / D With temperature T Service life t The fitting function takes the form of formula (2): （2）； In the formula: A , B 1, B 2 represents the fitting parameters. T Absolute temperature t This refers to the length of service.

6. The method for simulation modeling of rubber seals considering permanent deformation compensation according to claim 5, characterized in that, For different media conditions, the fitting parameters in formula (2) are as follows: Gas operating conditions: A =0.2113, B 1 = 12.5315 B 2=7716, goodness of fit R²=0.91; Oil-oil operating condition: A =0.2136, B 1 = 11.0775 B 2=7137, goodness of fit R²=0.73; Water-oil working condition: A =0.2218, B 1 = 17.3339 B 2=8374, goodness of fit R²=0.80; Water and air conditions: A =0.2210, B 1 = 19.5198 B 2=9126, goodness of fit R²=0.90; Oil and gas operating conditions: A =0.2137, B 1 = 9.9902, B =2=6554, goodness of fit R²=0.

70.

7. The method for simulating and modeling rubber seals based on permanent deformation compensation according to claim 1, characterized in that, In step S3, the extracted sealing contact area ratio parameter is the ratio of the actual contact area to the nominal contact area.

8. The method for simulating and modeling rubber seals based on permanent deformation compensation according to claim 1, characterized in that, In step S3, the analytical expression for the leakage probability is: ； in, k 1 = 21.97 ， This represents the actual contact area. This refers to the nominal contact area.

9. The simulation modeling method for rubber seals based on permanent deformation compensation according to claim 8, characterized in that, In step S3, the preset judgment criterion is: when the leakage probability... P leakage If the value is ≥0.95, the seal is determined to have a probability of leakage; otherwise, the specific leakage probability value is output.

10. The simulation modeling method for rubber seals based on permanent deformation compensation according to claim 1, characterized in that, The rubber seal is an O-ring, suitable for sealing systems in hydropower stations, petrochemical plants, and shipbuilding equipment.

Images