92 results about "Computational mechanics" patented technology

The invention relates to the technical field of computational mechanics and geological disaster simulation, and discloses a full-process numerical simulation method suitable for simulating continuous medium damage to discontinuous medium movement. The invention provides a novel chain-type disaster simulation method for solving the problems that calculation efficiency and large deformation precision are difficult to consider at the same time and the evolution process of a slope from a continuum to a fragmented body cannot be dynamically reflected in the prior art. The core of the method is that a dynamic criterion is set based on a unit real-time damage variable, a failed finite element (FEM) unit is adaptively converted into smoothed particle hydrodynamics (SPH) particles, and the physical state of the particles is accurately mapped; and a virtual particle coupling algorithm is adopted to process a two-domain interface, so that bidirectional transmission of mechanical parameters is realized. The method is executed circularly, the whole process of'continuous deformation-fragmentation flow 'of the side slope is simulated dynamically, high-precision integrated simulation of a'seismic source-propagation-response-movement' disaster chain in the same frame is achieved, and an efficient tool is provided for disaster risk assessment and prevention and control.

The invention discloses a multi-field coupling concrete dam parameter intelligent inversion method, system and device based on deep learning and a medium, and belongs to the technical field of dam engineering safety monitoring and computational mechanics crossing. The geometry and boundary modeling unit is used for uniformly expressing a dam body-dam foundation geometric region and a boundary type; simultaneously representing a state field and a parameter field by utilizing PINN inversion, calculating a physical residual error, and implementing parameter physical feasibility constraint and prior regularization; the organization optimizer, residual drive sampling and weight self-adaption are used for training scheduling, and result review, visualization and data export are carried out through the verification and output unit. According to the method, end-to-end parameter field inversion is realized on the premise of not depending on frequent mesh generation, the calculation cost of multiple forward modeling is reduced, the non-uniqueness of inversion is inhibited, and the method has higher recognition capability for engineering concerned areas such as permeation channels and degraded zones.

The invention belongs to the technical field of computational mechanics and artificial intelligence crossing in civil engineering, and discloses a foundation pit retaining wall deformation prediction method based on a physical information neural network, and the method comprises the following steps: 1, building a retaining wall physical information neural network agent model based on multi-source data fusion; 2, building a displacement-force inversion model of the retaining wall; and 3, carrying out iterative coupling solution on the support system. According to the method, mixed training is carried out by fusing physical rules and field measured data, so that the physical information neural network agent model follows a basic mechanics principle and can also fit specific engineering practice, the prediction result is more reliable, and compared with traditional finite element analysis, the efficiency is improved by several orders of magnitude; therefore, parameter optimization and real-time safety evaluation based on a large amount of calculation become possible.

The invention discloses a finite element-based carbon steel corrosion modeling method, which relates to the technical field of computational mechanics and finite element analysis, and comprises the following steps: S1, acquiring real-time data from an environment sensor, acquiring an initial parameter set containing temperature, salinity and material attributes, and determining initial corrosion rate distribution; s3, if the deviation value exceeds a preset threshold value, a parameter coupling relation is processed through a genetic algorithm, and an optimized parameter set is obtained and used for correcting model input; according to the finite element-based carbon steel corrosion modeling method, the precision and reliability of corrosion prediction in a complex environment are remarkably improved, and a scientific basis is provided for carbon steel structure maintenance.

The invention belongs to the technical field of crossing of artificial intelligence and computational mechanics, and particularly relates to an artificial intelligence-based MPM simulation method and system for a mold bag sand cofferdam, a three-dimensional simulation system of the mold bag sand cofferdam is constructed in combination with an MPM method, and the penetration state of sand and geotechnical cloth is detected in an accelerated manner by utilizing a penetration tracking algorithm and spatial Hash. An artificial intelligence algorithm is introduced, a reward function is designed with minimum penetration depth and penetration times as targets, penetration prediction of MPM particles is driven, the speed and stress of the MPM particles are further dynamically regulated and controlled on the basis of a predicted contact model and the artificial intelligence algorithm, penetration is prevented, numerical stability is ensured, and the prediction accuracy of the MPM particles is improved. The calculation precision is improved by adopting implicit integration and gradient descent methods; and bidirectional coupling of the sand and the geotechnical cloth is realized through a visualization module. The problems that in a traditional method, coupling simulation of sandy soil and a flexible structure is low in efficiency, easy to penetrate, unstable in numerical value and the like are solved, and the accuracy and reliability of simulation analysis of the mold bag sand cofferdam structure under the complex working condition are remarkably improved.

The invention discloses a surface material point embedding method for a multi-body contact problem, and belongs to the technical field of computational mechanics and numerical simulation. The method comprises the following steps: 1, performing discretization and parameter setting on a physical model of a multi-body contact system; 2, calculating a contact force; 3, transmitting the physical quantity of the material points of the contact body to background grid nodes to obtain node mass, momentum and force; 4, updating the grids; 5, updating the speed, the displacement, the stress and the density of the material point; 6, resetting the physical quantity of the background grid; seventhly, repeating the first step, the fourth step, the fifth step and the sixth step, sequentially updating the physical quantities of the material points of other contact bodies and the positions and the speeds of the slave nodes, and distributing background grid node memories of the areas where the corresponding contact bodies are located; and 8, after all contact bodies are updated, updating the current time to obtain a dynamic simulation result. According to the method, the boundary of the object can be accurately represented, complicated problems such as impact and multi-body contact can be effectively simulated, and meanwhile, the method has good stability and relatively low memory consumption.

The invention relates to the technical field of computational mechanics, in particular to a DDDA-based circle-circle contact interface non-occlusion sliding numerical simulation method, which comprises the following steps: constructing a non-occlusion joint model; establishing a circular unit-non-occlusion joint initial correlation; establishing a DDDA total balance equation; applying a non-occlusion contact constitutive to the associated circle unit contact; and solving the model to obtain the non-occlusion sliding process of the circle-circle contact interface. According to the invention, the non-occlusion sliding process of the circle-circle contact interface is realized, so that the sliding simulation result of the structural surface better conforms to the real mechanical response of the structural surface, and a more convincing numerical analysis tool is provided for the sliding process analysis of slope stability, the slip instability prediction of underground engineering surrounding rock along a weak surface and other key engineering problems. And the method has great engineering application value.

The present invention relates to the technical field of the intersection of computational mechanics and deep learning, and specifically to a deep learning method that combines a physical information neural network with a finite element method; specifically, partial differential equations are embedded into a neural network, comprising the following steps: S1, meshing the solution domain, and calculating the stiffness matrix and the node load array; S2, collecting observation data and boundary conditions; S3, using the energy represented by the stiffness matrix and the fitting error of the observation data as a loss function to train the neural network; S4, executing an optimization algorithm to optimize the neural network; S5, and repeating the training and optimization steps until convergence; the deep finite element method proposed in the present invention combines a physical information neural network with a finite element method, and exhibits significant advantages and beneficial effects in terms of the fusion of physical knowledge and observation data, modeling process and loss function design, model training stability and simplicity, computational efficiency and accuracy, pre- and post-processing technology and transfer learning capability, and applicability.

The invention belongs to the technical field of composite material structural mechanics and computational mechanics, and particularly relates to a composite material fan blade composite crack propagation simulation method, program and equipment and a storage medium. According to the method, key rotation and a failure judgment criterion based on a key deflection angle are introduced into the conventional state near-field dynamic model, so that the crack propagation behavior of the composite material under the shear or mixed mode load can be more accurately captured and simulated. According to the flat plate mapping method, the improved PD model suitable for the flat plate is successfully expanded and applied to the fan blade airfoil section, the problem of complex geometric modeling is solved, and the flat plate mapping method is suitable for complex geometric structures. According to the method, a near-field dynamics method is adopted, crack paths do not need to be preset or grid reconstruction does not need to be carried out, and complex expansion behaviors such as spontaneous bifurcation and merging of cracks can be naturally simulated; the PD model can be used for meticulously describing the mechanical behaviors and crack evolution of each layer and interlayer in the laminated plate, including intra-layer cracking and interlayer layering.

The invention relates to the technical field of civil engineering and computational mechanics, in particular to a load and environmental corrosion coupled reinforced concrete structure analysis method. The method comprises the following steps: carrying out structural finite element modeling; carrying out mechanical field solving; the diffusion characteristics of the crack-driven corrosive medium are updated; solving a corrosive medium diffusion field; the material performance is corrected; carrying out double-time-domain self-adaptive solution; by means of the mode, the coupling effect of complex loads and multiple environmental factors in the reinforced concrete structure can be comprehensively and accurately simulated within the engineering acceptable calculation cost, and reliable prediction of the full life cycle performance is achieved.

