不连续面三维网络模型图形显示及处理.doc

1 不连续面三维网络模型图形显示及处理摘要: 通过不连续面三维网络模拟获得了在模拟空间范围内具一定概率分布特征的不连续面, 这些不连续面的形态以圆盘形式表示,空间位置由圆盘中心点坐标及产状表示, 大小由圆盘直径体现。关键词:不连续面三维网络模型图形显示处理通过不连续面三维网络模拟获得了在模拟空间范围内具一定概率分布特征的不连续面, 这些不连续面的形态以圆盘形式表示,空间位置由圆盘中心点坐标及产状表示,大小由圆盘直径体现。根据模拟的结果可以进一步进行岩体稳定分析、岩体水力学性质等研究,但首先要解决的问题就是三维网络模型的图形化,即可视性,没有图形,研究将难以深入。为了使得研究者能直观地看到不连续面的三维网络模型, 并能对网络模型进行切割产生平面网络, 进而对这些图形直接进行处理, 比如输出等。笔者用 C++ 语言在 2 Borland C++ 集成环境下编写了一系列计算机程序,使得三维网络模型的应用更加方便可行。 1 程序组成本模块由以下文件组成: .cpp , .exe , , .prj , , .dat , , , , , , , hzktext , , hzf24 , , 。在 Borlandc++ 的集成环境中经工程文件编译, 可执行文件,在 DOS 或 Windows 环境下运行。 2 程序编制依据 结构面三维图形显示 求圆盘上各点三维坐标 3 对于一个节理圆盘,沿其法向观测时其为一圆形,沿其走向观测时其为一线段, 除此之外, 从其他任何角度去观测, 节理圆盘都为椭圆。 C++ 语言所提供的绘椭圆函数 ellipse() 只能绘制长轴为X轴、短轴为Y轴的椭圆, 而无法实现任意形状椭圆的显示,因此还得从圆盘轨迹着手,用 moveto() 、 o() 等函数来实现。由于三维网络模拟产生的节理是以具有中心点坐标、产状和直径的圆盘表示, 所以圆盘上的每一点相对于模型来说都是确定的,都可以用一组三维坐标表达。求圆盘各点的坐标是出于这样的出发点:对照图 1, 对于一圆盘,很容易得到其倾向线,以倾向方向的半径 OA 顺时针( 或逆时针) 旋转一角度,得到另一半径 OB ,求出 B 点的坐标, 旋转角的取值由 0°到 360 °, 即得到该圆盘的轨迹。对于倾角为α、倾向为β、半径为 R 中心点坐标为( Xo,Yo,Zo )的节理圆盘,其投影到水平面上则为一椭圆, 该椭圆长半轴沿走向方向( 以下所说的走向均指服从右手法则, 长度为a= R, 短半轴沿倾向方向, 长度为b= Rcosa 。 4 如果平移并沿 Z 轴旋转坐标系,使 X 轴与圆盘在水平面投影的椭圆长轴一致且原心在椭圆中心点, 则椭圆上一点的坐标可表示为: X’=a cosq = Rc (2-1-1) Y’=b sinq = Rcosa sinq (2-1-2) q 为所求坐标点与原心连线与椭圆长轴的夹角(逆时针为正) ,也就是圆盘上所求坐标点与圆盘中心点连线与走向线的夹角。再把坐标变换回来,即得到圆盘上任意一点的三维坐标,分别为: X= Xo +X’ cosb +Y’ sinb (2-1-3) Y= Yo +X’ sinb +Y’ cosb (2-1-4) Z= Zo - Rsin α’(2-1-5) 5 α’为所求坐标点与圆盘中心点连线的倾角α’是这样确定的,圆盘与水平面的交线必为走向线,在图1 中, OA 为圆盘倾向方向半径, OB 为与 OA 夹角为θ的一半径, 两半径与水平面分别交于 A’点和 B’点, 圆盘中心点 O 在水平面上的垂直投影点为 O’,A’B’即为走向线,O’A’B’构成水平面。因为 OO ’与A’O’垂直、 OB ’与B’O’垂直,并且 OA ’垂直于 A’B’( 走向线与倾向线) ,所以有: OA ’= OO ’/sina (α为圆盘倾角) (2-1-6) OB ’= OO ’/sina ‘(a‘为 OB 与水平面夹角) (2-1-7) OA ’= OB ’ sinq (2-1-8) 即: sina ‘= sina sinq (2-1-9) 6 所以 Z= Zo - Rsina sinq (2-1-10) 求圆盘在视面上的平面坐标计算机屏幕是个平面,要把一个三维物体在屏幕上显示出来,首先要把该物体各点的三维坐标变换为平面坐标。相对于物体, 屏幕是个投影面, 该面位于观测者与被观测物体之间, 如图 2。对于表征物体的三维坐标系 O_XYZ , 设观测点的位置为 V, 其直角坐标为(X