岩石力学-实验报告

1、岩石力学 -实验报告岩石力学与工程实验报告一、实验目的1、熟悉运用岩石力学的 phase 软件；2、运用岩石力学的基本理论，来计算某地的地应力值。二、实验软件1、岩石力学 phase 软件；2、autocad2006;3、 matlab6.5 软件；4、microsoftoffice2003 软件 。三、实验方法与步骤 1、选取九龙河溪古水电站地质构造带作为实验基础，并用 运用 autocad 软件绘制将该地区的断层、节理等地质构造单元；2、在 phase 软件中导入已绘制各种边界（断裂边界、材料 边界、 boundry）；3、进行网格划分；4、定义材料，并将所计算的模型设置正确的材料颜色；5

2、、运用 matlab 软件进行数据处理和计算；5.1、已知理塘、雅江、呷巴、长河坝、乾宁的最大主应力 及最小主应力， 利用工程力学的力学计算方法， 将已知应力点的(T3、最大主应力方向转换成c x、cy、Tcy、wx .可得出如表 1 所示的的实验数据：地名理塘雅江呷巴长河坝乾宁 c x7.4025735.3528234.5533733.1198512.883026cy5.897427315.9671774085.1466269146.090149323.22697392Tyl.960520.760290.044860.425860.56961x 坐标 -16.2352-8.73521.739

3、37.3222-0.3815y 坐标 14.60414.60414.001413.072820.9622表格1:将c1、c3转化为cx、cy 的数据表5.2、运用 matlab软件编程，求出各个地区的v、入、a值令 e=e;v=v; 1=入；a=a; yanshil 的源程序：e=inpu（t '请输入 e 的值: '）;v=inpu（t '请输入 v 的值: ' ）;g=e/2* （ 1 +v） l=e*v/ （1+v）*（1-2*v）a=l+2*g 对于古生代到三叠 纪的变质分布有： e=12500mpa， 0.22 运行 matlab 程序： yanshi

4、1 请 输 入 e 的 值 ：12500 请 输 入 v 的 值 ：0.22g=5.1230e+003l=4.0252e+003a=1.4271e+004 即 求 得 理 塘 g=370.3704； l=864.1975； a=1.6049e+0035.3、 在利用autocad的测量距离方法，得出理塘、雅江、呷巴、长河坝、乾宁的坐标，求得的数据如表2：地名 e ( mpa ) 入g a x 坐标理塘125000.224025.17565122.9508214271.08-121764 雅 江 125000.224025.17565122.9508214271.08-65514 呷 巴 125

5、000.224025.17565122.9508214271.0813044.75 长 河 坝 125000.224025.17565122.9508214271.0854916.5 乾 宁 125000.224025.17565122.9508214271.08-2861.25表格 2：各个地区的 x，y 坐标5.4、建立 matlab 的矩阵模型，求出系数 a1， a2，a3， a4， a5， b1， b2， b3， b4， b5matlab 的矩阵模型如下：a= a 02*a*x0a*y0入02*入*y入*x入02*入*x0入*y0a 02*a*ya*x0102*yx102*x0y;b=

6、 ox; oy; Ty/g;y 坐标 10953010953010501198046157217a*x=b ；即可得如下的系数矩阵：a= 1 427 1 ，0，-3475388088 ，0，1 563 1 1 1 392 ，0，4025.2，0， 881760312，-490124453 14271，0，-1869900588，0，1563111392，0，4025.2，0，881760312，-26370695314271，0，372324682，0，1498604846，0，4025.2，0，845376529.2，5250792914271，0，1567426743，0，13992144

7、66，0，4025.2，0，789309518.4，22104989614271 ， 0 ， -81667225 ， 0 ， 2243636672 ， 0， 4025.2 ， 0 ，1265655712，-115173054025.2， 0， -980248906 ， 0 ， 440880156 ， 0 ， 14271 ， 0， 3126205260，-17376940444025.2， 0， -527413906 ， 0 ， 440880156 ， 0 ， 14271 ， 0，3126205260，-9349502944025.2，0，105015858，0，422688265，0，1427

8、1，0，2997209691，1861623414025.2，0，442099792，0，394654759，0，14271，0，2798625024，7837133724025.2，0，-23034610，0，632827856，0，14271，0，4487273343，-408336120，1，0，219060，-121764，1，0，-243528，0，10953040，1， 0，219060，-65514，1， 0，-131028，0，109530 0，1，0，210021，13044.8，1，0，26089.6， 0，105010.5 0，1， 0，196092，54916.5，1，

9、0，109833，0， 98046 0，1，0，314433，-2861.3，1，0，-5722.6， 0，157216.5;b=-5.89743;-5.96718;-5.14663;-6.09015;-3.22697;-7.4026;-5.3 528;-4.5534;-3.1199;-2.883;0.000382689;0.000148407;0.000008756;0.000083127; 0.000111185;5.5 、利用以上模型来求解，从中任意选取 10 组可 求 a1，a2，a3，a4， a5 和 b1，b2，b3，b4，b5 的值分别如下：a1=-0.0007 ， a2=0 ，

10、 a3=-1e-10 ， a4=-1e-09 ， a5=3e-09 ， b1=0.00013，b2=-0.0003，b3=-2e-09，b4=6.5e-10，b5=1.7e-095.6、 根据以上的系数 a1，a2，a3， a4， a5， b1， b2，b3，b4，b5 可 将研究区域的不同坐标值找出，利用以下式子求出cx,ay,T<y值：a*a1+2 o(*a3+ ay*a5+ ?*b2+2 ?y*b4+ ?x*b5= ox?*a1+2 ?x*a3+ 入 y*a5+ a*b2+2 ay*b4+ ax*b5= cyb1+2x*b3+y*b5+a2+2y*a4+x*a5= t xy/g 求

11、得的实验数据见表 3：xcyycxtxy-270030.1-30000026.1019082723.443972-3.3845051-248837-30 000026.0187855322.919628-3.1933399-236123.4-30000025.96892 11422.60508-3.0786621-196566.8-30000025.8137738721.6264-2.7 218553-184185.2-30000025.7652113321.320064-2.6101715-15396 5.8-30000025.6466860720.572399-2.3375878-1237

12、46.4-30000025. 5281608119.824733-2.0650041-93526.95-30000025.4096355519.0 77067-1.7924204-60145.15-30000025.278706918.25116-1.491311 5-26763.34-30000025.1477782617.425253-1.19020256618.4692-3 0000025.0168496116.599346-0.889093645416.512-30000024.864 6776515.639435-0.539129473948.703-30000024.7527699

13、514.933 514-0.2817648104011.98-30000024.634857114.189711-0.0105895 134075.26-30000024.5169442513.4459090.26058577164138.53-3 0000024.399031412.7021060.53176105175764.95-30000024.3534 307912.4144540.63663307217176.63-30000024.1910077211.3898 781.01017268258588.32-30000024.0285846510.3653021.3837122 9

14、-300000267798.0381-0.9138556976.8671631-3.0855215-3000002 35596.07620.6249763837.8493492-3.1178096-300000203394.114 32.1638084638.8315352-3.1500977-300000158059.48854.330209 90210.214278-3.1955536-300000124366.99785.94026973911.241 926-3.2293362-30000090674.507127.55032957612.269575-3.2631 188-30000

15、056982.016449.16038941313.297223-3.2969014-30000 042538.324659.85060873513.737767-3.3113838-30000010072.42 38211.4020536414.728004-3.3439365-300000-22393.47712.95349 85415.71824-3.3764892-300000-49054.759214.2275586916.53143 1-3.4032218-300000-75716.041315.5016188417.344622-3.429954 4-300000-113096.

16、70117.2879248618.484762-3.4674351-300000-1 50477.36119.0742308819.624903-3.5049157-300000-187858.0212 0.860536920.765043-3.5423964-300000-225238.6822.6468429221.905184-3.579877-300000-262619.3424.4331489523.045324-3.617 3577-300000-30000026.2194549724.185465-3.6548383表格3:不同坐标的应力值 5.7、在phase中设定边界应力 值

17、导入所求的模型，即可得到所需的实验模型。 模型如图1所示:图1 :实验模型图四、实验成果由以上的模型，在phase软件中经过计算，可建立如下图所 示的成果：图 2: maximumshearstrain 图图 3: strengthfactor 图五、实验中遇到的问题及心得体会本题是岩石力学的基本实验之一，旨在通过学习了一些岩石力学的基本知识，来实地计算某地的地应力。 本题以九龙河溪古 水电站地质构造为实验数据的基础，来对该地区的地应力值进行了定量的计算。实验的过程中，由于是第一次接触岩石力学 phase 软件，也遇到了一定的问题。具体来说可以概括如下：1、工程力学基础知识学得不是特别的扎实，

18、致使在已知了理塘、雅江、呷巴、长河坝、乾宁的最大主应力及最小主应力后，处理实脸数16花了一定的时间，利用公式x y2x y2) 2x2花费了一定的时间，此题可用应力圆的方法来计算，使问题得到一定的简化2、由于知识的局限性， 所学的matlab知识也不是特别的熟 练。此题本来也可以 excel 来解决问题，求解线性方程组也有一 定的方便， 但关于这方面的知识自己学得不是特别的精通， 致使 用了 matlab 软件来解决。本题在一定的程度上也反映 matlab 的 局限性，如在求解的过程中，可从中选取任意 15 个线性方程组 中的 10 组数据来解决需要求的系数 a1，a2，a3，a4， a5，b

19、1， b2, b3, b4, b5。但matlab不能够显示数据过于小的数据，因 此在运行的过程中，会遇到显示 -inf 的情况，此方式就表示了在 matlab 的情况下，不能正常的显示数据。如果采用 excel 来解决 问题，或许不会遇到这样类似的问题。3、实验过程中，对计算模型的思路不清晰。致使在做题的过程中，在得出一些实验数据之后， 不能进行下一步的正确操作。 只有在等到老师讲解了一些相关的步骤之后， 才会有一定的思路。 同是， 对模型也不是特别的清楚， 由于 phase 软件全部是英文的 命令，所以在操作的过程中，也有了一定的难度。在实验时，我 不仅遇到了一些难题，而且还转了一定的圈子。比如在添加在 phase 软件中设定边