Deformation determination method and system for analysis object with complex geometric characteristics
Technical Field
The invention relates to the field of mechanical analysis, in particular to a deformation determination method and system for an analysis object with complex geometric characteristics.
Background
In mechanical analysis, some analysis objects have very complex geometric features, such as texture on the trunk of a tree, criss-cross joints in rock, faults, trabeculae inside bones, etc.; in addition, the geometric characteristics of some analysis objects are changed in the analysis process, such as new boundaries generated by crack propagation, tunnel faces continuously advancing in tunnel excavation, dynamic boundaries generated by underground water level changes and the like.
The mechanical analysis of such objects is mainly represented by numerical methods, in particular finite element methods. There are two major difficulties with finite element mechanical analysis of such objects: firstly, the existing numerical analysis methods generally need a CAD model of an analysis object, and the CAD model for establishing the analysis objects with complex geometric features has low automation degree and huge workload. Secondly, the high-quality computational grid required by the numerical analysis generated based on the CAD model has the difficulties of low automation degree and huge workload.
Disclosure of Invention
The invention aims to provide a deformation determining method and a deformation determining system for an analysis object with complex geometric characteristics, so as to improve the automation degree of mechanical analysis on a solid with complex geometric characteristics.
In order to achieve the purpose, the invention provides the following scheme:
a method of deformation determination of an analysis object having complex geometric features, comprising:
acquiring a geometric model of an analysis object;
constructing a mesh on the geometric model;
establishing a plurality of manifold units, a plurality of mathematical covers and a plurality of physical covers by using a numerical manifold method according to the grid;
acquiring a manifold unit containing a boundary from the manifold unit;
dividing the manifold unit containing the boundary to obtain a plurality of integral voxels;
determining a voxel value of an integral voxel containing a boundary as a first voxel value; determining a voxel value of an integral voxel not containing a boundary as a second voxel value; marking the physical coverage corresponding to the manifold unit with the boundary as the physical coverage to be updated;
dividing the physical coverage to be updated according to an integral voxel corresponding to a second voxel value in the physical coverage to be updated, and determining a connected domain contained in the physical coverage to be updated;
obtaining updated physical coverage by adopting a numerical manifold method according to the connected domain;
determining a global displacement function by adopting a numerical manifold method according to the mathematical coverage corresponding to the manifold unit with the boundary and the updated physical coverage;
and determining deformation information of the analysis object according to the global displacement function.
Optionally, the dividing the manifold unit including the boundary to obtain an integral voxel specifically includes:
and dividing the manifold unit with the boundary into integral voxels according to the resolution.
Optionally, the obtaining the updated physical coverage by using a numerical manifold method according to the connected domain specifically includes:
determining the physical coverage after updating according to the connected domain and the set integral voxel by adopting a numerical manifold method; and the set integral voxel is an integral voxel corresponding to the first voxel value and adjacent to the connected domain.
Optionally, the determining deformation information of the analysis object according to the global displacement function specifically includes:
determining the displacement distribution of the analysis object by utilizing a minimum potential energy principle according to the global displacement function;
according to the displacement distribution of the analysis object, determining the strain and stress distribution of the analysis object by using a continuous medium mechanics theory;
and determining deformation information of the analysis object according to the strain and stress distribution of the analysis object.
A deformation determination system with complex geometric feature analysis objects, comprising:
a geometric model obtaining module for obtaining a geometric model of the analysis object;
a construction module for constructing a mesh on the geometric model;
the system comprises an establishing module, a calculating module and a processing module, wherein the establishing module is used for establishing a plurality of manifold units, a plurality of mathematical covers and a plurality of physical covers by using a numerical manifold method;
the manifold unit acquisition module is used for acquiring a manifold unit containing a boundary from the manifold unit;
the dividing module is used for dividing the manifold unit containing the boundary to obtain a plurality of integral voxels;
a determining module for determining a voxel value of an integral voxel containing a boundary as a first voxel value; determining a voxel value of an integral voxel not containing a boundary as a second voxel value; marking the physical coverage corresponding to the manifold unit with the boundary as the physical coverage to be updated;
a connected component determining module, configured to divide the physical coverage to be updated according to an integral voxel corresponding to a second voxel value in the physical coverage to be updated, and determine a connected component included in the physical coverage to be updated;
the updating module is used for obtaining updated physical coverage by adopting a numerical manifold method according to the connected domain;
a global displacement function determining module, configured to determine a global displacement function by using a numerical manifold method according to the mathematical coverage corresponding to the manifold unit with the boundary and the updated physical coverage;
and the deformation information determining module is used for determining the deformation information of the analysis object according to the global displacement function.
Optionally, the dividing module specifically includes:
and the dividing unit is used for dividing the manifold unit containing the boundary into integral voxels according to the resolution.
Optionally, the update module specifically includes:
the updating unit is used for determining the updated physical coverage according to the connected domain and the set integral voxel by adopting a numerical manifold method; and the set integral voxel is an integral voxel corresponding to the first voxel value and adjacent to the connected domain.
Optionally, the deformation information determining module specifically includes:
the displacement distribution determining unit is used for determining the displacement distribution of the analysis object by utilizing a minimum potential energy principle according to the global displacement function;
the strain determining unit is used for determining the strain and stress distribution of the analysis object by utilizing a continuous medium mechanics theory according to the displacement distribution of the analysis object;
and the deformation information determining unit is used for determining the deformation information of the analysis object according to the strain and stress distribution of the analysis object.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects:
the invention provides a deformation determination method and a deformation determination system for an analysis object with complex geometric characteristics, wherein a plurality of integral voxels are obtained by dividing a manifold unit containing a boundary; determining a voxel value of an integral voxel containing a boundary as a first voxel value; determining a voxel value of an integral voxel not containing a boundary as a second voxel value; marking the physical coverage corresponding to the manifold unit containing the boundary as the physical coverage to be updated; dividing the physical coverage to be updated according to an integral voxel corresponding to a second voxel value in the physical coverage to be updated, and determining a connected domain contained in the physical coverage to be updated; the updated physical coverage is obtained by adopting a numerical manifold method according to the connected domain, and the calculation process avoids the explicit recording of line segments, so that the geometric parameters of the analysis object can be input through a CAD model and also can be input through a digital image. Through digital image input, the workflow with low automation degree and huge workload for establishing a CAD model for an analysis object with complex geometric characteristics can be avoided. The physical coverage is established based on the integral voxel value, and the method is simultaneously suitable for two-dimensional and three-dimensional and has high calculation efficiency. Compared with the prior art, the method is simpler to implement and is more suitable for mechanical analysis of objects with complex geometric characteristics.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings needed to be used in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings without inventive exercise.
FIG. 1 is a flow chart of a method for determining deformation of an object having complex geometric features according to the present invention;
FIG. 2 is a schematic diagram of a geometric model of a deformation determination method for an analysis object with complex geometric features according to the present invention;
FIG. 3 is a schematic diagram of a grid in accordance with the present invention;
FIG. 4 is a schematic diagram of the physical coverage of the present invention;
FIG. 5 is a schematic diagram of a geometric model after a mesh is constructed by the method for determining deformation of an object with complex geometric features according to the present invention;
FIG. 6 is a schematic diagram of the deformation determination method integral voxels of the analysis object with complex geometric features according to the present invention;
FIG. 7 is a schematic diagram of the deformation determination method for an analysis object with complex geometric features according to the present invention, integrating voxel setting voxel values;
FIG. 8 is a schematic diagram of a connected domain of a deformation determination method for an analysis object with complex geometric features according to the present invention;
FIG. 9 is a schematic diagram of a deformation determination system for analysis of objects with complex geometric features according to the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The invention aims to provide a deformation determination method and a deformation determination system for an analysis object with complex geometric characteristics, so as to improve the automation degree of mechanical analysis on complex solids.
In order to make the aforementioned objects, features and advantages of the present invention comprehensible, embodiments accompanied with figures are described in further detail below.
The numerical manifold method refers to a theory of reducing the difficulty of constructing an interpolation grid of a complex boundary problem and a dynamic boundary problem by subdividing coverage, and different fields refer to the theory by adopting different names, such as manifold elements, a unit decomposition method, a virtual point method, an expanded finite element method, a limited coverage method and the like.
As shown in fig. 1, the method for determining deformation of an analysis object with complex geometric features provided by the present invention includes:
step 101: a geometric model of the analysis object is obtained. As shown in fig. 2, the geometric model is a CAD model or a digital image. The digital image can be obtained by CT scanning or camera shooting, and can be any complex geometric model such as a tree, a human skeleton and the like.
Step 102: as shown in fig. 3, a grid for numerical manifold method calculation is constructed with a grid cell size of h. The computational grid may be established by network partitioning software, such as Gmsh; or by simple computer programming.
Step 103: a plurality of manifold cells, a plurality of mathematical overlays and a plurality of physical overlays are established from the grid using a numerical manifold method. A weight function is established on the mathematical overlay, and a local displacement function is established on the physical overlay. As shown in fig. 4, fig. 4(a) is a physical overlay centered on point 1, and fig. 4(b) is a physical overlay centered on point 2.
Step 104: as shown in fig. 5, the computational mesh is overlaid on the analysis object, and a manifold cell including a boundary is obtained from the manifold cells.
Step 105: and dividing the manifold unit containing the boundary to obtain a plurality of integral voxels.
Dividing the manifold unit containing the boundary to obtain an integral voxel, and specifically comprises the following steps:
and dividing the manifold unit with the boundary into integral voxels according to the resolution. As shown in fig. 6, in the present embodiment, a resolution of 10 × 10 is taken as an example.
Step 106: as shown in fig. 7, determining a voxel value of an integral voxel containing a boundary as a first voxel value; determining a voxel value of an integral voxel not containing a boundary as a second voxel value; and marking the physical coverage corresponding to the manifold unit with the boundary as the physical coverage to be updated. In this embodiment, the first voxel value is set to 0, and the second voxel value is set to 1.
Step 107: and dividing the physical coverage to be updated according to the integral voxel corresponding to the second voxel value in the physical coverage to be updated, and determining a connected domain contained in the physical coverage to be updated. Fig. 8(a) shows 3 connected domains generated by physical coverage centered on point 1, and fig. 8(b) shows 2 connected domains generated by physical coverage centered on point 2, and the connected domains are distinguished by different gradations.
Step 108: and obtaining the updated physical coverage by adopting a numerical manifold method according to the connected domain.
Obtaining the updated physical coverage by adopting a numerical manifold method according to the connected domain, specifically comprising:
determining the physical coverage after updating according to the connected domain and the set integral voxel by adopting a numerical manifold method; and the set integral voxel is an integral voxel corresponding to the first voxel value and adjacent to the connected domain. The integral voxel with the voxel value of 0 is classified into a physical coverage adjacent to the integral voxel. And regarding each connected domain as a physical coverage, wherein the number and the shape of the connected domains are the number and the shape of the updated physical coverage. According to the theory of the numerical manifold method, a local displacement function is given to the increased physical coverage, and the local displacement functions of all the physical coverage can be obtained.
Step 109: and determining a global displacement function by adopting a numerical manifold method according to the mathematical coverage corresponding to the manifold unit with the boundary and the updated physical coverage. According to the numerical manifold method, the global displacement function is uh(x)=w(x)ul(x) In that respect The argument x is the coordinate of an arbitrary point in the analysis target. w (x) is a weight function defined over the mathematical overlay. u. ofl(x) To define a local displacement function over physical coverage, where l denotes local meaning. u. ofh(x) For the global displacement function, h is the size of the grid cell used for the numerical manifold method.
Step 110: and determining deformation information of the analysis object according to the global displacement function.
Determining deformation information of the analysis object according to the global displacement function, specifically including:
and determining the displacement distribution of the analysis object by utilizing a minimum potential energy principle according to the global displacement function.
And determining the strain and stress distribution of the analysis object by utilizing a continuous medium mechanics theory according to the displacement distribution of the analysis object.
And determining deformation information of the analysis object according to the strain and stress distribution of the analysis object. The deformation information and the combination strength criterion judge the area of the analysis object which is easy to damage, thereby realizing the mechanical analysis of the analysis object.
As shown in fig. 9, a deformation determination system with a complex geometric feature analysis object includes:
a geometric model obtaining module 701, configured to obtain a geometric model of the analysis object.
A building module 702 configured to build a mesh on the geometric model.
A building module 703, configured to build a plurality of manifold units, a plurality of mathematical overlays and a plurality of physical overlays according to the grid by using a numerical manifold method.
A manifold unit obtaining module 704, configured to obtain a manifold unit including a boundary from the manifold unit.
A dividing module 705, configured to divide the manifold unit including the boundary to obtain a plurality of integral voxels.
A determining module 706 for determining a voxel value of an integral voxel containing a boundary as a first voxel value; determining a voxel value of an integral voxel not containing a boundary as a second voxel value; and marking the physical coverage corresponding to the manifold unit with the boundary as the physical coverage to be updated.
A connected component determining module 707, configured to divide the physical coverage to be updated according to an integral voxel corresponding to a second voxel value in the physical coverage to be updated, and determine a connected component included in the physical coverage to be updated.
An updating module 708, configured to obtain an updated physical coverage by using a numerical manifold method according to the connected domain.
A global displacement function determining module 709, configured to determine a global displacement function by using a numerical manifold method according to the mathematical coverage corresponding to the manifold unit with the boundary and the updated physical coverage.
And a deformation information determining module 710, configured to determine deformation information of the analysis object according to the global displacement function.
The dividing module 705 specifically includes:
and the dividing unit is used for dividing the manifold unit containing the boundary into integral voxels according to the resolution.
The updating module 708 specifically includes:
the updating unit is used for determining the updated physical coverage according to the connected domain and the set integral voxel by adopting a numerical manifold method; and the set integral voxel is an integral voxel corresponding to the first voxel value and adjacent to the connected domain.
The deformation information determining module 710 specifically includes:
and the displacement distribution determining unit is used for determining the displacement distribution of the analysis object by utilizing a minimum potential energy principle according to the global displacement function.
And the strain determining unit is used for determining the strain and stress distribution of the analysis object by utilizing a continuous medium mechanics theory according to the displacement distribution of the analysis object.
And the deformation information determining unit is used for determining the deformation information of the analysis object according to the strain and stress distribution of the analysis object.
The method and the device provided by the invention can be used for directly carrying out mechanical analysis based on the digital image, and conveniently simulating the complex geometric characteristics of the object. The new method has strong theoretical value and great engineering application value. The method takes solving the displacement field of an analysis object as an example, but the problems of solving temperature, seepage and the like can be solved by adopting the physical coverage generation method; the invention takes a square integral voxel as an example, and other integral voxels can be used. The method and the system provided by the invention are not only suitable for two-dimensional space, but also suitable for three-dimensional space.
The embodiments in the present description are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments are referred to each other. For the system disclosed by the embodiment, the description is relatively simple because the system corresponds to the method disclosed by the embodiment, and the relevant points can be referred to the method part for description.
The principles and embodiments of the present invention have been described herein using specific examples, which are provided only to help understand the method and the core concept of the present invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, the specific embodiments and the application range may be changed. In view of the above, the present disclosure should not be construed as limiting the invention.