Nonlinear Finite Element Modeling Method for Creep Analysis of Rubber Bearings

Abstract

The present invention discloses a non-linear finite element modeling method for creep analysis of rubber bearings, including: establishing a non-linear finite element model of a rubber bearing and performing stiffness and damping simulation calculations; conducting stiffness and damping tests on the rubber bearing and performing static calibration analysis on the non-linear finite element model of the rubber bearing according to the test results; when the static calibration analysis meets the first preset requirement, based on the non-linear finite element model of the rubber bearing that meets the preset requirement, establishing a thermosetting coupling non-linear model of the rubber bearing and performing steady-state and transient thermodynamics simulation analyses; conducting stress and temperature tests on the rubber bearing and performing dynamic calibration analysis on the non-linear finite element model of the rubber bearing according to the test results; when the dynamic calibration analysis meets the second preset requirement, setting boundary conditions and load conditions and performing creep simulation calculations on the rubber bearing. The present invention improves the creep simulation accuracy of rubber bearings and the practicability and operability of finite element simulation for creep calculation of rubber bearings.

Description

Technical Field

[0001] The present invention relates to the technical field of finite element simulation, and in particular to the finite element simulation of rubber bearings. Specifically, it is a non-linear finite element modeling method for creep analysis of rubber bearings. Background Art

[0002] The intermediate support of a vehicle drive shaft belongs to a kind of rubber bearing, which mainly bears the centrifugal force load of the drive shaft and plays a role in attenuating the bending vibration and torsional vibration of the drive shaft. It is a key system to improve the comfort performance of the vehicle. However, rubber is a viscoelastic material with highly non-linear mechanical properties. When rubber bearings work in an alternating load environment for a long time, creep phenomena often occur. Once creep occurs in the rubber bearing of the vehicle drive shaft, it will not only affect the vibration reduction and noise reduction effect of the rubber bearing, but also amplify the vibration of the drive shaft, making the comfort performance of the vehicle worse. Therefore, creep analysis and control of rubber bearings are the basis for reducing the vibration of the drive shaft and improving the comfort of the vehicle.

[0003] At present, the creep analysis of rubber bearings mainly adopts the finite element simulation method, that is, using ABAQUS software to establish a non-linear finite element model of the rubber bearing and setting creep loads for creep variable simulation calculation. The disadvantages of this method include: (1) lack of methods for correcting and statically calibrating the parameters of the rubber constitutive model; (2) lack of dynamic calibration methods for rubber bearings; (3) no consideration of the influence of vibration and temperature on creep. The above factors lead to large errors in the existing simulation methods, and the simulation results cannot truly reflect the creep characteristics of rubber bearings.

[0004] For example, a method for predicting the creep characteristics of a rubber isolator disclosed in the Chinese patent with the publication number CN112507595A on March 16, 2021, includes: establishing a geometric model of the rubber isolator; meshing the rubber part of the rubber isolator; adopting a hyperelastic-nonlinear viscoelastic superposition constitutive model and using a constitutive model parameter identification method to complete the material property setting; constraining the six degrees of freedom of all nodes of the metal outer tube of the rubber isolator; associating all nodes of the metal inner tube of the rubber isolator with the central node of the rubber isolator, and applying a force with a constant magnitude along the central node to complete the boundary condition setting; result post-processing: recording the change of the displacement of the central node of the rubber isolator over time to obtain the creep variable curve of the central node, and evaluating the creep characteristics of the rubber isolator according to the creep variable curve of the central node of the rubber isolator. This patent application predicts the creep characteristics of rubber isolators with different structures under different load conditions by means of finite elements, but its practicability and operability are relatively low.

[0005] In order to solve the shortcomings of existing finite element simulation methods, improve the creep simulation accuracy of rubber bearings, and enhance the practicability and operability of the finite element simulation method for rubber bearing creep calculation, the present invention proposes a non-linear finite element modeling method for rubber bearing creep analysis. Summary of the Invention

[0006] Aiming at the deficiencies of the prior art, the present invention provides a non-linear finite element modeling method for rubber bearing creep analysis to improve the creep simulation accuracy of rubber bearings and enhance the practicability and operability of the finite element simulation method for rubber bearing creep calculation.

[0007] The present invention is realized through the following technical solutions:

[0008] A non-linear finite element modeling method for rubber bearing creep analysis includes:

[0009] Establish a non-linear finite element model of a rubber bearing and conduct stiffness and damping simulation calculations;

[0010] Conduct stiffness and damping tests on the rubber bearing and perform static calibration analysis on the non-linear finite element model of the rubber bearing according to the test results;

[0011] When the static calibration analysis meets the first preset requirement, based on the non-linear finite element model of the rubber bearing that meets the preset requirement, establish a thermo-solid coupling non-linear model of the rubber bearing and conduct steady-state and transient thermodynamics simulation analyses;

[0012] Conduct stress and temperature tests on the rubber bearing and perform dynamic calibration analysis on the non-linear finite element model of the rubber bearing according to the test results;

[0013] When the dynamic calibration analysis meets the second preset requirement, set boundary conditions and load conditions and conduct rubber bearing creep simulation calculations.

[0014] The above technical solution first conducts stiffness and damping simulation calculations on the established non-linear finite element model of the rubber bearing, then conducts stiffness and damping tests on the rubber bearing, and calibrates the established non-linear finite element model of the rubber bearing according to the test results so that the coincidence degree between the simulation results of the calibrated model and the test results reaches the preset requirement; then, based on the calibrated non-linear finite element model of the rubber bearing, establish a thermo-solid coupling non-linear model of the rubber bearing and conduct steady-state and transient thermodynamics simulation analyses and stress and temperature tests on the rubber bearing. Similarly, calibrate the established thermo-solid coupling non-linear model of the rubber bearing according to the test results so that the coincidence degree between the simulation results of the calibrated model and the test results reaches the preset requirement; when both the non-linear finite element model of the rubber bearing and the thermo-solid coupling non-linear model of the rubber bearing based on this model meet the preset requirements, through boundary condition and load condition settings, conduct rubber bearing creep simulation calculations.

[0015] The above technical solution combines the static calibration, dynamic calibration of the rubber bearing, and the influence of vibration and temperature on creep, improving the creep simulation accuracy of the rubber bearing, and having practicability and operability.

[0016] As a further technical solution, a non-linear finite element model of the rubber bearing is established, and stiffness and damping simulation calculations are carried out, further including:

[0017] Establish a finite element mesh model of the bearing, a finite element mesh model of the rubber part, and a finite element mesh model of the bracket;

[0018] Establish the contact relationship between the bearing roller and the inner and outer rings of the bearing, and set the friction coefficient;

[0019] Establish a fixed constraint between the outer ring of the bearing and the rubber part;

[0020] Establish the contact relationship between the bracket and the rubber part, and set the friction coefficient;

[0021] Set the parameters of the Mooney-Rivlin constitutive model of the rubber part;

[0022] Establish a node at the center point of the bearing, and use the rigid connection element rb2 to connect the node at the center point of the bearing to the nodes on the inner surface of the bearing inner ring, and establish fixed constraints at the positions of two bolts on the bracket;

[0023] Apply a forced displacement load to the node at the center point of the bearing, and carry out the simulation calculations of the Z-direction stiffness and Y-direction stiffness of the rubber bearing;

[0024] Apply an impulse displacement load to the node at the center point of the bearing, and carry out the simulation calculations of the Z-direction damping and Y-direction damping of the rubber bearing.

[0025] As a further technical solution, carry out the stiffness and damping tests of the rubber bearing, and conduct static calibration analysis on the non-linear finite element model of the rubber bearing according to the test results, further including:

[0026] Test the Y-direction and Z-direction stiffness, Y-direction and Z-direction damping of the rubber bearing;

[0027] Establish a parametric model of the Mooney-Rivlin constitutive model, and set the rubber material property parameters c10, c01, and c11 as design variables;

[0028] Set the Y-direction and Z-direction stiffness test data, Y-direction and Z-direction damping test data of the rubber bearing as the objective function;

[0029] Determine the values of c10, c01, and c11 through optimization simulation calculations;

[0030] Modify the Mooney-Rivlin constitutive model of the rubber parts according to the optimization results of c10, c01 and c11;

[0031] Perform static simulation calculations on the non-linear finite element model of the rubber bearing after static calibration to verify the accuracy of the model.

[0032] As a further technical solution, establish a thermosetting coupling non-linear model of the rubber bearing and conduct steady-state and transient thermodynamics simulation analysis, further including:

[0033] Set the element type of the non-linear finite element model of the rubber bearing to solid70;

[0034] Set the heat generation rate and heat conduction coefficient on the contact surfaces between the bearing rollers and the inner and outer rings of the bearing, set the heat conduction coefficient on the connection surface between the outer ring of the bearing and the rubber parts, and set the heat conduction coefficient on the contact surface between the bracket and the rubber parts;

[0035] Apply the time-varying forces Fy and Fz at the node position of the center point of the bearing. Fy and Fz represent the forces in the y-direction and z-direction respectively;

[0036] Perform simulation calculations on frictional heat generation and heat transfer.

[0037] As a further technical solution, conduct stress and temperature tests on the rubber bearing, and perform dynamic calibration analysis on the non-linear finite element model of the rubber bearing according to the test results, further including:

[0038] Install strain gauges on the rubber parts to measure the stress test data of the rubber parts;

[0039] Arrange temperature sensors on the rubber parts to measure the temperature test data of the rubber parts;

[0040] Set the friction coefficients on the contact surfaces between the bearing rollers and the inner and outer rings of the bearing, and the friction coefficient on the contact surface between the bracket and the rubber parts as design variables;

[0041] Set the stress test data and temperature test data of the rubber bearing as the objective function;

[0042] Determine the values of the friction coefficients on each contact surface through optimization simulation calculations;

[0043] Modify the thermosetting coupling non-linear model of the rubber bearing according to the optimization results of the friction coefficients on each contact surface;

[0044] Perform thermodynamics simulation calculations on the non-linear finite element model of the rubber bearing after dynamic calibration to verify the accuracy of the model.

[0045] As a further technical solution, set boundary conditions and load conditions, and conduct creep simulation calculations on the rubber bearing, further including:

[0046] Apply the time-domain data of Fy and Fz at the center point node position of the non-linear finite element model of the rubber bearing after static calibration.

[0047] Apply the temperature field simulation results of the thermo-solid coupling model of the rubber bearing to the non-linear finite element model of the rubber bearing.

[0048] Define the creep model and coefficients of the rubber material, set the tolerance of the creep strain, and perform the display simulation solution calculation.

[0049] Plot the displacement change curve of a certain point on the rubber part to check the creep amount within the specified time.

[0050] As a further technical solution, when the static calibration analysis does not meet the first preset requirement, reset the Mooney-Rivlin constitutive model parameters of the rubber part and re-perform the simulation calculation of the stiffness and damping of the rubber bearing until the model after static calibration meets the preset requirement.

[0051] As a further technical solution, when the dynamic calibration analysis does not meet the second preset requirement, re-establish the thermo-solid coupling non-linear model of the rubber bearing and perform the steady-state and transient thermodynamics simulation analysis until the model after dynamic calibration meets the preset requirement.

[0052] As a further technical solution, the first preset requirement and the second preset requirement refer to the requirement of the coincidence degree between the model simulation results and the test results.

[0053] Compared with the prior art, the beneficial effects of the present invention are as follows:

[0054] The present invention combines finite element simulation and test technology, and proposes a non-linear finite element modeling method for creep analysis of rubber bearings. The coincidence degree between the stiffness and damping simulation data calculated by the non-linear finite element model of the rubber bearing established by this method and the test data is greater than 90%, and the coincidence degree between the thermo-solid coupling stress and temperature simulation data and the test data is greater than 90%. Finally, the creep amount of the rubber bearing is calculated. The present invention synthesizes the static calibration, dynamic calibration of the rubber bearing and the influence of vibration and temperature on creep, improves the creep simulation accuracy of the rubber bearing, and has practicability and operability. Description of the Drawings

[0055] Figure 1 It is a flow chart of the non-linear finite element modeling method for creep analysis of rubber bearings according to an embodiment of the present invention.

[0056] Figure 2 It is a schematic structural diagram of a rubber bearing according to an embodiment of the present invention.

[0057] Figure 3Comparison diagram of simulation data and test data of the Y-direction stiffness of the rubber bearing according to an embodiment of the present invention.

[0058] Figure 4 Comparison diagram of simulation data and test data of the Z-direction stiffness of the rubber bearing according to an embodiment of the present invention.

[0059] Figure 5 Comparison diagram of simulation data and test data of the Y-direction damping of the rubber bearing according to an embodiment of the present invention.

[0060] Figure 6 Comparison diagram of simulation data and test data of the Z-direction damping of the rubber bearing according to an embodiment of the present invention.

[0061] Figure 7 Comparison diagram of simulation data and test data of the stress of the rubber bearing according to an embodiment of the present invention.

[0062] Figure 8 Comparison diagram of simulation data and test data of the temperature of the rubber bearing according to an embodiment of the present invention.

[0063] Figure 9 Creep simulation curve diagram of the rubber bearing according to an embodiment of the present invention.

[0064] Description of the drawings: 1. Bearing bracket; 2. Rubber part; 3. Bearing; 4. Temperature sensor; 5. Strain gauge. Detailed implementation manners

[0065] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the drawings. Obviously, the described embodiments are only a part of the embodiments of the present invention, rather than all of the embodiments. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present invention without creative efforts shall fall within the scope of protection of the present invention.

[0066] In the description of the present invention, it should be noted that, unless otherwise clearly defined and limited, the terms "installation", "connection", and "connection" should be understood in a broad sense. For example, it can be a fixed connection, a detachable connection, or an integral connection; it can be a mechanical connection, an electrical connection, or a communication with each other; it can be directly connected, or indirectly connected through an intermediate medium, and it can be the internal connection of two components or the interaction relationship between two components. For those of ordinary skill in the art, the specific meanings of the above terms in the present invention can be understood according to specific situations.

[0067] See Figure 1 And Figure 2 , the present invention provides a non-linear finite element modeling method for creep analysis of rubber bearings, including the following steps:

[0068] (1) Establish a non-linear finite element model of the rubber bearing and conduct stiffness and damping simulation calculations

[0069] Establish a finite element mesh model for bracket 1, a finite element mesh model for rubber part 2, and a finite element mesh model for bearing 3. Establish the contact relationship between the bearing rollers and the inner and outer rings of the bearing and set the friction coefficient. Establish a fixed constraint between the outer ring of the bearing and the rubber part. Establish the contact relationship between the bracket and the rubber part and set the friction coefficient. Set the parameters of the Mooney-Rivlin constitutive model for the rubber part. Establish a node (Nc) at the center point of the bearing, and use rb2 to connect the center point Nc of the bearing to the nodes on the inner surface of the inner ring of the bearing, and establish fixed constraints at the positions of two bolts (A and B) on the bracket. Apply a forced displacement load at the center point Nc of the bearing to conduct simulation calculations for the Z-direction stiffness and Y-direction stiffness of the rubber bearing. Apply an impulse displacement load at the node at the center point of the bearing to conduct simulation calculations for the Z-direction damping and Y-direction damping of the rubber bearing.

[0070] (2) Conduct stiffness and damping tests on the rubber bearing and conduct static calibration analysis of the finite element model of the rubber bearing

[0071] Test the Y-direction and Z-direction stiffness, and Y-direction and Z-direction damping of the rubber bearing. Establish a parametric model of the Mooney-Rivlin constitutive model and set c10, c01, and c11 as design variables. Set the test data of the Y-direction and Z-direction stiffness and the Y-direction and Z-direction damping of the rubber bearing as the objective function. Determine the values of c10, c01, and c11 through optimized simulation calculations. Modify the Mooney-Rivlin constitutive model of the rubber part according to the optimization results of c10, c01, and c11. Conduct static simulation calculations on the finite element model of the rubber bearing after static calibration to verify the accuracy of the model.

[0072] The comparison results between the simulation data and the test data of the Y-direction stiffness of the rubber bearing are as Figure 3 shown. The average value of the simulation data is 74.3 N / mm, and the average value of the test data is 76.5 N / mm. The degree of coincidence between the two is 97.1%.

[0073] The comparison results between the simulation data and the test data of the Z-direction stiffness of the rubber bearing are as Figure 4 shown. The average value of the simulation data is 112.6 N / mm, and the average value of the test data is 115.4 N / mm. The degree of coincidence between the two is 97.6%.

[0074] The comparison results between the simulation data and the test data of the Y-direction damping coefficient of the rubber bearing are as Figure 5 shown. The average value of the simulation data is 0.3, and the average value of the test data is 0.32. The degree of coincidence between the two is 93.8%.

[0075] The comparison results of the simulation data and test data of the Z-direction damping coefficient of the rubber bearing are as follows Figure 6 shown. The average value of the simulation data is 0.38, and the average value of the test data is 0.41. The degree of coincidence between the two is 92.7%.

[0076] Through comparative analysis, it can be obtained that the degree of coincidence between the stiffness and damping simulation data and test data of the rubber bearing finite element model after static calibration is greater than 90%, indicating that the established rubber bearing finite element model is accurate.

[0077] If the degree of coincidence between the simulation data and the test data does not meet the preset requirements, reset the parameters of the Mooney-Rivlin constitutive model of the rubber parts, and repeat the above steps until the model meets the preset requirements.

[0078] (3) Establish a thermosetting coupling nonlinear model of the rubber bearing, and conduct steady-state and transient thermodynamics simulation analysis

[0079] Set the element type of the rubber bearing finite element model to solid70. Set the heat generation rate and heat conduction coefficient on the contact surfaces between the bearing rollers and the inner and outer rings of the bearing, set the heat conduction coefficient on the connection surface between the outer ring of the bearing and the rubber parts, and set the heat conduction coefficient on the contact surface between the bracket and the rubber parts. Apply the time-varying forces Fy and Fz at the node position of the center point of the bearing; conduct simulation calculations of frictional heat generation and heat transfer.

[0080] (4) Conduct stress and temperature tests on the rubber bearing, and conduct dynamic calibration analysis of the thermosetting coupling nonlinear model of the rubber bearing

[0081] Install strain gauges at position C on the rubber parts to measure the stress test data of the rubber parts. Arrange temperature sensors at position D on the rubber parts to measure the temperature test data of the rubber parts. Set the friction coefficients on the contact surfaces between the bearing rollers and the inner and outer rings of the bearing, and the friction coefficient on the contact surface between the bracket and the rubber parts as design variables. Set the stress test data and temperature test data of the rubber bearing as the objective function. Determine the values of the friction coefficients on each contact surface through optimized simulation calculations. Modify the thermosetting coupling nonlinear model of the rubber bearing according to the optimization results of the friction coefficients on each contact surface. Conduct thermodynamics simulation calculations on the dynamically calibrated rubber bearing finite element model to verify the accuracy of the model.

[0082] The comparison results of the simulation data and test data of the stress at position C of the rubber bearing are as follows Figure 7 shown. The average value of the simulation data is 227.6 MPa, and the average value of the test data is 249.8 MPa. The degree of coincidence between the two is 91.1%.

[0083] The comparison results of the simulation data and test data of the temperature at position D of the rubber bearing are as follows Figure 8As shown, the average value of the simulation data is 22.5 °C, and the average value of the test data is 23.8 °C. The degree of coincidence between the two is 94.5%.

[0084] Through comparative analysis, it can be obtained that the degree of coincidence between the stress and temperature simulation data and the test data of the thermo-solid coupling model of the rubber bearing after dynamic calibration is greater than 90%, indicating that the established finite element model of the rubber bearing is accurate.

[0085] If the degree of coincidence between the simulation data and the test data does not meet the preset requirements, then re-establish the thermo-solid coupling nonlinear model of the rubber bearing, and repeat the foregoing steps until the model meets the preset requirements.

[0086] (5) Set boundary conditions and load conditions, and perform creep simulation calculation of the rubber bearing

[0087] Apply the time-domain data of Fy and Fz at the center point Nc of the finite element model of the rubber bearing after static calibration; apply the temperature field simulation results of the thermo-solid coupling model of the rubber bearing to the finite element model of the rubber bearing; define the creep model and coefficient of the rubber material, set the tolerance of the creep strain, and perform display simulation solution calculation; draw the displacement change curve of the rubber part, as Figure 9 shown. After 100 s, the creep amount at point C is 2.68×10-6 mm, and the creep amount at point D is 2.41×10-6 mm.

[0088] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, and are not intended to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that they can still modify the technical solutions described in the foregoing embodiments, or perform equivalent replacements for some or all of the technical features; and these modifications or replacements do not cause the essence of the corresponding technical solutions to deviate from the technical solutions of the embodiments of the present invention.

Claims

1. A non-linear finite element modeling method for rubber bearing creep analysis, characterized in that, Including: Establish a non - linear finite element model of the rubber bearing and conduct stiffness and damping simulation calculations, specifically including: establishing a finite element mesh model of the bearing, a finite element mesh model of the rubber part, and a finite element mesh model of the bracket; establishing the contact relationship between the bearing rollers and the inner and outer rings of the bearing and setting the friction coefficient; establishing a fixed constraint between the outer ring of the bearing and the rubber part; establishing the contact relationship between the bracket and the rubber part and setting the friction coefficient; setting the Mooney - Rivlin constitutive model parameters of the rubber part; establishing a node at the center point of the bearing and connecting the node at the center point of the bearing to the nodes on the inner surface of the bearing inner ring using rb2, and establishing fixed constraints at the positions of two bolts on the bracket; applying a forced displacement load on the node at the center point of the bearing to conduct simulation calculations of the Z - direction stiffness and Y - direction stiffness of the rubber bearing; applying an impulse displacement load on the node at the center point of the bearing to conduct simulation calculations of the Z - direction damping and Y - direction damping of the rubber bearing; Conduct stiffness and damping tests on the rubber bearing and conduct static calibration analysis on the non - linear finite element model of the rubber bearing according to the test results; When the static calibration analysis meets the first preset requirement, based on the non - linear finite element model of the rubber bearing that meets the preset requirements, establish a thermo - solid coupling non - linear model of the rubber bearing and conduct steady - state and transient thermodynamics simulation analysis; Conduct stress and temperature tests on the rubber bearing and conduct dynamic calibration analysis on the non - linear finite element model of the rubber bearing according to the test results; When the dynamic calibration analysis meets the second preset requirement, set boundary conditions and load conditions and conduct creep simulation calculations of the rubber bearing; 2. The non-linear finite element modeling method for rubber bearing creep analysis according to claim 1, wherein Conduct stiffness and damping tests on the rubber bearing and conduct static calibration analysis on the non - linear finite element model of the rubber bearing according to the test results, further including: Testing the Y - direction and Z - direction stiffness, Y - direction and Z - direction damping of the rubber bearing; Establishing a parametric model of the Mooney - Rivlin constitutive model and setting c10, c01, and c11 as design variables; Setting the Y - direction and Z - direction stiffness test data and Y - direction and Z - direction damping test data of the rubber bearing as the objective function; Determining the values of c10, c01, and c11 through optimization simulation calculations; Modifying the Mooney - Rivlin constitutive model of the rubber part according to the optimization results of c10, c01, and c11; Conducting static - mechanics simulation calculations on the static - calibrated non - linear finite element model of the rubber bearing to verify the accuracy of the model.

