Bilinear constitutive parameter engineering establishment method, system, equipment and medium

Abstract

The invention belongs to the field of aero-engines, relates to a bilinear constitutive curve design technology, and discloses a bilinear constitutive parameter engineering establishment method, system and equipment and a medium. The method comprises the steps that a part made of a to-be-calibrated material is adopted, and bilinear constitutive parameter initial values are obtained according to tensile mechanical property data; according to the bilinear constitutive parameter initial value, elastic-plastic finite element analysis is conducted on the part, and first engineering stress and first engineering strain of the part under the working load and engineering residual strain of the part after the working load is unloaded are obtained; carrying out a working load loading and unloading test to obtain the real residual strain of the part; and evaluating and correcting the initial value of the bilinear constitutive parameter through the engineering residual strain and the real residual strain. According to the method, the constitutive curve capable of meeting engineering use is obtained through the tensile mechanical property data, data support is provided for strength design of engine parts, and the method has the advantages of being few in needed parameters, easy to operate and high in engineering operability.

Description

Technical Field

[0001] The invention belongs to the field of aeroengines, relates to a bilinear constitutive curve design technology, and specifically relates to a bilinear constitutive parameter engineering establishment method, system, equipment and medium. Background Art

[0002] Aircraft engines have many parts and many types of materials. Obtaining the mechanical properties, that is, the linear constitutive curve, of each material requires a lot of manpower, material resources, and time costs. In addition, with the gradual improvement of material processing, it is very difficult to obtain comprehensive mechanical properties data for each material in a short period of time.

[0003] At present, the material manual provides relatively comprehensive mechanical property data for commonly used materials, which can basically meet the use requirements. However, there are few performance data for most materials, especially non-critical materials. Usually, only tensile mechanical property data such as tensile strength, material yield strength, elastic modulus and cross-sectional shrinkage are provided. However, when designing parts, it is necessary to conduct elastic-plastic analysis on the parts to determine the residual deformation of the parts. Only tensile mechanical property data cannot meet its design requirements.

[0004] In view of this, it becomes more urgent to use the existing basic mechanical properties data to obtain constitutive parameters that meet engineering needs. Summary of the invention

[0005] In order to solve the technical problem that some materials in the existing material manuals only have tensile mechanical properties data, which does not meet the requirements of part elastic-plastic analysis, the present invention discloses a bilinear constitutive parameter engineering establishment method, which includes the following steps: S1. For a part made of a material to be calibrated, according to the elastic modulus, material section shrinkage, material yield strength and material tensile strength of the material to be calibrated, the initial values ​​of the bilinear constitutive parameters of the material to be calibrated are obtained, wherein the initial values ​​of the bilinear constitutive parameters include the initial value of the yield stress and the initial value of the tangent modulus; S2. performing elastic-plastic finite element analysis on the part according to the initial values ​​of the bilinear constitutive parameters to obtain the first engineering stress and the first engineering strain of the part under the working load, and the engineering residual strain of the part after the working load is unloaded; S3, performing a working load loading and unloading test on the part to obtain a true residual strain of the part; S4. By using the engineering residual strain and the true residual strain, the initial value of the bilinear constitutive parameter is evaluated and corrected according to the strain threshold to obtain the bilinear constitutive parameter.

[0006] Furthermore, in step S1, the initial values ​​of the bilinear constitutive parameters of the material to be calibrated are obtained according to the elastic modulus, the material cross-sectional shrinkage, the material yield strength and the material tensile strength of the material to be calibrated, including: S11, given the residual deformation of the part after load unloading, calculating the engineering stress and engineering strain of the deformation point corresponding to the residual deformation by using the material yield strength; S12, calculating the true stress and true strain of the deformation point through the engineering stress and the engineering strain; S13, obtaining the first true strain and the first true stress of the engineering fracture point according to the material cross-sectional shrinkage rate and the material tensile strength; S14. Obtain initial values ​​of bilinear constitutive parameters according to the elastic modulus, the true stress, the true strain, the first true strain and the first true stress.

