岩土工程数值分析问题的研究.doc

岩土工程数值分析问题的研究 　　【摘 要】首先，介绍了笔者对我国岩土工程数值分析现状的调查结果及分析，然后，分析了采用连续介质力学分析岩土工程问题的关键，并讨论分析了岩土本构理论发展现状，提出对岩土本构理论发展方向的思考，最后对数值分析在岩土工程分析中的地位作了分析。分析表明，岩土工程数值分析结果是岩土工程师在岩土工程分析过程中进行综合判断的重要依据之一；采用连续介质力学模型求解岩土工程问题的关键是如何建立岩土的工程实用本构方程；建立多个工程实用本构方程结合积累大量工程经验才能促使数值方法在岩土工程中由用于定性分析转变到定量分析。 　　【关键词】岩土工程；连续介质力学；数值分析；本构理论；本构模型 　　1 引 言 　　本文首先介绍笔者对我国岩土工程数值分析现状的调查结果，然后就岩土工程分析中的关键问题，如何发展岩土本构理论和数值分析在岩土工程分析中的地位这 3 个问题提出粗浅的看法。 　　2岩土工程分析中的关键问题 　　岩土工程分析中人们常常将用简化的物理模型去描述复杂的工程问题，再将其转化为数学问题并用数学方法求解。一个很典型的例子是，饱和软黏土地基大面积堆载作用下的沉降问题被简化为Terzaghi 一维固结物理模型，再转化为 Terzaghi 固结方程求解。采用连续介质力学模型求解工程问题一般包括 　　下述方程：①运动微分方程式（包括动力和静力分析两大类）；②几何方程（包括小应变分析和大应变分析两大类）；③本构方程（即力学本构方程）。对一具体工程问题，根据具体的边界条件和初始条件求解上述方程即可得到解答，对复杂的工程问题，一般需采用数值分析法求解。对不同的工程问题采用连续介质力学模型求解，所用的运动微分方程式和几何方程是相同的，不同的是本构方程、边界条件和初始条件。当材料为线性弹性体，本构方程为广义虎克定律。将岩土材料视为多相体，采用连续介质力学模型分析岩土工程问题一般包括下述方程[2]：①运动微分方程式（包括动力和静力分析两大类）；②总应力= 有效应力+ 孔隙压力（有效应力原理）；③连续方程（总体积变化为各相体积变化之和）；④几何方程，包括小应变分析和大应变分析两大类；⑤本构方程，即力学和渗流本构方程。 　　将多相体与单相体比较，基本方程多了 2 个，即有效应力原理和连续方程，且本构方程中多了渗流本构方程。对不同的岩土工程问题，基本方程中运动微分方程式、有效应力原理、连续方程和几何方程的表达式是相同的，不同的是本构方程。对一具体岩土工程问题，根据具体的边界条件和初始条件求解上述方程即可得到解答，一般需采用数值分析法求解。从上面分析可知，采用连续介质力学模型分析不同的岩土工程问题时，不同的是本构模型、边界条件和初始条件。对一个具体的岩土工程问题，边界条件和初始条件是容易确定的，而岩土的应力-应变关系十分复杂，采用的本构模型及参数对计算结果影响极大。 　　采用连续介质力学模型分析岩土工程问题一般需采用数值分析法求解，有限单元法对各种边界条件和初始条件，采用的各类本构方程都有较大的适应性。土的应力-应变关系十分复杂，自 Roscoe 和他的学生建立剑桥模型至今已近半个世纪，理论上已提出数百个本构方程，但得到工程应用认可的极少，或者说还没有。从这个角度讲，采用连续介质力学模型求解岩土工程问题的关键问题是如何建立岩土材料的工程实用本构方程。 　　3 如何发展岩土本构理论的思考 　　Janbu 认为，反映作用与效应之间的关系称为本构关系，力学中的虎克定律、电学中的欧姆定律、渗流学中的达西定律等反映的都是最简单的本构关系。岩土是自然、历史的产物，具有下述特性：土体性质区域性强，即使同一场地同一层土，沿深度和水平方向变化也很复杂；岩土体中的初始应力场复杂且难以测定；土是多相体，一般由固相、液相和气相三相组成，土体中的三相有时很难区分，而且处不同状态时，土的三相之间可以相互转化。土中水的状态又十分复杂；土体具有结构性，与土的矿物成分、形成历史、应力历史和环境条件等因素有关，十分复杂；土的强度、变形和渗透特性测定困难。岩土的应力-应变关系与应力路径、加荷速率、应力水平、成分、结构、状态等有关，土还具有剪胀性、各向异性等，因此，岩土体的本构关系十分复杂。至今人们建立的土体的本构模型类别有弹性模型、刚塑性模型、非线性弹性模型、弹塑性模型、黏弹性模型、黏弹塑性模型、边界面模型、内时模型、多重屈服面模型、损伤模型、结构性模型等等。已建立的本构模型多达数百个，但得到工程师认可的极少，或者说还没有。怎么走出困境？这是我们必须面对的难题。笔者认为，对土体本构模型研究应分为两大类，科学型模型的研究和工程实用性模型的研究。科学型模型重在揭示、反映某些特殊规律，如土的剪胀性、主应力轴旋转的影响等。该类模型也不能求全面，一个模型能反映一个