CN108090313B - Complex rock fracture model modeling and identifying method - Google Patents
Complex rock fracture model modeling and identifying method Download PDFInfo
- Publication number
- CN108090313B CN108090313B CN201810111651.6A CN201810111651A CN108090313B CN 108090313 B CN108090313 B CN 108090313B CN 201810111651 A CN201810111651 A CN 201810111651A CN 108090313 B CN108090313 B CN 108090313B
- Authority
- CN
- China
- Prior art keywords
- model
- ansys
- software
- proe
- fracture
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- Software Systems (AREA)
- General Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Computer Graphics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention relates to a complex rock fracture model modeling and recognition method based on ANSYS and introduced into Flac3D by ProE. Due to the fact that ANSYS has the functions of convenient and fast high-quality grid division, the Flac3D software has very strong geotechnical engineering numerical calculation capacity, however, the solid modeling capacity of the two types of software is weak, and a complex model cannot be established. According to the invention, the excellent modeling capability of ProE is adopted at the early stage to model the complex model, the powerful grid processor of ANSYS is used for grid division at the middle stage, and the Flac3D is used for simulation calculation at the later stage, so that advantages are complementary, and advantages are made up for the disadvantages, and error-free link among ProE, ANSYS and Flac3D is realized, thereby not only avoiding that a large amount of time and energy are consumed by research technicians for researching the complex FISH language in Flac3D, but also avoiding that the modeling and grid division of the complex model in the geotechnical engineering problem cannot be realized.
Description
Technical Field
The invention relates to the field of geotechnical engineering simulation research, in particular to a complex rock fracture model modeling and identifying method which is based on ANSYS and is introduced into Flac3D by ProE.
Background
With the rapid economic development of China in recent years, the geotechnical engineering problems related to complex geological environments are endless in the process of various large-scale infrastructure construction. With the continuous progress of the geotechnical field, the problems of deeper excavation depth, larger excavation area, more complex construction process and the like are faced in foundation pit engineering, tunnel engineering and the like. For example, Mu Bo Niger gold mine in south Africa has been excavated to under 4000 meters underground, and the deepest depth of a diversion tunnel of a silk screen hydropower station in China reaches 2500 meters. In this case, the geotechnical engineering problem is more demanding on the expertise of designers and builders.
With the rapid development of computer technology, various numerical calculation simulation methods are rapidly developed, and are obviously improved in terms of calculated amount, calculated speed and calculated precision. Of these, ANSYS and Flac3D are the superior software versions. ANSYS is large-scale general finite element software developed by American ANSYS company, and can be interfaced with most Computer Aided Design (CAD) software to realize data sharing and exchange. The ANSYS program provides the function of grid division of the CAD model with convenient use and high quality. The method comprises four grid division methods: extended partitioning, shadow partitioning, free partitioning, and adaptive partitioning. The function of the free grid divider of the ANSYS program is very powerful, the complex model can be directly divided, and the trouble caused by mismatching of grids of all parts when a user divides and assembles all the parts is avoided. The self-adaptive mesh partition is that after a solid model with boundary conditions is generated, a user instructs a program to automatically generate a finite element mesh, analyze and estimate the discrete error of the mesh, then redefine the size of the mesh, and analyze and estimate the discrete error of the mesh again until the error is lower than a user-defined value or the user-defined solving times are reached. And through the interface between ANSYS and Flac3D, model information in ANSYS can be converted in the form of a file that Flac3D can import. However, ANSYS also has the defects of itself, the entity modeling is simpler, and some complex models in geotechnical engineering are difficult to establish through a pretreatment module of ANSYS.
Flac3D is an excellent piece of geotechnical computing software developed by ITASCA, a company in the United states. The display lagrangian algorithm and the hybrid-discrete partitioning technique employed by Flac3D can very accurately simulate the material's satellite failure and flow. An "explicit solution" scheme is employed. Thus, an explicit solution takes almost the same time to solve the nonlinear stress-strain relationship as a linear constitutive relationship, while an implicit solution takes longer to solve the nonlinear problem. Moreover, the rigidity matrix does not need to be stored, which means that the multi-unit structure can be solved by adopting a memory with medium capacity; simulating a large deformation problem consumes almost no more computing time than a small deformation problem because there is no stiffness matrix to be modified. The Flac3D software operation interface is simple, many commands need to be written in the built-in FISH language, and the designer without programming background has great difficulty. And Flac3D is extremely difficult to construct complex models and to partition meshes.
The ProE operating software is CAD/CAM/CAE integrated three-dimensional software under the American Parameter Technology Corporation (PTC) flag. Pro/Engineer is accepted and popularized by the industry as a new standard in the field of current world mechanical CAD/CAE/CAM, is one of the mainstream CAD/CAM/CAE software nowadays, and particularly occupies an important position in the field of domestic product design. Due to the fact that the mechanical product model is more elaborate and complex compared to the geotechnical model. Therefore, through the powerful modeling function of the ProE, complex three-dimensional models in geotechnical engineering can be established, and the ProE and ANSYS have a perfect interface and can be directly imported into ANSYS software for grid division.
Therefore, it is necessary to combine the advantages of the three types of software and to eliminate the disadvantages of each software. The method realizes the establishment of the complex ProE complex geometric model at the early stage, the division of ANSYS complex grids at the middle stage and the analysis and calculation of Flac3D at the later stage, and solves the problem of the establishment and calculation of the Flac3D geotechnical engineering complex model for engineering technicians.
Disclosure of Invention
In order to achieve the purpose, the invention provides a complex model identification method for importing a ProE into Flac3D based on ANSYS, realizes a method for carrying out simulation calculation by importing a ProE three-dimensional model into Flac3D, and solves the problem of building calculation of a Flac3D geotechnical engineering complex model.
The technical scheme adopted by the invention is as follows: a complex rock fracture model modeling and identifying method comprises the following specific steps:
s1, modeling the complex rock fracture model by utilizing ProE software
(1) Extracting key part control point parameters according to a complex three-dimensional model of geotechnical engineering to be established, determining the long axis, the short axis and the thickness of the long, wide, high and elliptical fractures of the complex three-dimensional model, drawing a two-dimensional view of the model by using a draft function of ProE, and then obtaining the three-dimensional model by using a part module through stretching operation;
(2) introducing a rock three-dimensional model in the component-design of the ProE, and adding two positioning planes according to the coordinates of the key position control points;
(3) introducing an elliptical three-dimensional fracture, constraining the elliptical fracture model and a positioning plane through a constraint function, further fixing the position of the elliptical fracture model in a rock body, and radially carrying out radial operation on the elliptical fracture model through a radial function to obtain a three-dimensional solid model containing multiple fractures;
(4) after the three-dimensional entity model is built, converting node information and unit information in the model into field forms through an IGES standard, and outputting the fields meeting the IGES standard into an igs format file, namely outputting the model into the igs format file;
s2, carrying out gridding division on the model established by ProE software by using ANSYS software
(1) Importing an igs format file output by a three-dimensional model established by ProE software through an Import module Import function in ANSYS, and further importing node information and unit information of the model into an ANSYS preprocessing module; the ANASYS preprocessing module identifies field information in the igs file output by the ProE file, and further generates a three-dimensional model in ANSYS;
(2) carrying out meshing on the rock three-dimensional model by using a meshing tool of ANSYS, wherein the meshing tool influences the meshing mode, and divides the rock three-dimensional model into regular hexahedral meshes;
(3) respectively carrying out module marking on the rock mass and the elliptic fracture in the model, giving different material numbers to different modules, and carrying out data conversion on the material numbers to obtain a conversion file which is recognizable by Flac3D software with a suffix of lis;
s3, calculating by using a Flac3D software line model
Importing a model into Flac3D according to a path through an external grid input interface (Import from ANSYS) under a File menu (File) in Flac3D by a File with a suffix name of. lis output by ANSYS software, displaying a rock mass and an elliptical fracture grid, endowing material properties, boundary conditions and load conditions to each grid set through a simple FISH language, setting corresponding working condition conditions, and then calculating.
Further, in the modeling of the three-dimensional solid model by using the proE software, the model which needs several components to be assembled in combination is assembled by an assembly module in the ProE.
Further, the model established by the ProE software by utilizing ANSYS software is subjected to module marking on the rock mass and the elliptic fracture in the grid division, and the density attributes of the rock mass and the elliptic fracture are assigned by giving a module to a material in ANSYS without giving an actual value. By this material-giving method, the different modules are well distinguished for better property-giving in FLAC 3D.
In the invention, as ANSYS has a convenient and fast high-quality meshing function, the Flac3D software has very strong geotechnical engineering numerical calculation capability, but the solid modeling capability of the two types of software is weak, some complex models in the aspect of geotechnical engineering cannot be established, and the meshing function of the Flac3D software is single. Therefore, the complex model is modeled by adopting the excellent modeling capability of ProE in the early stage, the powerful grid processor of ANSYS is used for grid division in the middle stage, and the Flac3D is used for simulation calculation in the later stage, so that advantages are complemented, and the advantages are made up for the deficiencies.
Compared with the prior art, the invention has the beneficial effects that: the three-dimensional entity model is modeled by using the ProE software, so that the positions of fractures in rocks, the dip angles of the fractures, the rock bridge angles and other different combination modes are conveniently positioned, the modeling complexity of ANASYS or other software is greatly simplified, the model is not required to be cut by Boolean operation in ANSYS software, and the efficiency is improved; by the complex model establishing and identifying method, error-free link among ProE, ANSYS and Flac3D can be realized, so that the problem that a research technician consumes a large amount of time and energy for researching the complex FISH language in Flac3D is avoided, and the problem that the complex model in the geotechnical engineering problem cannot be modeled and the grid division is avoided.
Drawings
In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings needed to be used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments described in the present application, and other drawings can be obtained by those skilled in the art without creative efforts.
FIG. 1 is a flow chart of the present invention;
FIG. 2 is a model diagram of the ProE software for building a three-dimensional solid model;
FIG. 3 is a model diagram of ANSYS software gridding a model established by ProE software;
fig. 4 is a model diagram introduced into Flac 3D.
Detailed Description
The invention is further illustrated by the following examples.
Example 1
In the research of the project, a numerical simulation experiment of prefabricating complex multi-fracture in rock is developed, and the flow is shown in figure 1:
s1, modeling the complex rock fracture model by utilizing ProE software
(1) Determining the length, width and height of a rock model to be 60mm, 60mm and 60mm respectively, setting the major axis of an elliptical fracture to be 10mm and the minor axis to be 5mm, setting the thickness of the fracture to be 1mm, setting the included angle between the fracture and the horizontal direction to be 45 degrees, arranging 27 fractures in the cube model, and setting the distance between the fractures to be 15 mm; drawing a two-dimensional view of the model by using a sketch function of ProE, and then obtaining a three-dimensional model by using a part module through stretching operation;
(2) introducing a rock three-dimensional model in the component-design of the ProE, and adding two positioning planes according to the coordinates of the key position control points;
(3) then, introducing an elliptical three-dimensional fracture, constraining the elliptical fracture model and the positioning plane through a constraint function, and further fixing the position of the elliptical fracture model in the rock mass, wherein the established three-dimensional solid model is shown in figure 2;
(4) the model established by ProE software is converted into a field form by passing through an IGES standard, the node information and the unit information in the model are output as an igs format file, and the import module of ANSYS imports the node information and the unit information of the model into the ANSYS preprocessing module by identifying the fields under the IGES standard.
S2, carrying out gridding division on the model established by ProE software by using ANSYS software
(1) The ANASYS preprocessing module identifies field information in the igs file output by the ProE file, and further generates a three-dimensional model in ANSYS, as shown in FIG. 3;
(2) carrying out meshing on the rock three-dimensional model by using a meshing tool of ANSYS, wherein the meshing tool influences the meshing mode, and divides the rock three-dimensional model into regular hexahedral meshes;
(3) respectively carrying out module marking on the rock mass and the elliptical fracture in the model: assigning values to density attributes of the rock mass and the fracture through a material assigning module in ANSYS, and only distinguishing 1 and 2 without assigning actual values; giving different material numbers to different modules, and performing data conversion on the material numbers to obtain a conversion file identifiable by Flac3D software with a suffix of lis format;
s3, calculating by using a Flac3D software line model
Importing a model into Flac3D by a path through a File with a suffix name of. lis output by ANSYS software through an external grid input interface (Import from ANSYS) under a File menu (File) in Flac3D, and displaying a rock mass and an elliptic fracture grid as shown in FIG. 4; and (3) endowing material attributes, boundary conditions and load conditions to each grid set through a simple FISH language, setting corresponding working condition conditions, and then calculating. Through calculation, the text modeling method is adopted, the modeling steps are simplified, the modeling time is saved, and the text meshing and grouping method is adopted, so that the meshing efficiency is improved.
The above description is only a preferred embodiment of the present invention, but the scope of the present invention is not limited thereto, and the structure, connection, measurement parameters and the like of the components may be changed, and any person skilled in the art can substitute or change the technical solution and the inventive concept of the present invention within the technical scope of the present invention.
Claims (3)
1. A complex rock fracture model modeling and identifying method is characterized by comprising the following specific steps:
s1, modeling the complex rock fracture model by utilizing ProE software
(1) Extracting key part control point parameters according to a complex three-dimensional model of geotechnical engineering to be established, determining the long axis, the short axis and the thickness of the long, wide, high and elliptical fractures of the complex three-dimensional model, drawing a two-dimensional view of the model by using a draft function of ProE, and then obtaining the three-dimensional model by using a part module through stretching operation;
(2) introducing a rock three-dimensional model in the component-design of the ProE, and adding two positioning planes according to the coordinates of the key position control points;
(3) introducing an elliptical three-dimensional fracture, constraining the elliptical fracture model and a positioning plane through a constraint function, further fixing the position of the elliptical fracture model in a rock body, and radially carrying out radial operation on the elliptical fracture model through a radial function to obtain a three-dimensional solid model containing multiple fractures;
(4) after the three-dimensional entity model is built, converting node information and unit information in the model into field forms through an IGES standard, and outputting the fields meeting the IGES standard into an igs format file, namely outputting the model into the igs format file;
s2, carrying out gridding division on the model established by ProE software by using ANSYS software
(1) Importing an igs format file output by a three-dimensional model established by ProE software through an Import module imort function in ANSYS, and further importing node information and unit information of the model into an ANSYS preprocessing module; the ANASYS preprocessing module identifies field information in the igs file output by the ProE file, and further generates a three-dimensional model in ANSYS;
(2) carrying out meshing on the rock three-dimensional model by using a meshing tool of ANSYS, wherein the meshing tool influences the meshing mode, and divides the rock three-dimensional model into regular hexahedral meshes;
(3) respectively carrying out module marking on the rock mass and the elliptic fracture in the model, giving different material numbers to different modules, and carrying out data conversion on the material numbers to obtain a conversion file which is recognizable by Flac3D software with a suffix of lis;
s3, calculating by using a Flac3D software line model
Importing a model into Flac3D according to a path by a File with a suffix name of. lis output by ANSYS software through an external grid input interface Import from ANSYS under a File menu File in Flac3D, displaying a rock mass and an elliptic fracture grid, endowing material properties, boundary conditions and load conditions to each grid set through FISH language, setting corresponding working condition, and then calculating.
2. The modeling and identifying method for the complex rock fracture model according to claim 1, wherein in the modeling of the three-dimensional solid model by using the proE software, the model which requires the combination and assembly of several components is combined and assembled by an assembly module in the proE.
3. The complex rock fracture model modeling and identification method as claimed in claim 1, wherein the model established by the ProE software through ANSYS software is subjected to module marking on the rock mass and the elliptical fracture respectively, and the density attributes of the rock mass and the fracture are assigned through a material assignment module in ANSYS without assigning actual values.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810111651.6A CN108090313B (en) | 2018-02-05 | 2018-02-05 | Complex rock fracture model modeling and identifying method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810111651.6A CN108090313B (en) | 2018-02-05 | 2018-02-05 | Complex rock fracture model modeling and identifying method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108090313A CN108090313A (en) | 2018-05-29 |
| CN108090313B true CN108090313B (en) | 2021-03-19 |
Family
ID=62193750
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810111651.6A Active CN108090313B (en) | 2018-02-05 | 2018-02-05 | Complex rock fracture model modeling and identifying method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108090313B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109242965B (en) * | 2018-07-02 | 2019-08-30 | 中国科学院武汉岩土力学研究所 | A kind of model conversion method that 3DEC is imported by Midas based on c# |
| CN110322562A (en) * | 2019-07-02 | 2019-10-11 | 中国矿业大学(北京) | A kind of discrete element 3DEC model reconstruction method based on X-CT |
| CN110779795B (en) * | 2019-11-04 | 2022-05-10 | 中国地质大学（武汉） | Method for determining size of geomechanical modeling grid unit of fractured reservoir |
| CN113128095B (en) * | 2021-04-27 | 2022-04-22 | 中国科学院武汉岩土力学研究所 | Method and device for constructing three-dimensional grid model |
| CN113239600B (en) * | 2021-07-09 | 2021-09-17 | 成都理工大学 | Method for constructing two-dimensional random network model of complex rock mass |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107180135A (en) * | 2017-05-25 | 2017-09-19 | 西南交通大学 | A kind of cave type karst rock mass numerical model construction method |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080101770A1 (en) * | 2003-09-04 | 2008-05-01 | Digital Networks North America, Inc. | Method and Apparatus for Remotely Controlling a Receiver According to Content and User Selection |
| CN102184285B (en) * | 2011-05-03 | 2013-01-02 | 三一集团有限公司 | Digital prototyping design integrated development method and system |
| CN102750407B (en) * | 2011-11-16 | 2015-04-15 | 中国电子科技集团公司第三十八研究所 | Automatic processing and converting method for computer-aided design (CAD) /computer-aided engineering (CAE) data |
| CN103699737A (en) * | 2013-12-20 | 2014-04-02 | 北京瑞风协同科技股份有限公司 | Multi-discipline simulation software data exchange method |
| CN104361211B (en) * | 2014-10-24 | 2018-07-27 | 中冶长天国际工程有限责任公司 | The structure and methods for using them of Statistical Damage Constitutive Model For Rock |
| CN105927217B (en) * | 2016-04-29 | 2019-06-25 | 中国矿业大学 | A kind of nearly total rock up-protective layer mining Design method in seam mining |
| CN106227930A (en) * | 2016-07-18 | 2016-12-14 | 东北大学 | A kind of method of model identification being imported Flac3D by Midas based on Matlab |
-
2018
- 2018-02-05 CN CN201810111651.6A patent/CN108090313B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107180135A (en) * | 2017-05-25 | 2017-09-19 | 西南交通大学 | A kind of cave type karst rock mass numerical model construction method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108090313A (en) | 2018-05-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108090313B (en) | Complex rock fracture model modeling and identifying method | |
| CN1885295B (en) | Building integrated circuits using logical units | |
| US6823501B1 (en) | Method of generating the padring layout design using automation | |
| US7117469B1 (en) | Method of optimizing placement and routing of edge logic in padring layout design | |
| CN106227930A (en) | A kind of method of model identification being imported Flac3D by Midas based on Matlab | |
| US10180996B2 (en) | Multi-component computational fluid dynamics simulations | |
| CN101826016A (en) | Visual modeling and code skeleton generating method for supporting design of multinuclear parallel program | |
| CN107145651B (en) | ABAQUS three-dimensional infinite element boundary rapid modeling method based on INP file | |
| CN109858161A (en) | A kind of Abaqus Meshing Method based on Midas modeling and Matlab conversion | |
| CN103164228B (en) | A kind of generation method of field programmable gate array program and device | |
| CN109840878A (en) | It is a kind of based on SystemC towards GPU parameter management method | |
| US20100275168A1 (en) | Design method of semiconductor integrated circuit device and program | |
| CN112270024A (en) | BIM technology-based tunnel rapid modeling method | |
| CN115437528A (en) | Chemical process simulator based on digital twins and construction method thereof | |
| CN107644139B (en) | Attribute mapping method from CAD model to CAE model | |
| CN115984503B (en) | Geological profile generation method, system, electronic equipment and medium | |
| CN113011072A (en) | Discrete element complex model identification method based on MIDAS-PFC3D | |
| CN105631135A (en) | Rapid modeling analysis system based on ABAQUS foundation pit excavation and application system thereof | |
| Essawy et al. | Elemental graph data model: a semantic and topological representation of building elements | |
| CN106354913A (en) | Parametric design platform for helical gear drive system and application method thereof | |
| CN110704916A (en) | BIM technology-based large complex three-dimensional geological model grid coarsening method | |
| Mao et al. | Modeling method of FLAC3D based on RHINO-KUBRIX and deformation analysis of tunnel lining structure | |
| CN116541938B (en) | Canal engineering building modeling method, system, electronic equipment and medium | |
| CN112528515B (en) | FLAC based on hybrid programming technology 3D Complex model identification method | |
| CN114595610B (en) | Method for creating tunnel finite element model based on Revit secondary development and Ansys command stream |
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 |