A contact surface dynamic mechanical analysis method, apparatus, medium, and product

Abstract

This application relates to the field of dynamic mechanical analysis technology, providing a method, device, medium, and product for dynamic mechanical analysis of contact surfaces. The method includes: acquiring contact surface state information; acquiring sliding work state variables based on the contact surface state information and tangential stress information; acquiring dynamic shear parameter information based on the sliding work state variables; acquiring contact surface stage information based on the contact surface state information and dynamic shear parameter information; and updating tangential stress information based on the contact surface stage information, contact surface state information, and dynamic shear parameter information. By determining the sliding work state variables, dynamic shear parameter information, contact surface stage information, and updating tangential stress information, a closed-loop feedback mechanism of "energy dissipation - parameter degradation - stress update" is achieved for the contact surface under cyclic loading, improving the analytical accuracy of the dynamic mechanical response of the contact surface and meeting the need for precise analysis of structural mechanical stability in dynamic engineering scenarios.

Description

Technical Field

[0001] This application relates to the field of dynamic mechanical analysis technology, and in particular to a method, equipment, medium and product for dynamic mechanical analysis of contact surfaces. Background Technology

[0002] Contact surfaces or discontinuous interfaces are ubiquitous between various structures and materials, such as friction interfaces in mechanical connectors, construction joints in concrete structures, interlaminar bonding surfaces in composite materials, and faults and joints in rock masses. The mechanical properties of these interfaces are crucial to the stability and reliability of the entire structure. Under cyclic loading, contact surfaces may repeatedly open, slip, and close, leading to a continuous deterioration of shear parameters such as friction coefficient and cohesion, which may result in structural or component failure, damage accumulation, or even catastrophic damage. Therefore, accurately simulating and analyzing the mechanical response of contact surfaces under dynamic cyclic loading is an important prerequisite for conducting safety assessments and optimizing the design of engineering systems.

[0003] However, existing contact surface analysis models mainly focus on the mechanical response under static or monotonic loads. The simplified assumptions make it difficult for the models to accurately characterize the stress changes of the structure under cyclic loads, and cannot meet the needs of accurate analysis of structural mechanical stability in dynamic engineering scenarios.

[0004] Therefore, how to accurately analyze the dynamic mechanical response of the contact surface under cyclic loading is an urgent problem to be solved. Summary of the Invention

[0005] This application provides a method, equipment, medium, and product for dynamic mechanical analysis of contact surfaces, in order to solve the shortcomings of existing technologies that cannot accurately analyze the dynamic mechanical response of contact surfaces under cyclic loading.

[0006] This application provides a method for dynamic mechanical analysis of contact surfaces, including: Obtain contact surface state information, which includes the tangential displacement of the contact surface; Based on the contact surface state information and tangential stress information, the sliding work state variable is obtained, which characterizes the energy dissipation accumulated during sliding of the contact surface. Based on the sliding work state variable, obtain the dynamic shear parameter information; Based on the contact surface state information and the dynamic shear parameter information, the contact surface stage information is obtained, which indicates whether the contact surface is in the elastic stage or the plastic stage. The tangential stress information is updated based on the contact surface stage information, the contact surface state information, and the dynamic shear parameter information; The tangential stress information represents the stress condition of the contact surface.

[0007] According to the contact surface dynamic mechanical analysis method provided in this application, the step of obtaining the sliding work state variable based on the contact surface state information and tangential stress information includes: Calculate the energy dissipation value at the current time step based on the contact surface state information and the tangential stress information; Based on the sliding work state variable of the previous time step and the energy dissipation value of the current time step, obtain the sliding work state variable of the current time step; The initial value of the sliding work state variable is zero.

[0008] According to the contact surface dynamic mechanical analysis method provided in this application, the expression for obtaining the sliding work state variable is as follows: in, For the current time step The sliding work state variable; For the previous time step The sliding work state variable; For time step; The tangential stress at the contact surface; This refers to the tangential plastic slip. The critical energy threshold required for significant degradation of the contact surface; This represents the energy dissipation value at the current time step.

[0009] According to the dynamic mechanical analysis method for contact surfaces provided in this application, the dynamic shear parameter information includes the dynamic friction coefficient and dynamic cohesion; obtaining the dynamic shear parameter information based on the sliding work state variable includes: Based on the sliding work state variable and the friction coefficient decay law, the dynamic friction coefficient at the current time step is calculated and obtained. Based on the sliding work state variable and the cohesion softening law, the dynamic cohesion at the current time step is calculated and obtained.

[0010] According to the dynamic mechanical analysis method for contact surfaces provided in this application, the expression for obtaining the dynamic friction coefficient is as follows: in, The dynamic friction coefficient is mentioned above; The peak friction coefficient; The residual friction coefficient; The attenuation coefficient ( ); The sliding work state variable; And / or, the expression for obtaining the dynamic cohesion is: in, This refers to the dynamic cohesion; For initial cohesion, The critical state quantity for the complete loss of cohesion. Let be the sliding work state variable.

[0011] According to the dynamic mechanical analysis method for contact surfaces provided in this application, the step of obtaining contact surface stage information based on the contact surface state information and the dynamic shear parameter information includes: Based on the contact surface state information, determine the total tangential displacement increment and normal stress at the current time step; Based on the total tangential displacement increment and the tangential stress information of the previous time step, the predicted tangential stress is obtained; Based on the predicted tangential stress, the dynamic shear parameter information, and the normal stress, the contact surface stage information is obtained by calculation using the yield function.

[0012] According to the dynamic mechanical analysis method for contact surfaces provided in this application, the expression for obtaining the predicted tangential stress is as follows: in, For the predicted tangential stress; This refers to the tangential stress information from the previous time step; The tangential elastic stiffness of the contact surface; This is the total tangential displacement increment; And / or, the expression for the yield function is: in, For the predicted tangential stress; The dynamic friction coefficient is mentioned above; This refers to the dynamic cohesion; The normal stress; like If the value is less than or equal to zero, the elastic stage is used as the contact surface stage information; like If the value is greater than zero, the plastic stage is used as the contact surface stage information.

[0013] According to the dynamic mechanical analysis method for contact surfaces provided in this application, the step of updating the tangential stress information based on the contact surface stage information, the contact surface state information, and the dynamic shear parameter information includes: If the contact surface stage information is an elastic stage, then: The tangential stress information is updated based on the predicted tangential stress; If the contact surface stage information is the plastic stage, then: Based on the predicted tangential stress, the normal stress, and the dynamic shear parameters from the previous time step, the tangential plastic slip increment is obtained. The tangential stress information is updated based on the predicted tangential stress, the normal stress, and the dynamic shear parameters from the previous time step. The sliding work state variable at the current time step is corrected based on the tangential plastic slip increment, the updated tangential stress information, and the sliding work state variable from the previous time step.

[0014] According to the dynamic mechanical analysis method for contact surfaces provided in this application, the expression for obtaining the tangential plastic slip increment is as follows: in, This is the plastic slip increment; For the predicted tangential stress; The dynamic friction coefficient mentioned in the previous time step; This refers to the dynamic cohesion described in the previous time step; The normal stress; The tangential elastic stiffness of the contact surface; And / or, the expression for updating the tangential stress information based on the predicted tangential stress, the normal stress, and the dynamic shear parameter information of the previous time step is: in, The updated tangential stress information; For the predicted tangential stress; The dynamic friction coefficient mentioned in the previous time step; This refers to the dynamic cohesion described in the previous time step; The normal stress; It is a symbolic function; And / or, the expression for the sliding work state variable at the current time step that is modified is: in, The corrected sliding work state variable; The sliding work state variable from the previous time step; The updated tangential stress information; This is the plastic slip increment; The critical energy threshold required for significant degradation of the contact surface.

[0015] This application also provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement a contact surface dynamic mechanical analysis method as described above.

[0016] This application also provides a non-transitory computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements a dynamic mechanical analysis method for contact surfaces as described above.

[0017] This application also provides a computer program product, including a computer program that, when executed by a processor, implements a contact surface dynamic mechanical analysis method as described above.

[0018] This application provides a method, device, medium, and product for dynamic mechanical analysis of contact surfaces, which has at least the following beneficial effects: By acquiring contact surface state information and combining it with tangential stress information, a sliding work state variable is obtained, which characterizes the energy dissipation accumulated during contact surface sliding under cyclic loading. Based on the sliding work state variable, dynamic shear parameter information is obtained, enabling the shear parameters to be adjusted in real time according to the actual accumulated damage of the contact surface. Since the mechanical response of the contact surface differs between the elastic and plastic stages, contact surface stage information is obtained based on the contact surface state information and dynamic shear parameter information, clearly determining whether the contact surface is in the elastic or plastic stage. Finally, based on the contact surface stage information, an update strategy corresponding to the elastic and plastic stages is adopted, and the tangential stress information is updated by combining the contact surface state information and dynamic shear parameter information. In this way, by determining the sliding work state variable, dynamic shear parameter information, contact surface stage information, and updating tangential stress information, a closed-loop feedback mechanism of "energy dissipation-parameter degradation-stress update" is realized for the contact surface under cyclic loading. This significantly improves the analytical accuracy of the dynamic mechanical response of the contact surface under cyclic loading, which is beneficial to meeting the need for accurate analysis of structural mechanical stability in dynamic engineering scenarios. Attached Figure Description

[0019] To more clearly illustrate the technical solutions in this application or related technologies, the drawings used in the description of the embodiments or related technologies will be briefly introduced below. Obviously, the drawings described below are some embodiments of this application. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0020] Figure 1 This is one of the flowcharts of a dynamic mechanical analysis method for contact surfaces provided in this application.

[0021] Figure 2 This is the second flowchart of a dynamic mechanical analysis method for contact surfaces provided in this application.

[0022] Figure 3 This is a schematic diagram of the structure of the electronic device provided in this application. Detailed Implementation

[0023] To make the objectives, technical solutions, and advantages of this application clearer, the technical solutions of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some, not all, of the embodiments of this application. Based on the embodiments of this application, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this application.

[0024] Under real-world working conditions, such as cyclic loading scenarios involving earthquakes, wind loads, mechanical vibrations, vehicle traffic, or fatigue loads, the contact surfaces in a structure undergo repeated opening, sliding, and closing behaviors. This is a nonlinear dynamic mechanical response process, which typically leads to a continuous deterioration of shear parameters such as the friction coefficient and cohesion at the contact surfaces, potentially causing structural or component failure. Accurate simulation and analysis of the mechanical response of contact surfaces under dynamic cyclic loading can provide a foundation for engineering safety assessments.

[0025] However, in the dynamic mechanical analysis of the contact surface, it is usually assumed that the shear parameters are constant, making it difficult to accurately characterize the coupled process of shear parameter degradation under cyclic loading, which in turn affects the stress response, and the stress inversely affecting shear parameter degradation. The lack of consideration for this coupled process directly reduces the accuracy of the dynamic mechanical analysis, leading to insufficient precision.

[0026] To solve the above problems, the following will be combined with Figure 1 This application describes a method for dynamic mechanical analysis of contact surfaces, including: S100: Obtain contact surface state information, the contact surface state information including the tangential displacement of the contact surface; S200: Based on the contact surface state information and tangential stress information, obtain the sliding work state variable, which represents the energy dissipation accumulated during the sliding of the contact surface; S300: Obtain dynamic shear parameter information based on the sliding work state variable; S400: Based on the contact surface state information and the dynamic shear parameter information, obtain the contact surface stage information, wherein the contact surface stage information indicates whether the contact surface is in the elastic stage or the plastic stage; S500: Update the tangential stress information based on the contact surface stage information, the contact surface state information, and the dynamic shear parameter information; The tangential stress information represents the stress condition of the contact surface.

[0027] By acquiring contact surface state information and combining it with tangential stress information, a sliding work state variable is obtained, which characterizes the energy dissipation accumulated during contact surface sliding under cyclic loading. Based on the sliding work state variable, dynamic shear parameter information is obtained, enabling the shear parameters to be adjusted in real time according to the actual accumulated damage of the contact surface. Since the mechanical response of the contact surface differs between the elastic and plastic stages, contact surface stage information is obtained based on the contact surface state information and dynamic shear parameter information to clearly determine whether the contact surface is in the elastic or plastic stage. Finally, based on the contact surface stage information, update strategies corresponding to the elastic and plastic stages are adopted, and the tangential stress information is updated by combining the contact surface state information and dynamic shear parameter information.

[0028] In this way, by determining the sliding work state variable, dynamic shear parameter information, contact surface stage information, and updating tangential stress information, a closed-loop feedback mechanism of "energy dissipation-parameter degradation-stress update" is realized for the contact surface under cyclic loading. This significantly improves the analytical accuracy of the dynamic mechanical response of the contact surface under cyclic loading, which is beneficial to meeting the need for accurate analysis of structural mechanical stability in dynamic engineering scenarios.

[0029] It should be noted that the deterioration of shear parameters is related to energy dissipation under cyclic loading. Therefore, by constructing a sliding work state variable to reflect the energy dissipation accumulated by the sliding of the contact surface during the cyclic loading process, and then determining the dynamic shear parameter information based on the sliding work state variable, the dynamic changes of shear parameters under cyclic loading can be reflected, which is beneficial to improving the accuracy of dynamic mechanical analysis.

[0030] Understandably, updating tangential stress information is a closed-loop feedback process. After updating the tangential stress information, a new sliding work state variable is determined based on the contact surface state information and the updated tangential stress information. This process then repeats again, continuously cycling through the "energy dissipation - parameter degradation - stress update" cycle. This clearly reflects stress changes and enables dynamic mechanical analysis of the contact surface. The tangential stress information can be updated according to a preset update time step; that is, the tangential stress information is updated every preset time interval.

[0031] It should be further explained that the contact surface dynamic mechanical analysis method provided in this application can be integrated into the contact model, that is, implemented by the contact model. The contact model can simulate the response changes of the contact surface in the structure under different conditions.

[0032] In some embodiments of this application, the contact surface state information includes the tangential displacement of the contact surface, and may also include information such as the normal stress of the contact surface. The contact surface state information can be obtained by actual detection, simulation, extrapolation based on historical data, etc.

[0033] In some embodiments of the contact surface dynamic mechanical analysis method of this application, step S200 includes: Calculate the energy dissipation value at the current time step based on the contact surface state information and the tangential stress information; Based on the sliding work state variable of the previous time step and the energy dissipation value of the current time step, obtain the sliding work state variable of the current time step; The initial value of the sliding work state variable is zero.

[0034] Based on the contact surface state information and the tangential stress information, the energy dissipation value at the current time step is calculated. Combined with the sliding work state variable from the previous time step, the sliding work state variable at the current time step is obtained to dynamically reflect the energy dissipation accumulated during contact surface sliding. The initial value of the sliding work state variable is zero. By superimposing the energy dissipation value at the current time step with the accumulated value from the previous time step, the cumulative quantification of contact surface sliding energy dissipation is achieved as time steps progress. This allows the state variable to accurately characterize the sliding work accumulated on the contact surface under cyclic loading, thereby reflecting the degree of deterioration of shear parameters and ensuring the accuracy of subsequent acquisition of dynamic friction coefficient and dynamic cohesion.

[0035] Compared to assessing the degree of shear parameter degradation by simply calculating the number of cycles and sliding distance, this application calculates and accumulates the energy dissipation at each time step, thus providing a more solid physical basis for assessing the degree of shear parameter degradation and more accurately reflecting the degradation process of the actual structure.

[0036] In some embodiments of the contact surface dynamic mechanical analysis method of this application, the expression for obtaining the sliding work state variable is: in, For the current time step The sliding work state variable; For the previous time step The sliding work state variable; For time step; The tangential stress at the contact surface; This refers to the tangential plastic slip. The critical energy threshold required for significant degradation of the contact surface; This represents the energy dissipation value at the current time step.

[0037] Current time step Sliding work state variable It is based on the previous time step Sliding work state variable Add the current time step Internal tangential stress With tangential plastic slip The integral (i.e., the energy dissipation value) ) and critical energy threshold The ratio is determined by introducing a critical energy threshold. As a normalization factor, the sliding work state variable It not only quantifies the accumulated energy dissipation, but also correlates it with the critical deterioration criteria, providing dimensionless parameters with physical meaning for the evolution of subsequent shear parameters, which is conducive to obtaining the dynamic friction coefficient and dynamic cohesion more accurately.

[0038] It should be noted that the critical energy threshold The critical energy threshold can be determined through structural experiments. The introduction of this method can transform pure energy dissipation values ​​into indicators representing the degree of degradation.

[0039] In some embodiments of the dynamic mechanical analysis method for contact surfaces in this application, the dynamic shear parameter information includes the dynamic friction coefficient and dynamic cohesion; S300 includes: Based on the sliding work state variable and the friction coefficient decay law, the dynamic friction coefficient at the current time step is calculated and obtained. Based on the sliding work state variable and the cohesion softening law, the dynamic cohesion at the current time step is calculated and obtained.

[0040] Dynamic shear parameters include the dynamic friction coefficient and dynamic cohesion. Based on the sliding work state variable, and according to the friction coefficient decay law and the cohesion softening law, the dynamic friction coefficient and dynamic cohesion at the current time step are calculated. The friction coefficient and cohesion are adjusted in real time and adaptively based on the energy dissipation accumulated at the contact surface, thus reflecting the impact of contact surface deterioration on shear parameters.

[0041] Compared to assuming that the shear parameters are fixed as constants or only change with time, this application directly links the dynamic friction coefficient and dynamic cohesion to the sliding work state variable with clear physical meaning. This can reflect the real deterioration process of the shear parameters of the contact surface, and thus make the changes in the dynamic friction coefficient and dynamic cohesion more consistent with the actual mechanical response.

[0042] In some embodiments of this application, the friction coefficient decay law and the cohesion softening law can adopt different functional forms, which can be determined according to specific structure, materials and other factors.

[0043] In some embodiments of the contact surface dynamic mechanical analysis method of this application, the expression for obtaining the dynamic friction coefficient is as follows: in, The dynamic friction coefficient is mentioned above; The peak friction coefficient; The residual friction coefficient; The attenuation coefficient ( ); Let be the sliding work state variable.

[0044] Dynamic friction coefficient The friction coefficient is determined based on the exponential decay function that reflects the decay law of the friction coefficient. The friction coefficient starts from the initial peak friction coefficient. Gradually decays to the residual friction coefficient The attenuation rate is determined by the attenuation coefficient. Control, the entire process is controlled by the sliding work state variable. drive.

[0045] Peak friction coefficient With residual friction coefficient It can be determined by testing the actual structure.

[0046] In some embodiments of the dynamic mechanical analysis method for contact surfaces in this application, the expression for obtaining the dynamic cohesion is: in, This refers to the dynamic cohesion; For initial cohesion, The critical state quantity for the complete loss of cohesion. Let be the sliding work state variable.

[0047] Dynamic cohesion The cohesion is determined by the linear softening function of the reaction cohesion softening path, starting from the initial cohesion. With the sliding work state variable It decreases linearly with the increase, until... Critical state quantity Completely lost at that time.

[0048] Initial cohesion With critical state quantity It can be determined by testing the actual structure.

[0049] The dynamic friction coefficient and dynamic cohesion are determined by the above calculations, which enables an accurate quantitative description of the nonlinear deterioration of the shear parameters of the contact surface under dynamic loads. This allows for the simulation of the dynamic changes in shear parameters of the contact surface from undamaged to fully worn.

[0050] Besides the exponential decay function, the friction coefficient decay law and the cohesion softening law, besides the linear softening function, can be described using other mathematical functions to depict the dynamic changes of friction coefficient and cohesion, depending on the actual application scenario.

[0051] In some embodiments of the contact surface dynamic mechanical analysis method of this application, step S400 includes: Based on the contact surface state information, determine the total tangential displacement increment and normal stress at the current time step; Based on the total tangential displacement increment and the tangential stress information of the previous time step, the predicted tangential stress is obtained; Based on the predicted tangential stress, the dynamic shear parameter information, and the normal stress, the contact surface stage information is obtained by calculation using the yield function.

[0052] First, based on the contact surface state information, the total tangential displacement increment and normal stress at the current time step are determined. Then, based on the total tangential displacement increment and the tangential stress information from the previous time step, assuming no plastic deformation in the current time step, the predicted tangential stress is obtained. Based on the predicted tangential stress, dynamic shear parameters, and normal stress, the contact surface stage information is obtained through yield function calculation to determine whether the contact surface is in the elastic or plastic stage.

[0053] Therefore, by comparing the stress obtained based on elastic prediction with the dynamic shear parameters that are currently changing, it is possible to accurately determine whether plastic sliding has occurred at the contact surface, and thus determine the stage state of the contact surface at the current time step, providing a basis for subsequent stress updates.

[0054] It should be noted that the yield function can determine when a material or structure transitions from an elastic state to a plastic state. Assuming no plastic deformation occurs in the current time step (i.e., assuming elastic deformation), the predicted tangential stress, after being substituted into the yield function, will be used as updated tangential stress information if the material is determined to be in an elastic state. However, if it is determined to be in a plastic state, the predicted tangential stress is inaccurate and requires correction. Therefore, after obtaining the contact surface stage information, an update strategy corresponding to the stage state is adopted based on whether the contact surface is in an elastic or plastic stage at the current time step to ensure the accuracy of the updated tangential stress information, thereby ensuring the accuracy of the dynamic mechanical analysis.

[0055] In some embodiments of this application, the yield function corresponding to the Mohr-Coulomb yield criterion can be used, and its specific expression is as follows: in, For the yield function value, when At that time, the contact surface is in an elastic state, when When this occurs, the contact surface slides or cracks, indicating that it is in a plastic state; It is the tangential stress; The coefficient of friction; The normal stress at the contact surface; For cohesion.

[0056] It should be emphasized that, unlike the method of using a fixed friction coefficient and cohesion, this application calculates the yield by substituting dynamic shear parameters, namely dynamic friction coefficient and dynamic cohesion, into the above yield function. This makes the yield also dynamically evolving, which can more realistically reflect the process by which the contact surface gradually reduces its shear capacity due to damage accumulation under cyclic loading, and thus more easily enters the plastic state.

[0057] In addition to the functional form based on the Mohr-Coulomb yield criterion, the yield function can also adopt the functional form of the Drucker-Prager criterion or other plastic yield criteria when considering different stress states or material properties.

[0058] In some embodiments of the contact surface dynamic mechanical analysis method of this application, the expression for obtaining the predicted tangential stress is: in, For the predicted tangential stress; This refers to the tangential stress information from the previous time step; The tangential elastic stiffness of the contact surface; This represents the total tangential displacement increment.

[0059] Predicting tangential stress This prediction is based on the assumption that the current time step is in the elastic stage, i.e., without plastic deformation. Specifically, it involves combining the tangential stress information from the previous time step with the tangential elastic stiffness from the contact surface. and total tangential displacement increment The determined elastic stress increments are added together to predict the tangential stress state at the current time step, i.e., to predict the tangential stress.

[0060] In some embodiments of this application, the tangential elastic stiffness of the contact surface The tangential elastic stiffness of the contact surface can be determined through experiments on actual structures. It can be a constant, or a variable related to normal stress, contact area, etc.

[0061] In some embodiments of the contact surface dynamic mechanical analysis method of this application, the expression of the yield function is: in, For the predicted tangential stress; The dynamic friction coefficient is mentioned above; This refers to the dynamic cohesion; The normal stress; like If the value is less than or equal to zero, the elastic stage is used as the contact surface stage information; like If the value is greater than zero, the plastic stage is used as the contact surface stage information.

[0062] Yield function Based on the Mohr-Coulomb criterion, the predicted tangential stress will be... With the dynamic friction coefficient Normal stress and dynamic cohesion The shear strength is obtained through comparison and calculation. The yield function value is... When the yield function value is less than or equal to zero, it is in the elastic stage. A value greater than zero indicates the plastic stage. This provides a clear and quantitative standard for determining the stage of the contact surface, ensuring the accuracy of judging whether the contact surface is in the elastic or plastic stage during dynamic mechanical analysis.

[0063] refer to Figure 2 In some embodiments of the contact surface dynamic mechanical analysis method of this application, step S500 includes: If the contact surface stage information is an elastic stage, then: S510: Update the tangential stress information based on the predicted tangential stress; If the contact surface stage information is the plastic stage, then: S521: Based on the predicted tangential stress, the normal stress, and the dynamic shear parameter information of the previous time step, obtain the tangential plastic slip increment; S522: Update the tangential stress information based on the predicted tangential stress, the normal stress, and the dynamic shear parameter information of the previous time step; S523: Based on the tangential plastic slip increment, the updated tangential stress information, and the sliding work state variable of the previous time step, correct the sliding work state variable of the current time step.

[0064] In this application, the updating of tangential stress information adopts a stress updating algorithm of "elastic prediction-plastic correction", specifically as follows: If the contact surface stage information is in the elastic stage, then the tangential stress information is updated based on the predicted tangential stress calculated under the assumption of no plastic deformation, as expressed by: in, This is the updated tangential stress information; For the predicted tangential stress; This refers to the tangential stress information from the previous time step; The tangential elastic stiffness of the contact surface; This represents the total tangential displacement increment.

[0065] If the contact surface stage information indicates a plastic stage, then plastic correction is performed. The correction steps include: First, based on the predicted tangential stress, normal stress, and dynamic shear parameters from the previous time step, the tangential plastic slip increment is obtained. The plastic slip increment represents the sliding displacement exceeding the elastic limit, and its calculation expression is as follows: in, This is the plastic slip increment; For the predicted tangential stress; The dynamic friction coefficient mentioned in the previous time step; This refers to the dynamic cohesion described in the previous time step; The normal stress; This refers to the tangential elastic stiffness of the contact surface.

[0066] Then, based on the predicted tangential stress, normal stress, and dynamic shear parameters from the previous time step, the tangential stress information is updated to "correct" the stress to the yield surface, ensuring that it meets the yield criterion. The calculation expression is: Right now: in, The updated tangential stress information; For the predicted tangential stress; The dynamic friction coefficient mentioned in the previous time step; This refers to the dynamic cohesion described in the previous time step; The normal stress; It is a symbolic function.

[0067] Updated tangential stress information This ensures that the tangential stress equals the current dynamic shear strength, and is achieved through the sign function. This ensures that the direction of the tangential stress is consistent with the direction of the predicted stress, thereby strictly satisfying the yield criterion.

[0068] Finally, based on the tangential plastic slip increment, the updated tangential stress information, and the sliding work state variable from the previous time step, the sliding work state variable for the current time step is corrected to ensure its accuracy. The calculation expression is as follows: in, The corrected sliding work state variable; The sliding work state variable from the previous time step; The updated tangential stress information; This is the plastic slip increment; The critical energy threshold required for significant degradation of the contact surface.

[0069] Through the aforementioned stress update algorithm of "elastic prediction-plastic correction," the tangential stress information can be accurately updated through corresponding update steps, regardless of whether the contact surface is in the elastic or plastic stage. Simultaneously, this update mechanism ensures that the stress state remains on the yield surface during plastic slip. By correcting the sliding work state variable, the energy dissipated by plasticity is fed back into damage accumulation, realizing the update transfer of "stress update-energy dissipation." This, in turn, achieves a closed-loop feedback mechanism of "energy dissipation-parameter degradation-stress update" for the contact surface under cyclic loading.

