CN112487610A - Deformation determination method and system for analysis object with complex geometric characteristics - Google Patents

Deformation determination method and system for analysis object with complex geometric characteristics Download PDF

Info

Publication number
CN112487610A
CN112487610A CN202011236013.0A CN202011236013A CN112487610A CN 112487610 A CN112487610 A CN 112487610A CN 202011236013 A CN202011236013 A CN 202011236013A CN 112487610 A CN112487610 A CN 112487610A
Authority
CN
China
Prior art keywords
manifold
analysis object
integral
voxel
boundary
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
Application number
CN202011236013.0A
Other languages
Chinese (zh)
Other versions
CN112487610B (en
Inventor
武杰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hebei University of Technology
Original Assignee
Hebei University of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hebei University of Technology filed Critical Hebei University of Technology
Priority to CN202011236013.0A priority Critical patent/CN112487610B/en
Publication of CN112487610A publication Critical patent/CN112487610A/en
Application granted granted Critical
Publication of CN112487610B publication Critical patent/CN112487610B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开一种具有复杂几何特征分析对象的形变确定方法及系统,涉及力学分析领域,确定方法包括:获取分析对象的几何模型;建立计算网格;根据网格,利用数值流形方法建立多个流形单元、多个数学覆盖和多个物理覆盖;将含有边界的流形单元进行划分,得到多个积分体素;确定不含有边界的积分体素的体素值为第二体素值;并将含有边界的流形单元对应的物理覆盖标记为待更新的物理覆盖;通过连通域更新物理覆盖;根据含有边界的流形单元对应的数学覆盖和更新后的物理覆盖,采用数值流形方法,确定全局位移函数;根据全局位移函数确定分析对象的形变信息。本发明提供的方法及系统能够提高对复杂固体进行力学分析的自动化程度。

Figure 202011236013

The invention discloses a deformation determination method and system for an analysis object with complex geometric features, and relates to the field of mechanical analysis. The determination method includes: acquiring a geometric model of the analysis object; establishing a calculation grid; manifold elements, multiple mathematical covers and multiple physical covers; divide the manifold elements with boundaries to obtain multiple integral voxels; determine the voxel value of the integral voxels without boundaries as the second voxel value ; mark the physical coverage corresponding to the manifold element containing the boundary as the physical coverage to be updated; update the physical coverage through the connected domain; adopt the numerical manifold according to the mathematical coverage corresponding to the manifold element containing the boundary and the updated physical coverage method, the global displacement function is determined; the deformation information of the analysis object is determined according to the global displacement function. The method and system provided by the present invention can improve the automation degree of mechanical analysis of complex solids.

Figure 202011236013

Description

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.

Claims (8)

1.一种具有复杂几何特征分析对象的形变确定方法,其特征在于,包括:1. a deformation determination method with complex geometric feature analysis object, is characterized in that, comprises: 获取分析对象的几何模型;Obtain the geometric model of the analysis object; 在所述几何模型上构建网格;constructing a mesh on the geometric model; 根据所述网格,利用数值流形方法建立多个流形单元、多个数学覆盖和多个物理覆盖;establishing a plurality of manifold elements, a plurality of mathematical covers and a plurality of physical covers using a numerical manifold method according to the grid; 从所述流形单元中获取含有边界的流形单元;obtaining a manifold unit containing a boundary from the manifold unit; 将所述含有边界的流形单元进行划分,得到多个积分体素;dividing the manifold element containing the boundary to obtain a plurality of integral voxels; 确定含有边界的积分体素的体素值为第一体素值;确定不含有边界的积分体素的体素值为第二体素值;并将所述含有边界的流形单元对应的物理覆盖标记为待更新的物理覆盖;Determine the voxel value of the integral voxel with the boundary as the first voxel value; determine the voxel value of the integral voxel without the boundary as the second voxel value; Override physical overrides marked as pending update; 根据所述待更新的物理覆盖中第二体素值对应的积分体素将所述待更新的物理覆盖进行划分,确定所述待更新的物理覆盖包含的连通域;Divide 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 determine a connected domain included in the physical coverage to be updated; 根据所述连通域采用数值流形方法得到更新后的物理覆盖;Using the numerical manifold method to obtain the updated physical coverage according to the connected domain; 根据所述含有边界的流形单元对应的数学覆盖和所述更新后的物理覆盖,采用数值流形方法,确定全局位移函数;According to the mathematical coverage corresponding to the manifold element containing the boundary and the updated physical coverage, a numerical manifold method is used to determine a global displacement function; 根据所述全局位移函数确定所述分析对象的形变信息。Deformation information of the analysis object is determined according to the global displacement function. 2.根据权利要求1所述的具有复杂几何特征分析对象的形变确定方法,其特征在于,所述将所述含有边界的流形单元进行划分,得到多个积分体素,具体包括:2 . The method for determining the deformation of an object with complex geometric features according to claim 1 , wherein the dividing the manifold element containing the boundary to obtain a plurality of integral voxels, specifically comprising: 2 . 将所述含有边界的流形单元按照解析度划分为积分体素。The manifold element containing the boundary is divided into integral voxels according to the resolution. 3.根据权利要求1所述的具有复杂几何特征分析对象的形变确定方法,其特征在于,所述根据所述连通域采用数值流形方法得到更新后的物理覆盖,具体包括:3. The method for determining the deformation of an object with complex geometric features according to claim 1, wherein the method for obtaining the updated physical coverage by using a numerical manifold method according to the connected domain specifically comprises: 采用数值流形方法,根据所述连通域和设定积分体素确定为更新后的物理覆盖;所述设定积分体素为第一体素值对应的且与所述连通域相邻的积分体素。Using the numerical manifold method, the updated physical coverage is determined according to the connected domain and the set integral voxel; the set integral voxel is the integral corresponding to the first voxel value and adjacent to the connected domain voxels. 4.根据权利要求1所述的具有复杂几何特征分析对象的形变确定方法,其特征在于,所述根据所述全局位移函数确定所述分析对象的形变信息,具体包括:4 . The method for determining the deformation of an analysis object with complex geometric features according to claim 1 , wherein the determining the deformation information of the analysis object according to the global displacement function specifically includes: 5 . 根据所述全局位移函数,利用最小势能原理,确定所述分析对象的位移分布;According to the global displacement function, using the principle of minimum potential energy, determine the displacement distribution of the analysis object; 根据所述分析对象的位移分布,利用连续介质力学理论确定所述分析对象的应变和应力分布;According to the displacement distribution of the analysis object, use the continuum mechanics theory to determine the strain and stress distribution of the analysis object; 根据所述分析对象的应变和应力分布确定分析对象的形变信息。Deformation information of the analysis object is determined according to the strain and stress distribution of the analysis object. 5.一种具有复杂几何特征分析对象的形变确定系统,其特征在于,包括:5. A deformation determination system with a complex geometric feature analysis object, characterized in that, comprising: 几何模型获取模块,用于获取分析对象的几何模型;The geometric model acquisition module is used to acquire the geometric model of the analysis object; 构建模块,用于在所述几何模型上构建网格;a building block for building a mesh on the geometric model; 建立模块,用于利用数值流形方法建立多个流形单元、多个数学覆盖和多个物理覆盖;Build modules for building multiple manifold elements, multiple mathematical covers and multiple physical covers using numerical manifold methods; 流形单元获取模块,用于从所述流形单元中获取含有边界的流形单元;a manifold unit acquiring module, configured to acquire a manifold unit containing a boundary from the manifold unit; 划分模块,用于将所述含有边界的流形单元进行划分,得到多个积分体素;a dividing module, used for dividing the manifold element containing the boundary to obtain a plurality of integral voxels; 确定模块,用于确定含有边界的积分体素的体素值为第一体素值;确定不含有边界的积分体素的体素值为第二体素值;并将所述含有边界的流形单元对应的物理覆盖标记为待更新的物理覆盖;A determination module, configured to determine the voxel value of the integral voxel with the boundary as the first voxel value; determine the voxel value of the integral voxel without the boundary as the second voxel value; The physical coverage corresponding to the shape unit is marked as the physical coverage to be updated; 连通域确定模块,用于根据所述待更新的物理覆盖中第二体素值对应的积分体素将所述待更新的物理覆盖进行划分,确定所述待更新的物理覆盖包含的连通域;a connected domain determination module, configured to divide 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 determine the connected domain included in the physical coverage to be updated; 更新模块,用于根据所述连通域采用数值流形方法得到更新后的物理覆盖;an update module for obtaining an updated physical coverage by using a numerical manifold method according to the connected domain; 全局位移函数确定模块,用于根据所述含有边界的流形单元对应的数学覆盖和所述更新后的物理覆盖,采用数值流形方法,确定全局位移函数;a global displacement function determination module for determining a global displacement function by using a numerical manifold method according to the mathematical coverage corresponding to the manifold element containing the boundary and the updated physical coverage; 形变信息确定模块,用于根据所述全局位移函数确定所述分析对象的形变信息。A deformation information determination module, configured to determine the deformation information of the analysis object according to the global displacement function. 6.根据权利要求5所述的具有复杂几何特征分析对象的形变确定系统,其特征在于,所述划分模块,具体包括:6. The deformation determination system with a complex geometric feature analysis object according to claim 5, wherein the division module specifically comprises: 划分单元,用于将所述含有边界的流形单元按照解析度划分为积分体素。The dividing unit is used for dividing the manifold unit containing the boundary into integral voxels according to the resolution. 7.根据权利要求5所述的具有复杂几何特征分析对象的形变确定系统,其特征在于,所述更新模块,具体包括:7. The deformation determination system with complex geometric feature analysis object according to claim 5, wherein the update module specifically comprises: 更新单元,用于采用数值流形方法,根据所述连通域和设定积分体素确定为更新后的物理覆盖;所述设定积分体素为第一体素值对应的且与所述连通域相邻的积分体素。an update unit for using a numerical manifold method to determine the updated physical coverage according to the connected domain and the set integral voxel; the set integral voxel corresponds to the first voxel value and is connected to the Domain-adjacent integral voxels. 8.根据权利要求5所述的具有复杂几何特征分析对象的形变确定系统,其特征在于,所述形变信息确定模块,具体包括:8. The deformation determination system with a complex geometric feature analysis object according to claim 5, wherein the deformation information determination module specifically comprises: 位移分布确定单元,用于根据所述全局位移函数,利用最小势能原理,确定所述分析对象的位移分布;a displacement distribution determining unit, configured to determine the displacement distribution of the analysis object by using the principle of minimum potential energy according to the global displacement function; 应变确定单元,用于根据所述分析对象的位移分布,利用连续介质力学理论确定所述分析对象的应变和应力分布;a strain determination unit, configured to determine the strain and stress distribution of the analysis object by using the continuum mechanics theory according to the displacement distribution of the analysis object; 形变信息确定单元,用于根据所述分析对象的应变和应力分布确定分析对象的形变信息。The deformation information determination unit is used for determining the deformation information of the analysis object according to the strain and stress distribution of the analysis object.
CN202011236013.0A 2020-11-09 2020-11-09 Deformation determination method and system for analytical objects with complex geometric features Active CN112487610B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011236013.0A CN112487610B (en) 2020-11-09 2020-11-09 Deformation determination method and system for analytical objects with complex geometric features

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011236013.0A CN112487610B (en) 2020-11-09 2020-11-09 Deformation determination method and system for analytical objects with complex geometric features

Publications (2)

Publication Number Publication Date
CN112487610A true CN112487610A (en) 2021-03-12
CN112487610B CN112487610B (en) 2021-10-08

Family

ID=74928844

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011236013.0A Active CN112487610B (en) 2020-11-09 2020-11-09 Deformation determination method and system for analytical objects with complex geometric features

Country Status (1)

Country Link
CN (1) CN112487610B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113077460A (en) * 2021-04-28 2021-07-06 河北工业大学 Method and device for determining deformation amount of test piece and storage medium

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060089803A1 (en) * 2002-12-27 2006-04-27 Riken Method and device for numberical analysis of flow field of non-compressive viscous fluid, directly using v-cad data
CN1828672A (en) * 2006-03-28 2006-09-06 中山大学 A Geometric Data Subdivision Method Based on Hexahedral Interpolation Volume Subdivision
CN1912267A (en) * 2005-08-09 2007-02-14 同济大学 Beam unit counting method with thickness
CN101283378A (en) * 2005-08-22 2008-10-08 皇家飞利浦电子股份有限公司 Triangulation method of physical object surface
CN102054296A (en) * 2011-01-20 2011-05-11 西北大学 Grid deformation method based on local rigidity
US20120095738A1 (en) * 2009-06-25 2012-04-19 Asahi Glass Company Limited Calculation method for physical value, numerical analysis method, calculation program for physical value, numerical analysis program, calculation device for physical value, and numerical analysis device
CN105302974A (en) * 2015-11-06 2016-02-03 北京航空航天大学 Real-time cutting simulation method of flexible object on the basis of finite element and time-variant modal analysis
CN106611442A (en) * 2016-11-21 2017-05-03 吴怀宇 Detail shape-preserving 3D processing method and shape-preserving mesh processing system
CN107292953A (en) * 2017-06-23 2017-10-24 深圳市唯特视科技有限公司 A kind of Facial reconstruction method based on unsupervised autocoder
CN109493329A (en) * 2018-11-02 2019-03-19 河北工业大学 Digital Image Correlation Method based on local mesh reflnement

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060089803A1 (en) * 2002-12-27 2006-04-27 Riken Method and device for numberical analysis of flow field of non-compressive viscous fluid, directly using v-cad data
CN1912267A (en) * 2005-08-09 2007-02-14 同济大学 Beam unit counting method with thickness
CN101283378A (en) * 2005-08-22 2008-10-08 皇家飞利浦电子股份有限公司 Triangulation method of physical object surface
CN1828672A (en) * 2006-03-28 2006-09-06 中山大学 A Geometric Data Subdivision Method Based on Hexahedral Interpolation Volume Subdivision
US20120095738A1 (en) * 2009-06-25 2012-04-19 Asahi Glass Company Limited Calculation method for physical value, numerical analysis method, calculation program for physical value, numerical analysis program, calculation device for physical value, and numerical analysis device
CN102054296A (en) * 2011-01-20 2011-05-11 西北大学 Grid deformation method based on local rigidity
CN105302974A (en) * 2015-11-06 2016-02-03 北京航空航天大学 Real-time cutting simulation method of flexible object on the basis of finite element and time-variant modal analysis
CN106611442A (en) * 2016-11-21 2017-05-03 吴怀宇 Detail shape-preserving 3D processing method and shape-preserving mesh processing system
CN107292953A (en) * 2017-06-23 2017-10-24 深圳市唯特视科技有限公司 A kind of Facial reconstruction method based on unsupervised autocoder
CN109493329A (en) * 2018-11-02 2019-03-19 河北工业大学 Digital Image Correlation Method based on local mesh reflnement

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
NIU LING 等: "CFD Simulation Calculation of Inlet Manifold Modification Design", 《IEEE》 *
蔡永昌 等: "流形方法的矩形覆盖系统及其全自动生成算法", 《重庆大学学报(自然科学版)》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113077460A (en) * 2021-04-28 2021-07-06 河北工业大学 Method and device for determining deformation amount of test piece and storage medium

Also Published As

Publication number Publication date
CN112487610B (en) 2021-10-08

Similar Documents

Publication Publication Date Title
CN108052776B (en) Flood model construction and dynamic display method based on BIM and three-dimensional live-action model
Ito et al. Octree‐based reasonable‐quality hexahedral mesh generation using a new set of refinement templates
Fryazinov et al. Multi-scale space-variant FRep cellular structures
CN111243090B (en) Earthwork volume calculating method and system
CN104809362A (en) Method for preparing speckle pattern including heterogeneous deformation based on numerical calculation
CN108108563B (en) Cross-platform hydrodynamics modeling method based on BIM and GIS technology
CN112487610B (en) Deformation determination method and system for analytical objects with complex geometric features
CN109785443B (en) Three-dimensional model simplification method for large ocean engineering equipment
CN110111237B (en) Three-dimensional model coordinate conversion method based on Brep model
CN109933877B (en) Algebraic Multigrid 3D Variational Data Assimilation Method
Anderson et al. Surface and hypersurface meshing techniques for space–time finite element methods
CN111797447B (en) Complex foundation concrete dam pore channel reinforcement method
Ragia et al. 3D reconstruction and visualization of alternatives for restoration of historic buildings: a new approach
CN113297663B (en) BIM technology-based park water flow visualization analysis method
CN107391820A (en) Tunnel excavation supporting gridless routing analysis method
Liu et al. An improved Poisson surface reconstruction algorithm based on the boundary constraints
Agrawal et al. Conversion of a GIS surface data directly to a 3D STL part for terrain modeling
CN110704916B (en) BIM technology-based large complex three-dimensional geological model grid coarsening method
Guo et al. Bézier Segmentation of T-spline Solids in Parametric Domain.
Ripetskyi et al. Analysis of stress-strain state of the metal plate based on discrete data of displacement values
CN117036629B (en) Satellite joint production DSM method, device, system and computer readable storage medium
CN116051782B (en) Data processing and reconstruction modeling method, device and storage medium based on orthogonal grid curve interpolation
CN108765570A (en) A kind of implementation method and system that three-dimensional visualization is merged with numerical simulation
Chen et al. 3D Visualization of electromagnetic environment
Wu et al. Discontinuous mechanical analysis of manifold element strain of rock slope based on open source Gempy

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