Magnesium alloy micro-arc oxidation film thermal coupling simulation method and system based on marc

Abstract

The application discloses a magnesium alloy micro-arc oxidation film thermal coupling simulation method and system based on Marc, belongs to the technical field of magnesium alloy thermal coupling simulation, and has the advantages that by constructing a magnesium alloy micro-arc oxidation film finite element model, adapting to the form of a Marc solving framework, and then introducing a pulsed ion beam and customized stress response rules, not only the metallurgical bonding state of the raw material base surface is accurately restored, but also the material nonlinear mechanical behavior under high strain rate and wide temperature range in a real environment is compatible, the mechanical response process of the film layer under strong transient thermal shock working conditions is truly reproduced, a transient heat source function of the energy distribution of a quantitative space is constructed, and the dynamic solving of the coupling thermal stress field is realized, the linkage evolution of the temperature field, the stress field and the material phase change ablation behavior is realized, the reliability of the numerical simulation result is improved, and quantitative numerical support is provided for magnesium alloy micro-arc oxidation film ion beam modification process optimization.

Description

Technical Field

[0001] This invention relates to the field of thermo-mechanical coupling simulation technology for magnesium alloys, specifically to a thermo-mechanical coupling simulation method and system for magnesium alloy micro-arc oxidation films based on Marc. Background Technology

[0002] Existing finite element method (FEM) simulations can reflect the temperature and stress field distributions of magnesium alloy micro-arc oxide films to some extent under high-energy ion beam irradiation. However, current heat source models typically employ idealized Gaussian distributions or uniform heat flow assumptions, failing to adequately consider the transient energy deposition characteristics and spatial non-uniformity of high-intensity pulsed ion beams at the nanosecond scale, particularly the energy superposition effect and local peak variations during scanning, leading to deviations between the heat input process and the actual physical process. Secondly, in terms of simulation methods, current techniques generally employ sequential coupling to handle thermo-mechanical problems, meaning the temperature and stress fields are unidirectionally transferred, neglecting the influence of thermal stress reaction on material thermophysical parameters and local deformation behavior. Therefore, it is difficult to accurately reflect the multi-field coupling effects under strong transient thermal shock conditions. In summary, current numerical simulations struggle to simultaneously and accurately describe the ion beam energy deposition process, transient heat conduction process, and strongly coupled thermo-mechanical response process, thus failing to meet the engineering requirements for predicting the structural evolution and failure behavior of magnesium alloy micro-arc oxide films under high-energy irradiation.

[0003] Chinese Patent Publication No. CN110543744A discloses a method for constructing the thermo-mechanical coupling constitutive relation of zirconium alloys under neutron irradiation. Based on traditional crystal plasticity theory, it considers the influence of irradiation defects and temperature effects on the macroscopic mechanical behavior of zirconium alloys from a physical mechanism perspective. Combining the self-consistent elastoviscoplastic theory, a cross-scale theoretical system from single crystal to polycrystalline is established. Using a self-developed Matlab calculation program, the model parameters can be calibrated simply, quickly, and accurately. The numerical results of the model agree well with experimental data under different external loading conditions. However, it fails to fully consider the transient energy deposition characteristics and spatial non-uniformity of high-intensity pulsed ion beams at the nanosecond scale, resulting in significant discrepancies between the simulation and actual operating conditions, and thus low reliability of the simulation results. Summary of the Invention

[0004] This invention addresses the problem of low reliability in numerical simulation results caused by the inability of existing thermo-mechanical coupling simulation methods for magnesium alloy micro-arc oxidation films to simultaneously and accurately describe the ion beam energy deposition process, transient heat conduction process, and strong thermo-mechanical coupling response process. It provides a Marc-based thermo-mechanical coupling simulation method and system for magnesium alloy micro-arc oxidation films. By constructing a finite element model of the magnesium alloy micro-arc oxidation film and adapting it to the Marc solution framework, and introducing pulsed ion beams and customized stress response rules, it not only accurately reproduces the metallurgical bonding state of the raw material substrate but also accommodates the nonlinear mechanical behavior of materials under high strain rates and wide temperature ranges in real-world environments. It realistically reproduces the mechanical response process of the film layer under strong transient thermal shock conditions. Furthermore, by constructing a transient heat source function that quantifies the spatial energy distribution and dynamically solving the coupled thermal stress field, it achieves the linked evolution of the temperature field, stress field, and material phase transformation ablation behavior, thereby improving the reliability of the numerical simulation results and providing quantitative numerical support for optimizing the ion beam modification process of magnesium alloy micro-arc oxidation films.

[0005] In a first aspect, one technical solution provided in this embodiment of the invention is: a thermo-mechanical coupling simulation method for magnesium alloy micro-arc oxidation films based on Marc, comprising the following steps: S1. Using a magnesium alloy substrate as the bottom layer and a micro-arc oxide film as the surface layer, a finite element model of the magnesium alloy micro-arc oxide film was established in Marc software. S2. A pulsed ion beam is applied to the finite element model of the magnesium alloy micro-arc oxidation film, and stress response rules are set to obtain a thermo-mechanical coupling model. S3. Quantify the energy distribution characteristics of the pulsed ion beam in space in the thermo-coupling model to obtain the transient heat source function, and associate the transient heat source function with the thermo-coupling model to obtain the target simulation model; S4. Start the target simulation model and solve the thermal stress field of the finite element model of the magnesium alloy micro-arc oxide film. Based on the numerical distribution of the thermal stress field changing with time and combined with the computational domain boundary update mechanism, obtain the target simulation numerical of the thermo-mechanical coupling of the magnesium alloy micro-arc oxide film.

[0006] In this scheme, a finite element model of a magnesium alloy micro-arc oxidation film is constructed and adapted to the Marc solver framework to accurately reproduce the metallurgical bonding state of the film-substrate interface, ensuring the continuity of heat conduction and the consistency of mechanical load transfer, and providing a reliable geometric basis for subsequent multiphysics coupling calculations. By introducing pulsed ion beam irradiation load and customized stress response rules, the simplification assumption of traditional homogeneous continuous media is broken, and the nonlinear mechanical behavior of materials under high strain rates and wide temperature ranges is accommodated, thus realistically reproducing the mechanical response process of the film under strong transient thermal shock conditions, improving the physical reality of thermal stress calculations. Furthermore, by constructing a transient... The form of a state-state heat source function and its deep association with the model compensates for the shortcomings of traditional simplified heat source methods in describing the time-domain characteristics of nanosecond-level pulses and the superposition effect of multi-pulse energy. This makes the energy deposition process more closely resemble the actual ion beam irradiation physical process, thus significantly improving the accuracy of subsequent transient temperature field solutions. By coupling the dynamic solution of the thermal stress field with a real-time update mechanism of the computational domain boundary, the linkage evolution of temperature field, stress field and material phase transformation ablation behavior is realized. It can simultaneously predict the peak temperature of the film, stress concentration region and structural failure process, improving the reliability of numerical simulation results and providing quantitative numerical support for the optimization of ion beam modification process of magnesium alloy micro-arc oxidation film.

