毕业设计_堆石坝应力分析数值模拟GEOstudio.docx

坝体应力变形分析模型建立和参数取值应力变形模型是在渗透模型的基础上做了部分修改。地基本为无限边界，为了避免地基范围过小，而影响计算结果的准确性，将地基沿深度方向延伸了45m，上下游方向延伸45m。对部分尖锐三角形网格，通过绘制自由点调整集合形状等方式进行修正。为模拟分层填筑过程将，坝体分为9层，层高约为3m。网格边长为3m，结点1452个，单元1391个。土的刚度E一般是围压或者上覆土层应力的函数，邓肯—张模型建议用式（5-1）描述它们之间的关系。 （5-1）式中为围压的初始切线模量；为加载模量数；为大气压（作为标准化的参数使用）；为围压；n为定义围压对初始模量的影响指数。切线模量Et是偏应力（）和围压的函数，关系式为（5-2）（5-2）Et为切线模量；为土的内摩擦角；c为土的粘聚力；为双曲线的渐近线和最大剪切应力的比值，通常在0.75和1之间；为最大主应力；为最小主应力。各种材料的部分参数取值为表5-3、表5-4、表5-5：表5-3 材料部分参数取值参数坝壳过渡区沥青混凝土灰岩岩基泥岩岩基容重19kN/m321.5 kN/m324 kN/m323 kN/m323 kN/m3模型弹塑性弹塑性双曲线线弹性线弹性泊松比0.310.350.460.240.305弹性模量Ei方程80MPaEi方程23224MPa2205MPa表5-4 坝壳料Ei方程Ei（kPa）26084369584842463979755418504593233100493σ1（kPa）2050100200300400500600表5-5沥青混凝土Ei方程Ei（kPa）2050100200300400500600700σ1（kPa）51402664658073698093109943119218126955133651139591边界条件设定，岩基两侧设定位移控制，水平向位移为0；岩基底部是位移控制边界，竖向位移为0；上下游设置静水压力边界条件，上下游坝壳材料透水性相当强，在坝坡面的水压力是内外平衡的，而岩基的透水性相对于坝壳相当小，考虑这种实际受荷情况，将静水压力荷载施加在基岩、心墙上。而不是坝坡面上，这样的模型能够更为贴近真实情况。有限元计算成果分析在施工期，填筑完毕时的水平位移等值线分布图（图5-3），上下游变形较为图5-3 施工结束坝体水平位移等值线分布图对称，上游变形方向沿横坐标负方向，下游位移方向沿着横坐标正方向。在心墙竖直面上，大约在1/3处位移量最大。与工程经验较为吻合。最大位水平向移量约为0.025m。施工期，坝体竖直向位移等值线分布（图5-4），垂向位移关于心墙有对称性。在中心位置的垂向位移量相对上下部位和两侧坝体较大，地基一定范围内有微小的变形，离坝基较远处竖向位移量为0。这从侧面反映在模型建立时，坝基的范围足够大，没有因为区域较小而影响有限元计算的准确性。坝体上部分，心墙竖向位移较同等高程的坝壳更大，这是由于其刚度较小的缘故，计算结果较为合理，位移等值线较为光滑，不均均匀沉降问题不是很突出。坝体中下部分心墙和坝壳竖向位移很接近，沉降较为均匀。竖向位移最大为0.1039m，即为施工期的沉降量，约为坝体高程的0.3%。图5-4 施工结束坝体竖向位移等值线分布图 图5-5 施工结束坝体竖向总应力等值线分布图施工结束竖向总应力等值线分布（图5-5），沿竖直方向从上到下，竖向总应力不断递增，施工期的竖向总应力主要是自重应力引起，忽略了施工荷载，这样的整体应力分布情况与实际相一致。竖向总应力等值线在地基区域较为光滑平顺，在心墙、过渡料、坝壳料接触处有一定的波动，但是总体来说应力过渡还是较为平顺的。拱效应并不明显，过渡料有效的起到了协调心墙和坝壳材料变形的作用，一定程度上减小了应力差异可能导致的裂缝产生的可能。施工结束后的第三主应力等值线分布（图5-6）在-271kPa到616kPa范围变动，其中出现拉应力的部位主要有以下一些部位，坝趾和坝踵浅层以及坝基小范围内有拉应力的存在。可以考虑在上下游增设压重。沥青混凝土心墙没有拉应力的出现。心墙仅有50cm的厚度，心墙抗裂的关键问题就是避免拉应力的出现。在施工期计算中心墙没有受拉的应力状态。图5-6 施工结束坝体最小主应力等值线分布图正常蓄水位工况的应力变形分析，是在施工期模型的基础上，添加上浸润线，施加水荷载。浸润线的添加，计算式可以区分浮容重和天然容重，水荷载施加在较坝壳料更为不透水的岩基和心墙上，而不是坝面上。坝壳的透水性较好，坝面上的水压力处于平衡状态。蓄水到正常水位后坝体水平位移等值线分布图（5-7），由于蓄水的影响，上、下游坝壳的变位对称性被打破，同等高程上游的水平位移数值上较下游小。以最大位移为例-0.025m，下游为0.02m。上游蓄水后，水的浸润作用，土体的容重由天然容重转化为浮容重，而下游浸