CN117313471A - Door machine stress field inversion method and system and electronic equipment - Google Patents
Door machine stress field inversion method and system and electronic equipment Download PDFInfo
- Publication number
- CN117313471A CN117313471A CN202311272894.5A CN202311272894A CN117313471A CN 117313471 A CN117313471 A CN 117313471A CN 202311272894 A CN202311272894 A CN 202311272894A CN 117313471 A CN117313471 A CN 117313471A
- Authority
- CN
- China
- Prior art keywords
- stress
- gantry crane
- value
- initial load
- door machine
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 60
- 238000004458 analytical method Methods 0.000 claims abstract description 57
- 238000006073 displacement reaction Methods 0.000 claims abstract description 30
- 210000004027 cell Anatomy 0.000 claims description 35
- 239000011159 matrix material Substances 0.000 claims description 19
- 230000006870 function Effects 0.000 claims description 13
- 238000004088 simulation Methods 0.000 claims description 12
- 239000007787 solid Substances 0.000 claims description 9
- 210000003888 boundary cell Anatomy 0.000 claims description 7
- 238000001514 detection method Methods 0.000 claims description 6
- 238000005259 measurement Methods 0.000 claims description 6
- 239000011800 void material Substances 0.000 claims description 3
- 230000008569 process Effects 0.000 abstract description 14
- 238000012544 monitoring process Methods 0.000 abstract description 8
- 238000004364 calculation method Methods 0.000 description 13
- 238000010586 diagram Methods 0.000 description 11
- 238000004590 computer program Methods 0.000 description 6
- 230000036541 health Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000007689 inspection Methods 0.000 description 3
- 230000010354 integration Effects 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000012916 structural analysis Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000013506 data mapping Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 230000005489 elastic deformation Effects 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR 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; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR 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)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The invention provides a gantry crane stress field inversion method, a gantry crane stress field inversion system and electronic equipment, which belong to the field of gantry crane monitoring, and the method comprises the following steps: acquiring a three-dimensional model of a door machine, and determining the displacement freedom degree and initial load of the door machine; taking the displacement freedom degree as a constraint condition, taking an initial load as a boundary condition, carrying out finite element analysis on stress distribution on a gantry crane based on a gantry crane three-dimensional model, and determining all stress extreme points on the gantry crane; the stress extreme points include: stress maximum point and stress minimum point; and obtaining a stress value at a stress extreme point on the door machine, accumulating the measured stress value and a stress initial value determined based on an initial load to obtain an actual stress value at the stress extreme point, then carrying out interpolation inversion on the actual stress value based on a door machine three-dimensional model by adopting an equal geometry analysis method based on NURBS, and simulating to obtain the stress field distribution of the door machine. The invention realizes the comprehensiveness and real-time performance of stress monitoring in the running process of the door machine.
Description
Technical Field
The invention belongs to the field of door machine monitoring, and particularly relates to a door machine stress field inversion method, a door machine stress field inversion system and electronic equipment.
Background
Hydropower station gantry cranes (gantry cranes) are hoisting equipment for opening and closing various gates of hydropower stations, and in the current daily gantry crane running state monitoring and health state inspection process, various defects exist, such as: the monitoring of the states of the structure, the mechanism and the components in the running process of the crane is not comprehensive and sufficient; through manual regular inspection, the efficiency is low, the inspection can not comprehensively grasp the running state of the equipment in real time, and the health state of each part of the equipment can not be comprehensively and real-time estimated.
Disclosure of Invention
Aiming at the defects of the prior art, the invention aims to provide a stress field inversion method, a stress field inversion system and electronic equipment of a gantry crane, and aims to solve the problem that the state of health of each part of the gantry crane cannot be comprehensively and real-timely estimated in the prior art.
To achieve the above object, in a first aspect, the present invention provides a gantry crane stress field inversion method, including:
acquiring a three-dimensional model of a door machine, and determining the displacement freedom degree and initial load of the door machine;
taking the displacement freedom degree as a constraint condition, taking an initial load as a boundary condition, carrying out finite element analysis on stress distribution on a gantry crane based on the gantry crane three-dimensional model, and determining all stress extreme points on the gantry crane; the stress extreme points include: stress maximum point and stress minimum point; the initial load is used for determining initial stress field distribution on the door machine;
and then interpolating and inverting the actual stress value by adopting an isogeometric analysis method based on Non-uniform rational B-Splines (NURBS) based on a gantry crane three-dimensional model to simulate to obtain the stress field distribution of the gantry crane.
Specifically, after the initial load is determined, the initial load and the initial stress value at the extreme point are reversely calculated by comparing the deformation value of the finite element model at the extreme point of the stress of the gantry crane under different loads with the actually measured deformation value.
In addition, the stress value at the stress extreme point can be measured by the stress sensor.
After the stress field distribution of the gantry crane is obtained through simulation, the health state of each part of the gantry crane, such as data of strain, vibration or displacement, can be further estimated by combining the stress field distribution.
In one possible implementation, the initial load is determined by:
presetting an initial load value, carrying out finite element analysis on stress distribution on a gantry crane based on the preset initial load value, the displacement freedom degree and a gantry crane three-dimensional model, and preliminarily determining a stress extreme point on the gantry crane;
acquiring an actual deformation value after a preset load is added at a preliminarily determined stress extreme point on the door machine, comparing the actual deformation value with a deformation value obtained after the preset load is added by finite element analysis simulation, and if the deformation values are inconsistent, adjusting the preset initial load value until the deformation values are consistent;
and taking the initial load value corresponding to the consistent deformation value as the final determined initial load.
And then, based on the initial load, calculating to obtain initial distribution of stress fields on the gantry crane in a finite element model, namely, determining stress values of all points.
The stress value of each extreme point under the initial load with the same deformation value is the initial stress value of each extreme point. In the elastic deformation section of the structure, the actual stress value of each extreme point=the measured value of the stress sensor+the initial value.
In one possible implementation, the displacement degree of freedom is a displacement constraint of the gantry crane, and is determined according to parameters of the gantry crane.
In one possible implementation, the interpolation inversion is performed on the measured stress value based on the NURBS-based isogeometric analysis method adopted by the gantry crane three-dimensional model, including:
establishing a geometric model by adopting NURBS basis functions based on a gantry crane three-dimensional model;
dispersing a physical field corresponding to an actual stress value at a stress measurement point based on a NURBS basis function in the established geometric model so as to perform isogeometric analysis, and performing interpolation inversion to obtain the stress of each point on the gantry crane three-dimensional model; the discrete equilibrium equation is determined by the stiffness matrix, gantry crane displacement vector, and external load.
In one possible implementation, NURBS surface interpolation is performed using NURBS basis functions of second order or more.
In one possible implementation, the geometric model includes: a hole unit, a entity unit and a boundary unit;
for void cells, the stiffness matrix K e The method comprises the following steps:
K e ＝K solid *ρ e
wherein K is solid Representing the stiffness matrix, ρ, of the physical unit e Is the rigidity ratio of the cavity unit to the entity unit.
In one possible implementation, when each cell in the geometric model is identified, at least one additional detection point is introduced into the entity cells and/or the cavity cells near the boundary cells, and if the domain in which the at least one additional detection point is located is different from the domain in which the cell is located, the corresponding cell is re-divided into boundary cells.
In a second aspect, the present invention provides a gantry crane stress field inversion system comprising:
the gantry crane model modeling unit is used for establishing a three-dimensional model of the gantry crane and determining the displacement freedom degree and initial load of the gantry crane;
the extreme point analysis unit is used for taking the displacement freedom degree as a constraint condition, taking an initial load as a boundary condition, carrying out finite element analysis on stress distribution on the gantry crane based on the gantry crane three-dimensional model, and determining all stress extreme points on the gantry crane; the stress extreme points include: stress maximum point and stress minimum point; the initial load is used for determining initial stress field distribution on the door machine;
and the stress field inversion unit is used for acquiring the stress value at the stress extreme point on the measured gantry crane, accumulating the measured stress value and the stress initial value determined based on the initial load to obtain the actual stress value at the stress extreme point, and then carrying out interpolation inversion on the actual stress value based on the gantry crane three-dimensional model by adopting an isogeometric analysis method based on non-uniform rational B-spline NURBS, so as to obtain the stress field distribution of the gantry crane through simulation.
In one possible implementation manner, the gantry crane model obtaining unit presets an initial load value, performs finite element analysis on stress distribution on the gantry crane based on the preset initial load value, the displacement freedom degree and the gantry crane three-dimensional model, and preliminarily determines a stress extreme point on the gantry crane; acquiring an actual deformation value after a preset load is added at a preliminarily determined stress extreme point on the door machine, comparing the actual deformation value with a deformation value obtained after the preset load is added by finite element analysis simulation, and if the deformation values are inconsistent, adjusting the preset initial load value until the deformation values are consistent; and taking the initial load value corresponding to the consistent deformation value as the final determined initial load, and then calculating to obtain the initial stress distribution of the gantry crane based on the initial load, so as to determine the stress extreme point on the gantry crane based on the simulation of the initial stress field distribution during finite element analysis.
In a third aspect, the present invention provides an electronic device comprising: at least one memory for storing a program; at least one processor for executing a memory-stored program, which when executed is adapted to carry out the method described in the first aspect or any one of the possible implementations of the first aspect.
In a fourth aspect, the present invention provides a computer readable storage medium storing a computer program which, when run on a processor, causes the processor to perform the method described in the first aspect or any one of the possible implementations of the first aspect.
In a fifth aspect, the invention provides a computer program product which, when run on a processor, causes the processor to perform the method described in the first aspect or any one of the possible implementations of the first aspect.
In general, the above technical solutions conceived by the present invention have the following beneficial effects compared with the prior art:
the invention provides a gantry crane stress field inversion method, a gantry crane stress field inversion system and electronic equipment, which take a gantry crane as a research object entity, establish a high-precision three-dimensional simulation physical field model of the gantry crane, establish a data mapping relation between the physical entity and a virtual model on the basis of key positions in the actual working process of the gantry crane and real-time monitoring data on the mechanism, and establish a digital twin body of the gantry crane. And determining stress extreme points on the gantry crane through finite element simulation, measuring the stress at the extreme points, performing curved surface interpolation inversion on the measured stress through a NURBS isogeometric analysis method, determining stress field distribution on the gantry crane, realizing comprehensiveness and instantaneity of stress monitoring in the running process of the gantry crane, further analyzing physical field states such as strain, vibration and the like of the gantry crane by combining on-line monitoring results, and realizing visualization of the running state of equipment.
The invention provides a gantry crane stress field inversion method, a gantry crane stress field inversion system and electronic equipment, wherein the inversion method based on the three-dimensional modeling stress field only needs to sample work on an engineering site, so that the work load on the site is effectively reduced, the sampling of a plurality of measuring points can be realized, and the fine measurement of a plurality of measuring points of stress is realized; the method can establish a three-dimensional numerical model with engineering large scale and single or complex structure, and invert the stress distribution characteristics of an inversion research area through parallel calculation; visual and visualized overall stress distribution conditions can be obtained, and information such as vertical stress, maximum horizontal main stress, minimum horizontal main stress values and directions of different levels, different inclined planes and different structural positions of the model can be effectively extracted and checked.
Drawings
FIG. 1 is a flow chart of a gantry crane stress field inversion method provided by an embodiment of the invention;
FIG. 2 is a schematic diagram of the construction of a three-dimensional model of a gantry crane provided by an embodiment of the present invention;
FIG. 3 is a stress distribution diagram of deformation ratio 1 in finite element analysis of a gantry crane three-dimensional model provided by the embodiment of the invention;
FIG. 4 is a stress distribution diagram of 1000 ten thousand deformation proportion in finite element analysis of a gantry crane three-dimensional model provided by the embodiment of the invention;
FIG. 5 is a graph of the spatial discrete contrast of finite element analysis and NURBS unit analysis provided by an embodiment of the present invention;
FIG. 6 is a schematic diagram of an isogeometric analysis and extended finite element technique combined with an embodiment of the present invention;
FIG. 7 is a schematic diagram of a partial cell re-identification scheme for a false identification case according to an embodiment of the present invention;
FIG. 8 is a flow chart of a gantry crane stress field inversion algorithm provided by an embodiment of the invention;
fig. 9 is a diagram of a gantry crane stress field inversion system according to an embodiment of the present invention.
Detailed Description
The present invention will be described in further detail with reference to the drawings and examples, in order to make the objects, technical solutions and advantages of the present invention more apparent. It should be understood that the specific embodiments described herein are for purposes of illustration only and are not intended to limit the scope of the invention.
The term "and/or" herein is an association relationship describing an associated object, and means that there may be three relationships, for example, a and/or B may mean: a exists alone, A and B exist together, and B exists alone. The symbol "/" herein indicates that the associated object is or is a relationship, e.g., A/B indicates A or B.
In embodiments of the invention, words such as "exemplary" or "such as" are used to mean serving as an example, instance, or illustration. Any embodiment or design described herein as "exemplary" or "e.g." in an embodiment should not be taken as preferred or advantageous over other embodiments or designs. Rather, the use of words such as "exemplary" or "such as" is intended to present related concepts in a concrete fashion.
Next, the technical scheme provided in the embodiment of the present invention is described.
FIG. 1 is a flow chart of a gantry crane stress field inversion method provided by an embodiment of the invention; as shown in fig. 1, the method comprises the following steps:
s101, acquiring a three-dimensional model of a door machine, determining the displacement freedom degree and initial load of the door machine, and determining initial stress field distribution corresponding to the initial load;
s102, taking the displacement freedom degree as a constraint condition, taking an initial load as a boundary condition, carrying out finite element analysis on stress distribution on a gantry crane based on the gantry crane three-dimensional model, and determining all stress extreme points on the gantry crane; the stress extreme points include: stress maximum point and stress minimum point;
s103, obtaining a stress value at a stress extreme point on the measured gantry crane, accumulating the measured stress value and a stress initial value determined based on an initial load to obtain an actual stress value at the stress extreme point, then carrying out interpolation inversion on the actual stress value based on a gantry crane three-dimensional model by adopting an isogeometric analysis method based on non-uniform rational B-spline NURBS, and simulating to obtain the stress field distribution of the gantry crane.
In a specific embodiment, the gantry crane stress field inversion method provided by the invention comprises the following steps:
1) Constructing a three-dimensional model of a gantry crane: drawing a 3D geometric model of the gantry crane according to a two-dimensional drawing, firstly constructing a three-dimensional part model, then assembling the three-dimensional part model into a three-dimensional integral model, combining the assembled models into a single model for facilitating finite element analysis, and the model diagram is shown in FIG. 2.
2) Constructing a three-dimensional finite element model of a gantry crane: according to the portal crane model report, a three-dimensional finite element grid model of the portal crane model is established, and each physical and mechanical parameter, construction type, scale, quantity and distribution condition in the model research range are determined, so that displacement constraint, load condition and the like are determined.
3) Determining a measuring point of a part to be measured: and carrying out preliminary finite element analysis on the three-dimensional geometric model, and determining the positions and the number of the points to be measured approximately according to the positions where the stress is concentrated in the analysis result. Fig. 3 is a stress distribution diagram when the deformation ratio is 1 in the finite element analysis, and fig. 4 is a stress distribution diagram when the deformation ratio is 1000 ten thousand in the finite element analysis.
The stress extreme point is the stress maximum point and the stress minimum point.
4) The three-coordinate measuring machine is used for measuring the actual deformation value of the measuring point of the part to be measured: setting a plurality of measuring points, wherein the positions of the measuring points correspond to the measuring points selected by the finite element in the step 3. The stress data sigma of the normal direction of the position of the measuring point is directly measured or converted by using a stress sensor, wherein: σ is a matrix of nn×1 (deformed matrix of nn×1 obtained by the point to be measured).
5) Comparing the finite element model analysis distribution result in the step 3 with the measured data at the corresponding position in the step 4 to obtain measured data, and judging whether the distribution of the comparison result is consistent; if the measured points are consistent with the measured points, the step 6 is entered, otherwise, the measured points and the positions in the step 3 are corrected.
The finite element analysis needs to be performed based on an initial stress value, the initial stress value is determined by an initial load, and the initial load of the gantry crane is determined in the correction process. In the subsequent stress interpolation inversion process, the actual stress value is equal to the initial stress value plus the measured stress value, and the portal crane stress field distribution closer to the actual situation can be obtained based on the inversion of the actual stress value.
6) And carrying out high-dimensional NURBS curved surface interpolation on the measuring points based on NURBS interpolation inversion calculation on the model, wherein in the inversion calculation process, the measuring point data is taken as the node data as the basis to obtain the stress value of the surface points of the three-dimensional model of the gantry crane, and then obtaining the physical stress field distribution of the whole gantry crane.
Specifically, the existing stress-strain analysis is generally based on finite element analysis, but a finite element method is adopted for structural analysis, and an analysis model approximates a discrete three-dimensional geometric model to a finite element grid, so that the analysis is disjointed from the geometric model, and errors exist in a calculation model; the continuity between units is low, such as a C0 continuous displacement field and a discontinuous stress field between linear units, which affect analysis accuracy; when using higher order units, the computational efficiency is greatly reduced.
Therefore, the invention aims to use an equivalent geometric analysis method based on NURBS, and the structural analysis is constructed on the geometric model by adopting the same NURBS basis function to build the geometric model and the calculation model, so that the limitations of low-order continuity of the finite element inter-unit shape function, geometric model errors and the like can be improved. Isogeometric analysis uses control point information and NURBS basis functions to discrete the physical field x (coordinates, displacement, load, etc.):
x(ξ,η)＝∑N D (ξ,η)x D (1)
where (ζ, η) represents the parameter coordinates, x D Is the unknown coefficient at the D-th control point, N D (ζ, η) is a two-dimensional NURBS basis function. Fig. 5 is a spatial discrete result for a geometric block comprising 3 x 3 cells under different schemes. In fig. 5 (a), when a linear unit (p=1) is employed, the NURBS unit is equivalent to a lagrangian unit, and 16 control points coincide with nodes. For the secondary unit (p=2), the biquadratic lagrangian unit is adopted in (b) of fig. 5, the number of nodes rises to 49, and the control point number of the biquadratic NURBS unit in (c) of fig. 5 is only 25. This illustrates that the degree of freedom increase in NURBS unit analysis is relatively small when using higher order unit calculations.
Since higher order sequential cells facilitate smooth representation of the physical field, providing higher computational accuracy, NURBS cells are typically employed twice and above. For inter-cell continuity, the inter-cell continuity of biquadratic NURBS employed in fig. 5 (C) is C1, while the inter-cell continuity of biquadratic lagrangian employed in fig. 5 (b) is always C0, so that the analysis accuracy of NURBS cells of the same order is higher.
For the discrete equilibrium equation f=ku, K is the stiffness matrix, and the displacement vector u and the external load f are associated with the control point. Cell stiffness matrix K constituting a summary matrix K e The calculation formula is as follows:
wherein B represents a strain-displacement matrix and E represents an elastic matrix. The physical area of the cell is denoted omega e Corresponding regions in NURBS parameter space { ζ, η }, areIn the integrated parameter space->The corresponding area in (a) is->Through Jacobian matrix J 2 And J 1 The conversion of the integrated parameter space into NURBS parameter space and then into physical space is realized.
For the two-dimensional plane stress problem, the element composition of the B matrix is:
wherein N is i (i＝1,2,…,n c ) Is represented by n c NURBS surface basis function of each control point, jacobian matrix J 1 The method comprises the following steps:
gaussian integration region to NURBS parameter space [ ζ ] i ,ξ i+1 )×[η j ,η j+1 ) Is converted into:
thus, jacobian matrix J 2 The method comprises the following steps:
in the isogeometric analysis process, NURBS unit division is fixed, and a density-based 'substitute material' method is adopted to ensure non-singularity of an equation set. And (3) giving a weak material with extremely small elastic modulus to the cavity area, wherein in the numerical calculation process, the unit stiffness matrix expression is as follows:
wherein K is solid Representing the physical cell stiffness matrix ρ e Is the relative density of the units.
FIG. 6 is a diagram of an isogeometric analysis combined with extended finite element techniques provided by an embodiment of the present invention; when the isogeometric analysis is performed on the continuum structure in the isogeometric analysis process, as shown in fig. 6C, the rectangular structural domain is constructed by a two-dimensional NURBS basis function and is used as a background grid, wherein the dotted line represents the isogeometric analysis grid. The unit types of a in fig. 6 are divided into three types according to the duty ratio of the entity field part in the unit: 1) hole cells, 2) border cells, 3) solid cells. The quick cell identification scheme needs to identify whether each vertex (4/8 vertices in the case of 2D/3D) of the cell is in the physical domain. A cell is identified as a solid/hole cell if all vertices of the cell are in the solid or hole domain. In the rest of the cases, this cell will be identified as a boundary cell. As shown in fig. 6B, by the above identification scheme, it can be quickly identified that E1 is a hole unit, E2 is a boundary unit, and E3 is a physical unit.
As shown in fig. 7, the quick identification scheme is simple and easy to operate, but when the size of the holes/entities in the cell is smaller than the cell size, there may be a false identification. To address this problem, additional detection is introduced in the entity/hole cells near the boundary cells to ensure the correctness of cell identification. If all additional detection nodes remain within the entity/hole domain, then the element is still identified as an entity/hole element, otherwise the element will be re-identified as a boundary element. And after re-identification, local cell refinement is introduced, and the simulated boundary of the boundary cell is obtained. Notably, local cell refinement is only used to generate simulated boundaries, facilitating subsequent numerical integration, and does not increase the scale of the finite element system.
A flowchart of specific steps of the above process is shown in fig. 8.
Considering that the door machine model has a large number of units and the calculation amount of the high-order NURBS units is large, the parallel computer is considered to carry out NURBS interpolation calculation on the three-dimensional model in the self-grinding simulation software to obtain the stress field distribution of the whole model, the calculated stress magnitude and direction of each measuring point and the corresponding point of the inversion simulation result are compared, the stress inversion effect is verified, the boundary load is adjusted to enable the inversion value of each measuring point to be continuously approximate to the actual measurement value, and when the error between the inversion value of each measuring point and the actual measurement value is within 10%, the inversion result is the final result.
7) And checking the vertical stress, the maximum horizontal main stress, the minimum horizontal main stress contour map and the calculation result map data at different positions according to the calculation result, and analyzing the influence of factors such as the overall stress field distribution, the structure and the like on the stress field distribution, wherein the factors comprise the overall vertical main stress, the maximum horizontal main stress and the minimum horizontal main stress contour map of the gantry crane three-dimensional model.
The inversion method based on the three-dimensional modeling stress field only needs to sample the engineering site, so that the site workload is effectively reduced, sampling of a plurality of measuring points can be realized, and the fine measurement of the stress multiple measuring points is realized; the method can establish a three-dimensional numerical model with engineering large scale and single or complex structure, and invert the stress distribution characteristics of an inversion research area through parallel calculation; visual and visualized overall stress distribution conditions can be obtained, and information such as vertical stress, maximum horizontal main stress, minimum horizontal main stress values and directions of different levels, different inclined planes and different structural positions of the model can be effectively extracted and checked.
FIG. 9 is a diagram of a gantry crane stress field inversion system architecture provided by an embodiment of the present invention; as shown in fig. 9, includes:
a gantry crane model modeling unit 910, configured to build a three-dimensional model of the gantry crane, and determine a displacement degree of freedom and an initial load of the gantry crane;
the extremum point analysis unit 920 is configured to perform finite element analysis on stress distribution on the gantry crane based on the gantry crane three-dimensional model with the displacement degree of freedom as a constraint condition and the initial load as a boundary condition, and determine all stress extremum points on the gantry crane; the stress extreme points include: stress maximum point and stress minimum point;
the stress field inversion unit 930 is configured to obtain a stress value at a stress extreme point on the gantry crane, accumulate the measured stress value with a stress initial value determined based on an initial load to obtain an actual stress value at the stress extreme point, and then perform interpolation inversion on the actual stress value based on the gantry crane three-dimensional model by adopting an isogeometric analysis method based on a non-uniform rational B-spline NURBS, so as to obtain the stress field distribution of the gantry crane through simulation.
It should be understood that, the system is used to execute the method in the foregoing embodiment, and the corresponding program element in the system performs the principle and technical effects similar to those described in the foregoing method, and the working process of the system may refer to the corresponding process in the foregoing method, which is not repeated herein.
Based on the method in the above embodiment, the embodiment of the invention provides an electronic device. The apparatus may include: at least one memory for storing programs and at least one processor for executing the programs stored by the memory. Wherein the processor is adapted to perform the method described in the above embodiments when the program stored in the memory is executed.
Based on the method in the above embodiment, the embodiment of the present invention provides a computer-readable storage medium storing a computer program, which when executed on a processor, causes the processor to perform the method in the above embodiment.
Based on the method in the above embodiments, an embodiment of the present invention provides a computer program product, which when run on a processor causes the processor to perform the method in the above embodiments.
It is to be appreciated that the processor in embodiments of the invention may be a central processing unit (centralprocessing unit, CPU), other general purpose processor, digital signal processor (digital signalprocessor, DSP), application specific integrated circuit (application specific integrated circuit, ASIC), field programmable gate array (field programmable gate array, FPGA) or other programmable logic device, transistor logic device, hardware components, or any combination thereof. The general purpose processor may be a microprocessor, but in the alternative, it may be any conventional processor.
The steps of the method in the embodiment of the present invention may be implemented by hardware, or may be implemented by executing software instructions by a processor. The software instructions may be comprised of corresponding software modules that may be stored in random access memory (random access memory, RAM), flash memory, read-only memory (ROM), programmable ROM (PROM), erasable programmable PROM (EPROM), electrically erasable programmable EPROM (EEPROM), registers, hard disk, a removable disk, a CD-ROM, or any other form of storage medium known in the art. An exemplary storage medium is coupled to the processor such the processor can read information from, and write information to, the storage medium. In the alternative, the storage medium may be integral to the processor. The processor and the storage medium may reside in an ASIC.
In the above embodiments, it may be implemented in whole or in part by software, hardware, firmware, or any combination thereof. When implemented in software, may be implemented in whole or in part in the form of a computer program product. The computer program product includes one or more computer instructions. When loaded and executed on a computer, produces a flow or function in accordance with embodiments of the present invention, in whole or in part. The computer may be a general purpose computer, a special purpose computer, a computer network, or other programmable apparatus. The computer instructions may be stored in or transmitted across a computer-readable storage medium. The computer instructions may be transmitted from one website, computer, server, or data center to another website, computer, server, or data center by a wired (e.g., coaxial cable, fiber optic, digital Subscriber Line (DSL)), or wireless (e.g., infrared, wireless, microwave, etc.). The computer readable storage medium may be any available medium that can be accessed by a computer or a data storage device such as a server, data center, etc. that contains an integration of one or more available media. The usable medium may be a magnetic medium (e.g., a floppy disk, a hard disk, a magnetic tape), an optical medium (e.g., a DVD), or a semiconductor medium (e.g., a Solid State Disk (SSD)), or the like.
It will be appreciated that the various numerical numbers referred to in the embodiments of the present invention are merely for ease of description and are not intended to limit the scope of the embodiments of the present invention.
It will be readily appreciated by those skilled in the art that the foregoing description is merely a preferred embodiment of the invention and is not intended to limit the invention, but any modifications, equivalents, improvements or alternatives falling within the spirit and principles of the invention are intended to be included within the scope of the invention.
Claims (10)
1. A method for inverting a stress field of a door machine, comprising:
acquiring a three-dimensional model of a door machine, and determining the displacement freedom degree and initial load of the door machine;
taking the displacement freedom degree as a constraint condition, taking an initial load as a boundary condition, carrying out finite element analysis on stress distribution on a gantry crane based on the gantry crane three-dimensional model, and determining all stress extreme points on the gantry crane; the stress extreme points include: stress maximum point and stress minimum point; the initial load is used for determining initial stress field distribution on the door machine;
and obtaining a stress value at a stress extreme point on the door machine, accumulating the measured stress value and a stress initial value determined based on an initial load to obtain an actual stress value at the stress extreme point, then carrying out interpolation inversion on the actual stress value based on a door machine three-dimensional model by adopting an isogeometric analysis method based on non-uniform rational B-spline NURBS, and simulating to obtain the stress field distribution of the door machine.
2. The method of claim 1, wherein the initial load is determined by:
presetting an initial load value, carrying out finite element analysis on stress distribution on a gantry crane based on the preset initial load value, the displacement freedom degree and a gantry crane three-dimensional model, and preliminarily determining a stress extreme point on the gantry crane;
acquiring an actual deformation value after a preset load is added at a preliminarily determined stress extreme point on the door machine, comparing the actual deformation value with a deformation value obtained after the preset load is added by finite element analysis simulation, and if the deformation values are inconsistent, adjusting the preset initial load value until the deformation values are consistent;
and taking the initial load value corresponding to the consistent deformation value as the final determined initial load.
3. The method of claim 1, wherein the degree of freedom in displacement is a displacement constraint of a gantry crane, and is determined based on parameters of the gantry crane.
4. A method according to any one of claims 1 to 3, wherein interpolating and inverting the measured stress values based on NURBS-based isogeometric analysis method employed by the gantry crane three-dimensional model comprises:
establishing a geometric model by adopting NURBS basis functions based on a gantry crane three-dimensional model;
dispersing a physical field corresponding to an actual stress value at a stress measurement point based on a NURBS basis function in the established geometric model so as to perform isogeometric analysis, and performing interpolation inversion to obtain the stress of each point on the gantry crane three-dimensional model; the discrete equilibrium equation is determined by the stiffness matrix, gantry crane displacement vector, and external load.
5. The method of claim 4 wherein the NURBS surface interpolation is performed using NURBS basis functions of second order or more.
6. The method of claim 4, wherein the geometric model comprises: a hole unit, a entity unit and a boundary unit;
for void cells, the stiffness matrix K e The method comprises the following steps:
K e ＝K solid *ρ e
wherein K is solid Representing the stiffness matrix, ρ, of the physical unit e Is the rigidity ratio of the cavity unit to the entity unit.
7. The method according to claim 6, wherein at least one additional detection point is introduced in the physical cells and/or the void cells in the vicinity of the boundary cells when identifying each cell in the geometric model, and if the domain in which the at least one additional detection point is located is different from the domain in which the cell is located, the corresponding cell is re-divided into boundary cells.
8. A gantry crane stress field inversion system, comprising:
the gantry crane model modeling unit is used for establishing a three-dimensional model of the gantry crane and determining the displacement freedom degree and initial load of the gantry crane;
the extreme point analysis unit is used for taking the displacement freedom degree as a constraint condition, taking an initial load as a boundary condition, carrying out finite element analysis on stress distribution on the gantry crane based on the gantry crane three-dimensional model, and determining all stress extreme points on the gantry crane; the stress extreme points include: stress maximum point and stress minimum point; the initial load is used for determining initial stress field distribution on the door machine;
and the stress field inversion unit is used for acquiring the stress value at the stress extreme point on the measured gantry crane, accumulating the measured stress value and the stress initial value determined based on the initial load to obtain the actual stress value at the stress extreme point, and then carrying out interpolation inversion on the actual stress value based on the gantry crane three-dimensional model by adopting an isogeometric analysis method based on non-uniform rational B-spline NURBS, so as to obtain the stress field distribution of the gantry crane through simulation.
9. The system according to claim 8, wherein the gantry crane model acquisition unit presets an initial load value, performs finite element analysis on stress distribution on the gantry crane based on the preset initial load value, the displacement freedom degree and the gantry crane three-dimensional model, and preliminarily determines a stress extreme point on the gantry crane; acquiring an actual deformation value after a preset load is added at a preliminarily determined stress extreme point on the door machine, comparing the actual deformation value with a deformation value obtained after the preset load is added by finite element analysis simulation, and if the deformation values are inconsistent, adjusting the preset initial load value until the deformation values are consistent; and taking the initial load value corresponding to the consistent deformation value as the final determined initial load.
10. An electronic device, comprising:
at least one memory for storing a program;
at least one processor for executing the memory-stored program, which processor is adapted to perform the method according to any of claims 1-7, when the memory-stored program is executed.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311272894.5A CN117313471A (en) | 2023-09-27 | 2023-09-27 | Door machine stress field inversion method and system and electronic equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311272894.5A CN117313471A (en) | 2023-09-27 | 2023-09-27 | Door machine stress field inversion method and system and electronic equipment |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN117313471A true CN117313471A (en) | 2023-12-29 |
Family
ID=89284353
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311272894.5A Pending CN117313471A (en) | 2023-09-27 | 2023-09-27 | Door machine stress field inversion method and system and electronic equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117313471A (en) |
-
2023
- 2023-09-27 CN CN202311272894.5A patent/CN117313471A/en active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110442936B (en) | Equipment fault diagnosis method, device and system based on digital twin model | |
| CN109060279B (en) | Error analysis method for measuring bridge deflection by tilt angle sensor | |
| Rezchikov et al. | Mathematical models for estimating the degree of influence of major factors on performance and accuracy of coordinate measuring machines | |
| Franke et al. | Validation of OpenFOAM 1.6. x with the German VDI guideline for obstacle resolving micro-scale models | |
| CN106529055A (en) | Model updating method based on strain modal shape correlation | |
| WO2021208079A1 (en) | Method and apparatus for obtaining power battery life data, computer device, and medium | |
| CN112904139B (en) | High-voltage switch cabinet partial discharge positioning method and system considering temperature field change | |
| CN114547951B (en) | Dam state prediction method and system based on data assimilation | |
| JP2019537079A (en) | How to build stochastic models for large-scale renewable energy data | |
| CN111884207A (en) | Power grid topological structure visualization method, system and medium based on electrical distance | |
| CN113835060A (en) | Power transformer online state monitoring method and system based on digital twinning | |
| CN115455793A (en) | High-rise structure complex component stress analysis method based on multi-scale model correction | |
| Müller et al. | Comparison of different measures for the single point uncertainty in industrial X-ray computed tomography | |
| CN116821637A (en) | Building steel structure data processing method based on data twinning technology | |
| Sun et al. | Uncertainty calibration and quantification of surrogate model for estimating the machining distortion of thin-walled parts | |
| CN113420370B (en) | Intensity design load acquisition method for highly static indefinite structure | |
| CN113158589B (en) | Simulation model calibration method and device of battery management system | |
| Mehrad et al. | Prediction of surface reconstruction uncertainties for freeform surface inspection | |
| CN113808132A (en) | Three-dimensional network model quality detection method and device and computer equipment | |
| Müller et al. | Implementation of parameterized work piece deviations and measurement uncertainties into performant meta-models for an improved tolerance specification | |
| CN111507591B (en) | Power system state determination method, device, computer medium and storage medium | |
| CN117195610A (en) | Slope monitoring and early warning method and device, electronic equipment and readable storage medium | |
| CN116227045B (en) | Local stress strain field construction method and system for structural test piece | |
| CN117313471A (en) | Door machine stress field inversion method and system and electronic equipment | |
| CN115688571A (en) | Method and device for monitoring running state of transformer, electronic equipment and storage medium |
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 |