高面板堆石坝非线性数值分析快速求解法的研究
发布时间:2018-03-22 10:36
本文选题:高面板坝 切入点:水布垭 出处:《大连理工大学》2014年硕士论文 论文类型:学位论文
【摘要】:随着我国水利水电事业的飞速发展,混凝土面板堆石坝的高度被不断地刷新。现已建成的水布垭面板堆石坝高度为233m,是世界上最高的面板堆石坝,300m级的面板堆石坝目前也在规划研究中,而且在建和拟建的混凝土面板堆石坝坝址多坐落在地形、地质条件复杂的地方。随着坝高的逐渐增加,筑坝规模的增大,坝体的受力和变形问题变得越来越突出,更好的分析高面板坝的受力和变形对坝体的设计和施工显得至关重要。 高面板坝的数值分析是个相当复杂的非线性有限元问题,采用有限元法求解时,由于高面板坝的结构规模增大,模型网格数随之增多,为了更好地进行数值计算,如果采取加密措施,网格数量会更大,给计算带来了很大的困难。以前通用的直接法无论在内存占用还是求解效率上已经难以满足计算的需求,为了节约内存、提高计算速度和精度,需要寻求新的计算方法,这对高面板坝的数值模拟不仅具有重大的理论意义,还具有深远的工程实用意义。为了解决这一问题,本文主要做了如下三个方面的工作: 1.研究了预处理矩阵的方法、SSOR-PCG法及其改进的迭代格式、总体刚度矩阵的存储技术,然后结合土石坝有限元静力计算的需要,提出了采用行压缩稀疏格式存储总体刚度矩阵,SSOR-PCG法改进迭代格式求解大型线性方程组的快速求解法及一套快速求解法的编程实现方案,使用Visual Fortran语言编写了基于快速求解法的土石坝三维有限元静力计算程序。 2.选择规则面板堆石坝作为算例,采用不同密度的网格剖分方法,分别得到两种不同网格数量的有限元模型,然后分别采用快速求解法和直接法对规则面板坝模型进行有限元分析,验证了快速求解法应用到三维有限元计算程序的正确性,同时比较两种求解方法的计算效率并说明了快速求解法应用到大规模有限元计算的可行性。 3.以水布垭混凝土面板堆石坝为工程实例,建立了三维有限元计算模型,采用本文编写的程序对水布垭面板堆石坝进行了有限元计算,得到不同工况下的受力和变形结果,并对面板和坝体的应力和位移结果进行了分析。
[Abstract]:With the rapid development of water conservancy and hydropower industry in China, The height of concrete face rockfill dam has been constantly refreshed. The height of Shuibuya concrete face rockfill dam is 233m, which is the highest face rockfill dam in the world. Moreover, the sites of concrete face rockfill dams under construction and proposed construction are mostly located in places where the topography and geological conditions are complex. With the increasing of dam height and dam construction scale, the stress and deformation of dam body become more and more prominent. It is very important to analyze the stress and deformation of high-faced dam for the design and construction of dam body. The numerical analysis of high face slab dam is a very complicated nonlinear finite element problem. When the finite element method is used to solve the problem, the grid number of the model increases because of the increase of the structure scale of the high face slab dam, in order to carry out the numerical calculation better. If encryption measures are taken, the number of grids will be larger, which will bring great difficulties to the computation. In order to save memory, it is difficult to satisfy the needs of computation by using the common direct method in memory occupation and efficiency. To improve the calculation speed and precision, we need to find a new calculation method, which is not only of great theoretical significance, but also of far-reaching engineering practical significance for the numerical simulation of high-faced dam. The main work of this paper is as follows:. 1. SSOR-PCG method and its improved iterative format, the storage technology of total stiffness matrix, and the need of finite element static calculation of earth-rock dam are studied. This paper presents an improved iterative method for solving large linear equations by using row compression sparse format to store the total stiffness matrix / SSOR-PCG method. The 3D finite element static calculation program of earth-rock dam based on fast solution method is compiled with Visual Fortran language. 2. Taking regular face rockfill dam as an example, two finite element models with different mesh numbers are obtained by using different density mesh generation method. Then the fast solution method and the direct method are used to analyze the regular face slab dam model respectively, which verifies the correctness of the fast solution method applied to the 3D finite element calculation program. At the same time, the computational efficiency of the two methods is compared and the feasibility of applying the fast solution method to large-scale finite element calculation is illustrated. 3. Taking Shuibuya concrete face rockfill dam as an engineering example, a three-dimensional finite element calculation model is established. The finite element calculation of Shuibuya concrete face rockfill dam is carried out by using the program compiled in this paper, and the results of stress and deformation under different working conditions are obtained. The stress and displacement of face slab and dam are analyzed.
【学位授予单位】:大连理工大学
【学位级别】:硕士
【学位授予年份】:2014
【分类号】:TV641.43
【参考文献】
相关期刊论文 前10条
1 王忠雷;赵国群;马新武;;三维有限元刚度矩阵的压缩存储算法[J];材料科学与工艺;2012年02期
2 包劲青;杨强;刘耀儒;;适宜PCG方法的水工结构有限元刚度矩阵形成、组集和存储的新方法[J];华北水利水电学院学报;2009年04期
3 韩林;张子明;倪志强;;应用改进算法的对称逐步超松弛预处理共轭梯度法进行大体积混凝土仿真计算[J];河海大学学报(自然科学版);2010年03期
4 林绍忠;用预处理共轭梯度法求解有限元方程组及程序设计[J];河海大学学报;1998年03期
5 张永杰;孙秦;;有限元线性方程组改进的全稀疏存贮方法[J];航空计算技术;2006年02期
6 张永杰;孙秦;李江海;;大型稀疏线性方程组的改进ICCG方法[J];计算物理;2007年05期
7 陈璞,孙树立,袁明武;有限元分析快速解法[J];力学学报;2002年02期
8 朱俊高,殷宗泽;高土石坝混凝土防渗墙弹塑性应力变形分析[J];水利学报;1997年07期
9 林绍忠;对称逐步超松弛预处理共轭梯度法的改进迭代格式[J];数值计算与计算机应用;1997年04期
10 李守德,俞洪良;Goodman接触面单元的修正与探讨[J];岩石力学与工程学报;2004年15期
,本文编号:1648280
本文链接:https://www.wllwen.com/kejilunwen/shuiwenshuili/1648280.html