基于几何和拓扑关系的桥台三维参数化设计

基于几何和拓扑关系的桥台三维参数化设计

目前，三维书法可视化技术广泛应用于桥梁信息的各个方面，如桥梁设计、桥梁结构的分析和计算、方案比较的可视化分析和评估、方案的竞争效果的展示、图纸设计结果的可视化检验等。相比之下，对桥梁建设后的管理、检验和维护等三维图形可视化应用的研究和成果很少。桥梁种类虽然繁多,但其结构组成和空间分布规律性很强,因此,采用三维参数化建模成为构建三维可视化系统的必然选择。在桥梁三维建模的几何拓扑关系组织上,也有参考机械设计领域的成果。如在机械领域,产品设计围绕工业化生产,主要以三维几何模型为主,为研究零件间复杂的关系,研究者建立了基于拓扑树和元胞模型的零件族模型。在数据结构研究上,文献以扫描法和参数化特征造型为桥梁模型的主要造型方法,半边数据结构为模型的数据结构。文献采用基于知识的参数化建模方法对桥梁的参数化建模进行研究,提出了桥梁设计知识表达模型,知识获取方法,桥梁结构拼装算法,桥梁智能变动算法。文献结合斜拉桥结构的特点,采用参数化方法对斜拉桥进行三维建模,实现了斜拉桥三维可视化。本文对三维几何体参数化建模的过程进行了深入探讨,并结合桥梁结构特点,对桥梁组件的三维参数化建模及桥梁整体拼接进行了研究,实现了9种桥梁结构组合的三维可视化。1参数方程的描述参数化设计(ParametricDesign)是指将影响设计结果的各种要素进行归类,将这些要素变成某个函数的变量,并设计函数体中变量间的关系,则可以通过改变函数变量的不同取值,来得到不同的函数结果,或者说通过改变算法,能够获得相同函数所涵盖的不同设计方案。参数化设计可用如下公式来描述:其中,F:(f1,f2,…,fn)是一系列参数方程,D:(d1,d2,…,dn)是F函数的变量,表示几何尺寸间的约束关系,X:(x1,x2,…,xn)是F函数的变量,表示获取的几何特征点的坐标,包括结果值。对于三维几何体的参数化建模,参数化设计要素主要涉及几何体尺寸描述参数,尺寸参数映射为空间顶点坐标,几何体组成关系描述,几何体间空间拓扑关系描述等。这些参数描述和拓扑关系可以模块的形式嵌入到程序中。由于本文所关注的主要是用于桥梁检测维护和桥梁外观表现的建模,因此,采用了基于边界表示法(boundaryrepresentation,B-reps)的实体建模方法。尽管桥梁构件种类繁多,造型各异,但可以搜集到全面的桥梁各组件横截面的几何信息和纵向的跨径拓扑信息,因此,采用的是基于特征参数法的截面设计来进行桥梁各组件的三维参数化设计,同时基于桥梁纵向的跨径拓扑约束建立桥梁的整体拼接关系。2空间分布的数化建模本文主要针对空心板、T梁、小箱梁、连续板、连续等高箱梁、柱式墩、盖梁柱式墩、轻型桥台和重力式桥台这9种主要组件结构进行三维参数化建模,其模型主要是通过描述横截面参数进行纵向拉伸来建模。桥梁组件参数化建模过程如下:1)从桥梁设计图中获得组件横截面和纵向拉伸几何尺寸和空间分布拓扑关系描述变量;2)建立模型的局部坐标系,根据组件的外围几何分布,确定组成组件外表面多边形及多边形共线、线段共点的关系,从而完成将组件外表面尺寸变量映射为描述组件的顶点空间坐标的过程。此时每个空间顶点的3个坐标值完全由几何尺寸变量的表达式表示。在此过程中合理地设置每个组件的局部坐标系原点位置对于由组件拼接成为部件,由部件拼接为整个桥梁的便捷性影响很大;3)确定每个外表面多边形的绘制形状和组成该多边形的空间顶点;4)利用几何体的多边形绘制函数进行各外表面的绘制,并通过颜色设置或纹理映射等增加物体的真实感。2.1桥梁次级参数说明由于桥梁组件截面类型较多,每一类的截面都有其特有的从设计参数到空间三维坐标的映射关系,依据参数化造型理论,需要确立的特征参数主要有以下3类:形状特征参数、约束参数、层次结构参数。1)首先通过形状特征参数描述桥梁各组件的外观几何尺寸。其中又分主要参数和次级参数,主要参数是在桥梁各组件的结构计算中,依据规范及各组件的主尺寸等参数列出绘制各组件所需要的尺寸和参数。次级参数一般是指通过一定的约束条件由主要参数确定的参数。如:重力桥台中帽梁伸出前墙及伸出左右侧墙的尺寸等。2)然后再通过定位约束参数以及层次结构参数进行桥梁各组件的拼接以及整桥拼接关系的描述。定位约束参数体现在2个层次关系上,一个是组件中各个几何体间的位置关系,如在盖梁柱式墩中墩柱位于盖梁下部的位置描述;一个是在使用组件拼接桥梁时各个组件的位置关系,如第一块主梁位于左侧轻型桥台的位置。层次结构参数主要是描述一座桥梁的组成层次关系。由于桥梁的组件种类众多,形状参数各不相同,需要对参数进行合理归类。根据一般桥梁的构件图,桥梁建模参数可以归纳为以下5个部分:(1)桥梁整体布局建模基本参数:包括孔跨信息、各组件空间位置拓扑信息等。(2)主梁建模参数:包括主梁横截面类型、截面几何尺寸以及主梁的片数等。(3)桥墩建模参数:根据桥墩类型,获取相应的几何参数。按照便于建模的原则,将桥墩模型划分成若干部分,根据各部分空间位置拓扑关系拼接成桥墩整体。(4)桥台建模参数:包括不同桥台类型各组件的几何参数,各组件的空间位置拓扑关系参数。(5)桥面、护栏和其他部分的建模参数。2.2t梁建模过程对空心板、简支变连续小箱梁、连续板、T梁以及等高连续箱梁5种主梁的建模采用基于参数化草图产生三维实体模型的方法,主要以主梁横断面图和相关特征参数为基础。下部结构包括轻型桥台、重力桥台盖梁柱式墩和柱式墩。这些结构组成复杂,均是在对其各组成部分建模的基础上,再根据其空间结构位置分布关系进行拼接。如盖梁柱式墩由盖梁、墩柱、支座、垫石等组件构成,柱式墩主要由方形墩或圆形墩柱、支座、垫石等组件构成。这些组件尽管形状和组成各异,截面参数描述均不相同,但其建模过程相同,下面以T梁为例来介绍其建模过程。1)各部分的参数描述如下:T梁高度H,跨中肋板宽度bF,中梁顶宽bTB,马蹄宽度bD,边梁外侧半宽bT,马蹄过渡高度bD,悬臂加腋宽度bF,跨中马蹄高度hD,悬臂端部厚度hY,支点马蹄高度hB,悬臂根部厚度hR,如图1左侧所示。2)定义其局部坐标中心在O点,通过其尺寸约束关系求解梁体各顶点在三维坐标系中的坐标,各顶点坐标求解时是按照各面之间的连接关系和空间顶点位置关系进行映射。3)根据横断面两边顶点的对应关系,求解后面各顶点的三维坐标。由于主梁横断面的前后面都是对称的,所以后面的各顶点坐标都在前面的各对应点的坐标的x轴上加上length,其中length为主梁长度。4)建立网格顶点数组GridVertices。5)定义多边形顶点组成数组frontBase来决定GridVertices中的顶点在不同面中的组织方式。6)将frontBase加入到几何体节点frontGeometry下,再将frontGeometry节点加入到绘制T梁的叶节点Geode下,T梁的一面绘制完,其它各面的绘制过程与此类似,这样绘制出来的各个面就形成T梁的三维模型,T梁的三维模型如图2(c)所示。根据上述理论,绘制出的5种主梁三维模型如图2所示。2.3算法设计与分析在桥梁各组件建模过程中,有一些特殊的图元不能用上述方法进行建模,需要设计特殊的算法来进行建模。如:墩柱中具有圆弧倒角的方墩、圆柱墩中过渡加宽部分、桥台的近似锥面护坡以及主梁绘制中的凹多边形面片等。2.3.1建立假设点坐标墩柱主要分圆形墩柱和方形墩柱两种,而方形墩柱又分为无倒角、直角倒角和圆弧倒角3种。在绘制圆形墩柱时,可使用OpenSceneGraph中预定义几何体,调用绘制圆柱体函数即可。某些圆形墩柱中间有过渡加宽部分,可采用以下方法进行建模,如图3所示,加宽段上圆柱的半径为r1,下圆柱的半径为r2:1)在上下圆面上将圆弧等分为若干段,比如将每1/4圆弧平分成4段,每段圆弧所对应的圆心角为θ/4,其中θ=π/2。2)定义以图中O点为此部分局部坐标原点,计算得到插入点A的坐标(r1*cos(θ/4),-(r1*sin(θ/4)),0),同理可得到B点的坐标为(r1*cos(θ/2),-(r1*sin(θ/2),)0,C点的坐标为(r1*cos(3θ/4),-(r1*sin(3θ/4),)0,其它插入点的坐标也由此计算得到。3)再按上下对应的每4个点绘制小四边形,由此逼近圆锥效果,如图中以ABFE4点绘制四边形,依次类推,得到圆形墩柱三维模型如图4所示。2.3.2基于弧ab的近似诉讼桥台的组件非常多,这些组件中大部分的三维参数化建模与前面介绍的组件建模方法一致。但左右侧近似锥面护坡较为特殊,由于护坡在总体形状上接近于1/4圆锥体结构,但其锥底面的2个弧半径往往是不等长的,直接调用圆锥函数绘制出的效果很不理想(底面半径是等长的),如图6所示。为了能精细、准确地绘制护坡,本文设计了如下算法可完成护坡的建模。如图5所示(以左侧近似锥面护坡为例),A点为台前棱台底部端点(即台前棱台底部长度为L1),B点为左侧护坡底部端点(即左侧护坡底部长度为L2);求解由A、B形成的弧,并将弧AB进行等分,由O点与弧上相邻2点形成三角形来近似逼近锥面。其中弧AB的中心取线段AB的中垂线与x轴和y轴交点中距离最小的交点。对于L1大于L2如图5(a)所示,则弧AB的中心取在线段AB中垂线与x轴的交点C处;对于L1小于L2如图5(b)所示,则弧AB的中心取在线段AB中垂线与y轴的交点C处。整个左侧近似锥面护坡的参数由锥高(台后路基高度)zH、台前棱台底部长度L1及左侧护坡底部长度L2决定。1)由AB两点坐标可算出线段AB的方程为:y=-(L2/L1)*x+L2;2)由此可计算出线段AB的中垂线方程为:y-(L2/2)=(L1/L2)*(x-L1/2);3)当L1>L2时,如图5(a)所示,取此中垂线与x轴的交点为圆弧的圆心,故计算与x轴的交点坐标C(x1,y1,0),令y1=0,得到x1=L1/2-(L2*L2)/(2*L1);4)此时,得到所需绘制圆弧的半径r=L1-x1;5)计算弧AB的弧心角。由向量CA与向量CB夹角公式可得,其中θ为向量CA与向量CB的夹角。6)将圆弧平均分成n段,则第i点的坐标为(r*cos(θ/n*i)+x1,r*sin(θ/n*i),0),再以这些点中每相邻两点与O(0,0,zH)绘制三角形面片,则可得到逼近椎面。在上述第3步中,若L1小于等于L2时,如图5(b)所示,取其中垂线与Y轴的交点为圆弧圆心,计算与Y轴的交点坐标C(x2,y2,0),令x2=0,得到y2=L2/2-(L1*L1)/(2*L2);此时得到所需绘制圆弧的半径r=L2-y2;同样使用向量间夹角公式算出此时的夹角θ。绘制出的护坡三维模型如图7所示。3这座桥是连接3.1桥梁局部坐标系规则桥梁布局结构如图8所示,桥梁中心线是桥梁拼接中一个重要的参考位置,布孔线是两孔主梁的接缝线,也是本文中桥墩的默认位置。沿桥梁中心线纵桥向从左往右拼接桥梁,并定义桥梁中心线与世界坐标系的X轴重合。依据桥梁各组件空间位置拓扑关系进行整桥拼接,在桥梁一般的结构中,其上部结构即主梁是由跨径组成来描述的,每一跨又由结构相似的几孔组成,桥梁各组件的空间位置描述如下:左侧桥台,其局部坐标系原点与全桥的世界坐标系原点重合;第一孔主梁的左端位于左侧桥台的帽梁局部坐标系原点上;桥墩位于第一跨主梁与第二跨主梁之间下方;依次,从左往右,所有跨数主梁及桥墩拼接完,再绘制其右侧桥台,最后一孔主梁的最右端同样位于右侧桥台的帽梁局部坐标系原点。3.2培育几种独立的催化剂桥梁的5种主梁中空心板、简支变连续小箱梁、T梁是由若干块板梁横向拼接而成。本文采用如下算法完成了空心板的拼接,简支变连续小箱梁、T梁的拼接与此相同。各组件的局部坐标系原点位于左侧上部中心位置,而所拼接的空心板在将其拼接到整个桥梁上时,应确保其组件的局部坐标系原点位于板梁的左侧上端部中心位置,因此,当空心板块数是奇数或偶数时,其局部坐标系原点位置是有所差异的。由Num块空心板组成一孔的板梁拼接算法如下:1)判断板梁的块数Num是奇数还是偶数,若为偶数,则桥梁中心线位于Num/2号板与(Num/2+1)号板之间;若为奇数,则桥梁中心线位于(Num+1)/2号板的中心线上。2)进行循环控制(i=0;i<Num-2;i++,i为循环控制变量)以便绘制每一块中间板。若循环未结束则继续;若循环结束则进入6)。3)首先绘制第Num/2号中板,若Num为偶数,则其中心位置为pos(0,-bT/2,0);若Num为奇数,则其中心位置为pos(0,0,0)。bT为中间板或边板横截面上部长度。4)若Num为偶数,则循环绘制第Num/2号板的右边的第j块板梁,其中心位置为pos(0,-(bT/2+bT*j,0);若Num为奇数,则循环绘制第Num/2号板的右边的j块板梁,其中心位置为5)若Num为偶数,则循环绘制第Num/2号板的左边的k-1块板梁,其中心位置为pos(0,(bT/2+bT*(k-1),0);若Num为奇数,则循环绘制第Num/2号板的左边的k块板梁,其中心位置为pos(0,bT*k,0)。k++。6)若中间板绘制结束,则可绘制右边板。若Num为偶数,则其中心位置为pos(0,-((j+1)*bT+bT/2),0);若Num为奇数,则其中心位置为pos(0,-(j+1)*bT,0)。7)绘制左边板。若Num为偶数,则其中心位置pos(0,(k*bT+bT/2),0);若Num为奇数,则其中心位置pos(0,((k+1)*bT),0)。根据上述算法完成的空心板和小箱梁的拼接,如图9所示。3.3绘制右两侧桥台本文是基于桥梁整体纵向的位置分布约束建立桥梁的整体拼接关系。全桥拼接的算法如下:1)按照桥梁编号读取数据库中存储各桥梁组件的几何参数数据,并解析桥梁跨径组成信息,分别得到构成桥梁桥台数据的结构体向量、桥梁上部结构向量和桥梁下部结构向量;2)首先判断左侧桥台类型,若存储轻型桥台数据的结构体向量数组大于0,则调用绘制轻型桥台函数,若存储重力式桥台数据的结构体向量数组大于0,则调用绘制重力式桥台函数。左侧桥台局部坐标系原点定于帽梁上部中心位置,在全桥三维场景(世界坐标系)中位置为(-L/2,0,0),其中L是桥梁总长度;3)循环遍历上部结构向量数组,绘制每一跨。通过判断上部结构的类型,调用不同的上部结构绘制函数,如当前上部结构为空心板,则调