Method and device for identifying J-C constitutive parameters of manganese-aluminum bronze

Abstract

The invention belongs to the technical field of material modeling and simulation, and discloses an identification method and device for J-C constitutive parameters of manganese-aluminum bronze, and the method comprises the steps: obtaining a real stress-strain curve of a material through a quasi-static compression experiment, and carrying out a thermal-mechanical coupling compression experiment under the conditions of different temperatures and strain rates, extracting stress-strain response in a corresponding state, and estimating an initial parameter set of the model; constructing stress-strain response verification data through a thermal-mechanical coupling compression experiment; constructing a parameter optimization space based on the initial parameter set, and performing iterative optimization based on a fusion optimization algorithm of a genetic algorithm, a grey wolf algorithm and a particle swarm algorithm; the fitness function is weighted based on a variable weight on the basis of an average relative error of the experimental data curve and the fitting data curve, and the variable weight is determined according to the fitting deviation of the yield point and the plastic section. The J-C constitutive model obtained by the invention can accurately simulate the mechanical behavior of manganese-aluminum bronze.

Description

TECHNICAL FIELD

[0001] The present application belongs to the technical field of material modeling and simulation, and particularly relates to a method and device for identifying J-C constitutive parameters of manganese aluminum bronze. BACKGROUND

[0002] The statements in this section merely provide background information related to the present application and do not necessarily constitute prior art.

[0003] The J-C constitutive model refers to a constitutive relationship model for describing the stress-strain response of metal materials under different strain rates, temperatures and other load conditions, which can characterize the mechanical behaviors of materials such as yield, work hardening and thermal softening, and is one of the core models used in material modeling and simulation in the fields of aerospace, shipbuilding, weapons and the like. At present, for the parameter inversion in the model, the conventional technical implementation manner is to complete a series of mechanical property tests under different strain rates and different temperature conditions through a large number of grouped tensile and compressive experiments, obtain the stress-strain data of the material, and then estimate the five core parameters of the J-C constitutive model by using stepwise linear fitting or least squares method to realize the preliminary acquisition of the J-C constitutive model parameters and provide basic parameter support for material mechanical behavior simulation. However, this method relies on a large number of experiments to obtain parameters, needs to cover different strain rate and temperature combination conditions, and the experimental process needs to be controlled with high precision, resulting in a significant increase in experimental cost and a substantial extension of the parameter acquisition period. At the same time, the fitting result is extremely sensitive to the accuracy of the experimental data, and experimental errors are easily directly transmitted to the parameter results, resulting in poor technical and economic performance, and the parameter accuracy is limited by the experimental conditions, which is difficult to meet the efficient and high-precision engineering application requirements.

[0004] Another type of existing technology is to obtain material macro stress-strain response data under typical loading conditions, and perform parameter inversion based on genetic algorithm and other optimization algorithms to replace the traditional large number of grouped experimental fitting method, which can reduce the experimental cost to a certain extent and improve the parameter solving efficiency. However, there is a significant nonlinear coupling relationship between the five parameters of the J-C constitutive model, and especially for manganese aluminum bronze alloys containing α, β and κ multi-phase states, the strain rate dependence and temperature sensitivity are stronger, and the parameter coupling effect is more prominent. Based on the existing optimization algorithm for optimization, it is easy to cause the optimization period to be too long or to fall into local optimization. SUMMARY

[0005] To overcome the shortcomings of the above-mentioned prior art, the present application provides a method and device for identifying J-C constitutive parameters of manganese aluminum bronze, so that the constitutive model can more accurately simulate the mechanical behavior of manganese aluminum bronze.

[0006] To achieve the above-mentioned purpose, one or more embodiments of the present application provide a method for identifying J-C constitutive parameters of manganese aluminum bronze, comprising the following steps: The real stress-strain curve of the material is obtained through a quasi-static compression experiment, and a thermal-mechanical coupling compression experiment is performed under different temperature and strain rate conditions to extract the stress-strain response in the corresponding state, and estimate the initial parameter set of the model; A thermal-mechanical coupling compression experiment is performed under a set temperature and strain rate condition to extract the stress-strain response in the corresponding state as verification data; A parameter optimization space is constructed based on the initial parameter set, and an iterative optimization is performed based on a fusion optimization algorithm of a genetic algorithm, a grey wolf algorithm and a particle swarm algorithm to obtain an optimal parameter set; the fitness function used in the optimization process is weighted based on variable weights on the basis of the average relative error of the experimental data curve and the fitted data curve, and the variable weights are determined according to the fitting deviation at the yield point and the plastic section.

[0007] In some embodiments, the process of estimating the initial parameter set of the model comprises: The strain hardening term is separated from the constitutive model, a compression experiment is performed under quasi-static and room temperature conditions to obtain a first stress-strain curve, and the yield strength, strain hardening modulus and strain hardening index are fitted based on the strain hardening term; The strain hardening term and the strain rate hardening term are separated from the constitutive model, a compression experiment is performed under high strain rate and room temperature conditions to obtain a second stress-strain curve, the yield strain and yield stress of the material are determined, and the strain rate hardening index is calculated; A compression experiment is performed under high strain rate and high temperature conditions to obtain a third stress-strain curve, the yield strain and yield stress of the material are determined, and the thermal softening index is calculated.

