73 results about "Phase field models" patented technology

The invention discloses a fatigue life simulation evaluation method for a lightweight aluminum alloy material of a new energy automobile, and relates to the technical field of material life evaluation. A microstructure image is collected, coupling features are extracted through machine learning, and heterogeneous data fusion and enhancement are completed; generating a topological optimization structure based on a GAN, introducing a VPSC model to describe anisotropy according to a stress gradient dynamic grid, and constructing a dynamic finite element model; fusing vehicle driving data, predicting a load by using LSTM, performing VMD decomposition and environment correction, and realizing space-time correlation load spectrum reconstruction; a phase field model is used in a microcosmic mode, cracks are tracked in a macroscopic mode through XFEM, damage parameters are transmitted in a bidirectional coupling mode, and multi-physics field coupling simulation is carried out; fusing simulation and test data by adopting Bayesian reasoning, calculating life probability distribution, and correcting parameters when errors exceed the limit; according to the method, the fatigue life prediction error is finally reduced, the time consumption of single simulation is reduced, full-life-cycle evaluation and visual early warning are realized, an efficient scheme is provided for lightweight design, and industrial technology upgrading is promoted.

The invention provides a numerical method for immiscible solution drop swallowing behavior in limited shear flow, and relates to the technical field of immiscible solution drop swallowing behavior, a simulation task based on a phase field model and a Navier-Stokes equation is established, a computational domain is discretized through a finite element method, a physical field is initialized, and in each time step, the physical field of the immiscible solution drop swallowing behavior in limited shear flow is calculated. Updated velocity field and phase field data are obtained through coupling solution, the interface curvature change rate and the flow field vorticity intensity are comprehensively calculated, a global index representing the dynamic intensity of the system is generated, the index is compared with a preset threshold value, and the dynamic intensity of the system is calculated. When the index is higher than a threshold value, a smaller step length and a higher grid resolution are adopted, otherwise, a larger step length and a lower resolution are adopted, finally, solver configuration is automatically reconstructed according to the selected parameters, solving is continued under a new adjusting strategy, a self-adaptive numerical simulation cycle is formed, and a new time step is obtained. And completing simulation and outputting complete field data of the liquid drop swallowing process.

The invention discloses a polar diaphragm optimization method based on a high-fidelity electrochemical-phase field model, which comprises the following steps: constructing a high-fidelity electrochemical-phase field model considering a diaphragm substrate and a polar coating to describe the growth behavior of lithium dendrites in the electrodeposition process of a lithium battery; collecting physical property parameters of the diaphragm matrix and the polar coating, substituting the physical property parameters into the model, establishing a computational domain and performing grid division; carrying out numerical solution on the model by utilizing a finite difference method to obtain growth morphology data of the lithium dendrites on the surface of the anode; and performing quantitative characterization on dendritic crystal morphology uniformity based on the growth morphology data to obtain dendritic crystal indexes, adjusting diaphragm physical property parameters according to the indexes, and repeating simulation iteration to optimize the diaphragm structure. By introducing multi-phase field sequence parameters and comprehensively describing the influence of the diaphragm matrix and the polar coating on lithium ion transmission, electrochemical reaction and interface characteristics, the high-precision simulation and optimization of the polar diaphragm design are realized, and the growth of lithium dendrites can be effectively inhibited.

The invention provides a phase field model mass conservation correction method and system for two-phase channel turbulence, and relates to the technical field of two-phase flow numerical calculation, sequence parameters are adopted to describe a current-carrying phase and a dispersed phase in a flow field, the value of a current-carrying phase sequence parameter is 0, and the value of a dispersed phase sequence parameter is 1; comprising the following steps: (1) setting a sequence parameter target total mass M0 according to a mass conservation requirement, the M0 being a volume fraction of a full-field sequence parameter phi i0 at a simulation initial moment, (2) after simulation is started, calculating and updating a sequence parameter phi i (t) of each grid point of the full field at each time step, and (3) determining the total amount delta M (t) of the sequence parameters needing to be corrected according to the difference value between the real-time volume fraction M (t) and the target total mass M0. The method has the advantages that on the basis that an existing solver structure is kept, only the order parameters are operated, the implementation process is simplified, and a mass conservation algorithm which is efficient and easy to implement is achieved.

The invention discloses a reactor nuclear material cross-scale simulation method coupling microstructure simulation and mechanical-thermal performance calculation. The method comprises the following steps: a, planning a most unit sampling point through a most unit method; b, simulating through a mesoscale phase field model to obtain force thermal parameters; c, the force and heat parameters are input into a finite element model to be simulated to obtain force and heat performance; d, the mechanical and thermal performance is input into the phase field model for simulation to obtain the structure morphology of the target nuclear material; e, obtaining a macroscopic scale-mesoscopic scale coupling model by circulating the step c and the step d; f, obtaining a strong coupling model through physical time axis control; and g, constructing a database of the relationship between the mechanical and thermal parameters and the mechanical and thermal performance through a data-driven neural network algorithm. According to the method, strong coupling is achieved through parameter transmission of the mesoscale phase field model and the macroscale finite element, the macroscale and the mesoscale are connected through the most element method, the parameter transmission precision can be adjusted according to convergence, and novel reactor material mechanical and thermal performance prediction can be achieved.

The invention relates to the technical field of numerical simulation of additive manufacturing nickel-based superalloy, and discloses a temperature field and phase field coupling simulation method of SLM forming nickel-based superalloy, which comprises the following steps: establishing a three-dimensional model of a forming part and a substrate for processing to obtain a temperature field basic model; simulating an SLM forming process on the temperature field basic model according to preset parameters to obtain temperature field data, and importing the temperature field data into a phase field simulation module to serve as predefined temperature field data; coupling the gamma ' / gamma mismatch degree and temperature correlation data to the basic phase field model to obtain a gamma' phase change model, correcting the gamma 'phase change model in combination with predefined temperature field data, and adding a nucleation model based on a Poisson seed method; and based on the gamma'phase change model and the nucleation model, growth and coarsening behaviors of the gamma 'phase in the SLM material increasing process are obtained, and a simulation report is formed. According to the method, the vacancy of gamma'phase structure evolution of the existing nickel-based superalloy with high (Al + Ti) content under the SLM complex thermal history condition is solved, and important help is provided for regulating and controlling the structure to obtain the required mechanical property.

The invention discloses a rock material tension-shear failure simulation method and equipment based on a dynamic dual-phase field model. The method comprises the following steps: establishing a stress tensor physical model corresponding to a rock material constitutive relationship; respectively introducing a tension phase field and a shear phase field, and establishing a dynamic double-phase field model corresponding to the rock tension-shear failure physical process; deducing an energy functional of the dynamic double-phase field model; establishing a crack surface local coordinate system, and judging a crack contact state according to a frictional contact phase field theory; switching a stress tensor decomposition mode according to the contact state, and calculating a phase field evolution driving force; rock material tension-shear failure simulation based on the dynamic two-phase field model is achieved, and rock fracture form evolution data are output and used for evaluating rock slope stability.

The invention belongs to the technical field of fracture mechanics calculation, and particularly relates to a hybrid unit beam structure dynamic fracture phase field calculation method. The method comprises the following steps: determining material parameters required by the dynamic fracture phase field calculation of the hybrid unit beam structure; establishing a hybrid unit beam structure dynamic fracture phase field mathematical model according to the Lagrange action quantity principle, wherein the mathematical model comprises a displacement field control equation and a phase field evolution control equation; numerical discretization is carried out on the mathematical model, and a hybrid unit beam structure discrete model is established; and adopting a Newmark method in combination with an interlacing algorithm to carry out numerical solution on the mathematical model after numerical discretization until the displacement field and the phase field both meet convergence conditions, and obtaining a dynamic fracture phase field evolution result of the hybrid unit beam structure. The problem that an existing fracture phase field model cannot accurately simulate crack evolution of the beam structure on the thickness under the dynamic load is solved, and meanwhile the calculation cost for simulating the dynamic fracture of the beam structure is remarkably reduced.