3. The non-linear finite element modeling method for rubber bearing creep analysis according to claim 2, wherein Establishing a thermo - solid coupling non - linear model of the rubber bearing and conducting steady - state and transient thermodynamics simulation analysis, further including: Setting the element type of the non - linear finite element model of the rubber bearing to solid70; Setting the heat generation rate and heat conduction coefficient on the contact surfaces between the bearing rollers and the inner and outer rings of the bearing, setting the heat conduction coefficient on the connection surface between the outer ring of the bearing and the rubber part, and setting the heat conduction coefficient on the contact surface between the bracket and the rubber part; Applying time - varying forces Fy and Fz at the position of the node at the center point of the bearing; Conducting simulation calculations of frictional heat generation and heat transfer.

4. The non-linear finite element modeling method for rubber bearing creep analysis according to claim 3, wherein Conducting stress and temperature tests on the rubber bearing and conduct dynamic calibration analysis on the non - linear finite element model of the rubber bearing according to the test results, further including: Install strain gauges on the rubber parts to measure the stress test data of the rubber parts; Arrange temperature sensors on the rubber parts to measure the temperature test data of the rubber parts; Set the friction coefficients on the contact surfaces between the bearing rollers and the inner and outer rings of the bearing, and between the bracket and the rubber parts as design variables; Set the stress test data and temperature test data of the rubber bearing as the objective function; Determine the values of the friction coefficients on each contact surface through optimized simulation calculations; Modify the thermo-solid coupling nonlinear model of the rubber bearing according to the optimization results of the friction coefficients on each contact surface; Conduct thermodynamic simulation calculations on the nonlinear finite element model of the rubber bearing after dynamic calibration to verify the accuracy of the model.

5. The non-linear finite element modeling method for rubber bearing creep analysis according to claim 4, characterized in that Set boundary conditions and load conditions, and conduct creep simulation calculations on the rubber bearing, further including: Apply the time-domain data of Fy and Fz at the center point node position of the nonlinear finite element model of the rubber bearing after static calibration; Apply the temperature field simulation results of the thermo-solid coupling model of the rubber bearing to the nonlinear finite element model of the rubber bearing; Define the creep model and coefficients of the rubber material, set the tolerance of the creep strain, and conduct display simulation solution calculations; Plot the displacement change curve of a certain point on the rubber part to check the creep amount within the specified time.

6. The non-linear finite element modeling method for rubber bearing creep analysis according to claim 1, characterized in that, When the static calibration analysis does not meet the first preset requirement, reset the parameters of the Mooney-Rivlin constitutive model of the rubber part and re-conduct the stiffness and damping simulation calculations of the rubber bearing until the model meets the preset requirement after static calibration.

7. The non-linear finite element modeling method for rubber bearing creep analysis according to claim 2, wherein When the dynamic calibration analysis does not meet the second preset requirement, re-establish the thermo-solid coupling nonlinear model of the rubber bearing and conduct steady-state and transient thermodynamic simulation analyses until the model meets the preset requirement after dynamic calibration.

8. The non-linear finite element modeling method for rubber bearing creep analysis according to any one of claims 6 and 7, characterized in that The first preset requirement and the second preset requirement refer to the requirement for the degree of coincidence between the model simulation results and the test results.

Images