A phase field fracture simulation method, device and equipment for viscoplastic material

Abstract

The present invention relates to the field of viscoplastic material stress strain and failure damage simulation, and in particular discloses a phase field fracture simulation method, device, equipment and computer-readable storage medium for viscoplastic materials, which receives initial deformation field data and material loading history data of the tested viscoplastic material; loads the initial deformation field step by step according to the material loading history data to obtain the final field cumulative deformation energy density field data; wherein, after each load step calculation is completed, the corresponding stress strain field and phase field are generated, and after the corresponding cumulative strain energy density field is determined by the stress strain field, the phase field is corrected according to the cumulative strain energy density field to obtain the fracture phase field; and the next load step is calculated according to the material loading history data, the stress strain field and the fracture phase field. The present invention is closer to the actual situation of viscoplastic materials bearing loads and improves simulation accuracy.

Description

Technical Field

[0001] The present invention relates to the field of viscoplastic material stress strain and failure simulation, and in particular to a phase field fracture simulation method, device, equipment and computer-readable storage medium for viscoplastic materials. Background Art

[0002] The service environment of the hot end components of aviation gas turbine engines is very harsh. The mechanical response of the turbine disk and turbine blades under high temperature is very complex, which will produce nonlinear plastic deformation behaviors related to cyclic hardening, softening, creep, relaxation, etc., which will affect the crack initiation and expansion of the turbine disk material. At present, the calculation methods of material failure provided by commercial finite element software (including CZM, XFEM, node release technology, etc.) use relatively simple constitutive models and cannot describe the crack expansion behavior of the hot end component materials of aviation engines under the influence of complex nonlinear deformation behaviors such as cyclic hardening, softening, creep, relaxation, etc. The traditional continuum mechanics method can only perform static analysis on materials containing defects / microcracks, and cannot restore the real dynamic expansion process of cracks.

[0003] Therefore, it is necessary to develop a dynamic simulation method for crack propagation that can take into account the rate-dependent complex plastic deformation behaviors of materials such as cyclic hardening, softening, creep, and relaxation, and embed it into existing commercial finite element software to achieve engineering applications and calculations, in order to make up for the current inability to restore the true dynamic propagation process of cracks in viscoplastic materials and the low simulation accuracy. Summary of the invention

[0004] The purpose of the present invention is to provide a phase field fracture simulation method, device, equipment and computer-readable storage medium for viscoplastic materials to solve the problem that the existing technology cannot restore the real dynamic expansion process of cracks in viscoplastic materials and the simulation accuracy is low.

[0005] In order to solve the above technical problems, the present invention provides a phase field fracture simulation method for viscoplastic materials, comprising:

[0006] Receiving initial deformation field data and material loading history data of the tested viscoplastic material;

[0007] The initial deformation field is loaded step by step according to the material loading history data to obtain the final field cumulative deformation energy density field data; wherein, after each load step calculation is completed, the corresponding stress-strain field and phase field are generated, and after the corresponding cumulative strain energy density field is determined by the stress-strain field, the phase field is corrected according to the cumulative strain energy density field to obtain the fracture phase field; and the next load step is calculated according to the material loading history data, the stress-strain field and the fracture phase field.

[0008] Optionally, in the phase field fracture simulation method for viscoplastic materials, the calculation of each load step is completed to generate a corresponding stress-strain field and a phase field, and determining a corresponding cumulative strain energy density field through the stress-strain field includes:

[0009] Complete the calculation of the current load step and generate the corresponding stress-strain field and phase field;

[0010] Judging whether the stress-strain field and the phase field both meet preset convergence conditions based on the virtual displacement principle;

[0011] When the stress-strain field and the phase field meet a preset convergence condition, a corresponding cumulative strain energy density field is determined according to the stress-strain field.

[0012] Optionally, in the phase field fracture simulation method for viscoplastic materials, the phase field is corrected according to the accumulated strain energy density field to obtain a fracture phase field; and the calculation of the next load step is performed according to the material loading history data, the stress-strain field and the fracture phase field, including:

[0013] Correcting the phase field according to the accumulated strain energy density field to obtain a fracture phase field;

[0014] Judging whether the fracture phase field meets a preset convergence condition based on the virtual displacement principle;

[0015] When the fracture phase field meets the preset convergence condition, the calculation of the next load step is performed according to the material loading history data, the stress-strain field and the fracture phase field.

[0016] Optionally, in the phase field fracture simulation method of viscoplastic materials, after obtaining the final field accumulated deformation energy density field data, the method further includes:

[0017] A visualization result is obtained according to the final-field accumulated deformation energy density field data.

[0018] Optionally, in the phase field fracture simulation method for viscoplastic materials, the initial deformation field data is model data of mesh encryption in a potential crack region.

[0019] A phase field fracture simulation device for viscoplastic materials, comprising:

[0020] A receiving module, used for receiving the initial deformation field data and material loading history data of the tested viscoplastic material;

[0021] A calculation module is used to gradually load the initial deformation field with load steps according to the material loading history data to obtain the final field cumulative deformation energy density field data; wherein, after completing the calculation of each load step, the corresponding stress-strain field and phase field are generated, and after the corresponding cumulative strain energy density field is determined by the stress-strain field, the phase field is corrected according to the cumulative strain energy density field to obtain the fracture phase field; and the next load step is calculated according to the material loading history data, the stress-strain field and the fracture phase field.

[0022] Optionally, in the phase field fracture simulation device for viscoplastic materials, the calculation module includes:

[0023] The stress-strain calculation unit is used to complete the calculation of the current load step and generate the corresponding stress-strain field and phase field;

[0024] A calculation convergence checking unit, used for judging whether the stress-strain field and the phase field both meet preset convergence conditions based on the virtual displacement principle;

[0025] The cumulative strain energy calculation unit is used to determine the corresponding cumulative strain energy density field according to the stress strain field when the stress strain field and the phase field meet the preset convergence condition.

[0026] Optionally, in the phase field fracture simulation device for viscoplastic materials, the calculation module includes:

[0027] A correction unit, used for correcting the phase field according to the accumulated strain energy density field to obtain a fracture phase field;

[0028] A calibration convergence unit, used to determine whether the fracture phase field meets a preset convergence condition based on the virtual displacement principle;

[0029] The secondary load step calculation unit is used to calculate the next load step according to the material loading history data, the stress-strain field and the fracture phase field when the fracture phase field meets the preset convergence condition.

[0030] A phase field fracture simulation device for viscoplastic materials, comprising:

[0031] Memory for storing computer programs;

[0032] A processor is used to implement the steps of any one of the above-mentioned methods for simulating phase-field fracture of viscoplastic materials when executing the computer program.

[0033] A computer-readable storage medium having a computer program stored thereon, wherein the computer program, when executed by a processor, implements the steps of any of the above-mentioned phase-field fracture simulation methods for viscoplastic materials.

[0034] The phase field fracture simulation method of viscoplastic material provided by the present invention receives initial deformation field data and material loading history data of the tested viscoplastic material; loads the initial deformation field step by step according to the material loading history data to obtain final field cumulative deformation energy density field data; wherein, after each calculation of a load step is completed, a corresponding stress-strain field and a phase field are generated, and after the corresponding cumulative strain energy density field is determined by the stress-strain field, the phase field is corrected according to the cumulative strain energy density field to obtain a fracture phase field; and the calculation of the next load step is performed according to the material loading history data, the stress-strain field and the fracture phase field.

[0035] On the basis of the traditional phase field simulation based on load steps, the present invention further uses the cumulative strain energy density field corresponding to the latest stress-strain field to correct the phase field after the stress-strain field corresponding to each load step is calculated, and obtains the fracture phase field of the load step. This simulates the coupling effect of the initiation of microcracks inside the material and the complex nonlinear deformation behavior of the material itself as the load is applied after the viscoplastic material is subjected to complex external loads such as fatigue and creep, and its influence on the calculation of the stress-strain field and damage failure process of the material, which is closer to the actual situation of the viscoplastic material bearing the load, improves the simulation accuracy, and thus indicates the optimization and adjustment direction for related preparation and process work. The present invention also provides a phase field fracture simulation device, equipment and computer-readable storage medium for viscoplastic materials with the above-mentioned beneficial effects. BRIEF DESCRIPTION OF THE DRAWINGS

[0036] In order to more clearly illustrate the embodiments of the present invention or the technical solutions of the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0037] Figure 1 A flowchart of a specific implementation of the phase field fracture simulation method for viscoplastic materials provided by the present invention;

[0038] Figure 2 A structural diagram of a specific implementation of the phase-field fracture simulation device for viscoplastic materials provided by the present invention. DETAILED DESCRIPTION

[0039] In order to enable those skilled in the art to better understand the scheme of the present invention, the present invention is further described in detail below in conjunction with the accompanying drawings and specific implementation methods. Obviously, the described embodiments are only part of the embodiments of the present invention, rather than all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative work are within the scope of protection of the present invention.

[0040] The core of the present invention is to provide a phase field fracture simulation method for viscoplastic materials, and a flow chart of a specific implementation method thereof is shown as follows: Figure 1 As shown, it is called specific implementation mode 1, including:

[0041] S101: receiving initial deformation field data and material loading history data of a tested viscoplastic material.

[0042] The initial deformation field data includes a geometrically determined model, given material mechanical / physical property parameters and boundary conditions; the material loading history data includes N load step data, namely, first load step data, second load step data, ... Nth load step data, where N is a positive integer greater than 1.

[0043] S102: The initial deformation field is loaded step by step according to the material loading history data to obtain the final field cumulative deformation energy density field data; wherein, after each load step calculation is completed, the corresponding stress-strain field and phase field are generated, and after the corresponding cumulative strain energy density field is determined by the stress-strain field, the phase field is corrected according to the cumulative strain energy density field to obtain the fracture phase field; and the next load step is calculated according to the material loading history data, the stress-strain field and the fracture phase field.

[0044] It should be noted that the final field cumulative deformation energy density field data in the present invention is calculated step by step through multiple load steps included in the material loading history data. Specifically, this step is disassembled and explained, including:

[0045] Calculating the corresponding first stress-strain field and first phase field according to the initial deformation field and the first load step data;

[0046] Determine a corresponding first accumulated strain energy density field according to the first stress-strain field;

[0047] The first phase field is corrected according to the first accumulated strain energy density field to obtain a first fracture phase field, thereby completing the calculation related to the first load step;

[0048] Enter the second load step related calculation, calculate the corresponding second stress strain field and second phase field according to the first fracture phase field, the first stress strain field and the second load step data, and repeat this cycle, taking the stress strain field and fracture phase field obtained in each previous load step as the calculation basis of the stress strain field of the next load step, until the Nth load step is completed, and the final field cumulative deformation energy density field data can be obtained according to the Nth cumulative strain energy density field obtained at last.

[0049] As a specific implementation manner, the calculation of each load step is completed, a corresponding stress-strain field is generated, and the corresponding cumulative strain energy density field is determined through the stress-strain field, including:

[0050] A1: Complete the calculation of the current load step and generate the corresponding stress-strain field and phase field.

[0051] A2: Based on the virtual displacement principle, it is determined whether the stress-strain field and the phase field both meet the preset convergence conditions.

[0052] Of course, in addition to the virtual displacement principle, convergence may also be determined by other methods, such as determination based on the minimum potential energy principle, which is not limited in the present invention.

[0053] A2: When the stress-strain field and the phase field meet a preset convergence condition, a corresponding cumulative strain energy density field is determined according to the stress-strain field.

[0054] In addition, the phase field is corrected according to the accumulated strain energy density field to obtain a fracture phase field; and the next load step is calculated according to the material loading history data, the stress strain field and the fracture phase field, including:

[0055] B1: The phase field is corrected according to the accumulated strain energy density field to obtain the fracture phase field.

[0056] B2: Based on the virtual displacement principle, it is determined whether the fracture phase field meets the preset convergence condition.

[0057] Of course, in addition to the virtual displacement principle, convergence may also be determined by other methods, such as determination based on the minimum potential energy principle, which is not limited in the present invention.

[0058] B3: When the fracture phase field meets the preset convergence condition, the next load step is calculated according to the material loading history data, the stress-strain field and the fracture phase field.

[0059] In the above two specific implementations, a step of checking whether the calculation results have converged is added to avoid a calculation error in a certain step causing all subsequent calculations to be unobjectionable, resulting in a waste of computing power and reduced computing efficiency.

[0060] As a preferred implementation, after obtaining the final field accumulated deformation energy density field data, the method further includes:

[0061] According to the final field cumulative deformation energy density field data, a visualization result is obtained. The visualization result has a stronger reference significance and a stronger guiding role for the engineering evaluation of viscoplastic materials.

[0062] In addition, the initial deformation field data is model data for mesh encryption in the potential crack area.

[0063] The specific potential crack area can be selected by the staff according to the actual situation. Increasing the mesh density of the interface can significantly improve the calculation accuracy.

[0064] The phase field fracture simulation method of viscoplastic material provided by the present invention receives initial deformation field data and material loading history data of the tested viscoplastic material; loads the initial deformation field step by step according to the material loading history data to obtain final field cumulative deformation energy density field data; wherein, after each calculation of a load step is completed, a corresponding stress-strain field and a phase field are generated, and after the corresponding cumulative strain energy density field is determined by the stress-strain field, the phase field is corrected according to the cumulative strain energy density field to obtain a fracture phase field; and the calculation of the next load step is performed according to the material loading history data, the stress-strain field and the fracture phase field. On the basis of the traditional phase field simulation based on load steps, after the stress-strain field corresponding to each load step is calculated, the phase field is corrected using the accumulated strain energy density field corresponding to the latest stress-strain field to obtain the fracture phase field of the load step. This simulates the coupling effect of the initiation of microcracks inside the material and the complex nonlinear deformation behavior of the material itself as the load is applied after the viscoplastic material is subjected to complex external loads such as fatigue and creep, and its influence on the calculation of the stress-strain field and damage failure process of the material. It is closer to the actual situation of the viscoplastic material bearing the load, improves the simulation accuracy, and thus indicates the optimization and adjustment direction for related preparation and process work.

[0065] The following uses the ABAQUS software platform as an example to illustrate the implementation of the phase field fracture simulation method for viscoplastic materials, including:

[0066] 1) Finite element model establishment

[0067] a) Establish a finite element analysis model containing cracks / defective structures, and perform mesh encryption in areas where cracks may appear or extend to improve the accuracy of the calculation;

[0068] b) Using node nesting technology in the directly established finite element model, a group of puppet units without stiffness are established for subsequent phase field program calculations;

[0069] c) Apply boundary conditions of temperature, displacement and stress to the finite element model according to the actual load borne by the component.

[0070] 2) Material stress and strain calculation

[0071] a) Using Fortran programming language and referring to the requirements of ABAQUS platform for user subroutine secondary development interface (ABAQUS software is a general finite element software developed by Dassault), compile the user material subroutine plug-in (UMAT) for viscoplastic material stress and strain calculation;

[0072] b) The finite element software inputs the deformation field of the structure into UMAT and calculates the stress-strain field at each integration point based on the material loading history;

[0073] c) Return the stress and strain field of the integration point to the ABAQUS finite element software and calculate whether the structural stress and strain field meets the convergence conditions by the virtual displacement principle;

[0074] d) If the convergence condition is met, the next step of calculation is carried out; otherwise, the program subdivides the load step and returns to step 2-a for recalculation;

[0075] e) Based on the stress field and strain field of the integration point, calculate the cumulative strain energy density of the integration point, and return the result to ABAQUS.

[0076] 3) Crack propagation phase field iteration analysis

[0077] a) Use Fortran programming language and refer to the requirements of ABAQUS platform for user subroutine secondary development interface to compile material phase field calculation user subroutine plug-in (UEL);

[0078] b) The finite element software inputs the cumulative deformation energy density field calculated by UMAT into UEL. UEL corrects the stress result according to the phase field value of each integral point and assembles a new stiffness matrix;

[0079] c) Return the new stiffness matrix to ABAQUS software, and judge whether the stress-strain field and phase field are converged based on the principle of minimum potential energy: if converged, store the relevant variables and proceed to the next step of calculation; if not converged, the program will subdivide the load steps and return to step 2-a for recalculation.

[0080] 4) Visual post-processing of calculation results

[0081] a) Use internal variables to store and transfer the calculation results of each integration point UEL to UMAT;

[0082] b) Create separate user-defined variables in UMAT to store the calculated phase field values;

[0083] c) Use the ABAUS software post-processing module to interpolate and render the UEL calculation results to obtain the phase field cloud map of crack extension.

[0084] The following is an application example of the above-mentioned viscoplastic phase field simulation method, using UMAT and UEL program blocks to calculate the stress, strain and phase field information of the unit. At the same time, in order to ensure the convergence of the calculation, before each program block ends and enters the next integral point calculation, the convergence of the program needs to be judged. If the convergence requirements are not met, the load step needs to be subdivided and recalculated.

[0085] The above method is used to calculate the mechanical response of a certain equiaxed nickel-based superalloy under creep-fatigue load. The test results are in good agreement with the simulation results, indicating that the accuracy of the viscoplastic phase field fracture calculation program is good. Through these examples, the effectiveness of the program plug-in is further verified, indicating that the program can more accurately reflect the long-term and short-term mechanical properties of the material, and has certain engineering significance.

[0086] This viscoplastic phase field fracture calculation program can be combined with general finite element software to calculate the initiation and propagation of cracks in nickel-based high-temperature alloys under complex loads. In order to demonstrate the calculation performance of the program, this paper calculates the propagation law of two-dimensional and three-dimensional cracks inside a certain high-temperature alloy under creep-fatigue interaction loads. The calculation results show that the calculation results of crack propagation of this program are consistent with the experimental results, with good accuracy and calculation efficiency. These examples show that the program plug-in has good calculation accuracy and can be used to calculate the initiation and propagation behavior of cracks in actual engineering structures. The developed program plug-in has certain significance.

[0087] The phase field fracture simulation device for viscoplastic materials provided by an embodiment of the present invention is introduced below. The phase field fracture simulation device for viscoplastic materials described below and the phase field fracture simulation method for viscoplastic materials described above can be referred to each other.

[0088] Figure 2 The structural block diagram of the phase field fracture simulation device for viscoplastic materials provided by the embodiment of the present invention is shown in FIG. Figure 2 Phase field fracture simulation setup for viscoplastic materials can include:

[0089] The receiving module 100 is used to receive the initial deformation field data and material loading history data of the tested viscoplastic material;

[0090] The calculation module 200 is used to gradually load the initial deformation field with load steps according to the material loading history data to obtain the final field cumulative deformation energy density field data; wherein, after each load step calculation is completed, the corresponding stress-strain field and phase field are generated, and after the corresponding cumulative strain energy density field is determined by the stress-strain field, the phase field is corrected according to the cumulative strain energy density field to obtain the fracture phase field; and the next load step calculation is performed according to the material loading history data, the stress-strain field and the fracture phase field.

[0091] As a preferred implementation, the calculation module 200 includes:

[0092] The stress-strain calculation unit is used to complete the calculation of the current load step and generate the corresponding stress-strain field and phase field;

[0093] A calculation convergence checking unit, used for judging whether the stress-strain field and the phase field both meet preset convergence conditions based on the virtual displacement principle;

[0094] The cumulative strain energy calculation unit is used to determine the corresponding cumulative strain energy density field according to the stress strain field when the stress strain field and the phase field meet the preset convergence condition.

[0095] As a preferred implementation, the calculation module 200 includes:

[0096] A correction unit, used for correcting the phase field according to the accumulated strain energy density field to obtain a fracture phase field;

[0097] A calibration convergence unit, used to determine whether the fracture phase field meets a preset convergence condition based on the virtual displacement principle;

[0098] The secondary load step calculation unit is used to calculate the next load step according to the material loading history data, the stress-strain field and the fracture phase field when the fracture phase field meets the preset convergence condition.

[0099] As a preferred implementation, the calculation module 200 further includes:

[0100] The visualization unit is used to obtain a visualization result according to the final field accumulated deformation energy density field data.

[0101] The phase field fracture simulation method of viscoplastic material provided by the present invention comprises a receiving module 100 for receiving initial deformation field data and material loading history data of the tested viscoplastic material; a calculating module 200 for gradually loading the initial deformation field with load steps according to the material loading history data to obtain final field cumulative deformation energy density field data; wherein, after each calculation of a load step is completed, a corresponding stress-strain field and a phase field are generated, and after the corresponding cumulative strain energy density field is determined by the stress-strain field, the phase field is corrected according to the cumulative strain energy density field to obtain a fracture phase field; and the calculation of the next load step is performed according to the material loading history data, the stress-strain field and the fracture phase field. On the basis of the traditional phase field simulation based on load steps, after the stress-strain field corresponding to each load step is calculated, the phase field is corrected using the accumulated strain energy density field corresponding to the latest stress-strain field to obtain the fracture phase field of the load step. This simulates the coupling effect of the initiation of microcracks inside the material and the complex nonlinear deformation behavior of the material itself as the load is applied after the viscoplastic material is subjected to complex external loads such as fatigue and creep, and its influence on the calculation of the stress-strain field and damage failure process of the material. It is closer to the actual situation of the viscoplastic material bearing the load, improves the simulation accuracy, and thus indicates the optimization and adjustment direction for related preparation and process work.

[0102] The phase field fracture simulation device for viscoplastic materials of this embodiment is used to implement the aforementioned phase field fracture simulation method for viscoplastic materials. Therefore, the specific implementation methods of the phase field fracture simulation device for viscoplastic materials can be seen in the embodiment part of the phase field fracture simulation method for viscoplastic materials in the previous text. For example, the receiving module 100 and the computing module 200 are respectively used to implement steps S101 and S102 in the aforementioned phase field fracture simulation method for viscoplastic materials. Therefore, its specific implementation methods can refer to the descriptions of the corresponding embodiments of each part and will not be repeated here.

[0103] A phase field fracture simulation device for viscoplastic materials, comprising:

[0104] Memory for storing computer programs;

[0105] A processor is used to implement the steps of the phase field fracture simulation method of viscoplastic materials as described in any one of the above when executing the computer program. The phase field fracture simulation method of viscoplastic materials provided by the present invention receives the initial deformation field data and material loading history data of the tested viscoplastic material; the initial deformation field is gradually loaded with load steps according to the material loading history data to obtain the final field cumulative deformation energy density field data; wherein, after each load step calculation is completed, the corresponding stress-strain field and phase field are generated, and after the corresponding cumulative strain energy density field is determined by the stress-strain field, the phase field is corrected according to the cumulative strain energy density field to obtain the fracture phase field; and the next load step calculation is performed according to the material loading history data, the stress-strain field and the fracture phase field. On the basis of the traditional phase field simulation based on load steps, after the stress-strain field corresponding to each load step is calculated, the phase field is corrected using the accumulated strain energy density field corresponding to the latest stress-strain field to obtain the fracture phase field of the load step. This simulates the coupling effect of the initiation of microcracks inside the material and the complex nonlinear deformation behavior of the material itself as the load is applied after the viscoplastic material is subjected to complex external loads such as fatigue and creep, and its influence on the calculation of the stress-strain field and damage failure process of the material. It is closer to the actual situation of the viscoplastic material bearing the load, improves the simulation accuracy, and thus indicates the optimization and adjustment direction for related preparation and process work.

[0106] A computer-readable storage medium having a computer program stored thereon, wherein when the computer program is executed by a processor, the steps of the phase field fracture simulation method for viscoplastic materials as described in any one of the above are implemented. The phase field fracture simulation method for viscoplastic materials provided by the present invention receives the initial deformation field data and material loading history data of the tested viscoplastic material; the initial deformation field is gradually loaded with load steps according to the material loading history data to obtain the final field cumulative deformation energy density field data; wherein, after each load step calculation is completed, the corresponding stress-strain field and phase field are generated, and after the corresponding cumulative strain energy density field is determined by the stress-strain field, the phase field is corrected according to the cumulative strain energy density field to obtain the fracture phase field; and the next load step calculation is performed according to the material loading history data, the stress-strain field and the fracture phase field. On the basis of the traditional phase field simulation based on load steps, after the stress-strain field corresponding to each load step is calculated, the phase field is corrected using the accumulated strain energy density field corresponding to the latest stress-strain field to obtain the fracture phase field of the load step. This simulates the coupling effect of the initiation of microcracks inside the material and the complex nonlinear deformation behavior of the material itself as the load is applied after the viscoplastic material is subjected to complex external loads such as fatigue and creep, and its influence on the calculation of the stress-strain field and damage failure process of the material. It is closer to the actual situation of the viscoplastic material bearing the load, improves the simulation accuracy, and thus indicates the optimization and adjustment direction for related preparation and process work.

[0107] In this specification, each embodiment is described in a progressive manner, and each embodiment focuses on the differences from other embodiments. The same or similar parts between the embodiments can be referred to each other. For the device disclosed in the embodiment, since it corresponds to the method disclosed in the embodiment, the description is relatively simple, and the relevant parts can be referred to the method part.

[0108] It should be noted that, in this specification, relational terms such as first and second, etc. are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Moreover, the terms "comprise", "include" or any other variants thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or device. In the absence of further restrictions, an element defined by the statement "comprises a ..." does not exclude the presence of other identical elements in the process, method, article or device including the element.

[0109] Professionals may further appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, computer software, or a combination of the two. In order to clearly illustrate the interchangeability of hardware and software, the composition and steps of each example have been generally described in the above description according to function. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professionals and technicians may use different methods to implement the described functions for each specific application, but such implementation should not be considered to be beyond the scope of the present invention.

[0110] The steps of the method or algorithm described in conjunction with the embodiments disclosed herein may be implemented directly using hardware, a software module executed by a processor, or a combination of the two. The software module may be placed in a random access memory (RAM), a memory, a read-only memory (ROM), an electrically programmable ROM, an electrically erasable programmable ROM, a register, a hard disk, a removable disk, a CD-ROM, or any other form of storage medium known in the art.

[0111] The above is a detailed introduction to the phase field fracture simulation method, device, equipment and computer-readable storage medium for viscoplastic materials provided by the present invention. Specific examples are used herein to illustrate the principles and implementation methods of the present invention. The description of the above embodiments is only used to help understand the method of the present invention and its core idea. It should be pointed out that for ordinary technicians in this technical field, without departing from the principles of the present invention, the present invention can also be improved and modified in a number of ways, and these improvements and modifications also fall within the scope of protection of the claims of the present invention.

