Carbon fiber composite material mechanical property adjusting method based on strain rate effect

Abstract

The invention discloses a method for adjusting mechanical properties of a carbon fiber composite material based on a strain rate effect, which comprises the following steps: acquiring static mechanical parameters of carbon fibers, strain rate sensitive parameters of a matrix material and dynamic shear characteristic parameters of a fiber-matrix interface, and constructing a fiber-matrix-interface ternary parameter set; based on the ternary parameter set, calculating the dynamic ultimate strength of the matrix by adopting a Cowper-Symonds model introduced with an interface shear correction term, and generating a corrected matrix dynamic performance parameter set covering a preset strain rate range; according to the corrected dynamic performance parameter set of the matrix, the performance of the composite material is calculated through a three-phase collaborative mixing rule, and a dynamic performance parameter set of the composite material is output; and comparing the dynamic performance parameter set of the composite material with multi-strain-rate experimental data, and constructing a strain rate-interface coupling performance model. According to the embodiment of the invention, the prediction precision of the dynamic mechanical property of the composite material can be improved.

Description

TECHNICAL FIELD

[0001] The application belongs to the technical field of composite materials, and particularly relates to a method for adjusting mechanical properties of carbon fiber composite materials based on strain rate effect. BACKGROUND

[0002] Carbon fiber composite materials are widely used in the fields of aerospace and automobiles due to their high specific strength and high specific modulus, but their mechanical properties are significantly affected by strain rate effect under dynamic load. Existing researches mainly focus on the strain rate sensitivity of fibers or matrix single components, and ignore the dynamic mechanical behavior of fiber-matrix interface and its regulation on overall performance. The traditional mixing method does not consider the coupling relationship between the dynamic shear properties of the interface phase and the strain rate, resulting in insufficient simulation accuracy under high strain rate conditions. In addition, the existing model is only applicable to homogeneous materials and is difficult to characterize the synergistic effect of the three phases of the composite material. SUMMARY

[0003] The application aims to provide a method for adjusting mechanical properties of carbon fiber composite materials based on strain rate effect, so as to solve the problems in the prior art and improve the prediction accuracy of dynamic mechanical properties of composite materials by constructing a ternary parameter set of fiber-matrix-interface and a strain rate-interface coupling model.

[0004] One embodiment of the present application provides a method for adjusting mechanical properties of carbon fiber composite materials based on strain rate effect, which comprises the following steps:

[0005] Obtaining static mechanical parameters of carbon fibers, strain rate sensitive parameters of matrix materials and dynamic shear property parameters of fiber-matrix interface, and constructing a ternary parameter set of fiber-matrix-interface, wherein the dynamic shear property parameters include interface shear strength and thickness coefficient under different strain rates;

[0006] Based on the ternary parameter set, the Cowper-Symonds model introducing an interface shear correction term is used to calculate the dynamic ultimate strength of the matrix, and the logarithmic relationship formula containing the interface thickness coefficient is used to derive the elastic modulus of the matrix, thereby generating a set of corrected dynamic performance parameters of the matrix covering a preset strain rate range;

[0007] According to the set of corrected dynamic performance parameters of the matrix, the performance of the composite material is calculated by a three-phase synergistic mixing rule, the composite material is regarded as a three-phase system of fibers, matrix and interface phase, the main and auxiliary direction elastic moduli, shear modulus and dynamic tensile, compression and shear strength are calculated by volume fraction weighting, and a set of dynamic performance parameters of the composite material is outputted;

[0008] The composite material dynamic performance parameter set is compared with the multi-strain rate experimental data, an interface parameter sensitivity analysis is performed to optimize the interface coefficient weight in the ternary parameter set, a strain rate-interface coupling performance model is constructed, and is used as the simulation input under the high strain rate loading condition.

[0009] Optionally, the static mechanical parameters of the carbon fiber, the strain rate sensitive parameters of the matrix material, and the dynamic shear characteristic parameters of the fiber-matrix interface are obtained, a ternary parameter set of the fiber-matrix-interface is constructed, wherein the dynamic shear characteristic parameters include the interface shear strength and thickness coefficient under different strain rates, and the ternary parameter set includes:

[0010] A constant temperature static tensile experimental device is used to test the carbon fiber single wire under a strain rate of 0.001 s -1 The elastic modulus, the ultimate strength and the elongation at break are recorded, the average value is obtained through three parallel experiments, and a carbon fiber static mechanical parameter table is generated.

[0011] A split Hopkinson pressure bar device is used to perform a dynamic compression experiment on the matrix material in a strain rate range of 0.001 s -1 to 1000 s -1 The yield strength and the elastic modulus change data under different strain rates are collected, the strain rate related coefficients C, the strain rate sensitive index p and the elastic modulus sensitive coefficient k in the Cowper-Symonds model parameters are fitted, a matrix strain rate sensitive parameter table is generated, and the strain rate related coefficients C, the strain rate sensitive index p and the elastic modulus sensitive coefficient k in the Cowper-Symonds model parameters are fitted.

[0012] A micro-column indentation method is used to perform a dynamic test on the fiber-matrix interface, different strain rates are simulated by controlling the loading rate of the indenter, the relationship curve between the indentation depth and the shear load is recorded, the interface shear strength is calculated, the interface micro morphology is observed by using a scanning electron microscope, the interface thickness coefficient is extracted, and an interface dynamic shear characteristic parameter table is generated.

[0013] The carbon fiber static mechanical parameter table, the matrix strain rate sensitive parameter table and the interface dynamic shear characteristic parameter table are associated according to the fiber matrix volume fraction, a mapping relationship between the parameters is established, and a ternary parameter set of the fiber-matrix-interface is constructed.

[0014] Optionally, based on the ternary parameter set, the Cowper-Symonds model is used to calculate the dynamic ultimate strength of the matrix by introducing an interface shear correction term, the matrix elastic modulus is derived by combining a logarithmic relationship formula containing the interface thickness coefficient, a corrected matrix dynamic performance parameter set covering a preset strain rate range is generated, and the corrected matrix dynamic performance parameter set includes:

[0015] The interface shear strength parameter is extracted from the ternary parameter set, the ratio of the interface shear strength parameter to the static ultimate strength of the matrix is calculated as an interface shear correction coefficient, the interface shear correction coefficient is dynamically adjusted according to a linear relationship as the strain rate increases, and an interface shear correction coefficient sequence is generated.

[0016] The interface shear correction coefficient sequence is substituted into the Cowper-Symonds model, the dynamic ultimate strength calculation result is adjusted by adding a correction term, the dynamic ultimate strength of the matrix at different strain rates is obtained, and a corrected matrix dynamic ultimate strength sequence is generated;

[0017] The interface thickness coefficient is extracted from the three-parameter set and introduced into the logarithmic relationship formula of the matrix elastic modulus as a multiplier, the change amplitude of the elastic modulus with the strain rate is adjusted through the thickness coefficient, the elastic modulus at different strain rates is calculated, and a corrected matrix elastic modulus sequence is generated;

[0018] The corrected matrix dynamic ultimate strength sequence and the elastic modulus sequence are integrated, and data is completed at a strain rate interval of 0.001 s -1 to 1000 s -1 , to generate a corrected matrix dynamic performance parameter set covering the preset strain rate range.

[0019] Optionally, the composite material performance is calculated according to the corrected matrix dynamic performance parameter set through a three-phase synergistic mixing rule, the composite material is regarded as a three-phase system of fiber, matrix and interface phase, the main and auxiliary direction elastic modulus, shear modulus and dynamic tensile, compression and shear strength are calculated according to the volume fraction weighting, and a composite material dynamic performance parameter set is output, including:

[0020] The interface phase volume of a single fiber is calculated according to the fiber diameter and the interface thickness coefficient, the three-phase volume fraction proportions of the fiber, the matrix and the interface phase are derived in combination with the fiber volume fraction and the matrix volume fraction, and a three-phase volume fraction distribution table is generated;

[0021] Based on the three-phase volume fraction distribution table, the main direction and the auxiliary direction elastic modulus are calculated by using a three-phase synergistic mixing rule, wherein the fiber and the interface phase use static elastic modulus, and the matrix uses corrected dynamic elastic modulus, the results are obtained by weighted summation, and a composite material elastic modulus parameter to be integrated is generated;

[0022] The shear modulus is calculated according to the same mixing rule, wherein the shear modulus of the fiber and the interface phase takes the static value, the shear modulus of the matrix takes the corrected dynamic value, and the composite material shear modulus parameter is generated by weighted calculation according to the three-phase volume fraction;

[0023] The dynamic tensile, compression and shear strength are calculated respectively, wherein the fiber and the interface phase take the static strength parameter, the matrix takes the corrected dynamic strength parameter, the three-phase volume fraction is weighted and summed, and the composite material elastic modulus parameter and the composite material shear modulus parameter to be integrated are integrated, to obtain the composite material dynamic performance parameter set.

[0024] Optionally, the set of composite material dynamic performance parameters is compared with the multi-strain rate experimental data, the interface coefficient weight in the ternary parameter set is optimized through interface parameter sensitivity analysis, a strain rate-interface coupling performance model is constructed as the simulation input under high strain rate loading conditions, including:

[0025] The set of composite material dynamic performance parameters is compared with the experimental test data under different strain rates, the strength and modulus error values corresponding to each strain rate are calculated, and a performance error sequence is generated;

[0026] The performance error sequence is subjected to sensitivity analysis, the interface shear strength and thickness coefficient are adjusted separately through the control variable method, the change amplitude of the error value is observed, the influence weight of the two parameters on the error is determined, and an interface parameter sensitivity weight table is generated;

[0027] According to the sensitivity weight table, the interface shear strength and thickness coefficient in the ternary parameter set are adjusted according to the error minimization principle, and the iteration optimization is performed until the error value is lower than the preset threshold, and the optimized ternary parameter set is generated;

[0028] Based on the optimized ternary parameter set, the performance parameters of the composite material in the full strain rate range are recalculated, and a strain rate-interface coupling performance model is constructed as the simulation input under high strain rate loading conditions.

