CN104952108B - A kind of CT inversely changes the grid model optimization method of modeling technique - Google Patents
A kind of CT inversely changes the grid model optimization method of modeling technique Download PDFInfo
- Publication number
- CN104952108B CN104952108B CN201510256713.9A CN201510256713A CN104952108B CN 104952108 B CN104952108 B CN 104952108B CN 201510256713 A CN201510256713 A CN 201510256713A CN 104952108 B CN104952108 B CN 104952108B
- Authority
- CN
- China
- Prior art keywords
- unit
- group
- inclusion material
- model
- constituent unit
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Abstract
The invention discloses a kind of CT inversely changes the grid model optimization method of modeling technique, for solving to carry out in three-dimensional grid model process of reconstruction using CT scan picture, as grid component decision criteria limits the problem for causing inclusion material constituent content in grid model to be significantly lower than real content, initial mesh model is imported by FLAC3D software by interface routine, then according to the volume content of each component counted before stress and strain model, by the Fish Programming with Pascal Language in FLAC3D software, the inclusion material constituent unit quantity in initial mesh model is compensated.Reach under conditions of the distribution of inclusion material Component Space is not changed as far as possible, inclusion material constituent unit quantity in model is compensated, the deviation occurred in reverseization modeling process is reduced, the purpose of more accurate numerical model is provided to follow-up numerical simulation work.
Description
Technical field
The invention belongs to mathematic calculation field, is related to grid model optimization method, more particularly to a kind of CT inversely changes
The grid model optimization method of modeling technique.
Background technology
ICT (Industrial Computed Tomography) technology is industrial computed tomography or claims industry meter
Calculation machine Tomography, as the Dynamic Non-Destruction Measurement that a class is advanced, industrial production and scientific research field application increasingly
Extensively, combined with corresponding professional software using CT scan image and three-dimensional grid object model can be built, wherein used MIMICS
During software carries out reconstructing three-dimensional model, different components can be divided to natural material according to the difference of gray scale pre-value.Logical
Cross and 3 MATIC softwares are combined, further three-dimensional model gridding division can be carried out to model, and imports finite element software
Carry out numerical computations.
Yet with the restriction of grid decision criteria, the content for causing mineral to be mingled with component is occurred in that after grid component
Significantly reduce phenomenon.By reducing unit size, dividing elements precision can be improved, but while can also greatly increase numerical value
The amount of calculation of simulation.Three-dimensional grid model in such as Fig. 1 is the standard specimen model of diameter 50mm, height 100mm, contains
357565 tetrahedron elements, unit are referred to：Grid in three-dimensional grid model.The cell-average length of side is about 2mm, if by unit
Size reduction half, appointment generate more than 2860000 unit of 23 × 357565 ≈, and unit number is excessive, and computing hardware is required relatively
Height, while calculating speed can be had a strong impact on.Based on above-mentioned situation, in the urgent need to one kind is simple and reliable, calculate model is not affected
Under speed conditions, the CT that can effectively improve model counting accuracy inversely changes the grid model optimization method of modeling technique, effectively
Solve to be significantly lower than real content as grid component decision criteria limits inclusion material constituent content in the grid model for causing
Problem.
Content of the invention
Initial mesh model is imported FLAC3D software by interface routine by the present invention, is then counted according to before stress and strain model
Each component volume content, by the Fish Programming with Pascal Language in FLAC3D software, for being mingled with material in initial mesh model
Material constituent unit quantity is compensated.Reach under conditions of the distribution of inclusion material Component Space is not changed, in model as far as possible
Inclusion material constituent unit quantity is compensated, and reduces the deviation occurred in reverseization modeling process, to follow-up Numerical-Mode
Intend the purpose that work provides more accurate numerical model.To achieve these goals, the present invention is adopted the following technical scheme that：
A kind of CT inversely changes the grid model optimization method of modeling technique, it is characterised in that comprise the steps：
(1) assume when element number enough, and when between different units, size is more or less the same, inclusion material component
The quantity of unit is equal to its volume ratio with the ratio of total element number.
(2) according to the inclusion material constituent content in the archetype being grouped into by inclusion material component and stock group
With the difference of inclusion material constituent content in grid model, calculating needs in grid model to choose P stock component
Unit compensation is in inclusion material constituent unit.
(3) in FLAC3D software, Traversal Unit id, if unit component is judged to inclusion material constituent unit and is named as
1st group, then with inclusion material constituent unit centre coordinate as the centre of sphere, a spherical determinating area is generated as W millimeter with radius, will
The stock constituent unit that coordinate center falls in this region is defined as the 3rd group, while counting its total number for Q.
(4) N=P/Q, Traversal Unit id is made, if unit component is judged to the 3rd group, a random number M0 ＜ M ＜ 1 is generated,
If M<N, then unit be newly defined as " add " group, if M>=N, then unit be newly defined as the 2nd group.Then by the 1st group and
" add " group merges, and generates inclusion material constituent unit, and is named as group1, and remaining unit is defined as stock component list
Unit, is named as group2, finally completes model optimization.
Compared with existing correlation technique, the invention has the advantages that：
The invention provides a kind of new mathematic calculation, solves and inversely changes in modeling process in CT, due to grid
The restriction of decision criteria, causes the content of inclusion material component to occur in that the problem of obvious reduction after grid component.To the greatest extent
Under conditions of amount does not change the distribution of inclusion material Component Space, model is optimized, the method can be in size of mesh opening relatively
Greatly, under conditions of number of grid is limited, more accurate mathematical calculation model is provided.
Description of the drawings
Fig. 1 is the flow chart for carrying out three-dimensional grid model foundation based on CT scan picture；
Fig. 2 is the two-dimensional representation before the model of the present invention is optimized using computational methods；
Fig. 3 is the two-dimensional representation after the model of the present invention is optimized using computational methods；
Fig. 4 is to carry out three-dimensionalreconstruction model and original three-dimensional model and initial mesh model after model optimization using the present invention
Contrast effect figure；
Fig. 5 is model components hum pattern before and after stress and strain model.
In figure:1- grid model, 2- inclusion material constituent unit, 3- stock constituent unit, 4- archetype, 5- are pressed from both sides
Miscellaneous material component, 6- stock component, 7-FLAC3D software.
Specific embodiment
In conjunction with shown in Fig. 1, Fig. 2, Fig. 3, Fig. 4, a kind of CT inversely changes the grid model optimization method of modeling technique, including net
Lattice model 1, inclusion material constituent unit 2, stock constituent unit 3, archetype 4, inclusion material component 5, stock
Component 6, FLAC3D software 7, first-selection are wherein pressed from both sides according to inclusion material component 5 and stock component 6 in archetype 4
The difference of 5 content of inclusion material component in 5 content of miscellaneous material component and grid model 1, calculating needs choosing in grid model 1
Take P stock constituent unit 3 to compensate in inclusion material constituent unit 2；Then pass through interface routine by initial mesh mould
Type 1 imports FLAC3D software 7, using the Fish Programming with Pascal Language in FLAC3D software 7, for being mingled with initial mesh model 1
2 quantity of material component unit is compensated.
The standard coal test specimen model for taking certain ore deposit diameter 50mm, height 100mm is test specimen, illustrates that the present invention is specifically walked
Rapid as follows：
(1) assume when element number enough, and when between different units, size is more or less the same, inclusion material component
The quantity of unit 2 is equal to its volume ratio with the ratio of total element number, and this has obtained preferable checking in Figure 5.
(2) according to the inclusion material component 5 in the archetype 4 being made up of inclusion material component 5 and stock component 6
The difference of 5 content of inclusion material component in content and grid model 1, calculating needs in grid model 1 to choose 26962 bases
This material component unit 3 is compensated in inclusion material constituent unit 2.
(3) in FLAC3D software 7, Traversal Unit id, if unit component is judged to inclusion material constituent unit 2 and names
For the 1st group, then with 2 centre coordinate of inclusion material constituent unit as the centre of sphere, a spherical determinating area is generated as 1mm with radius,
The stock constituent unit 3 that coordinate center is fallen in this region is defined as the 3rd group, totally 43946 units.
(4) N=26962/43947=0.6135, Traversal Unit id is made, if unit component is judged to the 3rd group, generates one
Random number M (0,1), if M is ＜ N, unit is newly defined as " add " group, if M is ＞=N, unit is newly defined as the 2nd group.So
Afterwards the 1st group and " add " group are merged, inclusion material constituent unit 2 is generated, and group1 is named as, remaining unit is defined as base
This material component unit 3, is named as group2, finally completes model optimization.
Certainly, only presently preferred embodiments of the present invention described above, the present invention are not limited to enumerate above-described embodiment, should
When explanation, any those of ordinary skill in the art are under the teaching of this specification, all equivalents for being made, bright
Aobvious variant, all falls within the essential scope of this specification, ought to be protected by the present invention.
Claims (1)
1. a kind of CT inversely changes the grid model optimization method of modeling technique, it is characterised in that comprise the steps：
(1) assume when element number enough, and when between different units, size is more or less the same, inclusion material constituent unit
The ratio of quantity and total element number be equal to the volume of inclusion material constituent unit and the ratio of the volume of total unit；
(2) according to the inclusion material constituent content in the archetype being grouped into by inclusion material component and stock group and net
The difference of inclusion material constituent content in lattice model, calculating needs in grid model to choose P stock constituent unit
Compensate in inclusion material constituent unit；P is the number for the stock constituent unit for compensating, and P is natural number；
(3) in FLAC3D software, Traversal Unit id, if unit component is judged to inclusion material constituent unit and is named as the 1st
Group, then with inclusion material constituent unit centre coordinate as the centre of sphere, generate a spherical determinating area with radius as W millimeter, will sit
The stock constituent unit that mark center falls in this region is defined as the 3rd group, while counting its total number for Q；Q is the 3rd
The number of the included stock constituent unit of group, Q are natural number；
(4) N=P/Q, N are made by number and the 3rd group of stock group for including of the stock constituent unit for compensating
The ratio of the number of subdivision；Traversal Unit id, if unit component is judged to the 3rd group, generates random number M (0,1), if M<
N, then unit be newly defined as " add " group, if M>=N, then unit be newly defined as the 2nd group.Then by the 1st group and " add " group
Merge, inclusion material constituent unit is generated, and group1 is named as, remaining unit is defined as stock constituent unit, name
For group2, model optimization is finally completed；" add " group is the supplement of the 1st group of inclusion material constituent unit；2nd group i.e. by base
The group that this material component unit is formed.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510256713.9A CN104952108B (en) | 2015-05-20 | 2015-05-20 | A kind of CT inversely changes the grid model optimization method of modeling technique |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510256713.9A CN104952108B (en) | 2015-05-20 | 2015-05-20 | A kind of CT inversely changes the grid model optimization method of modeling technique |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104952108A CN104952108A (en) | 2015-09-30 |
| CN104952108B true CN104952108B (en) | 2017-03-08 |
Family
ID=54166733
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510256713.9A Expired - Fee Related CN104952108B (en) | 2015-05-20 | 2015-05-20 | A kind of CT inversely changes the grid model optimization method of modeling technique |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104952108B (en) |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2004043444A1 (en) * | 2002-11-06 | 2004-05-27 | Mount Sinai School Of Medicine | Treatment of amyotrophic lateral sclerosis with nimesulide |
| CN102768699B (en) * | 2012-06-14 | 2016-08-24 | 西安交通大学 | Method based on CT image Accurate Reconstruction dissimilar materials microcosmic FEM mesh |
| CN103822865B (en) * | 2014-03-20 | 2016-05-04 | 中国石油大学(华东) | A kind of high-resolution three-dimension digital cores modeling method |
-
2015
- 2015-05-20 CN CN201510256713.9A patent/CN104952108B/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| CN104952108A (en) | 2015-09-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zheng et al. | GPU-based parallel algorithm for particle contact detection and its application in self-compacting concrete flow simulations | |
| CN104866652A (en) | Finite element simulation method for shot-peening strengthening deformation based on ABAQUS | |
| CN103150758A (en) | Method for highly fine three-dimensional modeling of complex earth surface | |
| CN106960070A (en) | A kind of seepage simulation method that coal body is reconstructed based on finite element discretization member CT | |
| CN107908913B (en) | Earth power digital-analog method based on parallel computer | |
| CN104008251B (en) | A kind of Optimize magnetic circult method for designing of the Magnetic Flux Leakage Testing System based on orthogonal test | |
| CN108984829B (en) | Calculation method and system for stacking process of rock-fill concrete rock-fill body | |
| CN106886685A (en) | A kind of interface loading force prediction method of meter and three-dimensional fractal | |
| Du et al. | FEM-DEM coupling analysis for solid granule medium forming new technology | |
| CN108090313A (en) | A kind of Complex Rock model of fissuration models recognition methods | |
| CN104197932B (en) | Praxeology-based position tracking method | |
| CN105136623A (en) | Potential energy change based method for quantitatively characterizing packing segregation state of particles after falling | |
| Wu et al. | Effect of aggregate morphology characteristics on the voidage of aggregate loose packing based on 3D discrete element method | |
| CN104952108B (en) | A kind of CT inversely changes the grid model optimization method of modeling technique | |
| CN104973215A (en) | Global dynamic virtual reality carrying method based on precision database | |
| CN112668087B (en) | Ballasted track modeling analysis method and system based on finite element and discrete element coupling | |
| CN112184905B (en) | Three-dimensional modeling method, device and equipment for uneven roadway | |
| CN104794745A (en) | 3D and isogeometric mixed unit modeling method of rifling barrel | |
| CN109033566B (en) | Vehicle load numerical simulation method based on finite element-infinite element model | |
| CN103678732B (en) | The method that hydraulic engineering improves water environment proper scale is differentiated based on simple curve bathmometry | |
| Wang et al. | Three-dimensional cellular automata based particle flow simulations of mechanical properties of talus deposit | |
| Xia | Analysis of Machinery Manufacture and Production from the Metal Mechanics Perspective | |
| Zhang et al. | A Design Method and Application of Backfilling Hydraulic Supports in Reforming Top-Coal Caving Hydraulic Supports | |
| Liu et al. | Three-Dimensional Visual Simulation for Concrete Dam Construction Process | |
| CN106864677A (en) | A kind of design of liquid container saddle seat, manufacture and installation method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170308 Termination date: 20180520 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |