CN-120951703-B - Three-dimensional viscoelastic constitutive finite element implementation method and device, electronic equipment and storage medium

CN120951703BCN 120951703 BCN120951703 BCN 120951703BCN-120951703-B

Abstract

The application discloses a three-dimensional viscoelastic constitutive finite element implementation method, a device, electronic equipment and a storage medium, wherein the method and the device are applied to the electronic equipment, and particularly, a plurality of parameters of a high polymer material or a composite material transmitted by a main program are calculated based on the finite element, a plurality of sub-programs constructed in advance are called to calculate the plurality of parameters, and a simulation calculation result of the high polymer material or the composite material is obtained based on a pseudo strain-strain increment matrix. According to the technical scheme, a user can accurately predict the rigidity degradation and life prolonging behavior of the high polymer material or the composite material without repeated experiments, so that the problem that manpower and material resources are required to be consumed in actual engineering is solved.

Inventors

ZHAO YICONG
LIU YAJUAN
LI ZHONGLIN
LIN XINRU

Assignees

苏州数算软云科技有限公司

Dates

Publication Date: 20260512
Application Date: 20251015

Claims (6)

1. The three-dimensional viscoelastic constitutive finite element implementation method is applied to electronic equipment and is characterized by comprising the following steps of: Calculating a plurality of parameters of the main program transmitted into the high polymer material or the composite material based on the finite element; invoking a plurality of pre-built subroutines to calculate the parameters, and obtaining simulation calculation results of the high polymer material or the composite material based on a pseudo strain-strain increment matrix, wherein the parameters comprise time step, stress component and strain increment; The multiple subroutines comprise a pseudo strain increment subroutine, a tangential stiffness and damage variable computation subroutine, a pseudo strain update subroutine and a stress update subroutine, the multiple subroutines constructed in advance are called to calculate the multiple parameters, and a simulation calculation result of the high polymer material or the composite material is obtained based on a pseudo strain-strain increment matrix, and the method comprises the following steps: invoking the pseudo-strain updating subroutine, and calculating the strain increment by adopting implicit integration to obtain the pseudo-strain increment; Invoking the stress updating subroutine, calculating the pseudo strain increment by implicit integration to obtain a stress increment, and calculating based on the strain increment, the pseudo strain increment and the stress increment to obtain a consistent tangential stiffness matrix in a three-dimensional form; Invoking the pseudo strain increment subroutine to realize the calculation of the relation between the pseudo strain increment coefficient and the real strain increment, and obtaining the pseudo strain-strain increment matrix; And calling the tangent stiffness and damage variable computation operator program, completing the stress update in the incremental steps according to the time step and the stress component and in combination with the constructed incremental constitutive equation, and returning the tangent stiffness matrix of the incremental constitutive equation to the finite element computation main program so that the finite element computation main program obtains the stress and the strain increment through the Newton iteration technology.
2. The method of claim 1, wherein the simulation calculation results include stress strain and damage variables.
3. The three-dimensional viscoelastic constitutive finite element implementation device is applied to electronic equipment, and is characterized in that the three-dimensional viscoelastic constitutive finite element implementation device comprises: The data transmission module is configured to calculate a plurality of parameters of the main program transmitted into the high polymer material or the composite material based on the finite element; The subprogram calling module is configured to call a plurality of subprograms constructed in advance to calculate the plurality of parameters, and obtain a simulation calculation result of the high polymer material or the composite material based on a pseudo strain-strain increment matrix, wherein the plurality of parameters comprise time step, stress component and strain increment; the plurality of subroutines includes a pseudo strain increase subroutine, a tangential stiffness and injury variable calculation subroutine, a pseudo strain update subroutine, and a stress update subroutine, and the subroutine call module includes: The first calling unit is configured to call the pseudo-strain updating subprogram, and the calculation of the strain increment is realized by adopting implicit integration to obtain the pseudo-strain increment; The second calling unit is coordinated and called with the stress updating subprogram, calculates the pseudo strain increment by adopting implicit integration to obtain a stress increment, and calculates the stress increment based on the strain increment, the pseudo strain increment and the stress increment to obtain a consistent tangential stiffness matrix in a three-dimensional form; the third calling unit is configured to call the pseudo strain increment subroutine to calculate the relation between the pseudo strain increment coefficient and the real strain increment, so as to obtain the pseudo strain-strain increment matrix; And the fourth calling unit is configured to call the tangential stiffness and damage variable computation operator program, finish the stress updating in the incremental step according to the time step and the stress component and in combination with the constructed incremental constitutive equation, and return the tangential stiffness matrix of the incremental constitutive equation to the finite element computation main program so that the finite element computation main program obtains the stress and the strain increment through the Newton iteration technology.
4. The three-dimensional viscoelastic constitutive finite element implementation of claim 3, wherein the simulation calculation results comprise stress strain and damage variables.
5. An electronic device comprising at least one processor and a memory coupled to the processor, wherein: the memory is used for storing a computer program or instructions; The processor is configured to execute the computer program or instructions to cause the electronic device to implement the three-dimensional viscoelastic constitutive finite element implementation method according to any one of claims 1 to 2.
6. A computer-readable storage medium for use in an electronic device, wherein the storage medium carries one or more computer programs executable by the electronic device to enable the electronic device to implement the three-dimensional viscoelastic constitutive finite element implementation method according to any one of claims 1-2.

Description

Three-dimensional viscoelastic constitutive finite element implementation method and device, electronic equipment and storage medium Technical Field The application relates to the technical field of material simulation calculation, in particular to a three-dimensional viscoelastic constitutive finite element implementation method, a three-dimensional viscoelastic constitutive finite element implementation device, electronic equipment and a storage medium. Background In mechanical analysis of high molecular materials and composite materials, the effects of viscoelastic, damaging, aging time and prestrain often coexist, and the mechanical response of the materials is commonly affected. Existing commercial finite element tools typically provide only simple viscoelastic or damage models, and it is difficult to simulate the co-action under these effect couplings simultaneously. This results in a polymeric material or composite whose stiffness degradation and life-prolonging behavior cannot be accurately predicted when subjected to long-term loading or multiple cyclic loading. The inventor of the present application finds that no simulation method under the coupling influence of various factors exists in the prior art, so that in actual engineering, engineers either ignore performance degradation caused by damage accumulation, or need to perform a large number of experiments and repeated experiments to estimate the service life of materials, and consume a large amount of manpower and material resources. Disclosure of Invention In view of the above, the present application provides a three-dimensional viscoelastic constitutive finite element implementation method, device, electronic equipment and storage medium, which are used for implementing accurate prediction of rigidity degradation and life-prolonging behavior of a polymer material or a composite material, so as to solve the problem that human and material resources are required to be consumed in actual engineering. In order to achieve the above object, the following solutions have been proposed: the three-dimensional viscoelastic constitutive finite element implementation method is applied to electronic equipment and comprises the following steps: Calculating a plurality of parameters of the main program transmitted into the high polymer material or the composite material based on the finite element; And calling a plurality of pre-constructed subroutines to calculate the parameters, and obtaining a simulation calculation result of the high polymer material or the composite material based on the pseudo strain-strain increment matrix. Optionally, the simulation calculation result includes stress strain and damage variable. Optionally, the plurality of parameters includes a time step, a stress component, and a strain delta. Optionally, the plurality of subroutines includes a pseudo-strain increase subroutine, a tangential stiffness and damage variable calculation subroutine, a pseudo-strain update subroutine, and a stress update subroutine. Optionally, the step of calling a plurality of pre-built subroutines to calculate the plurality of parameters and obtaining a simulation calculation result of the polymer material or the composite material based on a pseudo strain-strain increment matrix includes the steps of: invoking the pseudo-strain updating subroutine, and calculating the strain increment by adopting implicit integration to obtain the pseudo-strain increment; Invoking the stress updating subroutine, calculating the pseudo strain increment by implicit integration to obtain a stress increment, and calculating based on the strain increment, the pseudo strain increment and the stress increment to obtain a consistent tangential stiffness matrix in a three-dimensional form; invoking the pseudo strain increment subroutine to realize the calculation of the relation between the pseudo strain increment coefficient and the real strain increment; And calling the tangent stiffness and damage variable computation operator program, completing the stress update in the incremental steps according to the time step and the stress component and in combination with the constructed incremental constitutive equation, and returning the tangent stiffness matrix of the incremental constitutive equation to the finite element computation main program so that the finite element computation main program obtains the stress and the strain increment through the Newton iteration technology. A three-dimensional viscoelastic constitutive finite element implementation device applied to electronic equipment, the three-dimensional viscoelastic constitutive finite element implementation device comprising: The data transmission module is configured to calculate a plurality of parameters of the main program transmitted into the high polymer material or the composite material based on the finite element; and the subprogram calling module is configured to call a plurality of subprograms construct