为实现岩石试样蠕变全过程的准确模拟,并从细观角度探究蠕变过程中微裂隙的发生和发展规律,在二维颗粒流程序(PFC2D)中开发出具有黏弹塑性特征的西原体流变接触本构模型,进一步提出包含两种非定常元件的非定常西原体模型,推导了模型本构关系和蠕变方程.在PFC2D中调用自定义西原体流变模型,通过参数调试,获得与真实试样具有相同强度特性的数值试样.以室内单轴压缩蠕变试验数据为基础,在Matlab中对模型非定常参数进行拟合反演分析.在此基础上,进行单轴压缩蠕变试验的模拟,计算过程中分别采用定常和非定常两种模型,并对微裂隙进行监测.对比分析结果表明:定常模型仅适用于衰减和稳定蠕变阶段;非定常模型也可用于描述加速蠕变阶段,从而准确模拟蠕变全过程;加速蠕变阶段主要是由微裂隙的加速发展而产生,加速蠕变将导致试样剪切破坏

4.5数值微分 先看一个实例: 已知20世纪美国人口的统计数据为(单位百万) 年份

9-1数值计算的基本思想 9-2流动与传热的数值计算 9-3 Saints22D软件简介

4.1 Newton-Cotes-公式 4.2复合求积法 4.3 Romberg算法 4.4* Gauss求积法 4.5数值微分

本研究应用Spalding数值计算方法,移植顶吹转炉成功的数模与计算软件,加入电磁场洛仑兹力,在理论分析的基础上研制了30t底电极直流电弧炉熔池电磁流场温度场数学模型,开发了计算软件,计算了流场、温度场

前言 二十一世纪是科学技术高速发展的信息时代,而计算机是处理 信息的主要工具,因此,人们已经认识到,计算机知识已成为 人类当代文化的一个重要组成部分。 计算机科学技术以惊人的速度向前发展,它的广泛应用已从传 统的数值计算领域发展到各种非数值计算领域。在非数值计算 领域里,数据处理的对象已从简单的数值发展到一般的符号, 进而发展到具有一定结构的数据。在这里,面临的主要问题是 :针对每一种新的应用领域的处理对象,如何选择合适的数据 表示(构构),如何有效地组织计算机存贮,并在此基础上又 如何有效地实现对象之间的\运算关系

综合前几节的内容我们知道 梯形公式,Simpson公式, Cotes公式的代数精度分别为 1次,3次和5次 复合梯形、复合 Simpson、复合 Cotes公式的收敛阶分别为 2阶、4阶和6阶 无论从代数精度还是收敛速度,复合梯形公式都是较差的 有没有办法改善梯形公式呢?

3.7数据拟合(最小二乘法) 实例:考察某种纤维的强度与其拉伸倍数的关系,下表 是实际测定的24个纤维样品的强度与相应的拉伸倍数 是记录:

当积分区间[a,b]的长度较大,而节点个数n+1固定时 直接使用 Newton-Cotes-公式的余项将会较大 而如果增加节点个数,即n+1增加时 公式的舍入误差又很难得到控制 为了提高公式的精度,又使算法简单易行往往使用复合方法 即将积分区间[a,b]分成若干个子区间 然后在每个小区间上使用低阶 Newton-Cotes-公式 最后将每个小区间上的积分的近似值相加

§5.3常微分方程组和高阶 微分方程的数值解法简介 一、常微分方程组的数值解法 下列包含多个一阶常微分方程的初值问题