A composite material damage analysis method, device, equipment and storage medium

Abstract

The application discloses a composite material damage analysis method, device, equipment and storage medium. The method comprises the following steps: calculating the equivalent strain of the axial direction of the fiber material based on the axial strain tensor; calculating the equivalent strain of the transverse direction of the matrix material based on other strain tensors except the transverse strain tensor; calculating the axial damage variable of the fiber material based on the equivalent strain of the axial direction of the fiber material; calculating the transverse damage variable of the matrix material based on the equivalent strain of the transverse direction of the matrix material; and calculating the stress value representing the damage degree of the composite material based on the axial damage variable of the fiber material and the transverse damage variable of the matrix material. The method is suitable for damage analysis of the composite material.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of material analysis, and in particular to a composite material damage analysis method, device, equipment and storage medium. BACKGROUND

[0002] The fracture phenomenon of solid materials is derived from the evolution of small-scale micro-cracks, micro-cracks and micro-cavities in space. When these defects develop to a certain extent under the driving of deformation, macroscopic visible fracture occurs, so it is necessary to analyze the state of these small-scale defects. In the prior art, the state of small-scale defects is converted into a damage state variable through statistics, and the evolution dynamics of the damage state variable is used to represent the fracture failure process of the material.

[0003] With the rapid development of the material field, composite materials are widely used due to their advantages. Composite materials are a typical anisotropic material because their stiffness and strength in the axial and transverse directions are different, so the existing isotropic continuous damage mechanics model cannot be applied to the damage analysis of composite materials. SUMMARY

[0004] Based on the above problems, the present application provides a composite material damage analysis method, device, equipment and storage medium, which is suitable for damage analysis of composite materials.

[0005] The embodiments of the present application disclose the following technical solutions:

[0006] In a first aspect, the present application provides a composite material damage analysis method, which comprises:

[0007] Based on the axial strain tensor, the equivalent strain in the axial direction of the fiber material is calculated; based on other strain tensors except the transverse strain tensor, the equivalent strain in the transverse direction of the matrix material is calculated;

[0008] Based on the equivalent strain in the axial direction of the fiber material, the axial damage variable of the fiber material is calculated; based on the equivalent strain in the transverse direction of the matrix material, the transverse damage variable of the matrix material is calculated;

[0009] Based on the axial damage variable of the fiber material and the transverse damage variable of the matrix material, a stress value representing the damage degree of the composite material is calculated.

[0010] Optionally, the calculation of the axial damage variable of the fiber material based on the equivalent strain in the axial direction of the fiber material; and the calculation of the transverse damage variable of the matrix material based on the equivalent strain in the transverse direction of the matrix material, comprises:

[0011] Based on the equivalent strain in the axial direction of the fiber material, the axial damage threshold of the fiber material is obtained;

[0012] calculating an axial damage variable of the fiber material based on the axial damage threshold value of the fiber material, the damage threshold value of the damage initial state and the damage threshold value of the complete damage state;

[0013] obtaining a transverse damage threshold value of the matrix material based on the equivalent strain of the transverse direction of the matrix material;

[0014] calculating a transverse damage variable of the matrix material based on the transverse damage threshold value of the matrix material, the damage threshold value of the damage initial state and the damage threshold value of the complete damage state.

[0015] Optionally, the stress value representing the damage degree of the composite material is calculated based on the axial damage variable of the fiber material and the transverse damage variable of the matrix material, comprising:

[0016] calculating a stiffness degradation of the composite material based on the axial damage variable of the fiber material and the transverse damage variable of the matrix material;

[0017] calculating the stress value representing the damage degree of the composite material based on the ratio of the stiffness degradation of the composite material to the strain tensor other than the transverse strain tensor.

[0018] Optionally, the equivalent strain of the transverse direction of the matrix material is calculated based on the strain tensor other than the transverse strain tensor, comprising:

[0019] calculating a plurality of eigen-strain values of the strain tensor based on the strain tensor other than the transverse strain tensor;

[0020] calculating the equivalent strain of the transverse direction of the matrix material according to the plurality of eigen-strain values.

[0021] Optionally, the axial damage threshold value of the fiber material is obtained based on the equivalent strain of the axial direction of the fiber material, and the transverse damage threshold value of the matrix material is obtained based on the equivalent strain of the transverse direction of the matrix material, comprising:

