Three-dimensional viscoelastic constitutive finite element implementation method, device, electronic equipment and storage medium
By employing the three-dimensional viscoelastic constitutive finite element method and utilizing pseudo-strain-strain increment matrices and subroutines, we can accurately predict the stiffness degradation and lifetime extension of polymeric or composite materials. This solves the problem of requiring extensive experiments in existing technologies and improves computational efficiency and accuracy.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- SUZHOU DIGITAL SOFT CLOUD TECH CO LTD
- Filing Date
- 2025-10-15
- Publication Date
- 2026-05-12
AI Technical Summary
Existing commercial finite element tools are unable to simultaneously simulate the stiffness degradation and life extension behavior of polymer or composite materials under long-term loads or multiple cyclic loads, resulting in the need for a large number of experiments and repeated tests in engineering to estimate the material life.
The three-dimensional viscoelastic constitutive finite element method is adopted. By using the pseudo-strain-strain increment matrix and combining pseudo-strain update, tangent stiffness and damage variable calculation subroutines, the simulation calculation of polymer materials or composite materials can be realized, and the stiffness degradation and life extension behavior can be accurately predicted.
It can accurately predict the stiffness degradation and life extension behavior of polymer materials or composite materials without repeated experiments, saving human and material resources.
Smart Images
Figure CN120951703B_ABST
Abstract
Description
Technical Field
[0001] This application relates to the field of materials simulation and calculation technology, and more specifically, to a method, apparatus, electronic device, and storage medium for realizing three-dimensional viscoelastic constitutive finite element method. Background Technology
[0002] In the mechanical analysis of polymer materials and composites, viscoelastic effects, damage effects, aging time, and pre-strain effects often coexist and jointly influence the mechanical response of the material. Existing commercial finite element analysis tools typically only provide simple viscoelastic or damage models, making it difficult to simultaneously simulate the combined effects of these coupled factors. This results in the inability to accurately predict the stiffness degradation and life extension behavior of polymer materials or composites under long-term loading or multiple cyclic loading.
[0003] The inventors of this application have discovered in the prior art that there is currently no simulation method for such coupled effects of multiple factors. As a result, in actual engineering, engineers either have to ignore the performance degradation caused by damage accumulation or have to conduct a large number of experiments and repeated tests to estimate the material life, which requires a lot of human and material resources. Summary of the Invention
[0004] In view of this, this application provides a three-dimensional viscoelastic constitutive finite element method, apparatus, electronic device and storage medium for accurately predicting the stiffness degradation and life extension behavior of polymer materials or composite materials, so as to solve the problem of consuming human and material resources in actual engineering.
[0005] To achieve the above objectives, the following solution is proposed:
[0006] A method for realizing a three-dimensional viscoelastic constitutive finite element model, applied to electronic devices, includes the following steps:
[0007] The main program for finite element analysis inputs multiple parameters of the polymer or composite material.
[0008] Multiple pre-built subroutines are invoked to calculate the multiple parameters, and simulation calculation results for the polymer material or composite material are obtained based on the pseudo-strain-strain increment matrix.
[0009] Optionally, the simulation results include stress, strain, and damage variables.
[0010] Optionally, the plurality of parameters include time step, stress components, and strain increment.
[0011] Optionally, the plurality of subroutines include a pseudo-strain increment subroutine, a tangent stiffness and damage variable calculation subroutine, a pseudo-strain update subroutine, and a stress update subroutine.
[0012] Optionally, the step of calling multiple pre-built subroutines to calculate the multiple parameters and obtaining the simulation calculation results for the polymer material or composite material based on the pseudo-strain-strain increment matrix includes the following steps:
[0013] The pseudo-strain update subroutine is invoked, and implicit integration is used to calculate the strain increment to obtain the pseudo-strain increment;
[0014] The stress update subroutine is called to calculate the pseudo-strain increment using implicit integration to obtain the stress increment. Based on the strain increment, the pseudo-strain increment, and the stress increment, a three-dimensional uniform tangent stiffness matrix is obtained.
[0015] The pseudo-strain increment subroutine is invoked to calculate the relationship between the pseudo-strain increment coefficient and the true strain increment.
[0016] The subroutine for calculating tangent stiffness and damage variables is invoked. Based on the time step, the stress components, and the constructed incremental constitutive equation, the stress is updated within the incremental step. The tangent stiffness matrix of the incremental constitutive equation is returned to the finite element calculation main program so that the finite element calculation main program can obtain the stress and strain increments through Newton's iteration technique.
[0017] A three-dimensional viscoelastic constitutive finite element method (VFE) implementation device, applied to electronic devices, comprising:
[0018] The data input module is configured to input multiple parameters of polymer materials or composite materials based on the finite element calculation main program;
[0019] The subroutine calling module is configured to call multiple pre-built subroutines to calculate the multiple parameters and obtain the simulation calculation results of the polymer material or composite material based on the pseudo-strain-strain increment matrix.
[0020] Optionally, the simulation results include stress, strain, and damage variables; the multiple parameters include time step, stress components, and strain increment; and the multiple subroutines include a pseudo-strain increment subroutine, a tangent stiffness and damage variable calculation subroutine, a pseudo-strain update subroutine, and a stress update subroutine.
[0021] Optionally, the subroutine calling module includes:
[0022] The first calling unit is configured to call the pseudo-strain update subroutine, and use implicit integration to calculate the strain increment to obtain the pseudo-strain increment;
[0023] The second calling unit calls the stress update subroutine, uses implicit integration to calculate the pseudo-strain increment to obtain the stress increment, and calculates a three-dimensional uniform tangent stiffness matrix based on the strain increment, the pseudo-strain increment and the stress increment.
[0024] The third calling unit is configured to call the pseudo-strain increment subroutine to calculate the relationship between the pseudo-strain increment coefficient and the real strain increment, and obtain the pseudo-strain-strain increment matrix.
[0025] The fourth calling unit is configured to call the tangent stiffness and damage variable calculation subroutine, complete the stress update within the incremental step according to the time step, the stress components, and the constructed incremental constitutive equation, and return the tangent stiffness matrix of the incremental constitutive equation to the finite element calculation main program so that the finite element calculation main program can obtain the stress and strain increments through Newton's iteration technique.
[0026] An electronic device includes at least one processor and a memory connected to the processor, wherein:
[0027] The memory is used to store computer programs or instructions;
[0028] The processor is used to execute the computer program or instructions to enable the electronic device to implement the three-dimensional viscoelastic constitutive finite element method as described above.
[0029] A computer-readable storage medium is used in an electronic device, the storage medium carrying one or more computer programs that can be executed by the electronic device, thereby enabling the electronic device to implement the three-dimensional viscoelastic constitutive finite element method as described above.
[0030] As can be seen from the above technical solution, this application discloses a method, apparatus, electronic device, and storage medium for realizing three-dimensional viscoelastic constitutive finite element methods. This method and apparatus are applied to electronic devices, specifically by using multiple parameters of a polymer or composite material input into the finite element calculation main program; calling multiple pre-built subroutines to calculate these parameters; and obtaining simulation results for the polymer or composite material based on the pseudo-strain-strain increment matrix. The technical solution of this application allows users to accurately predict the stiffness degradation and life extension behavior of polymer or composite materials without repeated experiments, thus solving the problem of consuming significant human and material resources in practical engineering. Attached Figure Description
[0031] To more clearly illustrate the technical solutions in the embodiments of this application or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of this application. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.
[0032] Figure 1 This is a flowchart illustrating a three-dimensional viscoelastic constitutive finite element method according to an embodiment of this application;
[0033] Figure 2 This is a block diagram of a three-dimensional viscoelastic constitutive finite element realization device according to an embodiment of this application;
[0034] Figure 3 This is a block diagram of another three-dimensional viscoelastic constitutive finite element implementation method according to an embodiment of this application;
[0035] Figure 4 This is a block diagram of an electronic device according to an embodiment of this application. Detailed Implementation
[0036] The technical solutions of the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this application, and not all embodiments. Based on the embodiments of this application, all other embodiments obtained by those of ordinary skill in the art without creative effort are within the scope of protection of this application.
[0037] Existing finite element simulations, which consider viscoelastic damage evolution, struggle to accurately account for pre-strain and aging time, simplifying the process by avoiding direct calculations and failing to meet the demands of optimized design. Therefore, the technical solution proposed in this application comprehensively considers various influencing factors of materials, resulting in newer and more accurate simulation algorithms and models. Through theoretical calculations and simulation comparisons of damage accumulation caused by multiple factors, the material properties are explored in depth.
[0038] This method couples viscoelastic behavior with the synergistic effects of aging time, pre-strain, and damage growth, employing a novel constitutive relation expression using pseudo-strain as a transition variable. First, considering the total stress, the constitutive model can be written in the following form:
[0039] (1);
[0040] To decouple the viscoelastic history effect from damage, this method uses pseudo-strain as an internal state variable. Pseudo-strain is used to characterize the equivalent deformation of the material if viscoelastic relaxation were not present. Specifically, the instantaneous elastic response and time-dependent response are separated, and the pseudo-strain and damage deformation are defined accordingly.
[0041] Definition of pseudo-strain:
[0042] (2);
[0043] Definition of damage variables:
[0044] (3)
[0045] In order to integrate the above constitutive model into the finite element solution framework, it is necessary to derive its incremental stress-strain relationship, i.e., the uniform tangent stiffness matrix.
[0046] In the derivation, considering that the damage D is also a function of strain, the chain rule is applied to expand the differentials of each term. After simplification, the incremental form of the one-dimensional tangent stiffness matrix constitutive relation is obtained.
[0047] In the one-dimensional case, the tangent stiffness matrix can be written as:
[0048] (4);
[0049] For a given one-dimensional constitutive expression, the parameters a, b, and K are obtained from the following relationship, where the external adjustment is: t a ( / day) and Ep (dimensionless) represent aging time and pre-strain rate, respectively, while other parameters are constants.
[0050]
[0051]
[0052] (5);
[0053] For small deformation problems, we can obtain an incremental form that is independent of k:
[0054] (6);
[0055] Construct the elasticity matrix:
[0056] (7)
[0057] The three-dimensional incremental form of this constitutive model can be written as:
[0058] (8);
[0059] in k 1 ，k 2 =1,2…，n comp Aging time and pre-strain need to be considered in order to calculate various parameters.
[0060] Viscoelastic materials, such as polymers and composites, exhibit a "memory effect" in response to stress. Directly integrating over historical strains results in extremely high computational costs. This application significantly improves numerical computation efficiency by deriving an integral incremental form of the viscoelastic constitutive equation.
[0061] First, by using mathematical transformations such as integration by parts, the pseudo-strain is... Represented as historical physical strain The convolution integral form clearly reflects the dependence of viscoelasticity on past strain.
[0062] Then, the total time interval is divided into several small increment steps by incremental discretization of time. △t Assuming the strain change is approximately linear within each step, the above integral formula is discretized, leading to the following derivation: ;
[0063] in For the strain increment at step k+1, This is the pseudo-strain value from the previous step.
[0064] By using recursive updates, only the pseudo-strain from the previous step needs to be retained. This can be based on the current strain increment. Calculate the pseudo-strain increment This leads to a new pseudo-strain: ;
[0065] In this way, internal variables can be updated efficiently without storing the complete historical strain sequence. Under the assumption of small deformation, it can be further simplified to a unified incremental update expression that is independent of the number of steps, indicating that the pseudo-strain increment of each step is only related to the data of the current and previous incremental steps, and is independent of the specific step number, which greatly simplifies the code implementation.
[0066] When implemented using computer equipment, the new strain increment can be adjusted at each iteration ( Recursive update ( ), thereby obtaining new Finally, the spurious strain is brought back into the damage equation to update the material stiffness coefficient, and thus the stress.
[0067] By employing the aforementioned recursive integral algorithm, the full-history memory effect of viscoelasticity is condensed into iterative updates of a small number of internal variables, which not only ensures simulation accuracy but also significantly improves computational efficiency. Based on the above prerequisites, this application provides the following specific implementation method:
[0068] Figure 1 This is a flowchart illustrating a three-dimensional viscoelastic constitutive finite element method according to an embodiment of this application.
[0069] like Figure 1 As shown, the three-dimensional viscoelastic constitutive finite element method provided in this embodiment is applied to electronic devices, and specifically includes the following steps:
[0070] S1. Input multiple parameters of polymer materials or composite materials into the main program based on finite element calculation.
[0071] The main finite element calculation program can be commercial finite element software, such as ANSYS, and can accept multiple parameters including but not limited to time step, stress components, and strain increments. It can also include multiple user-defined variables, including but not limited to: aging time, pre-strain, Poisson's ratio, relaxation times and relaxation modulus at each stage, pseudo-strain and its increment, and various material damage parameters.
[0072] S2. Call multiple pre-built subroutines to calculate multiple parameters and obtain simulation calculation results for polymer materials or composite materials.
[0073] The simulation results here include, but are not limited to, stress, strain, and damage variables. Several pre-built subroutines include a pseudo-strain increment subroutine, a tangent stiffness and damage variable calculation subroutine, a pseudo-strain update subroutine, and a stress update subroutine. The specific simulation calculations include the following:
[0074] First, the pseudo-strain update subroutine is called, and implicit integration is used to calculate the strain increment to obtain the pseudo-strain increment;
[0075] Then, the stress update subroutine is called to calculate the pseudo-strain increment using implicit integration to obtain the stress increment. Based on the strain increment, pseudo-strain increment, and stress increment, a three-dimensional uniform tangent stiffness matrix is obtained.
[0076] Next, the pseudo-strain increment subroutine is called to calculate the relationship between the pseudo-strain increment coefficient and the real strain increment; that is, by using the pseudo-strain to satisfy the historical integral relationship, and with the pseudo-strain as an intermediate variable, the pseudo-strain-strain increment matrix under the action of historical integral is established.
[0077] Finally, the subroutine for calculating tangent stiffness and damage variables is called. Based on the time step, stress components, and the constructed incremental constitutive equation, the stress is updated within the incremental step. The tangent stiffness matrix of the incremental constitutive equation is returned to the main program so that the main program can obtain the stress and strain increments through Newton's iteration technique.
[0078] As can be seen from the above technical solution, this embodiment provides a three-dimensional viscoelastic constitutive finite element method. This method is applied to electronic devices, specifically by using multiple parameters of polymeric materials or composite materials input into the finite element calculation main program; calling multiple pre-built subroutines to calculate multiple parameters; and obtaining simulation calculation results of polymeric materials or composite materials based on the pseudo-strain-strain increment matrix. The technical solution of this application allows users to accurately predict the stiffness degradation and life extension behavior of polymeric materials or composite materials without repeated experiments, thereby solving the problem of consuming human and material resources in actual engineering.
[0079] The flowcharts and block diagrams in the accompanying drawings illustrate the architecture, functionality, and operation of possible implementations of systems, methods, and computer program products according to various embodiments of this disclosure. In this regard, each block in a flowchart or block diagram may represent a module, segment, or portion of code containing one or more executable instructions for implementing a specified logical function. It should also be noted that in some alternative implementations, the functions indicated in the blocks may occur in a different order than those indicated in the drawings. For example, two consecutively indicated blocks may actually be executed substantially in parallel, and they may sometimes be executed in reverse order, depending on the functions involved. It should also be noted that each block in the block diagrams and / or flowcharts, and combinations of blocks in the block diagrams and / or flowcharts, can be implemented using a dedicated hardware-based system that performs the specified function or operation, or using a combination of dedicated hardware and computer instructions.
[0080] Although the operations are described in a specific order, this should not be construed as requiring these operations to be performed in the specific order shown or in a sequential order. In certain environments, multitasking and parallel processing may be advantageous.
[0081] It should be understood that the steps described in the method embodiments of this disclosure may be performed in different orders and / or in parallel. Furthermore, the method embodiments may include additional steps and / or omit the steps shown. The scope of this disclosure is not limited in this respect.
[0082] Computer program code for performing the operations of this disclosure can be written in one or more programming languages or a combination thereof, including but not limited to object-oriented programming languages such as Java, Smalltalk, and C++, as well as conventional procedural programming languages such as C or similar languages. The program code can be executed entirely on the user's computer, partially on the user's computer, as a standalone software package, partially on the user's computer and partially on a remote computer, or entirely on a remote computer or server. In cases involving remote computers, the remote computer can be connected to the user's computer via any type of network—including a local area network (LAN) or a wide area network (WAN)—or can be connected to an external computer.
[0083] Figure 2 This is a block diagram of a three-dimensional viscoelastic constitutive finite element implementation device according to an embodiment of this application.
[0084] like Figure 2 As shown, the three-dimensional viscoelastic constitutive finite element method provided in this embodiment is applied to electronic devices, specifically including a data input module 10 and a subroutine call module 20.
[0085] The data input module is used to input multiple parameters of polymer materials or composite materials into the main program based on the finite element calculation.
[0086] The main finite element calculation program can be commercial finite element software, such as ANSYS, and can accept multiple parameters including but not limited to time step, stress components, and strain increments. It can also include multiple user-defined variables, including but not limited to: aging time, pre-strain, Poisson's ratio, relaxation times and relaxation modulus at each stage, pseudo-strain and its increment, and various material damage parameters.
[0087] The subroutine calling module is used to call multiple pre-built subroutines to calculate multiple parameters and obtain simulation calculation results for polymer materials or composite materials.
[0088] The simulation results here include, but are not limited to, stress, strain, and damage variables. Several pre-built subroutines include a pseudo-strain increment subroutine, a tangent stiffness and damage variable calculation subroutine, a pseudo-strain update subroutine, and a stress update subroutine. This module specifically includes a first calling unit 21, a second calling unit 22, a third calling unit 23, and a fourth calling unit 24, as follows... Figure 3 As shown.
[0089] The first calling unit is used to call the pseudo-strain update subroutine, which uses implicit integration to calculate the strain increment and obtain the pseudo-strain increment.
[0090] The second calling unit is used to call the stress update subroutine, which uses implicit integration to calculate the pseudo-strain increment to obtain the stress increment, and calculates the three-dimensional uniform tangent stiffness matrix based on the strain increment, pseudo-strain increment and stress increment.
[0091] The third calling unit is used to call the pseudo-strain increment subroutine to calculate the relationship between the pseudo-strain increment coefficient and the real strain increment; that is, by using pseudo-strain to satisfy the historical integral relationship, and with pseudo-strain as an intermediate variable, a pseudo-strain-strain increment matrix under the action of historical integral is established.
[0092] The fourth calling unit is used to call the subroutine for calculating tangent stiffness and damage variables. Based on the time step, stress components, and the constructed incremental constitutive equation, it completes the stress update within the incremental step and returns the tangent stiffness matrix of the incremental constitutive equation to the main program, so that the main program can obtain the stress and strain increments through Newton's iteration technique.
[0093] As can be seen from the above technical solution, this embodiment provides a three-dimensional viscoelastic constitutive finite element realization device. This device is applied to electronic devices, specifically based on multiple parameters of polymeric materials or composite materials input into the finite element calculation main program; it calls multiple pre-built subroutines to calculate multiple parameters, and obtains the simulation calculation results of polymeric materials or composite materials based on the pseudo-strain-strain increment matrix. The technical solution of this application allows users to accurately predict the stiffness degradation and life extension behavior of polymeric materials or composite materials without repeated experiments, thereby solving the problem of consuming human and material resources in actual engineering.
[0094] The units described in the embodiments of this disclosure can be implemented in software or in hardware. The name of a unit does not necessarily limit the unit itself; for example, the first acquisition unit can also be described as "a unit that acquires at least two Internet Protocol addresses".
[0095] The functions described above in this document can be performed at least in part by one or more hardware logic components. For example, exemplary types of hardware logic components that can be used, without limitation, include: field-programmable gate arrays (FPGAs), application-specific integrated circuits (ASICs), application-specific standard products (ASSPs), system-on-a-chip (SoCs), complex programmable logic devices (CPLDs), and so on.
[0096] Figure 4 This is a block diagram of an electronic device according to an embodiment of this application.
[0097] The following is for reference. Figure 4This document illustrates a structural diagram suitable for implementing the electronic device in the embodiments of this disclosure. The terminal device in the embodiments of this disclosure may include, but is not limited to, mobile terminals such as mobile phones, laptops, digital broadcast receivers, PDAs (personal digital assistants), PADs (tablet computers), PMPs (portable multimedia players), in-vehicle terminals (e.g., in-vehicle navigation terminals), and fixed terminals such as digital TVs and desktop computers. This electronic device is merely an example and should not impose any limitations on the functionality and scope of use of the embodiments of this disclosure.
[0098] The electronic device may include a processing unit (e.g., a central processing unit, a graphics processing unit, etc.) 401, which can perform various appropriate actions and processes according to a program stored in a read-only memory (ROM) 402 or a program loaded from an input device 406 into a random access memory (RAM) 403. The RAM also stores various programs and data required for the operation of the electronic device. The processing unit, ROM, and RAM are interconnected via a bus 404. An input / output (I / O) interface 405 is also connected to the bus 404.
[0099] Typically, the following devices can be connected to the I / O interface: input devices including, for example, touchscreens, touchpads, keyboards, mice, cameras, microphones, accelerometers, gyroscopes, etc.; output devices 407 including, for example, liquid crystal displays (LCDs), speakers, vibrators, etc.; storage devices 408 including, for example, magnetic tapes, hard disks, etc.; and communication devices 409. Communication device 409 allows the electronic device to communicate wirelessly or wiredly with other devices to exchange data. Although electronic devices with various devices are shown in the figures, it should be understood that it is not required to implement or possess all of the devices shown. More or fewer devices may be implemented or possessed alternatively.
[0100] This application provides an embodiment of a computer-readable storage medium.
[0101] The aforementioned computer-readable storage medium is applied to an electronic device and carries one or more computer programs. When these programs are executed by the electronic device, the device calculates multiple parameters of the polymer or composite material input into the main program using finite element analysis; it then calls multiple pre-built subroutines to calculate these parameters and obtains simulation results for the polymer or composite material based on the pseudo-strain-strain increment matrix. This technical solution allows users to accurately predict the stiffness degradation and life extension behavior of polymer or composite materials without repeated experiments, thus solving the problem of consuming significant human and material resources in practical engineering.
[0102] It should be noted that the computer-readable medium described above in this disclosure can be a computer-readable signal medium or a computer-readable storage medium, or any combination of the two. A computer-readable storage medium can be, for example,—but not limited to—an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any combination thereof. More specific examples of a computer-readable storage medium may include, but are not limited to: an electrical connection having one or more wires, a portable computer disk, a hard disk, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), optical fiber, portable compact disk read-only memory (CD-ROM), optical storage device, magnetic storage device, or any suitable combination thereof.
[0103] In this disclosure, a computer-readable storage medium can be any tangible medium that contains or stores a program that can be used by or in connection with an instruction execution system, apparatus, or device. In this disclosure, a computer-readable signal medium can include a data signal propagated in baseband or as part of a carrier wave, carrying computer-readable program code. Such propagated data signals can take various forms, including but not limited to electromagnetic signals, optical signals, or any suitable combination thereof. A computer-readable signal medium can also be any computer-readable medium other than a computer-readable storage medium, which can send, propagate, or transmit a program for use by or in connection with an instruction execution system, apparatus, or device. The program code contained on the computer-readable medium can be transmitted using any suitable medium, including but not limited to: wires, optical fibers, RF (radio frequency), etc., or any suitable combination thereof.
[0104] The various embodiments in this specification are described in a progressive manner, with each embodiment focusing on the differences from other embodiments. The same or similar parts between the various embodiments can be referred to each other.
[0105] Although preferred embodiments of the present invention have been described, those skilled in the art, upon learning the basic inventive concept, can make other changes and modifications to these embodiments. Therefore, the appended claims are intended to be interpreted as including the preferred embodiments as well as all changes and modifications falling within the scope of the embodiments of the present invention.
[0106] Finally, it should be noted that in this document, relational terms such as "first" and "second" are used only to distinguish one entity or operation from another, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Furthermore, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or terminal device that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or terminal device. Without further limitations, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, method, article, or terminal device that includes said element.
[0107] The technical solution provided by the present invention has been described in detail above. Specific examples have been used to illustrate the principle and implementation of the present invention. The description of the above embodiments is only for the purpose of helping to understand the method and core idea of the present invention. At the same time, for those skilled in the art, there will be changes in the specific implementation and application scope based on the idea of the present invention. Therefore, the content of this specification should not be construed as a limitation of the present invention.
Claims
1. A three-dimensional viscoelastic constitutive finite element method for use in electronic devices, characterized in that: The method for realizing the three-dimensional viscoelastic constitutive finite element method includes the following steps: The main program for finite element analysis inputs multiple parameters of the polymer or composite material. Multiple pre-built subroutines are invoked to calculate the multiple parameters, and simulation calculation results for the polymer material or composite material are obtained based on the pseudo-strain-strain increment matrix. The multiple parameters include time step, stress components and strain increment. The multiple subroutines include a pseudo-strain increment subroutine, a tangent stiffness and damage variable calculation subroutine, a pseudo-strain update subroutine, and a stress update subroutine. The step of calling these pre-built subroutines to calculate the multiple parameters and obtaining simulation results for the polymer material or composite material based on the pseudo-strain-strain increment matrix includes the following steps: The pseudo-strain update subroutine is invoked, and implicit integration is used to calculate the strain increment to obtain the pseudo-strain increment; The stress update subroutine is called to calculate the pseudo-strain increment using implicit integration to obtain the stress increment. Based on the strain increment, the pseudo-strain increment, and the stress increment, a three-dimensional uniform tangent stiffness matrix is obtained. The pseudo-strain increment subroutine is invoked to calculate the relationship between the pseudo-strain increment coefficient and the real strain increment, thereby obtaining the pseudo-strain-strain increment matrix. The subroutine for calculating tangent stiffness and damage variables is invoked. Based on the time step, the stress components, and the constructed incremental constitutive equation, the stress is updated within the incremental step. The tangent stiffness matrix of the incremental constitutive equation is returned to the finite element calculation main program so that the finite element calculation main program can obtain the stress and strain increments through Newton's iteration technique.
2. The three-dimensional viscoelastic constitutive finite element method as described in claim 1, characterized in that, The simulation results include stress, strain, and damage variables.
3. A three-dimensional viscoelastic constitutive finite element realization device, applied to electronic equipment, characterized in that, The three-dimensional viscoelastic constitutive finite element realization device includes: The data input module is configured to input multiple parameters of polymer materials or composite materials based on the finite element calculation main program; The subroutine calling module is configured to call multiple pre-built subroutines to calculate the multiple parameters and obtain the simulation calculation results of the polymer material or composite material based on the pseudo-strain-strain increment matrix. The multiple parameters include time step, stress components and strain increment. The multiple subroutines include a pseudo-strain increment subroutine, a tangent stiffness and damage variable calculation subroutine, a pseudo-strain update subroutine, and a stress update subroutine. The subroutine calling module includes: The first calling unit is configured to call the pseudo-strain update subroutine, and use implicit integration to calculate the strain increment to obtain the pseudo-strain increment; The second calling unit calls the stress update subroutine, uses implicit integration to calculate the pseudo-strain increment to obtain the stress increment, and calculates a three-dimensional uniform tangent stiffness matrix based on the strain increment, the pseudo-strain increment and the stress increment. The third calling unit is configured to call the pseudo-strain increment subroutine to calculate the relationship between the pseudo-strain increment coefficient and the real strain increment, and obtain the pseudo-strain-strain increment matrix. The fourth calling unit is configured to call the tangent stiffness and damage variable calculation subroutine, complete the stress update within the incremental step according to the time step, the stress components, and the constructed incremental constitutive equation, and return the tangent stiffness matrix of the incremental constitutive equation to the finite element calculation main program so that the finite element calculation main program can obtain the stress and strain increments through Newton's iteration technique.
4. The three-dimensional viscoelastic constitutive finite element realization device as described in claim 3, characterized in that, The simulation results include stress, strain, and damage variables.
5. An electronic device, characterized in that, The electronic device includes at least one processor and a memory connected to the processor, wherein: The memory is used to store computer programs or instructions; The processor is used to execute the computer program or instructions to enable the electronic device to implement the three-dimensional viscoelastic constitutive finite element method as described in any one of claims 1 to 2.
6. A computer-readable storage medium for use in electronic devices, characterized in that, The storage medium carries one or more computer programs that can be executed by the electronic device, thereby enabling the electronic device to implement the three-dimensional viscoelastic constitutive finite element method as described in any one of claims 1 to 2.