140 results about "Finite element method" patented technology

The invention relates to the technical field of hydrogen production, in particular to a bipolar plate runner structure-mass transfer efficiency parameter simulation optimization method. Comprising the steps that a flow channel structure is designed, specifically, a bipolar plate flow channel is constructed through a fractal tree-shaped network topological structure, the flow channel is composed of multiple stages of branch channels, the section size of each stage of branch channel is decreased progressively according to a self-similarity proportion, and a three-dimensional flow channel network with fractal dimensions is formed; multi-physics field coupling modeling: establishing a multi-physics field coupling model including fluid flow, electrochemical reaction and heat and mass transfer based on computational fluid mechanics and a finite element method, wherein the model considers a turbulence effect in a flow channel, interface impedance of a porous medium diffusion layer and reaction kinetic parameters at the same time; parameter sensitivity analysis: screening key parameters influencing mass transfer efficiency through orthogonal test design, including flow channel fractal dimension, branch angle, porosity and surface wettability parameters, and constructing a parameter response curved surface; the flow channel structure of the bipolar plate can be accurately optimized, and the mass transfer efficiency is effectively improved.

The invention provides a high earth and rockfill dam digital twinborn model construction method, and belongs to the technical field of high earth and rockfill dam health monitoring, and the method comprises the steps: constructing a measuring point monitoring model based on sensor monitoring data; constructing a mapping model, and constructing a mapping model dynamic monitoring loss function for adversarial training by adopting a domain adaptive method; inputting predicted values of deformation monitoring effect quantity and osmotic pressure monitoring quantity of the measuring point monitoring model on the seepage process and the deformation process of the high earth and rockfill dam into the mapping model to obtain a digital twinborn model, and constructing an extreme working condition digital twinborn model based on the digital twinborn model; knowledge graph reasoning is combined with a dam slope stability coefficient to realize interaction of an extreme working condition digital twinning model and a solid model, and the high earth and rockfill dam deformation-seepage digital twinning framework is proposed for the first time on the basis of deep learning and finite element methods so as to realize a more economical, efficient and reliable monitoring scheme. And technical support is provided for health monitoring of the high earth and rockfill dam structure in the future.

The invention discloses an ice lake feedback monitoring and early warning method based on artificial intelligence and multi-source data fusion, and belongs to the technical field of natural disaster monitoring, and the method specifically comprises the steps: extracting ice lake state feature vectors of multiple spatio-temporal data through cross-modal self-supervised learning, and constructing a three-dimensional virtual ice lake model; based on the feature vectors, constructing a glacier three-dimensional stress field model by using a finite element method, constructing a seepage channel graph model by using a graph neural network, and fusing the features of the glacier three-dimensional stress field model and the seepage channel graph model to form mechanical-seepage coupling state vectors; constructing a pulse neural network to simulate glacier fracture extension and seepage mutation pulse events; constructing a heterogeneous graph federated architecture and outputting a federated weight matrix in combination with a dynamic weighted aggregation strategy; the federal learning weight and the three-dimensional model output result are fused, the outburst posterior probability is calculated through a Bayesian neural network, and a fourth-level early warning signal is generated; according to the invention, the multi-dimensional characterization and coupling process simulation of the state of the ice lake is realized, and the early warning accuracy and timeliness are improved.

The invention discloses an engine turbine shaft sequential optimization design method based on ensemble learning, and the method comprises the steps: carrying out the parametric modeling of an aero-engine turbine shaft through employing a finite element method, and determining a design variable; determining an optimized objective function and constraint conditions based on the turbine shaft structure; constructing a training point set by using the target function and the constraint condition, and establishing an ensemble learning agent model by using the training point set through an ensemble learning algorithm; executing a sequential optimization design process based on the ensemble learning agent model; searching a current optimal sequence design point by utilizing an intelligent algorithm and a sequential updating criterion, and judging whether the sequence design point meets a threshold value requirement or not; if the current optimal sequence design point meets the threshold requirement, outputting the current optimal sequence design point as a design result; otherwise, adding the sequence design points as new training points into the training point set to update the training point set, updating the ensemble learning agent model by using the updated training point set, and iteratively calculating the sequence design points.

The invention provides a generative and interpretable artificial intelligence-based layering design criterion mining method. The method mainly comprises the following steps of: (1) generating a training data set through a finite element method; (2) establishing and training a generative model; (3) establishing and training a symbol regression model; and (4) summarizing a layering design criterion. According to the method, through a data driving thought, a comprehensive composite material layer design rule aiming at a specific design target is mined from big data meeting design requirements.

The invention discloses a fluid distribution dynamic regulation and control method and device based on a shale reservoir microstructure, and relates to the technical field of oil and gas exploitation. The method comprises the steps that microstructure parameters of a target shale reservoir are obtained through a core experiment, and a multi-scale model is constructed based on the microstructure parameters; determining control equations of different physical fields of the multi-scale model; coupling the control equations of different physical fields, and resolving through a finite element method; performing a numerical simulation experiment based on the multi-scale model and the resolving result to determine regulation and control parameters; the regulation and control parameters are input into the multi-scale optimization model to optimize the regulation and control parameters, and a regulation and control strategy is obtained. The problems that an existing displacement experiment is single, the complex interaction process under the actual mining environment is not comprehensively considered, and the guiding effect on mining work is limited are solved. Business personnel can be guided to formulate a development scheme more scientifically and accurately, the mining efficiency and the success rate are improved, the development cost can be reduced, and meanwhile economic benefits are improved.

The invention discloses a parallel computing simulation method for thermal exploitation of a buried pipe well group in a geothermal reservoir, which comprises the following steps of: dispersing a grid volume unit of a buried pipe by adopting a parallel finite element unit method, and simplifying the buried pipe and circulating fluid in the pipe into a one-dimensional linear finite element unit; constructing a plurality of pipeline flow models arranged in an array, and combining the in-pipe one-dimensional linear units and the out-pipe three-dimensional volume units to obtain a well group type pipeline flow geometric model; a parallel finite element unit method is adopted to be combined with a non-isothermal pipeline flow physical field, an underground water seepage physical field and a porous medium heat transfer physical field to conduct multi-physical-field coupling, and a well group type pipeline flow numerical model is obtained; fluid flow and convective heat transfer in the buried pipe and pipe wall heat conduction are solved in a combined mode; and performing factor sensitivity analysis by changing different thermophysical parameters of the rock-soil body and heat exchange characteristic factors of the buried pipe well group. The invention provides a new method for researching fluid flowing and heat transfer behaviors of the buried pipe well group, and provides a new idea for optimal design and efficient simulation of thermal exploitation of the buried pipe well group in a geothermal reservoir.

The invention discloses a spiral acoustic black hole dynamic vibration absorption calculation method based on a reduced multi-body system transfer matrix method. Based on a reduced multi-body system transfer matrix method, a dynamic model of a combined system formed after a main structure is additionally provided with a spiral acoustic black hole vibration absorber is established, and therefore the characteristic frequency and steady-state response of the dynamic model are obtained through solving. Compared with a finite element method, the method has the advantages of high calculation precision and high calculation speed.

The invention relates to the technical field of multi-physics coupling, in particular to a multi-physics field coupling parameterized model establishment method, which comprises the following steps of: establishing a basic equation of each physics field through a physical law and a mathematical principle, and embedding coupling parameters in the equation to describe an interaction relationship among the multiple physics fields to form a coupling equation set; further adopting a finite element method or a finite difference method to solve the model, improving the calculation precision of an interface area through an adaptive grid encryption technology, and combining sensitivity analysis, a gradient descent method or a genetic algorithm to carry out iterative optimization on key parameters; the accuracy of the model is verified through comparison of an experiment and a simulation result, and if a significant difference exists, nonlinear coupling parameters are dynamically adjusted, so that calibration and reliability improvement of the model are realized. Through the reasonable method and steps, an accurate and reliable parameterized model with practical application value is established, so that the complex multi-physics field coupling problem is better understood and solved, and the calculation efficiency is improved.

The invention discloses a digital twinborn model construction and life prediction method for pressure vessel defect data, and belongs to the technical field of pressure vessel safety assessment, and the method comprises the following steps: a defect topological graph construction and feature extraction step: converting defect space distribution into a graph structure, and extracting features by adopting a graph attention network; a crack propagation physical constraint modeling step: predicting a defect propagation behavior by adopting a physical information constraint neural network; an uncertainty quantification and reliability evaluation step: calculating a confidence interval and a failure probability of the residual life by adopting a Monte Carlo Dropout method; the model parameters are dynamically adjusted through an online learning mechanism, high-precision real-time prediction of defect evolution and reliability evaluation of the residual life are achieved, and the calculation efficiency is improved by more than two orders of magnitude compared with a traditional finite element method.

The invention relates to a finite element method for calculating temperature distribution of a steel strip in a continuous hot galvanizing process, and belongs to the technical field of hot galvanizing methods. According to the technical scheme, the method comprises the steps that relevant sizes and technological parameters of a continuous hot galvanizing production line are obtained; model establishment and grid division are conducted, a three-dimensional model is established according to the related size, the model comprises molten zinc and a steel belt, and the molten zinc and the steel belt are set to be solid models; and setting boundary conditions, performing heat exchange calculation, and processing a calculation result. The method has the beneficial effects that the method is used for analyzing the temperature distribution of the steel strip in the hot galvanizing process and analyzing the flow field and the temperature field of molten zinc; and the change of working condition parameters such as the temperature of the steel brought into the zinc pot, the moving speed of the steel strip, the thickness and width of the steel strip, the temperature and depth of molten zinc and the like is accurately responded, so that the optimization of the continuous hot galvanizing process is guided.

The invention provides a concrete deformation prediction method and device based on a multi-field coupling model, and relates to the technical field of concrete structure engineering and multi-physics field coupling computation.The method comprises the steps that main variables of a hydration reaction field, a temperature field and a humidity field are discretized through a finite element method, and the discretized main variables are obtained; establishing a control equation and carrying out time stepping solution; calculating the evolution of a hydration reaction field by using a Newton-Raphson iteration method, and solving temperature and humidity distribution based on boundary conditions and a numerical integration method; by setting convergence error tolerance, whether calculation reaches balance or not is judged, and finally temperature deformation and humidity deformation distribution in the concrete is obtained. According to the method, accurate deformation prediction of the concrete structure in the hydration reaction process is achieved, error accumulation caused by single-field calculation in a traditional method is avoided, and calculation precision and engineering applicability are improved.

The present invention discloses a multi-objective thermal-fluid topological design optimization method for the cooling element of a precision gear grinding machine tool, comprising the following steps: Step 1: Define the design domain: Expand the cooling element into a plane, define a two-dimensional design domain representing the cooling element's geometric structure, and discretize it using the finite element method; Step 2: Numerical modeling: Model the cooling element as a porous medium; Under the conditions of laminar incompressible flow, the fluid dynamics are governed by the dimensionless forms of the continuity equation and the momentum equation; The fluid-solid coupled heat transfer model is governed by the dimensionless energy conservation equation; Step 3: Construct an objective function: The optimization objective is defined as a weighted function of heat transfer and fluid dissipated power, and a fluid-solid heat transfer topological optimization model is constructed; Step 4: Solve the objective function; Step 5: Restore the shape and structure of the cooling element. The present invention also discloses a precision gear grinding machine tool and a cooling element.

The invention provides a mining electric shovel working device health detection method based on multidisciplinary digital twinning, and the method comprises the steps: obtaining attitude data through employing a Latin hypercube sampling method; performing electric shovel simulation by using discrete element and finite element methods, and further obtaining a data set required for neural network training; training and testing the LSTM neural network and the RBF neural network by using the data set, and constructing a mining load prediction model and a mining stress prediction model; constructing a digital twinborn model by using a tilt angle sensor, a stay wire encoder and the virtual model; a laser radar is used for scanning the terrain, the material mass is calculated and obtained, the load and stress of data are sequentially predicted through a neural network model, the stress is transmitted to a digital twinborn model, and real-time monitoring is achieved in a color rendering mode. According to the mining electric shovel, the stress change of the working device is displayed in real time in the working process of the mining electric shovel, and the early warning capacity of the mining electric shovel when excavation exceeds the stress safety range in the complex working environment is improved.

The invention provides a method and a system for predicting the frictional wear residual size of a transmission line connection fitting, and the method comprises the steps: calculating the maximum wind speed and swing angle of the connection fitting through the data of a meteorological database and an engineering design data database, and simulating the wear process of the transmission line fitting through a finite element method. The method comprises the steps of obtaining a change relation of the abrasion loss of the connection hardware fitting along with the abrasion frequency, determining the average swing frequency of the hardware fitting according to the current operation and maintenance result of the connection hardware fitting, and predicting the remaining size of the hardware fitting after n years on the basis of the assumption that the connection hardware fitting operates according to the same swing frequency. The method and the system can scientifically and reasonably evaluate the service life of the power transmission line connection hardware fitting, and avoid careless mistakes caused by the fact that the worn hardware fitting is found only through sense organs of people in a traditional manual inspection mode.

The invention relates to the technical field of urban flood control and drainage, and discloses an urban inland inundation rainfall flood data assimilation method and system and a storage medium. The method comprises the steps of obtaining rainfall, flow and surface elevation data from multi-source hydrological observation data, and obtaining a multi-source data set after preprocessing; hydrological parameters are extracted based on the data set, and the influence weight of the hydrological parameters on the submerging range is analyzed and determined in combination with parameter sensitivity; updating the key hydrological parameters through a filtering algorithm to obtain a refined hydrological parameter set; a terrain gradient analysis and prediction algorithm is combined to establish a prediction model of a submerging boundary change rate, and an inversion time interval is dynamically adjusted; the surface elevation is corrected in real time through a finite element method, and the submerging prediction precision is improved; and further optimizing the hydrological parameters according to the real-time prediction result, and generating a drainage scheduling decision instruction. The method has high accuracy and real-time performance, scientific and effective decision support can be provided for an urban flood control and drainage system, and the emergency response capability is remarkably improved.

A method for simulating heat storage of a buried pipe well group in a geothermal reservoir adopts staggered arrangement to replace traditional sequential arrangement, so that the distance between buried pipes is more reasonable, heat interference is minimized, the volume shape of the buried pipes is simplified into a linear unit structure, pipeline volume units are dispersed by adopting a finite element method based on various theories, and the heat storage efficiency of the buried pipe well group is improved. The calculated amount is reduced; and fluid heat transfer characteristics in the pipeline are accurately captured. Then setting soil initial temperature field distribution and initial parameters of fluid in a buried pipeline unit, solving a pipeline unit model by adopting a proper numerical calculation method, and performing real-time monitoring to ensure that a calculation result is stable and accurate; factor sensitivity analysis is carried out by changing in-well parameters and operation parameters, the influence of rock-soil body thermophysical parameters and buried pipe heat exchange characteristic factors on the geothermal reservoir thermal exploitation efficiency is comprehensively researched, and a scientific basis is provided for optimal design and operation of a buried pipe well group. The invention further comprises a system, equipment and a storage medium for implementing the method.

The invention provides a freezing reinforcement project evaluation system and method based on hydrothermal coupling numerical model simulation, and relates to the technical field of computers and frozen soil engineering.According to the freezing reinforcement project evaluation system and method based on hydrothermal coupling numerical model simulation, a water flow equation and a heat transfer equation are established, and corresponding boundary and initial conditions are set; dividing a frozen soil reinforcement area into finite elements, discretizing an equation by utilizing a finite element method, iterating a global equation through a numerical solution method, simulating moisture and temperature changes in a freezing process, extracting key indexes, generating a freezing reinforcement energy efficiency index, comparing the freezing reinforcement energy efficiency index with a set efficiency interval, and evaluating a reinforcement engineering effect. And corresponding optimization suggestions are provided. According to the method, effective evaluation and optimization of the freezing reinforcement project are realized by establishing the hydrothermal coupling numerical model.

The invention relates to the technical field of hydrographic survey, in particular to a riverway water flow calculation method, which comprises the following steps: acquiring topographic data, hydrological data and roughness data, and constructing a one-dimensional riverway hydrodynamic numerical model; the one-dimensional river channel hydrodynamic numerical model is based on a discontinuous finite element method and a Runge-Kutta time discrete shallow water equation set, and a numerical flux is solved by combining an essential oscillation-free ENO slope limiter and an HLL Riemann solver; aiming at the confluence situation of tributaries in different intervals, the interval flow is calculated by utilizing an iterative method, including initially assuming that the tributary flow is zero, simulating non-tributary outflow by using a one-dimensional river hydrodynamic numerical model, determining the difference between the simulated non-tributary outflow and actual site flow measurement, and adjusting the tributary flow; and iteration is repeated until the error between the simulated outlet flow and the actual outlet flow reaches the preset precision. By constructing a one-dimensional river water power numerical model based on a Runge-Kutta time discrete shallow water equation set, high-precision simulation of river water flow calculation under the condition of branch flow convergence in different intervals is realized.

The invention relates to a method for calculating the power-on and power-off delay time performance of an electromagnetic valve, and the method comprises the following steps: obtaining the mass of a valve element of the electromagnetic valve, the spring stiffness and coil resistance data; establishing an electromagnetic valve equivalent circuit model, and listing a current equation; establishing an electromagnetic valve magnetic circuit mathematical model, and listing a magnetic field equation set; establishing an electromagnetic valve flow field mathematical model, and listing a flow field equation set; establishing a kinematic model of the electromagnetic valve, and listing a kinematic equation; using a finite element method to obtain a valve core electromagnetic force-displacement-current data table and a coil inductance-displacement-current data table; a finite element method is used for obtaining a valve element aerodynamic force-displacement data table; and obtaining the power-on and power-off delay time through data analysis processing and solving of the partial differential equation set. According to the method, prediction of the opening and closing delay response time of the electromagnetic valve is more accurate, visual adjustment in the simulation process is achieved, and the time for obtaining the response time of the electromagnetic valve and the test cost are greatly reduced while the preset precision requirement is met.

The invention discloses a dynamic behavior analysis method and system for kiln high-temperature synthesis, and relates to the technical field of high-temperature material synthesis monitoring, and the method comprises the steps: obtaining time sequence data through a high-frequency temperature sensor, a reaction intermediate detector and a flue gas component monitor, and carrying out the feature extraction; the micro-scale adopts molecular dynamics simulation to obtain local reaction dynamic behaviors, and the macro-scale simulates kiln temperature distribution based on a finite element method; fusing the feedback relationship between the local reaction rate and the temperature distribution by using multi-scale coupling; dividing the high-temperature synthesis process of the kiln into a plurality of continuous recognizable states, judging whether a state jump condition is met or not, and performing multi-objective optimization iteration to obtain a dynamic behavior analysis result and phase change stage judgment. According to the method, dynamic feedback of the reaction rate and the temperature field is realized, key inflection points and phase change characteristics can be accurately captured, the fitting precision of kinetic parameters is remarkably improved, and the robustness and the adaptive capacity of the model are enhanced.