The invention relates to a universal two-phase flow phase field model mesoscopic design method, device and equipment and a storage medium. The method comprises the steps that a two-phase flow control equation based on a phase field model is established; presetting an equilibrium state of the Boltzmann distribution function, calculating a low-order moment, and associating the low-order moment to a macroscopic quantity; performing analysis by using Chapman-Enskog to obtain an integral constraint condition required by a source item; utilizing a Hermite expansion method to obtain an expression of a source item; and outputting a mesoscopic design result. According to the universal two-phase flow phase field model mesoscopic design method provided by the invention, when the phase field model is used for two-phase flow simulation in the future, a mesoscopic format can be quickly obtained even if a new phase field equation is encountered, so that the mesoscopic method can be directly used for simulating the two-phase flow, and mesoscopic design does not need to be carried out again.

The application discloses a preparation method for regulating and controlling ferroelectric texture ceramic grain size by using a phase field method simulation technology, and comprises the following steps: preparing a series of microcrystalline powders with different size distributions; setting an initial phase field structure meeting the particle size distribution; establishing a phase field model for grain orientation growth; simulating and analoging the grain orientation growth process; preparing a film; preparing a green body with microcrystalline parallel distribution; and preparing ferroelectric texture ceramics with target grain size distribution. The application fully utilizes the computer simulation intelligent prediction guidance technology, does not need repeated trial and error experiments, saves time and effort, has good product reproducibility and stability, easily realizes accurate control of the ferroelectric texture ceramic grain size distribution, and is easy to be industrialized; the application effectively controls the occurrence of grain abnormal growth, maintains high grain orientation, and promotes the densification of the ceramic; and the application can not only prepare high-texture-density ferroelectric ceramics, but also regulate and control the grain size.

The application discloses a microstructure evolution prediction method in a selective laser sintering additive manufacturing process, according to the microstructure evolution characteristics of the selective laser sintering additive manufacturing, starting from the free energy density functional of the system, the time and space evolution expressions of temperature, conservative order parameter and non-conservative order parameter are derived in thermodynamic self-consistency, the Cahn-Hilliard equation containing temperature gradient driven diffusion and the determination method of model parameters are given, and numerical calculation is realized through finite element programming. Combined with powder bed laying and finite element calculation, the non-isothermal phase field model can simulate the sintering neck, grain boundary migration and other microstructure evolution in the selective laser sintering additive manufacturing process, simulate the influence of the huge temperature gradient in the selective laser sintering additive manufacturing, optimize the process parameters of the selective laser sintering additive manufacturing, reveal the influence law of the process parameters of the selective laser sintering additive manufacturing on the microstructure evolution and the final microstructure, and improve the design efficiency of the selective laser sintering additive manufacturing.

The invention discloses an integral transformation fracture phase field model implementation method of a random-order polynomial cohesion law, and belongs to the technical field of phase field model finite element analysis. The method comprises the following steps: inputting a parameter equation of a polynomial cohesion law of any order and a target order into a subprogram; deducing a cohesion eigenfunction of a random-order polynomial cohesion law by using a subprogram to determine an analytical expression of a random-order polynomial energy density eigenfunction, solving a target-order polynomial energy density eigenfunction expression and first two-order derivatives thereof, and determining an expansion expression at a zero point by using Maclaurin series; and inputting the target-order polynomial energy density eigenfunction expression, the first two-order derivatives thereof and the expansion expression at the zero point into a phase field model solver to obtain a calculation result of the target-order polynomial integral transformation fracture phase field model. According to the method, the simulation calculation of the integral transformation fracture phase field model of the target-order polynomial cohesion law can be realized under a universal finite element framework, and the accuracy and the efficiency are improved.

