IFC标准的装配式建筑空间自组织建模方法论文
摘要：基于国际IFC标准，研究装配式建筑构件的空间实时定位算法，将微机电系统（Micro-Electro-MechanicalSystem，MEMS）传感器应用于建筑业，对微机电系统传感器进行研究和开发，借助MEMS-IMU记录装配式构件在安装过程中的空间位置和姿态变化数据，结合已有的BIM构件库，研究基于IFC标准的建筑信息模型自动生成算法，以参数化获得最终IFC标准数据模型，为装配式建筑竣工模型的自动生成奠定基础。
关键词：建筑信息模型；IFC标准；MEMS传感器；空间定位
随着土木建筑工程项目的规模越来越大，建筑造型也越来越复杂，不同专业、各方人员协同困难，越来越多的工程项目期望使用或者正在使用BIM技术去集成、整合并分析建筑全生命周期的各种信息，加强信息共享，以便于对整个工程进行管理。同时，随着装配式建筑在国内的推广，越来越多的工程项目选择使用装配式技术来提高整体生产效率，提升工程质量。目前，国内外许多学者在装配式建筑和BIM技术结合方面做了不少研究[1-3]。目前大部分模型自动生成研究都是基于三维激光扫描系统及图像重建技术的逆向工程[4-10]。综合来说，激光扫描及图像重建技术在建筑业中的研究及应用还处于初级阶段，虽然许多研究已取得阶段性进展，但生成满足竣工交付要求的三维模型仍面临着很多技术挑战和实施困难。本文提出基于IFC标准的装配式构件空间定位方法，通过研究IFC标准中构件信息表达方法，结合构件位置及姿态捕捉算法，可生成符合IFC标准的建筑信息模型。
1基于IFC标准的构件空间定位算法
1.1符合国际标准的构件表达建筑项目构件种类繁多，在模型建立初期就应该对构件命名规则进行细化约定，规范项目参与人员对项目的设计、修改等行为，提高数据交互效率，保证建筑信息模型数据质量。建筑构件分类可借鉴北美地区广泛使用的OmniClass[11]标准Table23-Products中的分类方法，并在此基础上进行扩展，以对自组织生成模型中的构件进行命名。构件命名可根据装配式构件库中的构件基本信息按照预定规则自动生成并转为Unicode码，添加IFC文件实体属性Name字段。IFC标准采用EXPRESS语言描述[12]，并定义其包含的所有数据信息，唐春凤等[13]阐述了IFC文件的一般结构和EXPRESS语言。IFC标准技术架构分为4层，由上至下分别为：领域层、共享层、核心层、资源层。定义于核心层中的IfcRoot实体直接或者间接派生出资源层之外定义的每个实体。IFC模型中的三种基本实体类型(IfcObjectDenition，IfcPropertyDenition，IfcRelationship)都是由IfcRoot派生而来，实体间的关系见文献[14]Express-g图。三维几何建模常用到的构件在IFC标准中对应的实体都由实体IfcProduct派生或其子类派生，IfcProduct是对与几何或空间环境相关的任何对象的抽象表达，其子类通常设有形状表达和项目结构所涉及的对象坐标。1.2构件的空间定位目前国内外对装配式建筑安装过程中构件定位的研究基本处于起步阶段，本研究需要测量构件的实时位置及姿态，再根据其几何外观等参数，在远端计算机上实时显示及监测安装进度，进一步可以实现信息化的进度、材料、设备管理等内容，提升管理水平，提高工程质量，最后生成能用于后期运维的建筑信息模型，所以需要实时跟踪并记录构件在安装过程中的各项数据，考虑到施工现场条件限制以及信号遮挡等问题，本研究采用航位推算法，利用基于微机电系统的惯性测量单元(InertialMeasurementUnit，IMU)对装配式构件进行位置及姿态的追踪。IMU大多用在需要进行运动控制的设备，如汽车和机器人上，也被用在对姿态进行精密位移推算的场合，如潜艇、飞机、导弹和航天器的惯性导航设备等。相较于在现场使用GPS及微波定位等方法，其优点是体积小，抗冲击，可靠性高，寿命长，成本低，重量轻[15]，且不受施工场地使用环境限制，抗干扰能力强。但其也存在测量误差会随着时间的推移而增大的不足，所以需要通过一定的算法(如KalmanFiltering算法)来消除对应的误差，以得到满足系统需求的准确数据。1.2.1构件位置获取一个IMU一般包含有三轴加速度计和三轴陀螺仪，加速度计用来检测物体三个独立轴向的加速度数据，陀螺仪用来测量物体角速度数据。由于IMU只能采集到原始的加速度和角加速度信号，而不能直接得到构件的位置和姿态数据，所以需要利用相关算法处理器原始信号，以得到所需的位置和姿态信息。对于物体的加速度信号，可以通过时域积分，将加速度值a对时间t积分，同时给定初始速度，可以得到加速度计的速度函数：v(t)=∫0ta(t)dt=v*(t)+v0(1)将式(1)再次对时间t积分，同时给定初始位移，得到加速度计在局部坐标系下的三轴位移：r(t)=∫0tv(t)dt=r*(t)+r0(2)式(1)、(2)中：a(t)为加速度计原始信号，v*(t)为a(t)的原函数，v0为初始速度，r*(t)为v(t)的原函数，r0为初始位移。由于所测量的构件不能被看作空间中的一个质点，而加速度计测量的加速度数据仅仅是固连在自身上的坐标系中的数据，所以测得的数据并不是世界坐标系中的数据，这就需要进行进一步的坐标变换处理。此时就需要借助陀螺仪记录的方向参数。1.2.2构件姿态获取陀螺仪的使用和加速度计类似，它通过测量力矩计算角速率，通过角速率积分得到角度变化。一般建模过程中都会设置一个世界坐标系(WorldCoordinateSystem，WCS)，需要求得的数据为构件在世界坐标系下的绝对姿态，而固连在IMU上的坐标系可视为局部坐标系，假定局部坐标系和世界坐标系的初始位置重合，从世界坐标系到局部坐标系的变换可以用欧拉旋转或者四元数旋转等方式表达，为表达更加方便直观，这里选择用欧拉角表达变换过程，不妨设旋转次序为x-y-z，三个欧拉角为α、β、γ，则从世界坐标系Pw到局部坐标系Pι的变换为则从局部坐标系数据求世界坐标系数据只需求式(3)的逆变换，即：Pw=C-1(α)C-1(β)C-1(γ)Pι(5)不难看出沿轴旋转变换矩阵的逆矩阵即为绕坐标轴旋转一个相反的角度，即可得式(5)的等价公式：Pw=C(－α)C(－β)C(－γ)Pι(6)由于篇幅有限且误差处理及补偿算法较为复杂，在此不做论述。至此已经求得某个构件在世界坐标系的空间位置及姿态。1.2.3符合IFC标准的构件空间定位表达获取了构件空间位置及姿态数据，下一步就是将这些数据用符合IFC标准的语句表达出来。在IFC标准中，构件位置通过IfcObjectPlacement实体表达，它是定义对象坐标系的一种抽象父类，对于每个有形状表达的产品都需要提供IfcObjectPlacement。构件坐标表达形式有相对坐标、绝对坐标和网格三种，本研究采用相对坐标表达，由IfcProject表达项目的世界坐标系统，其下层的IfcSite包含了该场地的单一地理参考点的定义(使用WGS84坐标系统定义的经度、纬度及海拔)，IfcBuilding、IfcBuildingStorey等位置均以其上层坐标系作为参考坐标系，以此表达自身的局部坐标系位置。采用此种表达方法也更容易和前文所述的IMU采集到的数据结合。由于IFC标准规定每个IFC文件有且仅有一个IfcProject实体，而可以包含多个IfcSite等实体，所以每个构件实体需要通过IFC位置表达语句层层嵌套，最终关联到IfcSite实体的坐标系。此过程中最重要即获取当前局部坐标系在上层坐标系中的表达。由于IFC表达中一个坐标系需要原点坐标、Z轴和X轴向量确定，其中原点可以通过IMU数据解算，以得到现有构件坐标系在上层坐标系中的坐标原点(r1r2r3)，Z轴和X轴参考方向需要根据IMU数据解算得到的姿态转角进行变换得到。默认局部坐标系的Z轴坐标为(001)T，X轴坐标为(001)T，则局部坐标系Z轴、X轴在上层坐标系中的向量表达即为()zzzC123=-aRTSSSSSSSVXWWWWWWW()()001C-bC-cRTSSSSSSSVXWWWWWWW()xxxC123=-}RTSSSSSSSVXWWWWWWW()()100C-iC-{RTSSSSSSSVXWWWWWWW(7)则该构件的局部坐标系表达即为IFCCARTESIANPOINT((r1，r2，r3))；IFCDIRECTION((z1，z2，z3))；IFCDIRECTION((x1，x2，x3))。
2基于IFC标准的空间自组织建模
2.1自组织建模整体流程基于国际IFC标准，以NMBIM软件(上海交通大学BIM研究中心自主研发协同平台)为基础平台，使用VisualC++开发装配式建筑构件自组织建模软件。整个自组织建模软件分为测量单元，传输单元，模型生成单元三大部分，各个单元之间协同工作流程见图1，主要步骤如下：(1)装配式构件吊装前，通过BIM构件库[16]获取构件几何信息、材料信息、属性信息等基本信息，同时将这些基本信息读入缓存；(2)将测量单元贴附在构件表面特定位置并将测量单元初始化，开始构件吊装；(3)模型生成单元根据接收到的数据在建模平台实时显示安装进程；(4)某一构件安装完成，根据最终位置及构件信息，在数据库中生成安装完成的构件信息相对应的信息记录，同时拆下构件上的测量装置，进行下一个构件的吊装；(5)所有构件安装完成，保存安装过程中的数据文件至数据库，根据选项生成对应的IFC物理文件。2.2IFC模型文件的自动生成IFC文件分为文件头Header和数据Data部分，输出IFC文件前先根据预定义信息初始化文件头部分，记录与整个交换文件结构有关的信息。接下来向数据缓冲区输出个人和组织信息定义、归属历史等所有IFC文件都需要含有的通用实体信息等，紧接着根据外部数据资料实例化IfcProject实体，确定其相关属性及其几何环境表达、项目单位定义等信息。实例化完成一些必要的'实体后，创建IfcSite、IfcBuilding、IfcBuildingStorey及其它Element实体的实例，同时在各个实例构造函数中做好相关实体信息关联，同时根据逻辑判断，在生成相关实体实例的时候要创建相关实体间的关系实体实例。在完成所有实体数据信息实例的创建后，对容纳有所有实例的容器进行遍历，最后一步即关联一个文件流对象，将上述容器内的所有实例通过算法转化为文本，输出到磁盘文件中，最后调用IFC文件校验工具，生成校验日志。2.3案例验证根据前文所述的模型生成流程，在实验室验证整个系统的可行性。本系统可以自定义项目信息或者根据已有IFC模型数据解析项目信息，在原有IFC模型基础上继续添加新的装配式构件。本例采用已有的IFC模型数据，对该模型右上角的构件进行模拟吊装。在实验室中，将IMU测量装置粘贴在构件下表面，并对测量装置进行初始化，开始装配式构件的吊装，在实验室对一根模拟柱子进行移动，在移动的过程中，构件定位装置会一直记录加速度信号和角度信号数据，完成吊装后对IMU元件采集到的数据进行解析，得到解析后的位置和姿态数据。(限于篇幅，图中只展示原始加速度信号图像，未展示时域积分后的速度和位移信号图像)和最终数据解析结果。之后依据构件库中相应构件的数据信息。最后进行下一个构件的安装，待所有构件安装完毕，导出IFC模型数据(本例只添加一根柱子作为示意)。在构件树中选中新添加的名为test_column的柱子，可以看到其几何信息、材料信息、位置及其它属性信息。误差基本控制在7%以内。IMU惯性单元存在漂移问题(本方案采用的IMU惯性单元的零偏为1.0mg，零偏稳定性为±1.5mg/℃)，若采用精度更高的IMU惯性单元，误差将进一步减小。
3结论
本文通过对微机电传感器及IFC标准的研究，提出基于IFC标准的装配式构件空间定位算法，对MEMS传感器追踪到的构件移动过程数据进行自动化处理，借助数值积分算法，将加速度和角加速度信号进行时域积分，得到构件移动的三轴位移和三轴转角，实现构件的实时定位。然后结合已有BIM构件库，研究和开发了基于IFC标准的建筑信息模型生成软件，为装配式建筑竣工模型的生成提供一种新的思路和方法。