112 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 discloses an online detection method for the sealing performance of a sealing element, particularly relates to the technical field of quality detection of industrial sealing elements, and is used for solving the problems of reference value deviation and misjudgment caused by neglecting the influence of mounting stress in an existing detection method. According to the method, a sealed contact interface phase field model is constructed by collecting installation stress data in real time, and crack network permeability is analyzed to predict a leakage diffusion path; high / low stress sensitive areas are divided according to the spatial correlation between the stress gradient and the leakage density, and differential signal acquisition is implemented; extracting a bar code entropy value through topological data analysis, generating a compensation priority sequence, and dynamically correcting a detection reference; and finally, based on the regional weight fusion compensation signal, outputting a high-precision pressure distribution and leakage flow trend result. Cross-field coupling analysis of installation stress and leakage signals is achieved, and the accuracy and reliability of sealing performance detection under complex working conditions are improved.

The invention provides a hydraulic fracture evolution simulation method and device and electronic equipment.The hydraulic fracture evolution simulation method comprises the steps that a damaged rock mass is divided into a plurality of different calculation areas, and a seepage control equation of each calculation area is determined; in combination with seepage anisotropy, determining a seepage tensor of porous medium seepage-stress coupling in each calculation area; associating the seepage control equations of all the calculation areas with phase fields to obtain a target seepage control equation; and determining the displacement, fluid pressure and phase field of the damaged rock mass based on the target seepage control equation. Through the deduced seepage-stress coupling phase field model, the strain energy is decomposed through spectral decomposition, it is guaranteed that crack propagation occurs under the energy generated by tensile strain, and the problem that an existing hydraulic crack research method is limited is solved.

The invention relates to a structure anti-seepage waterproof technology, in particular to a concrete structure anti-seepage performance analysis method considering a multi-physics field coupling effect, which comprises the following steps of: establishing a chemical-thermal-mechanical multi-field coupling model in a concrete hydration process by comprehensively considering various influence factors such as hydration heat, temperature, humidity and structure constraint and combining a phase field theory; solving is carried out through refined numerical simulation, and the concrete crack development and evolution process is obtained. The method has the advantages that multiple influence factors such as hydration heat, temperature, humidity and structural constraint are comprehensively considered by constructing the multi-factor coupled phase field model, the cracking risk of the concrete can be accurately estimated, and comprehensiveness and accuracy of impermeability performance evaluation are improved. By combining finite element simulation and based on a high-precision numerical simulation technology, the anti-permeability performance of the concrete structure under different working conditions can be predicted, and a scientific basis is provided for engineering design.

The invention provides a multi-scale simulation method for high-fidelity selective laser melting forming and phase field-lattice Boltzmann coupling, and belongs to the technical field of selective laser melting numerical simulation. The method comprises the specific steps that a selective laser melting powder laying process model is constructed based on a discrete element method, and a randomly stacked powder bed model is obtained; a high-fidelity selective laser melting forming thermal fluid dynamic model is established, and macroscopic temperature distribution and flow field distribution of a molten pool are simulated; establishing a phase field model-lattice Boltzmann coupling model; and extracting temperature gradient, solidification velocity and melt flow velocity data of local positions of the boundary of the steady-state molten pool, and simulating solidification structure evolution under flowing of molten metal in the molten pool by taking the data as initial conditions of the phase field-lattice Boltzmann coupling model. According to the method, through high-fidelity multi-scale simulation, the relation between the process parameters and the microstructure is established more comprehensively, and reliable numerical support is provided for selective laser melting process optimization and microstructure prediction.

The invention discloses a fracturing fracture extension simulation method based on a self-adaptive multi-scale phase field model, and relates to the field of oil and gas reservoir renovation, and the method is characterized in that reservoir physical property and rock mechanics parameters are obtained, and a reservoir geologic model is established according to the parameters; dividing the reservoir geometric model into a plurality of finite element grids, and assigning grid attributes of the finite element grids by using the parameters to obtain preset hydraulic fractures in the finite element grids and initial variable values on nodes of the finite element grids; according to a preset position relation between the hydraulic fracture and the finite element mesh, constructing a multi-scale adaptive mesh for dynamic expansion of the hydraulic fracture; establishing a coupling mathematical model of a stress field, a seepage field and a phase field of hydraulic fracturing, and constructing a numerical calculation expression of the coupling mathematical model according to the multi-scale adaptive grid; and based on the numerical calculation expression, solving the coupling mathematical model by taking parameters of reservoir physical properties and rock mechanics and boundary conditions of the multi-scale adaptive grid as constraints to obtain hydraulic fracturing fracture extension simulation results under different reservoir and construction conditions.

The invention discloses a fan blade structure damage analysis method and device based on a sub-model, and is used for solving the technical problem that the existing fan blade structure damage and sub-model analysis efficiency is relatively low. The method comprises the following steps: constructing a fan blade finite element model according to obtained fan blade design parameters; performing stress calculation on the fan blade finite element model, and determining a fluid load and a Mises stress value at the Gaussian integral point of each blade unit in the fan blade finite element model at the current time; according to the fluid load, the comparison result between the Mises stress value at the Gaussian integral point of each blade unit at the current time and a preset stress threshold value and the fan blade finite element model, dynamic construction of a target fine mesh entity unit sub-model is completed; and combining a phase field model and a Lagrange-cohesion model to perform structural damage analysis on the target fine grid entity unit sub-model to generate a structural damage analysis result of the target fan blade.

The invention relates to the technical field of metal solidification numerical simulation and calculation material science, and provides a continuous casting billet solidification process phase field simulation method and system based on dynamic temperature boundary conditions. The method comprises the following steps: S1, collecting physical property parameters and component data of steel, continuous casting machine crystallizer parameters and secondary cooling working condition parameters; s2, selecting a phase field method model, and determining initial parameters; s3, solving a phase field control equation, a solute field control equation and a temperature field control equation to obtain phase field variables evolved along with time, solute distribution data and temperature field data; when the temperature field control equation is solved, continuous casting process parameters are introduced into the phase field model, and the actual boundary heat transfer process is simulated by considering the combination of D, N and R thermal boundary conditions. According to the method, any combination and stable solution of the three types of thermal boundary conditions are realized, the adaptability and calculation efficiency of phase field simulation to the complex thermal boundary conditions in the solidification process are remarkably improved, and the analysis precision, the calculation efficiency and the engineering adaptability are considered.

The invention discloses a selective laser melting simulation method simultaneously considering molten pool and microstructure evolution, which mainly comprises establishment and calculation methods of powder laying and a powder melting and solidification process model, and comprises the following specific steps: setting the geometric dimension and the number of internal crystal grains of a substrate to generate a polycrystalline substrate; establishing a discrete element model, simulating a roll shaft powder spreading process, and obtaining particle information; introducing a sequence parameter, constructing a free energy density functional, and establishing a heat-flow-microstructure strong coupling phase field model; constructing a powder bed geometric model, distributing grain orientation, and dividing grids for the powder bed geometric structure; setting material attributes and boundary conditions, and solving the phase field model by a finite element method; and carrying out visualization processing on a calculation result to obtain a molten pool and a microstructure morphology. Evolution behaviors of the molten pool and the microstructure in the selective laser melting process can be simulated at the same time, alternate simulation and data transmission of the powder laying process and the laser scanning process are achieved, and therefore actual machining production is guided, and technological parameters are optimized.

The invention belongs to the technical field of ductile metal fracture analysis, provides an explicit cohesion phase field material point method, and provides a brand new numerical calculation method for metal dynamic fracture damage research. In the method, aiming at the problems of large deformation and high-speed impact in ductile metal spallation and penetration failure, a cohesion phase field material point method and a marble compliant contact algorithm are innovatively provided. The cohesion phase field model can well describe the damage softening and fracture behaviors of the tough metal. The marble compliant contact algorithm effectively suppresses the problems that the speed of impacting the contact surface is over-oscillated and the normal vector outside the contact surface is difficult to calculate. Then, a control equation of dynamic fracture of the ductile metal is deduced through a Lagrange equation, and the control equation of the coupling displacement and cohesion phase field model is discretized in a material point method framework; and finally, a phase field-displacement field interlaced solution strategy is adopted to stably and efficiently solve the problems of contact, large deformation and other strong nonlinear large-scale ductile fracture damage.

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 rock material shear failure and friction fracture simulation method, device and equipment, and the rock material shear failure and friction fracture simulation method comprises the steps: determining a fourth-order tensor matrix corresponding to the stress tensor of a rock material according to a friction fracture phase field model; establishing a crack surface local coordinate system, correcting a traditional degradation function to obtain a target degradation function, and introducing the target degradation function and a crack density function; in the friction fracture phase field model, historical field variables are set, and a friction fracture phase field control equation is determined; and based on the friction fracture phase field control equation, obtaining a stress field, a phase field finite element weak form, a residual error and a stiffness matrix. According to the invention, the traditional phase field model is optimized, the calculation precision of the phase field model is improved, and the simulation result is more reliable.

The invention discloses a blade vibration fatigue strength and three-dimensional fatigue crack propagation prediction method which comprises the following steps: carrying out uniaxial tensile test and fracture toughness measurement on a blade test piece material to obtain material mechanical parameters of the blade test piece material; performing a bending vibration fatigue test on the blade test piece to obtain the fatigue strength of the blade test piece; a three-dimensional fatigue fracture phase field model is constructed, and phase field simulation parameters are obtained through finite element calculation; based on material mechanical parameters and phase field simulation parameters, a blade three-dimensional simulation model is established by simulating a blade sample in a bending vibration fatigue test, a periodic pressure load is applied to a blade basin of the blade to simulate an aerodynamic load in the service process of the blade, the three-dimensional fatigue crack propagation behavior of the blade is calculated, and the vibration fatigue strength of the blade is predicted. According to the method, a three-dimensional fatigue fracture phase field model is established, numerical solution is performed by using a finite element method, the three-dimensional fatigue crack propagation behavior of the blade under the periodic vibration load can be effectively predicted, and the blade fatigue prediction accuracy is improved.

The invention relates to the field of metallurgical casting, in particular to a phase field dynamics method for simulating precipitation evolution of a gamma'phase in nickel-based superalloy under external stress. The method comprises the following steps: S1, acquiring thermal and dynamic information of a gamma phase and a gamma'phase by adopting a CALPHAD method according to target components of the nickel-based superalloy; s2, constructing a strain energy item containing plastic strain, and defining stress-free strain generated by lattice mismatch during gamma-to-gamma 'phase transformation and a gradient energy item capable of considering different slip system information; s3, establishing a phase field model, solving a phase field control equation at a two-phase interface by adopting a Kim-Kim-Suzuki interface model to obtain a sequence parameter result value, and investigating the influence of different forms of external stress loads and different plastic strains on gamma'phase precipitation growth; and S4, the gamma'phase precipitation evolution structure in the nickel-based high-temperature alloy is subjected to visualization treatment, and the influence rule of the external stress and the plastic strain on gamma 'phase precipitation evolution is clarified.

The invention discloses a design method and system for a steel sealing layer of a compressed air energy storage underground artificial reservoir based on an isogeometric phase field method, and the method comprises the steps: building a coupling model of air thermodynamics in the reservoir and heat conduction of a wall medium, and carrying out the simulation calculation of the changes of temperature and pressure in the reservoir in the process of inflation-gas storage-deflation-gas storage; establishing a phase field model for thermal coupling fatigue failure analysis of the steel sealing layer, the concrete lining and the surrounding rock; constructing an isogeometric four-order phase field method calculation model for the fatigue failure analysis of the block type steel sealing layer, and performing numerical simulation based on the influence of a block type layout mode, a welding mode, welding defects and welding residual stress of the steel sealing layer on the fatigue failure of the steel sealing layer; and carrying out statistical analysis on a numerical simulation result, and determining a layout method and a welding mode of the steel sealing layer. According to the method, a new thought is provided for design of the steel sealing layer of the compressed air energy storage underground artificial cave depot, and the method assists in solving the electricity abandoning problem of new energy through the CAES technology.

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 dendritic crystal growth and motion simulation method and system based on an Euler-Lagrangian method.The method comprises the steps that firstly, a phase field model is constructed based on convection flux in a solute field; solving liquid convection on an Euler background grid, and solving solid motion through a Lagrange grid; and updating the moved solid boundary through the conservative phase field model. According to the method, on the basis of simulation of the growth condition of the dendritic crystal, the Lagrange grid is introduced into the discrete interface to simulate the motion behavior of the dendritic crystal, the alloy growth morphology and convection motion are simulated at the same time, the problem that the speed between the grid and the fluid is not matched is solved, and the problem of topological change of a complex interface is solved.

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.