The invention discloses a solver error control development method and system based on isogeometric analysis, and relates to the technical field of computational mechanics and numerical analysis, and the method comprises the steps: building an NURBS primary function according to isogeometric analysis, carrying out the grid division, generating initial grid data, calculating the stress error index of each grid unit through employing a Zienkiewicz-Zhu method, and carrying out the calculation of the stress error index of each grid unit. Grid optimization parameters are generated, local encryption processing is carried out on grid units, the order of an NURBS primary function is improved by adopting a p-refinement method, and an optimized grid structure is obtained; and inputting the optimized grid structure into a solver for iterative calculation, and when a calculation error is smaller than a preset error threshold value, outputting a solving result. According to the method, the optimal balance between the calculation precision and the resource consumption is realized through the precise geometric expression and the self-adaptive grid division of the NURBS primary function in combination with the Zienkiewicz-Zhu posterior error estimation technology.

The invention provides a composite material fan blade bending fracture forecasting method and device based on near-field dynamics and a storage medium, and belongs to the technical field of composite material structural mechanics and computational mechanics, and the method comprises the steps that a dangerous fan blade section under a preset load is determined; fan blade curves are integrated, an initial plane grid is created, polynomial fitting is carried out on discrete points, plane bending is carried out on transition area smoothing, meanwhile, error control is carried out at the large number position of the blade end, fan blade parameterized modeling is achieved, and the laying layer fiber direction is set; boundary conditions are applied, meshless discretization is carried out based on near-field dynamics, a material point motion equation is obtained, a near-field neighborhood is constructed, a family of each material point is searched, and a family member list is formed. The method is applied to bending deformation simulation of a large structure such as a fan blade section and damage simulation caused by deformation, modeling of the laminated plate is simplified by taking a single-layer material point as a laying layer, and natural crack propagation simulation is realized in combination with critical bending angle damage characteristics.

The invention discloses a method and a device for predicting creep deformation of a nickel-based single-crystal high-temperature alloy under creep-oxidation interaction, and belongs to the technical field of material engineering and computational mechanics, and the method is based on microstructure distribution characteristics of the nickel-based single-crystal high-temperature alloy after a creep test in vacuum and atmospheric environments. According to the method, a dynamic growth multilayer geometric model considering oxidation influence is established, and a method capable of predicting alloy creep deformation in different environments is developed by combining a creep constitutive model considering microstructure evolution based on a test observation result. The method mainly comprises the following steps: providing a multilayer geometric model considering oxidation influence; establishing growth kinetic models of different layers; obtaining microstructure evolution laws at different depths; building a multi-scale creep deformation prediction framework considering microstructure evolution; and acquiring model parameters and predicting creep deformation. The invention provides a set of method capable of accurately predicting the creep deformation of the alloy in different environments.

The invention discloses a synchronous bidirectional coupling solving method for large deformation of a flexible slender beam under uniformly distributed pressure based on an accurate load stiffness matrix, and relates to the field of computational mechanics. According to the method, a local curve coordinate system and a uniformly distributed pressure vector function are constructed through an absolute node coordinate Euler-Bernoulli beam unit discrete beam structure, a unit external force vector is derived based on a virtual work principle, and a precise load stiffness matrix is obtained through a vector value function differential rule. The load stiffness matrix and the elastic tangent stiffness matrix are synchronously assembled in the Newton-Rafson iteration step, and synchronous bidirectional coupling solution of the load and the configuration is achieved. According to the method, the numerical value loading limit of a traditional method is broken through, error accumulation is avoided, the calculation stability and precision are improved, the large deformation and even post-buckling behavior of the flexible slender beam under the uniformly distributed pressure can be accurately analyzed, and the method is suitable for design optimization of related equipment in the fields of ocean engineering, mechanical engineering and the like.