[0022] comparing the equivalent strain of the axial direction of the fiber material with the size of the axial damage threshold value of the fiber material at the last time, and taking the larger value as the axial damage threshold value of the fiber material;

[0023] comparing the equivalent strain of the axial direction of the matrix material with the size of the transverse damage threshold value of the matrix material at the last time, and taking the larger value as the transverse damage threshold value of the matrix material.

[0024] In a second aspect, the application provides a composite material damage analysis device, which comprises:

[0025] The first computing module, the second computing module, and the third computing module;

[0026] The first computing module comprises a first computing unit and a second computing unit; the first computing unit is configured to calculate an equivalent strain in the axial direction of the fiber material based on the axial strain tensor; and the second computing unit is configured to calculate an equivalent strain in the transverse direction of the matrix material based on a strain tensor other than the transverse strain tensor.

[0027] The second computing module comprises a third computing unit and a fourth computing unit; the third computing unit is configured to calculate an axial damage variable of the fiber material based on the equivalent strain in the axial direction of the fiber material; and the fourth computing unit is configured to calculate a transverse damage variable of the matrix material based on the equivalent strain in the transverse direction of the matrix material.

[0028] The third computing module is configured to calculate a stress value representing a damage degree of the composite material based on the axial damage variable of the fiber material and the transverse damage variable of the matrix material.

[0029] Optionally, the second computing module is specifically configured to:

[0030] The third computing unit in the second computing module is specifically configured to obtain an axial damage threshold of the fiber material based on the equivalent strain in the axial direction of the fiber material, and calculate the axial damage variable of the fiber material based on the axial damage threshold of the fiber material, a damage threshold of a damage initial state, and a damage threshold of a complete damage state.

[0031] The fourth computing unit in the second computing module is specifically configured to obtain a transverse damage threshold of the matrix material based on the equivalent strain in the transverse direction of the matrix material, and calculate the transverse damage variable of the matrix material based on the transverse damage threshold of the matrix material, a damage threshold of a damage initial state, and a damage threshold of a complete damage state.

[0032] Optionally, the third computing module is specifically configured to:

[0033] Calculate a stiffness degradation of the composite material based on the axial damage variable of the fiber material and the transverse damage variable of the matrix material.

[0034] Calculate a stress value representing a damage degree of the composite material based on a ratio of the stiffness degradation of the composite material to a strain tensor other than the transverse strain tensor.

[0035] Optionally, the first computing module is specifically configured to:

[0036] Calculate a plurality of eigenstrain values of the strain tensor based on a strain tensor other than the transverse strain tensor.

[0037] According to the multiple intrinsic strain values, the equivalent strain of the matrix material in the transverse direction is calculated.

[0038] In a third aspect, the present application provides a computer device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, when the processor executes the computer program, the damage analysis method of the composite material in any one of the first aspect is realized.

[0039] In a fourth aspect, the present application provides a computer readable storage medium, the computer readable storage medium stores instructions, when the instructions run on the terminal device, the terminal device executes the damage analysis method of the composite material in any one of the first aspect.

[0040] The present application provides a damage analysis method of a composite material, first, based on the axial strain tensor, the equivalent strain of the fiber material in the axial direction is calculated; based on other strain tensors except the transverse strain tensor, the equivalent strain of the matrix material in the transverse direction is calculated; secondly, based on the equivalent strain of the fiber material in the axial direction, the axial damage variable of the fiber material is calculated; based on the equivalent strain of the matrix material in the transverse direction, the transverse damage variable of the matrix material is calculated; finally, based on the axial damage variable of the fiber material and the transverse damage variable of the matrix material, the stress value representing the damage degree of the composite material is calculated.

[0041] Different damage state variables are used to describe the damage state of the fiber material and the matrix material respectively, so that the damage evolution of the fiber material and the matrix material can be independently characterized, thereby providing numerical support for the design and optimization of the composite material structure. BRIEF DESCRIPTION OF DRAWINGS

[0042] In order to more clearly illustrate the technical solutions of the embodiments of the present application or the prior art, the drawings needed in the embodiment or prior art description will be briefly introduced below. Obviously, the drawings in the following description are only some embodiments of the present application, and those skilled in the art can also obtain other drawings according to these drawings without creative labor.

[0043] Figure 1 A flowchart of a damage analysis method of a composite material provided by an embodiment of the present application;

[0044] Figure 2 A flowchart of another damage analysis method of a composite material provided by an embodiment of the present application;

[0045] Figure 3 An example diagram of axial tension of a fiber material provided by an embodiment of the present application;

[0046] Figure 4 Fig. 1 is a structural schematic diagram of a damage analysis device for a material according to an embodiment of the present application. DETAILED DESCRIPTION

[0047] The terms "first", "second", and "third" and the like in the description and in the claims of the present application and the drawings are used for distinguishing between similar objects and not necessarily for describing a particular sequential or chronological order. It is to be understood that the terms so used are interchangeable under appropriate circumstances and are not intended to limit the application unless otherwise explicitly indicated by context.

[0048] In the embodiments of the present application, the words "exemplary" and "for example" are used to mean serving as an example, instance, or illustration, at 99 least with respect to the specific embodiments described and not necessarily with respect to other embodiments. In other words, the words "exemplary" and "for example" are used herein to mean "one example, among others." As used herein, "for example" means "for the purpose of example, or to serve as an example, instance, or illustration." Any embodiment or design described herein as "exemplary" or "for example" is not necessarily to be construed as preferred or advantageous over other embodiments or designs.

[0049] Traditional strength criterion analysis cannot characterize the stiffness degradation process of the material until it is destroyed, so a model of isotropic continuum damage mechanics is proposed. However, the model cannot handle the special case of composite materials with different strengths in different directions (isotropic). After research, the present application provides a continuum damage mechanics model based on isotropy, which analyzes the stiffness degradation of composite materials under the framework of finite element analysis method, and provides data calculation basis for composite structure design optimization.

[0050] In view of this, the present application provides a damage analysis method for a composite material. First, based on an axial strain tensor, an equivalent strain in an axial direction of a fiber material is calculated; based on other strain tensors except a transverse strain tensor, an equivalent strain in a transverse direction of a matrix material is calculated; second, based on the equivalent strain in the axial direction of the fiber material, an axial damage variable of the fiber material is calculated; based on the equivalent strain in the transverse direction of the matrix material, a transverse damage variable of the matrix material is calculated; finally, based on the axial damage variable of the fiber material and the transverse damage variable of the matrix material, a stress value representing a damage degree of the composite material is calculated.

[0051] In order for those skilled in the art to better understand the technical scheme of the present application, the technical scheme of the embodiments of the present application will be described clearly and completely below in conjunction with the drawings in the embodiments of the present application. Obviously, the described embodiments are only some of the embodiments of the present application, not all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative labor fall within the scope of protection of the present application.

[0052] Referring to Figure 1 Fig. 1 is a structural schematic diagram of a damage analysis device for a material according to an embodiment of the present application.

[0053] As Figure 1As shown, the method comprises:

[0054] S101: Based on the axial strain tensor, the equivalent strain of the axial direction of the fiber material is calculated; based on other strain tensors except the transverse strain tensor, the equivalent strain of the transverse direction of the matrix material is calculated.

[0055] The composite material in the present application can be a fiber-reinforced composite material.

[0056] The fiber-reinforced composite material is mainly composed of fiber material and matrix material, wherein the fiber material includes but is not limited to glass fiber, carbon fiber, boron fiber, aramid fiber, silicon carbide fiber, asbestos fiber, etc.; the matrix material is divided into two categories of metal and non-metal, the commonly used metal matrix is aluminum, magnesium, copper, titanium and their alloys, and the non-metal matrix mainly includes synthetic resin, rubber, ceramic, graphite, carbon, etc.

[0057] The axial direction of the fiber material can be parallel to the axial direction; the transverse direction of the matrix material can be the width direction.

[0058] Specifically, the equivalent strain of the axial direction of the fiber material can be calculated by the following formula:

[0059] E eqv,1 ≡<E 11 >0

[0060] Wherein ε eqv,1 represents the axial equivalent strain of the fiber material, and ε 11 represents the axial strain tensor, and here <x> C is defined as:

[0061]

[0062] Optionally, a plurality of eigen-strain values of the strain tensor are calculated based on a strain tensor other than the transverse strain tensor; and the equivalent strain of the matrix material in the transverse direction is calculated according to the plurality of eigen-strain values.

[0063] Specifically, the strain tensor other than the transverse strain tensor can be calculated by the following formula:

[0064]

