论文部分内容阅读
随着我国水利水电事业的飞速发展,混凝土面板堆石坝的高度被不断地刷新。现已建成的水布垭面板堆石坝高度为233m,是世界上最高的面板堆石坝,300m级的面板堆石坝目前也在规划研究中,而且在建和拟建的混凝土面板堆石坝坝址多坐落在地形、地质条件复杂的地方。随着坝高的逐渐增加,筑坝规模的增大,坝体的受力和变形问题变得越来越突出,更好的分析高面板坝的受力和变形对坝体的设计和施工显得至关重要。高面板坝的数值分析是个相当复杂的非线性有限元问题,采用有限元法求解时,由于高面板坝的结构规模增大,模型网格数随之增多,为了更好地进行数值计算,如果采取加密措施,网格数量会更大,给计算带来了很大的困难。以前通用的直接法无论在内存占用还是求解效率上已经难以满足计算的需求,为了节约内存、提高计算速度和精度,需要寻求新的计算方法,这对高面板坝的数值模拟不仅具有重大的理论意义,还具有深远的工程实用意义。为了解决这一问题,本文主要做了如下三个方面的工作:1.研究了预处理矩阵的方法、SSOR-PCG法及其改进的迭代格式、总体刚度矩阵的存储技术,然后结合土石坝有限元静力计算的需要,提出了采用行压缩稀疏格式存储总体刚度矩阵,SSOR-PCG法改进迭代格式求解大型线性方程组的快速求解法及一套快速求解法的编程实现方案,使用Visual Fortran语言编写了基于快速求解法的土石坝三维有限元静力计算程序。2.选择规则面板堆石坝作为算例,采用不同密度的网格剖分方法,分别得到两种不同网格数量的有限元模型,然后分别采用快速求解法和直接法对规则面板坝模型进行有限元分析,验证了快速求解法应用到三维有限元计算程序的正确性,同时比较两种求解方法的计算效率并说明了快速求解法应用到大规模有限元计算的可行性。3.以水布垭混凝土面板堆石坝为工程实例,建立了三维有限元计算模型,采用本文编写的程序对水布垭面板堆石坝进行了有限元计算,得到不同工况下的受力和变形结果,并对面板和坝体的应力和位移结果进行了分析。