CN115295094A - Two-dimensional numerical model modeling method for material containing holes - Google Patents
Two-dimensional numerical model modeling method for material containing holes Download PDFInfo
- Publication number
- CN115295094A CN115295094A CN202210796584.2A CN202210796584A CN115295094A CN 115295094 A CN115295094 A CN 115295094A CN 202210796584 A CN202210796584 A CN 202210796584A CN 115295094 A CN115295094 A CN 115295094A
- Authority
- CN
- China
- Prior art keywords
- holes
- hole
- dimensional
- parameter
- material containing
- 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
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C60/00—Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/20—Drawing from basic elements, e.g. lines or circles
- G06T11/203—Drawing of straight lines or curves
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Computing Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention provides a two-dimensional numerical model modeling method for a material containing holes, which adopts C + + programming language as a bottom code and Javascript programming language as a script language. The modeling work of holes with any number and any shapes can be realized, and the size parameters of the holes can be accurately set to obey a certain specific distribution function (such as uniform distribution, normal distribution and Weber distribution). And then obtaining a two-dimensional numerical computation grid by utilizing Gmsh grid division software. The invention provides a method for rapidly acquiring a two-dimensional numerical grid of a hole-containing material for scientific research personnel and engineering design personnel, accurately sets the number, shape characteristics and spatial positions of holes, and can be used in the fields of numerical calculation of the hole-containing material and the like. The modeling method not only can simply and quickly construct the two-dimensional numerical model according to the actual size parameters of the holes, but also can simply and quickly construct the two-dimensional numerical model of the material containing the holes according to certain distribution characteristics.
Description
Technical Field
The invention belongs to the technical field of hole-containing materials, relates to numerical calculation, and particularly relates to a two-dimensional numerical model modeling method of a hole-containing material.
Background
The materials containing holes (rock, concrete, etc.) are widely used in the fields of national defense engineering, capital construction engineering, mine engineering, etc. Rock can be regarded as a natural material consisting of mineral particles and cementitious substances under long-term geological effects, and relevant indoor observations indicate that a large number of micro-cavities are present inside the rock. Concrete is an artificial material consisting of cement and sand, and related indoor observation results also show that a large number of micro-pores exist in the concrete. Stress concentration phenomena are generated near the micro-holes under the action of external load, so that cracks are initiated and expanded, and a plurality of current research results show that the micro-holes in the material have important influence on mechanical constitutive characteristics and mechanical parameters of the material.
In order to study the influence of the micro-cavities on the mechanical constitutive characteristics and mechanical parameters of the materials, researchers began to conduct research work based on numerical simulation. The pretreatment in the numerical simulation is a key link, but a simple and effective method for obtaining the two-dimensional numerical model of the material containing the holes is not available at present. On one hand, quick batch modeling work cannot be carried out on a large number of holes, and on the other hand, the size parameters and the spatial position parameters of the holes cannot be simply and quickly set. Not only is a great deal of time consumed for pretreatment work, but also the constructed two-dimensional numerical model cannot accurately reflect the micropore characteristics of the material. Therefore, a simple and effective two-dimensional numerical model modeling method for the material containing the holes is needed to be established.
Disclosure of Invention
Aiming at the defects in the prior art, the invention aims to provide a two-dimensional numerical model modeling method for a material containing holes, and solve the technical problem that the pretreatment work of numerical calculation in the prior art consumes long time.
In order to solve the technical problems, the invention adopts the following technical scheme:
a two-dimensional numerical model modeling method of a material containing holes comprises the following steps:
step 1, acquiring the total number N of holes needing to be characterized in a modeling stage in an object to be researched based on an indoor observation means;
step 2, aiming at N holes in the object to be researched, dividing the holes into holes with five shapes of triangle, rectangle, polygon, circle and ellipse by adopting an approximate method according to the shape characteristics of the holes, and respectively obtaining the number N of the holes with each shape i ；
Step 3, respectively acquiring the size parameters of each hole aiming at five types of holes, namely triangular holes, rectangular holes, polygonal holes, circular holes and elliptical holes;
step 4, aiming at five types of triangular, rectangular, polygonal, circular and elliptical holes, if the accurate size parameter of each hole cannot be obtained, but the size parameter of each hole is known to obey one of a uniform distribution type, a normal distribution type and a Weber distribution type, the size parameter of each hole is obtained by setting the parameter of the distribution type;
in step 4, the size parameter L is calculated as:
for uniform distribution: l = L min +(L max -L min )α；
For weber distribution:
for normal distributions:
in the formula:
alpha is a random value uniformly distributed in the interval of [0,1 ];
L min is the lower limit value of the size parameter;
L max the upper limit value of the size parameter;
λ is a proportional parameter;
k is a shape parameter;
mu is a mean value;
sigma is a standard deviation;
step 5, respectively acquiring spatial position parameters of each hole aiming at five types of holes, namely triangular holes, rectangular holes, polygonal holes, circular holes and elliptical holes;
the spatial position parameters comprise coordinate values X in the X direction and coordinate values Y in the Y direction;
step 6, firstly, establishing a two-dimensional geometric model of the whole object to be researched when holes are not considered, wherein the outer contour line of the geometric model can be in any irregular shape or regular shape;
step 7, according to the quantity, the size parameters and the spatial position parameters of the holes obtained in the steps 2, 3, 4 and 5, adopting Javascript script scripting language and utilizing a circular iteration method to establish a geometric model of each hole in the two-dimensional geometric model of the whole object to be researched;
step 8, completing the modeling work of the two-dimensional geometric model of the material containing the holes after the geometric models of all the holes are established; outputting the geometric data into a file in a geo file format according to the format type;
and 9, importing the file in the geo file format into Gmsh mesh generation software, and obtaining a two-dimensional numerical model of the material containing the holes by utilizing the 2D mesh generation function.
The invention also has the following technical characteristics:
in the step 1, the indoor observation means comprises a capillary pressure curve method, a CT scanning technology, a scanning electron microscope method and a nuclear magnetic resonance technology.
In step 3, for triangular holes, the dimensional parameters include the length l of each side i (ii) a For a rectangular hole, the dimensional parameter includes length a i And width b i (ii) a For polygonal holes, the dimensional parameters include the length of each side, l i (ii) a For a circular hole, the dimensional parameter includes the radius r i For oval holes, the dimensional parameters include the major axis a i And a short half shaft b i 。
Step 5, if the spatial position parameter of each hole cannot be accurately obtained, setting the spatial position parameter in a randomization manner; the spatial position parameters of the holes in the whole model can be randomly distributed, and the holes in the local area of the model can also be randomly distributed;
the coordinate values x and y are calculated as:
in the formula:
alpha represents a random value uniformly distributed in the interval of [0,1 ];
X min is the lower limit value of the coordinate value in the X direction;
X max is the upper limit value of the coordinate value in the X direction;
Y min is the lower limit value of the coordinate value in the Y direction;
Y max the upper limit value of the Y-direction coordinate value.
In step 6, the regular shapes include triangle, rectangle, circle and ellipse.
In step 8, the format types include four major format types, namely Point, line, lineLoop and Surface.
Compared with the prior art, the invention has the following technical effects:
the modeling method not only can simply and quickly construct the two-dimensional numerical model according to the actual size parameters of the holes, but also can simply and quickly construct the two-dimensional numerical model of the material containing the holes according to certain distribution characteristics.
The modeling method can greatly shorten the time of a pretreatment part in the numerical calculation process, can also represent the complex hole distribution characteristics in a two-dimensional numerical model, and lays a foundation for researching the damage process of the hole-containing material under the action of external load.
The modeling method can be applied to the fields of numerical calculation of materials containing holes and the like. The invention provides a method for rapidly acquiring a two-dimensional numerical model of a material containing holes for scientific research personnel and engineering design personnel, and realizes accurate setting of the number, shape characteristics and spatial positions of the holes in the modeling process.
Drawings
FIG. 1 is a flow chart of a two-dimensional numerical model modeling method of a material containing pores according to the present invention.
FIG. 2 is a schematic representation of a two-dimensional model of the void-containing material of the present invention.
FIG. 3 is a schematic representation of a geometric model of a two-dimensional void-containing material of the present invention.
Fig. 4 is a schematic diagram of a mesh model of a two-dimensional material containing pores according to the present invention, wherein (a) in fig. 4 is a schematic diagram of an overall mesh model, and (b) in fig. 4 is a schematic diagram of a partial mesh model.
The present invention will be explained in further detail with reference to examples.
Detailed Description
It should be noted that all models and parameters in the present invention, if not specifically mentioned, all use the models and parameters known in the art.
Currently, for materials (such as rocks, concrete and the like) with a plurality of holes inside, there is no simple and effective method for quickly constructing holes inside a model, so that the pretreatment work of numerical calculation consumes a lot of time. In order to solve the problems, the invention provides a two-dimensional numerical model modeling method of a material containing holes,
referring to fig. 2, in order to establish a two-dimensional numerical model of a material containing holes, the invention utilizes a C + + programming language and a Javascript programming language to realize the modeling work of holes with any number and any shape by a parameterized modeling method.
The method adopts C + + programming language as bottom layer code, and Javascript programming language is script language. The modeling work of holes with any number and any shapes can be realized, and the size parameters of the holes can be accurately set to obey a certain specific distribution function (such as uniform distribution, normal distribution and Weber distribution). And then obtaining a two-dimensional numerical computation grid by utilizing Gmsh grid generation software. The invention provides a method for rapidly acquiring a two-dimensional numerical grid of a hole-containing material for scientific research personnel and engineering design personnel, accurately sets the number, shape characteristics and spatial positions of holes, and can be used in the fields of numerical calculation of the hole-containing material and the like.
The following embodiments are given as examples of the present invention, and it should be noted that the present invention is not limited to the following embodiments, and all equivalent changes based on the technical solutions of the present invention are included in the protection scope of the present invention.
Example (b):
this embodiment provides a two-dimensional numerical model modeling method for a material containing pores, as shown in fig. 1, the method includes the following steps:
step 1, acquiring the total number N of holes needing to be characterized in a modeling stage in an object to be researched based on an indoor observation means;
in the step 1, the indoor observation means comprises a capillary pressure curve method, a CT scanning technology, a scanning electron microscope method and a nuclear magnetic resonance technology. The indoor observation means are known and commonly used means.
Step 2, aiming at the N holes in the object to be researched, an approximate method is adopted according to the shape characteristics of the holes, as shown in fig. 2, the holes are divided into five types of holes, namely triangular, rectangular, polygonal, circular and elliptical holes, and the number N of the holes in each shape is respectively obtained i ；
Step 3, respectively acquiring the size parameters of each hole aiming at five types of holes, namely triangular holes, rectangular holes, polygonal holes, circular holes and elliptical holes;
in step 3, for triangular holes, the dimensional parameters include the length l of each side i (ii) a For a rectangular hole, the dimension parameter includes length a i And width b i (ii) a For polygonal holes, the dimensional parameters include the length of each side, l i (ii) a For a circular hole, the dimensional parameter includes the radius r i For an elliptical hole, the dimensional parameters include the major axis a i And a short half shaft b i 。
Step 4, aiming at five types of triangular, rectangular, polygonal, circular and elliptical holes, if the accurate size parameter of each hole cannot be obtained, but the size parameter of each hole is known to obey one of a uniform distribution type, a normal distribution type and a Weber distribution type, the size parameter of each hole is obtained by setting the parameter of the distribution type;
in step 4, for uniform distribution, a lower limit value L of the size parameter subject to uniform distribution is set min And an upper limit value L max (ii) a Aiming at the Weber distribution, a proportion parameter lambda and a shape parameter k when the size parameter obeys the Weber distribution are required to be set; setting a mean value mu and a standard deviation sigma when the size parameter obeys normal distribution aiming at the normal distribution;
in step 4, the size parameter L is calculated as:
for uniform distribution: l = L min +(L max -L min )α；
For weber distribution:
for normal distributions:
in the formula:
alpha is a random value uniformly distributed in the interval of [0,1 ];
L min the lower limit value of the size parameter;
L max the upper limit value of the size parameter;
λ is a proportional parameter;
k is a shape parameter;
mu is a mean value;
σ is the standard deviation.
Step 5, respectively acquiring spatial position parameters of each hole aiming at five types of holes, namely triangular holes, rectangular holes, polygonal holes, circular holes and elliptical holes;
the spatial position parameter comprises a coordinate value X in the X direction and a coordinate value Y in the Y direction.
Step 5, if the spatial position parameter of each hole cannot be accurately obtained, setting the spatial position parameter in a randomization manner; not only can be arranged in the whole model and the spatial position parameters of the holes are randomly distributed, but also can be arranged in the local area of the model and the holes are randomly distributed
The coordinate values x and y are calculated as:
in the formula:
alpha represents a random value uniformly distributed in the interval of [0,1 ];
X min is the lower limit value of the coordinate value in the X direction;
X max is the upper limit value of the coordinate value in the X direction;
Y min is the lower limit value of the coordinate value in the Y direction;
Y max the upper limit value of the Y-direction coordinate value.
Step 6, firstly, establishing a two-dimensional geometric model of the whole object to be researched when the hole is not considered, wherein the outer contour line of the geometric model can be in any irregular shape or regular shape;
in step 6, the regular shapes include triangle, rectangle, circle and ellipse.
In this step, the specific process of establishing the two-dimensional geometric model may be a known common specific process of establishing a two-dimensional geometric model.
Step 7, according to the number, size parameters and spatial position parameters of the holes obtained in the steps 2, 3, 4 and 5, adopting Javascript scripting language, and establishing a geometric model of each hole in the two-dimensional geometric model of the whole object to be researched by using a cyclic iteration method, as shown in FIG. 3;
in this step, the specific process of building the geometric model of the hole may be a known and commonly used process of building the geometric model of the hole.
Step 8, completing the two-dimensional geometric model modeling work of the material containing the holes after the geometric models of all the holes are established; outputting the geometric data into a file in a geo file format according to the format type;
in step 8, the format types include four major format types, namely Point, line, lineLoop and Surface.
And 9, importing the file in the geo file format into Gmsh mesh generation software, and obtaining a two-dimensional numerical model of the material containing the holes by utilizing a 2D mesh generation function, as shown in figure 4. Then, known and commonly used related numerical calculation software can be introduced to carry out numerical simulation work.
Claims (6)
1. A two-dimensional numerical model modeling method of a material containing holes is characterized by comprising the following steps:
step 1, acquiring the total number N of holes to be characterized in a modeling stage in an object to be researched based on an indoor observation means;
step 2, aiming at N holes in the object to be researched, dividing the holes into holes with five shapes of triangle, rectangle, polygon, circle and ellipse by adopting an approximate method according to the shape characteristics of the holes, and respectively obtaining the number N of the holes with each shape i ；
Step 3, respectively acquiring the size parameters of each hole aiming at five types of holes, namely triangular holes, rectangular holes, polygonal holes, circular holes and elliptical holes;
step 4, aiming at five types of triangular, rectangular, polygonal, circular and elliptical holes, if the accurate size parameter of each hole cannot be obtained, but the size parameter of each hole is known to obey one of a uniform distribution type, a normal distribution type and a Weber distribution type, the size parameter of each hole is obtained by setting the parameter of the distribution type;
in step 4, the size parameter L is calculated as:
for uniform distribution: l = L min +(L max -L min )α；
For weber distributions:
for normal distributions:
in the formula:
alpha is a random value which is uniformly distributed in the interval of [0,1 ];
L min the lower limit value of the size parameter;
L max the upper limit value of the size parameter;
lambda is a proportional parameter;
k is a shape parameter;
mu is a mean value;
σ is the standard deviation;
step 5, respectively acquiring spatial position parameters of each hole aiming at five types of holes, namely triangular holes, rectangular holes, polygonal holes, circular holes and elliptical holes;
the spatial position parameters comprise coordinate values X in the X direction and coordinate values Y in the Y direction;
step 6, firstly, establishing a two-dimensional geometric model of the whole object to be researched when the hole is not considered, wherein the outer contour line of the geometric model can be in any irregular shape or regular shape;
step 7, according to the number, the size parameters and the spatial position parameters of the holes obtained in the steps 2, 3, 4 and 5, adopting Javascript scripting language and utilizing a cyclic iteration method to establish a geometric model of each hole in the two-dimensional geometric model of the whole object to be researched;
step 8, completing the modeling work of the two-dimensional geometric model of the material containing the holes after the geometric models of all the holes are established; outputting the geometric data into a file in a geo file format according to the format type;
and 9, importing the file in the geo file format into Gmsh mesh generation software, and obtaining a two-dimensional numerical model of the material containing the holes by utilizing the 2D mesh generation function.
2. The method according to claim 1, wherein the indoor observation means in step 1 comprises capillary pressure curve method, CT scanning technique, scanning electron microscope method and nuclear magnetic resonance technique.
3. The method according to claim 1, wherein in step 3, the dimensional parameters include the length of each side/, for the triangular hole i (ii) a For a rectangular hole, the dimensional parameter includes length a i And width b i (ii) a For polygonal holes, the dimensional parameters include the length of each side, l i (ii) a For a circular hole, the dimensional parameter includes the radius r i For oval holes, the dimensional parameters include the major axis a i And a short semi-axis b i 。
4. The method according to claim 1, wherein in step 5, if the spatial location parameter of each hole cannot be obtained accurately, the spatial location parameter is set in a randomized manner; the spatial position parameters of the holes in the whole model can be randomly distributed, and the holes in the local area of the model can also be randomly distributed;
the coordinate values x and y are calculated as:
in the formula:
alpha represents a random value uniformly distributed in the interval of [0,1 ];
X min is the lower limit value of the coordinate value in the X direction;
X max is the upper limit value of the coordinate value in the X direction;
Y min is the lower limit value of the coordinate value in the Y direction;
Y max the upper limit value of the Y-direction coordinate value.
5. The method according to claim 1, wherein the regular shapes in step 6 include triangular, rectangular, circular and elliptical shapes.
6. The method for modeling a two-dimensional numerical model of a material containing pores according to claim 1, wherein in step 8, the format types include four types, namely Point, line, lineLoop and Surface format types.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210796584.2A CN115295094A (en) | 2022-07-06 | 2022-07-06 | Two-dimensional numerical model modeling method for material containing holes |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210796584.2A CN115295094A (en) | 2022-07-06 | 2022-07-06 | Two-dimensional numerical model modeling method for material containing holes |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115295094A true CN115295094A (en) | 2022-11-04 |
Family
ID=83822958
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210796584.2A Pending CN115295094A (en) | 2022-07-06 | 2022-07-06 | Two-dimensional numerical model modeling method for material containing holes |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115295094A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116525041A (en) * | 2023-04-28 | 2023-08-01 | 江南大学 | Modeling method and performance prediction method for metal porous structure |
-
2022
- 2022-07-06 CN CN202210796584.2A patent/CN115295094A/en active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116525041A (en) * | 2023-04-28 | 2023-08-01 | 江南大学 | Modeling method and performance prediction method for metal porous structure |
| CN116525041B (en) * | 2023-04-28 | 2024-01-26 | 江南大学 | Modeling method and performance prediction method for metal porous structure |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114462124B (en) | Method for establishing and numerically simulating three-dimensional multiphase mesoscopic model of concrete | |
| Ouyang et al. | 3D meso-scale modeling of concrete with a local background grid method | |
| CN115295094A (en) | Two-dimensional numerical model modeling method for material containing holes | |
| CN109977481B (en) | Method for creating concrete assembly template and member steel bar based on CATIA | |
| Singh et al. | CFD simulation of the wind field around pyramidal roofed single-story buildings | |
| CN109979009A (en) | A kind of cement particle material three-dimensional with functionally gradient carefully sees mathematical model reconstructing method | |
| CN103345560B (en) | Based on the hydrated cementitious emulation of microscopic image feature and the method for macroscopic properties prediction | |
| CN112116713A (en) | High-precision steel bar arrangement method for linear bearing platform type components | |
| CN107368660B (en) | A kind of analysis method of predetermined period lattice material yield surface | |
| CN104298874A (en) | Meso-structure simulation based concrete creeping property prediction method | |
| Xu et al. | 3D mesostructure generation of fully-graded concrete based on hierarchical point cloud and aggregate coarsening | |
| Saucedo-Mora et al. | Method for the explicit insertion of microstructure in Cellular Automata Finite Element (CAFE) models based on an irregular tetrahedral Finite Element mesh: Application in a multi-scale Finite Element Microstructure MEshfree framework (FEMME) | |
| CN107063961A (en) | A kind of characterizing method of hydrated cementitious three dimensional pore structures feature | |
| CN117116390A (en) | Discretization method suitable for complex computational domain of grid discrete particle model | |
| Pugnale et al. | Computational morphogenesis: Design of freeform surfaces | |
| Jin et al. | Mesoscopic finite element modeling of concrete considering geometric boundaries of actual aggregates | |
| CN117237562A (en) | Method for generating three-dimensional polyhedral random aggregate model | |
| Cai et al. | Intelligent building system for 3D construction of complex brick models | |
| CN113343510B (en) | Two-dimensional random fracture grid generation method | |
| CN114282331B (en) | Method for generating three-dimensional grid of natural river channel | |
| CN113343423B (en) | Random fracture network generation method based on intensity spatial variability | |
| Tang et al. | A positivity-preserving conservative semi-Lagrangian multi-moment global transport model on the cubed sphere | |
| Yuan et al. | Improved random aggregate model for numerical simulations of concrete engineering simulations of concrete engineering | |
| CN113627040A (en) | Heterogeneous slope stability analysis method | |
| CN113343474A (en) | Microscopic analysis method of broken stone reinforced polymer composite material |
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 |