[0065] Since the above matrix is a three-order matrix, the eigenvalues of the above matrix are calculated to obtain three eigen-strain values, which are respectively represented by and The equivalent strain of the matrix material in the transverse direction is calculated by the eigen-strain values of the strain tensor by the following formula:

[0066]

[0067] Through the above process, the equivalent strain of the fiber material in the axial direction and the equivalent strain of the matrix material in the transverse direction can be calculated respectively.

[0068] S102: Based on the equivalent strain of the fiber material in the axial direction, the axial damage variable of the fiber material is calculated; and based on the equivalent strain of the matrix material in the transverse direction, the transverse damage variable of the matrix material is calculated.

[0069] Specifically, based on the equivalent strain of the fiber material in the axial direction, the axial damage threshold of the fiber material is obtained; and based on the axial damage threshold of the fiber material, the damage threshold of the initial damage state and the damage threshold of the complete damage state, the axial damage variable of the fiber material is calculated.

[0070] Based on the equivalent strain of the matrix material in the transverse direction, the transverse damage threshold of the matrix material is obtained; and based on the transverse damage threshold of the matrix material, the damage threshold of the initial damage state and the damage threshold of the complete damage state, the transverse damage variable of the matrix material is calculated.

[0071] The damage threshold of the initial damage state and the damage threshold of the complete damage state of the fiber material and the matrix material can be set in advance according to conditions, that is, the damage threshold of the initial damage state and the damage threshold of the complete damage state of the fiber material and the matrix material are known quantities in the calculation process.

[0072] t n+1 the axial damage threshold of the fiber material at the moment 1，n+1 and the transverse damage threshold r of the matrix material 2，n+1 ε can be obtained by calculating S101 eqv，1 and ε eqv，2 With t n Axial damage threshold r of fiber material at time t 1，n and the transverse damage threshold r of the matrix material 2，n The larger of the two values ​​is selected as the axial damage threshold and the transverse damage threshold at the current moment. Specifically, this can be calculated using the following formula:

[0073]

[0074]

