面向三维表面网格模型的有限元网格编码与索引方法技术

本发明专利技术涉及一种面向三维网格模型的有限元网格编码与索引方法，本方法能够对内存或外存中存储的体数据表达的模型进行准确的输出，提供对于给定编码规则的一种简单快速的索引方式。其操作步骤为：先根据内存或外存中的体素信息获取体数据中每个体素上的点的位置信息，再根据其位置信息计算其对应的编码。经过实例验证，此方法是一个有效的方法。

【技术实现步骤摘要】

本专利技术涉及一种面向三维表面网格模型的有限元网格编码与其索引方法，是一种基于三维八叉树体素模型的有限元网格编码与其索引方法，并能够进行应内存外运算。
技术介绍
在对模型进行有限元分析时需要对求解域进行离散化，求解域的离散化是有限元的核心技术之一，对离散化的网格进行编码需要有方便快速的索引算法。根据单元的维数特征可以分为一维、二维、三维单元，本专利技术涉及的是对三维体素单元进行的编码与索引。对应于不同的分析需求，需要针对不同的情况将CAD模型转化为不同的网格模型，对于三维模型可以是四面体、五面体或六面体。每一个网格的每一个顶点需要有一个编码，并且每一个顶点会对应若干个网格，需要进行唯一的定义。八叉树是一种用于描述三维空间的树状结构，数中每一个节点对应八个子节点，通过定义每个节点是否为空能够表达出对应的六面体体素模型，体素模型结构清晰，对模型进行处理比较方便，是常用的一种模型表达方式。体素模型所在空间为一个立方体包围盒，包围盒中按照八叉树的结构被划分成了众多的小立方体，根据标记存在与否表达图形。对六面体体素进行有限元分析需要进行网格编码，根据分析需求有时需要输出其他诸如四面体的网格。使用八叉树进行高精度的模型表达时，模型的信息量大，通常都需要信息存储在外存中，在使用时从外存中读取信息进行模型重建。外存中需要存储位置、领域关系等信息，本专利技术中的有限元网格生成方法在使用外存存储模型时，需要从外存中读取位置信息。r>张志明(CN103678792A)设计了预应力罐体结构有限元模型网格划分方法，能够根据不同的罐体壁厚生成不同的网格划分。洪军(CN102332047A)设计了一种可实现六面体网格划分的螺栓有限元参数化建模方法。本专利技术给出的是一种根据已有体素模型进行编码与索引的方式。
技术实现思路
本专利技术的目的在于针对已有技术的不足，提供一种面向三维表面网网格模型的有限元网格编码与索引方法。索引方法简单，能够根据此方法生成六面体网格，也可以根据需求输出四面体、五面体的网格。能够避免对顶点定义的重复，有效避免几何错误。为了达到上述目的，本专利技术的构思是：从内存或外存中读入每个体素的位置信息，以此计算出其上每个顶点的位置信息，再根据每个点的位置信息计算出相应的编码，最终建立空间坐标到一维编码的映射关系，将有限元网格转化为序列输出。其中包括两步映射，首先是根据体素的位置与点在体素上的位置到点在全局的位置的映射，再是点的位置到一维编码的映射。一.空间上点的编码定义为了将空间坐标映射到一维编码，需要对空间上的点进行编码。将立方体包围盒的一个顶点定为直角坐标系原点，三条坐标轴分别与此顶点相邻的三条边重合。坐标原点的节点编码为0，首先沿着Y轴正方向依次递增编号，排列完了一列后从第二列第一个节点开始继续递增编号，一层完了后从第二层继续编号，直到所有节点都被编号。如附图1所示。二.节点在体素上的编号的定义每个体素上有八个点，需要对这八个点定义一个在体素上的局部编号。如附图2所示，将体素平移至坐标原点，原点上的编号为4，位于原点Y轴正向的编号为5，位于原点X轴正向的编号为6，XOY平面上剩余的一个为7；偏向于Z轴正方向的四个点分别为XOY上四个点的编号减4。三.节点全局位置信息的获取读取体素的位置信息，每个体素具有在三个坐标轴方向上的编号，节点的全局位置信息通过体素的位置信息与节点在体素上的上的相对位置获得。若点位于Z轴的正方向一侧，则点在Z轴方向上的编号为体素在Z轴方向上的编号加1，若位于负方向一侧则点在Z轴方向上的编号为体素在Z轴方向上的编号；节点在X,Y轴上的编号也用同样方法获得。四.节点编码的计算获得了节点的全局位置信息后(也就是分别在三个坐标轴的方向上是第几个节点)，可以根据简单的方法计算出全局的编码：Cp＝x·q+y+z*q*q(1)Cp——点的全局编码；x,y,z——在X,Y,Z轴的方向上分别是第几个节点；q——每一列上的节点数；其中q＝2n-1。n——划分次数。根据上述专利技术构思，本专利技术采用如下述技术方案：一种针对表面网格模型的有限元划分编码与索引方法，其特征特在于操作步骤如下：1)节点全局位置信息的获取；2)节点编码的计算。所述步骤1)节点全局信息的获取：读取体素的信息，全局信息从体素的位置信息与节点相对于体素的位置获取，根据空间上点的编码的定义与节点在体素上的编号的定义，节点的全局信息获取可以表达为(其中x,y,z为节点分别为X,Y,Z方向上的第几个节点，zv,yv,xv为所在体素分别为X,Y,Z方向上的第几个体素)：若点在体素上位于Z轴的正方向一侧，则z＝zv+1若位于负方向一侧z＝zv若点在体素上位于Y轴的正方向一侧，则y＝yv+1若位于负方向一侧y＝yv若点在体素上位于X轴的正方向一侧，则x＝xv+1若位于负方向一侧x＝xv；所述步骤2)节点编码的计算：根据定义的全局编码规则与获取的节点信息，根据式1能够方便计算出节点的编码。映射关系一共有两步(如附图3所示)，两步映射可以用一个分段函数表达出来： C p = ( x - 1 ) * q + y + z * q * q ( b = 0 ) ( x - 1 ) * q + y + 1 本文档来自技高网...

【技术保护点】
一种面向三维网格模型的有限元编码与索引方法，其特征在于操作步骤如下：1)节点全局位置信息的获取；2）节点编码的计算。

【技术特征摘要】
1.一种面向三维网格模型的有限元编码与索引方法，其特征在于操作步骤如下：1)节
点全局位置信息的获取；2）节点编码的计算。
2.根据权利要求1所述的面向三维网格模型的有限元网格编码与索引方法，其特征在于：
所述步骤1）节点全局信息的获取：读取体素的位置信息，全局信息从体素的位置信息与节
点相对于体素的位置获取，根据空间上点的编码的定义与节点在体素上的编号的定义，节
点的全局信息获取可以表达为（其中x,y,z为节点分别为X,Y,Z方向上的第几个体素，,
,为所在体素分别为X,Y,Z方向上的第几个体素）：
若点在体素上位于Z轴的正方向一侧，则z=+1
若位于负方向一侧z=若点...

【专利技术属性】
技术研发人员：沈隽晟，李婉婉，姚远，
申请(专利权)人：上海大学，
类型：发明
国别省市：上海;31