本文選題：非連續(xù)變形分析 + 地下洞室群　； 參考：《山東大學》2014年博士論文

【摘要】：隨著大型水利水電工程以及深部資源開發(fā)工程的飛速發(fā)展,巖石非連續(xù)大變形問題的理論與計算成為學術(shù)界研究的焦點。塊體系統(tǒng)非連續(xù)變形分析(DDA)是近三十年來才發(fā)展起來的一種新的數(shù)值計算方法,它在滿足彈性理論基本方程的條件下能夠反映出巖體變形的不連續(xù)性,既具有離散元法可計算塊體大位移的特點,又具有有限元理論基礎(chǔ)的嚴密性,具有廣泛的應用前景,是目前巖土工程分析中的研究熱點之一。非連續(xù)變形分析方法由于發(fā)展的較晚,在工程應用方面有諸多的局限性。本文對非連續(xù)變形分析方法進行了改進,考慮了位移收斂準則和錨桿的塑性變形。然后建立了非連續(xù)變形分析方法與有限元方法的耦合程序,討論了并行算法。最后結(jié)合四川巴底水電站項目對該方法進行了深入研究,主要研究了等效力學參數(shù)的確定和動態(tài)施工過程等。本文的主要工作概括如下： 1.對非連續(xù)變形分析方法進行了改進,加入了位移收斂判斷準則。然后采用非連續(xù)變形分析方法對雙裂隙巖塊進行單軸和雙軸壓縮模擬試驗,研究了裂隙角度和側(cè)向應力大小對巖塊特性的影響,得到了裂隙巖塊在這兩種加載試驗中的破壞過程和應力—應變曲線,同時也得到了巖塊中裂隙的起裂應力和巖塊的峰值強度。在雙軸壓縮模擬試驗中繪制了裂隙角度為45。巖塊在不同側(cè)向壓力下的強度包絡線,并根據(jù)包絡線方程式確定參數(shù)的c、φ的取值,同時根據(jù)等效連續(xù)介質(zhì)分析方法分析了某地下洞室的穩(wěn)定性。隨后,在程序中加入了錨桿的塑性變形力學模型,并采用改進后的程序模擬了劈裂試驗中含裂隙試塊的錨固效果,得到錨桿的塑性力學特性可以更好的描述錨桿的受力,更貼近工程實際情況。 2.應用有限元參數(shù)化設(shè)計語言APDL,開發(fā)了非連續(xù)變形分析方法程序與有限元程序圖形用戶界面的接口程序ANSYS2DDARF,可以生成BLAC文件,使前處理更加直觀方便。同時,在線性位移模式下,推導了非連續(xù)變形分析程序和有限元程序耦合方法的總體平衡方程,并且詳細推導了相應的的彈性子矩陣、慣性子矩陣、荷載子矩陣和接觸子矩陣等。非連續(xù)變形分析方法與有限元方法耦合后,每個塊體的自由度由網(wǎng)格節(jié)點的數(shù)量所確定,更加準確的描述了塊體的變形。最后采用耦合分析方法分析了某地下洞室圍巖的位移影響因素,得到側(cè)壓力系數(shù)、圍巖破碎面粘聚力、內(nèi)摩擦角、抗拉強度對圍巖的位移影響較大,同時頂拱為支護的關(guān)鍵位置,為工程提供了參考。 3.考慮了非連續(xù)變形分析方法基于OpenMP的并行算法。首先提出了剛度矩陣并行組裝策略,對非連續(xù)變形分析的剛度矩陣進行了并行化處理。然后對非連續(xù)變形分析的方程組計算采用改進的SSOR-PCG方法。最后通過算例表明改進后的并行算法可以使程序具有良好的加速比,明顯的提高了模型的計算速度。 4.在非連續(xù)變形分析方法的實際工程應用中,以四川巴底水電站地下洞室群為研究背景,首先數(shù)值模擬中的力學參數(shù)考慮了巖體中節(jié)理的隨機性和巖體力學參數(shù)的不確定性,通過現(xiàn)場實測數(shù)據(jù)和數(shù)值模擬的方法進行等效確定。然后研究了動態(tài)施工方法,突出了對開挖順序、分步開挖、分步支護過程進行模擬計算分析,并給出相應的解決方法,使程序可以更加真實的反映實際施工過程。 5.將非連續(xù)變形分析方法與有限元方法的耦合問題,動態(tài)施工過程,基于OpenMP并行問題,錨桿的塑性變形力學模型和位移收斂準則使用C++語言在VC++平臺中編制到計算程序中,使程序同時具備較好的分析塊體變形的功能和工程應用的功能。
[Abstract]:With the rapid development of large water conservancy and hydropower projects and deep resource development projects, the theory and calculation of rock discontinuous large deformation have become the focus of academic research. The discontinuous deformation analysis (DDA) of block system is a new numerical method developed in the last thirty years. It is satisfied with the basic equations of elastic theory. Under the conditions, it can reflect the discontinuity of rock mass deformation, not only the discrete element method can calculate the large displacement of the block body, but also has the tightness of the finite element theory basis. It has a wide application prospect. It is one of the hot topics in the analysis of geotechnical engineering. The non continuous deformation analysis method is applied in engineering because of the late development. There are many limitations in the surface. In this paper, the method of discontinuous deformation analysis is improved, the displacement convergence criterion and the plastic deformation of the bolt are considered. Then the coupling program of the discontinuous deformation analysis method and the finite element method is established, and the parallel algorithm is discussed. Finally, the method is deeply studied in combination with the project of Sichuan Ba Di hydropower station. The determination of equivalent mechanical parameters and dynamic construction process are mainly studied. The main works in this paper are summarized as follows:
1. the discontinuous deformation analysis method was improved, and the criterion of displacement convergence was added. Then the discontinuous deformation analysis method was used to simulate the double fracture rock mass by single axis and biaxial compression simulation, and the influence of the fracture angle and the lateral stress size on the rock mass characteristics was studied, and the fractured rock block was obtained in the two loading tests. The fracture stress and the peak strength of the rock mass in the rock mass are also obtained by the failure process and the stress strain curve. In the biaxial compression simulation test, the strength envelope of the fracture angle of 45. rock under different lateral pressures is drawn, and the value of the parameters of C and phi is determined according to the equation of the envelope line, and the equivalent continuum is based on the equivalent continuum. The stability of a certain underground cavern is analyzed by the qualitative analysis method. Then, the mechanical model of the plastic deformation of the bolt is added to the program, and the anchorage effect of the fracture test block is simulated by the improved program. The plastic mechanical properties of the bolt can be obtained to describe the stress of the bolt better and close to the actual situation of the engineering.
2. using the finite element parameterized design language APDL, the interface program ANSYS2DDARF of the discontinuous deformation analysis program and the finite element program graphical user interface is developed. The BLAC file can be generated to make the preprocessing more intuitive and convenient. At the same time, the discontinuous deformation analysis program and the finite element program coupling side are derived under the linear displacement mode. The general equilibrium equation of the method and the corresponding elastic submatrix, inertia submatrix, load submatrix and contact submatrix are derived in detail. After coupling the discontinuous deformation analysis method with the finite element method, the degree of freedom of each block is determined by the number of grid nodes, and the deformation of the block is described more accurately. Finally, the coupling is used. The analysis method analyses the influence factors of the displacement of the surrounding rock of an underground cavern, and obtains the side pressure coefficient, the cohesive force of the broken wall of the surrounding rock, the internal friction angle and the tensile strength of the surrounding rock, and the top arch is the key position of the support, which provides the reference for the project.
3. the parallel algorithm based on OpenMP for discontinuous deformation analysis is considered. First, the parallel assembly strategy of stiffness matrix is proposed, and the stiffness matrix of discontinuous deformation analysis is parallelized. Then the improved SSOR-PCG method is used for the equation group calculation of discontinuous deformation analysis. Finally, the improved parallel calculation is shown by an example. The method can make the program have a good speedup, and obviously increase the computation speed of the model.
4. in the practical engineering application of the discontinuous deformation analysis method, taking the underground cavern group of Sichuan Ba Di hydropower station as the research background, the mechanical parameters in the numerical simulation take into account the randomness of joints in the rock mass and the uncertainty of the mechanical parameters of the rock mass, and then make the equivalent determination by the field measured data and the numerical simulation method. The dynamic construction method is studied, and the simulation calculation and analysis of the sequence of excavation, step excavation and step support are carried out, and the corresponding solutions are given so that the program can be more truly reflected in the actual construction process.
5. the coupling problem of the discontinuous deformation analysis method and the finite element method, the dynamic construction process, based on the OpenMP parallel problem, the plastic deformation mechanics model and the displacement convergence criterion of the anchor, are compiled into the calculation program in the VC++ platform by the C++ language, so that the program has a good function of analyzing the deformation of the block and the work of the engineering application. Yes.
【學位授予單位】：山東大學
【學位級別】：博士
【學位授予年份】：2014
【分類號】：TV223