[0007] Preferably, in S1, a magnesium alloy substrate is used as the bottom layer and a micro-arc oxidation film is used as the surface layer. A finite element model of the magnesium alloy micro-arc oxidation film is established in Marc software, including the following steps: In the Marc software, select the two-dimensional axisymmetric analysis type and set the model geometric parameters based on the actual physical parameters of the magnesium alloy; Based on the model's geometric parameters, mutually fitting upper and lower closed geometric regions are constructed in the same coordinate system, with the geometric center of the upper and lower closed geometric regions serving as the origin of the coordinate system; the upper geometric region is a micro-arc oxide film, and the lower geometric region is a magnesium alloy substrate; A fixed displacement constraint is applied to the lower surface of the magnesium alloy substrate, and the upper surface of the micro-arc oxide film is set as a free thermal boundary. With the origin of the coordinate system as the center, the mesh size is increased sequentially to the edges of the upper and lower closed geometric regions to divide the upper and lower closed geometric regions into meshes and obtain the finite element model of the magnesium alloy micro-arc oxide film.

[0008] In this scheme, by constructing a finite element model using a two-dimensional axisymmetric analysis method, the three-dimensional ion beam irradiation problem can be equivalent to a two-dimensional axisymmetric computational problem. This significantly reduces the element size while ensuring computational accuracy, thereby effectively improving solution efficiency. By constructing mutually fitting double-layer closed geometric regions, the in-situ bonding structure between the micro-arc oxide film and the magnesium alloy substrate can be accurately reproduced, thus ensuring the continuity of interfacial heat conduction and the consistency of mechanical load transfer, which is more in line with the actual metallurgical bonding state. Furthermore, by adopting a gradient mesh generation method that increases from the center to the edge, computational accuracy can be ensured in the strong gradient region of the beam spot effect, and the number of far-field elements can be controlled, thereby balancing solution accuracy and computational efficiency. The corresponding boundary constraints are also highly matched with the actual service conditions, making the simulation process as close to the actual situation as possible.

[0009] Preferably, in S2, a pulsed ion beam is applied to the finite element model of the magnesium alloy micro-arc oxidation film, and stress response rules are set to obtain a thermo-mechanical coupling model, including the following steps: A pulsed ion beam was applied at the center of the free thermal boundary, and the material thermal parameters and basic mechanical parameters were defined for the finite element model of the magnesium alloy micro-arc oxide film. Based on the material's thermal parameters and fundamental mechanical parameters, the stress response induced by the application of a pulsed ion beam to the finite element model of a magnesium alloy micro-arc oxide film is quantified as a stress response function, expressed by the following formula: in, The initial strain rate, For strain attenuation rate, For plastic strain modulus, The yield stress of a magnesium alloy micro-arc oxide film under the influence of a pulsed ion beam is given by the finite element model. The overall stress of the finite element model of the micro-arc oxide film on magnesium alloy. This is the equivalent plastic strain rate; The stress response function is used as the stress response constraint for the finite element model of magnesium alloy micro-arc oxidation film to obtain a thermo-mechanical coupling model.

[0010] In this scheme, by defining the thermal and fundamental mechanical parameters of the material for the finite element model, complete material property support can be provided for subsequent transient heat conduction and thermal stress calculations, thereby adapting to the dynamic changes of material properties over a wide temperature range and ensuring the physical consistency of multi-physics field calculations. By quantifying the stress response function under pulsed ion beam action and embedding it as a constraint into the model, the high strain rate nonlinear mechanical behavior of the material under nanosecond-level strong transient thermal shock can be accurately described, breaking through the accuracy limitations of the traditional linear elastic assumption. Furthermore, the constructed thermo-mechanical coupling model can realize the linkage response of the temperature field and stress field, thereby accurately reflecting the interfacial stress concentration effect caused by the difference in membrane-substrate properties, and effectively improving the accuracy of predicting the mechanical behavior of the membrane under irradiation conditions.

[0011] Preferably, in S3, the energy distribution characteristics of the pulsed ion beam in space in the thermo-coupling model are quantified to obtain the transient heat source function, and the transient heat source function and the thermo-coupling model are correlated to obtain the target simulation model, including the following steps: The coordinates of the heat source center are determined based on the scanning path of the pulsed ion beam and the current time. The total thermal energy deposited at the heat source center coordinates is calculated as the transient heat source function, as expressed in the following formula: in, Let be the energy density of the pulsed ion beam, x and y be the horizontal and vertical distances to the center of the heat source, respectively, and q be the charge of a single ion. , They represent the incident energy as follows: The axial and radial energy distribution of the ion beam, where t is the current time; The target simulation model is obtained by using the transient heat source function as the transient heat source constraint of the thermo-mechanical coupling model.

[0012] In this scheme, the coordinates of the heat source center are dynamically determined based on the pulsed ion beam scanning path and time nodes. At the same time, the total deposition heat energy is quantified to construct a transient heat source function, thereby accurately restoring the non-uniform energy distribution of the ion beam space and the nanosecond-level transient action characteristics. This makes up for the defect that the traditional simplified heat source model has a large deviation from the actual irradiation physical process, thus improving the reliability of the simulation results. Furthermore, the transient heat source function is embedded as a constraint into the thermo-mechanical coupling model to realize the real-time linkage between energy input and material thermo-mechanical response. This ensures that the heat flow loading can be dynamically updated with the irradiation process, effectively improving the solution accuracy of the transient temperature field, and also providing a reliable temperature-driven basis for subsequent stress evolution and phase transformation ablation analysis.

[0013] Preferably, in S4, the target simulation model is started, and the thermal stress field of the finite element model of the magnesium alloy micro-arc oxide film is solved, including the following steps: The target simulation model was started, and the initial room temperature was set. The temperature of each coordinate point in the finite element model of the magnesium alloy micro-arc oxidation film after being affected by the pulsed ion beam was calculated as the temperature field. The formula is as follows: in, This represents the temperature at any coordinate point (x, y) at time t. Here, k is the thermal diffusivity of magnesium oxide, and k is the thermal conductivity of the magnesium alloy. It is the density of magnesium alloy. This refers to the specific heat capacity of magnesium alloys, which is preset according to actual conditions. The thermal stress is calculated based on the temperature at each coordinate point after being affected by the pulsed ion beam, and the formula is expressed as follows: in, Where is thermal stress, and E is the elastic modulus of the magnesium alloy. Indicates the thermal strain modulus. It is the initial room temperature; The thermal stress field is obtained by calculating the thermal stress at each coordinate point in the finite element model of the magnesium alloy micro-arc oxide film.

[0014] In this scheme, the transient temperature field after pulse irradiation is solved point by point with room temperature as the initial condition. This can accurately capture the rapid temperature rise process and strong temperature gradient distribution under nanosecond-level pulse action, thus completely restoring the spatiotemporal evolution law of ion beam energy deposition. Based on the real-time temperature calculation of each coordinate point, the corresponding thermal stress is calculated to achieve spatiotemporal synchronous mapping of temperature field and stress field, ensuring a one-to-one correspondence between thermal load and mechanical response. This solution method can accurately reflect the interfacial stress concentration effect caused by the difference in thermophysical properties of the film substrate, providing quantitative data support for subsequent phase transition determination and structural failure area identification, and effectively improving the engineering reference value of simulation results.

[0015] Preferably, in S4, the target simulation values ​​for the thermo-mechanical coupling of the magnesium alloy micro-arc oxidation film are obtained based on the numerical distribution of the thermal stress field changing with time and combined with the computational domain boundary update mechanism, including the following steps: A fixed time period is used as the simulation window, and the thermal stress field and temperature field are calculated in each simulation window. The computational domain boundary of the finite element model of magnesium alloy micro-arc oxide film is updated according to the computational domain boundary update mechanism; The simulation is repeated for the next simulation window based on the updated computational domain boundary, and the thermo-mechanical coupling values ​​of the thermal stress field and temperature field changing with time are used as the target simulation values.

[0016] This scheme employs a fixed-time-window, step-by-step iterative solution mode, which can adapt to the temporal characteristics of nanosecond-level transient irradiation by ion beams. This allows for precise control of the computational accuracy and temporal resolution at each stage, balancing solution efficiency with the integrity of the evolution process. By introducing a dynamic update mechanism for the computational domain boundary, it can simultaneously track the morphological evolution caused by material melting, vaporization, and ablation, while adjusting the computational range and heat source interaction area in real time. This overcomes the limitation of traditional fixed-domain simulations in being unable to characterize structural evolution. Furthermore, through iterative output of full-time-domain thermo-mechanical coupling values, it achieves dynamic linkage between the temperature field, stress field, and material damage, fully reflecting the evolution law of the membrane structure and providing more realistic quantitative support for membrane failure prediction and modification process optimization.

[0017] Preferably, in S4, the computational domain boundary update mechanism is to analyze the temperature field under the current simulation window. If there is a region with a temperature higher than the boiling point of magnesium alloy, it is considered that vaporization has occurred, and the corresponding region is deleted. The remaining region after deleting the corresponding region is used as the computational domain boundary of the finite element model of magnesium alloy micro-arc oxidation film in the next simulation window.

[0018] In this scheme, the mechanism of determining vaporization ablation based on boiling point temperature threshold and dynamically updating the computational domain boundary can break through the simulation limitations of traditional fixed computational domains and realistically reproduce the dynamic evolution process of vaporization removal of film materials under high-energy ion beam irradiation. This mechanism adjusts the computational range in real time through unit deletion and synchronously links the heat source loading with adaptive correction according to the ablation morphology. It can quantitatively predict the expansion depth and distribution law of the ablation region, making the simulation results more consistent with the actual modification process and providing more accurate numerical support for process parameter optimization.

[0019] Secondly, one technical solution provided in this embodiment of the invention is: a thermo-mechanical coupling simulation system for magnesium alloy micro-arc oxidation film based on Marc, including a physical model construction module, a pulse application module, an energy distribution quantization module, a thermal stress field solution module, a boundary update module, and a result output module; The physical model construction module uses a magnesium alloy substrate as the bottom layer and a micro-arc oxidation film as the surface layer, and uses Marc software to build a finite element model of the magnesium alloy micro-arc oxidation film. The pulse application module applies a pulsed ion beam to the finite element model of the magnesium alloy micro-arc oxide film and sets stress response rules to obtain a thermo-mechanical coupling model. The energy distribution quantification module quantifies the spatial energy distribution characteristics of the pulsed ion beam in the thermo-coupling model to obtain the transient heat source function, and associates the transient heat source function with the thermo-coupling model to obtain the target simulation model; The target simulation model is started, and the thermal stress field solution module solves the thermal stress field of the finite element model of the magnesium alloy micro-arc oxide film. Based on the numerical distribution of the thermal stress field changing with time, the boundary is updated by the computational domain boundary update mechanism through the boundary update module. The result output module outputs the target simulation values ​​of the thermo-coupling of the magnesium alloy micro-arc oxidation film.

[0020] In this solution, a corresponding system is constructed to implement the thermo-mechanical coupling simulation method for magnesium alloy micro-arc oxidation film, thereby realizing human-computer interaction and improving the user experience.

[0021] Thirdly, one technical solution provided in this embodiment of the invention is: a computer device, including a processor, a communication interface, a memory, and a communication bus, wherein the processor, the communication interface, and the memory communicate with each other through the communication bus; the memory is used to store computer programs; and the processor is used to implement the steps of the above-mentioned Marc-based magnesium alloy micro-arc oxidation film thermo-coupling simulation method when executing the program stored in the memory.

[0022] Fourthly, one technical solution provided in this embodiment of the invention is: a computer-readable storage medium storing a computer program, wherein the computer program, when executed by a processor, implements the steps of the above-mentioned Marc-based magnesium alloy micro-arc oxidation film thermo-coupling simulation method.

[0023] The beneficial effects of this invention are as follows: By constructing a finite element model of magnesium alloy micro-arc oxidation film and adapting it to the Marc solution framework, and then introducing pulsed ion beams and customized stress response rules, this invention not only accurately restores the metallurgical bonding state of the raw material substrate, but also accommodates the nonlinear mechanical behavior of materials under high strain rates and wide temperature ranges in real environments. It realistically reproduces the mechanical response process of the film layer under strong transient thermal shock conditions. Furthermore, by constructing a transient heat source function that quantifies the spatial energy distribution and dynamically solving the coupled thermal stress field, it realizes the linkage evolution of temperature field, stress field and material phase transformation ablation behavior, thereby improving the reliability of numerical simulation results and providing quantitative numerical support for the optimization of ion beam modification process of magnesium alloy micro-arc oxidation film.

[0024] The above description of the invention is merely an overview of the technical solution of the present invention. In order to better understand the technical means of the present invention and to implement it in accordance with the contents of the specification, and to make the above and other objects, features and advantages of the present invention more apparent and understandable, specific embodiments of the present invention are described below. Attached Figure Description

[0025] Other features, objects, and advantages of the invention will become more apparent from the following detailed description of non-limiting embodiments with reference to the accompanying drawings. The drawings are for illustrative purposes only and are not intended to limit the invention. Furthermore, the same reference numerals denote the same parts throughout the drawings.

[0026] Figure 1 This is a flowchart of the thermo-mechanical coupling simulation method for magnesium alloy micro-arc oxidation film based on Marc according to the present invention; Figure 2 The temperature field distribution on the surface of the oxide film under irradiation with ion beam energy densities of 1 J / cm² and 2 J / cm² is shown. Figure 3The temperature field distribution on the surface of the oxide film under irradiation with ion beam energy densities of 4 J / cm² and 5 J / cm² is shown. Figure 4 The figure shows the maximum temperature distribution at different depths of the beam spot center of the oxide film under ion beam irradiation with energy densities of 1 J / cm², 2 J / cm², 4 J / cm², and 5 J / cm². Figure 5 The graph shows the temperature variation at different depths (0-4 μm) of the beam spot center of the oxide film under ion beam irradiation with energy densities of 1 J / cm², 2 J / cm², 4 J / cm², and 5 J / cm². Figure 6 The images show the residual stress at different depths on the oxide film surface under irradiation with ion beam energy densities of 1 J / cm², 2 J / cm², 4 J / cm², and 5 J / cm². Figure 7 The total strain diagrams are shown for different depths of the oxide film surface under irradiation with ion beam energy densities of 1 J / cm², 2 J / cm², 4 J / cm², and 5 J / cm². Figure 8 The graph shows the change of thermal stress on the oxide film surface with temperature under irradiation with ion beam energy densities of 1 J / cm², 2 J / cm², 4 J / cm², and 5 J / cm². Figure 9 The graph shows the evolution of plastic strain over time on the oxide film beam spot under ion beam irradiation with energy densities of 1 J / cm², 2 J / cm², 4 J / cm², and 5 J / cm². Figure 10 This is a schematic diagram of the thermo-coupling simulation system for magnesium alloy micro-arc oxidation film based on Marc according to the present invention; Figure 11 A schematic diagram of the structure of the computer device provided by the present invention. Detailed Implementation

[0027] To make the objectives, technical solutions, and advantages of this invention clearer, the invention will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only one preferred embodiment of this invention and are only used to explain this invention. They do not limit the scope of protection of this invention. All other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this invention.

[0028] Before discussing the exemplary embodiments in more detail, it should be mentioned that some exemplary embodiments are described as processes or methods depicted as flowcharts. Although the flowcharts describe the operations (or steps) as sequential processes, many of the operations (or steps) can be performed in parallel, concurrently, or simultaneously. Furthermore, the order of the operations can be rearranged. The process can be terminated when its operation is completed, but it may also have additional steps not included in the figures; the process may correspond to a method, function, procedure, subroutine, subroutine, etc.

[0029] Example 1: To address the problem that existing thermo-mechanical coupling simulation methods for magnesium alloy micro-arc oxidation films cannot simultaneously and accurately describe the ion beam energy deposition process, transient heat conduction process, and strong thermo-mechanical coupling response process, resulting in low reliability of numerical simulation results, this example provides a Marc-based thermo-mechanical coupling simulation method for magnesium alloy micro-arc oxidation films, such as... Figure 1 As shown, it includes the following steps: S1: Using a magnesium alloy substrate as the bottom layer and a micro-arc oxidation film as the surface layer, a finite element model of the magnesium alloy micro-arc oxidation film was established in Marc software.

[0030] In this embodiment, a magnesium alloy substrate is used as the bottom layer and a micro-arc oxidation film is used as the surface layer. A finite element model of the magnesium alloy micro-arc oxidation film is established in Marc software, including the following steps: In the Marc software, select the two-dimensional axisymmetric analysis type and set the model geometric parameters based on the actual physical parameters of the magnesium alloy; Based on the model's geometric parameters, mutually fitting upper and lower closed geometric regions are constructed in the same coordinate system, with the geometric center of the upper and lower closed geometric regions serving as the origin of the coordinate system; the upper geometric region is a micro-arc oxide film, and the lower geometric region is a magnesium alloy substrate; A fixed displacement constraint is applied to the lower surface of the magnesium alloy substrate, and the upper surface of the micro-arc oxide film is set as a free thermal boundary. With the origin of the coordinate system as the center, the mesh size is increased sequentially to the edges of the upper and lower closed geometric regions to divide the upper and lower closed geometric regions into meshes and obtain the finite element model of the magnesium alloy micro-arc oxide film.

[0031] Specifically, a two-dimensional axisymmetric thermo-mechanical coupling computational model was established in the MSC Marc finite element software. The model establishment process was as follows: a magnesium alloy substrate was used as the bottom layer structure, and a micro-arc oxide film layer structure was constructed on its upper surface. The entire structure was set along the Y-axis as the axis of symmetry, thus equating the three-dimensional ion beam irradiation process to a two-dimensional axisymmetric computational problem. The upper surface of the model was defined as the ion beam interaction surface, i.e., the free thermal boundary, where subsequent heat flow would be applied. The bottom of the model was set as a fixed constraint to simulate the support and constraint state of the actual magnesium alloy implant in bone tissue. The left and right boundaries and the axis of symmetry of the model were set as symmetric boundary conditions. Local mesh refinement was performed in the spatial ion beam interaction center region to achieve higher computational accuracy in the areas of maximum temperature gradient and stress concentration.

[0032] Marc software (MSC Marc) is an advanced nonlinear finite element analysis software, mainly used for numerical simulation of complex engineering structures. This embodiment uses this software and constructs a finite element model using a two-dimensional axisymmetric analysis type. This can transform the three-dimensional ion beam irradiation problem into a two-dimensional axisymmetric computational problem, significantly compressing the element size while ensuring computational accuracy, thereby effectively improving solution efficiency. By constructing mutually fitting double-layer closed geometric regions, the in-situ bonding structure between the micro-arc oxide film and the magnesium alloy substrate can be accurately reproduced, thus ensuring the continuity of interface heat conduction and the consistency of mechanical load transfer, which is more in line with the actual metallurgical bonding state. Furthermore, the use of a gradient mesh generation method that increases from the center to the edge can ensure computational accuracy in the strong gradient region of the beam spot and control the number of far-field elements, thereby balancing solution accuracy and computational efficiency. Its corresponding boundary constraints are also highly matched with actual service conditions, making the simulation process as close to the actual situation as possible.

[0033] S2: A pulsed ion beam is applied to the finite element model of the magnesium alloy micro-arc oxidation film, and stress response rules are set to obtain a thermo-mechanical coupling model.

[0034] In this embodiment, a pulsed ion beam is applied to the finite element model of the magnesium alloy micro-arc oxidation film, and a stress response rule is set to obtain a thermo-mechanical coupling model, including the following steps: A pulsed ion beam was applied at the center of the free thermal boundary, and the material thermal parameters and basic mechanical parameters were defined for the finite element model of the magnesium alloy micro-arc oxide film. Based on the material's thermal parameters and fundamental mechanical parameters, the stress response induced by the application of a pulsed ion beam to the finite element model of a magnesium alloy micro-arc oxide film is quantified as a stress response function, expressed by the following formula: in, The initial strain rate, For strain attenuation rate, For plastic strain modulus, The yield stress of a magnesium alloy micro-arc oxide film under the influence of a pulsed ion beam is given by the finite element model. The overall stress of the finite element model of the micro-arc oxide film on magnesium alloy. This is the equivalent plastic strain rate; The stress response function is used as the stress response constraint for the finite element model of magnesium alloy micro-arc oxidation film to obtain a thermo-mechanical coupling model.

