CN115619961B - 基于矩阵表达式的三维块体切割与搜索方法 - Google Patents
基于矩阵表达式的三维块体切割与搜索方法 Download PDFInfo
- Publication number
- CN115619961B CN115619961B CN202210883050.3A CN202210883050A CN115619961B CN 115619961 B CN115619961 B CN 115619961B CN 202210883050 A CN202210883050 A CN 202210883050A CN 115619961 B CN115619961 B CN 115619961B
- Authority
- CN
- China
- Prior art keywords
- spatial
- polygons
- polygon
- matrix
- block
- 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
- 239000011159 matrix material Substances 0.000 title claims abstract description 47
- 238000000034 method Methods 0.000 title claims abstract description 14
- 239000013598 vector Substances 0.000 claims abstract description 7
- 239000011435 rock Substances 0.000 claims description 8
- 238000010586 diagram Methods 0.000 description 7
- 101001121408 Homo sapiens L-amino-acid oxidase Proteins 0.000 description 2
- 102100026388 L-amino-acid oxidase Human genes 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
Classifications
-
- 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/10—Constructive solid geometry [CSG] using solid primitives, e.g. cylinders, cubes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T19/00—Manipulating 3D models or images for computer graphics
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computer Graphics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Geometry (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Computer Hardware Design (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Processing Or Creating Images (AREA)
Abstract
本发明提供一种基于矩阵表达式的三维块体切割与搜索方法，即通过矢量描述建立初始信息关系矩阵，对所有空间多边形进行求交运算，将空间多边形进行离散，建立描述空间多边形之间所有拓扑关系的矩阵；根据拓扑关系矩阵清除所有不能形成块体边界的多边形，建立只包含能够形成块体边界的新的拓扑关系矩阵：通过搜索建立的新的拓扑关系矩阵中的空间多边形实现三维块体搜索。本发明能够实现任意形状的块体搜索，块体与边界范围具有统一的表达式，实现三维全空间的搜索。
Description
技术领域
本发明涉及一种三维块体切割与搜索方法，具体地说，涉及一种基于矩阵表达式的三维块体切割与搜索方法。本发明属于岩土工程中岩体切割与块体搜索技术领域。
背景技术
节理岩体是岩土工程经常遇到的介质，岩体内包含大量结构面(如断层、节裂隙、层面等)，这些结构面相互交叠且形成了具有特性的岩体结构，岩体结构决定了节理岩体的变形、强度、稳定等特征。因此，高效快速描述上述特征，实现节理岩体的高效切割与搜索是开展节理岩体相关研究的前提和基础。
发明内容
鉴于上述原因，本发明的目的是提供一种基于矩阵表达式的三维块体切割与搜索方法。
为实现上述目的，本发明采用以下技术方案：一种基于矩阵表达式的三维块体切割与搜索方法，它包括如下步骤：
S1：通过矢量描述建立初始信息关系矩阵，对所有空间多边形进行求交运算，将空间多边形进行离散，建立描述空间多边形之间所有拓扑关系的矩阵；
S2：根据拓扑关系矩阵清除所有不能形成块体边界的多边形，建立只包含能够形成块体边界的新的拓扑关系矩阵：
S3：通过搜索步骤S2.2建立的新的拓扑关系矩阵中的空间多边形实现三维块体搜索：
S3.1：找到新的空间多边形之间拓扑关系矩阵中第1个编号大于0的空间多边形，将该空间多边形的编号作为起始面号i，将面号乘以-1表示该面已经使用过并存储该空间多边形新的面号，同时将编号为i的多边形中所有线段号存入临时线段数组；
S3.2：找到临时线段数组中第一个线段编号k，通过拓扑关系矩阵找到与编号为i的空间多边形公用边界k的多边形，当多边形数超过两个时，按步骤S2.2中排序的下一个多边形为j，删除j中所有与临时线段数组相同的边界，将剩余边界即线段存储形成新的临时线段数组；
S3.3：重复步骤S3.2直至临时线段数组为空，一个块体搜索结束；
S3.4：重复步骤S3.1-步骤S3.3直至所有块体搜索完成。
在本发明较佳实施例中，所述步骤S1建立描述空间多边形之间所有拓扑关系矩阵的方法为：
S1.1：输入初始信息，该初始信息包括用于描述三维块体的节点和空间多边形信息，每个空间多边形由一系列逆时针排列的节点组成，为保证数据的严密性，节点数组中不包含未被空间多边形使用的节点；
S1.2：每个空间多边形之间相互求交，记录下所有交线，根据交线对空间多边形进行离散形成用于块体搜索的空间多边形，离散后的空间多边形只在多边形边界(即线段)上相交；
S1.3：以空间多边形边界线段为基础，建立描述空间多边形公用线段与空间多边形之间的拓扑关系的矩阵。
在本发明较佳实施例中，所述步骤S2建立只包含能够形成块体边界的新的拓扑关系矩阵的方法为：
S2.1：以线段与空间多边形间拓扑关系的矩阵为基础，删除仅被使用一次的线段，同时删除使用该线段的空间多边形；
S2.2：删除不能形成块体的节点和线段后，建立新的空间多边形之间的拓扑关系矩阵；当线段被超过两个空间多边形使用时，矩阵中多边形按其对应的空间向量逆时针排序。
本发明能够实现任意形状的块体切割，块体与边界范围具有统一的表达式，实现三维全空间的快速搜索。
附图说明
图1是本发明三维块体切割与搜索方法流程图；
图2A是本发明输入的初始信息节点编号示意图；
图2B是本发明输入的初始信息空间多边形编号示意图；
图3A是本发明空间多边形求交前节点编号示意图；
图3B是本发明空间多边形求交后多边形编号示意图；
图4A是本发明去除只被使用1次的节点后，重新编号的节点示意图；
图4B是本发明去除只被使用1次的节点后，重新编号的空间多边形；
图5是本发明共线段多边形根据空间位置关系排序示意图；
图6是本发明三维块体搜索全过程图。
具体实施方式
下面结合附图和实施例对本发明的结构及特征进行详细说明。需要说明的是，可以对此处公开的实施例做出各种修改，因此，说明书中公开的实施例不应该视为对本发明的限制，而仅是作为实施例的范例，其目的是使本发明的特征显而易见。
如图1所示，本发明公开的基于矩阵表达式的三维块体切割与搜索方法包括如下步骤：
S1：通过矢量描述建立初始信息关系矩阵，对所有空间多边形进行求交运算，将空间多边形进行离散，建立描述空间多边形之间所有拓扑关系的矩阵，具体步骤为：
S1.1：输入初始信息，该初始信息包括用于描述三维块体的节点和空间多边形信息，每个空间多边形由一系列逆时针排列的节点组成，为保证数据的严密性，节点数组中不包含未被空间多边形使用的节点，如图2A和图2B所示。
S1.2：每个空间多边形之间相互求交，记录下所有交线，根据交线对空间多边形进行离散形成用于块体搜索的空间多边形，离散后的空间多边形只在多边形边界(即线段)上相交，如图3A和图3B所示。
S1.3：以空间多边形边界线段为基础，建立描述空间多边形公用线段与空间多边形之间的拓扑关系的矩阵。
S2：根据拓扑关系矩阵清除所有不能形成块体边界的多边形，建立只包含能够形成块体边界的新的拓扑关系矩阵：
S2.1：以线段与空间多边形间拓扑关系的矩阵为基础，删除仅被使用一次的线段，同时删除使用该线段的空间多边形，如图4A、图4B所示。
S2.2：删除不能形成块体的节点和线段后，建立新的空间多边形之间的拓扑关系矩阵；当线段被超过两个空间多边形使用时，矩阵中多边形按其对应的空间向量逆时针排序，如图5所示。
S3：通过搜索步骤S2.2建立的新的拓扑关系矩阵中的空间多边形实现三维块体搜索：
S3.1：找到新的空间多边形之间拓扑关系矩阵中第1个编号大于0的空间多边形，将该空间多边形的编号作为起始面号i，将面号乘以-1表示该面已经使用过并存储该空间多边形新的面号，同时将编号为i的多边形中所有线段号存入临时线段数组；
S3.2：找到临时线段数组中第一个线段编号k，通过拓扑关系矩阵找到与编号为i的空间多边形公用边界k的多边形，当多边形数超过两个时，按步骤S2.2中排序的下一个多边形为j，删除j中所有与临时线段数组相同的边界，将剩余边界即线段存储形成新的临时线段数组。
S3.3：重复步骤S3.2直至临时线段数组为空，一个块体搜索结束；
S3.4：重复步骤S3.1-步骤S3.3直至所有块体搜索完成；搜索过程如图6所示。
本发明能够实现任意形状的块体切割，块体与边界范围具有统一的表达式，实现三维全空间的快速搜索。
最后应说明的是：以上所述的各实施例仅用于说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述实施例所记载的技术方案进行修改，或者对其中部分或全部技术特征进行等同替换；而这些修改或替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
Claims (1)
1.一种基于矩阵表达式的三维节理岩体块体切割与搜索方法，其特征在于：它包括如下步骤：
S1：通过矢量描述建立初始信息关系矩阵，对所有空间多边形进行求交运算，将空间多边形进行离散，建立描述空间多边形之间所有拓扑关系的矩阵：
S1.1：输入初始信息，该初始信息包括用于描述三维块体的节点和空间多边形信息，每个空间多边形由一系列逆时针排列的节点组成，为保证数据的严密性，节点数组中不包含未被空间多边形使用的节点；
S1.2：每个空间多边形之间相互求交，记录下所有交线，根据交线对空间多边形进行离散形成用于块体搜索的空间多边形，离散后的空间多边形只在多边形边界线段上相交；
S1.3：以空间多边形边界线段为基础，建立描述空间多边形公用线段与空间多边形之间的拓扑关系的矩阵；
S2：根据拓扑关系矩阵清除所有不能形成块体边界的空间多边形，建立只包含能够形成块体边界的新的拓扑关系矩阵：
S2.1：以线段与空间多边形间拓扑关系的矩阵为基础，删除仅被使用一次的线段，同时删除使用该线段的空间多边形；
S2.2：删除不能形成块体的节点和线段后，建立新的空间多边形之间的拓扑关系矩阵；当线段被超过两个空间多边形使用时，矩阵中多边形按其对应的空间向量逆时针排序；
S3：通过搜索步骤S2.2建立的新的拓扑关系矩阵中的空间多边形实现三维块体搜索：
S3.1：找到新的空间多边形之间拓扑关系矩阵中第1个编号大于0的空间多边形，将该空间多边形的编号作为起始面号i，将面号乘以-1表示该面已经使用过并存储该空间多边形新的面号，同时将编号为i的多边形中所有线段号存入临时线段数组；
S3.2：找到临时线段数组中第一个线段编号k，通过拓扑关系矩阵找到与编号为i的空间多边形公用边界k的多边形，当多边形数超过两个时，按步骤S2.2中排序的下一个多边形为j，删除j中所有与临时线段数组相同的边界，将剩余边界即线段存储形成新的临时线段数组；
S3.3：重复步骤S3.2直至临时线段数组为空，一个块体搜索结束；
S3.4：重复步骤S3.1-步骤S3.3直至所有块体搜索完成。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210883050.3A CN115619961B (zh) | 2022-07-26 | 2022-07-26 | 基于矩阵表达式的三维块体切割与搜索方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210883050.3A CN115619961B (zh) | 2022-07-26 | 2022-07-26 | 基于矩阵表达式的三维块体切割与搜索方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115619961A CN115619961A (zh) | 2023-01-17 |
| CN115619961B true CN115619961B (zh) | 2024-04-30 |
Family
ID=84857691
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210883050.3A Active CN115619961B (zh) | 2022-07-26 | 2022-07-26 | 基于矩阵表达式的三维块体切割与搜索方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115619961B (zh) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101936008A (zh) * | 2010-09-30 | 2011-01-05 | 东北大学 | 岩体边坡三维模型及块体滑落分析方法 |
| CN102609982A (zh) * | 2012-01-20 | 2012-07-25 | 北京石油化工学院 | 空间地质数据非结构化模式的拓扑发现方法 |
| CN103065319A (zh) * | 2012-12-31 | 2013-04-24 | 上海同岩土木工程科技有限公司 | 一种空间多连通域的封闭面自动搜索方法 |
| CN104715507A (zh) * | 2015-04-09 | 2015-06-17 | 武汉大学 | 一种基于曲面片的三维地理实体自动构建方法 |
| CN107301678A (zh) * | 2017-06-29 | 2017-10-27 | 中国水利水电科学研究院 | 一种地质块体的三维建模方法 |
| CN107393005A (zh) * | 2017-07-25 | 2017-11-24 | 中国科学院武汉岩土力学研究所 | 一种三维岩石块体系统建模方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101582173B (zh) * | 2009-06-24 | 2012-07-11 | 中国石油天然气集团公司 | 复杂地质构造块状模型构建方法 |
-
2022
- 2022-07-26 CN CN202210883050.3A patent/CN115619961B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101936008A (zh) * | 2010-09-30 | 2011-01-05 | 东北大学 | 岩体边坡三维模型及块体滑落分析方法 |
| CN102609982A (zh) * | 2012-01-20 | 2012-07-25 | 北京石油化工学院 | 空间地质数据非结构化模式的拓扑发现方法 |
| CN103065319A (zh) * | 2012-12-31 | 2013-04-24 | 上海同岩土木工程科技有限公司 | 一种空间多连通域的封闭面自动搜索方法 |
| CN104715507A (zh) * | 2015-04-09 | 2015-06-17 | 武汉大学 | 一种基于曲面片的三维地理实体自动构建方法 |
| CN107301678A (zh) * | 2017-06-29 | 2017-10-27 | 中国水利水电科学研究院 | 一种地质块体的三维建模方法 |
| CN107393005A (zh) * | 2017-07-25 | 2017-11-24 | 中国科学院武汉岩土力学研究所 | 一种三维岩石块体系统建模方法 |
Non-Patent Citations (2)
| Title |
|---|
| 基于自由面上环路对应节理无限切割的关键块体搜索;石露;李小春;白冰;;岩土力学;20120710(07);全文 * |
| 边坡及洞室岩体的全空间块体拓扑搜索研究;张奇华;邬爱清;;岩石力学与工程学报;20081015(10);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115619961A (zh) | 2023-01-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104679949B (zh) | 基于XML路网数据的Paramics路网构建方法 | |
| CN101499097B (zh) | 基于散列表的数据流频繁模式内存压缩存储方法 | |
| US10838825B2 (en) | Implementing snapshot sets for consistency groups of storage volumes | |
| CN105095237A (zh) | 用于生成非关系数据库的模式的方法和设备 | |
| Kejlberg-Rasmussen et al. | Faster worst case deterministic dynamic connectivity | |
| CN107678883A (zh) | 一种基于存储系统的集群恢复方法及设备 | |
| CN107798063A (zh) | 快照处理方法和快照处理装置 | |
| CN111984732B (zh) | 在区块链上实现去中心化检索的方法、节点及区块链网络 | |
| CN107609011A (zh) | 一种数据库记录的维护方法和装置 | |
| CN115619961B (zh) | 基于矩阵表达式的三维块体切割与搜索方法 | |
| CN110011830A (zh) | 基于流量数据的通讯拓扑信息建模方法 | |
| CN111506552B (zh) | 一种树状结构的动态数据库设计方法及系统 | |
| CN103530067B (zh) | 一种数据操作的方法和设备 | |
| US9732593B2 (en) | Systems, methods, and computer medium to optimize storage for hydrocarbon reservoir simulation | |
| CN106855890A (zh) | 一种实现高性能数据库最终一致性全文检索的方法 | |
| CN112488708B (zh) | 区块链账户关联性查询方法及虚假交易筛选方法 | |
| CN115619893B (zh) | 基于矩阵表达式的二维块体切割与搜索方法 | |
| CN104539298A (zh) | 一种Reed-Muller逻辑电路极性快速转换方法 | |
| Ledoux et al. | Flipping to robustly delete a vertex in a Delaunay tetrahedralization | |
| CN105302889A (zh) | 数据存储结构的转换方法及装置 | |
| CN112329312B (zh) | 一种三维渗流应力耦合内聚力单元的快速生成方法 | |
| CN104484368B (zh) | 一种基于MapReduce的大规模图数据关键词搜索方法 | |
| CN104199895B (zh) | 一种基于双体内存系统的大数据实时排序处理器 | |
| CN106227036B (zh) | 一种对称离散事件系统在线控制律重构方法 | |
| CN117312633B (zh) | 基于带hbm的fpga的动态极大团枚举装置及方法 |
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 | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20230828 Address after: 430010 No.1, Liuhe Road, Jiang'an District, Wuhan City, Hubei Province Applicant after: CHINA THREE GORGES Corp. Applicant after: China Three Gorges Construction Engineering (Group) Co.,Ltd. Applicant after: China Institute of Water Resources and Hydropower Research Address before: 430010 No.1, Liuhe Road, Jiang'an District, Wuhan City, Hubei Province Applicant before: CHINA THREE GORGES Corp. Applicant before: China Institute of Water Resources and Hydropower Research |
|
| TA01 | Transfer of patent application right | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |