Uncured fiber reinforced metal laminate forming simulation method, system, medium and equipment

Abstract

The invention relates to the field of finite element analysis and computer simulation, in particular to an uncured fiber reinforced metal laminate forming simulation method, system, medium and equipment, and the method comprises the following steps: preprocessing a key variable according to the input key variable of uncured fiber to obtain a fiber strain direction; obtaining stretching and shearing strain of the fiber based on the strain direction; calculating the current shear angle, shear modulus and Young modulus based on the tensile and shear strain of the fiber, and updating the nominal stress of the fiber; and performing wrinkling judgment according to the updated nominal stress of the fiber, and performing damage judgment and evolution after the wrinkling judgment. According to the method, the weaving angle, fracture damage and wrinkling conditions of the fiber fabric can be effectively predicted and controlled in the forming process.

Description

Technical Field

[0001] The present invention relates to the technical field of finite element analysis and computer simulation, and in particular to a method, system, medium and equipment for simulating the forming of uncured fiber reinforced metal laminates based on a VUMAT subroutine. Background Art

[0002] Fiber Metal Laminate (FML) is widely used in aerospace, automotive, and shipbuilding fields due to its excellent strength and lightweight properties. However, during the uncured forming process of FML, the braiding angle of the fiber fabric will change with stress and deformation, and this change directly affects the fiber direction and mechanical properties of the final forming. In addition, FML may have defects such as fiber breakage damage and wrinkling during the forming process, which will reduce material properties and affect structural safety. Therefore, accurately simulating the changes in braiding angles and predicting breakage damage and wrinkling are difficult problems that need to be solved in forming simulation technology. Traditional methods cannot predict fabric shear angle deformation, breakage, and wrinkling, resulting in inaccurate numerical simulations. Summary of the invention

[0003] In view of the above problems, the purpose of the present invention is to provide a method, system, medium and equipment for simulating the forming of uncured fiber reinforced metal laminates, which can effectively predict and control the braiding angle, fracture damage and wrinkling of fiber fabrics during the forming process.

[0004] To achieve the above-mentioned purpose, in the first aspect, the technical solution adopted by the present invention is: a method for simulating the forming of uncured fiber-reinforced metal laminates, which includes: according to the key variables of the input uncured fiber, the key variables are preprocessed to obtain the fiber strain direction, and the tensile and shear strains of the fiber are obtained based on the strain direction; based on the tensile and shear strains of the fiber, the current shear angle, shear modulus and Young's modulus are calculated, and the nominal stress of the fiber is updated; according to the updated nominal stress of the fiber, the damage criterion coefficient is calculated, wrinkling judgment is performed through the shear angle, and after the wrinkling judgment, damage judgment and evolution are performed through the damage criterion coefficient.

[0005] Furthermore, the key variables of the unsolidified fiber include: the main yarn direction, the polynomial coefficient of Young's modulus, the polynomial coefficient of shear behavior, the tensile strength of the fiber bundle, the compressive strength of the fiber bundle, the operating factor β for calculating the fiber damage variable; the critical braiding angle size γ for the set fabric wrinkling w .

[0006] Furthermore, the fiber strain direction is obtained after preprocessing the key variables, including:

[0007] The deformation gradient tensor and stretch tensor are established based on key variables to characterize the extension or compression characteristics of the fiber material under external loads;

[0008] Based on the deformation gradient tensor and the stretch tensor, the rotation tensor and the local coordinate tensor of the rotation tensor are calculated to reflect the spatial rotation of the fiber direction caused by the deformation, and the rotation tensor is converted into the local coordinate system;

[0009] Based on the rotation tensor and its local coordinate tensor, the current fiber orientation is determined in real time to capture the changes in fiber braiding angle during the forming process;

[0010] The transformation matrix is ​​generated through the relative relationship between the fiber direction and the local coordinate system to convert the global loading condition into the stress state in the fiber direction;

[0011] The fiber strain direction is calculated based on the transformation matrix to obtain the tensile and shear strains.

[0012] Furthermore, based on the tensile and shear strains of the fiber, the current shear angle, shear modulus and Young's modulus are calculated, and the nominal stress of the fiber is updated, including:

[0013] σ f，t+1 =σ f，t +dε f ·C

[0014] Among them, σ f，t+1 is the updated nominal stress of the fiber; C is the current stiffness matrix, which consists of the shear modulus and Young's modulus.

[0015] Furthermore, wrinkle determination is performed based on the shear angle, including:

[0016] Determine whether the absolute value of the shear angle exceeds the critical braiding angle γ for fabric wrinkling w If it exceeds, wrinkling is determined to have occurred, otherwise wrinkling has not occurred.

[0017] Furthermore, damage judgment and evolution are: based on the 3D-Hashin criterion, the damage of the material during the forming process is judged and evolved, including:

[0018] By calculating different damage modes in the fiber direction, the 3D-Hashin damage criterion coefficient f is calculated to determine the damage variable By damage variable Identify possible failure scenarios;

[0019] When the damage state is determined, the stress in the fiber direction, the matrix stress and the corresponding state variables are gradually updated to reflect the damage accumulation and performance degradation of the material during loading and simulate the damage evolution process.

[0020] Furthermore, the damage variable for:

[0021]

[0022] Where e is the natural logarithm.

[0023] In the second aspect, the technical solution adopted by the present invention is: a simulation system for forming an uncured fiber reinforced metal laminate, which includes: a preprocessing module, which preprocesses the key variables according to the input uncured fiber key variables to obtain the fiber strain direction, and obtains the tensile and shear strains of the fiber based on the strain direction; a stress calculation update module, which calculates the current shear angle, shear modulus and Young's modulus based on the tensile and shear strains of the fiber, and updates the nominal stress of the fiber; a judgment evolution module, which calculates the damage criterion coefficient according to the updated nominal stress of the fiber, performs wrinkling judgment through the shear angle, and performs damage judgment and evolution through the damage criterion coefficient after the wrinkling judgment.

[0024] In a third aspect, the technical solution adopted by the present invention is: a computer-readable storage medium storing one or more programs, wherein the one or more programs include instructions, and when the instructions are executed by a computing device, the computing device executes any one of the above methods.

[0025] In a fourth aspect, the technical solution adopted by the present invention is: a computing device, comprising: one or more processors, a memory and one or more programs, wherein the one or more programs are stored in the memory and configured to be executed by the one or more processors, and the one or more programs include instructions for executing any of the above methods.

[0026] The present invention adopts the above technical solution, which has the following advantages:

[0027] The present invention uses a powerful material model to accurately simulate the response of composite materials under different loading conditions by considering the anisotropy, damage evolution, shear effect and nonlinear behavior of large deformation of textile materials. Through the combination of polynomial fitting, damage evolution and fourth-order Runge-Kutta integration method, the model has high accuracy and robustness in theory, and has good computational efficiency and flexibility, which is suitable for the design and analysis of high-performance composite materials such as aerospace. BRIEF DESCRIPTION OF THE DRAWINGS

[0028] Figure 1 is an overall flow chart of the uncured fiber reinforced metal laminate forming simulation method in an embodiment of the present invention;

[0029] Figure 2 It is a detailed flow chart of the uncured fiber reinforced metal laminate forming simulation method in an embodiment of the present invention. DETAILED DESCRIPTION

[0030] In order to solve the problem that the prior art cannot predict the shear angle deformation, fracture and wrinkling of fabrics, resulting in inaccurate numerical simulation, the present invention provides a method, system, medium and equipment for simulating the forming of uncured fiber reinforced metal laminates, and uses Abaqus software to perform structural, thermal and fluid dynamics analysis. The focus is on material modeling, especially the simulation of nonlinear behavior under dynamic or impact conditions, and special attention is paid to the change of fabric braiding angle, fracture damage prediction and wrinkling phenomenon simulation during the uncured forming process, so as to improve the accuracy of the forming process and material utilization.

[0031] In order to make the purpose, technical solution and advantages of the embodiment of the present invention clearer, the technical solution of the embodiment of the present invention will be clearly and completely described below in conjunction with the drawings of the embodiment of the present invention. Obviously, the described embodiment is a part of the embodiment of the present invention, not all of the embodiments. Based on the described embodiment of the present invention, all other embodiments obtained by ordinary technicians in this field belong to the scope of protection of the present invention.

[0032] It should be noted that the terms used herein are only for describing specific embodiments and are not intended to limit the exemplary embodiments according to the present application. As used herein, unless the context clearly indicates otherwise, the singular form is also intended to include the plural form. In addition, it should be understood that when the terms "comprise" and / or "include" are used in this specification, it indicates the presence of features, steps, operations, devices, components and / or combinations thereof.

[0033] In one embodiment of the present invention, a method for simulating the forming of an uncured fiber-reinforced metal laminate is provided to effectively predict and control the braiding angle, fracture damage and wrinkling of the fiber fabric during the forming process. Figure 1 , Figure 2 As shown, the method comprises the following steps:

[0034] 1) According to the key variables of the input uncured fiber, the fiber strain direction is obtained after preprocessing the key variables, and the tensile and shear strains of the fiber are obtained based on the strain direction, providing accurate stress-strain data for subsequent fracture damage and wrinkling prediction;

[0035] 2) Based on the tensile and shear strains of the fiber, the current shear angle, shear modulus and Young's modulus are calculated, and the nominal stress of the fiber is updated;

[0036] 3) The damage criterion coefficient f is calculated based on the updated nominal stress of the fiber, and wrinkling is judged by the shear angle. After wrinkling, damage judgment and evolution are performed using the damage criterion coefficient f.

[0037] In step 1), the key variables of the uncured fiber include the main yarn directions (warp and weft), the polynomial coefficients A of the Young's modulus, n (up to fifth order), polynomial coefficients B for shear behavior n (up to fifth order), tensile strength of fiber bundle X T , the compressive strength of the fiber bundle X C , the operating factor β (≥2) for calculating the fiber damage variable; the critical braiding angle γ for setting fabric wrinkling w (radians), f 0,1 and f 0,2 are the initial direction vectors of the warp and weft yarns of the fabric, respectively.

[0038] Specifically, for example, the key variables entered are:

[0039] props(1)-props(2) = main yarn direction 1 (warp and weft);

[0040] props(4)-props(8) = coefficients of the polynomial defining Young's modulus (up to fourth order);

[0041] props(10)-props(15) = polynomial coefficients defining shear behavior (up to fifth order);

[0042] props(17)-props(18) = shear angle of wrinkle criterion;

[0043] props(20)-props(22)=value that determines damage.

[0044] The shear angle of the wrinkle criterion is used to set the critical weaving angle γ of the fabric wrinkle w ;

[0045] The operating factor β is calculated by determining the damage values ​​corresponding to the tensile strength XT of the fiber bundle, the compressive strength XC of the fiber bundle and the fiber damage variable in order.

[0046] In the above step 1), the fiber strain direction is obtained after preprocessing the key variables, including:

[0047] 1.1) Establish the deformation gradient tensor F and the tensile tensor U based on the key variables to characterize the stretching or compression characteristics of the fiber material under the applied load, which serves as the basis for subsequent calculations;

[0048] 1.2) Based on the deformation gradient tensor and the stretch tensor, calculate the rotation tensor R and the local coordinate tensor of the rotation tensor R To reflect the spatial rotation of the fiber direction caused by deformation and transform it into the local coordinate system, so as to analyze the deformation response of the fiber more accurately;

[0049] Specifically, the rotation tensor R is:

[0050] R=F·U -1

[0051] Local Coordinate Tensor for:

[0052]

[0053] in, is the initial value of the local coordinate tensor; α represents the rotation angle; i represents the current running node.

[0054] 1.3) Based on the rotation tensor R and its local coordinate tensor Determine the current fiber direction f in real time α , to capture the changes in fiber braiding angle during the forming process;

[0055] The current fiber direction is the direction of the warp yarn and / or the weft yarn.

[0056] Specifically, the current fiber direction f α for:

[0057]

[0058] 1.4) Through the fiber direction f α The relative relationship with the local coordinate system generates the transformation matrix T α , converting the global loading condition into the stress state in the fiber direction;

[0059] Transformation matrix T α for:

[0060]

[0061] 1.5) Based on the transformation matrix T α Calculate the strain increment dε in the fiber direction f , tensile and shear strains are obtained, providing accurate stress-strain data for subsequent fracture damage and wrinkling prediction.

[0062] The stress state in the fiber direction is the fiber direction strain increment dε corresponding to the warp and weft yarns f :

[0063] dε f =T′ α ·dε·T α

[0064] Where dε is the strain increment at the current node; T′ α Represents the transformation matrix T α The transposed matrix of .

[0065] In the above step 2), the current shear angle γ and shear modulus G are calculated based on the fiber direction strain increment of the fiber. ij and Young's modulus E i , update the nominal stress σ of the fiber f , specifically:

[0066] Among them, the current shear angle γ is:

[0067] γ=γ+dε f +dσ f

[0068] The four-node Runge-Kutta method is used for interpolation, and the shear modulus G is:

[0069]

[0070] The four-node Runge-Kutta method is used for interpolation, and the Young's modulus E j (j=1, 2) is:

[0071]

[0072] Update the nominal fiber stress σ f for:

[0073] σ f，t+1 =σ f，t +dε f ·C

[0074] Where t is the time; C is the current stiffness matrix (composed of shear modulus G and Young's modulus E j composition):

[0075]

[0076] In the above step 3), wrinkling is judged by the shear angle, specifically:

[0077] Determine whether the absolute value of the shear angle |γ| exceeds the critical weaving angle γ for fabric wrinkling w If it exceeds, wrinkling is determined to have occurred, otherwise wrinkling has not occurred.

[0078] Specific:

[0079]

[0080] In the formula, Indicates whether wrinkling occurs. Wrinkling occurs.

[0081] That is, whether the absolute value of the current fabric shear angle is greater than the given value, if so, 1 represents wrinkling; otherwise, It is 0 (that is, the state variable SDV7 for judging whether the fabric is wrinkled is assigned to 1), which means it is not wrinkled (the state variable SDV7 for judging whether the fabric is wrinkled is assigned to 0).

[0082] In the above step 3), damage judgment and evolution are: judging and evolving the damage of the material during the forming process based on the 3D-Hashin criterion.

[0083] Among them, the damage of the material during the forming process is judged and evolved based on the 3D-Hashin criterion, including the following steps:

[0084] 3.1) By calculating different damage modes in the fiber direction (such as fiber tensile fracture, fiber compression failure, etc.), the 3D-Hashin damage criterion coefficient f is calculated to determine the damage variable By damage variable Identify possible failure scenarios;

[0085] Specifically, the damage variable for:

[0086]

[0087] Where e is the natural logarithm; β represents the operating factor for the calculation of fiber damage variables.

[0088] If the damage variable If it is 1, it means that the material at that point is completely failed (that is, it is in a completely damaged state), and the corresponding unit is deleted (the state variable SDV12 used to determine whether to delete the unit is assigned to 0); otherwise, the damage variable A value other than 1 indicates no damage.

[0089] 3.2) When the damage state is determined, the true stress in the fiber direction and the corresponding state variables are gradually updated to reflect the damage accumulation and performance degradation of the material during loading and accurately simulate the damage evolution process.

[0090] The corresponding state variables include axial stress and shear stress along the yarn direction. For example, they include SDV1 to SDV4, where SDV1 is the axial stress along the yarn direction 1, SDV2 is the shear stress along the yarn direction 1, SDV3 is the axial stress along the yarn direction 2, and SDV4 is the shear stress along the yarn direction 1.

[0091] Update the true stress σ f，t :

[0092]

[0093] Update the stress σ at this point:

[0094] σ=T′ α ·σ f，t ·T α

[0095] The output variables are:

[0096] state(*,1) = axial stress along yarn direction 1;

[0097] state(*,2) = shear stress along yarn direction 1;

[0098] state(*,3) = axial stress along yarn direction 2;

[0099] state(*,4) = shear stress along yarn direction 2;

[0100] state(*,5) = shear angle in radians;

[0101] state(*,6) = shear angle in degrees;

[0102] state(*,7) = fold criterion;

[0103] state(*,8) = fiber strain along yarn direction 1;

[0104] state(*,9) = fiber strain along yarn direction 2;

[0105] state(*,10) = exponential damage evolution index along the yarn direction 1;

[0106] state(*,11) = exponential damage evolution index along yarn direction 2;

[0107] state(*,12)=deleted flag;

[0108] state(*,13) = nominal stress of axial stress along yarn direction 1;

[0109] state(*,14) = nominal stress of shear stress along yarn direction 1;

[0110] state(*,15) = nominal stress of axial stress along yarn direction 2;

[0111] state(*,16) = nominal stress of shear stress along yarn direction 2.

[0112] In one embodiment of the present invention, a system for simulating the forming of an uncured fiber reinforced metal laminate is provided, comprising:

[0113] A preprocessing module, according to the key variables of the input uncured fiber, preprocesses the key variables to obtain the fiber strain direction, and obtains the tensile and shear strains of the fiber based on the strain direction;

[0114] The stress calculation update module calculates the current shear angle, shear modulus and Young's modulus based on the tensile and shear strains of the fiber, and updates the nominal stress of the fiber;

[0115] The judgment evolution module calculates the damage criterion coefficient f according to the updated nominal stress of the fiber, makes wrinkle judgment through the shear angle, and after wrinkle judgment, makes damage judgment and evolution through the damage criterion coefficient.

[0116] In the above embodiment, the key variables of the uncured fiber include: the main yarn direction, the polynomial coefficient of Young's modulus, the polynomial coefficient of shear behavior, the tensile strength of the fiber bundle, the compressive strength of the fiber bundle, the operating factor β for calculating the fiber damage variable; the critical braiding angle γ for the fabric wrinkling w .

[0117] In the above embodiment, the fiber strain direction is obtained after preprocessing the key variables, including:

[0118] The deformation gradient tensor and stretch tensor are established based on key variables to characterize the extension or compression characteristics of the fiber material under external loads;

[0119] Based on the deformation gradient tensor and the stretch tensor, the rotation tensor and its local coordinate tensor are calculated to reflect the spatial rotation of the fiber direction caused by the deformation, and are converted into the local coordinate system;

[0120] Based on the rotation tensor and its local coordinate tensor, the current fiber orientation is determined in real time to capture the changes in fiber braiding angle during the forming process;

[0121] The transformation matrix is ​​generated through the relative relationship between the fiber direction and the local coordinate system to convert the global loading condition into the stress state in the fiber direction;

[0122] The fiber strain direction is calculated based on the transformation matrix to obtain the tensile and shear strains.

[0123] In the above embodiment, based on the tensile and shear strains of the fiber, the current shear angle, shear modulus and Young's modulus are calculated, and the nominal stress of the fiber is updated, including:

[0124] σ f，t+1 =σ f，t +dε f ·C

[0125] Among them, σ f，t+1 is the nominal stress of the updated fiber; C is the current stiffness matrix, which is composed of the shear modulus G and Young's modulus Ej composition.

[0126] In the above embodiment, wrinkle determination is performed based on the shear angle, including:

[0127] Determine whether the absolute value of the shear angle exceeds the critical braiding angle γ for fabric wrinkling w If it exceeds, wrinkling is determined to have occurred, otherwise wrinkling has not occurred.

[0128] In the above embodiment, damage judgment and evolution are: judging and evolving the damage of the material during the forming process based on the 3D-Hashin criterion.

[0129] In this embodiment, the damage of the material during the forming process is judged and evolved based on the 3D-Hashin criterion, including:

[0130] Identify possible failure scenarios by calculating different damage modes along the fiber directions;

[0131] When the damage state is determined, the stress in the fiber direction, the matrix stress and the corresponding state variables are gradually updated to reflect the damage accumulation and performance degradation of the material during loading and simulate the damage evolution process.

[0132] In this embodiment, the damage variable for:

[0133]

[0134] Where e is the natural logarithm.

[0135] The system provided in this embodiment is used to execute the above-mentioned method embodiments. Please refer to the above-mentioned embodiments for specific processes and detailed contents, which will not be repeated here.

[0136] A computing device provided in one embodiment of the present invention may be a terminal, which may include: a processor, a communications interface, a memory, a display screen, and an input device. The processor, the communications interface, and the memory communicate with each other through a communication bus. The processor is used to provide computing and control capabilities. The memory includes a non-volatile storage medium and an internal memory, wherein the non-volatile storage medium stores an operating system and a computer program, and the computer program is executed by the processor to implement the methods in the above embodiments; the internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The communication interface is used to communicate with an external terminal in a wired or wireless manner, and the wireless manner may be implemented through WIFI, a management network, NFC (near field communication) or other technologies. The display screen may be a liquid crystal display screen or an electronic ink display screen, and the input device may be a touch layer covered on the display screen, or a key, trackball or touchpad provided on the housing of the computing device, or an external keyboard, touchpad or mouse, etc. The processor may call the logic instructions in the memory.

[0137] In addition, the logic instructions in the above-mentioned memory can be implemented in the form of software functional units and can be stored in a computer-readable storage medium when sold or used as an independent product. Based on such an understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art or the part of the technical solution, can be embodied in the form of a software product, which is stored in a storage medium and includes several instructions for a computer device (which can be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the method described in each embodiment of the present invention. The aforementioned storage medium includes: various media that can store program codes, such as a USB flash drive, a mobile hard disk, a read-only memory (ROM), a random access memory (RAM), a magnetic disk or an optical disk.

[0138] In one embodiment of the present invention, a computer program product is provided, wherein the computer program product includes a computer program stored on a non-transitory computer-readable storage medium, wherein the computer program includes program instructions. When the program instructions are executed by a computer, the computer can execute the methods provided by the above-mentioned method embodiments.

[0139] In one embodiment of the present invention, a non-transitory computer-readable storage medium is provided, wherein the non-transitory computer-readable storage medium stores server instructions, and the computer instructions enable a computer to execute the methods provided in the above embodiments.

[0140] The above embodiment provides a computer-readable storage medium, whose implementation principle and technical effect are similar to those of the above method embodiment, and will not be repeated here.

[0141] The present application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0142] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0143] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0144] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for simulating the forming of an uncured fiber reinforced metal laminate, characterized in that: include: According to the key variables of the input uncured fiber, the strain direction of the fiber is obtained after preprocessing the key variables, and the tensile and shear strains of the fiber are obtained based on the strain direction; Based on the tensile and shear strains of the fiber, the current shear angle, shear modulus and Young's modulus are calculated, and the nominal stress of the fiber is updated; The damage criterion coefficient is calculated according to the updated nominal stress of the fiber, and wrinkling is judged by the shear angle. After wrinkling judgment, damage judgment and evolution are performed through the damage criterion coefficient.