Claims

1. A phase field fracture simulation method for viscoplastic materials, characterized in that: include: Receiving initial deformation field data and material loading history data of the tested viscoplastic material; The initial deformation field is loaded step by step according to the material loading history data to obtain the final field cumulative deformation energy density field data; wherein, after each load step calculation is completed, the corresponding stress-strain field and phase field are generated, and after the corresponding cumulative strain energy density field is determined by the stress-strain field, the phase field is corrected according to the cumulative strain energy density field to obtain the fracture phase field; and the next load step calculation is performed according to the material loading history data, the stress-strain field and the fracture phase field; The calculation of each load step is completed, the corresponding stress-strain field and phase field are generated, and the corresponding cumulative strain energy density field is determined by the stress-strain field, including: Complete the calculation of the current load step and generate the corresponding stress-strain field and phase field; Based on the virtual displacement principle, judging whether the stress-strain field and the phase field both meet the preset convergence conditions; When the stress-strain field and the phase field meet a preset convergence condition, determining a corresponding cumulative strain energy density field according to the stress-strain field; The step of correcting the phase field according to the accumulated strain energy density field to obtain a fracture phase field; and calculating the next load step according to the material loading history data, the stress-strain field and the fracture phase field comprises: Correcting the phase field according to the accumulated strain energy density field to obtain a fracture phase field; Judging whether the fracture phase field meets a preset convergence condition based on the virtual displacement principle; When the fracture phase field meets the preset convergence condition, the calculation of the next load step is performed according to the material loading history data, the stress-strain field and the fracture phase field.

2. The phase field fracture simulation method of viscoplastic material according to claim 1, characterized in that: After obtaining the final field accumulated deformation energy density field data, the method further includes: A visualization result is obtained according to the final-field accumulated deformation energy density field data.

3. The phase field fracture simulation method of viscoplastic material according to claim 1, characterized in that: The initial deformation field data are model data for mesh encryption in the potential crack region.

4. A phase field fracture simulation device for viscoplastic materials, characterized in that: include: A receiving module, used for receiving the initial deformation field data and material loading history data of the tested viscoplastic material; A calculation module, used for gradually loading the initial deformation field with load steps according to the material loading history data to obtain the final field cumulative deformation energy density field data; wherein, after completing the calculation of each load step, the corresponding stress-strain field and phase field are generated, and after the corresponding cumulative strain energy density field is determined by the stress-strain field, the phase field is corrected according to the cumulative strain energy density field to obtain the fracture phase field; and the calculation of the next load step is performed according to the material loading history data, the stress-strain field and the fracture phase field; The calculation module comprises: The stress-strain calculation unit is used to complete the calculation of the current load step and generate the corresponding stress-strain field and phase field; A calculation convergence checking unit, used for judging whether the stress-strain field and the phase field both meet preset convergence conditions based on the virtual displacement principle; a cumulative strain energy calculation unit, configured to determine a corresponding cumulative strain energy density field according to the stress strain field when the stress strain field and the phase field meet a preset convergence condition; The calculation module comprises: A correction unit, used for correcting the phase field according to the accumulated strain energy density field to obtain a fracture phase field; A calibration convergence unit, used to determine whether the fracture phase field meets a preset convergence condition based on the virtual displacement principle; The secondary load step calculation unit is used to calculate the next load step according to the material loading history data, the stress-strain field and the fracture phase field when the fracture phase field meets the preset convergence condition.

5. A phase field fracture simulation device for viscoplastic materials, characterized in that: include: Memory for storing computer programs; A processor is used to implement the steps of the phase field fracture simulation method for viscoplastic materials as claimed in any one of claims 1 to 3 when executing the computer program.

6. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the steps of the phase-field fracture simulation method for viscoplastic materials according to any one of claims 1 to 3 are implemented.

Images