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

【摘要】：隨著大型水利水電工程以及深部資源開(kāi)發(fā)工程的飛速發(fā)展,巖石非連續(xù)大變形問(wèn)題的理論與計(jì)算成為學(xué)術(shù)界研究的焦點(diǎn)。塊體系統(tǒng)非連續(xù)變形分析(DDA)是近三十年來(lái)才發(fā)展起來(lái)的一種新的數(shù)值計(jì)算方法,它在滿足彈性理論基本方程的條件下能夠反映出巖體變形的不連續(xù)性,既具有離散元法可計(jì)算塊體大位移的特點(diǎn),又具有有限元理論基礎(chǔ)的嚴(yán)密性,具有廣泛的應(yīng)用前景,是目前巖土工程分析中的研究熱點(diǎn)之一。非連續(xù)變形分析方法由于發(fā)展的較晚,在工程應(yīng)用方面有諸多的局限性。本文對(duì)非連續(xù)變形分析方法進(jìn)行了改進(jìn),考慮了位移收斂準(zhǔn)則和錨桿的塑性變形。然后建立了非連續(xù)變形分析方法與有限元方法的耦合程序,討論了并行算法。最后結(jié)合四川巴底水電站項(xiàng)目對(duì)該方法進(jìn)行了深入研究,主要研究了等效力學(xué)參數(shù)的確定和動(dòng)態(tài)施工過(guò)程等。本文的主要工作概括如下： 1.對(duì)非連續(xù)變形分析方法進(jìn)行了改進(jìn),加入了位移收斂判斷準(zhǔn)則。然后采用非連續(xù)變形分析方法對(duì)雙裂隙巖塊進(jìn)行單軸和雙軸壓縮模擬試驗(yàn),研究了裂隙角度和側(cè)向應(yīng)力大小對(duì)巖塊特性的影響,得到了裂隙巖塊在這兩種加載試驗(yàn)中的破壞過(guò)程和應(yīng)力—應(yīng)變曲線,同時(shí)也得到了巖塊中裂隙的起裂應(yīng)力和巖塊的峰值強(qiáng)度。在雙軸壓縮模擬試驗(yàn)中繪制了裂隙角度為45。巖塊在不同側(cè)向壓力下的強(qiáng)度包絡(luò)線,并根據(jù)包絡(luò)線方程式確定參數(shù)的c、φ的取值,同時(shí)根據(jù)等效連續(xù)介質(zhì)分析方法分析了某地下洞室的穩(wěn)定性。隨后,在程序中加入了錨桿的塑性變形力學(xué)模型,并采用改進(jìn)后的程序模擬了劈裂試驗(yàn)中含裂隙試塊的錨固效果,得到錨桿的塑性力學(xué)特性可以更好的描述錨桿的受力,更貼近工程實(shí)際情況。 2.應(yīng)用有限元參數(shù)化設(shè)計(jì)語(yǔ)言APDL,開(kāi)發(fā)了非連續(xù)變形分析方法程序與有限元程序圖形用戶界面的接口程序ANSYS2DDARF,可以生成BLAC文件,使前處理更加直觀方便。同時(shí),在線性位移模式下,推導(dǎo)了非連續(xù)變形分析程序和有限元程序耦合方法的總體平衡方程,并且詳細(xì)推導(dǎo)了相應(yīng)的的彈性子矩陣、慣性子矩陣、荷載子矩陣和接觸子矩陣等。非連續(xù)變形分析方法與有限元方法耦合后,每個(gè)塊體的自由度由網(wǎng)格節(jié)點(diǎn)的數(shù)量所確定,更加準(zhǔn)確的描述了塊體的變形。最后采用耦合分析方法分析了某地下洞室圍巖的位移影響因素,得到側(cè)壓力系數(shù)、圍巖破碎面粘聚力、內(nèi)摩擦角、抗拉強(qiáng)度對(duì)圍巖的位移影響較大,同時(shí)頂拱為支護(hù)的關(guān)鍵位置,為工程提供了參考。 3.考慮了非連續(xù)變形分析方法基于OpenMP的并行算法。首先提出了剛度矩陣并行組裝策略,對(duì)非連續(xù)變形分析的剛度矩陣進(jìn)行了并行化處理。然后對(duì)非連續(xù)變形分析的方程組計(jì)算采用改進(jìn)的SSOR-PCG方法。最后通過(guò)算例表明改進(jìn)后的并行算法可以使程序具有良好的加速比,明顯的提高了模型的計(jì)算速度。 4.在非連續(xù)變形分析方法的實(shí)際工程應(yīng)用中,以四川巴底水電站地下洞室群為研究背景,首先數(shù)值模擬中的力學(xué)參數(shù)考慮了巖體中節(jié)理的隨機(jī)性和巖體力學(xué)參數(shù)的不確定性,通過(guò)現(xiàn)場(chǎng)實(shí)測(cè)數(shù)據(jù)和數(shù)值模擬的方法進(jìn)行等效確定。然后研究了動(dòng)態(tài)施工方法,突出了對(duì)開(kāi)挖順序、分步開(kāi)挖、分步支護(hù)過(guò)程進(jìn)行模擬計(jì)算分析,并給出相應(yīng)的解決方法,使程序可以更加真實(shí)的反映實(shí)際施工過(guò)程。 5.將非連續(xù)變形分析方法與有限元方法的耦合問(wèn)題,動(dòng)態(tài)施工過(guò)程,基于OpenMP并行問(wèn)題,錨桿的塑性變形力學(xué)模型和位移收斂準(zhǔn)則使用C++語(yǔ)言在VC++平臺(tái)中編制到計(jì)算程序中,使程序同時(shí)具備較好的分析塊體變形的功能和工程應(yīng)用的功能。
[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.
【學(xué)位授予單位】：山東大學(xué)
【學(xué)位級(jí)別】：博士
【學(xué)位授予年份】：2014
【分類號(hào)】：TV223

【參考文獻(xiàn)】

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

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

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

3 林紹忠;祁勇峰;蘇海東;;基于矩陣特殊運(yùn)算的高階流形單元分析[J];長(zhǎng)江科學(xué)院院報(bào);2006年03期

4 郭培璽;林紹忠;;粗粒料顆粒隨機(jī)分布的數(shù)值模擬[J];長(zhǎng)江科學(xué)院院報(bào);2007年04期

5 陳健,盛謙,高鋒;構(gòu)皮灘工程高邊坡開(kāi)挖數(shù)值模擬分析及穩(wěn)定性評(píng)價(jià)[J];長(zhǎng)江科學(xué)院院報(bào);1996年S1期

6 李靜;陳健云;周晶;;改進(jìn)的SSOR-PCG迭代法在接觸問(wèn)題研究中應(yīng)用[J];大連理工大學(xué)學(xué)報(bào);2006年04期

7 殷坤龍,姜清輝,汪洋;滑坡運(yùn)動(dòng)過(guò)程仿真分析[J];地球科學(xué);2002年05期

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

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

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

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

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

，

本文編號(hào)：1989286