[0008] In some embodiments, the iterative optimization based on the fusion optimization algorithm of the genetic algorithm, the grey wolf algorithm and the particle swarm algorithm comprises: An initial population is randomly generated in a search interval, and each population individual corresponds to a set of parameters; A global coarse search is performed on the initial population based on the genetic algorithm to obtain a GA optimal individual, and a GA optimal individual set is generated; The GA optimal individual is used as an initial alpha wolf candidate, a directional convergence is performed on the GA optimal individual set based on the grey wolf algorithm to obtain a GWO optimal individual, and a GWO individual set is generated; The GWO optimal individual is used as the central initial value of the particle swarm algorithm, and part of the GA optimal individual set and the GWO optimal individual set are used as the search range for local optimization to obtain a global optimal parameter set.

[0009] In some embodiments, the fitness function used by the fusion optimization algorithm is:

[0010] Wherein, P represents a parameter combination to be optimized, The number of sample points representing experimental data, The weight of the tth iteration, the initial value is 1, The experimental data stress value of the ith sample point, k represents the optimization stage, GA, GWO and PSO are identified, and m is the starting point of the continuous sample point corresponding to the yield point and the plastic section, The stress value of the yield point and the plastic section of the fitting curve, The fine-tuning weight.

[0011] In some embodiments, in the grey wolf algorithm, the position updating formula is:

[0012] Wherein, The alpha, beta and delta wolf pack leaders are represented by, The current round of grey wolf position, The position of a randomly selected grey wolf in the current group, Used to control the approaching speed and direction, Used to control the disturbance and contraction range.

[0013] In some embodiments, the method further comprises: Substitute the optimal parameter set into the constitutive model to perform a simulation experiment, and obtain a simulation stress-strain curve; According to the average relative error of the simulation stress-strain curve and the real experiment curve, it is judged whether the precision requirement is met, if it is met, the optimal parameter set is output, if it is not met, the optimal parameter set is taken as the starting individual of the search, and the fusion optimization algorithm based on genetic algorithm, grey wolf algorithm and particle swarm optimization algorithm is iterated and optimized again.

[0014] The second aspect of the application provides a device for identifying J-C constitutive parameters of manganese aluminum bronze, comprising: An initial parameter acquisition module is configured to obtain the real stress-strain curve of the material through a quasi-static compression experiment, and perform a thermal-mechanical coupling compression experiment under different temperature and strain rate conditions, extract the stress-strain response under the corresponding state, and estimate the initial parameter set of the model; A verification data acquisition module is configured to perform a thermal-mechanical coupling compression experiment under a set temperature and strain rate condition, extract the stress-strain response under the corresponding state, and use it as verification data; The fusion optimization module is configured to construct a parameter optimization space based on the initial parameter set, and to perform iterative optimization based on a fusion optimization algorithm of a genetic algorithm, a grey wolf algorithm and a particle swarm algorithm to obtain an optimal parameter set; a fitness function used in the optimization process is weighted based on variable weights on the basis of the average relative error of the experimental data curve and the fitting data curve, and the variable weights are determined according to the fitting deviation at the yield point and the plastic section.

[0015] The third aspect of the present application provides an electronic device comprising a processor and a memory, wherein the memory has stored thereon computer instructions which, when executed by the processor, cause the electronic device to perform the method.

[0016] The fourth aspect of the present application provides a computer-readable storage medium having stored thereon a computer program which, when executed by a processor, implements the method.

[0017] The fifth aspect of the present application provides a computer program product comprising a computer program which, when executed by a processor, implements the method.

[0018] The one or more technical solutions above aim at the material properties of manganese aluminum bronze, introduce variable weights based on the fitting deviation at the yield point and the plastic section in the fitness function, and perform iterative optimization based on a fusion optimization algorithm of a genetic algorithm, a grey wolf algorithm and a particle swarm algorithm, which combines the advantages of multiple algorithms, can guarantee the global search range to avoid local optimum, can improve the convergence efficiency, can focus on the key mechanical property regions of manganese aluminum bronze such as the yield point and the plastic section in the optimization process, and improves the fitting precision, so that the optimized constitutive model can more accurately simulate the mechanical behavior of manganese aluminum bronze. BRIEF DESCRIPTION OF DRAWINGS

[0019] The drawings accompanying the specification of the present application form a part thereof, serve to provide further understanding of the present application, and together with the specification explain the exemplary embodiments of the present application and do not limit the present application in any manner.

[0020] Figure 1 A flowchart of the method for identifying the J-C constitutive parameters of manganese aluminum bronze in the embodiments of the present application; Figure 2 A flowchart of the fusion optimization algorithm for solving the constitutive model in the embodiments of the present application; Figure 3 A stress-strain curve of a compression experiment and a preliminary fitting curve thereof; Figure 4 A stress-strain curve at 6500s-1, 600 DEG C and a verification of the preliminary fitting curve thereof; Figure 5 An iteration result graph of the fusion optimization algorithm; Figure 6 Fitting curve result graph for the optimal parameters obtained by the fusion algorithm; Figure 7 Stress-strain curve of compression experiment and its genetic algorithm fitting curve; Figure 8 Comparison chart of stress-strain curve of finite element simulation and stress-strain curve of experiment. DETAILED DESCRIPTION

[0021] Embodiments of the present application will be described in more detail by referring to the attached drawings. Although certain embodiments of the present application are shown in the drawings, it is understood that the present application can be implemented in various forms and should not be interpreted as being limited to the embodiments set forth herein, but rather these embodiments are provided to more thoroughly and completely understand the present application. It should be understood that the drawings and embodiments of the present application are for exemplary purposes only and are not intended to limit the scope of protection of the present application.

[0022] In the description of embodiments of the present application, the term "comprising" and similar terms thereof should be understood as open inclusion, i.e. "including but not limited to". The term "based on" should be understood as "at least partially based on".

[0023] In view of the defects of the prior art in the identification of J-C constitutive parameters of manganese aluminum bronze, the present application obtains data through quasi-static compression experiment combined with multi-temperature and multi-strain rate thermal-mechanical coupling compression experiment, avoids the high cost and long period problem of a large number of grouping experiments in traditional technology, and then adopts the fusion optimization strategy of genetic algorithm, grey wolf algorithm and particle swarm algorithm, balances global exploration and convergence speed by means of global search ability of genetic algorithm to avoid local optimum, directional convergence ability of grey wolf algorithm to compress parameter space, and local refining ability of particle swarm algorithm, breaks through the limitations of single algorithm in the parameter strong coupling scene of manganese aluminum bronze, and introduces variable weight based on yield point and plastic segment fitting deviation in the fitness function, so as to focus on the key mechanical region dominated by the multi-phase state of manganese aluminum bronze, solve the problem of low fitting precision caused by non-target adaptation of core characteristics in the prior art, and finally realize efficient and high-precision identification of J-C parameters of manganese aluminum bronze.

[0024] The J-C constitutive model is as follows:

[0025] σ is the flow strain, ε is the plastic strain, is the strain rate, and T is the actual temperature. The material flow behavior is defined by five parameters A, B, C, n and m, which are the yield strength, strain hardening modulus, strain rate hardening coefficient, strain hardening index and thermal softening index, respectively. In addition, , T r and T mrespectively, the reference strain rate, room temperature and melting temperature, which are taken as 0.001 s -1 , 20℃, 1400℃ for the manganese aluminum bronze alloy.

[0026] For the J-C constitutive model described above, one or more embodiments of the present application provide a J-C constitutive model parameter identification method for manganese aluminum bronze, as shown in Figure 1 and Figure 2 , comprising the following steps: S101, obtaining the true stress-strain curve of the material through the quasi-static compression experiment, and performing the thermal-mechanical coupling compression experiment under different temperature and strain rate conditions to extract the stress-strain response under the corresponding state, and estimate the initial parameter set of the model; S102, performing the thermal-mechanical coupling compression experiment under the set temperature and strain rate conditions, and extracting the stress-strain response under the corresponding state as the verification data; S103, constructing a parameter optimization space based on the initial parameter set, taking the error between the verification data and the fitting data as the fitness function, and iteratively optimizing based on the fusion optimization algorithm of genetic algorithm, grey wolf algorithm and particle swarm algorithm to obtain the optimal parameter set; the fitness function used in the optimization process is: based on the average relative error of the experimental data curve and the fitting data curve, the variable weight is weighted, and the variable weight is determined according to the fitting deviation of the yield point and the plastic section.

[0027] The above technical solution is aimed at the material characteristics of manganese aluminum bronze, introduces a variable weight based on the fitting deviation of the yield point and the plastic section in the fitness function, and iteratively optimizes based on the fusion optimization algorithm of genetic algorithm, grey wolf algorithm and particle swarm algorithm, which combines the advantages of multiple algorithms, can guarantee the global search range to avoid local optimum, and can improve the convergence efficiency, and can focus on the key mechanical property region of manganese aluminum bronze material such as yield point and plastic section in the optimization process, improve the fitting precision, and the obtained constitutive model can more accurately simulate the mechanical behavior of manganese aluminum bronze.

[0028] In addition, the initial parameter set is estimated by combining the quasi-static compression experiment with the thermal-mechanical coupling compression experiment under different temperatures and strain rates, and then the genetic-grey-wolf-particle swarm fusion optimization algorithm (GA-GWO-PSO) is introduced for global optimization, which can obtain constitutive parameters highly consistent with the measured curve under limited experimental conditions, overcome the defects of poor applicability and serious error accumulation of traditional item-by-item fitting method for complex multi-phase materials, and has the advantages of high precision, high efficiency and strong universality, which can be applied to the constitutive parameter inversion and simulation modeling of other high-strength copper-based and multi-phase alloys.

[0029] Step S101 specifically includes: S1011, separate the strain hardening term from the constitutive model, perform compression experiment under quasi-static and room temperature conditions to obtain a first stress-strain curve, and fit the yield strength, strain hardening modulus and strain hardening index based on the strain hardening term. It can be understood that under the condition of extremely low strain rate and room temperature, the strain rate term and the thermal softening term are approximately 1, and the contribution to the total stress can be ignored.

[0030] Specifically, the metal material to be tested is subjected to a Hopkinson pressure bar experiment under quasi-static conditions at a strain rate of 0.001 s -1 , room temperature, recorded as experiment 1, and the true stress-strain curve under this condition is obtained, and the constitutive model formula can be simplified as:

[0031] The yield strength A of the material is 315.98 MPa according to the stress-strain curve of experiment 1, and the strain hardening modulus B and the strain hardening index n are fitted by the least square method according to the simplified formula, recorded as A0, B0, n0; S1012, separate the strain hardening term and the strain rate hardening term from the constitutive model, perform compression experiment under high strain rate and room temperature conditions to obtain a second stress-strain curve, determine the yield strain and yield stress of the material, and calculate the strain rate hardening index. It can be understood that because the temperature is still room temperature, the thermal softening term is still 1. At this time, the difference in stress relative to experiment 1 is mainly caused by the strain rate effect. By comparing the stress values of experiment 2 and experiment 1 at the same strain and different strain rates, the strain rate hardening index C0 can be calculated.

[0032] Specifically, the metal material to be tested is subjected to a Hopkinson pressure bar experiment under quasi-static conditions at a strain rate of 4000 s -1 , temperature of 20℃, recorded as experiment 2, and the true stress-strain curve under this condition is obtained, and the constitutive model formula can be simplified as:

[0033] The yield strain of the material at this time is 0.0591, and the yield strength A is 739.77 MPa according to the stress-strain curve of experiment 2, and the strain rate hardening index C is calculated as 0.0185, recorded as C0. S1013, perform compression experiment under high strain rate and high temperature conditions to obtain a third stress-strain curve, determine the yield strain and yield stress of the material, and calculate the thermal softening index. It can be understood that by comparing the stress values of experiment 3 and experiment 2 at the same strain, similar strain rate and different temperatures, the thermal softening index m0 can be calculated.

[0034] Specifically, the metal material to be tested is subjected to a Hopkinson pressure bar experiment under quasi-static conditions at a strain rate of 4500 s -1The Hopkinson column test was conducted at a temperature of 600℃, designated Experiment 3, to obtain the actual stress-strain curves under these conditions. The constitutive model formula can then be written as:

[0035] Experiment 3 shows that the yield strain of the material is 0.02984 MPa and the yield strength A is 537.87 MPa. The calculated thermal softening index m is 1.9509, denoted as m0. Based on this, an initial parameter set P0=[A0,B0,C0,n0,m0] is obtained. This estimate is based on approximations obtained after decoupling the three effects, and may not be the globally optimal solution. In the following steps, a search interval is constructed based on this initial parameter set, and then the optimal parameter set is found through iterative optimization.

[0036] In step S102, a thermo-mechanical coupling compression experiment is conducted based on temperature and strain rate conditions different from those in the previous experiment, and the stress-strain response under the corresponding conditions is extracted as verification data.

[0037] For example, the metal material to be tested is subjected to a strain rate of 6500 s. -1 A Hopkinson compression test was conducted at 600℃, designated Experiment 4, to obtain the actual stress-strain curves under this condition. Using an independent set of data to evaluate the quality of the parameter set ensures that the identified parameters not only fit known data but also predict the material's behavior under unknown conditions, thereby enhancing the model's generalization and reliability.

[0038] As shown in Table 1, the JC constitutive parameters were initially fitted and solved using the first three sets of specific experimental data. The fitting curves for experiments 1-3 are shown in Table 1. Figure 3 As shown, they correspond to Figure 3 The fitting curves for Experiment 4 are as follows: (a, b, and c are given in the original text.) Figure 4 As shown, the initial constitutive parameters are then extended by ±30%, and the optimal parameters are obtained by fitting the parameters using a fusion optimization algorithm within the new interval. Finally, finite element simulation is performed using Abaqus to verify the results.

[0039] Table 1 Experimental Conditions

[0040] In step S103, the parameter values ​​in the initial parameter set are expanded to obtain the search space. The error between the validation data and the fitted data is used as the fitness function. The optimal parameter set is obtained by iterative optimization based on the fusion optimization algorithm of genetic algorithm (GA), gray wolf algorithm (GWO) and particle swarm optimization (PSO).

[0041] For example, let P0 = [A0, B0, C0, n0, m0]. Expanding P0 by ±30% yields the intervals A = [221.186, 410.774], B = [580.174, 1077.466], C = [0.01295, 0.02405], n = [0.28553, 0.53027], m = [1.36563, 2.53617], and the search space P = [A, B, C, n, m].