2. The uncured fiber reinforced metal laminate forming simulation method according to claim 1, characterized in that: The key variables of unsolidified fibers include: the main yarn direction, the polynomial coefficient of Young's modulus, the polynomial coefficient of shear behavior, the tensile strength of the fiber bundle, the compressive strength of the fiber bundle, the operating factor β for the calculation of the fiber damage variable; the critical braiding angle γ for the set fabric wrinkling w .

3. The uncured fiber reinforced metal laminate forming simulation method according to claim 1, characterized in that: The fiber strain direction is obtained after preprocessing the key variables, including: The deformation gradient tensor and stretch tensor are established based on key variables to characterize the extension or compression characteristics of the fiber material under external loads; Based on the deformation gradient tensor and the stretch tensor, the rotation tensor and the local coordinate tensor of the rotation tensor are calculated to reflect the spatial rotation of the fiber direction caused by the deformation, and the rotation tensor is converted into the local coordinate system; Based on the rotation tensor and its local coordinate tensor, the current fiber orientation is determined in real time to capture the changes in fiber braiding angle during the forming process; The transformation matrix is ​​generated through the relative relationship between the fiber direction and the local coordinate system to convert the global loading condition into the stress state in the fiber direction; The fiber strain direction is calculated based on the transformation matrix to obtain the tensile and shear strains.

4. The uncured fiber reinforced metal laminate forming simulation method according to claim 1, characterized in that: Based on the tensile and shear strains of the fiber, the current shear angle, shear modulus and Young's modulus are calculated, and the nominal stress of the fiber is updated, including: s f，t+1 =s f，t +dε f ·C Among them, σ f，t+1 is the nominal stress of the updated fiber; C is the current stiffness matrix, which consists of the shear modulus and Young's modulus.

5. The uncured fiber reinforced metal laminate forming simulation method according to claim 1, characterized in that: Wrinkle judgment by shear angle, including: Determine whether the absolute value of the shear angle exceeds the critical braiding angle γ for fabric wrinkling w If it exceeds, wrinkling is determined to have occurred, otherwise wrinkling has not occurred.

6. The uncured fiber reinforced metal laminate forming simulation method according to claim 1, characterized in that: Damage judgment and evolution: based on the 3D-Hashin criterion, the damage of the material during the forming process is judged and evolved, including: By calculating different damage modes in the fiber direction, the 3D-Hashin damage criterion coefficient f is calculated to determine the damage variable By damage variable Identify possible failure scenarios; When the damage state is determined, the stress in the fiber direction, the matrix stress and the corresponding state variables are gradually updated to reflect the damage accumulation and performance degradation of the material during loading and simulate the damage evolution process.

7. The uncured fiber reinforced metal laminate forming simulation method according to claim 6, characterized in that: Damage variables for: Where e is the natural logarithm.

8. An uncured fiber reinforced metal laminate forming simulation system, characterized in that: include: A preprocessing module, according to the key variables of the input uncured fiber, preprocesses the key variables to obtain the fiber strain direction, and obtains the tensile and shear strains of the fiber based on the strain direction; The stress calculation update module calculates the current shear angle, shear modulus and Young's modulus based on the tensile and shear strains of the fiber, and updates the nominal stress of the fiber; The judgment evolution module calculates the damage criterion coefficient according to the updated nominal stress of the fiber, makes wrinkle judgment through the shear angle, and after wrinkle judgment, makes damage judgment and evolution through the damage criterion coefficient.

9. A computer-readable storage medium storing one or more programs, characterized in that: The one or more programs include instructions, which, when executed by a computing device, cause the computing device to perform any one of the methods of claims 1 to 7.

10. A computing device, characterized in that: include: One or more processors, a memory, and one or more programs, wherein the one or more programs are stored in the memory and configured to be executed by the one or more processors, and the one or more programs include instructions for executing any one of the methods described in claims 1 to 7.

Images