数值模拟上机实习报告.DOCVIP

《数值模拟技术》课程报告 PAGE \* MERGEFORMAT 1 第一章 数值模拟技术研究文献综述 本门课详细地讲解了数值模拟的概念、原理、应用和几种重要的数值模拟方法的应用，主要包括有限差分法、有限元法和离散元法等。 一般而言，岩、土体处于三向受力状态，其破坏模式往往表现为压-剪破坏和拉伸破坏。要分析和预测岩、土体在外力作用下的变形、破坏，就需要对其变形、破坏情况进行较为直观地再现。岩土工程数值模拟正是从岩、土体的受力状态出发，来分析和预测岩、土体破坏情况的一种手段。 其基本原理是以典型试样的物理试验（室内试验或现场试验）获得的强度来表征整个地质体的岩、土体强度，以边界条件替代地质体周围所受的约束条件，借由本构关系表达岩、土体在外力作用下的应力-应变特性，最终了解、预测岩、土体变形破坏情况。它具有鲜明的时代特征，以计算机为实现平台，是信息化时代的产物。它当前的发展趋势是，通过与其它方法（如人工智能、人工生命科学、随机模拟、模糊数学、灰色理论以及分形理论等）交叉共生、相互耦合嫁接，以获得更广阔的发展空间。 从广义上来说，岩、土体的室内试验和原位试验也是一种模拟手段，可称之为物理模拟。之所以如此称谓，是因为它们也是为较真实地近似再现岩、土体在其所赋存的环境中所处的受力状态所采用的一种手段。从这个意义上来说，它与数值模拟的基本原理是相同的，因此，可以将数值模拟称为虚拟实验室模拟。所不同的是，数值模拟除可以进行常规尺寸模型的模拟外，还可以进行宏观和细观两个层面尺寸模型的模拟，而其输入的参数则需通过物理模拟来提供。因此，数值模拟是与物理模拟并行发展、相互补充和相互验证的试验系统。 相较于其它方法，数值模拟具有可重复和操作性强，费用低廉，不受模型尺寸控制，可视化程度高的优点，能有效延伸和扩展分析人员的认知范围，为分析人员洞悉岩、土体内部的破坏机理提供了强有力的可视化手段。 当然，作为一种分析方法，它也有自身的缺点，主要是易受制于岩、土体结构的描述和模型概化的准确性及合理性；受制于岩、土体物理试验模拟结果的准确性；受制于岩、土体本构关系与实际岩、土体力学响应特性拟合程度的高低。 下面对课堂所学的数值模拟方法进行总结： 1、有限差分方法(FDM)是计算机数值模拟最早采用的方法，至今仍被广泛运用。该方法将 求解域划分为差分网格，用有限个网格节点代替连续的求解域。有限差分法以Taylor级数展开等方法，把控制方程中的导数用网格节点上的函数值的差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法，数学概念直观，表达简单，是发展较早且比较成熟的数值方法。 对于有限差分格式，从格式的精度来划分，有一阶格式、二阶格式和高阶格式。从差分的空间形式来考虑，可分为中心格式和逆风格式。考虑时间因子的影响，差分格式还可以分为显格式、隐格式、显隐交替格式等。目前常见的差分格式，主要是上述几种形式的组合，不同的组合构成不同的差分格式。差分方法主要适用于有结构网格，网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。构造差分的方法有多种形式，目前主要采用的是泰勒级数展开方法。其基本的差分表达式主要有三种形式：一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等，其中前两种格式为一阶计算精度，后两种格式为二阶计算精度。通过对时间和空间这几种不同差分格式的组合，可以组合成不同的差分计算格式。 有限差分法的主要内容包括：如何根据问题的特点将定解区域作网格剖分；如何把原微分方程离散化为差分方程组以及如何解此代数方程组。此外为了保证计算过程的可行和计算结果的正确，还需从理论上分析差分方程组的性态，包括解的唯一性、存在性和差分格式的相容性、收敛性和稳定性。对于一个微分方程建立的各种差分格式，为了有实用意义，一个基本要求是它们能够任意逼近微分方程，这就是相容性要求。另外，一个差分格式是否有用，最终要看差分方程的精确解能否任意逼近微分方程的解，这就是收敛性的概念。此外，还有一个重要的概念必须考虑，即差分格式的稳定性。因为差分格式的计算过程是逐层推进的，在计算第n＋1层的近似值时要用到第n层的近似值，直到与初始值有关。前面各层若有舍入误差，必然影响到后面各层的值，如果误差的影响越来越大，以致差分格式的精确解的面貌完全被掩盖，这种格式是不稳定的，相反如果误差的传播是可以控制的，就认为格式是稳定的。只有在这种情形，差分格式在实际计算中的近似解才可能任意逼近差分方程的精确解。关于差分格式的构造一般有以下3种方法。最常用的方法是数值微分法，比如用差商代替微商等。另一方法叫积分插值法，因为在实际问题中得出的微分方程常常反映物理上的某种守恒原理，一般可以通过积分形式来表示。此外还可以用待定系数法构造一些精度较高的差