【參考文獻】

相關(guān)期刊論文 前10條

1 王建宙,于亞倫,唐春海;不連續(xù)變形分析及其在拆除爆破研究中的應用[J];爆破;2001年02期

2 林紹忠,蘇海東;大體積混凝土結(jié)構(gòu)仿真應力分析快速算法及應用[J];長江科學院院報;2003年06期

3 林紹忠;祁勇峰;蘇海東;;基于矩陣特殊運算的高階流形單元分析[J];長江科學院院報;2006年03期

4 郭培璽;林紹忠;;粗粒料顆粒隨機分布的數(shù)值模擬[J];長江科學院院報;2007年04期

5 陳健,盛謙,高鋒;構(gòu)皮灘工程高邊坡開挖數(shù)值模擬分析及穩(wěn)定性評價[J];長江科學院院報;1996年S1期

6 李靜;陳健云;周晶;;改進的SSOR-PCG迭代法在接觸問題研究中應用[J];大連理工大學學報;2006年04期

7 殷坤龍,姜清輝,汪洋;滑坡運動過程仿真分析[J];地球科學;2002年05期

8 黃潤秋;邊坡治理工程的數(shù)值模擬研究[J];地質(zhì)災害與環(huán)境保護;1996年01期

9 唐樹名,呂常新,鄧安福;預應力錨索加固路塹邊坡的數(shù)值流形分析[J];公路交通科技;2002年06期

10 朱維申,陳衛(wèi)忠,申晉;雁形裂紋擴展的模型試驗及斷裂力學機制研究[J];固體力學學報;1998年04期

相關(guān)博士學位論文 前1條

1 江巍;基于互補理論的非連續(xù)變形分析方法[D];中國科學院研究生院（武漢巖土力學研究所）;2009年

，

本文編號：1989286