三维空间分析模组实现了部分三维地理信息系统的功能，并进一步使三维地理信息系统迈向了虚拟现实。在系统的实现过程中利用IDL坐标系统的转换功能，成功地将外部数据可视化；实现了三维数据的实时获取，并解决了三维实体空间编辑与操作的实际问题，实现多视角多解析度的观察。

数字地形模型(DigitalTerrainModel，DTM)是地形表面形态属性信息的数字表达，是带有空间位置特徵和地形属性特徵的数字描述。数字地形模型中地形属性为高程时称为数字高程模型(DigitalElevationModel，DEM)。从数学角度讲，高程模型是高程Z关于平面坐标X、Y两个自变数的连续函式。数字高程模型(DEM)只是它的一个有限的离散的表示。在地理信息系统中，DEM模型是建立DTM的基础数据，其它的地形要素可由DEM直接或间接导出，称为“派生数据”，如坡度、坡向等。GIS的基础是空间资料库，三维地理空间定位和数字表达是地理信息系统的本质特徵，地形数据(如DEM等)作为空间资料库的某个特定结构的数据集合，或所有这些数据集合的总体，被包含在地理信息系统中，成为它的核心部分的实体。显然，对地形空间数据的真三维显示和在三维空间的查询与分析也是GIS的核心内容之一。研究使用互动式数据语言IDL(InteractiveDataLanguage)作为开发平台，开发一套基于DEM数据的三维显示、浏览和空间分析系统Analysis3D。该系统的实现使三维地形的展示更加迅速，空间分析功能在三维地形上实现，分析结果更加形象。

本系统的设计内容包括DEM数据读取、数据可视化、数据分析与计算以及数据的编辑与操作等。根据这些设计内容绘製了系统流程图(见图1)。

图 1 系统流程图

（1）DEM数据读取与显示

DEM档案结构描述使用ArcView导出的ASCIIRaster的数据档案，导出档案的扩展名为*.asc。用记事本打开导出的档案，如图2所示。前6行为头档案，这个头档案中包含了地理编码信息，后面是高程矩阵。头档案的意义如表1所示。

图2ArcView导出档案的头档案

DEM数据读取IDL有广泛的数据接口，对于ASCII档案可以使用“readf”函式读取，这一过程中需要构造两个变数，head=makearray(1，6，/string)用于构造字元型矩阵，用于装载头档案；然后使用“strjoin，strsplit，strtrim”等函式在字元之间添加空格并去除多余的空格，这样有利于单个数值的提取。最后使用“where”函式对数值进行準确的定位，通过上述操作可以成功地获取ASCII档案所带的地理编码信息。同样道理Data=makearray(ncols，nrows，/float)用于浮点型的高程矩阵，并通过“reverse”函式转换矩阵方向，得到正确的显示效果。

表1 ArcView导出档案头档案的意义

（2）三维编辑与操作

编辑功能是任何一款三维建模与分析软体所必须的功能，对于三维物体的编辑，IDL是通过矩阵运算来实现的，将原始矩阵与变换矩阵做简单的数学运算就可以轻鬆地实现三维实体的旋转、缩放与平移。

图3旋转,平移,缩放功能的实现

旋转、平移和缩放功能的实现使用Trackball对象的Update方法实现旋转、平移和缩放，bHaveTransform返回的是Update是否成功的布尔值(即是否旋转)，如果成功返回为1，反之为0。旋转功能使用TRANSFORM关键字返回变换矩阵，GetProperty获取当前模型的状态矩阵，SetProperty设定当前模型的矩阵为变换后矩阵。同样平移使用TRANSLATE关键字，缩放使用SCALE关键字。如图3所示为旋转、平移和缩放。

图4 滑鼠获取坐标值

读取空间三维坐标为了更直观地显示DEM的地理属性，就需要通过滑鼠的移动实时地获取DEM上任意一点的地理坐标，滑鼠的移动涉及到滑鼠的Motion事件，通过滑鼠事件控制可得到任一点坐标。归一化坐标还要用到萤幕坐标与数据坐标的转换。

程式通过PickData函式结合event.x，event.y实时获取各点的萤幕坐标，然后再将萤幕坐标与数据坐标相互转换，最终显示的是数据坐标。如图4所示框中显示的就是DEM的地理坐标。

结论与讨论Analysis3D的开发，实现了部分地理信息系统的空间分析功能，为将来进一步开发三维地理信息系统奠定了技术基础。利用IDL坐标系统的转换和读入外部数据的功能，成功地将外部数据读入并完成可视化。通过IDL编程实现了三维数据的实时获取，并解决了三维实体空间编辑与操作的实际问题，使三维地形的观察更随意，实现多视角、多解析度的观察。

根据相应的算法，实现空间分析与立体量测，所有过程均在三维空间中完成，使分析结果更形象。同时完成了三维空间中提取等高线，以及自动标注等高线功能，使这一工作的生产率大大提高。实现了二、三维互动式浏览，不但使只有专业人士才能解读的等高线信息变为三维地形，而且使二维和三维的显示互动，得到良好的视觉效果。这种三维的浏览与显示做到了“所见即所得”，成为进一步研究三维地理信息系统与虚拟现实之间关係桥樑。

