本发明涉及三维可视化领域和数值模拟领域，更具体地说，涉及一种三维可视化技术与数值模拟技术融合的实现方法与系统，主要应用于工程地质、灾害分析、物理模拟、工程建设、建筑设计、交通设计与施工等领域。
背景技术：
三维可视化技术：三维可视化技术是通过三维模型直观表达、理解目标物体的数字化技术，是物理实体的数字化映射。三维建模的一般过程包括：(1)获取研究对象或研究区域的各个层面的控制数据，如地表测绘数据、钻探获得的地下各个地层的上限和下限数据等；(2)利用各个层面的控制数据模拟研究对象各个分界曲面，如地层曲面、构造面等，得到研究对象的每个单元体的边界面，如各个地层单元的上下、左右和前后分界面，其中上、下两个面一般用曲面表达，左、右、前、后四个面一般用竖直的平面表达；(3)依据各个单元的边界曲面建立研究对象的每个实体单元，可采用多种三维数据结构进行实体单元的构建，如三棱柱结构、四面体结构、八叉树结构、b-rep结构、两种基本结构的混合结构等等；(4)各个单元实体的组合就是研究对象的实体三维模型。可以在三维模型上进行多种属性分析或空间分析，如工程开挖的工程量分析、滑坡等地质灾害的空间形态模拟等。
数值模拟技术：数值模拟技术是一种对研究对象的发展趋势进行动态模拟的技术，在工程建设、物理模拟、灾害分析等领域应用广泛。进行数据模拟分析，一般会把研究对象划分成一系列的网格，通过网格间物理力学参数的计算与传递计算，获得各个网格单元的运动趋势。如离散单元法，非连续变形分析(dda)方法、有限差分法等等。
三维可视化技术关注三维空间实体的构建、表达与属性-空间分析，但是三维模型的建立过程，是一个多次判断、多次处理的复杂而耗时的过程，在这个过程中，需要多次进行人机交互操作，而且，三维模型的各个实体单元之间联系密切，哪怕是一个单元的一个边界面上的一个控制点发生变化，对应的边界面、对应的实体单元、临近的一个或数个实体单元都要发生变化。这意味着甚至整个三维模型都要重建。因而，基于三维可视化技术对地质过程进行动态模拟面临严重的挑战。而数值模拟技术关注物理力学参数的计算、网格间参数的传递及过程的发展趋势模拟，但是在精细刻画三维实体以及属性-空间分析方面功能欠缺。两者具有很强的互补性。但是到目前为止，两者的有机结合还没有成功的先例。
技术实现要素：
本发明要解决的技术问题在于，针对现有技术中数值模拟技术还没有与三维可视化技术有机融合的技术缺陷，提供了一种三维可视化技术与数值模拟技术融合的实现方法与系统。
根据本发明的其中一方面，本发明解决其技术问题所采用的技术方案是：构造一种三维可视化技术与数值模拟技术融合的实现方法，包括如下步骤：
s1、在三维建模过程中获得研究对象每个实体单元的每个边界面的控制点；
s2、把每个边界面的控制点与数值模拟网格的节点进行对应分析，建立这些控制点与数值模拟网格节点的一一对应关系；
s3、通过数值模拟技术对所述数值模拟网格节点进行模拟；
s4、根据所述一一对应关系以及步骤s3得到的模拟结果，得到控制点在未来时刻的新的位置坐标；
s5、根据所述新的位置坐标建立未来时刻的三维模型。
优选地，在本发明的三维可视化技术与数值模拟技术融合的实现方法中，所述步骤s5具体为：
依据各个边界面上的控制点在未来不同时刻的坐标集，通过三维可视化中散点成面的方法，建立研究对象在未来不同时刻的三维模型。
优选地，在本发明的三维可视化技术与数值模拟技术融合的实现方法中，还包括步骤：通过获得的未来时刻各个控制点的三维坐标和未来时刻研究对象的三维模型，以未来某个时间点的三维模型为初始条件，计算未来另一时间节点的变化参数和进行发展趋势分析。
优选地，在本发明的三维可视化技术与数值模拟技术融合的实现方法中，所述步骤s5中建立三维模型具体为：得到研究对象在某个时间点的瞬态三维模型和/或在所述瞬态三维模型上基于三维可视化技术的属性-空间分析功能进行关键时间节点分析。
优选地，在本发明的三维可视化技术与数值模拟技术融合的实现方法中，所述步骤s5中建立三维模型具体为：基于步骤s4得到的未来时刻各个网格节点的三维坐标，得到未来任意时刻瞬间的真三维静态实体模型。
优选地，在本发明的三维可视化技术与数值模拟技术融合的实现方法中，所述步骤s5还包括步骤：通过未来不同瞬间的三维模型来实现研究对象的三维动态逼真模拟。
根据本发明的另一方面，本发明解决其技术问题，还提供一种三维可视化技术与数值模拟技术融合的实现系统，包括如下模块：
控制点获取模块，用于在三维建模过程中获得研究对象每个实体单元的每个边界面的控制点；
对应关系建立模块，用于把每个边界面的控制点与数值模拟网格的节点进行对应分析，建立这些控制点与数值模拟网格节点的一一对应关系；
数值模拟模块，用于通过数值模拟技术对所述数值模拟网格节点进行模拟；
数据处理模块，用于根据所述一一对应关系以及数值模拟模块得到的模拟结果，得到控制点在未来时刻的新的位置坐标；
模型建立模块，用于根据所述新的位置坐标建立未来时刻的三维模型。
优选地，在本发明的三维可视化技术与数值模拟技术融合的实现系统中，所述模型建立模块具体为用于：
依据各个边界面上的控制点在未来不同时刻的坐标集，通过三维可视化中散点成面的方法，建立研究对象在未来不同时刻的三维模型。
优选地，在本发明的三维可视化技术与数值模拟技术融合的实现系统中，所述模型建立模块中建立三维模型具体为：得到研究对象在某个时间点的瞬态三维模型和/或在所述瞬态三维模型上基于三维可视化技术的属性-空间分析功能进行关键时间节点分析。
优选地，在本发明的三维可视化技术与数值模拟技术融合的实现系统中，所述模型建立模块中建立三维模型具体为：基于数据处理模块得到的未来时刻各个网格节点的三维坐标，得到未来任意时刻瞬间的真三维静态实体模型。
本发明通过研究并确定对象各个实体单元的各个边界面上的控制点与数值模拟剖分的网格节点之间的一一对应关系，用网格节点代替边界控制点的方法实现了三维可视化与数值模拟的网格对接，从而使两者可以有效融合：在三维建模过程中获得研究对象每个实体单元的每个边界面的控制点，把每个边界面的控制点与数值模拟网格的节点进行对应分析，获得这些控制点与数值模拟网格节点的一一对应关系，通过数值模拟技术得到这些控制点在未来不同时刻的新的位置坐标。依据各个边界面上的控制点在未来不同时刻的坐标集，通过三维可视化中散点成面的方法，建立研究对象在未来不同时刻的三维模型。由于可以建立未来不同时刻的研究对象的三维模型，因而既可以得到研究对象在某个关键时间点的瞬态三维模型，也可以在此瞬态模型上基于三维可视化技术丰富的属性-空间分析功能进行关键时间节点分析。进而可以通过未来不同瞬间的三维模型来实现研究对象的三维动态逼真模拟。
本发明通过研究对象各个实体单元的各个边界面上的控制点与数值模拟剖分的网格节点之间一一对应关系的确定，用网格节点代替边界控制点的方法得到这些控制点在未来不同时刻的新的位置坐标，实现了三维可视化与数值模拟的网格对接，从而使两者可以有效地融合起来，并得到研究对象在未来不同时刻的三维模型，为研究对象的过程模拟提供更好的技术支持。
附图说明
下面将结合附图及实施例对本发明作进一步说明，附图中：
图1(a)是三维可视化的地层结构及边界面上的控制点示意图；图1(b)是图1(a)数值模拟划分的块体及网格示意图；
图2是本发明的三维可视化技术与数值模拟技术融合的实现方法的一优选实施例的流程图；
图3是本发明的三维可视化技术与数值模拟技术融合的实现方法的另一实施例的流程图；
图4是本发明的三维可视化技术与数值模拟技术融合的实现系统的模块组成图。
具体实施方式
为了对本发明的技术特征、目的和效果有更加清楚的理解，现对照附图详细说明本发明的具体实施方式。
下面以某个滑坡为例进行说明。为便于说明，以图1中的二维示意图代替三维进行说明，数值模拟软件以flac3d软件为例。其他软件的处理方法类似。
同时参考图2，本实施例的三维可视化技术与数值模拟技术融合的实现方法首先基于三维可视化软件建立滑坡各个地质单元的边界曲面，获得各个曲面上的控制点：图1(a)中有三个地层，地层i为粘土，ii为碎石土，iii是基岩。s为滑坡的滑面，位于基岩顶面。zk1、zk2等是钻孔。在三维可视化软件中建立地表、地层i、ii的底面曲面，滑面与地层ii的底面是同一个曲面。建立的数据依据包括地表测绘数据、钻探数据、地质剖面数据等，建面的方法是散点成面法。如zk1s、zk2s等钻孔与滑面s的交点和pms1、pms2等剖面pm与滑面s交线上的控制点共同构成了滑面s上的控制点，依据这些控制点(空间散点)，通过不规则三角网方法或grid格网方法或gocad中的dsi方法等等，都可完成滑面s的构建。在滑面s上标注zk1s等钻孔控制点和pms1等剖面上的控制点。其他面的构建方法类似。如图1(a)中的圆点和三角点。
其次，调用flac3d软件进行数值模拟：由于flac3d软件不能直接接收三维可视化格式的数据，我们通过ansys或其他软件中转一下，具体过程包括①三维可视化软件的曲面文件以通用的acis格式.sat存储，这种格式很多软件都可以打开，如autocad、rhino、catia、ansys等等；②把.sat文件导入到ansys中；③利用ansys-flac接口程序导入到flac3d中，进行块体及网格的划分，并准备计算。如图1(b)所示：地层i、ii、iii把研究对象划分为三个块体。块体划分好后，需要对每个块体进行网格剖分，一般会依据研究的关注点对不同块体进行疏密不同的网格剖分，如滑面附近的网格剖分较密。
然后确定计算前各个控制点的编号和坐标：计算之前要在地表、地层i底面和滑面上确定图1(b)中与钻孔控制点(zk1s点等)和剖面控制点(pms1点等)最靠近的块体网格节点的编号和对应坐标。以图1(b)右上角放大图中的zk1s和pms2点为例，由于最靠近zk1s点的网格节点是gi网格的右下角节点，可以把gi网格右下角的节点近似看做zk1s点，同理，把gj网格右上角的节点近似看做pms2点。其他各点的处理方法类似。这样可以获得每个曲面上每个散点(即每个曲面的控制点，包括钻孔控制点和剖面控制点等)的坐标及对应节点编号。如滑面上的zk1s点、zk2s点......等和pms1点、pms2点......等的节点编号及坐标。如果控制点不在某个节点附近，可以通过两个节点间的中点或等分点计算获得。这种用数值模拟网格节点来代替三维可视化曲面上各个控制点的方法，本发明称其为节点代替法。
再确定不同时步各个节点的新的坐标：在flac3d计算过程中，设置记录节点新坐标的时步间隔，进行计算。可以得到不同时步时散点对应节点的新坐标，即可以获得不同时步时钻孔控制点和剖面控制点的新的位置坐标；
再建立不同时步时的新曲面：这些控制点在不同时步时的新坐标就是三维可视化系统建立不同时步新模型的建面散点。把不同时步的散点坐标导入到三维可视化软件中，就可以建立不同时步时的三维模型，即得到未来不同时间瞬间的三维模型。
从而完成三维可视化技术与数值模拟技术的融合。
参考图3，其为本发明的三维可视化技术与数值模拟技术融合的实现方法的另一实施例的流程图。本实施与上述实施例的不同之处在于，在建立了不同时刻的三维模型后，可以以已经建立的某个时刻的三维模型为初始条件，采用相同的方法，计算未来另一时间节点的变化参数，从而对某个关键节点的参数与模型进行多次交叉验证。即步骤s1-步骤s5可以循环进行，直至得到理想的结果。
参考图4，其为本发明的三维可视化技术与数值模拟技术融合的实现系统的一优选实施例的原理框图，该三维可视化技术与数值模拟技术融合的实现系统，该系统包括顺次连接的控制点获取模块11、对应关系建立模块12、数值模拟模块13、数据处理模块14、模型建立模块15。控制点获取模块11用于在三维建模过程中获得研究对象每个实体单元的每个边界面的控制点；对应关系建立模块12用于把每个边界面的控制点与数值模拟网格的节点进行对应分析，建立这些控制点与数值模拟网格节点的一一对应关系；数值模拟模块13用于通过数值模拟技术对所述数值模拟网格节点进行模拟；数据处理模块14用于根据所述一一对应关系以及数值模拟模块得到的模拟结果，得到控制点在未来时刻的新的位置坐标；模型建立模块15用于根据所述新的位置坐标建立未来时刻的三维模型。
模型建立模块15具体为用于：
依据各个边界面上的控制点在未来不同时刻的坐标集，通过三维可视化中散点成面的方法，建立研究对象在未来不同时刻的三维模型。
模型建立模块15中建立三维模型具体为：得到研究对象在某个时间点的瞬态三维模型和/或在所述瞬态三维模型上基于三维可视化技术的属性-空间分析功能进行关键时间节点分析。
模型建立模块15中建立三维模型具体为：基于数据处理模块得到的未来时刻各个网格节点的三维坐标，得到未来任意时刻瞬间的真三维静态实体模型。
上面结合附图对本发明的实施例进行了描述，但是本发明并不局限于上述的具体实施方式，上述的具体实施方式仅仅是示意性的，而不是限制性的，本领域的普通技术人员在本发明的启示下，在不脱离本发明宗旨和权利要求所保护的范围情况下，还可做出很多形式，这些均属于本发明的保护之内。