132 results about "Discrete element method" patented technology

The invention discloses a mine semi-autogenous grinding process analog simulation method and system based on discrete elements, and relates to the field of analog simulation, and the method comprises the steps: according to a crushing characteristic label and a parameter mapping matrix, fusing Internet of Things sensor data, constructing a particle size distribution prediction model and an energy consumption optimization scheduling strategy, and obtaining an energy-saving operation plan; comparing the simulation data set with field production data, adjusting multi-phase coupling particle model parameters, generating a calibration parameter set, researching the influence law of semi-autogenous mill operation parameters by adopting a control variable method based on the calibration parameter set, obtaining a performance index database and an analysis report, and fusing multi-sensor data and laser radar data to obtain a semi-autogenous mill operation parameter analysis result. An ore grinding process sensing mechanism and optimization control logic are constructed, and a visual interface is generated; through multiphase modeling and data driving optimization, the simulation precision and the process suitability are comprehensively improved, and an efficient and economical solution is provided for the mine semi-autogenous grinding process.

The invention relates to the field of powder material manufacturing, and discloses a powder particle model parameter calibration method based on a discrete element method.The method comprises the steps that a static accumulation angle experiment is conducted, and the accumulation angle of micron-sized powder is obtained; performing particle size scaling according to a powder particle size distribution curve; discrete element intrinsic parameters are set according to powder characteristics and mold materials; setting to-be-calibrated parameters and ranges; performing PB design according to the to-be-calibrated parameters and the range; analyzing a test result and screening test parameters which have relatively obvious influence on the simplified and scaled powder accumulation angle; establishing a regression model and seeking an optimal combination by taking a relatively significant test parameter as an independent variable and taking an actually measured static accumulation angle as a target value; and finally determining the mineral particle model parameters based on the discrete element method through the established optimization combination verified by a numerical test. According to the method, the parameter calibration efficiency and the parameter calibration accuracy are improved, and an effective technical means is provided for key parameters of a powder particle numerical simulation experiment.

The invention discloses a direct coal liquefaction pump thermal wear prediction method considering a high temperature effect, and belongs to the technical field of fluid machinery and reliability engineering.The method comprises the steps that the high temperature effect is considered, a liquid phase control equation and a solid particle dynamic model are subjected to coupling iteration, and a solid-liquid two-phase flow numerical model is obtained; the method comprises the following steps: performing full flow field modeling on a fluid domain formed by flow passage components of a pump for direct coal liquefaction, and performing grid division on the fluid domain; constructing a heat-flow-solid multi-field coupling calculation system, and solving the solid-liquid two-phase flow numerical model by adopting a large eddy simulation method based on temperature correction and a discrete element method coupling method to obtain information of flow fields and particle fields at all positions in the pump; operation condition parameters are obtained, the thermal wear rate is obtained through a thermal wear rate calculation formula on the basis of flow field and particle field information of all positions in the pump, and the predicted wear amount of the flow passage component is obtained; and the prediction precision and applicability of the wear model are improved.

The invention provides a CFD-DEM coupling acceleration method, which is applied to the field of computational fluid dynamics and discrete element method coupling simulation and comprises a CFD module, a DEM module, a coupling controller and an optimization module. The DEM module predicts particle distribution based on LSTM to generate a probability distribution diagram; a gradient is calculated in combination with flow field data in a CFD module, and an irregular dynamic coupling region is divided by using a geometric deep learning model; constructing a GNN-based proxy model, and adaptively switching the proxy model and a traditional drag force model based on Reynolds number; the coupling controller monitors a calculation state, and triggers data synchronization of a hierarchical compression strategy when an asynchronous time step threshold value is reached; the optimization module monitors a momentum error after data synchronization, and if the momentum error exceeds a threshold value, a neural network parameter and a threshold value coefficient are fed back and adjusted; according to the method, the problems that a traditional method is low in calculation efficiency, rigid in region division and lack of error control are solved, the calculation efficiency is remarkably improved by 5-10 times, the data synchronization amount is reduced by 58-65%, and it is guaranteed that the precision error is smaller than or equal to 5%.

The invention belongs to the technical field of metal material mechanical property research, and discloses a gear shot peening tooth surface residual stress calculation method and system.The method comprises the steps that a shot speed calculation model in the gear shot peening process is built, shot peening process analog simulation is carried out, and dynamic information of collision between shots and a tooth surface is output; the projectile speed is converted from a machine tool coordinate system to a local coordinate system; establishing a sprayed target plate model and a projectile model, dividing grids, and assembling the grid-divided projectile model and the sprayed target plate model to construct a random multi-shot finite element geometric model; and performing finite element simulation calculation on the random multi-shot finite element geometric model, and outputting residual stress distribution of the gear. According to the method, computational fluid mechanics and a discrete element method are coupled, a high-precision calculation model of the shot speed in the gear shot blasting process is constructed, the forming mechanism and distribution characteristics of the tooth surface residual stress in the actual shot blasting process are effectively reflected, and the method is good in economical efficiency, high in simulation precision, high in process executability and closer to engineering practice.

The invention belongs to the field of geologic disaster prevention and control of geotechnical engineering, and provides a multi-scale numerical simulation method for a slope instability mechanism of a ductile shear zone, which comprises the following steps: acquiring physical and mechanical parameters of a rock-soil body material, and establishing a numerical model based on a particle flow discrete element method; a Johnson-Kendall-Robert contact model is adopted as a contact model between particles, the contact force between the particles is calculated, and meanwhile, a pressure interference parameter is defined to control the limit value of normal contact rigidity; performing dynamic modeling of the discrete element particle system based on a particle generator; static modeling of a discrete element particle system is carried out based on a Lagrange body continuous unit component; the acting force of discrete element particles is converted into equivalent grid node force through an equivalent node force conversion method, and surface attributes of all unit sets possibly involved in contact are defined in a solver.

The invention discloses a method for calibrating finite-discrete element input parameters of a water-containing rock, which comprises the following steps of: respectively carrying out a uniaxial tensile test, a uniaxial compression test and a triaxial compression test on a dry rock sample and rock samples with different water contents; obtaining an evolution equation of the tensile strength along with the moisture content, an evolution equation of the elastic modulus and the uniaxial compressive strength along with the moisture content, and an evolution equation of the cohesive force and the internal friction angle along with the moisture content; carrying out a lateral constraint water absorption swelling test on the dry rock to obtain a vertical expansion evolution curve; on the basis of the obtained results, the fracture energy, the humidity expansion coefficient and the humidity diffusion coefficient in the finite-discrete element method are calibrated; on the basis of the calibration parameters, finite-discrete element uniaxial and triaxial compression numerical simulation of the rock samples with different water contents is carried out, and simulation results are obtained; and when the error between the test result and the simulation result is within a preset threshold value, obtaining a final parameter calibration result. The invention provides a simple, accurate and efficient parameter calibration means for water-rock interaction mechanism research.

The invention relates to a karst area rock sample penetration damage simulation method and system based on discrete elements, and belongs to the field of geotechnical engineering. The invention discloses a karst area rock sample penetration damage simulation method based on discrete elements. The karst area rock sample penetration damage simulation method comprises the following steps: S1, constructing a fracture damage constitutive model reflecting a karst area rock discontinuous surface position; s2, inputting parameters and conditions based on a discrete element method and the fracture damage constitutive model reflecting the position of the rock discontinuous surface in the karst area, and calculating particle stress information; s3, calculating the particle speed and displacement according to the particle stress information to obtain a calculation result; s4, outputting the calculation result; wherein in the step S2, according to the position relation between the particles and the discontinuous surface, the damage coefficient of the discontinuous surface is introduced, the particle stress information is dynamically corrected, and the penetration mechanical response under different discontinuous surface positions is accurately simulated. According to the method, the position and the crushing damage degree of the discontinuous surface of the rock are considered, a new rock sample model is constructed, and more accurate simulation is realized.

The invention relates to the field of remote control, and discloses an intelligent bucket throwing operation method and system for remotely controlling a grab ship unloader, which are used for realizing intelligent bucket throwing operation for remotely controlling the grab ship unloader. The intelligent bucket throwing operation method for remotely controlling the grab ship unloader comprises the steps of generating a full-state expected trajectory containing position, speed and acceleration information, ensuring that the trajectory is smooth and accurate, constructing a nonlinear model prediction controller, integrating a nonlinear disturbance observer, effectively estimating and compensating wind disturbance and unmodeled dynamics, and enhancing the anti-jamming capability of the system. A discrete element method is utilized to simulate interaction of material particles, and in combination with bilge topographic data, spatial distribution after material throwing is accurately predicted, system stability is monitored in real time, and a control instruction is dynamically adjusted. The precision, the stability and the efficiency of bucket throwing operation of the grab ship unloader are remarkably improved, the operation difficulty and the labor cost are reduced, and the grab ship unloader can well adapt to complex working conditions.

The invention belongs to the technical field of slope engineering numerical simulation research, and discloses an engineering scale earth-rock mixed slope numerical modeling and critical dimension determination method and system, a terminal and a storage medium, and the method comprises the steps: generating geometric parameters of stones through employing a discrete element method, and building an earth-rock mixed slope finite element model according to the geometric parameters; the safety coefficient of the earth-rock mixed slope finite element model is calculated, the critical dimension of the earth-rock mixed slope is determined based on a hypothesis testing method, and a numerical result needed for calibrating the earth-rock mixed slope critical dimension calculation model is obtained; calibrating model parameters of the critical dimension calculation model based on the preset number of numerical simulation results; stone information of the soil-rock mixed slope in the engineering scale test is obtained; and substituting the obtained stone information into the soil-rock mixed slope critical dimension calculation model, and solving and outputting the engineering scale soil-rock mixed slope critical dimension. According to the method, the size of the representative unit of the soil-rock mixed slope under the engineering scale can be quickly and accurately determined.

This invention discloses a method and apparatus for simulating the motion of a projectile landslide. The method includes: acquiring rock mass parameters of the projectile rock mass; calculating the initial state parameters of the projectile rock mass during the landslide; analyzing the flight process of the projectile rock mass in the air using an aerodynamic analysis model based on the rock mass parameters and the initial state parameters to determine the state parameters of the projectile rock mass before impact with the ground; and analyzing the motion process of the projectile rock mass after impact and after impact using a discrete element method (DEM) model based on the state parameters before impact with the ground to obtain the final shape of the projectile rock mass. This embodiment, by simultaneously utilizing both an aerodynamic analysis model and a DEM model, can more realistically simulate the motion process of a projectile landslide.

The invention discloses a dynamic adding control method of calcium oxide in flue gas desulfurization, which relates to the technical field of environmental engineering, and comprises the following steps: constructing a computational fluid mechanics and discrete element method coupling simulation model, and introducing a humidity parameter field, a particle size parameter field and an agglomeration degree parameter field into the model, and generating a spatial distribution map and a blockage probability map of the calcium oxide particles in the conveying path, and generating a conveying path control priori according to the spatial distribution map and the blockage probability map. Through coupling simulation and multi-array sensing, accurate recognition of particle aggregation and targeted disturbance control are achieved, the effective adding efficiency of calcium oxide is improved, model parameters are optimized in combination with by-product quality feedback, dynamic closed loop and predictive control are achieved, the desulfurization efficiency is remarkably improved, resource waste is reduced, the by-product quality is stabilized, and the method is suitable for industrial production. And the intelligent level of the system is enhanced.

The invention discloses a CFD-DEM coupling and atomizing nozzle-based tobacco strip moisture regaining processing strength control method, which accurately controls a moisture regaining process through numerical simulation and process parameter optimization, and improves the product quality and the production efficiency. The method combines a bidirectional coupling theory of computational fluid dynamics (CFD) and a discrete element method (DEM) to simulate the movement, heat transfer and humidifying process of tobacco flakes in a moisture regaining cylinder, particularly, on the basis of a traditional CFD-DEM model, an atomizing nozzle model is innovatively introduced, a nozzle structure and an atomizing mechanism are optimized, uniform distribution of moisture is improved, and the moisture regaining effect is improved. The moisture content and temperature balance of the tobacco lamina are ensured, so that the moisture regaining effect is improved. Meanwhile, the invention provides a quantitative optimization method based on the processing strength (H), the relation between the processing strength and parameters such as temperature, moisture content and convective heat transfer coefficient is analyzed, the processing strength is kept in the optimal range by adjusting the parameters, and the stability of the moisture regaining process and the consistency of the quality of the tobacco lamina are ensured.

The invention provides a tailing particle classification method based on multi-dimensional morphological feature compression and clustering, and belongs to the technical field of geotechnical engineering, mineral processing and particle classification, and the method comprises the following steps: S1-S6, obtaining an STL file database of tailing particles; s7, calculating and extracting morphological parameters to form an original parameter data set; s8-S10: obtaining a three-dimensional morphological parameter set which retains the original morphological feature information; s11-S12, determining an optimal clustering number by adopting an elbow rule, a contour coefficient and Gap statistics, and performing particle morphology clustering grouping on the three-dimensional morphology parameters by adopting a K-means algorithm; and S13, carrying out stacking test simulation on the clustered particle categories by using a discrete element method, analyzing the stacking density and the average coordination number of different categories of particles, and verifying the physical significance of a classification result. The tailing particle classification method can systematically and efficiently integrate particle multi-dimensional form information, can directly associate the information with macroscopic physical and mechanical properties, and is data-driven, clear in mechanism and explainable in result.

This invention belongs to the field of oil and gas production enhancement and stimulation technology. It provides a FDEM simulation method and system for mineral-scale shale hydration-induced fracturing. The method includes the following steps: obtaining target reservoir stress, rock porosity, and permeability; preparing rock samples from the target stratigraphic layer and obtaining the mineral composition and content through X-ray diffraction experiments; constructing a polycrystalline geometric model of the minerals using random Voronoi mosaic technology; conducting core immersion hydration experiments to obtain the relationship between core water content and time; setting initial and boundary conditions for the model and establishing a numerical model of clay mineral hydration expansion-induced fracturing based on the finite discrete element method; and conducting numerical experiments on shale hydration-induced fracturing to obtain the micro-fracture propagation morphology induced by clay mineral hydration expansion. This invention considers the non-uniform mosaic distribution of multiple minerals in the reservoir rock, simulates the non-uniform stress distribution induced by clay mineral water absorption and expansion, and realizes numerical simulation of the mineral-scale shale hydration-induced fracturing process.

The application discloses a short-fiber reinforced ablative material bulk ablation simulation method, and belongs to the technical field of multifunctional composite materials and structure numerical simulation.The application comprises the following steps: constructing a discrete system composed of a series of spherical particle units by using a discrete element method; calculating the contact area, contact force and connecting force of the particle units by using a Hertz-Mindlin contact model and a parallel-bond model at each calculation time step; calculating the temperature of the particle units by using an energy conservation principle and a Fourier heat conduction law; calculating the displacement and angular displacement of the particle units by continuously integrating the momentum conservation equation twice; updating the velocity, displacement and diameter of each particle unit, and updating the porosity and pore pressure of the discrete system; iteratively solving and visualizing until a specific condition is met, and ending the simulation.The application characterizes the internal structure features of the short-fiber reinforced ablative material by using the discrete element method, so that the numerical simulation accuracy of the ablation process is improved, and the possibility of capturing the micro ablation mechanism is provided.

A kind of block discrete element-finite element combined modeling gas enrichment zone intelligent identification method, collects the stratum geological data and mine production data of abandoned mining area;Obtain stratum mechanical parameters and gas occurrence parameters;Establish a three-dimensional mine geological model in block discrete element UDEC;Carry out excavation, balance simulation of stoping section, and extract the stope fracture image after excavation;The stope fracture image is imported into image recognition software Coreldraw, and the coordinates of each stope fracture are obtained;A space grid model is established in MATLAB, and the permeability is given by using the relationship between the coordinates of the stope fracture and the space grid;The obtained permeability is imported into finite element software COMSOL in the form of interpolation method for gas flow simulation, and the gas emission and enrichment area are comprehensively determined. This method can accurately determine the enrichment area of underground gas in abandoned stope, and can provide scientific guidance for the layout of ground gas extraction well.