为了更真实地反映高层建筑的结构受力状况，一般设计单位都引进三维空间分析程式TBSA进行结构计算。我院于80年代后期引进该程式后，已用于许多工程的计算，在使用中遇到一些问题，积累一些经验，本文拟对有关问题加以探讨。

在高层建筑结构分析中，处理竖向荷载作用主要有三种计算模式：一是忽略竖向桿件的轴向变形；二是考虑竖向桿件的轴向变形，竖向荷载一次施加；三是考虑竖向桿件的轴向变形，模拟施工载入。

第一种计算模式忽略了竖向桿件的轴向变形，使墙与柱，柱与柱轴向变形差和梁刚度对结构的贡献不能体现，对计算结果产生如下影响：(1)轴压比大的构件轴力比实际偏大10%～20%，轴压比小的构件轴力比实际小10%～20%；(2)大部分梁的剪力比实际受力小；(3)樑柱的弯矩则是一些部位比实际受力大，一些部位比实际受力小。由此可见，忽略竖向桿件的轴向变形使部分樑柱墙设计偏于不安全。显然，这种载入模式是不合理的，已为大多高层分析程式所捨弃。

第二种计算模式对轴向变形估计偏大。第三种计算模式对轴向变形估计偏小，越接近顶层越明显。实际上构件的轴向变形应是介于两种计算模式之间，原因如下：(1)现实建筑施工载入与模拟施工载入仍有区别，模拟施工载入实际上是将全部荷载逐层全部加上，下层变形不受上层变形约束，且不影响上层；而在施工中，全部活载恆载中的内外装修荷载、大部分内隔墙荷载是在结构封顶后施加于结构的。(2)钢筋混凝土构件是弹塑性变形材料，房屋建成后的长期载入使受力状况又向一次载入靠拢。

因此，在高层建筑结构设计中，权宜之计是考虑竖向构件的轴向变形，并用第二、第三种模式进行结构计算，在两者之间取一插值进行构件配筋设计。

平面框架分析方法中，次梁与框架边梁是铰接的，而三维空间分析程式则考虑次梁与主梁弹性连线(如图5)，按其弹性刚度关係、位移协调来计算次梁端部弯矩和框架边梁的扭矩。两者的计算结果差异如下：(1)空间分析中由于边梁对次梁的约束，使次梁产生端部弯矩，框架边梁产生扭矩，平面框架分析方法中两者均为零；(2)由于端部弯矩，使P次梁空间分析中端部剪力较平面框架分析大。

图5 结构平面及框架边梁分析结果

实际上，钢筋混凝土构件在很小的扭矩下就会产生裂缝。根据试验，裂缝出现后，抗扭刚度会急剧下降，如图6。框架梁开裂后，其抗扭线刚度GI_T/L与次梁抗弯刚度EJ/m比值小，发生内力重分布(此时EJ亦可能减小，但辐度远较GI_T小)，框梁的扭矩减小，次梁的跨中弯矩相应增大。因此，其实际受力状况应介于平面框架分析与三维空间分析之间。

图6 扭矩- 扭转角试验曲线图

因此，运用三维空间分析软体时，应对框架边梁和次梁配筋进行人为干预，适当降低框架边梁扭矩，而次梁跨中弯矩应适当加大，否则会使梁偏于不安全。若偏于安全，可忽略次梁梁端弯矩，按简支计算次梁力。

在高层抗震结构中，结构延性越好，耗散的地震能量越多，抗震性能就越好。连梁是影响结构延性的主要构件之一，且连梁本身又是主要的抗震耗能构件。结构在承担地震前期地震脉冲冲击后，连梁首先屈服，它的弹塑性变形既可消耗地震能量，又可降低结构刚度，减小结构承担的地震能量，从而达到保护结构的目的。因此，要求连梁有良好的变形能力。为了使结构有良好的延性，连梁有良好的变形能力，就要求控制连梁高跨比，减小连梁刚度，採用严格的构造措施。因此，有人片面地认为连梁的高跨比越大，则结构的延性、连梁的耗能作用越小。实际上，不同部位的连梁对结构的贡献不同，且其耗能作用也不同。一般说来，连线不同薄壁柱的连梁贡献大，耗能作用大，如图7中的L-1。连线同一薄壁柱不同坪墙肢的连梁贡献小，耗能作用亦较小，如图7中的L-2。

图7 连梁平面

由于连梁是主要的弯曲耗能构件，为了使连梁如期屈服，就要求连梁实配弯矩筋不宜过大，且要控制其实测屈服强度与标準屈服强度的比值不大于1.25，以免使连梁受剪破坏，或者连梁未能如期屈服，结构薄弱部位发生转移，造成其它攀登的屈服破坏。

(1)高层建筑结构的内力介于模拟与不模拟施工载入间；

(2)框架边梁开裂后，抗扭线刚度下降，内力发生重分布；

(3)连线不同薄壁柱的连梁耗能作用大，对结构延性的贡献亦较大；

(4)在传力路径上三维空间分析较平面框架分析更合理；

(5)单片剪力墙平面外弯矩较大时，宜对其平面外配筋进行补算。