The invention discloses a calculation method of an integral transformation fracture phase field model for conformal interpolation reconstruction, and belongs to the technical field of finite element analysis of phase field models. The method comprises the following steps: inputting an interpolation coefficient obtained by performing segmented cubic Hermite interpolation reconstruction on a complex target curve into a subprogram; solving a general term expression of an interpolation segmentation energy density eigenfunction based on a continuity condition by utilizing a subprogram, and determining an interpolation segmentation energy density eigenfunction expression and first two-order derivatives thereof by combining an interpolation coefficient; and inputting the interpolation segmentation energy density eigenfunction expression, the first two-order derivatives of the interpolation segmentation energy density eigenfunction expression and the expansion at the inflection point into a solver for performing an integral transformation fracture phase field model to obtain a simulation result of the integral transformation fracture phase field model of the complex target curve. According to the method, the simulation calculation of the integral transformation fracture phase field model of the complex curve with the inflection point can be realized under a universal finite element framework, and the accuracy and the calculation efficiency are greatly improved.

The invention discloses a nuclear fuel particle fracture analogue simulation method based on a fracture phase field method, belongs to the technical field of nuclear fuel particle fracture simulation, and aims to solve the problems that the research efficiency is low and the application range of a numerical model is limited when a traditional method is used for a complex fracture condition. The method comprises the steps of building a phase field model, inputting material attribute parameters, setting boundary conditions of nuclear fuel particles, adding irreversibility constraints into an evolution equation of the phase field model, setting differentiated analogue simulation step lengths and outputting simulation results. Spatial distribution of phase field variables is used as characterization of a damaged area and a crack form, and real-time tracking of a crack tip position and an expansion path is not needed; for scenes such as multi-particle intersection cracks, branch cracks and irregular cracks which are difficult to process by a traditional method, accurate simulation can be realized only through phase field variable evolution without adjusting a grid or a calculation model.

The application belongs to but is not limited to the technical field of electronic core industry, and discloses a polycrystalline NCM electrode particle cracking and pulverization prediction method and system, which comprises the following steps: constructing an anisotropic polycrystalline fracture phase field model, including a grain boundary phase field equation, a solid mechanics equation, a Li + diffusion equation and a fracture phase field equation; constructing a two-dimensional geometry of a polycrystalline NCM electrode particle, setting required parameters of physical fields, carrying out grid division and setting initial conditions and boundary conditions; using a finite element simulation software to perform numerical solution, and obtaining cracking and pulverization distribution of the polycrystalline NCM electrode particle. + The application considers bidirectional coupling of Li diffusion and mechanical stress, cracking and pulverization and mechanical stress, and uses grain boundary phase field calculation to correct the critical energy release rate, so that the application can accurately reflect the mechanism of cracking and pulverization of the polycrystalline NCM electrode particle, and improve the accuracy of cracking and pulverization prediction.

A simulation method of microvoid evolution in electric field induced interconnection wire based on phase field method, including constructing a phase field model of microvoid evolution in electric field induced face-centered cubic metal {110} orientation "bamboo joint" interconnection wire based on coupled Cahn-Hilliard equation and Allen-Cahn equation; constructing phase field variables of microvoid and grain boundary; constructing distribution functions of electric conductivity and atomic migration rate varying with spatial coordinates; constructing system free energy function F , introducing a free energy function of four double potential well function to modify Cahn-Hilliard equation and Allen-Cahn equation, obtaining phase field control equation of the system; determining conservative phase field order parameter, non-conservative phase field order parameter and electrostatic field parameter, obtaining microvoid evolution in transgranular interconnection wire by coupling phase field and electrostatic field. A direct quantitative tool for the reliability of integrated circuit interconnection wire is proposed.

The invention provides a steel rail fatigue crack simulation method and system based on a fracture phase field fatigue model, and relates to the technical field of rail traffic structure safety and computational mechanics, and the method comprises the following steps: obtaining geometry, material and wheel load parameters of a steel rail; a three-dimensional steel rail finite element model is constructed according to the basic data, and local encryption grids are arranged in wheel-rail contact and other key areas; on the basis, establishing a fracture phase field model capable of describing natural initiation and expansion of cracks; a fatigue degradation model is further integrated, crack resistance attenuation caused by fatigue is simulated, and a phase-field fatigue coupling model is formed; and finally, simulating a moving wheel load through envelope loading, and performing efficient calculation by applying a cyclic transition strategy to realize accurate simulation of a fatigue crack initiation position, an expansion path and damage field distribution of the steel rail. According to the method, crack evolution can be naturally predicted, and the high-cycle fatigue analysis efficiency is remarkably improved.