The invention belongs to the field of geotechnical engineering, and relates to a method for determining particle strength distribution in a rockfill material humidifying process. According to the method, the change relation of the moisture content of the particles along with the humidifying time is tested in a particle soaking humidifying test, the probability distribution of the particle strength under different humidifying durations is given, and a probability distribution function of the particle strength under different moisture content conditions and a parameter determination method thereof are provided. The method lays a foundation for simulating the rockfill material humidifying process through a discrete element method, tracking particle structure evolution in the humidifying process, exploring the rockfill material humidifying deformation mechanism and finely analyzing the stress deformation state in the rockfill dam humidifying process. The method can be further applied to hydroelectric and wind-solar complementary power stations to determine reservoir water level operation schemes and provide technical support for increasing benefits of the complementary power stations.

The invention discloses a seed cultivation process optimization method based on a discrete element theory, and belongs to the technical field of seed cultivation. The method comprises the following steps: establishing a mechanical interaction model between soil particles and seeds according to an acquired physical characteristic parameter set by utilizing a discrete element method; acquiring a target physical characteristic parameter set, inputting the target physical characteristic parameter set into a mechanical interaction model, and performing discrete element analogue simulation output under a dynamic simulation condition to obtain a change curve of contact force along with seed growth displacement; and judging whether the contact force in the change curve meets a preset threshold value in each stage of seed growth, if not, continuously self-optimizing the target physical characteristic parameter set, and ending optimization until the judgment result is yes. According to the method, a seed optimization design scheme can be generated, improvement strategies for different soil types can be formulated, and a reliable basis is provided for agricultural machinery design optimization and seeding efficiency improvement.

The invention discloses a numerical simulation method for soil erosion caused by rainfall. The method comprises the following steps: establishing a numerical model and boundary conditions of a simulation area; modeling a soil particle system by adopting a discrete element method; modeling a wind and fluid-solid interaction system by adopting computational fluid dynamics; modeling rainfall by adopting a discrete element method; performing CFD-DEM coupling solution and erosion process simulation, calculating movement of all soil particles and raindrops under related acting force in each time step, and updating positions and speeds of the movement; according to the updated particle positions, the porosity in each CFD grid unit is calculated; based on the relative velocity of the particles and the fluid; fluid acting force borne by each particle is calculated and aggregated into a source item in a CFD grid unit; substituting the calculated porosity and the source item into the model, and solving the updated fluid velocity field and pressure field; key interaction judgment and implementation are carried out; iteratively circulating until the set total simulation time is reached; and outputting related data of the soil particles in the simulation process.

The application discloses a recycled aggregate crushing simulation and optimization method based on a DEM particle crushing algorithm, which is based on the Griffith strength theory and the discrete element method, realizes efficient and accurate simulation of particle crushing behavior through single-particle three-dimensional stress calculation, crushing criterion establishment, fracture surface construction based on a contact point origin and sub-particle attribute inheritance, and forms a method from parameter calibration to dosage quantification; the application quantifies the maximum dosage of recycled aggregate under different highway grades and base types by revealing the mesoscopic mechanism that recycled mortar is similar to natural crushed stone in mechanical properties and recycled bricks have a decisive control on the crushing value of the mixture, and provides a theoretical reference and practical guidance for the resource high-value utilization of building solid waste recycled aggregate in road bases.

