高填土公路箱涵垂直土压力有限元模拟分析.docVIP

高填土公路箱涵垂直土压力有限元模拟分析 　　摘要 利用大型有限元软件ANSYS，对高填土箱形涵洞进行施工过程模拟，分析涵洞土压力的变化规律，讨论沟槽宽度、边坡系数及涵洞截面几何尺寸对洞顶垂直土压力的影响。 　　关键词:箱形涵洞;土压力;有限元;模拟分析 　　中图分类号：U412.1文献标识码： A 　　 　　随着我国公路大力发展，各类涵洞广泛应用于公路工程中。为保证公路排水畅通，我国《公路设计规范》中要求每公里设置不少于四道涵洞，尤其我山西省地处黄土高原，大部分是山岭区，大填大挖工程量较大，其它领域如水利、市政建设及能源等领域也在应用涵洞进行排水。洞顶垂直土压力为涵洞的主要静荷载，尤其在高填方涵洞结构设计中起控制作用。 　　对于沟埋式涵洞垂直土压力的计算,上世纪初，Iowa 州立大学的A.Marston、M.G.Spangler 等人对此进行了一系列的研究。以后许多学者的相关研究都是在此基础上进行的。目前沟埋式涵洞土压力的计算方法仍然不完善，《水工建筑物荷载设计规范》中不推荐相应的计算公式。本文采用大型有限元软件ANSYS，对沟埋式涵洞土压力进行模拟计算，分析沟槽宽度、边坡系数及涵洞截面几何尺寸等对洞顶垂直土压力的影响。 　　1.有限元模拟 　　1.1 模型及参数 　　国内外管道施工技术规定，管周回填土的夯实要求达到其最大夯实密度的95％ ，而在一般填土条件下，管周土压力都处于非极限状态。土是非线性压缩体，其变形模量很大程度上随应力状态而变化，但当土体应力处于非极限状态时，线性与非线性两种理论对土体垂直压应力的分析结果十分接近。本文将涵洞周围土体和地基作为弹性介质。填土材料选用中密类碎石砂土，内摩擦角φ=30°,弹性模量Et=1.5×107Pa,容重γ=18kN/m3；地基弹性模量Ed=7.5×108Pa；对于钢筋混凝土涵洞，也作为弹性体，弹性模量E h=3×1010Pa，泊松比μ均取0.167。模拟单元选用平面四节点Plane42类型单元。 　　考虑到管土的共同作用，在土与涵洞以及土与地基之间设置接触对。与填土接触的涵洞表面和地基面作为目标面，目标单元选用Targe169；与涵洞表面和地基面接触的填土面作为接触面，用Conta171和Conta172作为接触单元。根据库伦土压力理论 , 填土与沟槽边壁、涵洞结构的接触单元摩擦系数分别取0.45和0.3。 　　考虑实际工程中的边界条件、涵洞结构及填土加载过程等特点，将计算模型作为平面应变问题。由于结构与荷载的对称性，取几何实体的一半进行计算，在对称轴上施加对称边界，计算范围左右两侧及地基底部边界约束所有方向位移。模拟时，取高h和宽D都为2m的箱涵为例。 　　1. 2施工过程模拟 　　在有限元整体计算中, 回填土体是一次性施加全部荷载，各层土体的重量由整个结构(计算域内所有单元) 承担；而在有限元施工模拟计算中,计算域荷载是逐渐增加的，两者在承载机理和形成的位移场上完全不同。土石坝、隧洞、边坡及涵洞等工程建设中具有明显的施工特点，荷载是随着工程进度逐级施加到结构物上的，它们的有限元分析宜采用模拟计算。本文使用ANSYS的荷载步多重计算功能和单元生死技术来实现土体分层填筑的施工过程模拟,洞顶填土高度H最大取到10D，分20层填筑。 　　2.结果及分析 　　2.1填土竖向位移分析 　　模拟计算时,填土压缩变形在填筑过程中瞬时完成，不考虑其固结等时间因素的影响。也就是说，第i层土填完后，其自重由自身及其以下各层已填筑的土体承担，它们的自重所引起的位移已经完成,对其上部填筑的土体位移毫无影响，即这部分土重不再引起后续填土的位移，最上部填土的沉降位移较小；模拟地基刚度相对于土体刚度大得多，在填土作用下，下边界沉降位移很小（最大位移约1.5mm），域内各计算点的沉降位移都由后续各填筑层土重对其产生的位移增量值累加而成。对于最下部土体，虽然所受上部土重相对最大，但其可压缩高度相对较小，位移较小；填筑高度中部附近的土体，相对上部和下部土体而言，所受土重和可变形土体高度都较大，累积沉降位移最大。 　　当洞顶填土高度较小时，回填土在自重作用下沉降，涵洞刚度大，其上土体沉降小；边壁的摩擦使其附近土体沉降也小；土体刚度小，胸腔土体的土柱高，压缩量大，其沉降位移大。此时，填土最大沉降位移点在胸腔上方，即涵洞的侧上方。填土高度较大时，槽壁与其附近土体的接触范围扩大，这部分土体所受槽壁向上的摩擦力越大，沉降位移相对较小；沟槽中心距边壁最远，所受边壁向上的摩擦最小。随填土高度的增大，最大沉降位移点由原来的洞顶侧上方逐渐转移到沟槽中心。以沟槽底宽B与涵洞宽D之比B/D=5，沟槽边壁系数m=0的沟埋式涵洞为例，图1为施工模拟计算时回填土体垂直位移随填土高度增加变化图，最大沉降位移