[0070] It is important to emphasize that this application not only corrects the stress in the plastic stage to satisfy the yield criterion, but more importantly, it corrects the sliding work state variable of the current time step. This means that the energy dissipated by each plastic sliding will be immediately reflected in the damage accumulation of the contact surface, thereby affecting the subsequent dynamic shear parameters and yield behavior.

[0071] It should be further emphasized that the contact surface dynamic mechanical analysis method provided in this application can be applied to contact models for simulating and analyzing contact surface-related structures to simulate the changes of the contact surface under different conditions, especially the dynamic mechanical response of the contact surface under cyclic loading. Based on the closed-loop feedback mechanism of "energy dissipation-parameter degradation-stress update", it can significantly improve the analysis accuracy of the contact surface dynamic mechanical response under cyclic loading, which is conducive to meeting the needs of accurate analysis of structural mechanical stability in dynamic engineering scenarios.

[0072] Figure 3 An example is a schematic diagram of the physical structure of an electronic device, such as... Figure 3As shown, the electronic device may include a processor 810, a communications interface 820, a memory 830, and a communication bus 840. The processor 810, communications interface 820, and memory 830 communicate with each other via the communication bus 840. The processor 810 can call logical instructions from the memory 830 to execute the aforementioned dynamic mechanical analysis method for contact surfaces.

[0073] Furthermore, the logical instructions in the aforementioned memory 830 can be implemented as software functional units and, when sold or used as independent products, can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of this application, in essence, or the part that contributes to the prior art, or a part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of this application. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0074] On the other hand, this application also provides a computer program product, which includes a computer program that can be stored on a non-transitory computer-readable storage medium. When the computer program is executed by a processor, the computer is able to execute a contact surface dynamic mechanical analysis method provided by the above methods.

[0075] In another aspect, this application also provides a non-transitory computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, is implemented to perform a dynamic mechanical analysis method for contact surfaces provided by the methods described above.

[0076] The device embodiments described above are merely illustrative. The units described as separate components may or may not be physically separate. The components shown as units may or may not be physical units; that is, they may be located in one place or distributed across multiple network units. Some or all of the modules can be selected to achieve the purpose of this embodiment according to actual needs. Those skilled in the art can understand and implement this without any creative effort.

[0077] Through the above description of the embodiments, those skilled in the art can clearly understand that each embodiment can be implemented by means of software plus necessary general-purpose hardware platforms, and of course, it can also be implemented by hardware. Based on this understanding, the above technical solutions, in essence or the part that contributes to the prior art, can be embodied in the form of a software product. This computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, magnetic disk, optical disk, etc., and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute the methods described in the various embodiments or some parts of the embodiments.

[0078] In the description of this application, it should be understood that the terms "first" and "second" are used for descriptive purposes only and should not be construed as indicating or implying relative importance or implicitly specifying the number of indicated technical features. Therefore, a feature defined as "first" or "second" may explicitly or implicitly include at least one of that feature. In the description of this application, "multiple" means at least two, such as two, three, etc., unless otherwise explicitly specified.

[0079] In the description of this specification, the references to terms such as "one embodiment," "some embodiments," "example," "specific example," or "some examples," etc., indicate that a specific feature, structure, material, or characteristic described in connection with that embodiment or example is included in at least one embodiment or example of this application. In this specification, the illustrative expressions of the above terms do not necessarily refer to the same embodiment or example. Furthermore, the specific features, structures, materials, or characteristics described may be combined in any suitable manner in one or more embodiments or examples. Moreover, without contradiction, those skilled in the art can combine and integrate the different embodiments or examples described in this specification, as well as the features of different embodiments or examples.

[0080] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of this application, and are not intended to limit them. Although this application has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some of the technical features. Such modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of this application.

Claims

1. A method for dynamic mechanical analysis of contact surfaces, characterized in that, include: Obtain contact surface state information, which includes the tangential displacement of the contact surface; Based on the contact surface state information and tangential stress information, the sliding work state variable is obtained, which characterizes the energy dissipation accumulated during sliding of the contact surface. Based on the sliding work state variable, obtain the dynamic shear parameter information; Based on the contact surface state information and the dynamic shear parameter information, the contact surface stage information is obtained, which indicates whether the contact surface is in the elastic stage or the plastic stage. The tangential stress information is updated based on the contact surface stage information, the contact surface state information, and the dynamic shear parameter information; The tangential stress information represents the stress condition of the contact surface.

2. The dynamic mechanical analysis method for contact surfaces according to claim 1, characterized in that, The step of obtaining the sliding work state variable based on the contact surface state information and tangential stress information includes: Calculate the energy dissipation value at the current time step based on the contact surface state information and the tangential stress information; Based on the sliding work state variable of the previous time step and the energy dissipation value of the current time step, obtain the sliding work state variable of the current time step; The initial value of the sliding work state variable is zero.