This invention belongs to the field of multi-scale computational mechanics technology for composite materials, and discloses a method and system for generating random fiber distribution RVEs considering random porosity in three-dimensional space. Existing methods do not incorporate the physical mechanism of porosity formation into the control logic of porosity distribution. The fiber random distribution RVE generation method and system proposed in this invention first plans the RVE model parameters; then, in a two-dimensional sketch, a fiber distribution satisfying periodic boundary conditions is generated through mesh layout and random perturbation algorithm; next, the density labels of the fibers are obtained through an improved density peak clustering algorithm, and the pore size is adaptively adjusted accordingly; subsequently, pores are randomly generated in three-dimensional space by spline curve rotation, ensuring periodic distribution; finally, the final RVE model is obtained through Boolean operations. This invention achieves a pore spatial distribution that conforms to actual physical laws, and can effectively generate fiber random distribution RVEs with controllable porosity, providing a method for parameterized modeling of microstructures in FRC calculations.

The invention belongs to the crossing field of computational mechanics and deep learning, and particularly discloses a heterogeneous composite material equivalent elastic modulus prediction method based on cross-scale differential-integral coupling, and the method comprises the steps: building a multi-scale elastic modulus prediction equation; constructing a cross-scale coupling differential equation set; performing data sampling on the multi-scale elastic modulus prediction equation to generate a training sample; designing a multi-channel neural network architecture; calculating an elastic modulus predicted value and an integral operator original function of each order of the defect field function through neural network forward propagation; through an automatic differential technology, a high-order derivative is accurately calculated, and a differential operator is obtained; driving the neural network to perform back propagation through a composite loss function to update network parameters; when the convergence condition is met, the output of the neural network is the solution of the multi-scale elastic modulus prediction equation. By constructing a macroscopic-microscopic dual-scale differential constraint system and a defect field integral error suppression framework, the bottleneck of traditional experimental test and numerical simulation in cross-scale mechanical analysis is effectively broken through.

This invention relates to the fields of computational mechanics and rock engineering, and particularly to a numerical simulation method for end-anchored rock mass engineering based on DDDA (Digital Directed Amplification and Analysis). The method includes: constructing a circular element discontinuous deformation analysis numerical model according to actual engineering conditions; performing discontinuous deformation analysis numerical calculations; analyzing the calculation results of the unanchored model to preliminarily determine the anchoring scheme; performing discontinuous deformation analysis numerical calculations again; determining whether the model contains unstable blocks, adjusting the anchoring scheme, and repeating the aforementioned steps until the model stabilizes, thus completing the design of the rock mass engineering anchoring scheme. This invention establishes a numerical simulation method capable of accurately simulating the interaction between anchor bolts (anchor cables) and rock mass, and integrating anchoring scheme design, stability assessment, and dynamic parameter optimization functions. It achieves accurate prediction of the deformation and failure behavior of anchored rock mass and optimizes anchoring design, featuring high simulation accuracy and strong practicality.

The invention provides a dam internal structural surface automatic identification method and device, and relates to the technical field of hydraulic engineering structure analysis and computational mechanics, and the method comprises the steps: reading a dam finite element grid, and defining boundary control parameters of a target structural surface position; the method comprises the following steps of: identifying a connecting surface simultaneously containing a first type of boundary nodes and a second type of boundary nodes in a model, and performing iterative expansion on the connecting surface by utilizing a topology propagation algorithm to form a continuous connecting surface set; positioning all unit surfaces formed by the first type of boundary nodes based on the set, and integrating the unit surfaces to obtain a target structural surface; and finally, outputting a multi-format grid file through coordinate transformation. According to the method and the device provided by the invention, the node boundary type is defined, the topology propagation algorithm is utilized to automatically identify and iteratively expand the connecting surface, the identification precision is improved, the processing efficiency is improved, the efficient and full-automatic identification of the internal structural surface of the dam is realized, and a core technical support is provided for structural analysis.

This invention belongs to the fields of geotechnical engineering and computational mechanics, and provides a finite element-peripheral dynamic coupling method and system for anchoring effects in underground engineering. The finite element-peripheral dynamic coupling method for anchoring effects in underground engineering includes: establishing a numerical model of the fractured surrounding rock in underground engineering based on periphery dynamics theory; constructing a corresponding anchor finite element model; for each anchor finite element model beam node, combining the numerical model of the fractured surrounding rock in underground engineering, establishing the coupling relationship between the beam node and the surrounding periphery dynamic material points; defining an interface key for each pair of beam nodes and their coupled periphery dynamic material points, and determining the mechanical behavior of the interface key; at each computation time step, performing force transfer and solving between the numerical model of the fractured surrounding rock in underground engineering and the anchor finite element model, updating the coupling relationship between the beam nodes and their coupled periphery dynamic material points, which can accurately simulate the anchoring effect in underground engineering.

The application discloses a simulation analysis method for an aircraft landing gear tire under heavy landing impact load, and relates to the technical field of aeronautical engineering and computational mechanics. First, in the tread and sidewall rubber area, the micro-heterogeneity of the rubber material in the tire is equivalent to the super-elastic mechanical behavior of the macro-continuous medium based on the homogenization theory. This process avoids the extremely dense mesh and complex constitutive relationship required for depicting the micro-features such as filler distribution and interface effect in the traditional model. Then, the rebar element representing force transmission under impact load is created in the pre-constructed rubber matrix grid, and the rebar element is coupled with the rubber matrix. The simplified coupled model can accurately capture the anisotropic stiffness, main force transmission path and reinforcement effect of the cord-rubber composite structure, so as to finally stably and accurately simulate the transient impact process of the tire under heavy landing impact load, which includes the triple nonlinearity of geometry, material and contact.

The application belongs to the field of computational mechanics, and discloses a scale order reduction near-field finite element calculation acceleration method and implementation method for structure fracture analysis.The acceleration method comprises the following steps: firstly, geometric modeling is performed on a target structure; secondly, the split finite element is recorded as a fine scale element and constitutes a fine scale grid; thirdly, based on a spatial neighborhood relationship, the fine scale element is collected into a coarse scale element and constitutes a coarse scale grid; fourthly, an equivalent shape function of the coarse scale element is obtained; fifthly, the overall stiffness matrix and the load vector of the coarse scale grid are calculated; sixthly, the overall node displacement of the coarse scale and fine scale grid is solved; seventhly, it is judged whether the current loading step or time step reaches the maximum loading step or time step; if yes, the displacement nephogram and equivalent damage nephogram of the target structure are output; if not, the key state of all fine scale near-field elements is updated; and finally, the implementation method is that the acceleration method is taken as a core calculation module to be implemented in different scenes.The application significantly improves the calculation efficiency of near-field dynamics in large-scale simulation analysis.

The invention relates to the technical field of underground engineering protection and computational mechanical simulation, and provides an underground structure explosion load dynamic response simulation and damage grading method based on ALE fluid-structure interaction, which comprises the following steps: step 1, unified domain modeling and working condition parameter optimization; 2, selecting a material model and a state equation; 3, initial field, boundary and detonation conditions are applied; step 4, displaying a dynamic ALE remapping solution; step 5, result extraction and multi-index fusion grading; step 6, performing multi-working-condition batch processing and sensitivity analysis; and step 7, engineering calibration of parameters and threshold values. According to the scheme, a three-dimensional ALE fluid-solid coupling simulation and damage grading technical system which can be used for engineering landing, is stable and efficient and can be used for quantitative evaluation is constructed by establishing an ALE fluid-solid coupling calculation framework integrating three domains of explosion gas, surrounding rock and underground structure.

The invention relates to the technical field of structural engineering and computational mechanics, in particular to a method and system for determining the axial pressure bearing capacity of a concrete filled steel tube at negative temperature, which responds to a bearing capacity calculation request, obtains working condition parameters of a target component and constructs a finite element simulation model for simulation analysis to obtain simulation bearing capacity data; the simulation bearing capacity data is compared with a standard theoretical value obtained through calculation according to a preset standard formula, a deviation coefficient is obtained, when the deviation coefficient is larger than a preset deviation threshold value, coefficient fitting is conducted through a preset size effect-negative temperature coupling correction model, a size effect correction coefficient and a negative temperature damage correction coefficient are obtained, and the size effect correction coefficient and the negative temperature damage correction coefficient are obtained. And based on the size effect correction coefficient and the negative temperature damage correction coefficient, the standard theoretical value is corrected, and the final axial compression bearing capacity is obtained. The problem of poor negative temperature adaptability caused by size effect and large-diameter deviation in the prior art is solved, and the structural risk caused by misjudgment of the bearing capacity is avoided.

This invention relates to the fields of computational mechanics and computer-aided engineering, specifically a mathematical model capable of predicting the forming shape of spherical bearings. The model takes the blank's geometric dimensions, mold parameters, and material elastic modulus as input. Based on the precise geometric relationships of the compression process, it establishes a system of equations, deriving the core relationship between the compression amount and the theoretical mean diameter, and solving for the theoretical shape without considering springback. By introducing the theory of compression-bending springback, it establishes springback geometric constraints and expressions for reverse stress. Based on the static equilibrium conditions under pure bending, it derives the explicit relationship between the radius after springback and the radius before springback, elastically correcting the theoretical shape to obtain a predicted shape closer to the actual forming result. The model's output can be compared and verified with finite element simulations and actual experimental data, exhibiting high prediction accuracy in the stable deformation stage with large compression amounts. This invention provides effective guidance for process parameter design and optimization, reducing the number of experiments and development cycles.

The application discloses a kind of composite material surface ablating dynamic mesh processing method and system, belong to the field of computational mechanics and heat transfer, including steps: S1, model initialization and parameter definition;S2, after determining model and parameter, coupling solving material surface ablation and thermal response;S3, after obtaining material surface ablation and thermal response, based on the element deletion criterion and execution of ablation depth;S4, after element deletion, execute dynamic boundary identification and update;S5, time step advances, repeat step S2 to step S4, until reach the preset total length of simulation.This application can not only accurately simulate the non-steady thermal response inside material, but also clearly and intuitively reproduce the macroscopic geometric shape change caused by ablation, and has superiority in handling complex moving boundary problems.

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

The application relates to the field of bridge engineering and computational mechanics, and discloses a nonlinear numerical simulation modeling method for a main cable of a suspension bridge. The method comprises the following steps: a finite element model of the main cable of the suspension bridge is established, and a group of internal state variables are introduced on integral points of cable elements of the finite element model, the internal state variables comprising a section topological configuration field, a friction interface state memory field, a slip activation field and a slip kinematics field; in an incremental step of finite element analysis, according to a preset coupling evolution control equation, the update values of the internal state variables driven by total strain increments are solved; then, the real stress tensor and the tangent stiffness matrix are calculated based on the updated internal state variables, and the calculation is returned to a main analysis program for balance iteration; through the calculation of each incremental step, the tracking and accumulation of the history-dependent response are realized. The application can capture the complex phenomena such as wear, slip and accumulation configuration rearrangement between steel wires, and improves the physical fidelity and calculation stability of numerical simulation.

The invention discloses a computational mechanical analysis system for evaluating bone growth induced by a bone healing scaffold, and relates to the technical field of bone healing evaluation, and the computational mechanical analysis system is specifically characterized in that under unified external load and time gating, quantitative indexes such as a regional average osteogenesis rate curve, an accumulated bone mass, time for reaching a preset healing threshold value, a cross section bridging rate and the like are used for evaluating bone growth induced by a bone healing scaffold; the induction osteogenesis effects of different stent structures are fairly compared, shape scoring and subjective weight are avoided, reproducibility and calibration are emphasized, the method can be used for preoperative structure screening and parameter sensitivity analysis and can also be in butt joint with follow-up images and experimental data, online calibration and treatment course optimization of individualized digital twinning are further supported, and the method is suitable for popularization and application. According to the method, the period boundary, the time-varying supply gating and the multi-physics field coupling of the RVE can be unified into the same calculation framework, and the defects of the prior art in microscopic geometric maintenance, time rhythm expression and biomechanical closed loop are effectively overcome.

The application discloses a reinforced concrete structure damage process simulation method based on an SPH method, belongs to the technical field of computational mechanics and structure engineering simulation, and comprises the following steps: firstly, removing the damage softening component from the total deformation gradient to obtain an effective physical state and a deformation gradient without interference; secondly, determining a plurality of mechanism criteria of tensile cracking, shear slip and compression plasticity exhaustion, and obtaining a comprehensive damage index through maximum value enveloping; mapping the damage index into a softening factor through a nonlinear square term, correcting the deformation gradient and dynamically adjusting the particle action strength; finally, performing stress gradual attenuation on critical damage particles, and removing the particles after meeting the step number condition. The application can consider the calculation efficiency and the macroscopic mechanical property under a single constitutive framework, accurately restores the material damage softening and the fracture process, effectively avoids numerical oscillation, and improves the precision and stability of the structure damage simulation.