[0035] Specifically, after the model is established, the thermal parameters and basic mechanical parameters of the magnesium alloy matrix and the micro-arc oxidation film material are defined, including density ρ, thermal conductivity k, and specific heat capacity C. p Parameters such as the elastic modulus E and the coefficient of thermal expansion α are updated and calculated using the wsklp subroutine to assess the mechanical behavior of the micro-arc oxidation film. By updating the yield stress function related to temperature and plastic strain, constitutive evolution parameters are provided for the finite element stress integration process. At each calculation time step, the system reads the strain rate, plastic strain, and temperature information of the current element and inputs it into the Deshpande-Evans constitutive model for calculation, thereby obtaining the stress response value of the current element. This calculation process yields the stress response function. The wsklp subroutine is a user-defined program based on Marc. Its implementation logic automates the reading of a series of parameters required for the stress response function calculation process, and then continuously calculates the overall stress of the magnesium alloy micro-arc oxidation film finite element model at different times based on the stress response function. This eliminates the need for manual operation, improving the automation level of the simulation process.

[0036] This embodiment defines the thermal and fundamental mechanical parameters of the material for the finite element model, providing complete material property support for subsequent transient heat conduction and thermal stress calculations. This adapts to the dynamic changes in material properties over a wide temperature range, ensuring physical consistency in multiphysics calculations. By quantifying the stress response function under pulsed ion beam action and embedding it as a constraint into the model, it can accurately describe the high strain rate nonlinear mechanical behavior of materials under nanosecond-level strong transient thermal shock, breaking through the accuracy limitations of traditional linear elastic assumptions. Furthermore, the constructed thermo-mechanical coupling model can achieve a linkage response between the temperature field and the stress field, thereby accurately reflecting the interfacial stress concentration effect caused by differences in film-substrate properties and effectively improving the accuracy of predicting the mechanical behavior of the film under irradiation conditions.

[0037] S3: Quantify the spatial energy distribution characteristics of the pulsed ion beam in the thermo-coupling model to obtain the transient heat source function, and associate the transient heat source function with the thermo-coupling model to obtain the target simulation model.

[0038] In this embodiment, the transient heat source function is obtained by quantifying the spatial energy distribution characteristics of the pulsed ion beam in the thermo-coupling model, and the target simulation model is obtained by associating the transient heat source function with the thermo-coupling model. This includes the following steps: The coordinates of the heat source center are determined based on the scanning path of the pulsed ion beam and the current time. The total thermal energy deposited at the heat source center coordinates is calculated as the transient heat source function, as expressed in the following formula: in, Let be the energy density of the pulsed ion beam, x and y be the horizontal and vertical distances to the center of the heat source, respectively, and q be the charge of a single ion. , They represent the incident energy as follows: The axial and radial energy distribution of the ion beam, where t is the current time; The target simulation model is obtained by using the transient heat source function as the transient heat source constraint of the thermo-mechanical coupling model.

[0039] Specifically, at the start of each finite element calculation time step, a custom flux subroutine is defined to first read the current time t and determine which ion beam pulse cycle is currently in; then, the location (x, y) of the heat source center is determined based on the ion beam scanning path. Since the interaction time between ions and target atoms is extremely short compared to the time required for heat transfer, energy is instantaneously deposited on the material surface, and this heat transfer process can be simulated as a two-dimensional heat conduction process. Assuming that the material density remains constant throughout the heating process and any phase transitions, the deposited thermal energy at the heat source center can be calculated using the transient heat source function, where the heat deposition in the x and y directions of the transient heat source function follows a Gaussian distribution: Where Q represents the heat at any location (x, y) and at any time t; Therefore, the total thermal energy deposited at any location can be expressed as: .

[0040] This embodiment dynamically determines the coordinates of the heat source center based on the pulsed ion beam scanning path and time nodes, and constructs a transient heat source function by quantifying the total deposition heat energy. This accurately restores the non-uniform energy distribution and nanosecond-level transient characteristics of the ion beam, making up for the large deviation between the traditional simplified heat source model and the actual irradiation physics process, thereby improving the reliability of the simulation results. Furthermore, the transient heat source function is embedded as a constraint into the thermo-mechanical coupling model to achieve real-time linkage between energy input and material thermo-mechanical response. This ensures that the heat flow loading can be dynamically updated with the irradiation process, effectively improving the accuracy of the transient temperature field solution and providing a reliable temperature-driven basis for subsequent stress evolution and phase transformation ablation analysis.

[0041] In this embodiment, the target simulation model is started, and the thermal stress field of the finite element model of the magnesium alloy micro-arc oxidation film is solved, including the following steps: The target simulation model was started, and the initial room temperature was set. The temperature of each coordinate point in the finite element model of the magnesium alloy micro-arc oxidation film after being affected by the pulsed ion beam was calculated as the temperature field. The formula is as follows: in, This represents the temperature at any coordinate point (x, y) at time t. Here, k is the thermal diffusivity of magnesium oxide, and k is the thermal conductivity of magnesium oxide. It is the density of magnesium oxide. This is the specific heat capacity of magnesium oxide, which is preset according to the actual situation; The thermal stress is calculated based on the temperature at each coordinate point after being affected by the pulsed ion beam, and the formula is expressed as follows: in, Where is thermal stress, and E is the elastic modulus of the magnesium alloy. Indicates the thermal strain modulus. It is the initial room temperature; The thermal stress field is obtained by calculating the thermal stress at each coordinate point in the finite element model of the magnesium alloy micro-arc oxide film.

[0042] This embodiment solves the transient temperature field after pulse irradiation point by point using room temperature as the initial condition. It can accurately capture the rapid temperature rise process and strong temperature gradient distribution under nanosecond-level pulse action, thus completely restoring the spatiotemporal evolution law of ion beam energy deposition. Based on the real-time temperature calculation of each coordinate point, the corresponding thermal stress is calculated to achieve spatiotemporal synchronous mapping of temperature field and stress field, ensuring a one-to-one correspondence between thermal load and mechanical response. This solution method can accurately reflect the interfacial stress concentration effect caused by the difference in thermophysical properties of film substrate, providing quantitative data support for subsequent phase transition determination and structural failure area identification, and effectively improving the engineering reference value of simulation results.

[0043] S4: Start the target simulation model and solve the thermal stress field of the finite element model of the magnesium alloy micro-arc oxide film. Based on the numerical distribution of the thermal stress field over time and combined with the computational domain boundary update mechanism, obtain the target simulation numerical values ​​of the thermo-mechanical coupling of the magnesium alloy micro-arc oxide film.