3. The dynamic mechanical analysis method for contact surfaces according to claim 2, characterized in that, The expression for obtaining the sliding work state variable is: in, For the current time step The sliding work state variable; For the previous time step The sliding work state variable; For time step; The tangential stress at the contact surface; This refers to the tangential plastic slip. The critical energy threshold required for significant degradation of the contact surface; This represents the energy dissipation value at the current time step.

4. The dynamic mechanical analysis method for contact surfaces according to claim 2, characterized in that, The dynamic shear parameter information includes the dynamic friction coefficient and dynamic cohesion; obtaining the dynamic shear parameter information based on the sliding work state variable includes: Based on the sliding work state variable and the friction coefficient decay law, the dynamic friction coefficient at the current time step is calculated and obtained. Based on the sliding work state variable and the cohesion softening law, the dynamic cohesion at the current time step is calculated and obtained.

5. The dynamic mechanical analysis method for contact surfaces according to claim 4, characterized in that, The expression for obtaining the dynamic friction coefficient is: in, The dynamic friction coefficient is mentioned above; The peak friction coefficient; The residual friction coefficient; The attenuation coefficient ( ); The sliding work state variable; And / or, the expression for obtaining the dynamic cohesion is: in, This refers to the dynamic cohesion; For initial cohesion, The critical state quantity for the complete loss of cohesion. Let be the sliding work state variable.

6. The dynamic mechanical analysis method for contact surfaces according to claim 4, characterized in that, The step of obtaining contact surface stage information based on the contact surface state information and the dynamic shear parameter information includes: Based on the contact surface state information, determine the total tangential displacement increment and normal stress at the current time step; Based on the total tangential displacement increment and the tangential stress information of the previous time step, the predicted tangential stress is obtained; Based on the predicted tangential stress, the dynamic shear parameter information, and the normal stress, the contact surface stage information is obtained by calculation using the yield function.

7. The dynamic mechanical analysis method for contact surfaces according to claim 6, characterized in that, The expression for obtaining the predicted tangential stress is as follows: in, For the predicted tangential stress; This refers to the tangential stress information from the previous time step; The tangential elastic stiffness of the contact surface; This is the total tangential displacement increment; And / or, the expression for the yield function is: in, For the predicted tangential stress; The dynamic friction coefficient is mentioned above; This refers to the dynamic cohesion; The normal stress; like If the value is less than or equal to zero, the elastic stage is used as the contact surface stage information; like If the value is greater than zero, the plastic stage is used as the contact surface stage information.

8. The dynamic mechanical analysis method for contact surfaces according to claim 6, characterized in that, The step of updating the tangential stress information based on the contact surface stage information, the contact surface state information, and the dynamic shear parameter information includes: If the contact surface stage information is an elastic stage, then: The tangential stress information is updated based on the predicted tangential stress; If the contact surface stage information is the plastic stage, then: Based on the predicted tangential stress, the normal stress, and the dynamic shear parameters from the previous time step, the tangential plastic slip increment is obtained. The tangential stress information is updated based on the predicted tangential stress, the normal stress, and the dynamic shear parameters from the previous time step. The sliding work state variable at the current time step is corrected based on the tangential plastic slip increment, the updated tangential stress information, and the sliding work state variable from the previous time step.

9. The dynamic mechanical analysis method for contact surfaces according to claim 8, characterized in that, The expression for obtaining the tangential plastic slip increment is: in, This is the plastic slip increment; For the predicted tangential stress; The dynamic friction coefficient mentioned in the previous time step; This refers to the dynamic cohesion described in the previous time step; The normal stress; The tangential elastic stiffness of the contact surface; And / or, the expression for updating the tangential stress information based on the predicted tangential stress, the normal stress, and the dynamic shear parameter information of the previous time step is: in, The updated tangential stress information; For the predicted tangential stress; The dynamic friction coefficient mentioned in the previous time step; This refers to the dynamic cohesion described in the previous time step; The normal stress; It is a symbolic function; And / or, the expression for the sliding work state variable at the current time step that is modified is: in, The corrected sliding work state variable; The sliding work state variable from the previous time step; The updated tangential stress information; This is the plastic slip increment; The critical energy threshold required for significant degradation of the contact surface.

10. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the computer program, it implements a dynamic mechanical analysis method for contact surfaces as described in any one of claims 1 to 9.

11. A non-transitory computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by the processor, it implements a dynamic mechanical analysis method for contact surfaces as described in any one of claims 1 to 9.

12. A computer program product, comprising a computer program, characterized in that, When the computer program is executed by the processor, it implements a dynamic mechanical analysis method for contact surfaces as described in any one of claims 1 to 9.

Images