CN119227469A - A structural mechanics reliability analysis method for container mobile laboratory - Google Patents
A structural mechanics reliability analysis method for container mobile laboratory Download PDFInfo
- Publication number
- CN119227469A CN119227469A CN202411750518.7A CN202411750518A CN119227469A CN 119227469 A CN119227469 A CN 119227469A CN 202411750518 A CN202411750518 A CN 202411750518A CN 119227469 A CN119227469 A CN 119227469A
- Authority
- CN
- China
- Prior art keywords
- container
- grid
- node
- shell
- unit
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
- G06T17/205—Re-meshing
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Computer Graphics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The invention belongs to the technical field of reliability analysis, and discloses a method for analyzing the structural mechanical reliability of a container mobile laboratory; the method comprises the steps of obtaining container load data and container basic data, establishing a container finite element model based on the container load data and the container basic data, obtaining a unit stress field and a deformation amplitude of a container based on the container finite element model, determining a failure mode of the container finite element model based on the unit stress field and the deformation amplitude, calculating reliability indexes of the failure mode, calculating the failure probability of the failure mode based on the failure mode and the reliability indexes corresponding to the failure mode, and sending the failure probability to a reliability analysis terminal, so that design quality, use safety and reliability of a container mobile laboratory are improved.
Description
Technical Field
The invention relates to the technical field of reliability analysis, in particular to a method for analyzing the structural mechanical reliability of a container mobile laboratory.
Background
The patent with the application publication number of CN115169020A discloses a truss structure system reliability analysis method based on an automatic updating model, which comprises the steps of establishing a truss finite element model of a truss structure to be analyzed, fitting to obtain a response surface equation according to rod stress output by the truss finite element model, establishing a function of each rod, establishing an optimization model according to the geometric meaning and the function of a reliability index, calculating the reliability index of each rod, obtaining a plurality of failure modes for representing failure of the truss structure according to the reliability index of each rod when determining that the reliability indexes of all the rods are converged by utilizing an iteration criterion, wherein each failure mode comprises a failure path, calculating the reliability index and the failure probability of each failure mode, and calculating the reliability index and the failure probability of the truss structure by utilizing PNET, so that the problem of an implicit function is solved, and the accuracy is higher.
The conventional analysis method in the field of container mobile laboratories still lacks a special analysis method for special structures, the conventional general method cannot fully consider special working conditions and constraint conditions of the containers, structural safety and use reliability of the containers are difficult to comprehensively evaluate, the conventional modeling method is too simplified, the actual state of the containers cannot be accurately described, the deviation between a model and the actual structure is caused, the analysis precision is not high, the conventional analysis only considers single or partial failure modes, potential failure risks of the containers under complex working conditions cannot be comprehensively evaluated, the analysis results lack comprehensiveness, meanwhile, uncertainty of parameters such as materials, geometry, loads and the like is ignored, the actual conditions cannot be truly reflected, the reliability evaluation results have deviation, in addition, the calculation efficiency is low, the period is long, the analysis results cannot be rapidly obtained, the engineering actual requirements are difficult to be met, and the accuracy and comprehensiveness of the analysis results are to be improved.
In view of the above, the present invention proposes a method for analyzing the mechanical reliability of a container mobile laboratory structure to solve the above-mentioned problems.
Disclosure of Invention
In order to overcome the defects in the prior art and achieve the purposes, the invention provides a method for analyzing the structural mechanics reliability of a container mobile laboratory, which comprises the following steps that S1, container load data and container base data are acquired, and a container finite element model is built based on the container load data and the container base data;
S2, acquiring a unit stress field and deformation amplitude of the container based on the finite element model of the container;
S3, determining a failure mode of the container finite element model based on the unit stress field and the deformation amplitude, and calculating a reliability index of the failure mode;
and S4, calculating the failure probability of the failure mode based on the failure mode and the reliability index corresponding to the failure mode, and sending the failure probability to the reliability analysis terminal.
Further, the container load data comprises static load and dynamic load, the container base data comprises geometric dimension and material attribute data of the container, wherein the geometric dimension comprises length, width and height, and the material attribute data comprises wall thickness and material type;
the static load comprises instrument dead weight load data and container dead weight load data, and the dynamic load comprises instrument inertia load data, instrument vibration load data and instrument impact load data.
Further, the method for building the finite element model of the container comprises the following steps:
establishing a three-dimensional geometric model of the container according to the geometric dimension of the container, and refining the three-dimensional geometric model based on the wall thickness of the container to obtain a refined geometric model;
Defining a thinned geometric model into corresponding material properties according to the material types, wherein the material properties comprise elastic modulus, poisson's ratio and density;
The method comprises the steps of dividing a shell unit of a thinned geometric model to generate a finite element grid model, applying instrument dead weight load data to corresponding instrument installation positions based on the finite element grid model, giving the container dead weight load data to the finite element grid model, converting instrument inertia load data into corresponding acceleration or force to be applied to the instrument installation positions, converting instrument vibration load data into vibration excitation to be applied to the corresponding instrument installation positions, converting instrument impact load data into impact force or acceleration time course to be applied to the corresponding instrument installation positions, and completing the establishment of the container finite element model.
Further, the method for refining the three-dimensional geometric model comprises the following steps:
dividing the three-dimensional geometric model into n grid cells initially, and calculating the average thickness of each grid cell, dividing the n grid cells into n1 new grid cells based on the average thickness of each grid cell, and recording as n1 subdivision grid cells;
selecting M2 subdivision grid cells from n1 subdivision grid cells as seed cells, and constructing a priority queue for each seed cell, wherein elements in the priority queue are neighborhood cells of the seed cell;
Defining a growing area, adding the neighborhood units into the growing area in sequence according to the sequence in the priority queue of each seed unit to obtain M2 area grids, reconstructing a topological relation according to the M2 area grids to finally obtain M subareas, and defining thickness values of each subarea to obtain a refined geometric model;
the manner of defining the thickness value includes:
Defining a thickness distribution function ;For the three-dimensional coordinates of the nodes on the three-dimensional geometric model, for the grid units in each partition, the corresponding thickness values are obtained in the corresponding thickness distribution functions according to the barycentric coordinates of the grid units, and the definition of the thickness values is completed;
Thickness distribution function Wherein, the method comprises the steps of,Is a constant coefficient of the number of the pieces of the material,For the characteristic dimensions of the container,In order to be the modulus of elasticity of the material,Is the moment of inertia of the cross section,In order to distribute the load,Is a bending moment load.
Further, the calculating the average thickness of each grid cell includes:
Defining the average on a three-dimensional geometric model Each node, for each nodeCalculate its and adjacent nodeNormal vector change rate of (2)Wherein, the method comprises the steps of, wherein,Is a nodeIs characterized by a normal vector of (c),Is a nodePresetting a change rate threshold value when the normal vector changes rateWhen the change rate threshold is larger than the change rate threshold, defining the nodeIs positioned in the corner area;
For defining nodes located in corner regions Giving it a thickness weightA thickness weight of 1 is given to the nodes which are not in the corner area;
wherein, the method comprises the steps of, Is a nodeThe wall thickness gradient values at the points,AndIn order to amplify the coefficient of the power,Is a nodeThe proportion of the volume of the grid unit;
wherein, the method comprises the steps of, Is a nodeThe volume fraction of the material to be processed is related,Is the firstThe total volume of the individual grid cells,As the average volume of the grid cells,Is a regularization parameter;
the acquisition mode of (1) comprises the following steps:
Finding a containing node Is provided in the form of a grid of cells, for each containing nodeConstructing a new small polyhedron so that nodesIs a vertex of the small polyhedron, calculates the sum of volumes of all the small polyhedrons as;
Wherein, the method comprises the steps of,For the wall thickness at the node j,Is the third node in the neighborhood of the node jThe wall thickness at the individual nodes is such that,Is the firstThe basis functions of the individual nodes are chosen,As a basis functionAt the nodeA gradient is provided at the point of the gradient,Is the firstThe weight coefficient of the individual node is determined,Is a nodeThe total number of nodes in the neighborhood;
wherein, the method comprises the steps of, Is a nodeAnd (3) withIs used for the distance of (a),To influence the radius parameter, thenAverage thickness of individual grid cells。
Further, the dividing the n grid cells into n1 new grid cells includes:
acquiring grid cells Neighbor cell set of (a)Neighborhood unit setComprising and grid cellsAll adjacent grid cells sharing nodes or edges are recorded as neighborhood cells;
Computing grid cells Maximum thickness difference between adjacent cellsWherein, the method comprises the steps of,Is a grid cellIs used for the production of a plastic film,Is a grid cellNeighborhood of cells of (a)An average thickness of (2);
preset thickness difference threshold interval If (1)At the position ofIn, then grid cellsOne time average is subdivided into two grid cells, ifGreater thanThen grid cellOne time average is thinned into four grid cells;
And combining the original grid cells with the grid cells obtained by average refinement to obtain n1 new grid cells.
Further, the method for obtaining the unit stress field and the deformation amplitude comprises the following steps:
for each shell element Based on the unit shape function and the material attribute, calculating to obtain a corresponding unit stiffness matrixAssembling the unit rigidity matrix of all the shell units to obtain the total rigidity matrix of the whole container finite element modelCell stiffness matrixIs to make the displacement strain functionIntegrating to obtain; Wherein, the method comprises the steps of, wherein, ,Is an elastic matrix; Is that Is to be used in the present invention,Is the derivative of the shape function,Is a coordinate jacobian matrix;
for each shell element Calculating corresponding unit quality matrixAnd assembled into a total mass matrix of the whole container finite element modelBased on the total mass matrixAnd node load, constructing and obtaining a total load matrix;
Solving linear equationsObtaining a node displacement matrixFor each shell elementAccording to a unit-shaped functionAnd node displacement matrixCalculating to obtain corresponding shell unit displacement fieldFor the displacement field of the shell unitPerforming differentiation operation to obtain corresponding shell element strain fieldShell unitCorresponding unit stress fieldThe deformation amplitude is the absolute value of the node displacement of the shell element in the node displacement matrix.
Further, the determining manner of the failure mode includes:
based on the cell stress field, calculate each shell cell Equivalent stress of (2);
Wherein, the method comprises the steps of,、AndThe components of normal stress of the shell element in three dimensions,、AndThe components of the shear stress of the shell units in the three-dimensional direction are respectively;
Presetting an equivalent stress threshold For each shell elementIf (if)Greater thanThe shell unit is thenAdding stress concentration area unit sets;
For stress concentration area unit setEach shell element inCalculating the yield allowance;
Wherein, the method comprises the steps of,Is a characteristic value of the material yield strength of the shell element,For the yield strength uncertainty correction factor,Is the material yield criterion index;
If it is Then define the shell elementIs a failure mode of yielding;
for stress concentration area unit set Each shell element inCalculating the fracture allowance;
Wherein, the method comprises the steps of,Correcting the coefficient for the ultimate strength uncertainty; For the strain triaxial correction coefficient, Is a material fracture criterion index; is a characteristic value of the material ultimate strength of the shell element;
Three-axis degree correction coefficient of strain Wherein, the method comprises the steps of,Is equivalent to Poisson's ratioThen define the shell elementFailure mode for fracture;
presetting a deformation threshold For all the shell unitsIf (if)Greater thanThe shell unit is thenAdding a set of deformation units;
To deform the unit setDividing to obtain K1 grid blocks, and dividing each grid blockAnalyzing the characteristic value to obtain a block stiffness matrix KP, and solving the characteristic value of KP;
If it isV is smaller than 0, then the grid blockBuckling, if the shell unitIf the grid blocks with buckling existing in the shell unit are larger than a preset buckling quantity threshold value, the shell unitIs a failure mode of buckling;
acquiring S-N curves of the material for each shell element Calculating the fatigue life cycle number according to the stress-strain cycle amplitude and the S-N curve;
If all areLess than a predetermined number of life cyclesAdding a fatigue crack initiation set PUY;
For each shell element in the PUY Establishing a crack propagation model, wherein the expression of the crack propagation model is as follows: wherein, the method comprises the steps of, For the size of the crack to be a size,For the number of cycles to be counted,AndIs a material constant which is a function of the material,Is the stress intensity factor range; Representing the increment of crack size per cycle;
Based on For crack sizePerforming cycle extension analysis whenReaching a preset critical dimensionWhen judging the corresponding shell unitPresetting an invalidation threshold for invalidation, if the number of circulation times when invalidation is judged to be smaller than the invalidation threshold, corresponding shell unitIs a failure mode of fatigue cracking;
and summarizing all the obtained failure modes to obtain the failure mode of the container finite element model.
Further, the calculation mode of the reliability index includes:
for each failure mode Establishing a corresponding limit state equationWherein, the method comprises the steps of,Is in failure modeIs characterized by a structural resistance of (a),Is in failure modeIs a structural action effect; is a substantially random variable vector;
For each substantially random variable Determining the statistical distribution type and calculating the statistical distribution parameters;
for each failure mode Is the limit state equation of (2)Calculating according to the statistical distribution parameters of the basic random variable XCorresponding mean valueAnd standard deviation;
Failure modeCorresponding reliability indexWherein, the method comprises the steps of,The coefficients are modified for the reliability index,And correcting the coefficient for the uncertainty of the reliability index.
Further, the obtaining manner of the failure probability includes:
Will fail mode Corresponding reliability indexStandardized to obtain standard normal distribution variableCalculating standard normal distribution variableCorresponding probability valueProbability valueI.e. failure modeCorresponding failure probability.
The invention relates to a technical effect and advantages of a container mobile laboratory structural mechanics reliability analysis method, which comprises the following steps:
The method comprehensively, accurately and reasonably analyzes the structural safety and the use reliability of the container mobile laboratory, comprehensively considers static and dynamic loads by establishing an accurate finite element model, greatly improves analysis precision by accurately simulating the change of the wall thickness of the container, particularly in a stress concentration area, comprehensively evaluates potential failure risks of the container by integrating various failure modes, enables analysis results to be more in accordance with actual conditions by introducing a reliability index calculation method based on a limit state equation, simultaneously can quickly obtain analysis results, remarkably improves analysis efficiency, further considers factors such as spatial distribution of material properties, stress concentration effect, fatigue crack propagation and the like, and further improves the comprehensive and accuracy of analysis.
Drawings
FIG. 1 is a schematic diagram of a method for analyzing the mechanical reliability of a container mobile laboratory structure;
Fig. 2 is a schematic diagram of a system for analyzing the mechanical reliability of a container mobile laboratory structure according to the present invention.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
Example 1
Referring to fig. 1, the method for analyzing the mechanical reliability of a container mobile laboratory structure according to the present embodiment includes:
s1, acquiring container load data and container base data, and establishing a container finite element model based on the container load data and the container base data;
S2, acquiring a unit stress field and deformation amplitude of the container based on the finite element model of the container;
S3, determining a failure mode of the container finite element model based on the unit stress field and the deformation amplitude, and calculating a reliability index of the failure mode;
and S4, calculating the failure probability of the failure mode based on the failure mode and the reliability index corresponding to the failure mode, and sending the failure probability to the reliability analysis terminal.
Further, the container load data comprises static load and dynamic load, the container base data comprises geometric dimension and material attribute data of the container, the geometric dimension comprises length, width and height, the material attribute data comprises wall thickness and material type, and the collection points of the wall thickness are enough at different positions of the wall plate;
the static load comprises instrument dead weight load data and container dead weight load data;
the dynamic load comprises instrument inertia load data, instrument vibration load data and instrument impact load data;
The instrument dead weight load data and the container dead weight load data refer to the gravity load of the instrument and the equipment installed in the container and the gravity load of the container.
The instrument inertial load data are inertial force load data generated by instrument equipment in the container during operation, and the instrument vibration load data are vibration load data generated by the instrument equipment in the container during operation, including load data caused by vibration frequency and vibration amplitude;
the instrument impact load data is impact load generated in the processes of starting, stopping and the like of instrument equipment in the container.
Further, the method for establishing the finite element model of the container comprises the following steps:
establishing a three-dimensional geometric model of the container according to the geometric dimension of the container, and refining the three-dimensional geometric model based on the wall thickness of the container to obtain a refined geometric model;
Defining the thinned geometric model into corresponding material attributes according to the material types, wherein the material attributes comprise elastic modulus, poisson ratio and density;
the method for refining the three-dimensional geometric model comprises the following steps:
initially dividing the three-dimensional geometric model into n grid cells, and calculating the average thickness of each grid cell;
Specifically, the three-dimensional geometric model is defined on average Each node, for each nodeCalculate its and adjacent nodeNormal vector change rate of (2)Wherein, the method comprises the steps of, wherein,Is a nodeIs characterized by a normal vector of (c),Is a nodePresetting a change rate threshold value when the normal vector changes rateWhen the change rate threshold is larger than the change rate threshold, defining the nodeIs positioned in the corner area;
For defining nodes located in corner regions Giving it a thickness weightA thickness weight of 1 is given to the nodes which are not in the corner area;
wherein, the method comprises the steps of, Is a nodeThe wall thickness gradient values at the points,AndIn order to amplify the coefficient of the power,Is a nodeProportion of the mesh cell volumeWherein, the method comprises the steps of,Is a nodeThe volume fraction of the material to be processed is related,Is the firstThe total volume of the individual grid cells (referring to the geometric volume of the grid cell itself),As the average volume of the grid cells,For regularization parameters, the value range is usually 0.01-0.1;
the acquisition mode of (1) comprises the following steps:
Finding a containing node Is provided in the form of a grid of cells, for each containing nodeA new small polyhedron (such as hexahedron) is constructed so that the nodesIs a vertex of the small polyhedron, calculates the sum of volumes of all the small polyhedrons as;
Wherein, the method comprises the steps of,For the wall thickness at the node j,Is the third node in the neighborhood of the node jThe wall thickness at the individual nodes is such that,Is the firstThe basis functions of the individual nodes (linear lagrangian basis functions are used),As a basis functionAt the nodeThe gradient at (calculated by means of numerical integration),Is the firstThe weight coefficient of the individual node is determined,Is a nodeTotal number of nodes within the neighborhood (domain obtained by predefined domain radius);
wherein, the method comprises the steps of, Is a nodeAnd (3) withIs used for the distance of (a),To influence the radius parameter (when calculating the wall thickness gradient at a node, wall thickness values at other nodes in the neighborhood of the node need to be considered, i.e. the radius range of the neighborhood is defined), thenAverage thickness of individual grid cells;
Dividing n grid cells into n1 new grid cells based on the average thickness of each grid cell, specifically, obtaining grid cellsNeighbor cell set of (a)Neighborhood unit setComprising and grid cellsAll adjacent grid cells sharing nodes or edges are recorded as neighborhood cells;
Computing grid cells Maximum thickness difference between adjacent cellsWherein, the method comprises the steps of,Is a grid cellIs used for the production of a plastic film,Is a grid cellNeighborhood of cells of (a)An average thickness of (2);
preset thickness difference threshold interval If (1)At the position ofIn, then consider the grid cellThe thickness of the region is greatly changed, and the grid cells are required to be arrangedOne time average is subdivided into two grid cells, ifGreater thanThen consider the grid cellThe thickness of the region is too varied, and the grid cells are required to be arrangedOne time average is thinned into four grid cells;
Combining the original grid cells with the grid cells obtained by average refinement to obtain n1 new grid cells, wherein the n1 new grid cells are marked as n1 subdivided grid cells;
Selecting M2 subdivision grid cells from n1 subdivision grid cells as seed cells, constructing a priority queue for each seed cell for determining the growth sequence of the seed cells, wherein elements in the priority queue are neighborhood cells of the seed cells, and similarity measurement (such as thickness difference, geometric distance and the like) is generally adopted as the basis of priority, wherein the higher the similarity is, the more preferentially the growth area is added;
defining a growing area, adding the neighborhood units into the growing area in sequence according to the sequence in the priority queue of each seed unit until a certain stopping condition is met (if the number of units in the area reaches a preset value, the similarity is lower than a threshold value and the like), obtaining M2 area grids, reconstructing a topological relation according to the M2 area grids, and finally obtaining M subareas;
Defining a thickness value for each partition to obtain a refined geometric model;
the manner in which the thickness values are defined includes:
Defining a thickness distribution function ;For the three-dimensional coordinates of the nodes on the three-dimensional geometric model, for the grid units in each partition, the corresponding thickness values are obtained in the corresponding thickness distribution functions according to the barycentric coordinates of the grid units, and the definition of the thickness values is completed;
Thickness distribution function Wherein, the method comprises the steps of,Is a constant coefficient of the number of the pieces of the material,For the characteristic dimensions of the container,In order to be the modulus of elasticity of the material,Is the moment of inertia of the section, is a function of the wall thickness,In order to distribute the load,Is a bending moment load;
It should be noted that the characteristic dimension of a container refers to the maximum value of the length, width or height of the container, the elastic modulus of a material represents the proportional relation between stress and strain of the material during elastic deformation, different materials (such as steel, aluminum alloy and the like) have different elastic modulus values, the section moment of inertia is a mechanical property parameter describing the geometric shape of a section, and for a container, the section moment of inertia is proportional to the third power of the height of the section;
the distributed load refers to the load size received by the unit length/area along the length, width and height directions of the container, the bending moment load refers to the bending moment load acting on the container, and the bending moment load can cause the container to bend and deform;
The method comprises the steps of dividing a thinned geometric model into shell units (dispersing the whole structural model into a plurality of shell units according to the form of the shell units) to generate a finite element grid model (comprising a plurality of shell units), applying instrument dead weight load data to corresponding instrument installation positions based on the finite element grid model, giving the container dead weight load data (which can be realized by defining gravity acceleration) to the finite element grid model, converting instrument inertia load data into corresponding acceleration or force to be applied to the instrument installation positions, converting instrument vibration load data into vibration excitation to be applied to the corresponding instrument installation positions, converting instrument impact load data into impact force or acceleration time course to be applied to the corresponding instrument installation positions, and completing the establishment of the container finite element model.
Further, the method for obtaining the unit stress field and the deformation amplitude comprises the following steps:
for each shell element Based on the unit shape function and the material attribute, calculating to obtain a corresponding unit stiffness matrixAssembling the unit rigidity matrix of all the shell units to obtain the total rigidity matrix of the whole container finite element modelCell stiffness matrixIs to make the displacement strain functionIntegrating to obtain; Wherein, the method comprises the steps of, wherein, ,Is an elastic matrix; Is that Is to be used in the present invention,Is the derivative of the form function (function of polynomial interpolation of the shell elements),Is a coordinate jacobian matrix;
for each shell element Calculating corresponding unit quality matrixAnd assembled into a total mass matrix of the whole container finite element model;
The mass matrix of the shell unit is obtained by integrating the material density and the volume of the shell unit according to the shape function of the shell unit, and finally, the mass matrix of the shell unit can be obtained by superposing and summing the mass matrices of all the shell units.
Based on a total mass matrixAnd node load (load applied to finite element model nodes), constructing to obtain a total load matrix;
Solving linear equationsObtaining a node displacement matrixThe node displacement matrix is a matrix formed by node displacement of the shell unit;
for each shell element According to a unit-shaped function(Also as a function of polynomial interpolation of the shell elements) and node displacement matricesCalculating to obtain corresponding shell unit displacement fieldFor the displacement field of the shell unitPerforming differentiation operation to obtain corresponding shell element strain fieldShell unitCorresponding unit stress field;
The deformation amplitude is the absolute value of the node displacement of the shell element in the node displacement matrix.
Further, the determining manner of the failure mode includes:
based on the cell stress field, calculate each shell cell Equivalent stress of (2);
Wherein, the method comprises the steps of,、AndThe components of normal stress of the shell element in three dimensions,、AndThe components of the shear stress of the shell element in the three-dimensional direction (all obtained based on the element stress field);
Presetting an equivalent stress threshold For each shell elementIf (if)Greater thanThe shell unit is thenAdding stress concentration area unit sets;
For stress concentration area unit setEach shell element inCalculating the yield allowance;
Wherein, the method comprises the steps of,Is a characteristic value of the material yield strength of the shell element,For the correction of the uncertainty of the yield strength, usually we take,For the material yield criterion index, generally 1.5 to 2.0 is taken;
If it is Then define the shell elementIs a failure mode of yielding;
for stress concentration area unit set Each shell element inCalculating the fracture allowance;
Wherein, the method comprises the steps of,A correction factor for the ultimate intensity uncertainty, typically less than 1; For the strain triaxial correction coefficient, Is a material fracture criterion index (the same range of values as the material yield criterion index); is a characteristic value of the material ultimate strength of the shell element;
Three-axis degree correction coefficient of strain Wherein, the method comprises the steps of,Is equivalent to Poisson's ratioThen define the shell elementFailure mode for fracture;
presetting a deformation threshold For all the shell unitsIf (if)Greater thanThe shell unit is thenAdding a set of deformation units;
To deform the unit setDividing to obtain K1 grid blocks, and dividing each grid blockAnalyzing the characteristic value, and calculating to obtain a block stiffness matrix KP;
specifically, an initial vector is selected, a triple diagonal matrix is constructed through matrix-vector multiplication, all eigenvalues and eigenvectors of the triple diagonal matrix are solved, and the eigenvalues of the original matrix are approximated by the eigenvectors;
solving the characteristic value of KP If (1)V is smaller than 0, then the grid blockBuckling, if the shell unitIf the grid blocks with buckling existing in the shell unit are larger than a preset buckling quantity threshold value, the shell unitIs a failure mode of buckling;
acquiring S-N curve (stress-cycle number curve) of the material for each shell element Calculating the fatigue life cycle number according to the stress-strain cycle amplitude and the S-N curveSpecifically, the periodic component and the amplitude of the stress-strain periodic amplitude are extracted through methods such as rain flow counting, and the corresponding fatigue life cycle number can be estimated by substituting the amplitude into an S-N curve.
The stress-strain cycle amplitude is obtained by calculating stress/strain time response of a unit under a load process according to finite element analysis, extracting cycle components from the response time by adopting a fractional counting method or a rain flow counting method and the like, and calculating to obtain the stress-strain cycle amplitude;
If it is Less than a predetermined number of life cyclesShell elementBelongs to the fatigue crack initiation dangerous area, allLess thanAdding a fatigue crack initiation set PUY;
For each shell element in the PUY Establishing a crack propagation model, wherein the expression of the crack propagation model is as follows: wherein, the method comprises the steps of, For the size of the crack to be a size,For the number of cycles (representing the number of cycles the structure undergoes under cyclic loading),AndIs a constant of materialReflecting the intrinsic crack growth rate of the material under specific circumstances,Reflecting the material pairIs sensitive to the degree of sensitivity of (c),Is the stress intensity factor range; Representing the increment of crack size per cycle;
Based on For crack sizes (finite element, boundary element, etc.)Performing cycle extension analysis whenReaching a preset critical dimensionWhen judging failure, presetting a failure threshold, and if judging that the cycle number of the failure is smaller than the failure threshold, the corresponding shell unitIs a failure mode of fatigue cracking;
and summarizing all the obtained failure modes to obtain the failure mode of the container finite element model.
Further, the calculation method of the reliability index includes:
for each failure mode determined in the container finite element model Establishing a corresponding limit state equationWherein, the method comprises the steps of,Is in failure modeIs characterized by a structural resistance of (a),Is in failure modeIs a structural action effect; Is a substantially random variable vector (including random variables such as load, material properties, geometry, etc.);
For each substantially random variable Determining the statistical distribution type (such as normal distribution, lognormal distribution and the like) and calculating the statistical distribution parameters (such as mean value, standard deviation and the like);
for each failure mode Is the limit state equation of (2)Calculating according to the statistical distribution parameters of the basic random variable XCorresponding mean valueAnd standard deviation;
Failure modeCorresponding reliability indexWherein, the method comprises the steps of,For the reliability index correction factor, usually 0.1-0.2,For the uncertainty correction coefficient of the reliability index, usually 0.1-0.3 is adopted; The larger the probability of occurrence of the corresponding failure mode is, the smaller the reliability is.
Further, the obtaining method of the failure probability includes:
Will fail mode Corresponding reliability indexNormalization (linear transformation obeys standard normal distribution) to obtain standard normal distribution variableCalculating standard normal distribution variableCorresponding probability valueProbability valueI.e. failure modeCorresponding failure probability.
The method and the device for analyzing the structural safety and the use reliability of the container mobile laboratory comprehensively, accurately and reasonably analyze the structural safety and the use reliability of the container mobile laboratory, comprehensively consider static and dynamic loads by establishing an accurate finite element model, greatly improve analysis precision by accurately simulating the change of the wall thickness of the container, particularly in a stress concentration area, comprehensively evaluate potential failure risks of the container by integrating multiple failure modes, enable an analysis result to be more in line with actual conditions by introducing a reliability index calculation method based on a limit state equation, and simultaneously quickly obtain the analysis result, remarkably improve analysis efficiency, further improve the comprehensiveness and the accuracy of analysis by considering factors such as spatial distribution of material properties, stress concentration effect, fatigue crack expansion and the like, and improve the design quality, the use safety and the reliability of the container mobile laboratory in general.
Example 2
Referring to fig. 2, the detailed description of the embodiment is not provided in the description of embodiment 1, and a system for analyzing the mechanical reliability of a container mobile laboratory structure is provided, which includes:
The model building module is used for acquiring container load data and container base data and building a container finite element model based on the container load data and the container base data;
the stress deformation fitting module is used for acquiring a unit stress field and deformation amplitude of the container based on the container finite element model;
The index output module is used for determining the failure mode of the container finite element model based on the unit stress field and the deformation amplitude and calculating the reliability index of the failure mode;
The failure judgment module calculates the failure probability of the failure mode based on the failure mode and the corresponding reliability index, and sends the failure probability to the reliability analysis terminal, and the modules are connected in a wired and/or wireless mode to realize data transmission among the modules.
Example 3
The embodiment discloses an electronic device, which comprises a memory, a processor and a computer program stored on the memory and capable of running on the processor, wherein the running mode of the method for analyzing the mechanical reliability of the container mobile laboratory structure is realized when the processor executes the computer program.
Since the electronic device described in this embodiment is an electronic device used to implement the method for analyzing mechanical reliability of a container mobile laboratory structure in this embodiment of the present application, based on the method for analyzing mechanical reliability of a container mobile laboratory structure described in this embodiment of the present application, a person skilled in the art can understand a specific implementation manner of the electronic device in this embodiment and various modifications thereof, so how to implement the method in this embodiment of the present application in this electronic device will not be described in detail herein. As long as the person skilled in the art implements the electronic device adopted by the method for analyzing the mechanical reliability of the container mobile laboratory structure in the embodiment of the application, the electronic device belongs to the scope of protection required by the application.
The above formulas are all formulas with dimensionality removed and numerical calculation, the formulas are formulas with the latest real situation obtained by software simulation through collecting a large amount of data, and preset parameters and threshold selection in the formulas are set by those skilled in the art according to the actual situation.
The above description is only a preferred embodiment of the present invention, and the protection scope of the present invention is not limited to the above examples, and all technical solutions belonging to the concept of the present invention belong to the protection scope of the present invention. It should be noted that modifications and adaptations to those skilled in the art without departing from the principles of the present invention are intended to be comprehended within the scope of the present invention.
Claims (10)
1. A method for analyzing the mechanical reliability of a container mobile laboratory structure, which is characterized by comprising the following steps:
s1, acquiring container load data and container base data, and establishing a container finite element model based on the container load data and the container base data;
S2, acquiring a unit stress field and deformation amplitude of the container based on the finite element model of the container;
S3, determining a failure mode of the container finite element model based on the unit stress field and the deformation amplitude, and calculating a reliability index of the failure mode;
and S4, calculating the failure probability of the failure mode based on the failure mode and the reliability index corresponding to the failure mode, and sending the failure probability to the reliability analysis terminal.
2. The method for analyzing the mechanical reliability of the structure of the container mobile laboratory according to claim 1, wherein the container load data comprises static load and dynamic load, the container base data comprises geometric dimension and material attribute data of the container, the geometric dimension comprises length, width and height, the material attribute data comprises wall thickness and material type;
the static load comprises instrument dead weight load data and container dead weight load data, and the dynamic load comprises instrument inertia load data, instrument vibration load data and instrument impact load data.
3. The method for analyzing the mechanical reliability of the container mobile laboratory structure according to claim 2, wherein the method for establishing the finite element model of the container comprises the following steps:
establishing a three-dimensional geometric model of the container according to the geometric dimension of the container, and refining the three-dimensional geometric model based on the wall thickness of the container to obtain a refined geometric model;
Defining a thinned geometric model into corresponding material properties according to the material types, wherein the material properties comprise elastic modulus, poisson's ratio and density;
The method comprises the steps of dividing a shell unit of a thinned geometric model to generate a finite element grid model, applying instrument dead weight load data to corresponding instrument installation positions based on the finite element grid model, giving the container dead weight load data to the finite element grid model, converting instrument inertia load data into corresponding acceleration or force to be applied to the instrument installation positions, converting instrument vibration load data into vibration excitation to be applied to the corresponding instrument installation positions, converting instrument impact load data into impact force or acceleration time course to be applied to the corresponding instrument installation positions, and completing the establishment of the container finite element model.
4. A method of analyzing mechanical reliability of a container mobile laboratory structure according to claim 3, wherein the method of refining the three-dimensional geometric model comprises:
dividing the three-dimensional geometric model into n grid cells initially, and calculating the average thickness of each grid cell, dividing the n grid cells into n1 new grid cells based on the average thickness of each grid cell, and recording as n1 subdivision grid cells;
selecting M2 subdivision grid cells from n1 subdivision grid cells as seed cells, and constructing a priority queue for each seed cell, wherein elements in the priority queue are neighborhood cells of the seed cell;
Defining a growing area, adding the neighborhood units into the growing area in sequence according to the sequence in the priority queue of each seed unit to obtain M2 area grids, reconstructing a topological relation according to the M2 area grids to finally obtain M subareas, and defining thickness values of each subarea to obtain a refined geometric model;
the manner of defining the thickness value includes:
Defining a thickness distribution function ;For the three-dimensional coordinates of the nodes on the three-dimensional geometric model, for the grid units in each partition, the corresponding thickness values are obtained in the corresponding thickness distribution functions according to the barycentric coordinates of the grid units, and the definition of the thickness values is completed;
Thickness distribution function Wherein, the method comprises the steps of,Is a constant coefficient of the number of the pieces of the material,For the characteristic dimensions of the container,In order to be the modulus of elasticity of the material,Is the moment of inertia of the cross section,In order to distribute the load,Is a bending moment load.
5. The method for analyzing mechanical reliability of a container mobile laboratory structure according to claim 4, wherein the calculating the average thickness of each grid cell comprises:
Defining the average on a three-dimensional geometric model Each node, for each nodeCalculate its and adjacent nodeNormal vector change rate of (2)Wherein, the method comprises the steps of, wherein,Is a nodeIs characterized by a normal vector of (c),Is a nodePresetting a change rate threshold value when the normal vector changes rateWhen the change rate threshold is larger than the change rate threshold, defining the nodeIs positioned in the corner area;
For defining nodes located in corner regions Giving it a thickness weightA thickness weight of 1 is given to the nodes which are not in the corner area;
wherein, the method comprises the steps of, Is a nodeThe wall thickness gradient values at the points,AndIn order to amplify the coefficient of the power,Is a nodeThe proportion of the volume of the grid unit;
wherein, the method comprises the steps of, Is a nodeThe volume fraction of the material to be processed is related,Is the firstThe total volume of the individual grid cells,As the average volume of the grid cells,Is a regularization parameter;
the acquisition mode of (1) comprises the following steps:
Finding a containing node Is provided in the form of a grid of cells, for each containing nodeConstructing a new small polyhedron so that nodesIs a vertex of the small polyhedron, calculates the sum of volumes of all the small polyhedrons as;
Wherein, the method comprises the steps of,For the wall thickness at the node j,Is the third node in the neighborhood of the node jThe wall thickness at the individual nodes is such that,Is the firstThe basis functions of the individual nodes are chosen,As a basis functionAt the nodeA gradient is provided at the point of the gradient,Is the firstThe weight coefficient of the individual node is determined,Is a nodeThe total number of nodes in the neighborhood;
wherein, the method comprises the steps of, Is a nodeAnd (3) withIs used for the distance of (a),To influence the radius parameter, thenAverage thickness of individual grid cells。
6. The method for analyzing the mechanical reliability of the container mobile laboratory structure according to claim 5, wherein the dividing the n grid cells into n1 new grid cells comprises:
acquiring grid cells Neighbor cell set of (a)Neighborhood unit setComprising and grid cellsAll adjacent grid cells sharing nodes or edges are recorded as neighborhood cells;
Computing grid cells Maximum thickness difference between adjacent cellsWherein, the method comprises the steps of,Is a grid cellIs used for the production of a plastic film,Is a grid cellNeighborhood of cells of (a)An average thickness of (2);
preset thickness difference threshold interval If (1)At the position ofIn, then grid cellsOne time average is subdivided into two grid cells, ifGreater thanThen grid cellOne time average is thinned into four grid cells;
And combining the original grid cells with the grid cells obtained by average refinement to obtain n1 new grid cells.
7. The method for analyzing mechanical reliability of a container mobile laboratory structure according to claim 6, wherein the means for acquiring the unit stress field and the deformation amplitude comprises:
for each shell element Based on the unit shape function and the material attribute, calculating to obtain a corresponding unit stiffness matrixAssembling the unit rigidity matrix of all the shell units to obtain the total rigidity matrix of the whole container finite element modelCell stiffness matrixIs to make the displacement strain functionIntegrating to obtain; Wherein, the method comprises the steps of, wherein, ,Is an elastic matrix; Is that Is to be used in the present invention,Is the derivative of the shape function,Is a coordinate jacobian matrix;
for each shell element Calculating corresponding unit quality matrixAnd assembled into a total mass matrix of the whole container finite element modelBased on the total mass matrixAnd node load, constructing and obtaining a total load matrix;
Solving linear equationsObtaining a node displacement matrixFor each shell elementAccording to a unit-shaped functionAnd node displacement matrixCalculating to obtain corresponding shell unit displacement fieldFor the displacement field of the shell unitPerforming differentiation operation to obtain corresponding shell element strain fieldShell unitCorresponding unit stress fieldThe deformation amplitude is the absolute value of the node displacement of the shell element in the node displacement matrix.
8. The method for analyzing mechanical reliability of a container mobile laboratory structure according to claim 7, wherein the determining manner of the failure mode comprises:
based on the cell stress field, calculate each shell cell Equivalent stress of (2);
Wherein, the method comprises the steps of,、AndThe components of normal stress of the shell element in three dimensions,、AndThe components of the shear stress of the shell units in the three-dimensional direction are respectively;
Presetting an equivalent stress threshold For each shell elementIf (if)Greater thanThe shell unit is thenAdding stress concentration area unit sets;
For stress concentration area unit setEach shell element inCalculating the yield allowance;
Wherein, the method comprises the steps of,Is a characteristic value of the material yield strength of the shell element,For the yield strength uncertainty correction factor,Is the material yield criterion index;
If it is Then define the shell elementIs a failure mode of yielding;
for stress concentration area unit set Each shell element inCalculating the fracture allowance;
Wherein, the method comprises the steps of,Correcting the coefficient for the ultimate strength uncertainty; For the strain triaxial correction coefficient, Is a material fracture criterion index; is a characteristic value of the material ultimate strength of the shell element;
Three-axis degree correction coefficient of strain Wherein, the method comprises the steps of,Is equivalent to Poisson's ratioThen define the shell elementFailure mode for fracture;
presetting a deformation threshold For all the shell unitsIf (if)Greater thanThe shell unit is thenAdding a set of deformation units;
To deform the unit setDividing to obtain K1 grid blocks, and dividing each grid blockAnalyzing the characteristic value to obtain a block stiffness matrix KP, and solving the characteristic value of KP;
If it isV is smaller than 0, then the grid blockBuckling, if the shell unitIf the grid blocks with buckling existing in the shell unit are larger than a preset buckling quantity threshold value, the shell unitIs a failure mode of buckling;
acquiring S-N curves of the material for each shell element Calculating the fatigue life cycle number according to the stress-strain cycle amplitude and the S-N curve;
If all areLess than a predetermined number of life cyclesAdding a fatigue crack initiation set PUY;
For each shell element in the PUY Establishing a crack propagation model, wherein the expression of the crack propagation model is as follows: wherein, the method comprises the steps of, For the size of the crack to be a size,For the number of cycles to be counted,AndIs a material constant which is a function of the material,Is the stress intensity factor range; Representing the increment of crack size per cycle;
Based on For crack sizePerforming cycle extension analysis whenReaching a preset critical dimensionWhen judging the corresponding shell unitPresetting an invalidation threshold for invalidation, if the number of circulation times when invalidation is judged to be smaller than the invalidation threshold, corresponding shell unitIs a failure mode of fatigue cracking;
and summarizing all the obtained failure modes to obtain the failure mode of the container finite element model.
9. The method for analyzing the mechanical reliability of the container mobile laboratory structure according to claim 8, wherein the calculation mode of the reliability index comprises the following steps:
for each failure mode Establishing a corresponding limit state equationWherein, the method comprises the steps of,Is in failure modeIs characterized by a structural resistance of (a),Is in failure modeIs a structural action effect; is a substantially random variable vector;
For each substantially random variable Determining the statistical distribution type and calculating the statistical distribution parameters;
for each failure mode Is the limit state equation of (2)Calculating according to the statistical distribution parameters of the basic random variable XCorresponding mean valueAnd standard deviation;
Failure modeCorresponding reliability indexWherein, the method comprises the steps of,The coefficients are modified for the reliability index,And correcting the coefficient for the uncertainty of the reliability index.
10. The method for analyzing mechanical reliability of a container mobile laboratory structure according to claim 9, wherein the obtaining manner of the failure probability comprises:
Will fail mode Corresponding reliability indexStandardized to obtain standard normal distribution variableCalculating standard normal distribution variableCorresponding probability valueProbability valueI.e. failure modeCorresponding failure probability.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202411750518.7A CN119227469B (en) | 2024-12-02 | 2024-12-02 | Structural mechanical reliability analysis method for container mobile laboratory |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202411750518.7A CN119227469B (en) | 2024-12-02 | 2024-12-02 | Structural mechanical reliability analysis method for container mobile laboratory |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN119227469A true CN119227469A (en) | 2024-12-31 |
| CN119227469B CN119227469B (en) | 2025-04-18 |
Family
ID=94040823
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202411750518.7A Active CN119227469B (en) | 2024-12-02 | 2024-12-02 | Structural mechanical reliability analysis method for container mobile laboratory |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN119227469B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN119962322A (en) * | 2025-04-09 | 2025-05-09 | 江西省水利科学院(江西省大坝安全管理中心、江西省水资源管理中心) | A safety analysis method for sluice structures under cyclic loading |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109362020A (en) * | 2018-09-30 | 2019-02-19 | 浙江中科电声研发中心 | A kind of Numerical Simulation Analysis method of speaker frame dynamic stiffness |
| US20190362041A1 (en) * | 2017-10-31 | 2019-11-28 | China University Of Mining And Technology | Reliability robust design method for multiple failure modes of ultra-d eep well hoisting container |
| CN110705163A (en) * | 2019-09-30 | 2020-01-17 | 北京航空航天大学 | A Reliability Analysis Method for Fatigue System of Composite Laminated Structure |
| CN111783327A (en) * | 2020-05-28 | 2020-10-16 | 西北工业大学 | Reliability analysis method of turbine blade under multi-field load based on support vector machine |
| CN115169020A (en) * | 2022-05-25 | 2022-10-11 | 石家庄铁道大学 | Reliability analysis method of truss structure system based on automatic update model |
| CN118886242A (en) * | 2024-06-26 | 2024-11-01 | 浙江工业大学 | An analytical method for calibrating the ultimate load of cracked structures under complex loads |
| CN119004602A (en) * | 2024-08-05 | 2024-11-22 | 重庆电讯职业学院 | Safety risk assessment method and system for assembled building |
-
2024
- 2024-12-02 CN CN202411750518.7A patent/CN119227469B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20190362041A1 (en) * | 2017-10-31 | 2019-11-28 | China University Of Mining And Technology | Reliability robust design method for multiple failure modes of ultra-d eep well hoisting container |
| CN109362020A (en) * | 2018-09-30 | 2019-02-19 | 浙江中科电声研发中心 | A kind of Numerical Simulation Analysis method of speaker frame dynamic stiffness |
| CN110705163A (en) * | 2019-09-30 | 2020-01-17 | 北京航空航天大学 | A Reliability Analysis Method for Fatigue System of Composite Laminated Structure |
| CN111783327A (en) * | 2020-05-28 | 2020-10-16 | 西北工业大学 | Reliability analysis method of turbine blade under multi-field load based on support vector machine |
| CN115169020A (en) * | 2022-05-25 | 2022-10-11 | 石家庄铁道大学 | Reliability analysis method of truss structure system based on automatic update model |
| CN118886242A (en) * | 2024-06-26 | 2024-11-01 | 浙江工业大学 | An analytical method for calibrating the ultimate load of cracked structures under complex loads |
| CN119004602A (en) * | 2024-08-05 | 2024-11-22 | 重庆电讯职业学院 | Safety risk assessment method and system for assembled building |
Non-Patent Citations (1)
| Title |
|---|
| 韩向阳;赵继广;段永胜;辛腾达;: "推进剂贮箱结构失效概率研究", 计算机仿真, no. 11, 15 November 2019 (2019-11-15) * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN119962322A (en) * | 2025-04-09 | 2025-05-09 | 江西省水利科学院(江西省大坝安全管理中心、江西省水资源管理中心) | A safety analysis method for sluice structures under cyclic loading |
Also Published As
| Publication number | Publication date |
|---|---|
| CN119227469B (en) | 2025-04-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Chowdhury et al. | High‐dimensional model representation for structural reliability analysis | |
| CN119227469B (en) | Structural mechanical reliability analysis method for container mobile laboratory | |
| EP2525296B1 (en) | Three-dimensional fluid simulation method | |
| CN119647209B (en) | A probe card simulation method and system based on reaction force analysis | |
| CN113987889A (en) | Method for fitting and predicting service life of welding structure ultrahigh cycle fatigue main S-N curve | |
| CN118586214B (en) | Furniture product tipping simulation method, system and equipment based on statics iteration | |
| Lv et al. | Seismic response characteristics of nuclear island structure at generic soil and rock sites | |
| CN116227045B (en) | Local stress strain field construction method and system for structural test piece | |
| CN118797841A (en) | A method, medium and system for determining the deformation degree of a large steel structure corridor after being suspended | |
| US6731996B1 (en) | Method of using isotropic plasticity to model the forming of anisotropic sheet | |
| CN118886599A (en) | Ancient building risk assessment method and device | |
| CN120257737A (en) | Parametric rapid modeling method, system, equipment and medium for foundation pit finite element model | |
| CN119047277B (en) | Submarine tunnel deformation prediction method under earthquake action | |
| CN119358084A (en) | Saturated porous media dynamic analysis method, device, storage medium and electronic equipment | |
| CN116484483B (en) | A reliability evaluation method for vibration isolation performance of nonlinear vibration isolators | |
| CN105404740B (en) | A kind of amendment sensitivity Structural Damage Identification based on Joint Index and limited information | |
| CN118690595A (en) | A crack growth analysis method based on sub-model method | |
| CN116910862A (en) | Method for rapidly and transitionally identifying damage degree of post-earthquake concrete arch dam structure | |
| CN113792461B (en) | A composite time-domain analysis method for dynamic response of engineering structures under extreme loads | |
| CN115859715A (en) | Method for calculating equivalent performance of periodic porous structure material | |
| CN111415346B (en) | A Gridding-based Damage Characterization Method for Particle Models | |
| Li et al. | Wind-resistant analysis for large-span single-layer reticulated shells | |
| CN114925580A (en) | A vibration assessment method for the foundation of a large-scale shaking table | |
| CN113722859B (en) | Uncertainty structure static response determination method based on convex polyhedron model | |
| CN119845469B (en) | Bridge construction structure stress real-time monitoring method and system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |