73 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 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 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 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.

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.

The invention relates to the technical field of rock material elastic modulus field and Poisson's ratio field prediction, in particular to a rock elastic modulus and Poisson's ratio field prediction method based on physical constraints, which comprises the following steps: constructing a rock mechanical sample data set; the data set comprises a rock elastic modulus field, a Poisson's ratio field, a strain field, a stress field and space coordinates, a physical information neural network is constructed, an encoder, a decoder and a physical calculation layer are constructed, and the physical calculation layer is used for calculating a mechanical response strain field; the method comprises the steps of constructing a multi-element composite loss function fused with physical constraints on the basis of a physical information neural network, performing iterative training on the physical information neural network on the basis of the multi-element composite loss function, and outputting a high-resolution elastic modulus field and a Poisson's ratio field of a target rock test piece. It is ensured that the predicted elastic modulus and Poisson's ratio parameter field meets the basic law of rock mechanics, and the credibility of the prediction result is improved.

The present application relates to the field of digital design and computational mechanics topology optimization of power engineering foundation, specifically to a concrete pole connection mechanical optimization analysis and design system, comprising: a micro image acquisition module: constructing microstructure image dataset; a heterogeneous space mapping module: constructing anisotropic permeability tensor field, thereby mapping the mechanical transmission path of solid materials to a digitally defined virtual fluid permeation field; a parallel flow field solving module: calculating to obtain a stress flux distribution thermodynamic map representing stress transmission strength; a topology inversion optimization module: reconstructing the geometric reinforcement parameters and material distribution of the connector based on the identification results; a design verification closed loop module: triggering a new round of simulation calculation until the preset convergence condition is met, outputting the final optimized design drawing; the present application solves the contradiction between microstructure complexity and engineering design timeliness, supports rapid verification and multiple optimization of design schemes.

The invention discloses a discrete element method for representing elastic-plastic fiber deformation, and belongs to the field of computational mechanics and numerical simulation. According to the method, to-be-analyzed elastic-plastic fibers are divided into ideal elastic-plastic fibers and damaged elastic-plastic fibers, for different types of elastic-plastic fibers, a bonded ball model is adopted to represent fiber dynamic characteristics, and a bending moment-bending angle analysis model of elastic-plastic bonds is established by integrating cross section stress of the bonds in the fibers; and further constructing a bending moment increment model at any bending angle, and simplifying by adopting a piecewise linear method to obtain a simplified model of the elastic-plastic key bending moment-bending angle relationship, so as to obtain the deformation characteristics of the to-be-analyzed elastic-plastic fiber. According to the discrete element model capable of representing the elastic-plastic bending deformation of the fibers, accurate description and piecewise linear simplification of nonlinear mechanical behaviors of the fiber materials are achieved, and the calculation efficiency is remarkably improved while the calculation precision is guaranteed.

The invention provides a composite material plate shell multi-scale simulation method based on clustering order reduction, and belongs to the technical field of composite material multi-scale computational mechanics, and the method comprises the steps: building a voxelization 3D model and an associated 2D plate shell discrete model; constructing a single Lippmann-Schwinger control equation on the basis of a plate-shell theory and a periodic boundary, solving a macroscopic equivalent resultant force / bending moment of the single Lippmann-Schwinger control equation by adopting FFT (Fast Fourier Transform), and grouping discrete points through a clustering algorithm to realize order reduction; the method comprises the following steps of: clustering a plurality of plate shell models, disassembling an original equation into a cluster-level discrete equation based on a clustering result, constructing a constitutive updating scheme based on a 3D model, solving the model based on a Newton-Raphson iteration scheme, and distributing a cluster result to pixels of the plate shell model according to a solving result to obtain a mechanical response of the model. According to the method, the theoretical model is established through physical reasonable hypothesis; effective unit condensation of a real 3D model is realized by constructing a voxelization-based finite prism discrete scheme; and the elastoplasticity and damage failure analysis of the composite material is realized by constructing a 2D plan-3D solid model alternate solution scheme.