[0042] Based on the characteristics of manganese-aluminum bronze—its high strain rate sensitivity, high-temperature softening properties, and yield characteristics sensitive to working conditions—a fitness function is established. Specifically, a weighted average is applied based on the average relative error between the experimental and fitted data curves. The variable weights are determined according to the fitting deviation between the yield point and the plastic segment, with an initial value of 1. The function is defined as follows:

[0043] Where P = [A, B, C, n, m], This indicates the number of sample points in the experimental data. This represents the weight in the t-th iteration, with an initial value of 1. The stress-strain curve for Experiment 4 is shown. Let represent the experimental data for the i-th sample point, and k represent the optimization stage, identified by GA, GWO, and PSO. When k=GA, ... This represents the fitting curve corresponding to the parameter combination during the optimization process within the search space based on the genetic algorithm, when k=GWO. The fitting result of the gray wolf optimization algorithm is given when k=PSO. This is the fitting result of the particle swarm optimization algorithm. The critical transition region and the plastic segment where the yield point is located on the experimental curve both correspond to continuous data sample points, and these data sample points are generally continuous. The data sample points corresponding to the yield point and the plastic segment are mN. This represents the stress values ​​at the yield point and plastic segment of the fitted curve. The weights are then fine-tuned. The value can be 1.4.

[0044] In the fitness function above, variable weights are set on the basis of traditional average relative error to address the high strain rate sensitivity and high temperature softening characteristics of manganese aluminum bronze. The weights are driven by the fitting deviation of the yield point and plastic segment, and the variable weights are amplified in a directional manner to avoid overestimating stress in the high temperature phase transformation region. This better reflects the characteristics of high temperature softening and avoids the model overestimating stress in the high temperature phase transformation region. It also makes the weights automatically tilt towards high strain rate conditions without additional manual intervention.

[0045] The iterative optimization based on the fusion algorithm of genetic algorithm (GA), gray wolf algorithm (GWO) and particle swarm optimization (PSO) specifically includes: (1) Randomly generate an initial population within the search interval, with each individual in the population corresponding to a set of parameters. The total population size is 1000, which serves as the initial population for the GA population.

[0046] (2) Perform a global coarse search on the initial population based on the genetic algorithm (GA) to obtain the GA optimal individuals and generate a set of GA optimal individuals. For example, the set of GA optimal individuals retains the first 500 individuals.

[0047] Specifically, a global coarse search using a genetic algorithm is performed on the GA population. The population size is 1000, and the maximum number of iterations is 200. A tournament selection mechanism is used to select superior individuals for reproduction. Simulated binary crossover (SBX) is used to generate new individuals. Gaussian mutation is used to perturb individual genes with a certain probability to maintain population diversity. During the optimization process, the fitness convergence criterion is set to a relative error of less than 10. -4 To ensure computational accuracy, the best individual in each generation is retained for the next generation, and its output best individual is used as the initial solution for the next stage. At this point, the solution is obtained. Corresponding The optimal individual is then written into both the GWO and PSO populations to achieve information sharing. Leveraging the diversity advantage of the GA population, the solution space can be quickly covered, avoiding local optima traps.

[0048] (3) The best individual in GA is used as the initial α wolf candidate. The best individual in GA is then converged directionally using the gray wolf algorithm within the set of best individuals in GA to obtain the best individual in GWO and generate the set of GWO individuals. For example, the set of best individuals in GWO retains the first 150 individuals.

[0049] Specifically, the Grey Wolf Algorithm is used to achieve directional convergence of the GWO population, with the output of the GA stage as the basis. As the initial alpha wolf candidates for GWO, the new population was subjected to 200 iterations. Every 50 generations, the three individuals with the best fitness were reselected to determine the alpha, beta, and delta wolf pack leaders. The position of ω wolf is updated using the prey-encircling behavior of gray wolves. Let P(t) be the position of a randomly selected gray wolf in the current population. This is to prevent the population from excessively shrinking around individuals α, β, and δ, thus avoiding getting trapped in a local optimum. The solution obtained is... Corresponding The position update formula is:

[0050] in, Used to control the approximation speed and direction. Used to control the amount of disturbance and the range of contraction. , , where a is a constant that decreases linearly from 2 to 0, and r1 and r2 are random numbers between [0,1].

[0051] The Grey Wolf Algorithm, by simulating the cooperative search mechanism of the social hierarchy of grey wolves, enables the population to converge quickly to the region indicated by the GA optimal solution, compressing the optimization range and improving the quality of the solution. Furthermore, by introducing... Random individual perturbations enable the optimization process to guide the population towards the optimal parameter region based on the α, β, and δ leader wolves, while also preventing excessive aggregation of the population due to strong coupling of manganese aluminum bronze parameters, thus maintaining the diversity of population parameters and adapting to the characteristics of high strain rate strengthening and high temperature softening.

[0052] (4) Using the GWO optimal individual as the central initial value of the particle swarm algorithm, and using a portion of the GA optimal individual set and the GWO optimal individual set as the search range, perform local optimization to obtain the global optimal parameter set.

[0053] Specifically, the PSO population is locally refined using the particle swarm optimization algorithm, and the output of the GWO stage is... As the central initial value for PSO, an initial population of PSO is constructed and refined iteratively is performed using the standard formula of the particle swarm optimization algorithm. The adaptive inertia ω decreases linearly from 0.9 to 0.4, the weights of individual c1 and population experience c2 are both 1.6, and r1 and r2 are random numbers between [0,1]. The initial velocity of the particle is 0. This is the optimal solution in the initial population of PSO. This represents the currently calculated particle position in the PSO population. The solution obtained at this point is... Corresponding The standard formula for the particle swarm optimization algorithm is:

[0054] By iteratively updating particle velocity and position, and leveraging the strong local exploitation capabilities of PSO, the system searches within a local region near the GWO optimal solution to achieve precise parameter fine-tuning and obtain the globally optimal parameter set.

[0055] (5) The constitutive parameters are used in finite element simulation to conduct simulation experiments, and the simulated stress-strain curves are obtained. It is then determined whether the accuracy requirements are met, specifically whether the average relative error between the simulated curve and the actual experimental curve is less than 2%. If this requirement is met, then... Let be the constitutive parameter of JC, and if it is not satisfied, let Optimize again.

[0056] The iteration terminated when the average relative error was less than 2% during the third iteration, indicating that the JC constitutive model at this point was the optimal parameter combination model selected by the optimization algorithm. Figure 5 As shown, the objective function value reaches 0.0135, satisfying the model conditions. The fitted curve obtained with these optimal parameters is compared with the experimental curve, as shown below. Figure 6 As shown, it can accurately reflect the mechanical behavior of materials under different strain, strain rate, and temperature conditions. The values ​​of its parameters (such as initial yield stress, strain hardening coefficient, strain rate sensitivity coefficient, temperature softening coefficient, etc.) have achieved a high degree of fit to the experimental data through multiple generations of evolutionary iterations. The final JC constitutive parameters are A=312.77MPa, B=771.71MPa, n=0.38, C=0.0219, and m=1.892.

[0057] like Figure 7 As shown, the stress-strain curves of the plastic stage under four typical experimental conditions (different combinations of strain rate and temperature) were fitted and verified. Experiments 1-4 correspond to... Figure 7 a, b, c, and d in the figure represent the degree of fit (coefficient of determination) before and after optimization. The comparison results are shown in Table 2.

[0058] Table 2 Comparison of Fit Before and After Optimization

[0059] The results show that the fusion optimization algorithm can improve the overall fitting accuracy between the model prediction curve and the experimental curve under different thermo-mechanical loading conditions. Especially under high strain rate and high temperature conditions, the optimized model describes the stress response of the material in the plastic stage more accurately, and the coefficient of determination of the curve is higher. The maximum increase reached 0.0353, and the average increase reached 0.0169. For example... Figure 8 As shown, the obtained JC model parameters are in good agreement with the experimental stress-strain curves in the finite element simulation, and the obtained parameters can predict the mechanical behavior of the material. This optimization strategy significantly improves the stability and versatility of the J-C model parameters under complex thermo-mechanical coupling conditions, and provides a reliable parameter basis for accurate modeling of manganese-aluminum bronze-like multiphase alloys in high-performance cutting, impact, and high-temperature deformation simulations.

[0060] The above one or more embodiments obtain real stress-strain data under different temperatures and strain rates through quasi-static and high-strain-rate compression experiments to establish an initial J-C model parameter set. Then, a Genetic-Grey Wolf-Particle Swarm Optimization (GA-GWO-PSO) fusion optimization algorithm is introduced into the parameter space for global optimization. Combined with dynamic crossover mutation probability adjustment, elite retention, parallel cooperative evolution, adaptive error weighting, and a center-encircling mechanism, efficient identification of five parameters—yield strength, strain hardening modulus, strain rate sensitivity coefficient, strain hardening exponent, and thermal softening exponent—is achieved. Subsequently, through finite element simulation feedback correction and error re-optimization mechanisms, the physical consistency and prediction accuracy of the parameters are further improved. Through multi-intelligent algorithm collaborative optimization and adaptive feedback correction mechanisms, the high-dimensional coupling and convergence instability problems existing in the traditional parameter solution process are solved, significantly enhancing the model's predictive ability under strain rate, temperature, and microstructure effects.

[0061] Furthermore, one or more of the above embodiments can obtain constitutive parameters that are highly consistent with the measured curves under limited experimental conditions, overcoming the shortcomings of traditional term-by-term fitting methods, such as poor applicability to complex multiphase materials and severe error accumulation. They have the advantages of high precision, high efficiency and strong versatility, and can be extended to the field of constitutive modeling and thermo-mechanical coupling simulation calibration of titanium alloys, aluminum alloys and high-temperature structural materials, providing highly reliable theoretical and computational support for the design simulation of complex alloy materials in the fields of shipbuilding, aviation and marine engineering.

