数理方程在岩土工程中的应用.docVIP

数理方程在岩土工程中的应用 岩土工程是以实践和试验为基础的工程性学科，但近年来正在发展中的计算土力学，为岩土工程的发展和应用工程实践提供了便捷通道。现实中的工程问题是不能或很难用工程试验的方法来究的，怎样在试验前作较准确的预测，由于理论的发展远滞后于工程实践的应用需要，人们必须寻求新的路径：既能满足实践的定量需要，又尽可能的符合理论的定性要求。因此，发展出多种偏微分问题的处理方法，《数学物理方程与特殊函数》作为一门工具性的基础学科在计算土力学中显得尤为重要，在处理一些实际课题时，电子计算机已越来越成为一个重要的工具，要能有效地将数学物理方程遇特殊函数同电子计算机来解决实际工程问题，其先决条件是：(1)建立合理的数学物理模型。对决定岩土性质的重要变量及参数，通过大量偏微分方程及数学模型来描述；比较及优化各种模型，选定能符合实际工程的模型。(2)确定合理的数学物理方程的边界条件与初始条件。实际中的边界条件往往是复杂多变的，初始条件更是无法精确地确定，所以就存在“抓主忽次”的问题（即能真实地反映问题，又能简化方程，更能方便计算）。(3) 对相应的偏微分方程进行定性的研究。许多偏微分方程的非级数解的存在与否仍然备受争议，我们只要确定存在性、稳定性、适应性才能进行下一步的研究分析。 寻求或选择有效的求解方法，特别是数值的求解方法（即设幂级数为微分方程的解，确定系数即可，取满足精度要求的有限项进行计算）。 总结几点微分方程在岩土中的应用： 一、太沙基一维渗透固结微分方程 达西定律： 压缩定律： 有效应力原理： 单元体流入的水量： 流出水量： 水量的变化为： 单元土体由于压缩，空隙体积的改变而发生变化 土粒体积保持不变 则 压缩定律： 用体积压缩系数表示为： 土体压缩系数表示为 土体压缩的过程就是孔隙水压力和有效应力转化过程： 孔隙水压力 则 将和代入 得 令 （竖向固结系数） 这就是饱和土的一维渗透固结微分方程，也称为一维热传导公式。 一维渗透固结微分方程式的解答 根据一定的初始条件和边界条件，可以求解一维渗透固结微分方程，求得任一深度z在任一时间的孔隙水压力u的表达式。 初始条件和边界条件如下： 根据以上初始条件和边界条件采用分离变量法，可求得满足初始条件和边界条件的解答如下： 定解问题 按照分离变量法，令 解得： 太沙基一维渗透固结微分方程就是经过应用土力学中的基本知识推导出数理方程中的一维热传导方程，然后再根据边界条件用分离变量法进行求解，这样可以求解土的固结。 二、二维渗流的基本方程 取一微元土体，面积为，厚度为 单位时间流入这个微元体的水量为 单位时间流出这个微元体的水量为 对于饱和土流入微单元的水量必须等于流出的水量，即： 从而化简得到 根据达西定律，得： 和 代入 对于各向同性的均质土， 则最终化简为： 这就是饱和各向同性的均质土中二维渗流的基本微分方程。也是著名的拉普拉斯公式，通过求解一定边界条件下的拉普拉斯公式，即可求得该条件下的渗流场。 三、数理方程在土壤环境中的应用 随着城市经济的快速发展发展和人口的不断增加，人类活动对城市土壤质量的影响日显突出。对城市土壤地质异常的考证、评价，研究人类活动影响下城市地质环境的演变模式，日益成为人们关注的焦点。刘哲民应用单因子指数和尼梅罗综合污染指数对宝鸡土壤的重金属污染进行了评价。庞资估等也对复杂地形下瞬时污染源扩散机理进行了模拟与研究。对城市土壤重金属污染问题，首先结合单因子指数定义出归一化加权污染指数，对城市内不同区域的重金属污染进行分析，然后根据重金属元素的空间分布忽略一些次要因素，建立二维偏微分方程模型模拟污染物的扩散情况，模型用有限差分法迭代求出拉普拉斯数值解，实现对污染源位置的估计。 可得以下拉普拉斯方程： 在给定的一个区域里，假设最初的污染源集中于一点，因此可以给如下初值条件： ? 其中c（x,y,z）为某一点重金属元素的浓度，为重金属元素的传播系数，为降解系数。求解得： ??? 通过对各功能区污染程度分析，得出污染集中分布在海拔较低的人类活动密集区，同时也不考虑土壤类型对扩散影响，建上了重金属污染物在各向同性二维传播的偏微分方程模型，利用数理方程的算法，实现污染源点坐标估计。实验数据验证了所提出模型的可行性和有效性。 通过分析汞在土壤中的迁移规律可得汞在城市土壤表皮扩散主要受各方向上都存在着浓度梯度和水流的影响．由此本文以浓度梯度结合水流影响，首先建立偏微分方程；然后对数据进行回归得到扩散方程，从而确定污染源中心。利用偏微分方程建立了一种确定土壤表皮汞污染源的三维扩散偏微分方程模型．通过对数据进行回归得到扩散方程，从而确定污染源中心。这种利用数理方程的计算方法不但使污