本发明涉及计算机程序应用领域，尤其涉及一种基于bim技术的双侧壁导坑法隧道开挖单元划分方法。
背景技术：
在建筑信息管理bim技术中，通过对工程项目各项相关信息数据的采集、整理能够对工程项目进行总体的管理，在bim工程中，现有的隧道施工模型都是通过从设计模型或是手工进行划分若干的单元再进行施工组织的，这样的施工组织粗糙；目前的施工模型是直接遵从设计模型进行的，施工信息与设计模型没有有效地对接；而且在隧道进行开挖模拟过程中，不能有效的形成快捷挖掘方案。这就亟需本领域技术人员解决相应的技术问题。
技术实现要素：
本发明旨在至少解决现有技术中存在的技术问题，特别创新地提出了一种基于bim技术的双侧壁导坑法隧道开挖单元划分方法。
为了实现本发明的上述目的，本发明提供了一种基于bim技术的双侧壁导坑法隧道开挖单元划分方法，包括如下步骤：
s1，对隧道模型进行设计导入，配置隧道参数属性，按照隧道开挖工法的特点，建立开挖截面的约束条件，生成开挖截面轮廓，根据隧道中心线生成出开挖面，再与工程设计结果模型进行布尔运算，获得开挖结果模型；
s2，根据这些分层模型生成拆分单元，对这些分层模型进行单元属性配置和拓展，获得该开挖模型施工的bim模型。
所述的基于bim技术的双侧壁导坑法隧道开挖单元划分方法，优选的，所述s1包括：
s1-1，根据选择的隧道截面，计算截面的法向量vsecnormal；
s1-2，根据隧道中心线和隧道截面上任意一点a，计算该点在隧道中心线上的投影a1，从而获取隧道中心线在a1点的法平面和切线方向；再结合截面的法向量和a1点的切线方向向量，通过计算其夹角的余弦值确定a1的切线方向向量是否与预期反向，最终计算出隧道中心线在a1法平面的法向量vpronormal；
s1-3，规定以隧道基准开挖线方向为h方向，计算h方向向量，再与a1的切线方向求得隧道开挖的v方向，计算公式为v＝vpronormal*h，即向量v为投影面法向量vpronormal与方向向量h的叉乘；
s1-4,根据选择的隧道截面轮廓，获取施工工程数据中对象模型截面的最小包围盒l和包围盒长轴最大点maxp以及最小点minp；
s1-5，根据求得的包围盒长轴maxp和minp，结合投影面计算出截面中心点minpos。
所述的基于bim技术的双侧壁导坑法隧道开挖单元划分方法，优选的，所述s2包括：
s2-1，对每个开挖结果实体进行属性配置，包括：颜色、透明度、施工类型、人员、时间进行赋值；
s2-2，根据h和v方向和包围盒长轴长度，分别计算出隧道截面的宽和高，计算公式为：
包围盒长轴在h方向的投影向量：
vhprolen＝((((maxp–minp)％h)/(h.x*h.x+h.y*h.y+h.z*h.z))*h)；
隧道截面宽度：
hlen＝sqrt(vhprolen.x*vhprolen.x+vhprolen.y*vhprolen.y+vhprolen.z*vhprolen.z)；
其中h为s1-3求得的方向向量，h.x/h.y/h.z分别为向量h的x/y/z坐标，vhprolen.x/vhprolen.y/vhprolen.z分别为向量vhprolen的x/y/z坐标，sqrt为求平方根；
s2-3，根据输入的截面参数结合，计算出隧道截面开挖轮廓；通过扫略的方式沿着隧道中心线将隧道截面开挖轮廓扫略成开挖面，再将开挖面与开挖体按照先左右侧导坑后上下的顺序进行布尔运算，获得开挖结果模型。
综上所述，由于采用了上述技术方案，本发明的有益效果是：
本发明解决如何从设计院获得的设计模型转换为施工单位方便使用的高精度施工模型，即通过行业既定的规则或是约束条件，并对bim属性信息进行管理，自动划分便于实施的施工组织单元，实现项目施工组织的bim模型自动生成，提高工程项目施工管理的精度和准确性，通过本方法基于bim技术的双侧壁导坑法隧道开挖单元划分方法，提高项目执行效率。
本发明的附加方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解，其中：
图1是本发明工作流程图；
图2是本发明实施例示意图；
图3是本发明划分实施例示意图；
图4是本发明方形隧道划分实施例示意图；
图5是本发明菱形隧道划分实施例示意图；
图6是本发明总体示意图。
具体实施方式
下面详细描述本发明的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，仅用于解释本发明，而不能理解为对本发明的限制。
在本发明的描述中，需要理解的是，术语“纵向”、“横向”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。
其中bim(buildinginformationmodeling)技术为建筑信息建模。
其中在隧道数据模型导入到bim系统时，需要执行如下步骤：
ⅰ，直接读取工程设计模型，并对其模型数据组织结构进行转换解析；根据施工工程建设要求，形成设计施工工程基本数据信息，根据建设要求，读取导入不同格式设计结果模型；根据不同设计工具设计结果模型的方式，转换解析设计结果模型数据组织结构productstructure；按照施工特点及定位，直接读取各个不同模型设计工具设计结果模型，不依赖于不同模型设计工具。
ⅱ，按照模型数据组织结构特点，对每一个模型数据组织结构进行分类，根据分类，获取其每个分类的数据，其中包含模型组织结构上的图形数据，模型组织结构上的属性数据；根据施工工程模型数据组织结构要求，按照ⅰ的转换解析后的模型数据组织结构productstructure特点，对每一个模型数据组织结构进行分类，分类形式为装配product，零件part，体body；根据施工工程模型数据组织结构要求，对分类后的模型数据组织结构，分别获取其属性数据和图形数据；施工工程模型数据组织结构要求：装配product的子节点只能是装配product以及零件part节点，零件part的子节点只能是体body，其中装配product和零件part节点均仅代表施工工程模型树组织结构，而体body代表该施工工程模型数据下的一个几何图形数据；
ⅲ，根据获取到的模型数据组织结构，以及需要的新的模型数据组织结构，组织出新的需要的模型组织结构；
ⅳ，根据获取到的模型数据组织结构，以及对模型数据组织结构分类后获取到的对应模型属性数据和模型图形数据，按照直接加载和快速加载两种模式，分别对模型几何拓扑体进行相应处理；
根据施工工程模型数据组织结构要求，以及对设计结果模型数据组织结构分类后分别获取的对应模型属性数据和模型图形数据，针对施工工程模型数据导入的两种方式，直接加载和快速加载，分别处理；
根据施工工程模型数据要求，直接加载方式要求加载设计结果模型的几何实体数据，设计结果模型的几何实体数据由模型几何model和构造几何construction两部分组成，其中模型几何是指拓扑类指向的类，不包括具体的形状信息，构造几何是指模型几何中的声明，包含实际的形状信息；模型几何又称模型的拓扑结构，构造几何又称模型的解释结构，一个正方体的几何实体表示，1：其中模型几何拓扑是指，这个正方体几何中有多少块(正方体只有一个块)，在块中有多少面(正方体有六个面)，面上有多少边界边(12个边界边)，边界边上有多少点(八个点)，点和边之间怎样连接，2：其中构造几何指，具体的面，面上具体的边(由哪些点构成)，具体的边，具体的点坐标，构造几何是具体解释模型几何拓扑。模型实体数据中的模型几何model与构造几何construction的区别在与，前者不包含具体的形状数据信息，而后者包含实际的形状数据信息。
根据施工工程模型数据要求，直接加载设计结果模型数据需要进行精度和模型空间控制，加载设计结果模型数据的精度和模型空间控制公式如下：
设定a:绝对最小值(10e-6)
b:规格化最小值(10e-10)
c:曲线曲面的逼近精度(10e-3)
d:被认为是零的最大值(10e-11)
模型空间计算算法：
modelspace＝a/b＝10e-6/10e-10＝10e4；
根据施工工程模型数据要求，快速加载方式要求快速加载设计结果模型的图形数据，设计结果模型的图形数据是指，不包含有拓扑的纯图形数据，图形数据表示方法如下：
设定点：p1，p2，p3，p4，p5，p6
其中p1＝p4，p3＝p6，
点列表：plist＝{p1x，p1y，p1z，p2x，p2y，p2z，p3x，p3y，p3z，p4x，p4y，p4z，p5x，p5y，p5z，p6x，p6y，p6z}
面列表数据集：flist＝{pn1，p1，p2，p3，pn2，p5，p4，p6}
面face＝plist+flist
pn：点数量
的原始离散图形数据，表示图形为图3：
该图形，为一个矩形，离散后变为两个三角形，遵守右手定则，法向全部向外，优化过程为：
点去重：plist＝{p1x，p1y，p1z，p2x，p2y，p2z，p3x，p3y，p3z，p5x，p5y，p5z}，
面去重数据集：flist＝{pn1，p1，p2，p3，pn2，p5，p1，p3}
降低面片数量：例如两个四边形构成一个大的四边形，优化前为四个三角网格，优化后为两个三角网格，三角网格是计算机硬件渲染是需要的数据，三角网格数量越低，可提升渲染效率。
所述的基于bim技术的工程设计模型导入方法，优选的，所述s5包括：
根据施工工程模型数据要求，将s4中获取到的已经过处理的图形数据，以及属性数据，分别绑定到s3中重组的施工工程模型数据组织结构的对应节点上。
ⅴ，根据直接加载和快速加载两种模式，分别向模型数据组织结构各个分类中添加处理后的模型图形信息。
如图1和2所示，本发明提供了一种基于bim技术的双侧壁导坑法隧道开挖单元划分方法，包括如下步骤：
s1，对隧道模型进行设计导入，配置隧道参数属性，根据用户对挖掘隧道的需求获取相应的隧道参数，按照隧道开挖工法的特点进行预处理了，建立开挖截面的约束条件，生成拆分隧道模型轮廓，对开挖的拆分截面轮廓进行扫描，根据隧道中心线生成出开挖面，再与工程设计结果模型进行布尔运算，获得开挖的拆分隧道模型的结果；
s2，根据这些分层模型生成拆分单元，对这些分层模型进行单元属性配置和拓展，获得该开挖模型施工的bim模型。
优选的，所述s1包括：
s1-1，根据选择的隧道截面，计算截面的法向量vsecnormal；
s1-2，根据隧道中心线和隧道截面上任意一点a，计算该点在隧道中心线上的投影a1,从而获取隧道中心线在a1点的法平面和切线方向；再结合截面的法向量和a1点的切线方向向量，通过计算其夹角的余弦值确定a1的切线方向向量是否与预期反向,最终计算出隧道中心线在a1法平面的法向量vpronormal；
s1-3，规定以隧道基准开挖线方向为横轴h方向进行拆分，计算h方向向量，再与a1的切线方向求得隧道开挖的纵轴v方向进行拆分，计算公式为v＝vpronormal*h，即向量v为投影面法向量vpronormal与方向向量h的叉乘；
s1-4,根据选择的隧道截面轮廓，获取施工工程数据中对象模型截面的最小包围盒l和包围盒长轴最大点maxp以及包围盒长轴最小点minp；
s1-5，根据求得的包围盒长轴最大点maxp以及包围盒长轴最小点minp，结合投影面计算出截面中心点minpos。
如图3所示，优选的，所述s2包括：
s2-1，对每个开挖结果实体进行属性配置，包括：颜色、透明度、施工类型、人员、时间进行赋值；
s2-2，根据h和v方向和包围盒长轴长度，分别计算出隧道截面的宽和高，计算公式为：
包围盒长轴在横轴h方向的投影向量：
vhprolen＝((((maxp–minp)％h)/(h.x*h.x+h.y*h.y+h.z*h.z))*h)；
隧道截面宽度：
hlen＝sqrt(vhprolen.x*vhprolen.x+vhprolen.y*vhprolen.y+vhprolen.z*vhprolen.z)；
其中h为s1-3求得的方向向量，h.x/h.y/h.z分别为向量h的x/y/z坐标，vhprolen.x/vhprolen.y/vhprolen.z分别为向量vhprolen的x/y/z坐标，sqrt为求平方根；根据图3、图4和图5的具体实施例进行分析来看，首先采集竖向数据h，然后采集横向数据w，从而对于数据的采集能够有序进行，使工程数据更加富有调理，当然，根据工程施工的具体需要，针对其他实际情况采用先采集横向数据w，然后采集竖向数据h的方法，也能够起到很好的技术效果，实际操作中根据数据需求进行灵活调整。
如图6所示，在建筑信息管理bim技术中，尤其对于隧道施工工程的数据模型建构过程中为了实现对于施工工程数据与机器的相互识别认证，从而对于需要对于隧道数据模型进行模型的导入操作，从而导入到bim系统之中进行处理，将隧道模型根据不同的实际施工条件进行模型的拆分和规划，形成针对不同隧道模型的划分手段，从而提高隧道施工工程效率，将隧道模型划分完毕之后，通过pbs结构自动挂接方法，对隧道数据模型进行建构和数据名称统筹整理，将建构和数据名称统筹整理后的隧道数据模型进行导出操作，上述描述展现了pbs在整体bim中的工作环节，对于施工中的隧道模型划分具有重要的指导意义。
s2-3，根据输入的截面参数结合，计算出隧道截面开挖轮廓；通过扫略的方式沿着隧道中心线将隧道截面开挖轮廓扫略成开挖面，再将开挖面与开挖体按照先左右侧导坑后上下的顺序进行布尔运算，获得开挖结果模型。
如图3所示，对于传统隧道形状采用s2-3的步骤进行施工操作，绘制隧道数据模型；
在隧道数据模型进行划分之后，通过pbs进行自动挂接，具体步骤如下：
①，读入pbs数据；按照pbs编码、pbs描述、pbs分类的数据列，对pbs数据在excel中进行组织并读入系统；导入的pbs数据包含列标题，其中pbs编码是必须具备的数据列；若导入pbs数据除了pbs编码、描述、分类之外，还包含属性或工程量等信息，需将其添加到对应的pbs之后；对于数据列组织的顺序无任何要求，任意顺序数据列组织都支持导入并识别；支持单一pbs导入，也支持pbs工程量同时导入。
②，对数据列进行映射；对导入的pbs数据将pbs编码、pbs描述、pbs分类的数据列，进行一对一的映射；
③，录入pbs每一级的结构编码样例；pbs编码样例支持：工程项目、单位工程、分部工程、分项工程、单元工程、施工单元6个级别；
对于该pbs编码进行任意级别个数导入，
实际编码级别1＝工程项目+单位工程；
实际编码级别2＝工程项目+单位工程+分部工程；
实际编码级别3＝工程项目+单位工程+分部工程+分项工程；
实际编码级别4＝工程项目+单位工程+分部工程+分项工程+单元工程；
实际编码级别5＝工程项目+单位工程+分部工程+分项工程+单元工程+施工单元；
其中工程项目不填入编码样例，从而生成一个临时根节点，在预览功能中，对该根节点进行自动校正，将工程项目设置为根节点。
④，进行自动挂接；自动计算各级结构，自动组织子父节点关系；
从pbs数据池中提取指定的数据列，放入有效数据池中；
依次获取工程项目、单位工程、分部工程、分项工程、单元工程、施工单元六个级别样例编码；根据编码样例解析编码占位符，其算法为：
设pbs编码列当前编码为cur_code，pbs编码列下一个编码为next_code；
那么遍历pbs编码列，若cur_code_len！＝next_code_len则占位符replacechar为空字符串，这时pbs编码中每一级编码的位数都是不一样的，其中cur_code_len为前一个编码字符串长度，next_code_len为后一个编码字符串长度；
若cur_code_len＝＝next_code_len则为等位pbs编码，这时pbs编码每一级别的编码位数是一样的，对于每个级别不足的位数，采用某个特定占位符为补充，例如0；此时对pbs编码的每一个字符进行遍历：设：在同一个位置j，cur_code的字符为cur_code_c，next_code的字符为next_code_c；若cur_code_c！＝next_code_c，next_code_c即为占位符replacechar；解析pbs编码样例每个级别的编码位数；根据编码样例解析编码占位符，其算法为：若占位符replacechar为空串，则每级编码位数为实际编码样例位数；若占位符replacechar不是空串，则需要对每级编码样例剔除掉占位符的位数；根据占位符、各级位数对pbs编码从顶层到底层依次查找匹配，先查找父节点fnode，组织父节点fnode的属性及其他数据，并将父节点fnode的属性及其他数据与pbs编码绑定；查找父节点的下一级节点集合vchildren，并把vchildren挂接到fnode下；直到查找的子集合vchildren为空结束查找匹配；
⑤，pbs结构预览；对挂接好的pbs结构提供预览功能。
将pbs挂接完成的隧道数据模型进行导出操作，具体执行步骤如下：
㈠，选择工程施工模型导出格式，以及工程施工模型导出格式对应的版本；
㈡，根据选择的工程施工模型，遍历其工程施工模型数据组织结构，并对此工程施工模型数据组织结构进行分类，分类方式为product装配，part零件，body体；
㈢，根据获取到的工程施工模型数据组织结构，以及需要的新的模型数据组织结构，组织出新的需要的模型组织结构；
㈣，根据获取到的模型数据组织结构，以及对模型数据组织结构分类后获取到的对应模型属性数据和模型图形实体b-rep数据，分别填到对应的工程施工模型数据组织结构中；根据施工工程模型数据组织结构要求，以及对设计结果模型数据组织结构分类后分别获取的对应模型属性数据和模型图形b-rep数据，分别填充到对应工程施工模型数据组织结构各分类节点；根据施工工程模型数据要求，工程施工模型，设计结果模型的几何实体数据由模型几何model和构造几何construction两部分组成，其中模型几何是指拓扑类指向的类，不包括具体的形状信息，构造几何是指模型几何中的声明，包含实际的形状信息，模型几何和构造几何的关系如下：其中抽象几何对应模型几何，具体几何对应构造几何；根据施工工程模型数据要求，工程施工模型数据导出需要进行精度和模型空间控制，工程施工模型数据导出的精度和模型空间控制公式如下：
设定a：绝对最小值(10e-6)
b：规格化最小值(10e-10)
c：曲线曲面的逼近精度(10e-3)
d：被认为是零的最大值(10e-11)
模型空间计算算法：
modelspace＝a/b＝10e-6/10e-10＝10e4。
㈤，根据选择的工程施工模型导出格式以及导出格式对应版本，将获取到的工程施工模型数据组织结构以及各分类结构上的工程施工模型属性数据和图形b-rep数据导出。
如图4和5所示，其中对于方形隧道或者棱形隧道在具体执行操作中也进行了仿真实验，对于新式的隧道模型同样采用s3-4的步骤进行执行操作，其实在本发明的双侧壁导坑法隧道开挖单元划分的执行过程中，进行拆分的过程完全一致，最终也能够实现相应效果，针对将来根据不同的施工方案转换隧道形状，在采用该双侧壁导坑法隧道开挖单元划分时，一样实现快速的拆分模型，从而对实际的隧道施工工程起到显著的指导意义。
上述步骤采用的数据需要依次进行数据的获取操作，这样才能够准确的得到施工每一阶段需要的工程材料数据，能够保证工程项目预算精度。而且针对不同的数据其配置内容并不相同，这需要根据不断的试验才能够得到上述工作方法。
尽管已经示出和描述了本发明的实施例，本领域的普通技术人员可以理解：在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由权利要求及其等同物限定。