[0062] Based on the above method, one or more embodiments of the present invention also provide a device for identifying the constitutive parameters of manganese aluminum bronze JC, comprising: The initial parameter acquisition module is configured to obtain the true stress-strain curve of the material through quasi-static compression experiments, and to conduct thermo-mechanical coupling compression experiments under different temperatures and strain rates to extract the stress-strain response under the corresponding conditions and estimate the initial parameter set of the model. The verification data acquisition module is configured to conduct a thermo-mechanical coupling compression experiment under set temperature and strain rate conditions, and extract the stress-strain response under the corresponding conditions as verification data. The fusion optimization module is configured to construct a parameter optimization space based on the initial parameter set, and perform iterative optimization based on a fusion optimization algorithm of genetic algorithm, gray wolf algorithm and particle swarm algorithm to obtain the optimal parameter set. The fitness function used in the optimization process is: based on the average relative error between the experimental data curve and the fitted data curve, a weighted average is applied, and the variable weight is determined according to the fitting deviation at the yield point and the plastic segment.

[0063] One or more embodiments of the present invention also provide an electronic device that can be used to implement the method for identifying the constitutive parameters of manganese aluminum bronze (JC) described in the above embodiments. The electronic device includes one or more processors, one or more memories coupled to the processors, and a communication module coupled to the processors.

[0064] The memory may include one or more non-volatile memories and one or more volatile memories. Examples of non-volatile memories include, but are not limited to, at least one of the following: read-only memory (ROM), erasable programmable read-only memory (EPROM), flash memory, hard disk, compact disc (CD), digital video disc (DVD), or other magnetic and / or optical storage. Examples of volatile memories include, but are not limited to, at least one of the following: random access memory (RAM), or other volatile memories that do not persist during power-off periods. The computer program may be stored in the ROM. When the processor executes the computer program, it implements the above-described method for identifying the constitutive parameters of JC in manganese aluminum bronze.

[0065] In some embodiments, the program may be tangibly contained in a computer-readable medium, which may include in a device (such as in memory) or other storage device accessible by the device. The program may be loaded from the computer-readable medium into RAM for execution. The computer-readable medium may include any type of tangible non-volatile memory, such as ROM, EPROM, flash memory, hard disk, whereby the computer-readable storage medium stores a computer program that, when executed by a processor, implements the aforementioned method for identifying constitutive parameters of manganese aluminum bronze JC.

[0066] In the above embodiments, implementation can be achieved entirely or partially through software, hardware, firmware, or any combination thereof. When implemented using software, it can be implemented entirely or partially in the form of a computer program product. The computer program product includes one or more computer instructions. When these computer program instructions are loaded and executed on a server or terminal, they generate, in whole or in part, the processes or functions described in the embodiments of this application. The computer instructions can be stored in a computer-readable storage medium or transmitted from one computer-readable storage medium to another. For example, the computer instructions can be transmitted from one website, computer, server, or data center to another website, computer, server, or data center via wired (e.g., coaxial cable, fiber optic cable, digital subscriber line) or wireless (e.g., infrared, wireless, microwave, etc.) means. The computer-readable storage medium can be any available medium accessible to the server or terminal, or a data storage device such as a server or data center that integrates one or more available media. The available medium can be a magnetic medium (e.g., floppy disk, hard disk, and magnetic tape), an optical medium (e.g., digital video disk (DVD), etc.), or a semiconductor medium (e.g., solid-state drive).

[0067] Furthermore, although the operations are described in a specific order, this should be understood as requiring that such operations be performed in the specific order shown or in sequential order, or requiring that all illustrated operations be performed to achieve the desired result. In certain environments, multitasking and parallel processing may be advantageous. Similarly, although several specific implementation details are included in the above discussion, these should not be construed as limiting the scope of this application. Certain features described in the context of individual embodiments may also be implemented in combination in a single implementation. Conversely, various features described in the context of a single implementation may also be implemented individually or in any suitable sub-combination in multiple implementations.

[0068] Although the subject matter has been described using language specific to structural features and / or methodological logic, it should be understood that the subject matter defined in the appended claims is not necessarily limited to the specific features or actions described above. Rather, the specific features and actions described above are merely illustrative examples of implementing the claims.

Claims

1. A method for identifying constitutive parameters of JC in manganese aluminum bronze, characterized in that, Includes the following steps: The true stress-strain curve of the material is obtained by quasi-static compression experiments, and thermo-mechanical coupled compression experiments are carried out under different temperatures and strain rates. The stress-strain response under the corresponding conditions is extracted, and the initial parameter set of the model is estimated. A thermo-mechanical coupling compression experiment was conducted under set temperature and strain rate conditions, and the stress-strain response under the corresponding conditions was extracted as verification data. A parameter optimization space is constructed based on the initial parameter set. An iterative optimization algorithm based on the fusion of genetic algorithm, gray wolf algorithm and particle swarm algorithm is used to obtain the optimal parameter set. The fitness function used in the optimization process is: based on the average relative error between the experimental data curve and the fitted data curve, a weighted average is applied. The variable weight is determined based on the fitting deviation between the yield point and the plastic segment.

