热力检查室内力计算方法研究.docx

中文摘要摘要：热力检查室板的内力可采用解析和有限元数值方法。解析计算又可根据板的支承条件分为简支、固支和弹性固定三种方法。有限元数值方法可以针对检查室整体进行结构内力计算。热力检查室设计中，板的内力计算一直以来采用解析方法，即把一个检查室分成六块单板，查表计算其弯矩，然后按纯弯构件配置钢筋。这样做没有考虑六块板的协调变形，也难以考虑人孔和沟口等开孔对板内力的影响。导致设计结果通常过于保守，造成钢材的浪费。而有限元数值方法得到的结果偏小，如果采用会使得工程安全储备不足。针对这一问题，本文开展检查室内力的数值计算方法研究，提出能够用于设计的简化数值计算方法。本文的主要内容和结果如下：(D把检查室各块板视为四边简支板和固支板，查表计算单块板的弯矩。(2)采用ANSYS软件的Shen63单元建立检查室模型，计算得到检查室内力的数值解。(3)提出检查室计算的框架方法，计算检查室的内力，并与以上两种方法得到的结果进行对比分析；根据对比结论对框架方法的荷载分配进行调整，得到了更合理的计算结果；再采用长宽比为2:1的检查室板验证框架方法的合理性。(4)把框架方法应用于实际工程案例，验证框架方法的实用性。通过解析解、数值解以及框架方法得到的解之间的对比可知：框架方法得出的内力解整体处于解析解和数值解之间。板四周简支条件下得到的板跨中弯矩是板中内力的上限值，数值计算得到的板跨中弯矩是板跨中弯矩的下限值；由板四周固支得到板支座弯矩是板端弯矩的上限值，数值模拟得到的板支座弯矩是板支座弯矩的下限值。本文的研究成果，对热力检查室结构及与其类似的地下结构的研究和设计具有较好的参考价值。关键词：热力检查室；板内力计算方法；解析方法；有限元数值模拟分类号：TU924ABSTRACTABSTRACTABSTRACT:Theoreticalmethodandnumericalmethodareoftenusedincalculationofinternalforcesofinspectionwells.Accordingtoboundaryconditions,aslabcanbesimplysupported,fixedandelasticallyfixed.FEMnumericalmethodcanbeusedtocalculateinternalforcesforaninspectionasawhole.Indesignofheatingprojects,theoreticalmethodismainlyadoptedincalculationofinternalforcesofinspectionwells,inwhichaninspectionwellisdividedintosixslabsandcalculatedseparately.Thentheslabsarereinforcedaccordingtothemoment.Indoingso,thecoordinateddeformationofsixslabsisignored,anditisdifficulttoconsidertheinfluenceofanopeningoraholetointernalfbrces.Asaresult,thedesignisusuallytooconservative,causingmaterialwaste.Ontheotherhand,numericalsimulationusuallygivesresultswhicharetoosmall.Ifadopted,thedesigncannotguaranteeenoughsafetymargin.Tosolvethisproblem,thispaperfocusesonanewmethodwhichcanbeusedindesignofheatingprojects.Themaincontentsandconclusionsareasbelow:(I)TakingawellassixSlabswithfour-sidessimplysupportedandfixedrespectiVelyJntemalforcesarecalculatedthereotically.(2)Shell63elementisusedandtheinternalforcesofwellsarecomputedbyANSYSsoftware.(3)Aframemethodisproposedandinspectionwellsarecalculatedbythismethod.Comparisonisconductedovertheresultsbytheoreticalmethod,numericalSimulationandframemethod.Uponthecomparison,theloadsontheframearemodified,whichleadstomorerationalresults.Thentwowellswithlength-to-widthrationbeing2:1arecalculated,whichshowsthattheframemethodcanalsogivecorrectresults(4)Theframemethodisusedtoanengineeringcaseandverified.Accordingtothecomparisononresultsbytheoreticalmethod,numericalsimulationandframemethod,theresultsbyframemethodliebetweentheoreticalmethodandnumericalsimulation.Forthemomentatcenterofasimplysupportedslab,theoreticalmethodgivestheIargestvaluewhilenumericalsimulationgivesthesmallest.Forthemomentatfoursidesofafixedslab,theoreticalmethodgivesthelargestvaluewhilenumericalsimulationgivesthesmallest.TheresultsofthispaperprovideagoodreferencefordesignandfurtherresearchOfinspectionwellsandotherundergroundstructures.KEYWORDS:inspectionWelhcalculationmethodforinternalforcesofaslab;theoreticalmethod;FEMnumericalsimulationCLASSNO:TU924目录中文摘要iABSTRACTiii第一章绪论11.1 弓I言11.2 热力检查室结构设计原理及存在问题11.3 国内外研究现状31.4 本文研究的背景及意义41.5 本文研究的内容和方法5第二章检查室内力查表计算72.1 概述72.2 典型检查室计算数据72.3 查表计算检查室板内力7第三章检查室内力数值模拟计算153.1 有限元分析的基本原理153.2 ANSYS软件中Shell63单元介绍173.3 检查室内力三维数值计算20第四章检查室内力按框架方法计算234.1 检查室框架计算方法的提出234.2 框架方法的内力计算234.2.1 框架荷载导算234.2.2 框架内力计算264.2.3 三种计算方法结果对比374.3 调整后框架方法的内力计算394.3.1 框架荷载导算394.3.2 框架内力计算424.3.4 调整荷载系数504.3.5 三种计算方法结果对比514.4 改变矩形面为1:2时的对比计算524.4.1 由解析计算得到的内力524.4.2 由数值模拟得到的内力554.4.3 由框架得到的内力564.4.4 三种计算方法结果对比分析59第五章工程案例615.1 工程概况615.2 传感器的布置625.3 监测结果635.4 用框架方法计算工程案例675.4.1 检查室计算数据675.4.2 查表计算检查室板内力685.4.3 检查室内力三维数值计算695.4.4 检查室内力平面框架解析计算705.4.5 双筋矩形截面偏心受压内力计算735.5 计算结果与监测结果的对比74第六章主要结论和展望756.1 主要结论756.2 展望75参考文献77作者简历79独创性声明81学位论文数据集83第一章绪论1.1 引言随着我国经济的飞速发展，快速的城市化进程对土建行业的需要和要求随之大大增加，在城市市政工程中，工程设计人员要面临集中供热的问题，对于城市市政工程，目前主要有明挖法和喑挖法。在我国北方地区集中供热在冬季供暖中应用的很广泛。与其它众多供热方式相比，其主要的优势在于：卫生、节能、环保。为了节省城市使用空间、不影响城市的美观，城市中的集中供热系统往往采用的是地下敷设方式。却也带来了一些新的问题。埋设在地下的管道一旦发生诸如个别管线泄漏事故，如何及时采取有效措施进行处理确保其正常工作？目前最为经济、有效的方法是在管道分支处和套管补偿器、波形补偿器、疏水器、阀门、排气、放水和除污装置等易于出现故障的管道附件处，设置专门用于管道维护和检修之用的特殊构筑物一一热力检查室，对管路进行定期的维护和检修1。与工程建设不相适应的是，热力检查室的内力计算一直以来采用解析方法，即把一个检查室分成六块单板，查表计算其弯矩，然后按纯弯构件配置钢筋。这样做没有考虑六块板的协调变形，也难以考虑人孔和沟口等开孔对板内力的影响。导致设计结果通常过于保守，造成钢材的浪费。而有限元数值方法得到的结果偏小，如果采用会使得工程存在安全隐患。针对这一问题，地下热力结构设计指南评审会中提出检查室需要采用更合理的结构计算方法。为此，本文开展检查室板内力的数值计算方法研究，提出能够用于设计的简化数值计算方法，为广大设计人员提供具有实用性的参考。1.2 热力检查室结构设计原理及存在问题近年来，在我国城市所新建的集中供热系统中，热力检查室采用的多为矩形的截面形式。与其它的截面形式相比较，矩形截面的检查室主要具有施工简便，易于与周围其它地面构筑物连接等优点。但是一直以来，由于各种不同的原因，对于这种矩形检查室的结构设计，却始终没有一种比较好的专用设计方法。在现今的检查室计算中，常用到两种简化模型：单元板简化模型和整室简化模型。单元板简化模型，是将检查室六块板当做相互独立的板，按照四边简支和固支分别计算板中和板端处弯矩，这种简化模型不用分析板与板的相互作用，不用考虑板与板之间的协调变形，所以计算方便，但由于忽略了检查室的整体变形，会产生内力不连续，变形不协调等问题。工程中通常在板与板连接处加上一定的加固措施，抵消整体计算产生的不利影响，从而保证检查室结果的安全使用。整室简化模型，是将检查室作为一个整体分析其受力情况的一种简化模型。此模型最突出优点是：模型中考虑了板间内力的相互影响作用，避免了检查室内力及变形的不连续性和板间交接处内力不平衡问题。但是，简化模型计算中需要考虑的影响因素较多，求解难度很大，不便于实践中直接应用。所以直到今天，对矩形检查室结构的设计，仍沿用着矩形水池一类矩形小室结构计算中所通用的那五种计算方法，即检查室的古典弹性计算法、力学计算法、弯矩分配法、有限差分法和有限元法。古典弹性计算方法，是一种最基本的热力检查室单元板简化模型的计算方法。它是以结构的平衡条件及应变协调条件为依据，由板的弹性曲面微分方程，即板的拉格朗日(Iagrange)方程，结合板的一定边界条件，来进行单板的内力及变形计算。力学计算法，是一种在复杂结构问题求解中所普遍用到的计算方法。这种计算方法基本的思想是，通过各种方式的简化、转换，将一个复杂的结构体模型转化为一系列相对简单的梁单元，然后应用结构力学的基本计算理论，进行结构的具体计算。弯矩分配法属于一种改进了的力学计算法。基本设计思想是：在相邻板的共同棱边上，首先加以人为约束，使板的边缘固定，此时，板的边缘必产生固端弯矩。由于两相邻板不同而产生的不平衡力矩有使板边缘产生转动的倾向。当放松人为约束时，由于在不平衡力矩作用下棱边发生转动而使两相邻板产生变形和新的分配弯矩。此分配弯矩是板抵抗棱边转动所产生的弯矩，等于将不平衡力矩按相邻板的劲度比反向分配于各板的弯矩值。各板边缘的最终弯矩等于分配弯矩与固端弯矩的代数和。有限差分法对于板面上有限量的点的挠度用一个联立线性代数方程组代替拉格朗日四阶偏微分方程进行求解。求得这些挠度后，再通过与弹性理论相似的方法，利用各种内力与挠度间的转换关系，求得板的其它内力。在应用联立线性代数方程组替代板的拉格朗日四阶偏微分方程时，需要用到有限差分算子的倒数以及内力。有限元技术是一种既适合于检查室的单元板简化模型计算，又适合于检查室整室简化模型计算的数值方法。这种计算方法与有限差分法最大的不同，在于它是将单元板分成许多矩形，三角形或四边形的面积或单元，而不是网格，所研究的是分割出的小面积而不是网格。每个单元都具有己知的或可精确地近似的弯曲变位的特性。一般的分析方法是将荷载集中作用在割出的单元的角部（或结点）上，然后再在每个结点上（有时也在中间的一些点上）重新恢复斜率和变形的连续性，以满足平衡条件和边界条件。在当今有限元技术的发展阶段中，有许多技巧因素渗入了科学方法之内。选择单元的形状、单元的数量、具体的变位特征、单元的定位、单元的数量等都会影响解的最终结果。有限元技术采用较简便的方法即可满足几乎任何边界条件，所以它是一个有力的分析工具。目前已有很多专门用于有限元运算的软件，例如：ansys,abaqus等。这些有限元软件能帮助我们分析工程应用中复杂的结构体系。五种计算方法各自都有优点和不足，但在实际工程中还有需要进一步解决的问题。在实际工程中，设计人员将检查室的各块板单独按照荷载结构模型计算，实际上，组成检查室的各板之间并非相互独立，互不影响的，各单元板在承载的过程中，向与之相邻的板面进行传递，从而使得检查室所受外部荷载在检查室内部产生了重分配，这种把各板当做独立板的计算方法计算结果过于保守，会造成钢材的浪费。而有限元数值方法得到的结果偏小，如果采用会使得工程存在安全隐患。所以热力检查室板的内力计算方法有待于我们进一步探讨。1.3 国内外研究现状对于城市地下热力结构工程，目前国内外都有相关的工程实践，但具体的规范性的还很稀缺，多是借鉴相关专业的设计规范和设计经验。对于矩形检查室板间连接及其相互影响的研究是一个十分复杂的问题。在检查室中各块板的刚度不一，研究中需要考虑组成检查室各板件的材料性质、各板件间钢筋的连接形式以及各板面上荷载的作用情况等诸多难以确定因素的影响，所以一直以来国内外对此问题的研究都很少。工程中往往将各块板从检查室中简化出来单独计算，一块完整不开孔的板，在大多数支撑条件下，其弯矩可按建筑结构静力计算手册通过查表算出，但这一解析法的缺陷有二：不能考虑有些支撑条件；内力只有弯矩，没有杆件的轴力；而在实际检查室中，还要处理各板件连接的问题，可以通过对板间相交处的内力的近似重分配（或按板的线刚度分配，或按板的劲度分配），来解决相邻板面上弯矩的不平衡问题，以及要求在实际结构配筋中，通过一些构造上的规定来抵消板间作用对整室结构可能带来的一些不利影响2O目前对检查室的内力-一般采用解析方法查表计算。没有考虑六块板的协调变形；没有考虑杆件的轴力；对开孔的影响采用局部加强的办法进行处理。对地下检查室这样类似的构筑物，如地铁车站，除按单板解析计算外，有些设计单位开始采用平面框架有限元数值方法计算。复杂情况下，进行三维有限元计算以作设计参考。但是，地下检查室简化为平面框架，一般是在长宽比在2倍及以上是合适的，短于这个尺寸，平面框架上作用的荷载如何确定，需要深入研究。为了解决实际问题，有时规范会对实际情况做相应的简化处理，如在给水排水工程钢筋混凝土水池结构设计规程中规定3,当顶板为预制装配板搁置在池壁顶端而无其它连接措施时，顶板应视为简支于池壁，池壁顶端应视为自由端当预制顶板与池壁顶端有抗剪钢筋时，该节点视为被支承；当顶板与池壁为整体浇注并配置连接钢筋时，该节点视为弹性固定；当仅配置抗剪钢筋时，该节点视为钱支承；池壁与底板、条形基础或斗槽连接，可视池壁为固端支承；对于软基上的水池，应考虑地基变形的影响，宜按弹性固定计算；当池壁为双向受力时，相邻池壁间的连接应视为弹性固定。目前有限元数值计算方法应用广泛，其突出优势是能够考虑不同材料、不同几何尺寸（包括板的开孔）、不同支撑条件及不同荷载形式等复杂条件。因此，数值计算方法在实践中使用越来越多。但是，把数值计算的结果直接用于设计的却不多。一般认为其计算结果偏小，不能保证安全。热力事业在最近几十年得到蓬勃发展，与工程建设不相适应的是，目前的热力结构相关设计规范和计算方法无法满足需求，规范仅涉及到明挖的管沟与检查室，对暗挖及较大埋深结构没有明确规定，尤其对热力地下结构工程的设计，尚没有专业的设计规范和行业内专家认同的一些经验数据。目前热力地下结构设计还是参照城镇供热管网结构设计规范及地下铁道设计规范等相关规范，所以急需要编制地下热力工程的结构设计规范或指南，从而规范各种结构的计算模型，对计算方法进行分析，在实验和分析的基础上对各种计算参数进行合理的取值，使地下热力工程的结构更加安全、经济和合理。1.4 本文研究的背景及意义本文的研究是以北京市热力工程设计有限责任公司与北京交通大学合作项目热力检查室内力计算方法研究的一部分，其目的是为热力检查室板结构设计提供相应的技术依据。本文针对热力检查室的内力计算方法这一问题，提出了一种解决检查室板内力的新途径，经验证表明所计算出来的值在解析解和数值模拟解之间，通过工程实例也得到印证，证明提出方法的可行性，为检查室板内力计算提出了一种更合理的结构计算方法，可以为工程设计人员在设计中作为参考。本文对热力检查室结构设计方法的研究，有着较强的理论价值和现实意义。1.5 本文研究的内容和方法本文的主要内容包括：(1) 对典型检查室进行查表计算，得到各个检查室内力的解析结果。(2) 对典型检查室进行三维有限元数值模拟，得到各个检查室内力的数值模拟结果。(3) 提出检查室闭合框架计算方法，用此方法对典型检查室进行计算，得到各个检查室的内力。(4) 对以上三种方法得到的计算结果进行分析对比，论证方法的可行性。(5) 用本文提出的闭合框架计算内力方法与北京郑常庄热电厂配套热力管线工程(万寿路一西三环)的监测结果相比较，验证方法的可行性。本文中主要用到的研究方法(1) 在各典型检查室板的内力计算中，分别对各块板进行了解析计算。对于板中的内力依照四边简支的简支板计算，对于跨端的内力依照各板为四边固支板计算，计算方法参考结构静力计算手册。(2) 利用ANSYS软件对检查室三维整体进行了数值模拟。(3) 再用本文提出闭合框架方法算出各典型检查室的板内力，最后对三种计算结果进行对比分析。(4) 分析北京郑常庄热电厂配套热力管线工程(万寿路一西三环)的监测结果，与本文提出的方法计算结果进行对比分析，分析其合理性。第二章检查室内力查表计算本章主要包括如下内容：(1)选取用于对比计算的典型检查室，确定其尺寸和荷载。(2)根据建筑结构静力计算手册，把检查室板视为四边简支和固支板，分别查表计算单块板的内力。2. 2典型检查室计算数据选择的3个典型检查室进行对比计算，其几何尺寸及所承受的荷载如表2.1所示。表2.1检查室的计算数据检查室编号检查室内净空(m)顶板、底板、边境厚(m)顶板竖向土压力(kPa)墙顶侧向土压力(kPa)墙底侧向土压力(kPa)长宽高恒载活载合计恒载活载合计恒载活载合计18880.5114.3025.20139.5041.918.4050.31102.878.40111.27212880.5114.3025.20139.5041.918.4050.31102.878.40111.273880.5114.3025.20139.5041.91S.1(150.31133.358.40141.75混凝土强度等级C30,容重取25kNm30结构自重取结构体积X混凝土重度。底板反力计算如下：检查室1：P=(9×9×0.5×2+8×8×0.5×2+8×9×0.5×2)×25(9X9)+139.5=66.9753+139.5=206.48kPa检查室2:P=(13×9×0.5×2+8×8×0.5×2+8×13×0.5×2)×25(13×9)+139.5=60.90+139.5=200.40kPa检查室3:P=(9×9×0.5×2+8×12×0.5×2+12×9×0.5×2)X25/(9X9)+139.5=87.96+139.5=227.46kPa3. 3查表计算检查室板内力顶板计算时要加上混凝土的自重，底板计算时要减去混凝土的自重。顶底板混凝土自重为：25×0.5=12.5kPa0顶板的荷载为：139.50+12.5=152kPa。底板的荷的值为:检查室1:206.48-12.5=193.98kPa;检查室2:200.4072.5二187.90kPa;检查室3:227.46-12.5=214.96kPa°列出顶板、底板、侧墙计算用荷载如表2.2。表2.2检查室解析计算弯矩用荷载查编检室号检查室内净空(m)顶板、底板、边墙厚(In)顶板竖向土压力(kPa)墙顶侧向土压力(kPa)墙底侧向土压力(kl,a)底板反力(kPa)净长净宽净高18880.515250.31111.27193.98212880.515250.31111.27187.903881/0.515250.31141.75214.96(2)侧墙梯形荷载的处理侧墙所受荷载可分解成矩形均布荷载和三角形荷载，如图2.1所示。侧墙弯矩即为矩形均布荷载和三角形荷载作用下弯矩的线性叠加。111.27kpa50.31kpo60.96kpo图2.1梯形荷载的等效均布荷载最大弯矩在板中间，但三角形荷载最大弯矩离板中间有一点距离，故最大弯矩值不一定在板中间。对于长8.5m的简支梁，均布荷载50.3IkPa最大弯矩在梁中心(4.25m处)，弯矩值为454.36kNm;三角形荷载60.96kPa最大弯矩在4.9Im处，弯矩值为282.55kN-m,梁中间4.25m处弯矩值为275.27kNm。叠加起来，最大弯矩在4.5Im处，为452.66+279.83=732.49kN*而中间点弯矩为454.36+275.27=729.63kNm0729.63/732.49=99.61%,误差在1%以内。故可近似认为中点即为梁的最大值点，推广到板上，均布荷载和三角形荷载在板中间的叠加值近视认为是板的最大弯矩值。(3)四边简支板计算把检查室各块板当成简支板，查表计算跨中两个方向的弯矩。计算公式为M=kpL2式中M跨内或支座单位板宽内的弯矩，kN.mok弯矩系数，由查表得到。P板上的均布荷载，kPa。L计算跨度，取板尺寸LX和Ly中的较小值。建筑结构静力计算手册表中跨内弯矩系数系泊松比U=O的计算用表。当UWo时，跨内弯矩需进行如下换算：MM=%+5查表得3个检查室的弯矩系数，见表2.3。表2.3四边简支板的弯矩土卜算系数检查室编号检查室部位均布载荷三角形荷载MxMyMx不1顶板0.03680.0368-底板0.03680.0368-侧墙0.03680.03680.01840.01842顶板0.02860.0710底板0.02860.0710前后墙0.02860.07100.01430.0355左右墙0.03680.03680.01840.01843顶板0.03680.0368底板0.03680.0368侧墙0.07100.02860.03550.01433个检查室的弯矩计算结果见表2.4。表2.4四边简支板的弯矩检查室编号检查室板计算尺寸(m)计算结果(kNm)顶板底板前后墙左右墙MxMyMxMyMxMyMxMy18.5X8.5X8.5484.97484.97618.90618.90257.77257.77257.77257.77212.5×8.5×8.5470.03842.54581.041041.53249.83447.82257.77257.7738.5×8.5×12.5484.97484.97685.84685.84532.30296.95532.30296.95(4)四边固支板计算把检查室各块板当成固支板，查表计算跨中及固端的弯矩。当U二O时，四边固支板跨中计算系数如表2.5。表2.5四边固支板跨中的弯矩计算系数检查室编号检查室部位均布载荷三角形荷载MxMyMxMy1顶板0.01760.0176底板0.01760.0176侧墙0.01760.01760.00880.00882顶板0.01060.0331底板0.01060.0331前后墙0.01060.03310.00530.0165左右墙0.01760.01760.00880.00883顶板0.01760.0176-底板0.01760.0176侧墙0.03310.01060.01650.0053算出四边固支板跨中弯矩如表2.6o表2.6四边固支板跨中弯矩检查室编号检查室板计算尺寸(m)计算结果(kN-m)顶板底板前后墙左右墙MxMyMxMyMxMyMxMy18.5×8,5×8.5231.94231.94296.00296.00123.28123.28123.28123.28212.5X8.5X8.5189.11386.79233.77478.14100.47205.36123.28123.2838.5X8.5X12.5231.94231.94328.01328.01244.03119.41244.03119.41当0时，固端弯矩不变，不用修正。固端计算系数如表2.7。表2.7四边固支板固端弯矩计算系数检查室编号检查室部位均布载荷三角形荷载MxMyMxMylMy21顶板-0.0513-0.0513底板-0.0513-0.0513侧墙-0.0513-0.0513-0.0257-0.0334-0.01792顶板-0.0570-0.0747底板-0.0570-0.0747前后墙-0.0570-0.0747-0.0284-0.0457-0.0291左右墙-0.0513-0.0513-0.0257-0.0334-0.01793顶板-0.0513-0.0513一底板-0.0513-0.0513一-侧墙-0.0747-0.0570-0.0374-0.0418-0.0152算出四边固支板固端弯矩如表2.8o表2.8四边固支板固端弯矩编号检查室板计算尺寸(m)计算结果(kNm)顶板底板前后墙左右墙MxMyMxMyMxMylMy2MxMylMy218.5X8.5X8.5-563.38-563.38-718.97-718.97-299.66-333.58-186.47-299.66-333.58-186.47212.5×8.5×8.5-625.97-820.36-773.82-1014.11-332.27-472.81-271.53-299.66-333.58-186.4738.5×8.5X12.5-563.38-563.38-796.73-796.73-518.61-483.34-207.19-518.61483.34-207.19(5)检查室轴力的计算侧墙轴力为顶板传过来的压力和侧墙自重产生，而顶板传过来的压力按双向板的受力方式传力；顶底板轴力由侧向荷载产生，计算方式也把侧墙按双向板的受力方式传力，按受力面的中点处荷载值作为计算的平均值，作用示意图如图2.2。对于检查室一1)侧墙竖直方向顶部轴力为152X9X9/4/8.5:362.12kN/m底部轴力为362.12+0.5X8X25=462.12kN/m水平方向轴力为(50.31+111.27)/2×9X9/4/8.5=192.47kNm2)顶板左右方向轴力为(3/4X50.31+14×111.27)×9×948.5=156.16kNm前后方向与左右方向相同轴力为(3/4X50.31+1/4X111.27)×9×948.5=156.16kNm3)底板左右方向轴力为(1/4X50.31+34×111.27)×9X9/4/8.5=228.78kNm前后方向与左右方向相同轴力为(1/4X50.31+34×111.27)×9×9/4/8.5=228.78kNm对于检查室二1)前后侧墙竖直方向顶部轴力为152X(2×13-9)X9/4/12.5=465.12kN/m底部轴力为465.12+0.5×8×25=565.12kN/m水平方向轴力为(50.31+111.27)2×9×948.5=192.47kNm左右侧墙竖直方向顶部轴力为152X9X9/4/8.5=362.12kN/m底部轴力为362.12+0.5×8×25=462.12kN/m水平方向轴力为(50.31+111.27)2×9×948.5=192.47kNm2)顶板左右方向轴力为(3/4X50.31+14×111.27)×9×9/4/8.5=156.16kNm前后方向轴力为(3/4X50.31+14×111.27)×(2×13-9)×9/4/12.5=200.58kNm3)底板左右方向轴力为(1/4X50.31+34×111.27)×9×9/4/8.5=228.78kNm前后方向轴力为(1/4X50.31+3/4×111.27)×(2×13-9)X9/4/12.5=293.85kNm对于检查室三1)侧墙竖直方向顶部轴力为152X9X9/4/8.5=362.12kN/m底部轴力为362.12+0.5X12X25=512.12kN/m水平方向轴力为(50.31+141.75)2×(2×13-9)X9/4/12.5=293.85kNm2)顶板左右方向轴力为(0.827X50.31+0.173×141.75)×9×9/4/8.5=157.54kNm前后方向与左右方向相同轴力为(0.827×50.31+0.173×141.75)×9X9/4/8.5=157.54kNm3)底板左右方向轴力为(0.173X50.31+0.827×141.75)×9×9/4/8.5=300.01kNm前后方向与左右方向相同轴力为(0.173X50.31+0.827X141.75)×9×948.5=300.OlkNZm得到各板的轴力表2.91.1 9各检查室板的轴力值检查室编号检查室内净空(m)计算结果(kN/m)顶板底板前后墙左右墙左右向前后向左右向前后向水平向竖向水平向竖向18X8X8156.16156.16228.78228.78192.47412.12192.47412.12212×8×8156.16200.58228.78293.85192.47515.12192.47412.1238X8X12157.54157.54300.Ol300.Ol293.85437.12293.85437.12第三章检查室内力数值模拟计算1.2 有限元分析的基本原理有限单元法(FEM,FiniteElementMethOd)是随着电子计算机的发展而发展起来的一种很有效的数值方法。1960年，美国的CknlghR.W.在一篇论文中首次使用“有限元法”这个词。它的创立和应用在工程分析中具有重要的意义。在国外，20世纪70年代初有限单元法的基本的理论己基本成熟，商业的有限元分析软件开始陆续出现。在我国，60年代初，数学家冯康教授和他的研究组提出了一种以变分原理为基础的三角形近似法，并用偏微分方程求得了近似解，在严密的数学基础上证明了它的收敛性、稳定性和误差估计，编制了计算程序，并成功的应用于刘家峡水坝的应力分析等一些课题一一这个方法就是我们通常应用的有限单元法。有限单元法是在计算机技术和数值分析方法支持下发展起来的，为解决复杂的工程分析计算问题提供了有效地方法。有限单元法的主要优点是它应用范围相当广泛，几乎在所有的连续介质和场问题中都得到了应用，有限单元法的基本思想是将问题的求解步骤划分为一系列的单元，单元之间仅靠节点连接。单元内的节点