岩石力学数值试验实验报告

姓名：郑周立学号： 1108010103 班级：采矿111班指导教师：左宇军

同组人：郑周立、周义现、胡斌、朱红伟、高言、

王坤

实验名称：圆孔对岩石力学性质影响的数值加载

试验

2014年5月16日

圆孔对岩石力学性质影响的数值加载试验

一、实验目的:

1.通过对RFPA2D学习，知道RFPA2D基本使用方法。

2.了解RFPA2D模拟试验的条件和RFPA2D的基本功能。

3.通过操作端部效应对岩石力学性质影响的数值实验，了解每一步操作以及岩石破裂过程，最终完成实验得到结果。

二、实验原理:

RFPA-2D是一种基于有限元应力分析和统计损伤理论的材料破裂过程分析数值计算方法，是一个能够模拟材料渐进破裂直至失稳全过程的数值试验工具。

三、

1、试样尺寸： 100mm*51mm

2、基元数: 100*51

3、应力分析模式: 平面应变

4、圆孔：半径10mm

5、加载方式：单轴压缩

6、加载条件：竖向位移加载

7、均质度m=2

8、加载量：每步0.002mm

9、实验内容：

（1）、应力-应变曲线；

（2）、强度；

（3）、破坏模式

四、实验内容：

(一)、操作步骤：

第一步启动RFPA，新建模型建立存放的根目录

第二步划分网格，单击在弹出的窗口中设置模型的大小，单击确定第三步选择施加荷载模式...

（二）实验结果

弹性模量图

第1步

第4步（开始破坏）

第7步（开始横向破坏）

第32步（彻底破坏）

第200步

最大剪应力图第1步

第4步（开始破坏）

第33步（彻底破坏）

第200步

最大主应力图

第1步

第3步（开始破坏）

第32步（彻底破坏）

第200步

最小主应力图第1 步

第4 步（开始破坏）

第31 步（完全破坏）

第200步

应力/力-加载步曲线

位移/应变-加载步曲线

在本试验中破坏的模式有脆性破坏和剪切破坏。

五、实验心得

这次试验是端部效应对岩石性质影响的数值实验，在实验中尤其注重对垫板的性质不同而产生的滑动条件下，不同垫板的弹性模量对岩石力学的影响，垫板刚度和试样刚度是很重要的，不能忘记输入，在加载方式里，设置好加载条件以，加载量，缺一不可。每一步都需要认真对待才能很好的完成实验。

在这门课中以及实验中我收获很多，从书本学习到实际操作，从实际操作到获得结果。这个过程是不容易的，由于这个软件不是太熟悉，所以我不断操作，不断实验和请教同学，最后才得以完成，所以这个实验是不容易的。

在岩石力学教学中，岩石的拉、压、剪基本实验及岩石的破裂与失稳过程是一个重要的基本教学内容。由于岩石材料的非均匀性、非连续性，以及外载荷作用下微缺陷之间相互作用的复杂性，现有的解析方法尚缺少有效的手段对此过程进行研究，理论上很难对岩石的破裂与失稳过程做准确的描述。因此，目前有关岩石破裂与失稳过程的研究，仍然主要依赖于现场观测和实验物理实验。现场观测对工程而言是非常必要的，但由于这种方法受到现场条件、人力、物力和人力的限制，很难在教学中得到充分利用；物理实验虽直观，但有关岩石破裂过程现象的复杂性，以及实验室观测手段、经费等条件限制，通常的岩石力学教学很难通过大量的物理实验向学生直观演示各种岩石

变形、破坏的复杂现象。因此，数值试验方法可能补充常规的实验教

学，达到岩石力学辅导学习的目的。

我们的实验是端部效应对岩石力学性质影响的数值实验，本实验的难度在于多次输入垫板刚度/式样刚度的值不一样。垫板泊松比/式样泊松比也不同，是个难点。在加载步数中观察也是一个很重要的过程，每个人观察的结果和确定的过程也是不同所以造成实验结果大有不同。而出现这些问题之后，我们小组成员互相商量，问题出现在哪里怎样解决问题。

另外，在试验中最后我觉得团队合作也是个很重要的事情，缺乏了团队合作，就缺乏工作效率和准确性。所以团队合作在认识事情中都能彰显出重要意义和作用。

最后，我认为不管是此次实验还是做任何事，都应该一丝不苟，细心对待，以及坚持到底才能获得最终的回报。另外，在工作中也要保持细心认真的态度，才能更好解决工作中的问题和实际问题，让自己体现出价值。

岩石力学-硕士研究生课程报告-中南大学

硕士研究生课程报告题目顺层高边坡稳定性影响因素及工程灾害防治姓名曾义专业班级岩土13级任课教师阳军生张学民中南大学土木工程学院

引言近年来，随着铁路公路建设步伐加快，铁路公路等级不断提高，边坡防护建设工程中所遇到的岩土边坡安全稳定性问题也相应增多，并成为岩土工程中比较常见的技术难题。由于工程建设的需要，往往在一定程度上破坏或扰动原来较为稳定的岩土体而形成新的人工边坡，因而普遍存在着边坡稳定的问题需要解决。国家实施西部大开发战略以来，西部山区高等级公路得到迅速发展。在山区修建高等级公路不可避免会遇到大量的深挖高填路基，就目前建设的高速公路情况看:一般情况下，100km长的山区高等级公路，挖填方路基段落长度占路线总长度的60%以上。已建高速公路最高的填方已达到50多米,最高的挖方边坡高度已超过100m。尽管山区高等级公路的建设越来越倡导环境保护，尽量避免深挖高填，但路基作为公路的主要结构，其边坡稳定问题不可避免。在山区复杂多变的地质条件下建设高等级公路，其边坡稳定性问题必将受到人们的普遍关注，高边坡岩土安全状况直接关系到公路交通运输安全。虽然计算理论方法、地质探测技术、现代监测技术、边坡加固技术及施工技术不断的在进步，但顺层边坡稳定性问题和高边坡稳定性问题，时至今日依然是国内外学者研究的热点问题，并逐步涌现出许多的新的研究方向。 1、顺倾高边坡稳定性研究现状随着人类工程活动的发展，对边坡问题的研究也在不断深入，归纳前人对边坡问题的研究大致可分为以下几个阶段：人们对边坡稳定性的关注和研究最早是从滑坡现象开始的（张倬元等，2001）。19世纪末和20世纪初期，伴随着欧美资本主义国家的工业化而兴起的大规模土木工程建设（如修筑铁路、公路，露天采矿，天然建材开采等），出现了较多的人工边坡，诱发了大量滑坡和崩塌，造成了很大的损失。这时，人们才开始重视边坡失稳给人类造成的危害，并开始借用一般材料分析中的工程力学理论对滑坡进行半经验、半理论的研究。 20世纪50年代，我国学者引进苏联工程地质的体系，继承和发展了“地质历史分析”法，并将其应用于滑坡的分析和研究中，对边坡稳定性研究起到了推动作用（张倬元等，1994）。该阶段学者们着重边坡地质条件的描述和边坡类型的划分，采用工程地质类比法评价边坡稳定性。 20世纪60年代，世界上几起灾难性的边坡失稳事件的发生（如意大利的瓦依昂滑坡造成近3000人死亡和巨大的经济损失）（张倬元等，1994）,使人们逐渐认识到了结构面对边坡稳定性的控制作用以及边坡失稳的时效特征，初步形

数值分析实验报告1

ve=zeros(1,21); ve(n)=ess; r=roots(a+ve); -6-m s=max(abs(r-rr)) end end 利用符号函数：（思考题一）a=poly(1:20); y=poly2sym(a); rr=solve(y) n

很容易的得出对一个多次的代数多项式的其中某一项进行很小的扰动，对其多项式的根会有一定的扰动的，所以对于这类病态问题可以借助于MATLAB来进行问题的分析。学号：06450210 姓名：万轩实验二插值法

岩石力学研究进展报告

岩石力学研究新进展报告姓名：XXX 学号：XXXXXXXX 专业：岩土工程

岩石力学研究新进展报告 1 引言时光如白驹过隙，一学期的《XXXXX》课程在不知不觉间结课了。这一学期的学习，使我在岩石力学方面有了很大的启发，特别是分形理论在岩石力学中的应用令我神往。下面我对岩石力学研究的新进展做简要报告。岩石力学可以作为固体力学的一个新分支，用以研究岩石材料的力学性能和岩石工程的特殊设计方法。岩石力学经过近50年的发展，在土木工程、水利工程、采矿工程、石油工程、国防工程等领域都得到了广泛的应用，随着科学技术的进步，岩石力学涉及的领域会进一步扩大。岩石力学是一门内涵深，工程实践性强的发展中学科。岩石力学面对的是“数据有限”的问题，输入给模型的基本参数很难确定，而且没有多少对过程（特别是非线性工程）的演化提供信息的测试手段。另一方面，对岩体的破坏机体还不能准确的解释。岩石力学所涉及的力学问题是多场（应力场、温度场、渗流场、甚至还存在电磁场等)、多相（固、液、气）影响下的地质构造和工程构造相互作用的耦合问题。这就表明，工程岩体的变形破坏特征是极为复杂的，其大多数是高度非线性的。目前，岩石力学的许多数学模型是不准确和不完整的，可以广泛接受和适用的概化模型并不多。基于此，近年来，多种数值方法、细观力学、断裂与损伤力学、系统科学、分形理论、块体理论等在岩石力学中的应用以及各种人工智能、神经网络、遗传算法、进化算法、非确定性数学等域岩石力学的交叉学科的兴起，为我们提供了全新和有效的思维方式和研究方法，更能激发研究者的创新精神，这也为突破岩石力学的确定性研究方法提供了强有力的理论基础[1]。本报告主要对分形岩石力学、块体岩石力学、断裂与损伤岩石力学和岩石细观力学四部分的研究新进展做简要报告。由于时间和精力有限（最近导师安排的任务非常多，而且要准备英语和政治期末考试），每部分内容除第一大段的研究新进展综述外，只对近几年的三篇比较好的文献做分析说明，包括两篇中文学术论文和一篇外文学术论文，这12篇学术论文我都比较仔细的看了。以后若有机会和时间，我会在导师和各位老师同学的不吝赐教下，努力做岩石力学的创新性研究，届时会在文献综述部分查阅和介绍更多最新以及更优秀的文献。 2 分形岩石力学从古至今，岩石已成为人们熟知的工程材料，它是由矿物晶粒、胶结物质和大量各种不同阶次、不规则分布的裂隙、薄弱夹层等缺陷构成，是一种成分和结构高度复杂的孔隙体。岩石力学经过近50年的发展，人们尝试用各种数学力学方法研究和描述岩石复杂的自然结构性状和物理力学性质，提出了多种岩石力学分析和计算方法，为解决实际工程中的岩石力学问题创造了条件。19世纪70年代Mandelbrot创立分形几何学，提出了一种定量研究和描述自然界中极不规则且看似无序的复杂结构、现象或行为的新方法，从此分形几何学广泛地应用于自然科学研究的各个领域，并且在经济学等社会科学也有很巧妙的应用。19世纪80年代，分形几何学开始应用于岩石力学研究，开始形成分形岩石力学这一门新兴交叉学科。人们逐渐发现岩石力学领域中的分形现象相当普遍，不仅岩石的自然结构性状、缺陷几何形态、分布以及地质结构产状、断层几何形态、分布都观察到分形特征或分形结构，而且岩石体强度、变形、破断力学行为以及能量耗

数值分析实验报告

数值分析实验报告姓名：周茹学号： 912113850115 专业：数学与应用数学指导老师：李建良

线性方程组的数值实验一、课题名字：求解双对角线性方程组二、问题描述考虑一种特殊的对角线元素不为零的双对角线性方程组（以n=7为例） ?????????? ?????? ? ???? ?d a d a d a d a d a d a d 766 55 44 3 32 211??????????????????????x x x x x x x 7654321=?????????? ? ???????????b b b b b b b 7654321 写出一般的n （奇数）阶方程组程序（不要用消元法，因为不用它可以十分方便的解出这个方程组）。三、摘要本文提出解三对角矩阵的一种十分简便的方法——追赶法，该算法适用于任意三对角方程组的求解。四、引言对于一般给定的d Ax =，我们可以用高斯消去法求解。但是高斯消去法过程复杂繁琐。对于特殊的三对角矩阵，如果A 是不可约的弱对角占优矩阵，可以将A 分解为UL ，再运用追赶法求解。

五、计算公式（数学模型）对于形如????? ?? ????? ??? ?---b a c b a c b a c b n n n n n 111 2 2 2 11... ... ...的三对角矩阵UL A =，容易验证U 、L 具有如下形式： ??????? ????? ??? ?=u a u a u a u n n U ...... 3 3 22 1 ， ?? ????? ? ?? ??????=1 (1) 1132 1l l l L 比较UL A =两边元素，可以得到 ? ?? ??-== = l a b u u c l b u i i i i i i 111 i=2, 3, ... ,n 考虑三对角线系数矩阵的线性方程组 f Ax = 这里()T n x x x x ... 2 1 = ,()T n f f f f ... 2 1 = 令y Lx =，则有 f Uy = 于是有 ()?????-== --u y a f y u f y i i i i i 1 1 11 1 * i=2， 3, ... ,n 再根据y Lx =可得到

岩石力学试验报告-2010

长沙理工大学岩石力学试验报告年级班号姓名同组姓名实验日期月日理论课教师：指导教师签字：批阅教师签字：实验一实验二实验三实验四实验五实验六实验七

试验一、岩石单向抗压强度的测定一、试验的目的：测定岩石的单轴抗压强度Rc。当无侧限试样在纵向压力作用下出现压缩破坏时，单位面积上所承受的载荷称为岩石的单轴抗压强度，即试样破坏时的最大载荷与垂直于加载方向的截面积之比。本次试验主要测定天然状态下试样的单轴抗压强度。二、试样制备： 1、试料可用钻孔岩心或坑槽探中采取的岩块。在取料和试样制备过程中，不允许人为裂隙出现。 2、本次试验采用圆柱体作为标准试样，直径为5cm，允许变化范围为4.8～5.4cm，高度为10cm，允许变化范围为9.5～10.5cm。 3、对于非均质的粗粒结构岩石，或取样尺寸小于标准尺寸者，允许采用非标准试样，但高径之比宜为2.0～2.5。 4、制备试样时采用的冷却液，必须是洁净水，不许使用油液。 5、对于遇水崩解、溶解和干缩湿胀的岩石，应采用干法制样。 6、试样数量：每组须制备3个。 7、试样制备的精度。 (1)在试样整个高度上，直径误差不得超过0.3mm。 (2)两端面的不平行度，最大不超过0.05mm。 (3)端面应垂直于试样轴线，最大偏差不超过0.25。三、试样描述：试验前的描述，应包括如下内容： 1、岩石名称、颜色、结构、矿物成分、颗粒大小，风化程度，胶结物性质等特征。 2、节理裂隙的发育程度及其分布，并记述受载方向与层理、片理及节理裂隙之间的关系。 3、量测试样尺寸，检查试样加工精度，并记录试样加工过程中的缺陷。试件压坏后，应描述其破坏方式。若发现异常现象，应对其进行描述和解释。四、主要仪器设备：

数值计算实验报告

(此文档为word格式，下载后您可任意编辑修改！) 2012级6班###(学号)计算机数值方法实验报告成绩册姓名:宋元台学号：成绩：

数值计算方法与算法实验报告学期： 2014 至 2015 第 1 学期 2014年 12月1日课程名称: 数值计算方法与算法专业:信息与计算科学班级 12级5班实验编号： 1实验项目Neton插值多项式指导教师：孙峪怀姓名：宋元台学号：实验成绩：一、实验目的及要求实验目的：掌握Newton插值多项式的算法，理解Newton插值多项式构造过程中基函数的继承特点，掌握差商表的计算特点。实验要求： 1. 给出Newton插值算法 2. 用C语言实现算法二、实验内容三、实验步骤(该部分不够填写.请填写附页)

1.算法分析：下面用伪码描述Newton插值多项式的算法： Step1 输入插值节点数n，插值点序列{x(i),f(i)},i=1,2,……,n,要计算的插值点x. Step2 形成差商表 for i=0 to n for j=n to i f(j)=((f(j)-f(j-1)(x(j)-x(j-1-i)); Step3 置初始值temp=1,newton=f(0) Step4 for i=1 to n temp=(x-x(i-1))*temp*由temp(k)=(x-x(k-1))*temp(k-1)形成 (x-x(0).....(x-x(i-1)* Newton=newton+temp*f(i); Step5 输出f（x）的近似数值newton(x)=newton. 2.用C语言实现算法的程序代码 #includeMAX_N) { printf("the input n is larger than MAX_N,please redefine the MAX_N.\n"); return 1; } if(n<=0) { printf("please input a number between 1 and %d.\n",MAX_N); return 1; } printf("now input the (x_i,y_i)i=0,...%d\n",n); for(i=0;i<=n;i++) { printf("please input x(%d) y(%d)\n",i,i);

数值分析实验报告模板

数值分析实验报告模板篇一：数值分析实验报告(一)(完整) 数值分析实验报告 1 2 3 4 5 篇二：数值分析实验报告实验报告一题目：非线性方程求解摘要：非线性方程的解析解通常很难给出，因此线性方程的数值解法就尤为重要。本实验采用两种常见的求解方法二分法和Newton法及改进的Newton法。利用二分法求解给定非线性方程的根，在给定的范围内，假设f(x,y)在[a,b]上连续，f(a)xf(b) 直接影响迭代的次数甚至迭代的收敛与发散。即若x0 偏离所求根较远，Newton法可能发散的结论。并且本实验中还利用利用改进的Newton法求解同样的方程，且将结果与Newton法的结果比较分析。前言：（目的和意义）掌握二分法与Newton法的基本原理和应用。掌握二分法的原理，验证二分法，在选对有根区间的前提下，必是收

敛，但精度不够。熟悉Matlab语言编程，学习编程要点。体会Newton使用时的优点，和局部收敛性，而在初值选取不当时，会发散。数学原理：对于一个非线性方程的数值解法很多。在此介绍两种最常见的方法：二分法和Newton法。对于二分法，其数学实质就是说对于给定的待求解的方程f(x)，其在[a,b]上连续，f(a)f(b) Newton法通常预先要给出一个猜测初值x0，然后根据其迭代公式xk?1?xk?f(xk) f'(xk) 产生逼近解x*的迭代数列{xk}，这就是Newton法的思想。当x0接近x*时收敛很快，但是当x0选择不好时，可能会发散，因此初值的选取很重要。另外，若将该迭代公式改进为 xk?1?xk?rf(xk) 'f(xk) 其中r为要求的方程的根的重数，这就是改进的Newton 法，当求解已知重数的方程的根时，在同种条件下其收敛速度要比Newton法快的多。程序设计：本实验采用Matlab的M文件编写。其中待求解的方程写成function的方式，如下 function y=f(x);

岩石力学数值试验实验报告

岩石力学数值试验实验报告姓名：郑周立学号： 1108010103 班级：采矿111班指导教师：左宇军同组人：郑周立、周义现、胡斌、朱红伟、高言、王坤实验名称：圆孔对岩石力学性质影响的数值加载试验 2014年5月16日

圆孔对岩石力学性质影响的数值加载试验一、实验目的: 1.通过对RFPA2D学习，知道RFPA2D基本使用方法。 2.了解RFPA2D模拟试验的条件和RFPA2D的基本功能。 3.通过操作端部效应对岩石力学性质影响的数值实验，了解每一步操作以及岩石破裂过程，最终完成实验得到结果。二、实验原理: RFPA-2D是一种基于有限元应力分析和统计损伤理论的材料破裂过程分析数值计算方法，是一个能够模拟材料渐进破裂直至失稳全过程的数值试验工具。三、 1、试样尺寸： 100mm*51mm 2、基元数: 100*51 3、应力分析模式: 平面应变 4、圆孔：半径10mm 5、加载方式：单轴压缩 6、加载条件：竖向位移加载 7、均质度m=2 8、加载量：每步0.002mm

9、实验内容：（1）、应力-应变曲线；（2）、强度；（3）、破坏模式四、实验内容： (一)、操作步骤：第一步启动RFPA，新建模型建立存放的根目录第二步划分网格，单击在弹出的窗口中设置模型的大小，单击确定第三步选择施加荷载模式... （二）实验结果弹性模量图第1步

第4步（开始破坏）第7步（开始横向破坏）第32步（彻底破坏）第200步

最大剪应力图第1步

第4步（开始破坏）第33步（彻底破坏）第200步最大主应力图

数值分析实验报告1

实验一误差分析实验（病态问题）实验目的：算法有“优”与“劣”之分，问题也有“好”与“坏”之别。对数值方法的研究而言，所谓坏问题就是问题本身对扰动敏感者，反之属于好问题。通过本实验可获得一个初步体会。数值分析的大部分研究课题中，如线性代数方程组、矩阵特征值问题、非线性方程及方程组等都存在病态的问题。病态问题要通过研究和构造特殊的算法来解决，当然一般要付出一些代价（如耗用更多的机器时间、占用更多的存储空间等）。问题提出：考虑一个高次的代数多项式 )1.1() ()20()2)(1()(20 1∏=-=---=k k x x x x x p 显然该多项式的全部根为1,2,…,20共计20个，且每个根都是单重的。现考虑该多项式的一个扰动 )2.1(0 )(19=+x x p ε 其中ε是一个非常小的数。这相当于是对（）中19x 的系数作一个小的扰动。我们希望比较（）和（）根的差别，从而分析方程（）的解对扰动的敏感性。实验内容：为了实现方便，我们先介绍两个Matlab 函数：“roots ”和“poly ”。 roots(a)u = 其中若变量a 存储n+1维的向量，则该函数的输出u 为一个n 维的向量。设a 的元素依次为121,,,+n a a a ，则输出u 的各分量是多项式方程 01121=+++++-n n n n a x a x a x a 的全部根；而函数 poly(v)b =

的输出b 是一个n+1维变量，它是以n 维变量v 的各分量为根的多项式的系数。可见“roots ”和“poly ”是两个互逆的运算函数。 ;000000001.0=ess );21,1(zeros ve = ;)2(ess ve = ))20:1((ve poly roots + 上述简单的Matlab 程序便得到（）的全部根，程序中的“ess ”即是（）中的ε。实验要求：（1）选择充分小的ess ，反复进行上述实验，记录结果的变化并分析它们。如果扰动项的系数ε很小，我们自然感觉（）和（）的解应当相差很小。计算中你有什么出乎意料的发现表明有些解关于如此的扰动敏感性如何（2）将方程（）中的扰动项改成18x ε或其它形式，实验中又有怎样的现象出现（3）（选作部分）请从理论上分析产生这一问题的根源。注意我们可以将方程（）写成展开的形式， ) 3.1(0 ),(1920=+-= x x x p αα 同时将方程的解x 看成是系数α的函数，考察方程的某个解关于α的扰动是否敏感，与研究它关于α的导数的大小有何关系为什么你发现了什么现象，哪些根关于α的变化更敏感思考题一：（上述实验的改进）在上述实验中我们会发现用roots 函数求解多项式方程的精度不高，为此你可以考虑用符号函数solve 来提高解的精确度，这需要用到将多项式转换为符号多项式的函数poly2sym,函数的具体使用方法可参考Matlab 的帮助。

(完整版)哈工大-数值分析上机实验报告

实验报告一题目：非线性方程求解摘要：非线性方程的解析解通常很难给出，因此线性方程的数值解法就尤为重要。本实验采用两种常见的求解方法二分法和Newton法及改进的Newton法。前言：（目的和意义）掌握二分法与Newton法的基本原理和应用。数学原理：对于一个非线性方程的数值解法很多。在此介绍两种最常见的方法：二分法和Newton法。对于二分法，其数学实质就是说对于给定的待求解的方程f(x)，其在[a,b]上连续，f(a)f(b)<0，且f(x)在[a,b]内仅有一个实根x*，取区间中点c，若，则c恰为其根，否则根据f(a)f(c)<0是否成立判断根在区间[a,c]和[c,b]中的哪一个，从而得出新区间，仍称为[a,b]。重复运行计算，直至满足精度为止。这就是二分法的计算思想。

Newton法通常预先要给出一个猜测初值x0，然后根据其迭代公式产生逼近解x*的迭代数列{x k}，这就是Newton法的思想。当x0接近x*时收敛很快，但是当x0选择不好时，可能会发散，因此初值的选取很重要。另外，若将该迭代公式改进为其中r为要求的方程的根的重数，这就是改进的Newton法，当求解已知重数的方程的根时，在同种条件下其收敛速度要比Newton法快的多。程序设计：本实验采用Matlab的M文件编写。其中待求解的方程写成function的方式，如下 function y=f(x); y=-x*x-sin(x); 写成如上形式即可，下面给出主程序。二分法源程序： clear %%%给定求解区间 b=1.5; a=0;

%%%误差 R=1; k=0;%迭代次数初值 while (R>5e-6) ; c=(a+b)/2; if f12(a)*f12(c)>0; a=c; else b=c; end R=b-a;%求出误差 k=k+1; end x=c%给出解 Newton法及改进的Newton法源程序：clear %%%% 输入函数 f=input('请输入需要求解函数>>','s') %%%求解f(x)的导数 df=diff(f);

岩石力学试验报告

岩石力学实验指导书及实验报告班级姓名山东科技大学土建学院实验中心编

目录一、岩石比重的测定二、岩石含水率的测定三、岩石单轴抗压强度的测定四、岩石单轴抗拉强度的测定五、岩石凝聚力及内摩擦角的测定（抗剪强度试验）六、岩石变形参数的测定七、煤的坚固性系数的测定

实验一、岩石比重的测定岩石比重是指单位体积的岩石（不包括孔隙）在105～110o C 下烘至恒重的重量与同体积4o C 纯水重量的比值。一、仪器设备岩石粉碎机、瓷体或玛瑙体、孔径0.2或0.3毫米分样筛、天平（量0.001克）、烘箱、干燥器、沙浴、比重瓶。二、试验步骤 1、岩样制备：取有代表性的岩样300克左右，用机械粉碎，并全部通过孔径0.2（或0.3）毫米分样筛后待用。 2、将蒸馏水煮沸并冷却至室温取瓶颈与瓶塞相符的100毫升比重瓶，用蒸馏水洗净，注入三分之一的蒸馏水，擦干瓶的外表面。 3、取15g 岩样（称准到0.001克）得g 借助漏斗小心倒入盛有三分之一蒸馏水的比重瓶中，注意勿使岩样抛撒或粘在瓶颈上。 4、将盛有蒸馏水和岩样的比重瓶放在沙浴上煮沸后再继续煮1～1.5小时。 5、将煮沸后的比重瓶自然冷却至室温，然后注入蒸馏水，使液面与瓶塞刚好接触，注意不得留有气泡，擦干瓶的外表面，在天平上称重得g 1。 6、将岩样倒出，比重瓶洗净，最后用蒸馏水刷一遍，向比重瓶内注满蒸馏水，同样使液面与瓶塞刚好接触，不得留有气泡，擦干瓶的外表面，在天平上称重得g 2。三、结果：按下式计算： s d g g g g d 1 2-+= 式中：d ——岩石比重； g ——岩样重、克； g 1——比重瓶、岩样和蒸馏水合重、克； g 2——比重瓶和满瓶蒸馏水合重、克； d s ——室温下蒸馏水的比重、d s ≈1

数值分析实验报告

实验一、误差分析一、实验目的 1．通过上机编程，复习巩固以前所学程序设计语言及上机操作指令； 2．通过上机计算，了解误差、绝对误差、误差界、相对误差界的有关概念； 3．通过上机计算，了解舍入误差所引起的数值不稳定性。二．实验原理误差问题是数值分析的基础，又是数值分析中一个困难的课题。在实际计算中，如果选用了不同的算法，由于舍入误差的影响，将会得到截然不同的结果。因此，选取算法时注重分析舍入误差的影响，在实际计算中是十分重要的。同时，由于在数值求解过程中用有限的过程代替无限的过程会产生截断误差，因此算法的好坏会影响到数值结果的精度。三．实验内容对20,,2,1,0 =n ，计算定积分 ?+=10 5dx x x y n n . 算法1：利用递推公式 151--=n n y n y , 20,,2,1 =n , 取 ?≈-=+=1 00182322.05ln 6ln 51dx x y . 算法2：利用递推公式 n n y n y 51511-= - 1,,19,20 =n . 注意到 ???=≤+≤=10 10202010201051515611261dx x dx x x dx x , 取 008730.0)12611051(20120≈+≈y .: 四.实验程序及运行结果程序一： t=log(6)-log(5);

n=1; y(1)=t; for k=2:1:20 y(k)=1/k-5*y(k-1); n=n+1; end y y =0.0884 y =0.0581 y =0.0431 y =0.0346 y =0.0271 y =0.0313 y =-0.0134 y =0.1920 y =-0.8487 y =4.3436 y =-21.6268 y =108.2176 y =-541.0110 y =2.7051e+003 y =-1.3526e+004 y =6.7628e+004 y =-3.3814e+005 y =1.6907e+006 y =-8.4535e+006 y =4.2267e+007 程序2： y=zeros(20,1); n=1; y1=(1/105+1/126)/2;y(20)=y1; for k=20:-1:2 y(k-1)=1/(5*k)-(1/5)*y(k); n=n+1; end 运行结果：y = 0.0884 0.0580 0.0431 0.0343 0.0285 0.0212 0.0188 0.0169

中南大学ANSYS上机实验报告

ANSYS上机实验报告小组成员：郝梦迪、赵云、刘俊一、实验目的和要求本课程上机练习的目的是培养学生利用有限单元法的商业软件进行数值计算分析，重点是了解和熟悉ANSYS的操作界面和步骤，初步掌握利用ANSYS建立有限元模型，学习ANSYS分析实际工程问题的方法，并进行简单点后处理分析，识别和判断有限元分析结果的可靠性和准确性。二、实验设备和软件台式计算机，ANSYS10.0软件三、基本步骤 1）建立实际工程问题的计算模型。实际的工程问题往往很复杂，需要采用适当的模型在计算精度和计算规模之间取得平衡。常用的建模方法包括：利用几何、载荷的对称性简化模型，建立等效模型。 2）选择适当的分析单元，确定材料参数。侧重考虑一下几个方面：是否多物理耦合问题，是否存在大变形，是否需要网格重划分。 3）前处理（Preprocessing）。前处理的主要工作内容如下：建立几何模型（Geometric Modeling），单元划分（Meshing）与网格控制，给定约束（Constraint）和载荷（Load）。在多数有限元软件中，不能指定参数的物理单位。用户在建模时，要确定力、长度、质量及派生量的物理单位。在建立有限元模型时，最好使用统一的物理单位，这样做不容易弄错计算结果的物理单位。建议选用kg，N，m，sec；常采用kg，N，mm，sec。 4）求解（Solution）。选择求解方法，设定相应的计算参数，如计算步长、迭代次数等。 5）后处理（Postprocessing）。后处理的目的在于确定计算模型是否合理、计算结果是否合理、提取计算结果。可视化方法（等值线、等值面、色块图）显

数值分析2016上机实验报告

序言数值分析是计算数学的范畴，有时也称它为计算数学、计算方法、数值方法等，其研究对象是各种数学问题的数值方法的设计、分析及其有关的数学理论和具体实现的一门学科，它是一个数学分支。是科学与工程计算（科学计算）的理论支持。许多科学与工程实际问题（核武器的研制、导弹的发射、气象预报）的解决都离不开科学计算。目前，试验、理论、计算已成为人类进行科学活动的三大方法。数值分析是计算数学的一个主要部分,计算数学是数学科学的一个分支,它研究用计算机求解各种数学问题的数值计算方法及其理论与软件实现。现在面向数值分析问题的计算机软件有：C,C++,MATLAB,Python,Fortran等。 MATLAB是matrix laboratory的英文缩写，它是由美国Mathwork公司于1967年推出的适合用于不同规格计算机和各种操纵系统的数学软件包，现已发展成为一种功能强大的计算机语言，特别适合用于科学和工程计算。目前，MATLAB应用非常广泛，主要用于算法开发、数据可视化、数值计算和数据分析等，除具备卓越的数值计算能力外，它还提供了专业水平的符号计算，文字处理，可视化建模仿真和实时控制等功能。本实验报告使用了MATLAB软件。对不动点迭代，函数逼近（lagrange插值，三次样条插值，最小二乘拟合），追赶法求解矩阵的解，4RungeKutta方法求解，欧拉法及改进欧拉法等算法做了简单的计算模拟实践。并比较了各种算法的优劣性，得到了对数值分析这们学科良好的理解，对以后的科研数值分析能力有了极大的提高。

目录序言 (1) 问题一非线性方程数值解法 (3) 1.1 计算题目 (3) 1.2 迭代法分析 (3) 1.3计算结果分析及结论 (4) 问题二追赶法解三对角矩阵 (5) 2.1 问题 (5) 2.2 问题分析（追赶法） (6) 2.3 计算结果 (7) 问题三函数拟合 (7) 3.1 计算题目 (7) 3.2 题目分析 (7) 3.3 结果比较 (12) 问题四欧拉法解微分方程 (14) 4.1 计算题目 (14) 4.2.1 方程的准确解 (14) 4.2.2 Euler方法求解 (14) 4.2.3改进欧拉方法 (16) 问题五四阶龙格-库塔计算常微分方程初值问题 (17) 5.1 计算题目 (17) 5.2 四阶龙格-库塔方法分析 (18) 5.3 程序流程图 (18) 5.4 标准四阶Runge-Kutta法Matlab实现 (19) 5.5 计算结果及比较 (20) 问题六舍入误差观察 (22) 6.1 计算题目 (22) 6.2 计算结果 (22) 6.3 结论 (23) 7 总结 (24) 附录

岩层实验报告

中国矿业大学矿业工程学院实验报告

《岩层控制》实验报告实验一矿山岩体力学实验注：包括岩石抗拉、抗压、抗剪三个内容。岩石的抗拉强度试验一、实验目的与要求岩石在单轴拉伸载荷作用下达到破坏时所能承受的最大拉应力称为岩石的单轴抗拉强度。由于进行直接拉伸实验在准备试件方面要花费大量的人力、物力和时间，因此采用间接拉伸实验方法来测试岩石的抗拉强度。劈裂法是最基本的方法。二、实验仪器（1）钻石机或车床，锯石机，磨石机或磨床。（2）劈裂法实验夹具，或直径2.0mm钢丝数根。（3）游标卡尺（精度0.02mm），直角尺，水平检测台，百分表架和百分表。（4）材料实验机。三、实验原理图3-1显示的是在压应力作用下，沿圆盘直径y-y的应力分布图。在圆盘边缘处，沿y-y方向（σy）和垂直y-y（σx）方向均为压应力，而离开边缘后，沿y-y方向仍为压应力，但应力值比边缘处显著减少，并趋于平均化；垂直y-y方向变成拉应力。并在沿y-y的很长一段距离上呈均匀分布状态。虽然拉应力的值比压应力值低很多，但由于岩石的抗拉强度很低，所以试件还是由于x方向的拉应力而导致试件沿直径的劈裂破坏，破坏是从直径中心开始，然后向两端发展，反映了岩石的抗拉强度比抗压强度要低得多的事实。 χy r/R 0.5 -0.5x σyσx y 压缩拉伸应力值/MPa 160120804040 图3-1 劈裂实验应力分布示意图四、实验内容

(1) 了解试件的加工机具、检测机具，规程对精度的要求及检测方法； (2) 学会材料实验机的操作方法及拉压夹具的使用方法； (3) 学会间接测试岩石抗压强度及数据处理方法。五、实验步骤（1）测定前核对岩石名称和岩样编号，对试件颜色、颗粒、层理、裂隙、风化程度、含水状态机加工过程中出现的问题进行描述，并填入记录表1-1内。（2）检查试件加工精度，测量试件尺寸，填入记录表内。（3）选择材料实验机度盘时，一般应满足下式：0.2 P 0< P max <0.8P 0 （4）通过试件直径两端，沿轴线方向画两条互相平行的线作为加载基线。把试件放入夹具内，夹具上、下刀刃对准加载基线，用两侧夹持螺钉固定好试件，或用两根直径2.0mm 的钢丝放在加载基线上，钢丝间用橡皮筋固定。（5）把夹好试件的夹具或夹好钢丝的试件放入材料实验机的上、下承压板之间，使试件的中心线和材料实验机的中心线在一条直线上。（6）开动材料实验机，施加数百牛载荷后，松开夹具两侧夹持螺钉，然后以0.03~0.05MPa/s 的速度加载，直至试件破坏。（7）记录破坏载荷，对破坏后的试件进行摄影或描述。六、注意事项（1）记录试件的完整状态，（2）选择合适的材料实验机及合适的实验机度盘值，（3）夹具对试件的加载方向要与试件的轴线在一平面上，（4）选择合适的加载速率。七、数据处理表1-1 计算试件单向抗拉强度： R 1= 102?DL P π=5.98MPa 式中 R 1—试件的抗拉强度，MPa ； P —试件破坏载荷，kN; D —试件直径，cm; L —试件厚度，cm 。八、误差分析（1）试件自身各方面的影响；（2）系统误差；

岩体力学实验..

岩体力学实验一.实验目的岩石单轴压缩是指岩石在单轴压缩条件下的强度、变形和破坏特征。通过该实验掌握岩石单轴压缩实验方法，学会岩石单轴抗压强度、弹性模量、泊松比的计算方法；了解岩石单轴压缩过程的变形特征和破坏类型。二.实验设备、仪器和材料 1.钻石机、锯石机、磨石机； 2.游标卡尺，精度0.02mm； 3.直角尺、水平检测台、百分表及百分表架； 4.YE-600型液压材料试验机； 5.JN-16型静态电阻应变仪； 6.电阻应变片（BX-120型）； 7.胶结剂，清洁剂，脱脂棉，测试导线等。三.试样的规格、加工精度、数量及含水状态 1. 试样规格：采用直径为50 mm，高为100 mm的标准圆柱体，对于一些裂隙比较发育的试样，可采用50 mm×50 mm×100 mm的立方体，由于岩石松软不能制取标准试样时，可采用非标准试样，需在实验结果加以说明。 2. 加工精度： a 平行度：试样两端面的平行度偏差不得大于0.1mm。检测方法如图5-1所示，将试样放在水平检测台上，调整百分表的位置，使百分表触头紧贴试样表面，然后水平移动试样百分表指针的摆动幅度小于10格。 b 直径偏差：试样两端的直径偏差不得大于0.2 mm，用游标卡尺检查。 c 轴向偏差：试样的两端面应垂直于试样轴线。检测方法如图5-2所示，将试样放在水平检测台上，用直角尺紧贴试样垂直边，转动试样两者之间无明显

缝隙。 3.试样数量：每种状态下试样的数量一般不少于3个。 4.含水状态：采用自然状态，即试样制成后放在底部有水的干燥器内1～2 d ，以保持一定的湿度，但试样不得接触水面。四.电阻应变片的粘贴 1.阻值检查：要求电阻丝平直，间距均匀，无黄斑，电阻值一般选用120欧姆，测量片和补偿片的电阻差值不超过0.5Ω。 2.位置确定：纵向、横向电阻应变片粘贴在试样中部，纵向、横向应变片排列采用“┫”形，尽可能避开裂隙，节理等弱面。 3.粘贴工艺：试样表面清洗处理→涂胶→贴电阻应变片→固化处理→焊接导线→防潮处理。五.实验步骤 1. 测定前核对岩石名称和试样编号，并对岩石试样的颜色、颗粒、层理、裂隙、风化程度、含水状态等进行描述。 2. 检查试样加工精度。并测量试样尺寸，一般在试样中部两个互相垂直方向测量直径计算平均值。 3. 电阻应变仪接通电源并预热数分钟后, 连接测试导线，接线方式采用公 1—百分表 2-百分表架 3-试样 4水平检测台图5-1 试样平行度检测示意图 1—直角尺 2-试样 3- 水平检测台图5-2 试样轴向偏差度检测示意图图5-3 电阻应变片粘贴

数值分析实验报告

学生实验报告实验课程名称开课实验室学院年级专业班学生姓名学号开课时间至学年学期

if(A(m,k)~=0) if(m~=k) A([k m],:)=A([m k],:); %换行 end A(k+1:n, k:c)=A(k+1:n, k:c)-(A(k+1:n,k)/ A(k,k))*A(k, k:c); %消去end end x=zeros(length(b),1); %回代求解 x(n)=A(n,c)/A(n,n); for k=n-1:-1:1 x(k)=(A(k,c)-A(k,k+1:n)*x(k+1:n))/A(k,k); end y=x; format short;%设置为默认格式显示,显示5位（2）建立MATLAB界面利用MA TLAB的GUI建立如下界面求解线性方程组：详见程序。五、计算实例、数据、结果、分析下面我们对以上的结果进行测试，求解：

? ? ? ? ? ? ? ? ? ? ? ? - = ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? - - - - - - 7 2 5 10 13 9 14 4 4 3 2 1 13 12 4 3 3 10 2 4 3 2 1 x x x x 输入数据后点击和，得到如下结果：更改以上数据进行测试，求解如下方程组： 1 2 3 4 43211 34321 23431 12341 x x x x ?? ???? ?? ???? ?? ???? = ?? ???? - ?? ???? - ???? ?? 得到如下结果：

数值分析实验报告1

实验一误差分析实验1.1（病态问题）实验目的：算法有“优”与“劣”之分，问题也有“好”与“坏”之别。对数值方法的研究而言，所谓坏问题就是问题本身对扰动敏感者，反之属于好问题。通过本实验可获得一个初步体会。数值分析的大部分研究课题中，如线性代数方程组、矩阵特征值问题、非线性方程及方程组等都存在病态的问题。病态问题要通过研究和构造特殊的算法来解决，当然一般要付出一些代价（如耗用更多的机器时间、占用更多的存储空间等）。问题提出：考虑一个高次的代数多项式 )1.1() ()20()2)(1()(20 1∏=-=---=k k x x x x x p 显然该多项式的全部根为1,2,…,20共计20个，且每个根都是单重的。现考虑该多项式的一个扰动 )2.1(0 )(19=+x x p ε 其中ε是一个非常小的数。这相当于是对（1.1）中19x 的系数作一个小的扰动。我们希望比较（1.1）和（1.2）根的差别，从而分析方程（1.1）的解对扰动的敏感性。实验内容：为了实现方便，我们先介绍两个Matlab 函数：“roots ”和“poly ”。 roots(a)u = 其中若变量a 存储n+1维的向量，则该函数的输出u 为一个n 维的向量。设a 的元素依次为121,,,+n a a a ，则输出u 的各分量是多项式方程 01121=+++++-n n n n a x a x a x a 的全部根；而函数 poly(v)b = 的输出b 是一个n+1维变量，它是以n 维变量v 的各分量为根的多项式的系数。可见“roots ”和“poly ”是两个互逆的运算函数。 ;000000001.0=ess );21,1(zeros ve = ;)2(ess ve =

现场岩石力学试验报告模板

工程勘察：证书编号 45040Ⅲ -211-U 桂林漓江**水库枢纽工程现场岩石试验报告广西*******勘察设计研究院

核定：审查：校核：编写：试验：

1工作概况 (1) 2 现场混凝土与岩体抗剪（断）试验 (1) 2.1 抗剪(断)试验试样布置及地质条件 (1) 2.2 抗剪（断）试验试样制备情况 (2) 2.3 抗剪(断)试验方法 (2) 2.4 抗剪(断)试验成果整理方法 (3) 2.5 抗剪(断)试验破坏机理分析 (3) 2.6 抗剪断试验成果分析 (4) 3 现场岩体变形试验 (5) 3.1 岩体变形试验试样布置及地质条件 (7) 3.2 岩体变形试点制作 (7) 3.3 岩体变形试验方法 (7) 3.4 岩体变形试验成果整理 (7) 3.5 岩体变形试验成果分析 (8) 4 建议 (9)

1 工作概况桂林漓江**水库枢纽工程位于广西桂林市为漓江一级支流，距离桂林**km有等外公路从**至**村。该水库枢纽主要任务是调蓄讯期洪水水量，枯水期向漓江补水，并利用补水水能发电。拟建枢纽最大坝高约**m，正常高水位**m，总库容约为**万m3，通过引水隧洞到下游厂房发电，电站装机容量为**MW。坝址现场岩体力学试验于****日至*****日坝轴线左岸及坝轴线下游200m右岸进行现场混凝土与岩体抗剪（断）试验及现场岩体变形试验，共完成工作量见表1。表1 现场岩石试验工作量表试验数据采集和处理采用8098多功能岩土检测系统，该微机系统于1991年4月通过广西科学技术委员会的技术鉴定，开工前经广西计量测试研究所率定。各项技术指标均符合DLJ204-81，SLJ2-81《水利水电工程岩石试验规程》（试行），DL5006-92《水利水电工程岩石试验规程（补充部分）》。 2 现场混凝土与岩体抗剪（断）强度试验 2.1抗剪(断)试验试样布置及地质条件 a) 现场混凝土与岩体抗剪（断）试验在坝址区内进行，分别选强、弱风化泥质粉砂岩各12个点（即3组），详见表2。岩层产状一般为**?/NW∠**?,周围岩石为砂岩、泥岩互层。