岩石力学实验报告心得.docxVIP

岩石力学实验报告心得 　　岩石力学数值试验实验报告　　姓名：莫道兴学号：　　班级：采矿101班指导教师：左宇军老师　　同组人：许云飞辜巍　　实验名称：不同均质度对岩石力学性质影响的数据拉伸试验　　XX年5月16号　　一、实验目的　　1、通过对RFPA学习，知道RFPA基本使用方法；　　2、了解RFPA模拟试验的条件和RFPA的基本功能。　　二、实验原理　　RFPA是一种基于有限元应力分析和统计损伤理论的材料破裂过程分析数值计算方法，是一个能够模拟材料渐进破裂直至失稳全过程的数值试验工具。　　三、实验步骤如下;　　RFPA数值模型　　本模型拟采用尺寸为100X50mm的岩石试件模型。试件模型划分为100X50个单元。采用平面应力模型。逐渐个加载过程采用位移控制的直接拉伸加载方式，即通过在试件端部施加位移实现拉伸加载，每步加载位移量S=-。　　操作步骤　　第一步，启动RFPA，新建模型建立存放的根目录　　第二步，划分网格，单击　　单击确定在弹出的窗口中设置模型的大小，2D2D2D2D2D　　第三步，得到材料图形　　第四步，施加荷载　　第五步，单击求解控制信息，设置条件如下：　　岩石力学与工程　　实验报告　　一、实验目的　　1、熟悉运用岩石力学的phase软件；　　2、运用岩石力学的基本理论，来计算某地的地应力值。　　二、实验软件　　1、岩石力学phase软件；　　2、autoCADXX;　　3、matlab软件；　　4、microsoftofficeXX软件。　　三、实验方法与步骤　　1、选取九龙河溪古水电站地质构造带作为实验基础，并用运用autoCAD软件绘制将该地区的断层、节理等地质构造单元；　　2、在phase软件中导入已绘制各种边界；　　3、进行网格划分；　　4、定义材料，并将所计算的模型设置正确的材料颜色；　　5、运用matlab软件进行数据处理和计算；　　、已知理塘、雅江、呷巴、长河坝、乾宁的最大主应力及最小主应力，利用工程力学的力学计算方法，将已知应力点的σ1、σ3、最大主应力方向转换成σx、σy、τxy、τyx.可得出如表1所示的的实验数据：　　地名　　理塘　　雅江　　呷巴　　长河坝　　乾宁σx　　σy　　τxy　　x坐标---y坐标　　　　表格1：将σ1、σ3转化为σx、σy的数据表　　、运用matlab软件编程，求出各个地区的ν、λ、α值　　令E=E;v=ν；l=λ；a=α；　　Yanshi1的源程序：　　E=input(请输入E的值:);　　v=input(请输入v的值:);　　G=E/[2*(1+v)]　　l=E*v/[(1+v)*(1-2*v)]　　a=l+2*G　　对于⑤古生代到三叠纪的变质分布有：E=12500MPa，　　运行matlab程序：yanshi1　　请输入E的值:12500　　请输入v的值:　　G=+003　　l=+003　　a=+004　　即求得理塘G=；l=；a=+003　　、在利用autoCAD的测量距离方法，得出理塘、雅江、呷巴、长河坝、乾宁的坐标，求得的数据如表2：　　地名E(MPa)μλGαx坐标理塘12500　　-雅江12500　　-65514呷巴12500　　长河坝12500　　乾宁12500　　-表格2：各个地区的x,y坐标　　、建立matlab的矩阵模型，求出系数A1,A2,A3,A4,A5,B1,B2,B3,B4,B5　　Matlab的矩阵模型如下：　　A=[α02*α*X0a*Y0λ02*λ*Yλ*X　　λ02*λ*X0λ*Y0α02*α*Yα*X　　0102*YX102*X0Y];　　b=[σx;σy;τxy/G];y坐标　　A*x=b；　　即可得如下的系数矩阵：　　A=[14271,0,-,0,,0,,0,,-271,0,-,0,,0,,0,,-271,0,,0,,0,,0,,71,0,,0,,0,,0,,271,0,-,0,,0,,0,,-,0,-,0,,0,14271,0,,-,0,-,0,,0,14271,0,,-,0,,0,,0,14271,0,,,0,,0,,0,14271,0,,,0,-,0,,0,14271,0,,-,1,0,,-,1,0,-,0,,1,0,,-65514,1,0,-,0,,1,0,,,1,0,,0,0,1,0,,,1,0,,0,,1,0,,-,1,0,-,0,];b=[-;-;-;-;-;-;-;-;-;-;;;;;];　　、利用以上模型来求解，从中任意选取10组可求A1,A2,A3,A4,A5和B1,B2,B3,B4,B5的值分别如下：　　A1=-,A2=0,A3=-1E-10,A