[0044] In this embodiment, the target simulation values ​​for the thermo-mechanical coupling of the magnesium alloy micro-arc oxidation film are obtained based on the numerical distribution of the thermal stress field over time and combined with the computational domain boundary update mechanism, including the following steps: A fixed time period is used as the simulation window, and the thermal stress field and temperature field are calculated in each simulation window. The computational domain boundary of the finite element model of magnesium alloy micro-arc oxide film is updated according to the computational domain boundary update mechanism; The simulation is repeated for the next simulation window based on the updated computational domain boundary, and the thermo-mechanical coupling values ​​of the thermal stress field and temperature field changing with time are used as the target simulation values.

[0045] This embodiment employs a fixed-time-window, step-by-step iterative solution mode, which can adapt to the temporal characteristics of nanosecond-level transient irradiation by ion beams. This allows for precise control of the computational accuracy and temporal resolution at each stage, balancing solution efficiency with the integrity of the evolution process. By introducing a dynamic update mechanism for the computational domain boundary, it can simultaneously track the morphological evolution caused by material melting, vaporization, and ablation, while adjusting the computational range and heat source interaction area in real time. This overcomes the limitation of traditional fixed-domain simulations in being unable to characterize structural evolution. Furthermore, through the iterative output of full-time-domain thermo-mechanical coupling values, it achieves dynamic linkage between the temperature field, stress field, and material damage, fully reflecting the evolution law of the membrane structure and providing more realistic quantitative support for membrane failure prediction and modification process optimization.

[0046] In this embodiment, the computational domain boundary update mechanism analyzes the temperature field under the current simulation window. If there is a region with a temperature higher than the boiling point of magnesium alloy, it is considered that vaporization has occurred, and the corresponding region is deleted. The remaining region after deleting the corresponding region is used as the computational domain boundary of the finite element model of magnesium alloy micro-arc oxidation film in the next simulation window.

[0047] This embodiment uses a mechanism based on boiling point temperature threshold to determine vaporization ablation and dynamically update the computational domain boundary. This mechanism overcomes the limitations of traditional fixed computational domain simulations and realistically reproduces the dynamic evolution process of vaporization removal of film materials under high-energy ion beam irradiation. The mechanism adjusts the computational range in real time by deleting units and synchronously links the heat source loading with adaptive correction according to the ablation morphology. It can quantitatively predict the expansion depth and distribution pattern of the ablation region, making the simulation results more consistent with the actual modification process and providing more accurate numerical support for process parameter optimization.

[0048] To further illustrate the solution of this embodiment, the solution will be described through the following scenario: Scenario 1: Simulation of mild modification under low energy density pulsed ion beam irradiation: In scenario 1, the ion beam energy density is set to 1-2 J / cm². 2 The pulse width is 50-100 ns, the pulse period is 1-10 μs, the beam radius σ is set to 0-1 mm, and the number of irradiations is 1-10. A transient heat source function is loaded through the flux subroutine, so that the heat flux presents a low-intensity periodic input in time and a wide diffusion distribution in space.

[0049] Under these conditions, the temperature field distribution was obtained by solving the thermal stress field using the finite element method. The results show that the highest surface temperature of the micro-arc oxide film is lower than the melting point T of the material. m=3125K, therefore the computational domain boundary update mechanism was not triggered; only a thermoelastic deformation response occurred, and the computational domain boundary information was recorded. The recorded boundary information indicates that the ablation boundary did not change significantly, and the computational domain maintained a continuous structure.

[0050] Scenario 2: Simulation of densification modification under medium energy density pulsed ion beam irradiation: In scenario 2, the ion beam energy density is set to 3-4 J / cm², the pulse width to 100-150 ns, the pulse period to 800 ns, the beam spot radius σ to 0.3 mm, and the number of irradiations to 1-10. Under these conditions, the flux subroutine performs multi-pulse superimposed heat source loading, which significantly increases the heat flux density in the local area.

[0051] During the solution process of the thermal stress field, the local temperature reaches T. m =3125K, some units enter the melting state, but do not exceed the boiling point T. b =3873 K, thus triggering the computational domain boundary update mechanism, but no large-scale cell deletion occurred, and the computational domain boundary information was recorded. The records show that the melting boundary gradually expanded, but the ablation boundary remained basically stable.

[0052] Scenario 3: Simulation of ablation and failure under high-energy-density pulsed ion beam irradiation: In scenario 3, the ion beam energy density is set to 4-6 J / cm², the pulse width to 150-300 ns, the pulse period to 500 ns, the beam spot radius σ to 0.2 mm, and the number of irradiations to 1-10. Under these conditions, the flux subroutine executes high-intensity multi-pulse superimposed heat flux input, significantly increasing the local heat input rate.

[0053] During the solution of the thermal stress field, the local temperature rapidly exceeds the material's boiling point T. b =3873K, at which point the computational domain boundary update mechanism is triggered. Elements exceeding the threshold are ablated and removed, and the computational domain boundary is updated and recorded in real time. The ablation boundary expansion process is also recorded, showing that the ablation region gradually advances in depth over time, forming a significant material loss area. Simultaneously, due to the extremely high thermal gradient, the flux subroutine dynamically corrects the heat source's area of ​​action based on boundary changes, causing the heat source to shift with changes in surface morphology.

[0054] The thermo-mechanical coupling results of magnesium alloy micro-arc oxidation films in three scenarios with increasing ion beam energy density are as follows: Figure 2-8 As shown, Figure 2 The ion beam energy densities are 1 J / cm² and 2 J / cm², respectively. 2 Temperature field distribution diagram of oxide film surface under irradiation. Figure 3The temperature field distribution diagrams of the oxide film surface under ion beam energy densities of 4 J / cm² and 5 J / cm² are shown. It can be seen that as the ion beam energy density increases, the temperature of the outer layer decreases, indicating that the failure area becomes larger.

[0055] Figure 4 The figure shows the maximum temperature distribution at different depths of the beam spot center of the oxide film under ion beam irradiation with energy densities of 1 J / cm², 2 J / cm², 4 J / cm², and 5 J / cm². As can be seen from the figure, the higher the ion beam energy density, the lower the depth at which the temperature drop occurs.

[0056] Figure 5 The figures show the temperature changes at different depths (0-4 μm) of the magnesium oxide film at the beam spot center under ion beam energy densities of 1 J / cm², 2 J / cm², 4 J / cm², and 5 J / cm², respectively. Figures a, b, c, and d represent ion beam energy densities of 1 J / cm², 2 J / cm², 4 J / cm², and 5 J / cm². It can be seen that the surface temperature of magnesium oxide is not high at ion beam energy densities of 1 J / cm² and 2 J / cm², and no large-area ablation occurs; the main response is thermoelastic deformation. The surface temperature of magnesium oxide is higher at ion beam energy densities of 4 J / cm² and 5 J / cm², and some areas exhibit a melting state.

[0057] Figure 6 The residual stress at different depths on the oxide film surface under ion beam irradiation with energy densities of 1 J / cm², 2 J / cm², 4 J / cm², and 5 J / cm² is shown. Figure 7 The total strain at different depths on the oxide film surface is represented by [value]. Figure 8 This is a graph showing the change of thermal stress on the oxide film surface with temperature. Figure 9 The graph shows the evolution of plastic strain on the oxide film over time. These experimental results clearly and intuitively demonstrate the impact of different ion beam energy densities on the oxide film, providing quantitative numerical support for optimizing the ion beam modification process of magnesium alloy micro-arc oxide films.

[0058] Example 2: This example also provides a thermo-coupling simulation system for magnesium alloy micro-arc oxidation films based on Marc, such as... Figure 10 As shown, it includes a physical model building module, a pulse application module, an energy distribution quantization module, a thermal stress field solving module, a boundary update module, and a result output module; The physical model construction module uses a magnesium alloy substrate as the bottom layer and a micro-arc oxidation film as the surface layer, and uses Marc software to build a finite element model of the magnesium alloy micro-arc oxidation film. The pulse application module applies a pulsed ion beam to the finite element model of the magnesium alloy micro-arc oxide film and sets stress response rules to obtain a thermo-mechanical coupling model. The energy distribution quantification module quantifies the spatial energy distribution characteristics of the pulsed ion beam in the thermo-coupling model to obtain the transient heat source function, and associates the transient heat source function with the thermo-coupling model to obtain the target simulation model; The target simulation model is started, and the thermal stress field solution module solves the thermal stress field of the finite element model of the magnesium alloy micro-arc oxide film. Based on the numerical distribution of the thermal stress field changing with time, the boundary is updated by the computational domain boundary update mechanism through the boundary update module. The result output module outputs the target simulation values ​​of the thermo-coupling of the magnesium alloy micro-arc oxidation film.

[0059] This embodiment implements the thermo-coupling simulation method for magnesium alloy micro-arc oxidation film in this solution by constructing a corresponding system, thereby realizing human-computer interaction and improving the user experience.

[0060] The embodiments also provide a computer device, such as Figure 11 As shown, it includes a processor, a communication interface, a memory, and a communication bus, wherein the processor, the communication interface, and the memory communicate with each other through the communication bus; Memory is used to store computer programs; When the processor executes the program stored in memory, it implements a Marc-based thermal coupling simulation method for magnesium alloy micro-arc oxidation films.

[0061] The communication bus mentioned in the above electronic devices can be a Peripheral Component Interconnect (PCI) bus or an Extended Industry Standard Architecture (EISA) bus, etc. This communication bus can be divided into address bus, data bus, control bus, etc. For ease of illustration, only one thick line is used to represent it in the diagram, but this does not indicate that there is only one bus or one type of bus.

[0062] The communication interface is used for communication between the aforementioned electronic devices and other devices.

[0063] The memory may include random access memory (RAM) or non-volatile memory (NVM), such as at least one disk storage device. Optionally, the memory may also be at least one storage device located remotely from the aforementioned processor.

[0064] The processors mentioned above can be general-purpose processors, including central processing units (CPUs), network processors (NPs), etc.; they can also be digital signal processors (DSPs), application-specific integrated circuits (ASICs), field-programmable gate arrays (FPGAs), or other programmable logic devices, discrete gate or transistor logic devices, or discrete hardware components.

[0065] This application also provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements a Marc-based method for simulating the thermal coupling of magnesium alloy micro-arc oxidation films.

[0066] As can be seen from the above embodiments, it has at least the following substantial effects: (1) By constructing a finite element model of magnesium alloy micro-arc oxidation film and adapting it to the form of Marc solution framework, the present invention accurately restores the metallurgical bonding state of the film-substrate interface, ensuring the continuity of heat conduction and the consistency of mechanical load transmission, and providing a reliable geometric carrier for subsequent multi-physics field coupling calculations. (2) By introducing pulsed ion beam irradiation load and customized stress response rules, this invention breaks through the simplified assumption of traditional uniform continuous medium, and is compatible with the nonlinear mechanical behavior of materials under high strain rate and wide temperature range, thereby truly reproducing the mechanical response process of the film under strong transient thermal shock conditions and improving the physical reality of thermal stress calculation. (3) This invention makes up for the shortcomings of traditional simplified heat source methods in describing the time-domain characteristics of nanosecond pulses and the superposition effect of multi-pulse energy by constructing a transient heat source function to quantify the spatial energy distribution and relating it to the model depth. This makes the energy deposition process more consistent with the actual ion beam irradiation physical process, thereby significantly improving the accuracy of subsequent transient temperature field solution. (4) This invention realizes the linkage evolution of temperature field, stress field and material phase transformation ablation behavior by coupling dynamic solution of thermal stress field and real-time update mechanism of computational domain boundary. It can simultaneously predict the peak temperature of film, stress concentration area and structural failure process, improve the reliability of numerical simulation results, and provide quantitative numerical support for the optimization of ion beam modification process of magnesium alloy micro-arc oxidation film.

[0067] The specific embodiments described above are preferred embodiments of the present invention and are not intended to limit the specific scope of the present invention. The scope of the present invention includes, but is not limited to, these specific embodiments. All equivalent changes made in accordance with the shape and structure of the present invention are within the protection scope of the present invention.

Claims

1. A thermo-mechanical coupling simulation method for magnesium alloy micro-arc oxidation films based on Marc, characterized in that: Includes the following steps: S1. Using a magnesium alloy substrate as the bottom layer and a micro-arc oxide film as the surface layer, a finite element model of the magnesium alloy micro-arc oxide film was established in Marc software. S2. A pulsed ion beam is applied to the finite element model of the magnesium alloy micro-arc oxide film, and stress response rules are set to obtain a thermo-mechanical coupling model. S3. Quantify the energy distribution characteristics of the pulsed ion beam in space in the thermo-coupling model to obtain the transient heat source function, and associate the transient heat source function with the thermo-coupling model to obtain the target simulation model; S4. Start the target simulation model and solve the thermal stress field of the finite element model of the magnesium alloy micro-arc oxide film. Based on the numerical distribution of the thermal stress field changing with time and combined with the computational domain boundary update mechanism, obtain the target simulation numerical of the thermo-mechanical coupling of the magnesium alloy micro-arc oxide film.

2. The thermo-mechanical coupling simulation method for magnesium alloy micro-arc oxidation films based on Marc according to claim 1, characterized in that: In S1, a finite element model of the magnesium alloy micro-arc oxide film is established in Marc software, with the magnesium alloy substrate as the bottom layer and the micro-arc oxide film as the surface layer. The steps include: In the Marc software, select the two-dimensional axisymmetric analysis type and set the model geometric parameters based on the actual physical parameters of the magnesium alloy; Based on the model's geometric parameters, mutually fitting upper and lower closed geometric regions are constructed in the same coordinate system, with the geometric center of the upper and lower closed geometric regions serving as the origin of the coordinate system; the upper geometric region is a micro-arc oxide film, and the lower geometric region is a magnesium alloy substrate; A fixed displacement constraint is applied to the lower surface of the magnesium alloy substrate, and the upper surface of the micro-arc oxide film is set as a free thermal boundary. With the origin of the coordinate system as the center, the mesh size is increased sequentially to the edges of the upper and lower closed geometric regions to divide the upper and lower closed geometric regions into meshes and obtain the finite element model of the magnesium alloy micro-arc oxide film.