[0007] Furthermore, in step S14, obtaining the initial value of the bilinear constitutive parameter according to the elastic modulus, the true stress, the true strain, the first true strain and the first true stress includes: S141, obtaining an initial value of the tangent modulus of the material to be calibrated according to the true stress, the true strain, the first true strain and the first true stress; S142, according to the initial value of the tangent modulus and the elastic modulus The first true strain and the first true stress are used to obtain an initial value of the yield stress of the material to be calibrated.

[0008] Furthermore, in step S3, a working load loading and unloading test is performed on the part to obtain the real residual strain of the part, including: S31, measuring a first dimension of the part before loading a working load; S32, performing a working load loading and unloading test on the part, and measuring a second dimension of the part after the working load is unloaded; S33. Calculate the true residual strain of the part using the first dimension and the second dimension.

[0009] Further, in step S4, the bilinear constitutive parameter is obtained by evaluating and correcting the initial value of the bilinear constitutive parameter according to the strain threshold through the engineering residual strain and the true residual strain, including: S41, using the engineering residual strain and the true residual strain, according to the formula Calculate the ratio, where is the true residual strain, is the engineering residual strain; S42, when the ratio is less than or equal to a given strain threshold, judging that the initial value of the bilinear constitutive parameter is reasonable, and taking the initial value of the bilinear constitutive parameter as the final bilinear constitutive parameter; S43. When the ratio is greater than a given strain threshold, it is determined that the initial values ​​of the bilinear constitutive parameters are unreasonable, and the initial values ​​of the bilinear constitutive parameters are corrected.

[0010] Furthermore, in step S43, when the ratio is greater than a given strain threshold, it is determined that the initial value of the bilinear constitutive parameter is unreasonable, and the initial value of the bilinear constitutive parameter is corrected, including: S431, using the first engineering stress, the first engineering strain, the elastic modulus and the true residual strain to obtain a second true stress and a second true strain at a working load point of the part; S44. Calculate bilinear constitutive parameters based on the elastic modulus, the true stress, the true strain, the second true stress and the second true strain.

[0011] The embodiment of the present invention also provides a bilinear constitutive parameter engineering establishment system, which includes a bilinear constitutive parameter initial value calculation module, a simulation module, a loading and unloading test module, and an evaluation and correction module.

[0012] The bilinear constitutive parameter initial value calculation module is used to obtain the bilinear constitutive parameter initial values ​​of the material to be calibrated according to the elastic modulus, material section shrinkage, material yield strength and material tensile strength of the material to be calibrated, and the bilinear constitutive parameter initial values ​​include the yield stress initial value and the tangent modulus initial value; The simulation module is used to perform elastic-plastic finite element analysis on the part according to the initial value of the bilinear constitutive parameter, and obtain the first engineering stress and the first engineering strain of the part under the working load, and the engineering residual strain of the part after the working load is unloaded; The loading and unloading test module is used to perform a working load loading and unloading test on the part to obtain the real residual strain of the part; The evaluation and correction module is used to evaluate and correct the initial value of the bilinear constitutive parameter according to the strain threshold through the engineering residual strain and the true residual strain to obtain the bilinear constitutive parameter.

[0013] An embodiment of the present invention also provides a computer device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the computer program, it implements any of the above-mentioned bilinear constitutive parameter engineering establishment methods to solve the technical problem that some materials in the existing material manuals only have tensile mechanical property data, which does not meet the requirements of part elastic-plastic analysis.

[0014] An embodiment of the present invention also provides a computer-readable storage medium, which stores a computer program for executing any of the above-mentioned bilinear constitutive parameter engineering establishment methods to solve the technical problem that some materials in existing material manuals only have tensile mechanical properties data, which does not meet the requirements of part elastic-plastic analysis.

[0015] Compared with the prior art, the beneficial effects that can be achieved by at least one of the above technical solutions adopted in the embodiments of this specification include at least: the bilinear constitutive parameter engineering establishment method of the present invention first calculates an initial value of a bilinear constitutive parameter through the tensile mechanical property data in the material manual, and then performs the elastic-plastic loading and unloading of the working load through simulation technology and real experiments, respectively, and evaluates and corrects the initial value of the bilinear constitutive parameter to obtain the final bilinear constitutive parameter. The method of the present invention obtains a constitutive curve that can meet engineering use through tensile mechanical property data, provides data support for the strength design of engine parts, and the method has the advantages of few required parameters, easy operation and strong engineering operability. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0017] Figure 1 A flow chart of a method for establishing bilinear constitutive parameter engineering disclosed in an embodiment of the present invention; Figure 2 It is a schematic diagram of establishing a bilinear constitutive parameter curve through tensile mechanical property data of a material to be calibrated according to an embodiment of the present invention; Figure 3 The schematic diagram of the bilinear constitutive parameter engineering establishment system disclosed in the embodiment of the present invention; Among them, 301 is a bilinear constitutive parameter initial value calculation module; 302 is a simulation module; 303 is a loading and unloading test module; 304 is an evaluation and correction module. DETAILED DESCRIPTION

[0018] The embodiments of the present application are described in detail below with reference to the accompanying drawings.

[0019] The following describes the implementation methods of the present application through specific specific examples, and those skilled in the art can easily understand other advantages and effects of the present application from the content disclosed in this specification. Obviously, the described embodiments are only a part of the embodiments of the present application, rather than all the embodiments. The present application can also be implemented or applied through other different specific implementation methods, and the details in this specification can also be modified or changed in various ways based on different viewpoints and applications without departing from the spirit of the present application. It should be noted that the following embodiments and the features of the embodiments can be combined with each other in the absence of conflict. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in the field without making creative work belong to the scope of protection of the present application.

[0020] The embodiment of the present invention discloses a method for establishing bilinear constitutive parameters engineering, see Figure 1 and Figure 2 As shown, the method comprises the following steps: S1. For a part made of a material to be calibrated, the initial values ​​of the bilinear constitutive parameters of the material to be calibrated are obtained according to the elastic modulus E, the material cross-sectional shrinkage, the material yield strength and the material tensile strength of the material to be calibrated. The initial values ​​of the bilinear constitutive parameters include the initial value of the yield stress and the initial value of the tangent modulus. The tensile mechanical properties data of the material to be calibrated include the elastic modulus E, the material cross-sectional shrinkage, the material yield strength and the material tensile strength. , the material yield strength is and material tensile strength , all of which can be obtained from the materials manual.

[0021] S2. Performing elastic-plastic finite element analysis on the part according to the initial value of the bilinear constitutive parameter to obtain the first engineering stress of the part under the working load. and the first engineering strain , and the engineering residual strain of the part after unloading of the working load ; S3. Perform working load loading and unloading tests on the parts to obtain the true residual strain of the parts. ; S4, through the engineering residual strain and the true residual strain , the bilinear constitutive parameters are obtained by evaluating and correcting the initial values ​​of the bilinear constitutive parameters according to the strain threshold.

[0022] Furthermore, in step S1, according to the elastic modulus E and the cross-sectional shrinkage rate of the material to be calibrated , material yield strength and material tensile strength , obtaining the initial values ​​of the bilinear constitutive parameters of the material to be calibrated, including: S11. Given the residual deformation of the part after load unloading, calculate the engineering stress and engineering strain at the deformation point D corresponding to the residual deformation by using the material yield strength. .

[0023] In the specific implementation, see Figure 2 As shown in the figure, point D is the deformation point corresponding to the residual deformation of the material to be calibrated after unloading is 0.2% (that is, the length of OA is 0.2%), and the engineering stress at point D is for , engineering strain for , where the elastic modulus E is the slope of the straight line segment OC.

[0024] S12, calculating the true stress of the deformation point through the engineering stress and the engineering strain and true strain Specifically, the formula ; Formula for calculation.

[0025] S13, according to the material cross-sectional shrinkage rate and material tensile strength, obtain the first true strain of the engineering fracture point T and the first true stress Specifically, the formula and formula Perform calculations.

[0026] S14. Obtain initial values ​​of bilinear constitutive parameters according to the elastic modulus, the true stress, the true strain, the first true strain and the first true stress.

[0027] Furthermore, in step S14, obtaining the initial value of the bilinear constitutive parameter according to the elastic modulus, the true stress, the true strain, the first true strain and the first true stress includes: S141, obtaining an initial value of the tangent modulus of the material to be calibrated according to the true stress, the true strain, the first true strain and the first true stress Specifically, the formula Calculated.

[0028] S142, obtaining the initial value of the yield stress of the material to be calibrated according to the initial value of the tangent modulus, the elastic modulus, the first true stress and the first true strain Specifically, the initial value of the yield stress is also referred to as Figure 2The true stress at point G shown in can be expressed by the formula Calculated.

[0029] Furthermore, in the specific implementation of step S2, the material properties can usually be calculated using a bilinear isotropic hardening model. When performing simulation calculations, the initial values ​​of the bilinear constitutive parameters (yield stress and tangent modulus data) are input, and finite element calculations are first performed under the working load to obtain the first engineering stress of the part under the working load. The first engineering response Based on this calculation result, the working load is unloaded to obtain the engineering residual strain of the part after unloading. .

[0030] Furthermore, in step S3, a working load loading and unloading test is performed on the part to obtain the real residual strain of the part, including: S31, measuring the first dimension of the part before loading the working load ; S32, performing a working load loading and unloading test on the part, and measuring a second dimension of the part after unloading the working load ; S33, calculating the true residual strain of the part by using the first dimension and the second dimension Specifically, the formula Calculated.

[0031] Further, in step S4, the bilinear constitutive parameter is obtained by evaluating and correcting the initial value of the bilinear constitutive parameter according to the strain threshold through the engineering residual strain and the true residual strain, including: S41, using the engineering residual strain and the true residual strain, according to the formula Calculate the ratio.

[0032] S42. When the ratio is less than or equal to a given strain threshold, the initial value of the bilinear constitutive parameter is judged to be reasonable, and the initial value of the bilinear constitutive parameter is used as the final bilinear constitutive parameter. This parameter can be used to carry out elastic-plastic finite element calculation of the part under other loads, wherein the given strain threshold can be set to 10%.

[0033] S43. When the ratio is greater than a given strain threshold, it is determined that the initial values ​​of the bilinear constitutive parameters are unreasonable, and the initial values ​​of the bilinear constitutive parameters are corrected.

[0034] Furthermore, in step S43, when the ratio is greater than a given strain threshold, it is determined that the initial value of the bilinear constitutive parameter is unreasonable, and the initial value of the bilinear constitutive parameter is corrected, including: S431, using the first engineering stress, the first engineering strain, the elastic modulus and the true residual strain, correcting the first engineering stress of the part under the working load determined by the engineering fracture point T, to obtain the second true stress of the part at the working load point H and the second true strain Specifically, we can use the formula The second true strain is calculated as .

[0035] S44, calculating bilinear constitutive parameters according to the elastic modulus, the true stress, the true strain, the second true stress and the second true strain. Specifically, after obtaining the second true stress and the second true strain of the working load point H, they can be used to replace the first true strain and the first true stress in steps S141 and S142, and recalculate the bilinear constitutive parameters of the material to be calibrated.

[0036] Based on the same inventive concept, a bilinear constitutive parameter engineering establishment system is also provided in the embodiments of the present invention, as described in the following embodiments. Since the principle of solving the problem by the bilinear constitutive parameter engineering establishment system is similar to the bilinear constitutive parameter engineering establishment method shown in the above embodiments, the implementation of the bilinear constitutive parameter engineering establishment system can refer to the implementation of the bilinear constitutive parameter engineering establishment method, and the repeated parts will not be repeated. As used below, the term "unit" or "module" can be a combination of software and / or hardware that implements a predetermined function. Although the device described in the following embodiments is preferably implemented in software, the implementation of hardware, or a combination of software and hardware, is also possible and conceived.

[0037] Figure 3 is a structural block diagram of a bilinear constitutive parameter engineering establishment system disclosed in an embodiment of the present invention, such as Figure 3 As shown, the bilinear constitutive parameter engineering establishment system includes a bilinear constitutive parameter initial value calculation module 301, a simulation module 302, a loading and unloading test module 303 and an evaluation and correction module 304, and the structure is described below.