The invention discloses a flow-induced noise integration method, device and equipment based on a finite element method and a storage medium, and relates to the technical field of noise integration, and the method comprises the steps: obtaining an aerodynamic noise source based on a sound source type and a time domain flow field file; transmitting an aerodynamic noise source from the fluid grid to the acoustic grid based on the acoustic grid file and a data mapping algorithm to obtain a time domain sound source; converting the time domain sound source into a frequency domain sound source based on a time-frequency conversion strategy; and constructing a numerical model based on the frequency domain sound source, and solving the numerical model to obtain a sound propagation result. The problems that in a traditional mixing method, due to the fact that a plurality of independent modules exchange large-scale intermediate data through a file system, the process is split, reading and writing are time-consuming, and user operation is tedious are solved, integrated solving of aerodynamic noise simulation from sound source extraction to sound field calculation is achieved, user operation steps are simplified, and user experience is improved. Therefore, on the premise that the calculation precision is guaranteed, the overall execution efficiency and the automation level of the simulation process are remarkably improved.

The invention discloses a multi-physics field coupling numerical calculation method based on multiplexing shift Laplace preprocessing acceleration. In a multi-physical field coupling calculation process induced by a frequency domain electromagnetic field, a group of electromagnetic linear systems discretized by a finite element method which changes slowly needs to be solved numerically, and generally, time and memory are consumed very much. The method can effectively reduce the number of times of inversion operation of the preprocessing matrix, and significantly reduces the demand for computing resources. Different from a traditional algorithm using approximate inverse of a preprocessing matrix, the method provided by the invention constructs accurate inverse of the preprocessing matrix, and improves the stability and convergence rate of an iteration method, thereby ensuring high efficiency and high robustness of linear system solving. According to the method, the shift Laplace preprocessing technology is introduced into the multi-physics field coupling calculation for the first time, and the electromagnetic linear system sequence in the multi-physics field coupling simulation can be efficiently calculated.

The application relates to a seepage simulation method in a multi-scale fractured medium based on coupling finite elements, which comprises the following steps: 1. acquiring statistical distribution characteristics of related parameters of multi-scale fractures in a simulation domain, and constructing a multi-scale fracture geometric model; 2. respectively performing intersection line detection on each scale fracture in the multi-scale fracture geometric model; 3. simplifying the multi-scale fracture geometric model; 4. adopting coupling units to establish coupling relationships among the each scale fracture models; and 5. simulating a multi-scale discrete fracture-matrix model. The application is based on a coupling finite element method, and couples three scale fracture conceptual models, so that the difficulties in calculation precision, calculation cost and mesh partitioning of a single traditional model in solving multi-scale problems are avoided, and a feasible new way is provided for large-scale complex rock mass underground water flow analysis.

The present invention relates to the field of image processing technology, and more particularly to a method and system for visual inspection of small glass products based on the finite element method. The method comprises collecting images of glass products with defects; constructing a DLHPDE model and calculating weight values ​​using the DLHPDE model; improving the DLHPDE model using the finite element method and solving partial differential equations based on different weight values; and outputting a denoised image. The present invention utilizes the DLHPDE model in partial differential equation image processing and uses the finite element method to obtain a numerical solution to the model, thereby resolving the weak marginalization problem in the prior art during the denoising process of glass product image inspection. The finite element method has high computational accuracy and can approximate the boundaries of regions with complex geometric shapes, effectively improving the efficiency of finite element pre-processing.

The invention provides a medium loading resonant cavity analysis method based on a volume-surface integral equation, and belongs to the technical field of resonance characteristic numerical analysis. According to the method, an internal resonance generalized eigenvalue equation based on a body-surface integral equation and internal resonance factors representing internal resonance attributes of body and surface current characteristic modes are defined, so that accurate solution of a unique internal resonance mode is realized; the method avoids the problem that the traditional characteristic value method is difficult to directly solve the internal resonance mode and the mode field of the metal cavity with the aperture after the medium is loaded, is suitable for the loading medium and the metal cavity with any shape and placement relation, and achieves accurate modeling through the volume-surface mixed grid on the premise of guaranteeing the solving precision. Compared with a finite element method, the method has lower modeling complexity.