[0029] Another embodiment of the application provides a carbon fiber composite material mechanical property adjustment system based on strain rate effect, the system comprising:

[0030] An acquisition module is configured to acquire static mechanical parameters of carbon fibers, strain rate sensitive parameters of matrix materials, and dynamic shear characteristic parameters of fiber-matrix interfaces, and construct a ternary parameter set of fibers-matrix-interfaces, wherein the dynamic shear characteristic parameters include interface shear strength and thickness coefficients under different strain rates;

[0031] A derivation module is configured to calculate the dynamic ultimate strength of the matrix by using a Cowper-Symonds model with an introduced interface shear correction term based on the ternary parameter set, derive the elastic modulus of the matrix by combining a logarithmic relationship formula containing the interface thickness coefficient, and generate a corrected matrix dynamic performance parameter set covering a preset strain rate range;

[0032] A calculation module is configured to calculate the performance of the composite material by using a three-phase synergistic mixing rule according to the corrected matrix dynamic performance parameter set, treat the composite material as a three-phase system of fibers, matrixes, and interfaces, calculate the main and auxiliary direction elastic moduli, shear moduli, and dynamic tensile, compression, and shear strengths by volume fraction weighting, and output a set of composite material dynamic performance parameters;

[0033] A construction module is configured to compare the composite material dynamic performance parameter set with multi-strain rate experimental data, optimize the interface coefficient weight in the ternary parameter set through interface parameter sensitivity analysis, construct a strain rate-interface coupling performance model as a simulation input under a high strain rate loading condition.

[0034] A further embodiment of the present application provides a storage medium having a computer program stored therein, wherein the computer program is configured to execute the method described in any of the above when executed.

[0035] A further embodiment of the present application provides an electronic device comprising a memory having a computer program stored therein and a processor configured to execute the computer program to execute the method described in any of the above.

[0036] Compared with the prior art, the method for adjusting the mechanical properties of the carbon fiber composite material based on the strain rate effect provided by the present application acquires the static mechanical parameters of the carbon fiber, the strain rate sensitive parameters of the matrix material and the dynamic shear characteristic parameters of the fiber-matrix interface, and constructs a ternary parameter set of the fiber-matrix-interface; based on the ternary parameter set, the Cowper-Symonds model introducing an interface shear correction term is used to calculate the dynamic ultimate strength of the matrix, and a corrected matrix dynamic performance parameter set covering a preset strain rate range is generated; according to the corrected matrix dynamic performance parameter set, the performance of the composite material is calculated through a three-phase synergistic mixing rule, and a composite material dynamic performance parameter set is output; the composite material dynamic performance parameter set is compared with multi-strain rate experimental data, and a strain rate-interface coupling performance model is constructed, so that the prediction accuracy of the dynamic mechanical properties of the composite material can be improved by constructing the ternary parameter set of the fiber-matrix-interface and the strain rate-interface coupling model. BRIEF DESCRIPTION OF DRAWINGS

[0037] Figure 1 A hardware structure block diagram of a computer terminal of the method for adjusting the mechanical properties of the carbon fiber composite material based on the strain rate effect provided by the embodiment of the present application is provided.

[0038] Figure 2 A flowchart of the method for adjusting the mechanical properties of the carbon fiber composite material based on the strain rate effect provided by the embodiment of the present application is provided.

[0039] Figure 3 A structure diagram of the system for adjusting the mechanical properties of the carbon fiber composite material based on the strain rate effect provided by the embodiment of the present application is provided. DETAILED DESCRIPTION

[0040] The embodiments described below with reference to the drawings are exemplary and are only used to explain the present application, and cannot be explained as a limitation of the present application.

[0041] The embodiment of the present application first provides a carbon fiber composite material mechanical property adjustment method based on strain rate effect, which can be applied to electronic devices, such as computer terminals, specifically, common computers and the like.

[0042] The following will be described in detail taking the running on the computer terminal as an example. Figure 1 The hardware structure block diagram of the computer terminal provided by the carbon fiber composite material mechanical property adjustment method based on strain rate effect is shown in the embodiment of the present application. As shown in the figure, Figure 1 The computer device includes a processor, a memory and a network interface connected through a system bus, wherein the memory can include a non-volatile storage medium and an internal memory.

[0043] The non-volatile storage medium can store an operating system and a computer program. The computer program includes program instructions, which, when executed, can make the processor execute any kind of carbon fiber composite material mechanical property adjustment method based on strain rate effect.

[0044] The processor is used to provide computing and control capabilities to support the operation of the entire computer device.

[0045] The internal memory provides an environment for the running of the computer program in the non-volatile storage medium, which, when executed by the processor, can make the processor execute any kind of carbon fiber composite material mechanical property adjustment method based on strain rate effect.

[0046] The network interface is used for network communication, such as sending assigned tasks and the like. Those skilled in the art can understand that Figure 1 The structure shown in the figure is only a block diagram of part of the structure related to the scheme of the present application, and does not constitute a limitation on the computer device to which the scheme of the present application is applied. The specific computer device can include more or less components than those shown in the figure, or combine certain components, or have a different component arrangement.

[0047] It should be understood that the processor can be a central processing unit (CPU), and the processor can also be other general-purpose processors, digital signal processors (DSP), application specific integrated circuits (ASIC), field programmable gate arrays (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. Among them, the general-purpose processor can be a microprocessor or the processor can also be any conventional processor, etc.

[0048] Referring to Figure 2 Embodiments of the present application provide a method for adjusting the mechanical properties of carbon fiber composites based on the strain rate effect, which can include the following steps:

[0049] S201, obtain the static mechanical parameters of carbon fibers, the strain rate sensitive parameters of the matrix material, and the dynamic shear characteristic parameters of the fiber-matrix interface, and construct a ternary parameter set of the fiber-matrix interface, wherein the dynamic shear characteristic parameters include the interfacial shear strength and thickness coefficient under different strain rates;

[0050] Specifically, a constant temperature static tensile test device can be used to test the carbon fiber monofilament under a strain rate of 0.001 s -1 , record the elastic modulus, ultimate strength and elongation at break, take the average value through three parallel experiments, and generate a carbon fiber static mechanical parameter table;

[0051] Experimental device configuration and environmental control

[0052] The core of the constant temperature static tensile test device is a precision electronic universal testing machine (such as Instron 5967 type), equipped with an environmental box with a temperature control accuracy of ±0.5°C. Before testing, the carbon fiber monofilament (diameter 7 microns, length 50 millimeters) is fixed at both ends on a paper fixture with a special adhesive, and the fixture is clamped into a pneumatic clamp to avoid damage to the fiber caused by clamping force. The temperature of the environmental box is set to 23±0.5°C and the humidity is 50±5% RH to eliminate the influence of temperature and humidity on the resin-based fiber. The strain rate is strictly controlled at 0.001 s -1 (i.e. 0.001 times the initial gage length per second of tensile stretching), which is realized in displacement control mode through the servo system of the testing machine with a displacement resolution of 0.1 microns. The force sensor has a range of 5kN and an accuracy of 0.5%, and real-time load data is collected; the extensometer has a gage length of 25 millimeters and an accuracy of 0.5 microns, and the deformation amount within the gage length is recorded synchronously.

[0053] Test procedure and data acquisition

[0054] After starting the tensile program, the testing machine applies an axial tensile force at a constant rate. The load-displacement curve is displayed in real time on the control software (such as Bluehill Universal), and when the curve shows a clear inflection point (yield point), the ultimate strength (unit: megapascal MPa) is recorded, and the elastic modulus (unit: gigapascal GPa) is calculated by the slope of the linear segment (stress increment / strain increment). A high-speed camera (1000 frames / second) is triggered at the moment of fracture to capture the fracture position, and the elongation at break (fracture displacement / initial gage length x 100%) is calculated in combination with the extensometer data. The single test takes about 30 minutes, each group of carbon fiber samples is repeated three times independently, and a new sample is replaced each time to avoid fatigue effects.

[0055] Data processing and parameter table generation

[0056] The three experimental data are imported into a statistical analysis software (such as Minitab), and abnormal values (such as data deviation caused by clamping slip) are removed. The elastic modulus takes the arithmetic mean of the linear regression slopes of the three experiments; the ultimate strength takes the average of the three peak values; and the elongation at break is calculated according to the average of the three fracture displacements. Finally, a structured parameter table is generated, which contains three columns of key data:

[0057] Elastic modulus (E_f): for example, 230 GPa, representing the fiber's resistance to deformation;

[0058] Ultimate strength (σ_f_max): for example, 3.5 GPa, representing the fiber's maximum load-carrying capacity;

[0059] Elongation at break (ε_f): for example, 1.5%, representing the fiber's ductility.

[0060] Using a split Hopkinson pressure bar device, dynamic compression experiments are conducted on the matrix material at a strain rate range of 0.001 s -1 to 1000 s -1 The yield strength and elastic modulus change data at different strain rates are collected, and the strain rate related coefficients C, strain rate sensitive index p and elastic modulus sensitive coefficient k in the Cowper-Symonds model parameters are fitted to generate a matrix strain rate sensitive parameter table;

[0061] The split Hopkinson pressure bar (SHPB) device is composed of a launch system, an impact rod, an incident rod, a transmission rod and an absorption rod (all made of high-strength steel with a diameter of 20 millimeters). The matrix material (such as epoxy resin EPOLAM 5015) is made into a cylindrical sample (diameter 5 millimeters, height 5 millimeters). Before the experiment, the rod system is adjusted to have a coaxiality error of less than 0.02 millimeters. The low strain rate section (0.001 s -1 -1 s -1 ) uses hydraulic servo control to compress the sample at a constant rate; the high strain rate section (10 s -1 -1000 s -1 ) uses a gas gun to launch the impact rod, and the impact velocity is calibrated by a laser speedometer (for example, 20 meters per second corresponds to about 500 s -1 ). The sample is coated with lubricating grease at both ends to reduce friction, and strain gauges (gauge length 2 millimeters) are attached to the incident rod and the transmission rod, which are connected to a dynamic strain meter (such as Vishay 7000 series) to collect stress wave signals, with a sampling frequency of 10 megahertz.

[0062] The dynamic compression experiment is divided into groups according to the strain rate gradient: 0.001 s -1 (quasi-static), 0.1 s -11s -1 10s -1 100s -1 500s -1 1000s -1 There are 7 groups in total, each group repeated 5 times. The original signal is processed using the three-wave method: incident wave εi(t), reflected wave εi(t), and reflected wave εi(t). r (t), transmitted wave ε t Substitute (t) into the calculated specimen stress σ(t) and strain ε(t). Extract the yield strength (stress corresponding to 0.2% plastic strain, in MPa) and elastic modulus (slope of the initial linear segment, in GPa) from the stress-strain curve. For example, at 100 s... -1 Under these conditions, the yield strength of the epoxy resin increased from 50 MPa to 85 MPa, and the modulus increased from 3.5 GPa to 4.2 GPa. Response curves were plotted using data at different strain rates, with the horizontal axis representing the logarithmic strain rate. The vertical axis represents the strength ratio σ_d / σ_s (the ratio of dynamic yield strength to static yield strength) and the modulus ratio E_d / E_s.

[0063] Based on data fitting, the Cowper-Symonds model describes the relationship between dynamic yield strength and strain rate, expressed as follows: The strain rate correlation coefficient C (unit: s) -1 The strain rate sensitivity threshold reflects the material's strain rate sensitivity, and the strain rate sensitivity index p (dimensionless) controls the growth rate. The least squares method is used to fit the experimental data: for example, C = 80 s⁻¹. -1 This indicates that when the strain rate reaches 80s -1 The strength increases by a factor of two. The elastic modulus fits an independent model: The elastic modulus sensitivity coefficient k (dimensionless) characterizes the magnitude of modulus growth. The reference strain rate is used. The final matrix strain rate sensitive parameter table is generated, including C values ​​(e.g., 80.2 s). -1 p-value (e.g., 3.8), k-value (e.g., 0.15), and correlation coefficient R0. 2 (Must be greater than 0.98).

[0064] Dynamic testing was conducted at the fiber-matrix interface using the micropillar indentation method. Different strain rates were simulated by controlling the indenter loading rate, and the relationship curve between indentation depth and shear load was recorded. The interfacial shear strength was calculated, and the microstructure of the interface was observed using scanning electron microscopy. The interfacial thickness coefficient was extracted, and a table of dynamic shear characteristic parameters of the interface was generated.

[0065] Micro-pillar indentation method uses a nanoindenter (e.g. Hysitron TI 980) equipped with a wedge-shaped indenter (90° tip angle, 100 nm tip radius). The specimen is a carbon fiber / matrix micro-pillar array: first, a 5x5 micro-pillar array (10 μιη height) is cut by focused ion beam (FIB) on the cross-section of the composite, with each pillar containing a single fiber (centered) and symmetric matrix on both sides. The indenter is positioned with 0.1 μιη precision, and the indentation test is performed right above the fiber-matrix interface. Dynamic tests are achieved by controlling the indenter displacement rate: four levels of 0.01 μιη / s, 0.1 μιη / s, 1 μιη / s, 10 μιη / s are set, corresponding to local strain rates of 1 s -1 , 10 s -1 , 100 s -1 , 1000 s -1 . The load-displacement curve is recorded simultaneously (sampling frequency 100 Hz), with emphasis on monitoring the shear load component (induced by the asymmetric structure of the wedge-shaped indenter).

[0066] The interfacial shear strength is calculated based on the maximum shear load criterion: when the indentation depth reaches the center of the interface (about 5 μιη), the maximum load value F_max (unit μN) corresponding to the shear load drop point is divided by the shear area A_s (micro-pillar width x indentation depth, e.g. 5 μιη x 5 μιη = 25 μιη 2 ), and the shear strength τ = F_max / A_s (unit MPa) is obtained. 10 micro-pillars are tested at each strain rate, and the average value is taken after removing outliers. For example, at 1000 s -1 , the shear strength increases from 40 MPa in static to 65 MPa. Scanning electron microscopy (SEM) observation is performed after testing: using a field emission electron microscope (e.g. ZEISS Gemini 500) to observe the interface morphology in the indentation area at 5 kV voltage. The interfacial transition zone (ITZ) thickness is measured: 10 line scans are taken along the fiber radial direction, and the energy dispersive spectrometer (EDS) is used to analyze the element gradient (e.g. carbon element content change), and the element mutation interval is defined as the interface thickness δ (unit nm).

[0067] The dynamic test data are combined with the morphology analysis to establish the correlation. With the increase of strain rate, SEM observation finds that the interface damage mode changes from ductile debonding (matrix residue) to brittle fracture (clean interface), and the interface thickness δ decreases with the increase of strain rate (e.g. from 200 nm in static to 120 nm at 1000 s -1 ). Finally, the interface dynamic shear characteristic parameter table is generated, containing four columns of data: strain rate (s -1 ), average interface shear strength (MPa), standard deviation (MPa), average interface thickness coefficient δ (nm). The thickness coefficient is defined as δ / δ0 (δ0 is the static interface thickness), which is used for subsequent model correction.

[0068] The carbon fiber static mechanical parameter table, matrix strain rate sensitive parameter table, and interface dynamic shear characteristic parameter table are correlated according to the fiber-matrix volume fraction to establish the mapping relationship between the parameters and construct a ternary parameter set of fiber-matrix-interface.

[0069] Parameter correlation needs to be based on the microstructure of the composite material. First, determine the fiber volume fraction V_f (e.g., 60%) and matrix volume fraction V_m (e.g., 38%) of the target composite material. The interfacial phase volume fraction V_i is calculated using the formula V_i = V_f × [(d_f + 2δ)^2 / d_f^2 - 1] (d_f is the fiber diameter, δ is the interfacial thickness). For example, when d_f = 7 μm and δ = 0.2 μm, V_i ≈ 2% (of the total volume). Establish the logical correlation between three parameter tables: the carbon fiber static parameter table provides the inherent fiber properties (E_f, σ_f); the matrix strain rate sensitive parameter table provides the dynamic correction function (C, p, k); and the interfacial dynamic parameter table provides strain rate-related parameters. and

[0070] The mapping relationship is constructed by defining interface variables: in terms of strain rate. For public index (0.001s) -1 up to 1000s -1 For each strain rate node:

[0071] Calculate the dynamic yield strength using the Cowper-Symonds model from the matrix parameter table.

[0072] This can be obtained from the interpolation of the interface parameter table. Below and

[0073] The carbon fiber parameters remain static (the strain rate sensitivity of the carbon fiber is negligible);

[0074] Generate a two-dimensional relational matrix, where rows correspond to strain rates and columns include: V_f, V_m, Vi, E_f, σ_f, There are 12 parameters.

[0075] The final ternary parameter set is stored in a database, supporting strain rate queries.

[0076] S202, Based on the ternary parameter set, the dynamic ultimate strength of the matrix is ​​calculated using the Cowper-Symonds model with an interface shear correction term, and the elastic modulus of the matrix is ​​derived by combining the logarithmic relationship formula including the interface thickness coefficient, thereby generating a modified matrix dynamic performance parameter set covering a preset strain rate range.

[0077] Specifically, the interface shear strength parameter can be extracted from the ternary parameter set, and the ratio of the interface shear strength parameter to the static ultimate strength of the matrix is calculated as the interface shear correction coefficient, which is dynamically adjusted in a linear relationship with the increase of the strain rate, and a sequence of interface shear correction coefficients is generated.

[0078] Firstly, the fiber-matrix-interface ternary parameter set that has been constructed is called, which is stored in the form of a structured database, containing three sub-libraries of carbon fiber static mechanical parameters, matrix strain rate sensitive parameters, and interface dynamic shear characteristic parameters. The system locates the interface dynamic shear characteristic parameter table through data indexing, and extracts the key fields: the interface shear strength (ISS) values corresponding to different strain rate levels (such as 10s -1 , 100s -1 , 500s -1 , etc.). For example, the ISS value recorded at 100s -1 strain rate is 45 megapascals (MPa). At the same time, the matrix static ultimate strength (MSUS) is extracted from the matrix strain rate sensitive parameter table, which is the benchmark value of the matrix material failure strength measured at the reference strain rate 0.001s -1 (such as 80 MPa). The system performs calculation for each strain rate point: divide the ISS value at the current strain rate by the MSUS value to obtain the initial interface shear modification coefficient (IISMCo) corresponding to the strain rate. For example, at 100s -1 , IISMCo = 45 MPa / 80 MPa = 0.5625.

[0079] Due to the limited experimental data points (usually only covering several discrete strain rates), the system needs to construct a continuous function relationship. By analyzing the trend of IISMCo at different strain rates (for example, 0.40 at 10s -1 , 0.56 at 100s -1 , and 0.72 at 500s -1 ), it is found that it approximately presents a linear growth relationship. The system uses the least squares method to fit a straight line: let the strain rate be the independent variable (s -1 ), and IISMCo be the dependent variable, the fitting equation is where the slope a (e.g., 0.00065 s) reflects the growth rate of the correction factor with strain rate, and the intercept b (e.g., 0.395) represents the baseline correction value when approaching static state. This linear relationship is defined as the Dynamic Adjustment Rule (DAR). Based on the DAR, the system calculates the corresponding correction factor at fixed intervals (e.g., every 50 s -1 up to 1000 s -1 ) within a pre-defined continuous strain rate range (0.001 s -1 ). For example, at 300 s -1 , IISMCo = 0.00065 x 300 + 0.395 = 0.590.

[0080] The final generated Interface Shear Modification Coefficient Sequence (ISMCoSeq) is a two-dimensional array. The first column is the strain rate value in ascending order (e.g., [0.001, 50, 100, 150,..., 1000] s -1 ), and the second column is the corresponding correction factor calculated by the DAR (e.g., [0.395, 0.428, 0.460, 0.493,..., 0.995]). The sequence is stored in standard JSON or CSV format, with each data point containing values accurate to four decimal places. The system verifies the monotonicity of the sequence (ensuring that the coefficient strictly increases with strain rate) and the boundary rationality (e.g., the coefficient does not exceed 1.2 at 1000 s -1 ) through the data verification module. Abnormal values will trigger the experimental data review process.

[0081] Substitute the Interface Shear Modification Coefficient Sequence into the Cowper-Symonds model, adjust the dynamic ultimate strength calculation results by adding a correction term, obtain the matrix dynamic ultimate strength at different strain rates, and generate the corrected matrix dynamic ultimate strength sequence;

[0082] The Cowper-Symonds model (Cowper-Symonds Model, CSM) is a classic model for describing material dynamic strengthening, and its standard form is: where is the strain rate, and C (strain rate related coefficient) and p (strain rate sensitivity index) are parameters obtained from matrix dynamic compression experiments (e.g., C = 1200 s -1 , p = 4.2). The system first calls the CSM parameters in the ternary parameter set to calculate the uncorrected matrix dynamic ultimate strength (UMDUS). For example, at Time: UMDUS = 80 MPa x [1 + (200 / 1200)^(1 / 4.2)] = 102.3 MPa.

[0083] To embody the enhancement of interface effect on the strength of matrix, a modification term (MT) is added to CSM. The specific operation is: querying the interface shear modification coefficient k_iss corresponding to the current strain rate from ISMCoSeq (for example, k_iss = 0.525 at 200 s -1 ). Define the modification factor as (1 + k_iss), multiply it with the original CSM formula to get the modified model: Modified MDUS (MMDUS) = UMDUS x (1 + k_iss). Still take 200 s -1 as an example: MMDUS = 102.3 MPa x (1 + 0.525) = 156.0 MPa. The physical meaning of this modification term is: the higher the interface shear strength (the larger k_iss), the stronger the constraint on the matrix, making the matrix show higher equivalent bearing capacity. The system performs this calculation for each point in the preset strain rate sequence.

[0084] The generated Modified Matrix Dynamic Ultimate Strength Sequence (MMDUSSeq) has the same structure as ISMCoSeq, containing strain rate-strength value pairs. For example, at [50, 200, 1000] s -1 , the strength values are [112.8, 156.0, 254.7] MPa, respectively. The system draws the strength-strain rate curve through the visualization module and compares it with the unmodified CSM curve (the modified curve is shifted upward as a whole and the slope is increased), verifying the rationality of the interface modification. This sequence is written as the "ultimate strength" sub-table of the modified matrix dynamic performance parameter set.

[0085] The interface thickness coefficient is extracted from the ternary parameter set and introduced as a multiplier into the logarithmic relationship formula of the matrix elastic modulus. The thickness coefficient adjusts the amplitude of the change of the elastic modulus with the strain rate, and the elastic modulus at different strain rates is calculated to generate the modified matrix elastic modulus sequence.

[0086] Interface Thickness Coefficient (ITC) is extracted from the Interface Dynamic Shear Property Table. This coefficient is the ratio of the average interface thickness (e.g. 0.8 microns) to the fiber diameter (e.g. 7 microns) measured by scanning electron microscopy (ITC = 0.114), reflecting the relative size of the interface. Meanwhile, the Logarithmic Relation Formula (LRF) describing the strain rate effect on the matrix modulus is obtained from the Matrix Strain Rate Sensitivity Table: where E_sta is the static modulus (e.g. 3.2 Giga Pascals, GPa), k is the elastic modulus sensitivity coefficient (e.g. 0.18), is the reference strain rate (0.001 s -1 ).

[0087] ITC is introduced as a Modulation Multiplier (MM) into LRF, forming the modified formula: Modified Dynamic Elastic Modulus

[0088] where (1 + ITC) constitutes a modulus amplification modulation factor. For example, when ITC = 0.114, the coefficient becomes k_new = 0.18 x (1 + 0.114) = 0.2005. The physical meaning is that a thicker interface layer (larger ITC) can more effectively transfer the strain rate load and amplify the rate-sensitive effect of the matrix modulus. Take as an example for calculation:

[0089] Unmodified Modulus = 3.2 x [1 + 0.18 x log10(500 / 0.001)] = 3.2 x [1 + 0.18 x 5.699] ≈ 5.28 GPa;

[0090] Modified Modulus = 3.2 x [1 + 0.2005 x 5.699] ≈ 5.67 GPa (amplification of about 7.4%).

[0091] Iterate through the preset strain rate sequence, calculate the MDEM value at each point, and generate the Modified Matrix Elastic Modulus Sequence (MMEMSeq). For example, when the strain rate sequence is [1, 100, 1000] s -1 , the modulus is [3.42, 4.95, 6.01] GPa, respectively. The system performs two verifications: 1) Compare the baseline curve when ITC = 0 to ensure that the modified curve maintains a reasonable logarithmic growth trend; 2) Check that the modulus at extremely high strain rates (e.g. 1000 s -1 ) does not exceed the theoretical limit of the material (e.g. 10 GPa). The sequence data is stored in the "Elastic Modulus" sub-table of the Dynamic Performance Parameter Set.

[0092] The modified matrix dynamic ultimate strength sequence and the elastic modulus sequence are integrated, and data completion is performed at intervals of 0.001 s -1 to 1000 s -1 to generate a modified matrix dynamic performance parameter set covering the preset strain rate range.

[0093] The system has obtained two core sequences: MMDUSSeq (ultimate strength sequence) and MMEMSeq (elastic modulus sequence). Since the strain rate sampling points of the two sequences may not be completely consistent (for example, the strength sequence contains 50 s -1 points and the modulus sequence does not contain), data alignment needs to be performed. A linear interpolation algorithm is used: for the missing strain rate points of the modulus sequence between its adjacent known points and the completion value is calculated according to the formula . For example, at , the modulus values at 100 s -1 and 200 s -1 are used to obtain the result.

[0094] To achieve full range coverage, the system generates a uniform strain rate array (USRA) at fixed intervals (such as 10 s -1 ) between the minimum (0.001 s -1 ) and maximum (1000 s -1 ) strain rates. For each strain rate point in each USRA, perform: 1) if the point already exists in the original sequence, directly call the strength / modulus value; 2) if it does not exist, call the interpolation algorithm to complete the data. For example, at the 750 s -1 point (not in the original sequence), the strength values at 500 s -1 and 1000 s -1 are used to obtain MMDUS≈210.4 MPa. Finally, a complete sequence containing 101 data points (0.001 s -1 to 1000 s -1 every 10 s -1 ) is obtained.

[0095] The processed strength sequence and the modulus sequence are combined to construct a structured database - the modified matrix dynamic performance parameter set (MMDPD). The data set contains three columns:

[0096] Column 1: Strain Rate (SR), unit s-1 , precision 0.001;

[0097] Column 2: Modified Ultimate Strength (MUS), unit MPa, precision 0.1;

[0098] Column 3: Modified Elastic Modulus (MEM), unit GPa, precision 0.01.

[0099] For example, at 1000 s -1 Row: SR = 1000, MUS = 254.7, MEM = 6.01. The dataset is stored in HDF5 binary format, accompanied by metadata description (e.g. generation algorithm version, interpolation method). The system automatically generates a verification report, including data range statistics (e.g. modulus minimum 3.20 GPa @ 0.001 s -1 , maximum 6.01 GPa @ 1000 s -1 ) and abnormal point detection log, ensuring the reliability of subsequent hybrid rule calculations.

[0100] S203, according to the modified matrix dynamic performance parameter set, the performance of the composite material is calculated by the three-phase synergistic hybrid rule, the composite material is regarded as a three-phase system of fiber, matrix and interface phase, the main and auxiliary direction elastic modulus, shear modulus and dynamic tensile, compression and shear strength are calculated according to the volume fraction, and the dynamic performance parameter set of the composite material is output;

[0101] Specifically, the interface phase volume of a single fiber can be calculated according to the fiber diameter and interface thickness coefficient, and the three-phase volume fraction proportion of fiber, matrix and interface phase can be derived by combining the fiber volume fraction and the matrix volume fraction, so as to generate a three-phase volume fraction distribution table;

[0102] Microscopic calculation of interface phase volume

[0103] The average diameter (e.g. 7 microns, unit μm) and interface thickness coefficient (e.g. 0.15, dimensionless ratio) of carbon fiber are extracted from the three-parameter set. The interface thickness coefficient is calculated by dividing the actual thickness (e.g. 1.05 microns) obtained by scanning electron microscope (SEM) observation of the fiber-matrix interface region by the fiber radius (3.5 microns). The interface phase of a single fiber is modeled as a uniform cylindrical shell wrapping the fiber, and its volume calculation formula is: cylindrical shell volume = π × (fiber radius + interface thickness) 2 × unit length - π × (fiber radius) 2 × unit length.

[0104] For example, if the fiber radius is 3.5 microns and the interfacial thickness is 1.05 microns, the interfacial phase volume of a single fiber is π x (4.55 2 -3.5 2 ) x 10 -12 cubic meters. This calculation needs to go through all fiber size samples (e.g. diameter distribution 6.8-7.2 microns) and take the average as the interfacial phase volume of a single fiber (e.g. 2.3 x 10 -11 cubic meters per meter).

[0105] Integrated calculation of three-phase volume fractions

[0106] Based on the microstructure of the composite material, input the macroscopic fiber volume fraction (e.g. 60%, i.e. V_f = 0.6) and the matrix volume fraction (e.g. 38%, i.e. V_m = 0.38). Assuming that the total number of fibers in 1 cubic meter of composite material is N, then:

[0107] Total fiber volume = N x single fiber volume;

[0108] Total interfacial phase volume = N x single fiber interfacial phase volume;

[0109] Matrix volume = 1 - total fiber volume - total interfacial phase volume.

[0110] Solve the value of N by simultaneous equations, and finally get the three-phase volume fractions:

[0111] Interfacial phase volume fraction V_i = total interfacial phase volume / 1 cubic meter (e.g. V_i = 0.02, i.e. 2%);

[0112] Corrected fiber volume fraction V_f' = total fiber volume / 1 cubic meter (e.g. V_f' = 0.588);

[0113] Corrected matrix volume fraction V_m' = matrix volume / 1 cubic meter (e.g. V_m' = 0.392).

[0114] Generate a three-phase volume fraction distribution table containing V_f', V_m', V_i three columns of data and corresponding strain rate labels (e.g. 0.001 s -1 , 100 s -1 , etc.).

[0115] Dynamic correlation processing of volume fractions

[0116] Because the interfacial thickness coefficient changes with strain rate (the interface may be compressed and thinned at high strain rate), the volume fractions need to be calculated independently at different strain rates. For example, at 1000 s -1At a strain rate of 100 s-1, the interphase thickness coefficient can be reduced to 0.12, resulting in a V_i of 1.7%. The system generates an independent volume fraction distribution table for each strain rate node through interpolation algorithms (such as cubic spline interpolation), ensuring that the three-phase proportion is accurately associated with the strain rate in subsequent calculations.

[0117] Based on the three-phase volume fraction distribution table, the main direction and secondary direction elastic modulus are calculated using the three-phase cooperative mixing rule, where the fiber and interphase phase use the static elastic modulus, and the matrix uses the modified dynamic elastic modulus. The results are obtained by weighted summation to generate the composite elastic modulus parameters to be integrated;

[0118] Main direction elastic modulus (E1) calculation principle

[0119] The main direction (fiber axial direction) elastic modulus calculation uses the parallel mixing rule:

[0120] Fiber contribution: fiber static elastic modulus (such as 230 GPa) x modified fiber volume fraction V_f';

[0121] Interphase contribution: interphase static elastic modulus (such as 8 GPa, measured by nanoindentation) x interphase volume fraction V_i;

[0122] Matrix contribution: extract the matrix dynamic elastic modulus at the current strain rate (such as 3.2 GPa at a strain rate of 100 s -1 -1) from the modified matrix dynamic performance parameter set x the modified matrix volume fraction V_m';

[0123] Main direction modulus:

[0124] E1 = Σ (modulus of each phase x volume fraction) = E_f x V_f' + E_i x V_i + E_m_dyn x V_m'; For example: E1 = 230 x 0.588 + 8 x 0.02 + 3.2 x 0.392 ≈ 138.2 GPa.

[0125] Secondary direction elastic modulus (E2) calculation principle

[0126] The secondary direction (perpendicular to the fiber) elastic modulus uses the series mixing rule:

[0127] Sum of each phase compliance: compliance is the inverse of modulus;

[0128] Fiber compliance = 1 / fiber static modulus x V_f';

[0129] Interphase compliance = 1 / interphase static modulus x V_i;

[0130] Matrix compliance = 1 / matrix dynamic modulus x V_m';

[0131] E2 = 1 / (sum of compliances) = 1 / (S_f x V_f' + S_i x V_i + S_m_dyn x V_m');

[0132] Example: E2 = 1 / (1 / 230 x 0.588 + 1 / 8 x 0.02 + 1 / 3.2 x 0.392) = 5.8 GPa.

[0133] Strain rate dependent modulus parameter integration

[0134] Iterate through a pre-defined strain rate range (0.001 s -1 to 1000 s -1 ), repeat the above calculation at each strain rate node. Generate composite elastic modulus parameters to be integrated.

[0135] Calculate shear modulus according to the same mixing rule, where fiber and interphase shear modulus take static values, and matrix shear modulus takes modified dynamic values, weighted according to the volume fraction of the three phases, to generate composite shear modulus parameters;

[0136] Shear modulus (G12) calculation model

[0137] Shear modulus characterizes the material's ability to resist in-plane shear deformation, using the improved Halpin-Tsai model:

[0138] Fiber shear modulus G_f: static value (e.g. 15 GPa);

[0139] Interphase shear modulus G_i: static value (e.g. 3 GPa, obtained by microcolumn shear test);

[0140] Matrix shear modulus G_m_dyn: extracted from the modified matrix dynamic performance parameter set (e.g. 1.2 GPa at a strain rate of 500 s -1 );

[0141] Equivalent shear modulus G12 calculation formula: G12 = G_m_dyn x [(1 + ζηV_f') / (1 - ηV_f')].

[0142] Where: η = [(G_f / G_m_dyn) - 1] / [(G_f / G_m_dyn) + ζ]; ζ is the fiber shape factor (1.0 for long fibers);

[0143] Interphase modification: the interphase volume is considered as part of the matrix, and the equivalent matrix modulus G_m_eq = (G_i x V_i + G_m_dyn x V_m') / (V_i + V_m').

[0144] Dynamic shear modulus calculation example

[0145] At strain rate 200 s -1 For example:

[0146] The matrix dynamic shear modulus G_m_dyn = 1.5 Giga Pascals GPa;

[0147] The equivalent matrix modulus G_m_eq = (3x0.02 + 1.5x0.392) / (0.02 + 0.392) ≈ 1.62 Giga Pascals GPa;

[0148] Calculate η = [(15 / 1.62)-1] / [(15 / 1.62)+1] ≈ 0.82;

[0149] G12 = 1.62x[(1+1x0.82x0.588) / (1-0.82x0.588)] ≈ 4.3 Giga Pascals GPa.

[0150] Full strain rate range parameter generation

[0151] Take logarithmic coordinate points at 10 times intervals (0.001, 0.01, 0.1, 1, 10, 100, 1000 s -1 ), calculate G12 value point by point, and generate composite material shear modulus parameter table.

[0152] Calculate dynamic tensile, compression and shear strength respectively, wherein the fiber and interface phase take static strength parameters, the matrix takes modified dynamic strength parameters, and the three-phase volume fraction is weighted and summed, and the composite material elastic modulus parameters and the composite material shear modulus parameters to be integrated are integrated, to obtain the composite material dynamic performance parameter set.

[0153] Dynamic tensile strength (σ_t) calculation

[0154] Three-phase weighted model using maximum stress criterion:

[0155] Fiber tensile strength σ_f_t: static value (such as 4900 Mega Pascals MPa);

[0156] Interface phase tensile strength σ_i_t: static value (such as 80 Mega Pascals MPa, obtained by interface peeling test);

[0157] Matrix tensile strength σ_m_t_dyn: extracted from the modified matrix dynamic ultimate strength sequence (such as 120 Mega Pascals MPa at strain rate 50 s -1 );

[0158] Composite tensile strength: σ_t = σ_f_t x V_f' + σ_i_t x V_i + σ_m_t_dyn x V_m'.

[0159] For example: σ_t = 4900 x 0.588 + 80 x 0.02 + 120 x 0.392 ≈ 2940 MPa.

[0160] Dynamic compressive strength (σ_c) and shear strength (τ) calculation

[0161] Compression strength needs to consider the matrix support effect: σ_c = σ_f_c x V_f' + k x σ_m_c_dyn x V_m'.

[0162] Wherein:

[0163] σ_f_c: fiber compression strength static value (such as 4200 MPa) ;

[0164] σ_m_c_dyn: matrix dynamic compression strength (different from tensile strength, needs to be corrected separately).

[0165] k is the support coefficient (empirical value 1.2-1.5, determined by fiber arrangement)

[0166] Shear strength adopts interface dominant model: τ = τ_i x (V_i)^0.5 + τ_m_dyn x V_m'. Wherein:

[0167] τ_i: interface static shear strength (extracted from ternary parameter set, such as 60 MPa) ;

[0168] τ_m_dyn: matrix dynamic shear strength (such as 90 MPa at strain rate 200 s -1 ).

[0169] Multi-parameter integration and structured output

[0170] Calculate all strength parameters for each strain rate node: tensile strength σ_t; compression strength σ_c; shear strength τ.

[0171] Integrate elastic parameters: main / secondary direction modulus E1, E2; shear modulus G12.

[0172] Build composite material dynamic performance parameter set, including parameters: strain rate (s -1 ), E1 (GPa), E2 (GPa), G12 (GPa), σ_t (MPa), σ_c (MPa), τ (MPa). This parameter set can be directly input into finite element software (such as Abaqus) as material card (Material Card) for high strain rate impact simulation.

[0173] S204, compare the composite material dynamic performance parameter set with the multi-strain rate experimental data, optimize the interface coefficient weight in the ternary parameter set through interface parameter sensitivity analysis, construct a strain rate-interface coupling performance model as the simulation input under high strain rate loading conditions.

[0174] Specifically, the composite material dynamic performance parameter set can be compared with the experimental test data under different strain rates point by point, the strength and modulus error values corresponding to each strain rate are calculated, and a performance error sequence is generated.

[0175] After completing the theoretical calculation of the composite material dynamic performance parameter set, its accuracy needs to be verified through experimental data. The experimental data comes from standardized dynamic testing covering the target strain rate range (such as 0.001 to 1000 seconds per minute). Dynamic tensile, compression and shear experiments are carried out using a split Hopkinson pressure bar (SHPB) device:

[0176] Dynamic tensile test: the carbon fiber composite sample is clamped between the incident bar and the transmission bar, and the stress pulse is generated by driving the impact bar to impact the incident bar through high-pressure gas. The strain gauge attached to the bar body is used to record the pulse waveform, and the strain rate (SR) and dynamic tensile strength (DTS) of the sample are calculated based on one-dimensional stress wave theory.

[0177] Dynamic compression test: cylindrical samples are placed in the SHPB device, and dynamic compression strength (DCS) is obtained through the same principle.

[0178] Dynamic shear test: interlaminar shear samples with notches or thin-walled cylindrical samples are used to realize pure shear loading on the SHPB through special clamps, and dynamic shear strength (DSS) is measured.

[0179] At the same time, the deformation field of the sample is recorded synchronously through high-speed photography or digital image correlation technology (DIC), which assists in verifying the measured values of elastic modulus (EM) and shear modulus (SM). At least three repeated experiments are carried out under each strain rate, and the average value after removing abnormal values is taken as the final experimental data point.

[0180] The Composite Dynamic Performance Parameter Set (CDPPS) generated from the theoretical calculation contains the following key parameters: principal direction elastic modulus (E11), secondary direction elastic modulus (E22), shear modulus (G12), dynamic tensile strength (σ_t), dynamic compressive strength (σ_c), dynamic shear strength (τ_s). The theoretical prediction value is given for each parameter at strain rate intervals (e.g. 10 points in logarithmic intervals) from 0.001 sec-1 to 1000 sec-1.

[0181] The point-by-point comparison process is as follows:

[0182] Extract the experimental measurement value (Experimental Value, EV) and the theoretical prediction value (Theoretical Value, TV) at the same strain rate.

[0183] Calculate the absolute error value (Absolute Error Value, AEV) for each parameter:

[0184] Strength parameter error: AEV_σ = |EV_σ - TV_σ|;

[0185] Modulus parameter error: AEV_E = |EV_E - TV_E|.

[0186] Further calculate the relative error percentage (Relative Error Percentage, REP) to standardize the error: REP = (AEV / EV) x 100%.

[0187] For example, at a strain rate of 500 sec-1, if the experimental dynamic tensile strength is 1200 MPa and the theoretical prediction value is 1150 MPa, the absolute error is 50 MPa and the relative error is 4.17%. This process is repeated for all parameters at all strain rate points.

[0188] Arrange the error calculation results (including absolute error and relative error) of all parameters at each strain rate point in ascending order of strain rate to form the Performance Error Sequence (PES). This sequence uses a structured data storage format and contains the following fields:

[0189] Strain Rate (SR, unit: reciprocal second); Parameter Type (PT, such as E11, σ_t, etc.); Experimental Value (EV); Theoretical Value (TV); Absolute Error (AE); Relative Error (RE).

[0190] For example, at the strain rate point of 100 reciprocal second, the PES may record: "SR = 100 s -1 -1, PT = σ_t, EV = 1100 MPa, TV = 1050 MPa, AE = 50 MPa, RE = 4.55%". This sequence is the core input of subsequent sensitivity analysis and parameter optimization, and needs to ensure data integrity and traceability.

[0191] Sensitivity analysis is performed on the performance error sequence. By adjusting the interfacial shear strength and thickness coefficient separately through the controlled variable method, the change range of the error value is observed to determine the influence weight of the two parameters on the error, and an interfacial parameter sensitivity weight table is generated.

[0192] Sensitivity analysis (Sensitivity Analysis, SA) aims to quantify the influence degree of interfacial parameters (interfacial shear strength and interfacial thickness coefficient) on the theoretical prediction error. The controlled variable method (Controlled Variable Method, CVM) is used for operation:

[0193] Baseline parameter set: The initial three-parameter set (including fiber static parameters, matrix strain rate sensitive parameters, and interfacial dynamic shear characteristic parameters) is taken as the baseline (Baseline Set, BS).

[0194] Parameter perturbation rule: Interfacial shear strength (Interfacial Shear Strength, ISS): perturbed by ±10%, ±20% amplitude based on the baseline value, generating four perturbation values (such as ISS_+10%, ISS_+20%, ISS_-10%, ISS_-20%).

[0195] Interfacial thickness coefficient (Interfacial Thickness Coefficient, ITC): similarly perturbed by ±10%, ±20% amplitude, generating four perturbation values (ITC_+10% and the like).

[0196] Perturbation isolation principle: each time only one parameter is perturbed (such as adjusting only ISS), and the rest of the parameters (including ITC) remain unchanged at the baseline value, ensuring single source of influence.

[0197] For each perturbed parameter combination, recalculate the flow:

[0198] Using the perturbed ISS or ITC value, recalculate the Modified Matrix Dynamic Performance Parameter Set (MMDPS).

[0199] Based on the MMDPS, regenerate the Composite Dynamic Performance Parameter Set (New CDPPS) through the three-phase synergistic mixing rule.

[0200] Compare the New CDPPS with experimental data again point by point to generate a new Performance Error Sequence (New PES).

[0201] Key operation examples:

[0202] When ISS increases by 20% (ITC remains unchanged), the theoretical dynamic tensile strength at a strain rate of 1000 seconds increases from 1050 MPa to 1080 MPa, with an absolute error of 50 MPa and a relative error of 4.55% compared to the original experimental value of 1100 MPa, which decreases to 20 MPa and 1.82%, respectively.

[0203] When ITC decreases by 20% (ISS remains unchanged), the theoretical dynamic tensile strength at the same strain rate decreases to 1030 MPa, with an absolute error of 70 MPa and a relative error of 6.36% increase.

[0204] By statistically analyzing the error value change trend under different perturbations, calculate the Sensitivity Weight (SW):

[0205] Change amplitude calculation: for each strain rate point, record the relative error change ΔRE caused by parameter perturbation (e.g. ΔRE = RE_original - RE_new when ISS + 20%).

[0206] Take the average of ΔRE for all strain rate points to get the Mean Error Variation (MEV) for this perturbation amplitude.

[0207] Weight distribution: the Sensitivity Weight of ISS SW_ISS is proportional to MEV_ISS (the average MEV caused by ISS perturbation); the Sensitivity Weight of ITC SW_ITC is proportional to MEV_ITC (the average MEV caused by ITC perturbation).

[0208] Normalization: make SW_ISS + SW_ITC = 1.

[0209] An Interfacial Parameter Sensitivity Weight Table (IPSWT) is finally generated.

[0210] According to the sensitivity weight table, the interfacial shear strength and thickness coefficient in the ternary parameter set are adjusted according to the error minimization principle, and iteratively optimized until the error value is lower than the preset threshold, to generate an optimized ternary parameter set.

[0211] The iterative optimization process is based on the sensitivity weight table (IPSWT) and the error minimization objective:

[0212] Initial adjustment strategy:

[0213] A larger adjustment step (such as ±5%) is used for high-weight parameters (such as ISS, weight 0.65), and a smaller step (such as ±2%) is used for low-weight parameters (such as ITC, weight 0.35).

[0214] Direction determination mechanism:

[0215] If increasing ISS can reduce the error (as in the example above), then preferentially attempt to adjust ISS in the positive direction; if reducing ITC will increase the error, then avoid adjusting ITC in the negative direction.

[0216] First round adjustment example: initial ISS = 80 MPa, ITC = 0.15.

[0217] According to the sensitivity weight, ISS is increased by 5% to 84 MPa, and ITC remains unchanged.

[0218] Recalculate CDPPS and new error sequence, find that the average relative error is reduced from 8.2% to 7.0%.

[0219] Multiple iterations and convergence control:

[0220] Error evaluation criteria: take the mean relative error (MRE) as the core indicator, and set the preset threshold to 5% (which can be set according to actual needs).

[0221] Adaptive step size adjustment:

[0222] If the MRE decreases by more than 1% after this round of adjustment, then maintain the same adjustment direction and step size in the next round;

[0223] If the MRE decreases by less than 0.5%, then halve the step size (e.g. ISS adjustment step size from 5% to 2.5%);

[0224] If the MRE increases, then adjust the parameters in the opposite direction (e.g. change from increasing to decreasing ISS).

[0225] Loop termination condition:

[0226] Condition 1: MRE ≤ 5% (reached preset threshold);

[0227] Condition 2: MRE variation amplitude is less than 0.2% in three consecutive iterations (converged stably).

[0228] Generate Optimized Ternary Parameter Set (OTPS):

[0229] Parameter update: replace the corresponding values in the original ternary parameter set with the final converged ISS (91 MPa) and ITC (0.147), and keep the rest of the parameters (such as fiber static modulus, matrix Cowper-Symonds coefficient) unchanged.

[0230] Data verification: recalculate the composite performance in the full strain rate range using OTPS; compare the new predicted value with the experimental data to confirm that MRE is stable at 4.9% (lower than the 5% threshold).

[0231] Documented output: OTPS is stored in a structured database (such as JSON or XML format), containing fields:

[0232] Fiber parameter group (static modulus, strength, etc.);

[0233] Matrix parameter group (strain rate dependent coefficient C, sensitivity index p, etc.);

[0234] Interface parameter group (optimized ISS, ITC, and strain rate correlation table);

[0235] Optimization metadata (number of iterations, final MRE, convergence criterion).

[0236] This optimized parameter set is the basis for building a high-precision strain rate-interface coupling model.

[0237] Based on the optimized ternary parameter set, recalculate the performance parameters of the composite material in the full strain rate range, and build a strain rate-interface coupling performance model as the simulation input under high strain rate loading conditions.

[0238] Full strain rate performance calculation:

[0239] Input parameters: use all data in the optimized ternary parameter set (OTPS), especially the optimized interface shear strength (ISS = 91 MPa) and interface thickness coefficient (ITC = 0.147).

[0240] Calculation flow reproduction:

[0241] Step 1: Re-calculate the dynamic property parameter set of the matrix (including the interface-corrected Cowper-Symonds model and the logarithmic modulus formula) according to the interface parameters in the OTPS.

[0242] Step 2: Calculate the dynamic property parameter set of the composite (New CDPPS) based on the new matrix parameters through the three-phase synergistic mixing rule.

[0243] Data coverage: The strain rate range is extended to 0.001 second to 5000 seconds (covering higher strain rate conditions), with points taken at 0.1 logarithmic intervals (such as 0.001, 0.01, 0.1, 1, 10, 100, 1000, 5000 seconds), ensuring data density in key intervals.

[0244] Strain rate-interface coupling performance model construction:

[0245] Model structure definition:

[0246] Independent variable: Strain rate ( Unit: seconds);

[0247] Dependent variable: Composite performance parameters (E11, E22, G12, σ_t, σ_c, τ_s);

[0248] Implicit variable: Optimized interface parameters Both are functions of strain rate (provided by the dynamic correlation table in the OTPS).

[0249] Using a combination of piecewise functions and response surface fitting:

[0250] Low strain rate segment Performance changes slowly, fitted with linear or polynomial;

[0251] High strain rate segment Performance is significantly nonlinear, using a physically-based exponential form (such as ).

[0252] Each performance parameter is independently established as a sub-model, for example:

[0253] Dynamic tensile strength model:

[0254] Principal direction modulus model:

[0255] Model verification: Plot the theoretical data of New CDPPS and experimental data on the same strain rate-performance coordinate graph, and verify the model accuracy through the coefficient of determination (R 2 >0.98). ​

[0256] Implementation as high strain rate simulation input:

[0257] Model integration approach: Embed the constructed Strain Rate-Interface Coupled Performance Model (SR-ICPM) into the User Material Subroutine (UMAT / VUMAT) of commercial finite element software (e.g. Abaqus, LS-DYNA).

[0258] In the subroutine, according to the real-time strain rate of the current integration point (provided by the solver), the SR-ICPM is called to calculate the instantaneous material parameters (such as elastic modulus, yield strength).

[0259] Typical application scenarios:

[0260] Ballistic impact simulation: A projectile hits a carbon fiber composite armor plate at 1500 meters per second (m / s), with a local strain rate exceeding 3000 per second. The SR-ICPM automatically calls high strain rate parameters of 5000 per second to predict material response.

[0261] Turbine blade containment analysis: After the engine blade breaks, it hits the casing at high speed, with a strain rate of about 2000 per second. The model predicts the dynamic anti-penetration ability of the composite casing based on optimized interface parameters.

[0262] Engineering value:

[0263] By accurately reflecting the evolution of interface parameters (ISS, ITC) with strain rate, the simulation prediction reliability of high-speed collision, explosion impact and other scenarios is significantly improved.

[0264] The optimized ternary parameter set (OTPS) can be reused to the same type of composite material system, reducing the cost of experimental calibration.

[0265] Finally, the model is delivered in the form of dynamic link library (DLL) or source code, becoming a core simulation tool for composite material structure design under high strain rate conditions.

[0266] It can be seen that the static mechanical parameters of the carbon fiber, the strain rate sensitive parameters of the matrix material and the dynamic shear characteristic parameters of the fiber-matrix interface are acquired, a ternary parameter set of the fiber-matrix-interface is constructed, based on the ternary parameter set, the Cowper-Symonds model introducing an interface shear correction term is used to calculate the dynamic ultimate strength of the matrix, a corrected dynamic performance parameter set of the matrix covering a preset strain rate range is generated, according to the corrected dynamic performance parameter set of the matrix, the performance of the composite material is calculated through a three-phase synergistic mixing rule, a dynamic performance parameter set of the composite material is output, the dynamic performance parameter set of the composite material is compared with the multi-strain rate experimental data, a strain rate-interface coupling performance model is constructed, and thus the prediction accuracy of the dynamic mechanical performance of the composite material can be improved by constructing the ternary parameter set of the fiber-matrix-interface and the strain rate-interface coupling model.

[0267] Another embodiment of the present application provides a carbon fiber composite material mechanical performance adjustment system based on strain rate effect, referring to Figure 3 , the system can include:

[0268] The acquisition module 301 is configured to acquire the static mechanical parameters of the carbon fiber, the strain rate sensitive parameters of the matrix material and the dynamic shear characteristic parameters of the fiber-matrix interface, and construct a ternary parameter set of the fiber-matrix-interface, wherein the dynamic shear characteristic parameters include the interface shear strength and thickness coefficient under different strain rates.

[0269] The derivation module 302 is configured to, based on the ternary parameter set, use the Cowper-Symonds model introducing an interface shear correction term to calculate the dynamic ultimate strength of the matrix, derive the elastic modulus of the matrix in combination with a logarithmic relationship formula containing the interface thickness coefficient, and generate a corrected dynamic performance parameter set of the matrix covering a preset strain rate range.

[0270] The calculation module 303 is configured to, according to the corrected dynamic performance parameter set of the matrix, calculate the performance of the composite material through a three-phase synergistic mixing rule, regard the composite material as a three-phase system of fiber, matrix and interface, calculate the main and auxiliary direction elastic modulus, shear modulus and dynamic tensile, compression and shear strength according to the volume fraction, and output a dynamic performance parameter set of the composite material.

[0271] The construction module 304 is configured to compare the dynamic performance parameter set of the composite material with the multi-strain rate experimental data, optimize the interface coefficient weight in the ternary parameter set through interface parameter sensitivity analysis, construct a strain rate-interface coupling performance model, and use the model as the simulation input under the high strain rate loading condition.

[0272] The embodiment of the present application also provides a storage medium, and the storage medium stores a computer program, wherein the computer program is set to execute the steps in any of the method embodiments when running.

[0273] Specifically, in the embodiment, the storage medium can be configured to store a computer program for executing the following steps:

[0274] S201, obtain static mechanical parameters of carbon fibers, strain rate sensitive parameters of a matrix material, and dynamic shear characteristic parameters of a fiber-matrix interface, and construct a ternary parameter set of the fiber-matrix-interface, wherein the dynamic shear characteristic parameters include interface shear strength and thickness coefficients under different strain rates;

[0275] S202, based on the ternary parameter set, calculate the dynamic ultimate strength of the matrix by using a Cowper-Symonds model introducing an interface shear correction term, derive the elastic modulus of the matrix by combining a logarithmic relationship formula containing the interface thickness coefficient, and generate a corrected matrix dynamic performance parameter set covering a preset strain rate range;

[0276] S203, according to the corrected matrix dynamic performance parameter set, calculate the performance of the composite material by a three-phase synergistic mixing rule, regard the composite material as a three-phase system of fibers, matrixes, and interfaces, calculate the main and auxiliary direction elastic moduli, shear moduli, and dynamic tensile, compression, and shear strengths by volume fraction weighting, and output a composite material dynamic performance parameter set;

[0277] S204, compare the composite material dynamic performance parameter set with multi-strain rate experimental data, optimize the interface coefficient weight in the ternary parameter set through interface parameter sensitivity analysis, construct a strain rate-interface coupling performance model as a simulation input under a high strain rate loading condition.

[0278] The embodiment of the application further provides an electronic device including a memory and a processor, the memory stores a computer program, and the processor is configured to run the computer program to execute the steps in any one of the method embodiments.

[0279] Specifically, the electronic device can further include a transmission device and an input-output device, wherein the transmission device is connected with the processor, and the input-output device is connected with the processor.

[0280] Specifically, in the embodiment, the processor can be configured to execute the following steps by the computer program:

[0281] S201, obtain static mechanical parameters of carbon fibers, strain rate sensitive parameters of a matrix material, and dynamic shear characteristic parameters of a fiber-matrix interface, and construct a ternary parameter set of the fiber-matrix-interface, wherein the dynamic shear characteristic parameters include interface shear strength and thickness coefficients under different strain rates;

[0282] S202, based on the three parameter set, the Cowper-Symonds model introducing the interface shear correction term is used to calculate the dynamic ultimate strength of the matrix, the elastic modulus of the matrix is derived by combining the logarithmic relationship formula containing the interface thickness coefficient, and the corrected matrix dynamic performance parameter set covering the preset strain rate range is generated;

[0283] S203, according to the corrected matrix dynamic performance parameter set, the performance of the composite material is calculated by the three-phase synergistic mixing rule, the composite material is regarded as a three-phase system of fiber, matrix and interface phase, the main and auxiliary direction elastic modulus, shear modulus and dynamic tensile, compression and shear strength are calculated by volume fraction weighting, and the dynamic performance parameter set of the composite material is output;

[0284] S204, the dynamic performance parameter set of the composite material is compared with the multi-strain rate experimental data, the interface coefficient weight in the three parameter set is optimized through the interface parameter sensitivity analysis, the strain rate-interface coupling performance model is constructed, and is used as the simulation input under the high strain rate loading condition.

[0285] The above describes the structure, features and effects of the application in detail according to the embodiments shown in the drawings. The above description is only the preferred embodiments of the application, but the application is not limited by the drawings. Any change or modification within the scope of the application, or equivalent embodiments with equivalent changes, shall be within the protection scope of the application.

Claims

1. A method for adjusting the mechanical properties of carbon fiber composite materials based on strain rate effect, characterized in that, The method includes: The static mechanical parameters of carbon fiber, the strain rate sensitive parameters of the matrix material, and the dynamic shear characteristics of the fiber-matrix interface are obtained to construct a ternary parameter set of fiber-matrix-interface. The dynamic shear characteristics parameters include the interfacial shear strength and thickness coefficient under different strain rates. Based on the ternary parameter set, the dynamic ultimate strength of the matrix is ​​calculated using the Cowper-Symonds model with an interface shear correction term. The elastic modulus of the matrix is ​​derived by combining the logarithmic relationship formula including the interface thickness coefficient, and a modified matrix dynamic performance parameter set covering the preset strain rate range is generated. Based on the modified matrix dynamic performance parameter set, the composite material properties are calculated using the three-phase synergistic mixing rule. The composite material is regarded as a three-phase system of fiber, matrix, and interface phase. The elastic modulus, shear modulus, and dynamic tensile, compressive, and shear strength in the principal and secondary directions are calculated by weighting by volume fraction, and the composite material dynamic performance parameter set is output. The dynamic performance parameter set of the composite material is compared with the multi-strain rate experimental data. The interface coefficient weights in the ternary parameter set are optimized through interface parameter sensitivity analysis to construct a strain rate-interface coupled performance model, which serves as the simulation input under high strain rate loading conditions.

2. The method according to claim 1, characterized in that, The static mechanical parameters of the carbon fiber, the strain rate-sensitive parameters of the matrix material, and the dynamic shear characteristics of the fiber-matrix interface are obtained to construct a ternary parameter set for the fiber-matrix-interface. The dynamic shear characteristics include the interfacial shear strength and thickness coefficient at different strain rates, including: A constant-temperature static tensile testing apparatus was used, and the test was conducted at 0.001 s. -1 The carbon fiber monofilament was tested under strain rate, and the elastic modulus, ultimate strength and elongation at break were recorded. The average value of three parallel experiments was taken to generate a table of static mechanical parameters of carbon fiber. Using a split Hopkinson pressure bar device, at 0.001s -1 up to 1000s -1 Dynamic compression experiments were conducted on the matrix material within the strain rate range. Data on the changes in yield strength and elastic modulus under different strain rates were collected. The strain rate correlation coefficient C, strain rate sensitivity index p, and elastic modulus sensitivity coefficient k in the Cowper-Symonds model parameters were obtained by fitting, and a table of matrix strain rate sensitivity parameters was generated. Dynamic testing was conducted at the fiber-matrix interface using the micropillar indentation method. Different strain rates were simulated by controlling the indenter loading rate, and the relationship curve between indentation depth and shear load was recorded. The interfacial shear strength was calculated, and the microstructure of the interface was observed using scanning electron microscopy. The interfacial thickness coefficient was extracted, and a table of dynamic shear characteristic parameters of the interface was generated. The carbon fiber static mechanical parameter table, matrix strain rate sensitive parameter table, and interface dynamic shear characteristic parameter table are correlated according to the fiber-matrix volume fraction to establish the mapping relationship between the parameters and construct a ternary parameter set of fiber-matrix-interface.

3. The method according to claim 2, characterized in that, Based on the ternary parameter set, the dynamic ultimate strength of the matrix is ​​calculated using the Cowper-Symonds model with an interface shear correction term. The elastic modulus of the matrix is ​​derived using a logarithmic formula including the interface thickness coefficient, generating a corrected set of matrix dynamic performance parameters covering a preset strain rate range, including: The interfacial shear strength parameter is extracted from the ternary parameter set, and its ratio with the static ultimate strength of the matrix is ​​calculated as the interfacial shear correction coefficient. This coefficient is dynamically adjusted according to a linear relationship as the strain rate increases, generating an interfacial shear correction coefficient sequence. Substitute the interface shear correction coefficient sequence into the Cowper-Symonds model, adjust the dynamic limit strength calculation results by adding correction terms, obtain the matrix dynamic limit strength under different strain rates, and generate the corrected matrix dynamic limit strength sequence. The interface thickness coefficient is extracted from the ternary parameter set and used as a multiplier in the logarithmic relationship formula of the matrix elastic modulus. The change of elastic modulus with strain rate is adjusted by the thickness coefficient, and the elastic modulus under different strain rates is calculated to generate the corrected matrix elastic modulus sequence. The integrated and corrected matrix dynamic ultimate strength sequence and elastic modulus sequence are calculated at 0.001s. -1 up to 1000s -1 The strain rate intervals are used to complete the data and generate a set of corrected matrix dynamic performance parameters covering the preset strain rate range.

4. The method according to claim 3, characterized in that, The composite material properties are calculated based on the modified matrix dynamic performance parameter set using the three-phase synergistic mixing principle. The composite material is considered as a three-phase system of fiber, matrix, and interface phases. The elastic modulus, shear modulus, and dynamic tensile, compressive, and shear strength in the principal and secondary directions are calculated using volume fraction weighting. The resulting composite material dynamic performance parameter set includes: The interfacial phase volume of a single fiber is calculated based on the fiber diameter and interfacial thickness coefficient. Combining the fiber volume fraction and matrix volume fraction, the volume fraction ratios of the fiber, matrix, and interfacial phases are derived, and a three-phase volume fraction distribution table is generated. Based on the three-phase volume fraction distribution table, the elastic modulus of the principal and secondary directions is calculated using the three-phase synergistic mixing rule. The static elastic modulus of the fiber and interface phases is used, and the modified dynamic elastic modulus of the matrix is ​​used. The results are obtained by weighted summation, and the elastic modulus parameters of the composite material to be integrated are generated. The shear modulus is calculated according to the same mixing rule, where the shear modulus of the fiber and the interface phase is taken as the static value, and the shear modulus of the matrix is ​​taken as the corrected dynamic value. The shear modulus parameters of the composite material are generated by weighting the three phase volume fractions. Dynamic tensile, compressive, and shear strengths were calculated separately. Static strength parameters were used for the fiber and interfacial phases, while modified dynamic strength parameters were used for the matrix. The three phases were weighted and summed according to their volume fractions. The elastic modulus parameters and shear modulus parameters of the composite material to be integrated were also combined to obtain the set of dynamic performance parameters of the composite material.

5. The method according to claim 4, characterized in that, The process of comparing the dynamic performance parameter set of the composite material with multi-strain rate experimental data, optimizing the interface coefficient weights in the ternary parameter set through interface parameter sensitivity analysis, and constructing a strain rate-interface coupled performance model as the simulation input under high strain rate loading conditions includes: The dynamic performance parameter set of composite materials is compared point by point with the experimental test data at different strain rates. The strength and modulus error values ​​corresponding to each strain rate are calculated to generate a performance error sequence. Sensitivity analysis was performed on the performance error sequence. The interface shear strength and thickness coefficient were adjusted separately by the control variable method. The change range of the error value was observed, the influence weight of the two parameters on the error was determined, and an interface parameter sensitivity weight table was generated. Based on the sensitivity weight table, the interface shear strength and thickness coefficient in the ternary parameter set are adjusted according to the principle of minimizing error. The optimization is iterated until the error value is lower than the preset threshold, and the optimized ternary parameter set is generated. Based on the optimized ternary parameter set, the performance parameters of the composite material in the full strain rate range were recalculated, and a strain rate-interface coupled performance model was constructed as the simulation input under high strain rate loading conditions.

6. A system for adjusting the mechanical properties of carbon fiber composite materials based on strain rate effect, characterized in that, The system includes: The acquisition module is used to acquire the static mechanical parameters of carbon fiber, the strain rate sensitive parameters of the matrix material, and the dynamic shear characteristics of the fiber-matrix interface, and to construct a ternary parameter set of fiber-matrix-interface. The dynamic shear characteristics parameters include the interfacial shear strength and thickness coefficient under different strain rates. The derivation module is used to calculate the dynamic ultimate strength of the matrix based on the ternary parameter set, using the Cowper-Symonds model with an interface shear correction term, and derive the elastic modulus of the matrix by combining the logarithmic relationship formula including the interface thickness coefficient, and generate a modified matrix dynamic performance parameter set covering a preset strain rate range. The calculation module is used to calculate the properties of composite materials based on the modified matrix dynamic performance parameter set and the three-phase synergistic mixing rule. The composite material is regarded as a three-phase system of fiber, matrix and interface phase. The elastic modulus, shear modulus and dynamic tensile, compressive and shear strength in the main and secondary directions are calculated by volume fraction weighting, and the composite material dynamic performance parameter set is output. The module is used to compare the dynamic performance parameter set of the composite material with multi-strain rate experimental data, optimize the interface coefficient weights in the ternary parameter set through interface parameter sensitivity analysis, and construct a strain rate-interface coupled performance model as the simulation input under high strain rate loading conditions.

7. The system according to claim 6, characterized in that, The acquisition module is specifically used for: A constant-temperature static tensile testing apparatus was used, and the test was conducted at 0.001 s. -1 The carbon fiber monofilament was tested under strain rate, and the elastic modulus, ultimate strength and elongation at break were recorded. The average value of three parallel experiments was taken to generate a table of static mechanical parameters of carbon fiber. Using a split Hopkinson pressure bar device, at 0.001s -1 up to 1000s -1 Dynamic compression experiments were conducted on the matrix material within the strain rate range. Data on the changes in yield strength and elastic modulus under different strain rates were collected. The strain rate correlation coefficient C, strain rate sensitivity index p, and elastic modulus sensitivity coefficient k in the Cowper-Symonds model parameters were obtained by fitting, and a table of matrix strain rate sensitivity parameters was generated. Dynamic testing was conducted at the fiber-matrix interface using the micropillar indentation method. Different strain rates were simulated by controlling the indenter loading rate, and the relationship curve between indentation depth and shear load was recorded. The interfacial shear strength was calculated, and the microstructure of the interface was observed using scanning electron microscopy. The interfacial thickness coefficient was extracted, and a table of dynamic shear characteristic parameters of the interface was generated. The carbon fiber static mechanical parameter table, matrix strain rate sensitive parameter table, and interface dynamic shear characteristic parameter table are correlated according to the fiber-matrix volume fraction to establish the mapping relationship between the parameters and construct a ternary parameter set of fiber-matrix-interface.

8. The system according to claim 7, characterized in that, The derivation module is specifically used for: The interfacial shear strength parameter is extracted from the ternary parameter set, and its ratio with the static ultimate strength of the matrix is ​​calculated as the interfacial shear correction coefficient. This coefficient is dynamically adjusted according to a linear relationship as the strain rate increases, generating an interfacial shear correction coefficient sequence. Substitute the interface shear correction coefficient sequence into the Cowper-Symonds model, adjust the dynamic limit strength calculation results by adding correction terms, obtain the matrix dynamic limit strength under different strain rates, and generate the corrected matrix dynamic limit strength sequence. The interface thickness coefficient is extracted from the ternary parameter set and used as a multiplier in the logarithmic relationship formula of the matrix elastic modulus. The change of elastic modulus with strain rate is adjusted by the thickness coefficient, and the elastic modulus under different strain rates is calculated to generate the corrected matrix elastic modulus sequence. The integrated and corrected matrix dynamic ultimate strength sequence and elastic modulus sequence are calculated at 0.001s. -1 up to 1000s -1 The strain rate intervals are used to complete the data and generate a set of corrected matrix dynamic performance parameters covering the preset strain rate range.

9. A storage medium, characterized in that, The storage medium stores a computer program, wherein the computer program is configured to execute the method of any one of claims 1-5 when it is run.

10. An electronic device comprising a memory and a processor, characterized in that, The memory stores a computer program, and the processor is configured to run the computer program to perform the method of any one of claims 1-5.

Images