[0038] The bilinear constitutive parameter initial value calculation module 301 is used to obtain the bilinear constitutive parameter initial values ​​of the material to be calibrated according to the elastic modulus, material section shrinkage, material yield strength and material tensile strength of the material to be calibrated, and the bilinear constitutive parameter initial values ​​include the yield stress initial value and the tangent modulus initial value; The simulation module 302 is used to perform elastic-plastic finite element analysis on the part according to the initial value of the bilinear constitutive parameter, and obtain the first engineering stress and the first engineering strain of the part under the working load, and the engineering residual strain of the part after the working load is unloaded; The loading and unloading test module 303 is used to perform a working load loading and unloading test on the part to obtain the real residual strain of the part; The evaluation and correction module 304 is used to evaluate and correct the initial value of the bilinear constitutive parameter according to the engineering residual strain and the true residual strain according to the strain threshold to obtain the bilinear constitutive parameter.

[0039] The bilinear constitutive parameter engineering establishment method of the present invention first calculates an initial value of a bilinear constitutive parameter through tensile mechanical property data in a material manual, and then performs elastic-plastic loading and unloading of a working load through simulation technology and real tests, and evaluates and corrects the initial value of the bilinear constitutive parameter to obtain the final bilinear constitutive parameter. The method of the present invention obtains a constitutive curve that can meet engineering use through tensile mechanical property data, provides data support for the strength design of engine parts, and has the advantages of few required parameters, easy operation, and strong engineering operability.

[0040] In this embodiment, a computer device is provided, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the computer program, any of the above-mentioned bilinear constitutive parameter engineering establishment methods is implemented.

[0041] Specifically, the computer device may be a computer terminal, a server or a similar computing device.

[0042] In this embodiment, a computer-readable storage medium is provided, wherein the computer-readable storage medium stores a computer program for executing any of the above-mentioned bilinear constitutive parameter engineering establishment methods.

[0043] Specifically, computer-readable storage media include permanent and non-permanent, removable and non-removable media that can be used to store information by any method or technology. Information can be computer-readable instructions, data structures, program modules or other data. Examples of computer-readable storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technology, compact disk read-only memory (CD-ROM), digital versatile disk (DVD) or other optical storage, magnetic cassettes, tape disk storage or other magnetic storage devices or any other non-transmission media that can be used to store information that can be accessed by a computing device. As defined herein, computer-readable storage media does not include temporary computer-readable media (transitory media), such as modulated data signals and carrier waves.

[0044] Obviously, those skilled in the art should understand that the modules or steps of the above-mentioned embodiments of the present invention can be implemented by a general computing device, they can be concentrated on a single computing device, or distributed on a network composed of multiple computing devices, and optionally, they can be implemented by a program code executable by a computing device, so that they can be stored in a storage device and executed by the computing device, and in some cases, the steps shown or described can be executed in a different order from that here, or they can be made into individual integrated circuit modules, or multiple modules or steps therein can be made into a single integrated circuit module for implementation. In this way, the embodiments of the present invention are not limited to any specific combination of hardware and software.

[0045] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. For those skilled in the art, the embodiments of the present invention may have various modifications and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention shall be included in the protection scope of the present invention.

Claims

1. A method for establishing bilinear constitutive parameters engineering, characterized in that: include: A part made of a material to be calibrated is used to obtain initial values ​​of bilinear constitutive parameters of the material to be calibrated according to the elastic modulus, material cross-sectional shrinkage, material yield strength and material tensile strength of the material to be calibrated, wherein the initial values ​​of bilinear constitutive parameters include an initial value of yield stress and an initial value of tangent modulus; Performing elastic-plastic finite element analysis on the part according to the initial values ​​of the bilinear constitutive parameters to obtain the first engineering stress and the first engineering strain of the part under the working load, and the engineering residual strain of the part after the working load is unloaded; Performing working load loading and unloading tests on the part to obtain the true residual strain of the part; The bilinear constitutive parameters are obtained by evaluating and correcting the initial values ​​of the bilinear constitutive parameters according to the strain threshold through the engineering residual strain and the true residual strain.

2. The method for establishing bilinear constitutive parameters engineering according to claim 1, characterized in that: According to the elastic modulus, material section shrinkage, material yield strength and material tensile strength of the material to be calibrated, the initial value of the bilinear constitutive parameter of the material to be calibrated is obtained, including: Given the residual deformation of the part after load unloading, the engineering stress and engineering strain of the deformation point corresponding to the residual deformation are calculated through the material yield strength; Calculating the true stress and true strain of the deformation point through the engineering stress and the engineering strain; Obtaining the first true strain and the first true stress of the engineering fracture point according to the material cross-sectional shrinkage rate and the material tensile strength; Initial values ​​of bilinear constitutive parameters are obtained according to the elastic modulus, the true stress, the true strain, the first true strain and the first true stress.

3. The method for establishing bilinear constitutive parameters engineering according to claim 2, characterized in that: Obtaining initial values ​​of bilinear constitutive parameters according to the elastic modulus, the true stress, the true strain, the first true strain and the first true stress includes: Obtaining an initial value of the tangent modulus of the material to be calibrated according to the true stress, the true strain, the first true strain and the first true stress; The initial value of the yield stress of the material to be calibrated is obtained according to the initial value of the tangent modulus, the elastic modulus, the first true strain and the first true stress.

4. The method for establishing bilinear constitutive parameters engineering according to claim 1, characterized in that: Performing working load loading and unloading tests on the parts to obtain the true residual strain of the parts, including: Measuring a first dimension of the part before loading the working load; Performing working load and working load unloading tests on the part, and measuring a second dimension of the part after the working load is unloaded; The true residual strain of the part is calculated by using the first dimension and the second dimension.

5. The method for establishing bilinear constitutive parameters engineering according to claim 2, characterized in that: By using the engineering residual strain and the real residual strain, the initial value of the bilinear constitutive parameter is evaluated and corrected according to the strain threshold to obtain the bilinear constitutive parameter, including: Using the engineering residual strain and the true residual strain, according to the formula Calculate the ratio, where is the true residual strain, is the engineering residual strain; When the ratio is less than or equal to a given strain threshold, the initial value of the bilinear constitutive parameter is judged to be reasonable, and the initial value of the bilinear constitutive parameter is used as the final bilinear constitutive parameter; When the ratio is greater than a given strain threshold, it is determined that the initial value of the bilinear constitutive parameter is unreasonable, and the initial value of the bilinear constitutive parameter is corrected.

6. The method for establishing bilinear constitutive parameters engineering according to claim 5, characterized in that: When the ratio is greater than a given strain threshold, it is judged that the initial value of the bilinear constitutive parameter is unreasonable, and the initial value of the bilinear constitutive parameter is corrected, including: S431, using the first engineering stress, the first engineering strain, the elastic modulus and the true residual strain to obtain a second true stress and a second true strain at a working load point of the part; S44. Calculate bilinear constitutive parameters based on the elastic modulus, the true stress, the true strain, the second true stress and the second true strain.

7. A bilinear constitutive parameter engineering establishment system, characterized in that: include: A bilinear constitutive parameter initial value calculation module, wherein the bilinear constitutive parameter initial value calculation module is used to obtain the bilinear constitutive parameter initial values ​​of the material to be calibrated according to the elastic modulus, material cross-sectional shrinkage, material yield strength and material tensile strength of the material to be calibrated, wherein the bilinear constitutive parameter initial values ​​include the yield stress initial value and the tangent modulus initial value; A simulation module, wherein the simulation module is used to perform elastic-plastic finite element analysis on the part according to the initial value of the bilinear constitutive parameter, to obtain the first engineering stress and the first engineering strain of the part under the working load, and the engineering residual strain of the part after the working load is unloaded; A loading and unloading test module, wherein the loading and unloading test module is used to perform a working load loading and unloading test on the part to obtain a true residual strain of the part; An evaluation and correction module is used to evaluate and correct the initial value of the bilinear constitutive parameter according to the engineering residual strain and the true residual strain according to the strain threshold to obtain the bilinear constitutive parameter.

8. A computer device comprising a memory, a processor and a computer program stored in the memory and executable on the processor, characterized in that: When the processor executes the computer program, the bilinear constitutive parameter engineering establishment method according to any one of claims 1 to 6 is implemented.

9. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program for executing the bilinear constitutive parameter engineering establishment method according to any one of claims 1 to 6.

Images