This invention provides an information processing method that can shorten the computation time required to obtain results equivalent to those obtained by a Pareto front using only finite element method calculations in the optimization calculation of design parameters using a genetic algorithm with CAE analysis such as the finite element method. [Solution] An information processing method for optimizing design parameters that define an individual model representing a design object using a genetic algorithm, wherein the computer calculates the response value of the individual model in each generation using a first calculation method or a second calculation method which has a faster calculation speed than the first calculation method. When the computer is calculating the response value of the individual model using the second calculation method, it calculates an optimization progress index that indicates the progress of the design parameters, and determines whether the optimization of the design parameters has stalled based on the amount of change in the optimization progress index with respect to generational progress, and if it determines that the optimization has stalled, it calculates the response value of the individual model using the first calculation method.

The present application relates to the technical field of aircraft engine simulation design, and particularly relates to a method for estimating aircraft engine thrust based on flight parameters. The aircraft engine thrust estimation method provided by the present application estimates the aircraft engine thrust from the perspective of flight parameter data, in combination with aircraft engine design parameters, taking into account task requirements and adopting the finite element method. Compared with the traditional scheme based on data, i.e., relying on a large amount of flight parameter data and combining statistical rules to obtain the aircraft engine thrust spectrum, the aircraft engine thrust estimation formula derived based on the physical formula has more reliable physical meaning and physical constraints, the distortion degree of the obtained aircraft engine thrust spectrum caused by data sampling error is smaller, the trend and characteristics are more in line with the real rules, and the aircraft engine thrust spectrum has more guiding significance and prediction value, thereby providing strong theoretical support for aircraft engine thrust spectrum prediction.

This invention relates to the field of slope stability analysis technology, and in particular to a method, device, and storage medium for measuring the stability of unsaturated soil slopes. The unsaturated soil slope stability measurement method of this invention is based on the soil-water characteristic curve of unsaturated soil and uses the finite element method to solve the Richards equation to obtain the transient seepage field distribution of the unsaturated soil slope. This method is more consistent with reality than the seepage field assumed by empirical formulas in the past. While the SMP criterion proposed by Songgang-Zhongjing considers the influence of intermediate principal stress, the SMP criterion is only suitable for cohesive soil failure criteria. Therefore, this invention modifies the SMP criterion by considering the influence of cohesion, friction angle, and matrix suction, making it applicable to unsaturated cohesive soil slopes. Based on the calculation of the slope's self-weight stress, and using a modified SMP material model, we conduct stability analysis of unsaturated soil slopes using the strength reduction method, considering the influence of the intermediate principal stress of the soil, resulting in more reasonable results.

The application relates to a finite element numerical implementation method of an asphalt mixture fractional viscoelastic model, and belongs to the field of constitutive behaviors of asphalt mixtures. The application aims at the problem that the existing finite element method of the fractional viscoelastic model cannot consider both the calculation efficiency and the accuracy. The application comprises the following steps: obtaining a complex frequency domain modulus expression of the fractional viscoelastic model based on Laplace transformation; constructing a one-dimensional form of a stress-strain minimum-order differential expression of the asphalt mixture fractional viscoelastic model through time domain / complex frequency domain equivalent substitution; constructing a discrete form of a stress updating expression of the asphalt mixture fractional viscoelastic model based on non-classical method discrete fractional derivative terms in the stress-strain differential form; constructing a three-dimensional form of stress updating expressions and Jacobian matrix updating expressions of the fractional viscoelastic model based on common / abnormal Poisson ratio assumptions and an incremental iteration method; and writing a user subroutine UMAT for finite element numerical calculation. The application is used for finite element numerical implementation of the fractional viscoelastic model.