2. The method for identifying constitutive parameters of JC in manganese-aluminum bronze as described in claim 1, characterized in that, The process of estimating the initial parameter set of the model includes: The strain hardening term was isolated from the constitutive model, and a compression experiment was conducted under quasi-static, room temperature conditions to obtain the first stress-strain curve. The yield strength, strain hardening modulus and strain hardening exponent were obtained by fitting the strain hardening term. The strain hardening term and strain rate hardening term were separated from the constitutive model. Compression experiments were conducted under high strain rate and room temperature conditions to obtain the second stress-strain curve. The yield strain and yield stress of the material were determined, and the strain rate hardening index was calculated. Compression experiments were conducted under high strain rate and high temperature conditions to obtain the third stress-strain curve, determine the yield strain and yield stress of the material, and calculate the thermal softening index.

3. The method for identifying constitutive parameters of JC bronze as described in claim 1, characterized in that, Iterative optimization based on a fusion optimization algorithm of genetic algorithm, gray wolf algorithm and particle swarm optimization includes: An initial population is randomly generated within the search interval, and each individual in the population corresponds to a set of parameters. A global coarse search is performed on the initial population based on the genetic algorithm to obtain the optimal individual in the GA and generate a set of optimal individuals in the GA. The best individual in GA is used as the initial α wolf candidate. The gray wolf algorithm is used to perform directional convergence within the set of best individuals in GA to obtain the best individual in GWO and generate the set of GWO individuals. Using the best individual in the GWO as the central initial value of the particle swarm optimization algorithm, and taking a portion of the best individual set in the GA and the best individual set in the GWO as the search range, local optimization is performed to obtain the globally optimal parameter set.

4. The method for identifying constitutive parameters of JC bronze as described in claim 3, characterized in that, The fitness function used in the fusion optimization algorithm is: Where P represents the combination of parameters to be optimized. This indicates the number of sample points in the experimental data. This represents the weight in the t-th iteration, with an initial value of 1. Let represent the experimental stress value of the i-th sample point, k denote the optimization stage (identified by GA, GWO, and PSO), and m be the starting point of the continuous sample points corresponding to the yield point and the plastic segment. To fit the stress values ​​at the yield point and in the plastic segment of the curve, This is for fine-tuning the weights.

5. The method for identifying constitutive parameters of JC bronze as described in claim 3, characterized in that, In the Grey Wolf Algorithm, the position update formula is: in, These represent the α, β, and δ wolf pack leaders, respectively. Indicates the position of the gray wolf in the current round. This represents the position of a randomly selected gray wolf from the current group. Used to control the approximation speed and direction. Used to control the amount of disturbance and the range of contraction.

6. The method for identifying constitutive parameters of JC in manganese-aluminum bronze as described in claim 1, characterized in that, The method further includes: Substitute the optimal parameter set into the constitutive model, conduct simulation experiments, and obtain the simulated stress-strain curves; Based on the average relative error between the simulated stress-strain curve and the actual experimental curve, it is determined whether the accuracy requirement is met. If it is met, the optimal parameter set is output. If it is not met, the optimal parameter set is used as the starting individual for the search, and iterative optimization is performed again based on the fusion optimization algorithm of genetic algorithm, gray wolf algorithm and particle swarm algorithm.

7. A device for identifying constitutive parameters of JC in manganese-aluminum bronze, characterized in that, include: The initial parameter acquisition module is configured to obtain the true stress-strain curve of the material through quasi-static compression experiments, and to conduct thermo-mechanical coupling compression experiments under different temperatures and strain rates to extract the stress-strain response under the corresponding conditions and estimate the initial parameter set of the model. The verification data acquisition module is configured to conduct a thermo-mechanical coupling compression experiment under set temperature and strain rate conditions, and extract the stress-strain response under the corresponding conditions as verification data. The fusion optimization module is configured to construct a parameter optimization space based on the initial parameter set, and perform iterative optimization based on a fusion optimization algorithm of genetic algorithm, gray wolf algorithm and particle swarm algorithm to obtain the optimal parameter set. The fitness function used in the optimization process is: based on the average relative error between the experimental data curve and the fitted data curve, a weighted average is applied, and the variable weight is determined according to the fitting deviation at the yield point and the plastic segment.

8. An electronic device, characterized in that, It includes a processor and a memory, the memory storing computer instructions that, when executed by the processor, cause the electronic device to perform the method of any one of claims 1 to 6.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the program is executed by the processor, it implements the method according to any one of claims 1 to 6.

10. A computer program product, the computer program product comprising a computer program, characterized in that, When the computer program is executed by a processor, it implements the method as described in any one of claims 1 to 6.

Images