[0075] The axial damage threshold of the fiber material and the transverse damage threshold of the matrix material can be obtained through a piecewise function g. i (r i Perform the corresponding calculations, where g i (r i ) requires monotonicity and g i (r i )∈[0,1]:

[0076]

[0077] When i = 1, g i (r i ) represents the axial damage threshold of the fiber material, r i，0 The damage threshold r represents the initial stage of damage to the fiber material. i，1 This represents the damage threshold corresponding to the moment when the fiber material is in a completely damaged state.

[0078] When i = 2, g i (r i ) represents the axial damage threshold of the matrix material, r i，0 The damage threshold r represents the initial stage of damage to the matrix material. i，1 This represents the damage threshold corresponding to the moment when the matrix material is in a fully damaged state.

[0079] The axial damage variable of the fiber material and the transverse damage variable of the matrix material can be calculated based on the damage threshold of the fiber material and the damage threshold of the matrix material using the following formula:

[0080] ω1=g(r1)

[0081] ω2=g(r2)

[0082] The damage variable of the fiber material in the axial direction and the damage variable of the matrix material in the transverse direction can be calculated respectively by S102.

[0083] S103: calculating the stress value representing the damage degree of the composite material based on the axial damage variable of the fiber material and the transverse damage variable of the matrix material.

[0084] Specifically, the stiffness reduction of the composite material is calculated based on the axial damage variable of the fiber material and the transverse damage variable of the matrix material; and the stress value representing the damage degree of the composite material is calculated based on the ratio of the stiffness reduction of the composite material to the strain tensor other than the transverse strain tensor.

[0085] The damage variable of the fiber material in the axial direction and the damage variable of the matrix material in the transverse direction calculated by S102 can be used to calculate the stiffness reduction of the composite material, which is specifically described as follows:

[0086]

[0087] Wherein,

[0088]

[0089]

[0090] In the above formula, the 1 direction is the fiber axial direction, and the 2 and 3 directions are the fiber transverse directions. That is, in this model, the effective bearing area of different 1 directions and 2 and 3 directions is considered, which is reasonable for fiber-reinforced composite materials.

[0091] When an object is deformed due to external factors (force, humidity, temperature field change, etc.), internal forces are generated between different parts of the object to resist the action of the external factors and try to restore the object from the deformed position to the pre-deformed position. The internal force per unit area at a certain point of the cross section under consideration is called stress.

[0092] The stress can be calculated by the following formula:

[0093] σ = L s : s

[0094] First, the equivalent strain of the fiber material in the axial direction is calculated based on the axial strain tensor; the equivalent strain of the matrix material in the transverse direction is calculated based on the strain tensor other than the transverse strain tensor; second, the axial damage variable of the fiber material is calculated based on the equivalent strain of the fiber material in the axial direction; the transverse damage variable of the matrix material is calculated based on the equivalent strain of the matrix material in the transverse direction; finally, the stress value representing the damage degree of the composite material is calculated based on the axial damage variable of the fiber material and the transverse damage variable of the matrix material.

[0095] Different damage state variables are used to describe the damage state of the fiber material and the matrix material respectively, so that the damage evolution of the fiber material and the matrix material can be independently characterized, thereby providing numerical support for the design and optimization of the composite material structure.

[0096] Referring to Figure 2 , the figure is a flowchart of another composite material damage analysis method provided by the embodiment of the application.

[0097] As Figure 2 shown, the method comprises:

[0098] S201: based on the axial strain tensor, the equivalent strain of the axial direction of the fiber material is calculated; based on other strain tensors except the transverse strain tensor, the equivalent strain of the transverse direction of the matrix material is calculated.

[0099] Specifically, the equivalent strain of the axial direction of the fiber material can be calculated by the following formula:

[0100] E eqv，1 ≡<E 11 >0

[0101] Wherein ε eqv，1 represents the axial equivalent strain of the fiber material, ε 11 represents the axial strain tensor, and here <x> C is defined as:

[0102]

[0103] Specifically, other strain tensors except the transverse strain tensor can be calculated by the following formula:

[0104]

[0105] Since the above matrix is a three-order matrix, the eigenvalues of the above matrix are calculated, and three eigen-strain values are obtained, which are represented by and respectively. The equivalent strain of the matrix material in the transverse direction can be calculated by the following formula:

[0106]

[0107] The eigen-strain value with the largest absolute value among the three eigen-strain values is taken as the equivalent strain of the matrix material in the transverse direction.

[0108] S202: Based on the equivalent strain of the fiber material in the axial direction, the axial damage variable of the fiber material is calculated; and based on the equivalent strain of the matrix material in the transverse direction, the transverse damage variable of the matrix material is calculated.

[0109] S203: Based on the axial damage variable of the fiber material and the transverse damage variable of the matrix material, the stress value representing the damage degree of the composite material is calculated.

[0110] Referring to Figure 4 , which is a structural schematic diagram of a damage analysis device of a composite material provided by an embodiment of the present application.

[0111] As shown in Figure 4 , the device comprises:

[0112] a first calculation module 401, a second calculation module 402, and a third calculation module 403;

[0113] The first calculation module 401 comprises a first calculation unit 4011 and a second calculation unit 4012; the first calculation unit 4011 is configured to calculate the equivalent strain of the fiber material in the axial direction based on the axial strain tensor; and the second calculation unit 4012 is configured to calculate the equivalent strain of the matrix material in the transverse direction based on other strain tensors except the transverse strain tensor;

[0114] The second calculation module 402 comprises a third calculation unit 4013 and a fourth calculation unit 4014; the third calculation unit is configured to calculate an axial damage variable of the fiber material based on the equivalent strain of the axial direction of the fiber material; and the fourth calculation unit 4014 is configured to calculate a transverse damage variable of the matrix material based on the equivalent strain of the transverse direction of the matrix material.

[0115] The third calculation module 403 is configured to calculate a stress value representing the damage degree of the composite material based on the axial damage variable of the fiber material and the transverse damage variable of the matrix material.

[0116] Optionally, the third calculation unit 4013 in the second calculation module 402 is specifically configured to obtain an axial damage threshold of the fiber material based on the equivalent strain of the axial direction of the fiber material; and calculate the axial damage variable of the fiber material based on the axial damage threshold of the fiber material, a damage threshold of a damage initial state and a damage threshold of a complete damage state.

[0117] The fourth calculation unit 4014 in the second calculation module 402 is specifically configured to obtain a transverse damage threshold of the matrix material based on the equivalent strain of the transverse direction of the matrix material; and calculate the transverse damage variable of the matrix material based on the transverse damage threshold of the matrix material, the damage threshold of the damage initial state and the damage threshold of the complete damage state.

[0118] Optionally, the third calculation module 403 is specifically configured to:

[0119] calculate a stiffness degradation of the composite material based on the axial damage variable of the fiber material and the transverse damage variable of the matrix material.

[0120] calculate the stress value representing the damage degree of the composite material based on the ratio of the stiffness degradation of the composite material to a strain tensor other than the transverse strain tensor.

[0121] Optionally, the first calculation module 401 is specifically configured to:

[0122] calculate a plurality of eigenstrain values of the strain tensor based on the strain tensor other than the transverse strain tensor; and calculate the equivalent strain of the transverse direction of the matrix material according to the plurality of eigenstrain values.

[0123] Figure 3 For the calculation example corresponding to the fiber axial tension, it can be known from the figure that the above method can adopt different damage state variables to respectively describe the damage states of the fiber material and the matrix material, so that the damage evolutions of the fiber material and the matrix material can be independently represented, thereby providing numerical support for the design and optimization of the composite material structure.

[0124] The embodiment of the present application provides a computer readable storage medium, which stores a computer program, and the program is executed by a processor to realize the damage analysis method of the composite material.

[0125] In practical application, the computer readable storage medium can adopt any combination of one or more computer readable media. The computer readable medium can be a computer readable signal medium or a computer readable storage medium. The computer readable storage medium can be, for example but not limited to, an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, device or apparatus, or any combination of the above. More specific examples (non-exhaustive list) of the computer readable storage medium include an electrical connection having one or more wires, a portable computer diskette, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the above. In this embodiment, the computer readable storage medium can be any tangible medium containing or storing a program that can be used by or in connection with an instruction execution system, apparatus or device.

[0126] The computer readable signal medium can include a data signal propagating in baseband or propagating as a carrier wave in a propagated data signal, in which computer readable program code is embodied. Such propagated data signal can take a variety of forms, including but not limited to electro-magnetic, optical, or any suitable combination thereof. The computer readable signal medium can also be any computer readable medium that is not a storage medium and that can communicate, propagate or transport program for use by or in connection with an instruction execution system, apparatus, or device.

[0127] The program code contained in the computer readable medium can be transmitted by any suitable medium, including but not limited to wireless, wire line, optical cable, RF, etc., or any suitable combination thereof.

[0128] Computer program code for carrying out operations of the present application can be written in any combination of one or more programming languages, including an object oriented programming language such as Java, Smalltalk, C++ or the like, and conventional procedural programming languages, such as the "C" programming language or similar programming languages. The program code can execute entirely on the user's computer, partly on the user's computer, as a stand-alone software package, partly on the user's computer and partly on a remote computer or entirely on the remote computer or server. In the latter scenario, the remote computer can be connected to the user's computer through any type of network, including a local area network (LAN) or a wide area network (WAN), or the connection can be made to an external computer (for example, through the Internet using an Internet Service Provider).

[0129] It should be noted that each of the embodiments described in the specification of the present application adopts a progressive manner, and the same or similar parts between each embodiment can be referred to each other. Each embodiment focuses on the difference from other embodiments. Especially, the device is described simply because it is basically similar to the method embodiment. The relevant part can be referred to the part of the method embodiment. The device described above is only schematic and the various components illustrated in the device can be implemented or not implemented physically, that is, they can be located in one place or distributed on multiple network units. Part or all of the modules can be selected according to the actual needs to achieve the purpose of the embodiment. Those skilled in the art can understand and implement it without creative labor.

[0130] The above description is only one specific embodiment of the present application, but the protection scope of the present application is not limited to this. Any changes or replacements within the technical scope disclosed in the present application can be easily thought by those skilled in the art without creative labor, which should be covered in the protection scope of the present application. Therefore, the protection scope of the present application should be subject to the protection scope of the claims.< / x> < / x>

Claims

1. A damage analysis method for composite materials, characterized in that, The method includes: Based on the axial strain tensor, the equivalent strain in the axial direction of the fiber material is calculated; based on strain tensors other than the transverse strain tensor, the equivalent strain in the transverse direction of the matrix material is calculated. Based on the equivalent strain of the fiber material in the axial direction, the axial damage variable of the fiber material is calculated; based on the equivalent strain of the matrix material in the transverse direction, the transverse damage variable of the matrix material is calculated. Based on the axial damage variable of the fiber material and the transverse damage variable of the matrix material, the stress value characterizing the damage degree of the composite material is calculated. The step of calculating the equivalent strain in the transverse direction of the matrix material based on strain tensors other than the transverse strain tensor includes: calculating multiple intrinsic strain values ​​of the strain tensor based on strain tensors other than the transverse strain tensor; and calculating the equivalent strain in the transverse direction of the matrix material based on the multiple intrinsic strain values. Specifically, the calculation of axial damage variables of the fiber material based on the equivalent strain in the axial direction and the calculation of transverse damage variables of the matrix material based on the equivalent strain in the transverse direction include: obtaining an axial damage threshold of the fiber material based on the equivalent strain in the axial direction; calculating an axial damage variable of the fiber material based on the axial damage threshold, the damage threshold at the damage initiation state, and the damage threshold at the complete damage state; obtaining a transverse damage threshold of the matrix material based on the equivalent strain in the transverse direction; and calculating a transverse damage variable of the matrix material based on the transverse damage threshold, the damage threshold at the damage initiation state, and the damage threshold at the complete damage state. The step of calculating the stress value characterizing the degree of damage to the composite material based on the axial damage variable of the fiber material and the transverse damage variable of the matrix material includes: calculating the stiffness loss of the composite material based on the axial damage variable of the fiber material and the transverse damage variable of the matrix material; and calculating the stress value characterizing the degree of damage to the composite material based on the ratio of the stiffness loss of the composite material to other strain tensors besides the transverse strain tensor.

2. The method as described in claim 1, characterized in that, The axial damage threshold of the fiber material is obtained based on the equivalent strain in the axial direction of the fiber material, and the lateral damage threshold of the matrix material is obtained based on the equivalent strain in the lateral direction of the matrix material, including: The equivalent strain in the axial direction of the fiber material is compared with the magnitude of the axial damage threshold of the fiber material at the previous moment, and the larger value of the two is taken as the axial damage threshold of the fiber material. The equivalent strain in the axial direction of the matrix material is compared with the magnitude of the transverse damage threshold of the matrix material at the previous moment, and the larger value is taken as the transverse damage threshold of the matrix material.

3. A damage analysis apparatus for composite materials, the apparatus comprising: The first calculation module, the second calculation module, and the third calculation module; The first computing module includes a first computing unit and a second computing unit; The first calculation unit is used to calculate the equivalent strain of the fiber material in the axial direction based on the axial strain tensor; the second calculation unit is used to calculate the equivalent strain of the matrix material in the transverse direction based on strain tensors other than the transverse strain tensor. The second calculation module includes a third calculation unit and a fourth calculation unit; the third calculation unit is used to calculate the axial damage variable of the fiber material based on the equivalent strain in the axial direction of the fiber material. The fourth calculation unit is used to calculate the lateral damage variable of the matrix material based on the equivalent strain in the lateral direction of the matrix material. The third calculation module is used to calculate the stress value characterizing the degree of damage of the composite material based on the axial damage variable of the fiber material and the transverse damage variable of the matrix material. Specifically, the second calculation unit is used to calculate multiple intrinsic strain values ​​of the strain tensor based on strain tensors other than the transverse strain tensor; and to calculate the equivalent strain in the transverse direction of the matrix material based on the multiple intrinsic strain values. Specifically, the third calculation unit is used to obtain the axial damage threshold of the fiber material based on the equivalent strain in the axial direction of the fiber material; and to calculate the axial damage variable of the fiber material based on the axial damage threshold, the damage threshold at the damage initiation state, and the damage threshold at the complete damage state. The fourth calculation unit is specifically used to obtain the transverse damage threshold of the matrix material based on the equivalent strain in the transverse direction of the matrix material; and to calculate the transverse damage variable of the matrix material based on the transverse damage threshold, the damage threshold at the damage initiation state, and the damage threshold at the complete damage state. Specifically, the third calculation module is used to calculate the stiffness loss of the composite material based on the axial damage variable of the fiber material and the transverse damage variable of the matrix material; and to calculate the stress value characterizing the degree of damage of the composite material based on the ratio of the stiffness loss of the composite material to other strain tensors except the transverse strain tensor.

4. A computer device, characterized in that, include: A memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor, when executing the computer program, implements the damage analysis method for composite materials as described in claim 1 or 2.

5. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores instructions that, when executed on a terminal device, cause the terminal device to perform the damage analysis method for composite materials as described in claim 1 or 2.

Images