The application discloses a fracturing fracture extension simulation method based on an adaptive multi-scale phase field model and relates to the field of oil and gas reservoir reconstruction. The technical scheme is as follows: parameters of reservoir physical properties and rock mechanics are acquired, and a reservoir geological model is established according to the parameters; the reservoir geometric model is divided into a plurality of finite element grids, and the grid attributes of the finite element grids are valued by using the parameters, so that preset hydraulic fractures in the finite element grids and initial variable values on nodes of the finite element grids are obtained; a multi-scale adaptive grid for dynamic expansion of the hydraulic fractures is constructed according to the position relationship between the preset hydraulic fractures and the finite element grids; a coupling mathematical model of stress field, seepage field and phase field of hydraulic fracturing is established, a numerical calculation expression of the coupling mathematical model is constructed according to the multi-scale adaptive grid; and the coupling mathematical model is solved based on the numerical calculation expression, with the parameters of reservoir physical properties and rock mechanics and the boundary conditions of the multi-scale adaptive grid as constraints, so that the hydraulic fracturing fracture extension simulation results under different reservoirs and construction conditions are obtained.

The invention belongs to the technical field of microwave rock breaking, and provides a phase field simulation method for quasi-static hydrofracture of a thermoelastic porous medium under the condition of high-energy fracturing, which comprises the following steps of: firstly, acquiring parameters of a target reservoir and parameters of fluid in the target reservoir; secondly, microwave parameters are set on the basis of target reservoir parameters and fluid parameters in a target reservoir, and a phase field model of quasi-static hydraulic fracturing of the thermoelastic porous medium under the microwave-assisted fracturing condition is established; then, the target reservoir is irradiated with the microwave parameters, and the phase field model is used for simulating the crack evolution process during hydraulic fracturing; according to the method, the crack initiation difficulty of hydraulic fracturing can be reduced through microwave-assisted heating.

The invention provides a rock dynamic shear failure simulation method and equipment considering frictional contact. The method comprises the following steps: calculating a stress tensor and a four-order elastic modulus matrix according to an elastic constitutive relation of a rock material; establishing a crack surface local coordinate system based on a maximum shear stress criterion, and determining normal and tangential unit vectors; introducing a crack density function and a corresponding degradation function of the cohesion phase field model, and configuring to eliminate the influence of the phase field characteristic length; deducing an energy functional of a dynamic shear phase field model, judging a contact state according to normal strain of a crack surface, and deducing strong forms of a displacement field control equation and a phase field control equation; variational processing is conducted on the strong forms of the displacement field control equation and the phase field control equation, corresponding weak forms are obtained, rock dynamic shear failure simulation considering frictional contact is achieved based on the weak forms, and rock mass stability evaluation in geological disaster prevention and control is supported. The model is insensitive to the phase field characteristic length, and the simulation result is more reliable.

The application is a phase field model of complex phase structure transformation based on the theory of stacking fault, comprising the following steps: establishing a sublattice free energy model of complex phase structure according to the thermodynamic parameters of the alloy, and thermodynamically describing the alloy system; establishing a stress-strain model of phase transformation according to the stacking fault theory and the phase structure transformation in the alloy, and solving the elastic strain energy of the alloy system; establishing a phase field evolution equation related to the composition and order parameter; setting appropriate initial parameters, solving the phase field equation to obtain the evolution results of the composition field and the order parameter field of the alloy with time and space, and drawing visual images; and analyzing the microstructure evolution diagram of each precipitated phase of the alloy with time. The application provides a phase field model of complex phase structure transformation based on the theory of stacking fault, and the method can predict the process of phase transformation or decomposition of the alloy strengthening phase at high temperature for a long time.

This invention belongs to the field of electrical equipment technology and discloses a method for calculating the dynamic response of dry bushings under seismic loading and establishing a core crack damage model. After simulating the vibration response of transformer bushings under seismic conditions, this invention can obtain various important data, including acceleration response, displacement response, actual stress, deformation degree, and stress concentration. This allows for the adjustment and optimization of various structural and material parameters. Through this simulation, the various physical effects in 110kV resin-impregnated paper dry bushings can be predicted, and the service life of the material in actual scenarios can be accurately determined. This invention establishes two fracture phase-field models to simulate crack propagation in thin plates under different loading conditions. By establishing a phase-field model for isotropic plate fracture problems, numerical calculations are performed using finite element analysis software to study the crack propagation process of thin plates with initial cracks under tensile and shear loads, and the influence of various parameters on crack propagation is calculated.

The present disclosure relates to a two-phase interface area calculation method and system for simulating two-phase flow by a phase field model, the method comprising the following steps: S1, extracting the order parameter value at each lattice point from the phase field model simulation result, and calculating the order parameter gradient value at each lattice point; S2, calculating the bulk free energy density and the interface free energy density of each point in the two-phase system according to the order parameter value and the order parameter gradient value; S3, volume-integrating the calculated bulk free energy density and the interface free energy density to obtain the system free energy F(φ); S4, calculating the surface tension value σ of the two-phase system according to the parameter setting of the phase field model; and S5, calculating the interface surface area A according to the system free energy F(φ) and the surface tension value σ F The method of the present disclosure avoids the complex interface discretization process, thus avoiding the complex grid discretization algorithm and data structure storage mode related thereto, and the method is easy to implement parallel computing and can further improve the calculation efficiency.

This invention discloses a method for establishing a phase-field model of a polycrystalline phase relaxor ferroelectric, the method comprising the following steps: 1. Determining the free energy of the polycrystalline phase relaxor ferroelectric; 2. Determining the order parameter of the polycrystalline phase relaxor ferroelectric as spontaneous polarization. P i 3. Represent free energy as spontaneous polarization P i The invention employs several methods: 1) introducing ferroelectric-ferroelectric phase boundaries and random local electric fields; 2) solving the Ginzburg-Landau partial differential equations to obtain the spatial evolution of spontaneous polarization, thereby determining the domain structure of polycrystalline relaxed ferroelectrics and ultimately establishing a phase-field model for polycrystalline relaxed ferroelectrics. Based on a model ferroelectric system, this invention, by considering ferroelectric-ferroelectric phase boundaries and random local electric fields, reproduces the unique domain structure characteristics of polycrystalline relaxed ferroelectrics with the coexistence of multiple nanometer-sized ferroelectric phases, establishes a phase-field model for polycrystalline relaxed ferroelectrics, and improves the theoretical model of polycrystalline relaxed ferroelectrics, which is beneficial for guiding the design of novel high-performance ferroelectric materials.

The invention provides a microalloy steel multi-mode mechanical property prediction method and device and a readable storage medium, and relates to the technical field of steel performance prediction. The method comprises the following steps: constructing and verifying a phase field model based on thermodynamics and phase field dynamics of a microalloyed steel system to obtain a phase field simulation microstructure image of microalloyed steel; alloy component data and heat treatment process data on a production line are collected, and phase field simulation microscopic structure images and mechanical property indexes corresponding to the alloy component data and the heat treatment process data are obtained; performing standardized preprocessing and sample construction on the phase field simulation microscopic structure image, the alloy component data, the heat treatment process data and the mechanical property indexes, and performing data set division; establishing a multi-modal deep learning model containing an image branch and a numerical branch; carrying out training, evaluation and hyper-parameter optimization based on the training set and the test set, and determining an optimal multi-modal deep learning model; and determining a mechanical property prediction result of the to-be-tested sample based on the optimal multi-modal deep learning model.

The invention relates to the field of modeling methods, and particularly discloses an SAC305 lead-free solder viscoplastic corrosion fracture phase field modeling method, which comprises the following steps: determining SAC305 material viscoplastic constitutive model parameters based on stress-strain data of a uniaxial tensile experiment; carrying out a corrosion experiment on the SAC305 material, and determining a corrosion morphology evolution rule of the SAC305 material; and establishing a corrosion-fracture phase field model, and calibrating phase field model parameters. According to the scheme, constitutive model parameters of SAC305 are calibrated according to a uniaxial tensile experiment, corrosion boundary conditions of the model are confirmed according to a corrosion experiment result, and a corrosion-fracture phase field model is constructed, so that corrosion-induced strength attenuation, crack initiation and expansion paths can be accurately simulated; the method can be used for reliability evaluation, design optimization and service life prediction of SAC305 welding spots in electronic packaging, and provides a basis for engineering application of key interconnection parts in a corrosion environment.

The invention discloses a calculation method for adding an irregular sub-domain without a flux boundary condition in a computational domain, which comprises the following steps of: determining an irregular sub-domain of a substrate / impurity phase added in the computational domain, and drawing the irregular sub-domain into a substrate / impurity phase shape image; processing the image into data with a required model size; establishing a calculation equation of the calculation method coupled with the irregular sub-domain added with the flux-free boundary condition in the calculation domain; performing discrete difference on the calculation equation by adopting a finite difference method to obtain a discrete difference format equation; writing a phase field preposition program by applying a discrete difference format equation; importing the data into a phase field preposition program; taking an operation result of the phase field preposition program as an initial condition of a substrate / impurity phase, and importing the operation result into a constructed phase field model for describing a dendritic crystal growth process to obtain a calculation model coupled with the substrate / impurity phase with an irregular domain; a numerical solution of the calculation model is obtained through computer simulation. According to the invention, the operation efficiency and accuracy are improved.

The application discloses a multiphase flow interaction simulation method based on MPM and a phase field model, and comprises the following steps: step S1, storing mass parameters, phase field parameters and momentum parameters of Lagrange particles in an Euler grid; step S2, updating the phase field of the Euler grid according to the chemical potential to control the energy of the whole system; step S3, updating the momentum of the Euler grid according to gravity, elastic force and the surface tension of the updated equivalent phase field form; step S4, after the grid momentum is updated, processing the boundary on the Euler grid by using a slip boundary condition; step S5, after the phase field and the momentum are updated, converting the momentum parameters and the phase field parameters from the Euler grid variables to Lagrange variables; and step S6, updating the state of the Lagrange particles: updating the position and mass state of the Lagrange particles according to the information returned to the Lagrange particles.

The application discloses a multi-direction forging pass reduction rate distribution optimization method based on a crystal phase field simulation, relates to the technical field of metal material plastic working, and comprises the following steps: constructing a crystal phase field model based on a classical density functional theory and a phase field theory; determining a numerical solution method according to a dynamic equation in the crystal phase field model; simulating a multi-pass forging process based on the crystal phase field model according to forging simulation parameters and a forging reduction rate distribution scheme, and obtaining micro-characteristic data of the forging process through the numerical solution method; evaluating indexes according to the micro-characteristic data, and obtaining a pass reduction rate distribution scheme with the optimal fine-grain strengthening effect as the multi-direction forging optimization process scheme. The above method is based on the crystal phase field model and the numerical solution method, can predict the crystal refining effect of different forging reduction rate distribution schemes, obtain the pass reduction rate distribution scheme with the optimal fine-grain strengthening effect, and does not need to rely on actual trial forging, thereby providing a scientific basis for process optimization.

The application discloses a kind of multi-scale simulation methods of vacuum self-consumption electric arc smelting process ingot, it is related to metal smelting technical field, including: based on the actual working condition parameters and material thermal physical parameters of vacuum self-consumption electric arc smelting, the three-dimensional geometric model of smelting area is established and function area is divided;Multi-physics field coupling numerical model is constructed in macroscopic simulation software, and simulation calculation result is obtained;Simulation result is input to mesoscopic simulation software through standardized data interface, as the initial condition and boundary condition of mesoscopic scale simulation, based on phase field model and solute diffusion equation, the grain growth morphology, solute distribution and microstructure evolution process are simulated;Crystal growth rate, solidification rate and solute distribution result are fed back to multi-physics field coupling numerical model, and iterative calculation is carried out, to realize the closed loop coupling of macroscopic and mesoscopic simulation.The application realizes the accurate control of ingot microstructure and performance.