大型地下钢筋混凝土岔管结构优化分析2013113111055460.docVIP

大型地下钢筋混凝土岔管结构优化分析 1 地下岔管结构优化分析有限元网格自动生成方法 　　地下钢筋混凝土岔管结构是一种复杂的多重介质空间结构，采用有限元数值方法进行优化分析时，首先须对地下钢筋混凝土岔管结构进行离散化。复杂空间结构的三维有限元网格离散剖分工作是一项繁杂费时的工作，而且数据输入时容易出现人为差错。虽然国内有一些著名软件如ViziCAD，AutoCAD具有较好的交互式绘图功能，但是这些软件往往需要有限元使用者有一定的建模技巧，同样须花相当多的时间。当进行岔管结构优化分析时，每分析一个方案，都必须重新建模，对修改结构方案进行重新离散和重新组织数据，所以对优化参数的选择和优化结构的判断十分不便[2]。因此编制一个适应于各种不同类型的地下钢筋混凝土岔管有限元网格自动生成程序，对岔管进行优化分析是十分重要的。 　　地下岔管钢筋混凝土岔管，衬砌与围岩共同承载，因此网格剖分应包括衬砌和围岩。为了适用于各种地下岔管衬砌与围岩联合承载的三维有限元网格自动生成，本文利用映射法原理编制岔管结构网格的自动生成程序。只要用户给出结构轮廓的某些特征值，程序根据这些特征值，调用相关算法，即可进行网格的自动生成，不再需要用户计算单元信息(包括各种约束信息、锚杆信息、受载面信息).因此大大减小了前处理工作量，节省了时间，降低了出错率。使得地下岔管结构优化分析计算变得简捷而快速。程序生成的有限元数据可以通过图形显示程序显示出来，并可通过人机交互交式，对生成的网格进行修改、添加、删除，保证有限元分析数据的准确性。 　　地下钢筋混凝土岔管有限元网格自动剖分程序主要分3步进行。(1)根据钢筋混凝土岔管结构的几何特征和有限单元信息规定，首先对规则岔管结构进行离散。如图1(a)对称岔管，只需输入主、支管直径D1，D2；衬砌厚度d1,d2；主管和总管长度S1和S；岔管的半岔角θ、半锥角β和两支管的间距L，则可分别建立主管和支管的轮廓方程： (1) 　　式中：X，Y，Z为以岔管相贯线交点为圆点的三维坐标；C0，C1，C2为沿X轴线所取的不同段面常数。对于岔管段可根据岔角参数(半岔角θ、半锥角β)建立以腰梁和U梁为局部坐标的相关方程，通过数学变换则可求所需相贯线的坐标。再根据几何关系，则可推出岔管锥段斜截面的交线参数方程[3]。如图1(b)所示，任意倾角为α1的斜截面π与母线OD的交点为P，由图中的几何关系则可得到P点在截面上的向径OP和夹角α的表达式为： ρp=cos(β+α1)cos(β-α)/cos(β+α1-α)cosβ，tgα=(1-cosφ)sinβcosβ/cos2β+cosφsin2β (2) 式中：α为母线OD在对称面OAB上的投影OC与OA的夹角(∠AOC)；φ角为母线OD的旋转角，即从OA开始旋转到OD的角度∠AED的角度。 根据O点的坐标和精度要求控制，选取一定间距截面和不同的φ角，按照式(2)则计算出斜截面π上的任意交点P的坐标。对于卜型岔采用类似方法可求得岔管轮廓线尺寸和相贯线坐标。 图1 对称岔管结构参数 图2 对称岔和卜型岔衬砌单元 图3 对称岔整体有限元单元网格 (2)对岔管外围的围岩进行网格剖分。首先可根据已离散的岔管衬砌单元和节点，按其规则将围岩向外扩充3～4层单元，为此可采用对岔管衬砌外边界进行自动标识，而后将其向围岩拓展。然后根据岩体所要计算的边界范围以及网格离散的计算精度确定计算模型边界的单元编号和结点坐标。再采用插值方法对中间岩体网格进行逐步加密过渡，则可确定岔管外围的全部岩体单元网格。这一部分关键在于对岔管岔角处网格剖分，因为岔角处应力较为集中，往往是有限元分析的重点，要求网格较密，同时又要求与第一部分网格相协调。因此在对岔管外围规则部分岩体网格进行离散时，可采用管道截面交线逐渐由圆或椭圆过渡到矩形，从而使得岔管岔角处围岩形成为一个大的三棱柱，对于较为规则的三棱柱，可以较为容易地进行六面体网格剖分。 (3)对岔管分岔角处进行修圆处理。地下钢筋混凝土衬砌岔管分岔角处的应力集中现象十分明显，在开挖时，常常对分岔角处进行修圆处理。岔角修圆后过渡段不再是一锥管，但截面交线参数方程仍为椭圆。本文采用在分岔角处进行滚动圆滑处理，即根据修圆程度要求，选用一个合适半径的圆球，将岔管过渡段母线滚动圆滑，求得母线与斜截面的交点。然后根据修圆的控制范围选定不修圆的边界，再对修圆后的母线和不修圆的边界之间的网格节点坐标进行插值计算，则可得到较均匀的修圆岔角网格。采用上述方法，只需修改初始的几个输入参数，则可离散如图2、图3所示各种岔管参数网格。