This invention discloses a discrete element method, system, medium, and equipment for simulating equivalent crystalline materials, relating to the field of numerical simulation of rock mechanics. The method includes: converting granite images into binary images using a binarization method to identify various minerals within the binary images; determining the size data of various minerals in the binary images based on a particle size statistical algorithm; extracting micromechanical parameters using nanoindentation experiments combined with a trial-and-error method, and constructing multiple sets of equivalent crystalline models based on the mineral size data; dividing each set of equivalent crystalline models into multiple regions using a Brazilian splitting experiment; determining the optimal magnification for each region of each set of equivalent crystalline models by magnifying from non-critical regions to critical regions using a magnification algorithm, obtaining differentiated equivalent crystalline models; and performing Brazilian splitting numerical simulation using these differentiated equivalent crystalline models to obtain simulation results.

This invention discloses a numerical simulation method for the maintenance and construction of semi-flexible pavement based on dual-scale coupling, comprising: determining the dimensions of the milled and repaved section based on the original pavement structure obtained by three-dimensional ground-penetrating radar and falling-weight deflectometer; establishing a finite element model for the original pavement structure; simulating the semi-flexible layer after milling and repaving using discrete element method and coupling it into the finite element model to construct a dual-scale coupled model; performing simulation calculations on the compaction and molding of large-void parent asphalt mixture; conducting simulation derivation of self-leveling cement slurry filling based on the construction process and the dual-scale coupled model; establishing a discrete element numerical model for cement curing to characterize cement curing behavior; and evaluating and guiding the maintenance of the semi-flexible pavement based on the dual-scale coupled model. This invention realizes process control and performance evaluation during the construction period of semi-flexible pavement and guidance for maintenance and repair during service life, making up for the shortcomings of high testing costs and insufficient theoretical data for different construction conditions in semi-flexible pavement.

The invention relates to a concrete structure fault movement damage simulation method based on damage contact constitutive, which belongs to the technical field of geotechnical engineering numerical simulation, and has the technical key points that the method is based on a plastic damage constitutive theory proposed by A.D.Jefferson and is combined with a discrete element method to simulate the fault movement damage of a concrete structure. A contact model suitable for cracking evolution of the concrete lining structure under the fault shearing condition is constructed; the method comprises the following steps: establishing a contact surface force-displacement increment relationship and a failure criterion by introducing parameters such as normal rigidity, tangential rigidity, tensile strength, cohesive force, a friction angle and fracture energy, carrying out secondary development coding based on a 3DEC development port by utilizing a C + + language, and generating a constructed DLL module file; by constructing a typical engineering numerical model, the development process of concrete cracks under fault movement is simulated, and sensitivity analysis is carried out. The method can accurately capture a structural damage path and a stress concentration part, and is suitable for engineering application scenes such as tunnel crossing fault and high-stress area support optimization.

The invention relates to the field of rock mechanics and mine engineering, and discloses an anchor rod prestress applying simulation method, system and equipment and a storage medium. The method comprises the steps that surrounding rock excavation simulation is conducted on a surrounding rock numerical model, and after the surrounding rock numerical model reaches a balanced state, an anchor rod is activated and anchor rod parameters are input; prestress is applied to the surrounding rock numerical model, the state of the anchor rod is modified through a mixed finite element-discrete element method, and the anchor rod is divided into a shallow section and a deep section which are disconnected; oppositely stretching the shallow section and the deep section according to the pre-tightening force data until the numerical model of the surrounding rock reaches the balance state again; parameters of the anchor rod are modified, the shallow section and the deep section are bonded to obtain a complete anchor rod, and a prestress simulation result is generated through an anchor rod constitutive equation and a surrounding rock constitutive equation. According to the invention, the real process of forming the near-pressure and far-tension prestress in the surrounding rock by oppositely stretching the anchor rod can be reflected, and the interaction process between the pre-stressed anchor rod and the surrounding rock can be accurately captured.

The application discloses a kind of inversion methods of tunnel drill and blast method construction blasting load curve based on discrete element method, it is related to tunnel blasting technique numerical simulation technical field, specifically in solving the problem that discrete element method is applied to the design and optimization of tunnel blasting load in background technology, realize the significant problem of tunnel drill and blast method construction, by designing a kind of inversion methods of tunnel drill and blast method construction blasting load curve based on discrete element method, the method is combined with high-precision rock mechanics model and advanced optimization algorithm, from actual blasting result, the best blasting load parameter is accurately inverted, using the method of the application will significantly improve the safety and economy of tunnel drill and blast method construction, it has important practical significance and wide application prospect to tunnel engineering field.

This invention provides a tailings particle classification method based on multi-dimensional morphological feature compression and clustering, belonging to the fields of geotechnical engineering, mineral processing, and particle classification technology. The method includes: S1-S6: acquiring an STL file database of tailings particles; S7: calculating and extracting morphological parameters to form an original parameter dataset; S8-S10: obtaining a three-dimensional morphological parameter set retaining the original morphological feature information; S11-S12: determining the optimal number of clusters using the elbow rule, silhouette coefficient, and Gap statistic, and using the K-means algorithm to cluster the three-dimensional morphological parameters into particle morphology groups; S13: using the discrete element method to conduct a stacking test simulation on the clustered particle categories, analyzing the stacking density and average coordination number of different particle categories, and verifying the physical meaning of the classification results. This tailings particle classification method can systematically and efficiently integrate multi-dimensional particle morphological information and directly correlate this information with macroscopic physical and mechanical properties. It is data-driven, has a clear mechanism, and the results are interpretable.

The invention discloses a manned lunar vehicle elastic wheel traction performance rapid evaluation method and device, and the method comprises the steps: building a virtual lunar soil model based on a discrete element method, and obtaining the ground mechanics parameters of the interaction of lunar soil and wheels through a virtual pressing plate test and a virtual shear test; on the basis of the terramechanics parameters and the wheel parameters, the wheel soil acting force is calculated through an improved semi-empirical model, so that the traction performance of the wheels is evaluated; solving a wheel soil contact parameter by adopting an iteration mode, and evaluating the traction performance of the wheel according to the calculated wheel soil acting force; wherein the improved semi-empirical model divides the contact area of the wheel and the lunar soil into an elastic contact area with elastic deformation and a rigid contact area without deformation, the wheel soil acting force in the two areas is calculated respectively, and the wheel structure parameters and the lunar soil parameters are analyzed in a simulation mode according to the wheel structure parameters and the lunar soil parameters. And the traction performance such as hook traction force, rolling resistance and rolling efficiency of the elastic wheel can be quickly calculated.

The invention discloses a level set-based accelerated simulation method for rock-soil granular materials in any shape. Comprising the steps of establishing a level set particle model of any shape; obtaining particle data through the shape template; according to an extended multi-level Linked-cell method, constructing a potential contact pair list; calculating point-to-surface contact information on the vertex of the particle model; according to the normal contact force, tangential contact force and contact torque are calculated, and the particle motion state is updated; and outputting analog information. Compared with a traditional discrete element method, the high-performance level set discrete element algorithm for the particles in any shape can consider more extreme particle gradation distribution and complex particle shapes, has good calculation efficiency and can solve the problem of large-scale particle material simulation; and by establishing the shape template and the Verlet list, memory consumption and memory access delay are reduced, and the computing performance of the GPU can be fully utilized.

The application provides a coal mine underground hydraulic cavity forming process simulation optimization method and system, obtains target coal seam parameters, establishes a three-dimensional coal seam structure, performs model particle filling, and forms a coal numerical model; sets cavity forming parameters, establishes a model three-dimensional structure of the cavity forming system, generates abrasive particles and liquid particles according to a ratio, and sets boundary conditions; performs parameter calibration and balances the model among the particles of coal, water and abrasive; performs numerical simulation calculation of a hydraulic cavity forming process coupled with a discrete element method and computational fluid dynamics, simulates a hydraulic cavity forming operation process, and records the distribution of local stress, density and porosity in a target cavity forming area in the simulation process; according to multiple numerical simulation calculation results under different cavity forming parameters, the response relationship of pressure relief and permeability improvement effects to each hydraulic cavity forming parameter is analyzed, based on the analysis results, an optimization scheme of the hydraulic cavity forming parameters is determined, and a final hydraulic cavity forming process scheme is obtained. The application improves the accuracy and applicability of simulation.