水利工程论文-太浦河泵站泵房底板空间和平面有限元结构分析及其比较研究.doc

水利工程论文-太浦河泵站泵房底板空间和平面有限元结构分析及其比较研究摘要：软土地基上基础板的结构分析通常简化为平面问题进行，这种简化与实际情况存在较大差异。了解这种差异特性，对于正确评估简化计算的结果，保证设计的安全和经济具有重要意义。太浦河泵站泵房底板是典型的大型软土地基基础板。本文以该底板为计算对象，分别应用有限单元法进行了空间三维和简化的平面二维结构分析，分析中考虑了底板、基础、墩墙和顶板的共同作用。通过对计算结果的比较分析，找出差异上的一些规律，为类似结构的计算提供一定的参考。关键词：太浦河泵站软基底板结构分析有限元ANSYS程序3D/2DFEAandComparativeStudyofFoundationSlabofTaipuhePumpingStationAbstract:Structureanalysisoffoundationslabonsoftfoundationisusuallysimplifiedtoplaneproblem.Suchsimplificationhasquitemoredifferencetotheact.Tofindoutthisdifferenceissignificativetocorrectlyevaluateresultsandensuredesignsafetyaswellaseconomy.ThefoundationslabofTaipuhepumpingstationisatypicalhugefoundationslabonsoftfoundation.Inthisthesis,accurate3Dandsimplified2DFEAisdoneseparatelyontheslab,andslab-ground-upperstructureinteractionistakenintoaccountintheanalysis.Bycomparativeanalysisingtheresults,someconclusionisdrawn,whichcanbereferenceinsimilarstructureAnalysis.Keywords:TaipuhePumpingStationSoftfoundationSlabStructureanalysisFEAANSYSprogram1、问题的提出软土地基上水闸、泵站底板的结构计算通常是采用“截条法”进行，即沿着水流方向截取单位宽度的底板简化为横向（垂直水流方向）弹性基础梁或框架进行计算，实质上是把空间问题简化为平面应变问题进行计算。“截条法”在计算理论上存在着诸多缺陷：首先，水闸、泵站的底板在几何上是一个复杂的空间结构，荷载的分布往往也很复杂，几何特征和荷载特性均与弹性理论中的平面应变条件有较大的差别。其次，“截条法”中引进的不平衡剪力及其分配，其实只是为了满足结构的静力平衡而采取的一种“补偿”，实际上并不能反应结构真实的受力情况。再次，“截条法”在计算不平衡剪力的时候假设地基反力顺水流方向直线分布也不符合实际情况。因此，采用“截条法”进行底板结构分析，其计算结果与实际情况会产生较大的误差。这种误差与底板的几何形状和荷载分布有关。对于水闸底板，由于其形状和荷载相对比较匀称和简单，“截条法”的计算结果在某些条件下尚可以接受，而对于大型泵站底板，由于涉及到机电设备布置和流道的影响，其几何形状和荷载分布相当复杂，对于这样的空间体系，如果简化为平面问题分析，误差将更显著。按照空间问题进行底板、地基和上部结构的联合受力分析，可以比较精确地求得其应力和变形情况。随着计算机技术的发展和具备强大运算、前后处理功能的大型有限元分析系统出现，对这样的结构进行比较精确的分析成为了可能。由于这种分析方法工作量的巨大，在工程实践中，根据现行水闸和泵站设计规范的推荐，通常还是用“截条法”简化为平面问题分析，在这种情况下，如何正确评估这种简化计算所导致差异的，对于保证设计的安全和经济具有重要意义。这种评估应该包括以下几个内容：从定性的角度比较两种计算方法所求得的应力和变形的分布规律的异同；从定量的角度对应力和内力进行比较；比较平面计算中选择不同的计算断面对结果的影响等。本文以太浦河泵站泵房底板为计算对象，分别进行了空间三维和平面二维有限元分析，通过对这两种计算模型结果的比较，分析了它们之间的异同，对平面“截条法”的简化计算结果做出评价。2、工程概况和基本资料太浦河泵站工程位于江苏省吴江市东太湖边，工程的主要功能是在枯水年份抽取太湖水满足上海市的供水要求。泵站内安装6台单机流量为50m3/s的斜15。轴伸泵。泵房为堤身式块基型结构，分三块底板，每块底板平面尺寸为40.45m(顺水流方向)22.50m（垂直水流方向），底板平均厚度2m，廊道部位为1.55m。底板以下8m深度范围内地基采用水泥土搅拌桩复合地基。1.45m高程以下的底板、墩墙和流道体型参见图1。本文以完建期中跨底板单元为计算对象。3、空间有限元分析（1）建模和求解取中跨底板单元完建工况为计算对象，对1.45m高程以下包括顶板、墩墙、底板和地基进行整体空间三维有限元分析。地基深度取50m，水平方向各往外取22.5m。混凝土结构部分采用能够较好模拟弯曲变形的空间非协调单元，底板、流道和墩墙为C20混凝土，弹性模量为2.55104MPa，顶板为C25混凝土，弹性模量为2.55104MPa。地基部分采用空间协调单元，底板以下8m的按复合地基考虑，其变形模量为55MP，其余为天然地基，变形模量为10MP。将上部荷载导算到1.45m高程。计算底板单元的两侧作用有相邻底板的边荷载，按平均基底压力165kPa考虑。采用国际著名的有限元分析系统ANSYS进行三维有限元建模、分析和后处理。有限元模型见图1，水流方向为X，底板高度方向为Y，垂直水流方向为Z。图1三维有限元模型（2）计算结果分析底板竖向位移（图2）：进水侧小，出水侧大，在垂直水流方向则比较均匀。最大竖向位移发生在出口处，为100mm，廊道水泵支墩处的沉降为92mm，顶板在电机部位的沉降为97mm。整块底板的计算不均匀沉降差为14.5mm。顺水流方向应力（图3）：底板进水流道段上部为压应力区，平均值为0.27MPa，沿垂直水流方向贯穿底板上部，在墩墙下部压应力区范围比较大。底板进水流道段下部为拉应力区，平均值为0.11MPa，沿图垂直水流方向贯穿底板下部，流道中部拉应力区比较大。底板出水流道段除了墩墙下部出现几个不连续的压应力区外，基本上是拉应力区，拉应力平均值0.11MPa；底板廊道段几乎全部是是拉应力区，只是底部靠近边缘处出现小范围的压应力区。顺水流向的最大拉应力为0.50MPa，出现在三个部位：廊道上表面、进水流道出口段下表面、出水流道中段下表面。垂直水流方向应力（图4）：整个底板上面受拉，下面受压。廊道及其附近为较高应力区，最大拉应力0.54MPa出现在进水流道出口段和出水流道进口段上表面，最大压应力0.83MPa出现在进口流道后部的下表面。图2竖向位移等值云图图3顺水流向应力云图图4垂直水流向应力云图4、平面有限元分析（1）计算断面的选择根据泵房底板不同部位的支承形式、跨度和受力条件，进、出水流道按弹性地基上的框架、机组段按弹性地基上的梁计算。选择计算剖面的时候，为了比较不同的上部荷载对结果的影响，分别在进出流道段选择包含排架柱、混凝土承重墙等承受上部荷载的断面以及不承受上部荷载的断面。选择的6个计算断面为：C断面（X=4.0m，包含进水侧胸墙）、D断面（X=7.225m，无上部荷载）、E断面（X=10.45m，包含进水侧上部排架柱）、F断面（X=13.75m，包含水泵支墩）、I断面（X=31.45m，包含出水侧上部排架柱）、J断面（X=33.5m，无上部荷载）。各断面位置见图1。（2）计算荷载的确定假设基础反力沿水流方向直线分布，根据偏心受压公式求得底板反力。对于给定的计算断面，地基反力、自重和上部荷载是确定的，据此可以求出不平衡剪力。根据剪应力分布图形的面积比将该不平衡剪力分配给顶板、墩墙和底板，并作为分布力施加。底板、墩墙和顶板分担的不平衡剪力大概为15%，70%和15%左右。由于地基反力在顺水流向的分布相对比较均匀，而上部荷载分布往往具有很大的突变性，所以对于上部荷载值大的计算断面，不平衡剪力的方向是向上的，反之不平衡剪力的方向向下。（3）平面有限元计算为了便于与空间有限元计算结果的比较，同时避免梁截面大高跨比所引