岩石多场耦合作用的研究是当前研究的热点难点问题，为了更好的分析岩石在多场耦合作用条件下的作用机理，主要通过实验和数值模拟两方面进行研究。在总结国内外多场耦合微观–细观–宏观多尺度力学试验设备的改进和研发、数值模拟软件及耦合计算程序的开发等方面的研究现状的基础上，展望多场多相耦合作用下岩石力学实验设备和数值分析的研究方向。为了研究岩石多场耦合作用下的力学性能，通过改进和研发设计了不同物理场多场耦合试验系统，在开发试验设备的基础上引起和发展现代无损探测手段，比如实时CT（Computed tomography）扫描技术，电镜扫描技术（SEM）、核磁共振技术（NMRI）、X射线立体成像法、超声波技术等，既能无损检测到岩石的内部孔隙微细观结构及演化过程，也能得出岩石在温度?水流?应力?化学（THMC）多场耦合作用中各物理场的宏观关系，进一步从微细观和宏观相结合的角度得出岩石在多场耦合作用下的性能。随着计算机技术的进步，岩石多场耦合作用下的数值模拟软件及耦合计算程序的开发有了一定的发展，特别是TOUGHREACT与FLAC3D相结合的THMC四场耦合作用的数值模拟软件和数值仿真软件Comsol与Matlab对接的多场耦合计算程序的开发，为岩石多场耦合模拟的开展提供了技术支持

在微观力学行为分析的基础上，对90W合金宏观力学性能及其与微观结构因素（粘结相力学参数）之间的关系进行了计算机数值模拟研究.结果表明：钨合金性能与粘结相力学参数密切相关.随着粘结相弹性模量增加，合金的抗拉强度增加，但延伸率降低.当粘结相屈服强度800MPa时，合金抗拉强度随粘结相屈服强度增加而增大，在粘结相屈服度为800MPa时达到最大值.随粘结相抗拉强度增加，合金抗拉强度和延伸率均呈近似线性规律增加.合金延伸率对粘结相应变硬化模量极为敏感

提出了钨合金的组织结构模型.在此基础上通过有限元法对钨合金拉伸变形时组织结构中的应力分布及其变化规律进行了计算机数值模拟研究.结果表明：钨合金拉伸变形时拉伸方向最大正应力和最大Mises应力区域在钨颗粒相中，同时钨颗粒还将受到垂直于拉伸方向的压应力；最大剪切应力分布在粘结相中，随拉伸变形增加，粘结相最先进入塑性状态；由于粘结相的塑性变形，钨合金中应力不断重组，应力逐渐在钨颗粒中集中

第一章 概论 第二章 MATLAB的基础准备及入门 第三章 MATLAB的数值计算功能 第四章 MATLAB的符号计算功能 第五章 计算结果的可视化 第六章 MATLAB Notebook 第七章 SIMULINK动态仿真集成环境 第八章 MATLAB的程序设计

方程是在科学研究中不可缺少的工具 方程求解是科学计算中一个重要的研究对象 几百年前就已经找到但是对于更高次数 了代数方程中二次至的代数方程目前仍 五次方程的求解公式无有效的精确解法 对于无规律的非代数方程的求解也无精确解法 因此研究非线性方程的数值解法成为必然

第一章概论 计算流体力学在近二三十年中有了突飞猛进的发展,而且正 在以更快的速度前进。推动这一发展的原因一方面是实际问题 的需要,特别是宇航事业的需要另一方面是计算技术的飞速发 展和巨型计算机的出现。 计算流体力学是多种领域的交叉学科,它所涉及的学科有流 体力学、偏微分方程的数学理论、计算何、数值分析、计算机 科学等。它的发展促进了这些学科的进一步发展。最终体现计算 流体力学水平的是解决实际问题的能力

一、牛顿法及其几何意义 二、 收敛性及其收敛速度 三、 计算实例及其程序演示

第 1 章 基础准备及入门 第 2 章 符号计算 第 3 章 数值数组及向量化运算 第 4 章 数值计算 第 5 章 数据和函数的可视化 第 6 章 M 文件和函数句柄 第 7 章 SIMULINK 交互式仿真集成环境 第 8 章 图形用户界面（GUI） 第 9 章 Notebook

我们知道对矩阵进行一次初等变换,就相 当于用相应的初等矩阵去左乘原来的矩阵。 因此我们这个观点来考察Gaus消元法用 矩阵乘法来表示,即可得到求解线性方程 组的另一种直接法:矩阵的三角分解

定理7.3.1设矩阵A∈Rn,且非奇异,则一定存在正交矩 nxn 阵,上三角矩阵R,使 A=OR (7.3.2) 且当要求R的主对角元素均为正数时,则分解式(7.3.2)是唯一的。 证明存在性有矩阵A的非奇 Householder异性及变换矩 阵的性质(3)知,一定可构造n-1个H矩阵:H1,H2,…,Hn-1使 A+1=HA(k=1,2,…n-1)