3. The thermo-mechanical coupling simulation method for magnesium alloy micro-arc oxidation films based on Marc according to claim 2, characterized in that: In S2, a pulsed ion beam is applied to the finite element model of the magnesium alloy micro-arc oxidation film, and stress response rules are set to obtain a thermo-mechanical coupling model, including the following steps: A pulsed ion beam was applied at the center of the free thermal boundary, and the material thermal parameters and basic mechanical parameters were defined for the finite element model of the magnesium alloy micro-arc oxide film. Based on the material's thermal parameters and fundamental mechanical parameters, the stress response induced by the application of a pulsed ion beam to the finite element model of a magnesium alloy micro-arc oxide film is quantified as a stress response function, expressed by the following formula: in, The initial strain rate, For strain attenuation rate, For plastic strain modulus, The yield stress of a magnesium alloy micro-arc oxide film under the influence of a pulsed ion beam is given by the finite element model. The overall stress of the finite element model of the micro-arc oxide film on magnesium alloy. This is the equivalent plastic strain rate; The stress response function is used as the stress response constraint for the finite element model of magnesium alloy micro-arc oxidation film to obtain a thermo-mechanical coupling model.

4. The thermo-mechanical coupling simulation method for magnesium alloy micro-arc oxidation films based on Marc according to claim 1, characterized in that: In S3, the energy distribution characteristics of the pulsed ion beam in space in the thermo-coupling model are quantified to obtain the transient heat source function. The transient heat source function and the thermo-coupling model are then correlated to obtain the target simulation model, including the following steps: The coordinates of the heat source center are determined based on the scanning path of the pulsed ion beam and the current time. The total thermal energy deposited at the heat source center coordinates is calculated as the transient heat source function, as expressed in the following formula: in, Let be the energy density of the pulsed ion beam, x and y be the horizontal and vertical distances to the center of the heat source, respectively, and q be the charge of a single ion. , They represent the incident energy as follows: The axial and radial energy distribution of the ion beam, where t is the current time; The target simulation model is obtained by using the transient heat source function as the transient heat source constraint of the thermo-mechanical coupling model.

5. The thermo-mechanical coupling simulation method for magnesium alloy micro-arc oxidation films based on Marc according to claim 1, characterized in that: In S4, the target simulation model is started, and the thermal stress field of the finite element model of the magnesium alloy micro-arc oxide film is solved, including the following steps: The target simulation model was started, and the initial room temperature was set. The temperature of each coordinate point in the finite element model of the magnesium alloy micro-arc oxidation film after being affected by the pulsed ion beam was calculated as the temperature field. The formula is as follows: in, This represents the temperature at any coordinate point (x, y) at time t. Here, k is the thermal diffusivity of magnesium oxide, and k is the thermal conductivity of magnesium oxide. It is the density of magnesium oxide. This is the specific heat capacity of magnesium oxide, which is preset according to the actual situation; The thermal stress is calculated based on the temperature at each coordinate point after being affected by the pulsed ion beam, and the formula is expressed as follows: in, Where is thermal stress, and E is the elastic modulus of the magnesium alloy. Indicates the thermal strain modulus. It is the initial room temperature; The thermal stress field is obtained by calculating the thermal stress at each coordinate point in the finite element model of the magnesium alloy micro-arc oxide film.

6. The thermo-coupling simulation method for magnesium alloy micro-arc oxidation films based on Marc according to claim 5, characterized in that: In S4, the target simulation values ​​for the thermo-mechanical coupling of the magnesium alloy micro-arc oxide film are obtained based on the numerical distribution of the thermal stress field over time and combined with the computational domain boundary update mechanism. This includes the following steps: A fixed time period is used as the simulation window, and the thermal stress field and temperature field are calculated in each simulation window. The computational domain boundary of the finite element model of magnesium alloy micro-arc oxide film is updated according to the computational domain boundary update mechanism; The simulation is repeated for the next simulation window based on the updated computational domain boundary, and the thermo-mechanical coupling values ​​of the thermal stress field and temperature field changing with time are used as the target simulation values.

7. The thermo-mechanical coupling simulation method for magnesium alloy micro-arc oxidation films based on Marc according to claim 6, characterized in that: In S4, the computational domain boundary update mechanism analyzes the temperature field under the current simulation window. If there is a region with a temperature higher than the boiling point of magnesium alloy, it is considered that vaporization has occurred, and the corresponding region is deleted. The remaining region after deleting the corresponding region is used as the computational domain boundary of the finite element model of magnesium alloy micro-arc oxidation film in the next simulation window.

8. A Marc-based thermo-coupling simulation system for magnesium alloy micro-arc oxidation films, applicable to the Marc-based thermo-coupling simulation method for magnesium alloy micro-arc oxidation films according to any one of claims 1-7, characterized in that: It includes a physical model building module, a pulse application module, an energy distribution quantization module, a thermal stress field solving module, a boundary update module, and a result output module; The physical model construction module uses a magnesium alloy substrate as the bottom layer and a micro-arc oxidation film as the surface layer, and uses Marc software to build a finite element model of the magnesium alloy micro-arc oxidation film. The pulse application module applies a pulsed ion beam to the finite element model of the magnesium alloy micro-arc oxide film and sets stress response rules to obtain a thermo-mechanical coupling model. The energy distribution quantification module quantifies the spatial energy distribution characteristics of the pulsed ion beam in the thermo-coupling model to obtain the transient heat source function, and associates the transient heat source function with the thermo-coupling model to obtain the target simulation model; The target simulation model is started, and the thermal stress field solution module solves the thermal stress field of the finite element model of the magnesium alloy micro-arc oxide film. Based on the numerical distribution of the thermal stress field changing with time, the boundary is updated by the computational domain boundary update mechanism through the boundary update module. The result output module outputs the target simulation values ​​of the thermo-coupling of the magnesium alloy micro-arc oxidation film.

9. A computer device, characterized in that: It includes a processor, a communication interface, a memory, and a communication bus, wherein the processor, the communication interface, and the memory communicate with each other through the communication bus; the memory is used to store computer programs; the processor is used to execute the program stored in the memory to implement the steps of the Marc-based thermo-coupling simulation method for magnesium alloy micro-arc oxidation film as described in any one of claims 1-7.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program that, when executed by a processor, implements the steps of the Marc-based thermo-coupling simulation method for magnesium alloy micro-arc oxidation films as described in any one of claims 1-7.

Images