来源：中北大学 内容包括： 线性方程组的解法 矩阵特征值与特征向量的计算 非线性方程与方程组组的解法 插值与逼近 数值积分 微分方程的数值解法

为了研究岩石的Ⅰ型动态断裂韧度测试方法,用大直径(100 mm)分离式霍普金森压杆(SHPB)对边裂纹平台圆环(edge cracked flattened ring,ECFR)砂岩试样进行径向冲击实验。分别用普通应变片和高精度裂纹扩展计2种测试元件监测试样起裂时刻和测定裂纹传播速度,比较了它们的准确性和合理性。用实验-数值结合普适函数分析测定了砂岩的动态起裂韧度和动态扩展韧度,初步探讨了动态加载率影响动态起裂韧度的原因,以及裂纹扩展速度对动态扩展韧度的影响,对ECFR试样裂纹扩展路径弯折现象也进行了理论分析。研究结果表明:裂纹扩展计测定比应变片更为灵敏、准确、合理。加载率在(0.74~4.48)×104MPa·m1/2/s范围内动态起裂韧度随动态加载率的增大而增大,裂纹扩展速度在(0.24~0.34)cR范围内动态扩展韧度随裂纹扩展速度的增大而提高。

库水位下降的饱和-非饱和土坝渗流数值模拟，徐杨军，夏元友，基于非饱和土渗流的理论，借助Geo-slope的有限元数值模拟软件，对库水位下降过程中的饱和-非饱和土坝进行了数值模拟分析，分析结果�

为揭示天然气水合物试采工程中储层物理参数演化规律,建立了耦合化学分解和气、水两相渗流物理场的数学模型,并采用有限单元法进行求解,对不同降压幅度、初始绝对渗透率和含水饱和度等影响下的储层物理参数的时空演化规律进行研究。结果表明:井筒开始降压后,井周迅速形成压降漏斗,压降辐射范围随时间逐渐增大,且水合物分解范围和压降范围呈现明显的一致性;井筒降压后短时间内产气速率达到峰值,该峰值随降压幅度和初始绝对渗透率的增大而增大,随初始含水饱和度的增大而减小;水合物分解速率和渗透率增大速率与降压幅度和初始绝对渗透率成正相关,与初始含水饱和度成负相关。

为研究深立井基岩段地面预注浆浆液渗流规律,以淮南某矿改扩建工程第二副井基岩段地面预注浆工程作为研究对象,将浆液裂隙流等效为多孔介质渗流,建立考虑注浆压力与岩石变形耦合效应的流固耦合注浆模型,运用数值模拟方法,分析裂隙基岩段地面预注浆浆液扩散半径和注浆压力与岩石渗透系数之间的关系。结果表明在其他因素保持不变的情况下,多孔注浆注入浆液量和浆液的扩散半径随注浆压力和渗透系数的增加而逐渐增加。揭示深立井裂隙基岩段地面预注浆浆液渗流规律,并成功应用于工程实践,为今后类似工程注浆参数设计与优化提供了参考。

为实现具有防渗墙防渗措施的河道堤防二维渗流场真实区域的数值模拟,采用有限体积法的Fluent数值模拟软件,利用其多孔介质模型,求解堤防基础的渗流场,并对采取防渗墙的单一防渗措施条件下的渗流场进行分析,结果表明,随着防渗墙深度的增加,其对河道堤防防渗效果明显但不是无限制的,由此可为堤防渗流计算及工程所需的合理的防渗墙深度提供依据。

应用FLUENT软件建立了综放面采空区瓦斯渗流的三维数值模型,通过数值模拟对比分析了没有瓦斯抽采和有瓦斯尾巷抽排两种情况下的综放面采空区及其上覆岩层内的瓦斯渗流分布规律,通过现场验证,分析了瓦斯尾巷在实际应用中的效果,模拟结果和实际效果较吻合,使工作面达到了安全生产的要求。

转子碰摩系统的非线性数值分析，王璋奇，杨文刚，建立了基于瞬时碰摩力的偏心转子碰摩模型，该碰摩模型为分段线性系统，数值求解该类系统需采用接缝法，但是直接的数值方法无法准

目前,多数矿山采用FLAC3D进行边坡稳定数值模拟,然而在使用FLAC3D软件时,会发现建模比较困难,并浪费大量时间。结合工程实例,介绍CASS与ANSYS的FLAC3D非工作帮稳定性分析,进而分析3种软件在数值模拟中各自的作用,并分析边坡的稳定性。

通过对磨粒流加工原理和运动规律的分析,探讨喷油嘴微小孔磨粒流加工特性。以流体力学软件FLUENT为平台,采用固液两相Mixture湍流模型对磨粒流加工过程中磨粒流的流动形态进行数值模拟,模拟结果为磨粒流加工参数的优化选择提供了理论依据。