数值模拟预案.doc

第一部分：数值模拟技术研究文献综述 浅析数值模拟技术 1．引言 近年来，随着国大规模地进行“西部大开发”和“南水北调”等巨型工程，越来越多的岩土工程难题摆在我们面前，单纯依靠经验显然已不能有效指导工程问题的解决，迫切需要更强有力的分析手段来进行这些问题的研究和分析。自R.W. Clough 上世纪60年代末首次将有限元引入某土石坝的稳定性分析以来，数值模拟技术在岩土工程领域取得了巨大的进步，并成功解决了许多重大工程问题。特别是个人电脑的普及及计算性能的不断提高，使得分析人员在室内进行岩土工程数值模拟成为可能。在这样的背景下，数值模拟特别是三维数值模拟技术逐渐成为当前中国岩土工程研究和设计的主流方法之一，也使得岩土工程数值模拟技术成为当今高校和科研院所岩土工程专业学生学习的一个热点。 采用大型通用软件对工程进行数值模拟计算，在目前已成为项目科研、工程设计、风险评估等岩土类项目的必须，学习和掌握nsys、、等数值计算软件已成为学校、科研院所对从业人员的基本要求。数值模拟的基本原理 一般而言，岩、土体处于三向受力状态，其破坏模式往往表现为压-剪破坏和拉伸破坏。要分析和预测岩、土体在外力作用下的变形、破坏，就需要对其变形、破坏情况进行较为直观地再现。岩土工程数值模拟正是从岩、土体的受力状态出发，来分析和预测岩、土体破坏情况的一种手段。 其基本原理是以典型试样的物理试验（室内试验或现场试验）获得的强度来表征整个地质体的岩、土体强度，以边界条件替代地质体周围所受的约束条件，借由本构关系表达岩、土体在外力作用下的应力-应变特性，最终了解、预测岩、土体变形破坏情况。它具有鲜明的时代特征，以计算机为实现平台，是信息化时代的产物。通过与其它方法（如人工智能、人工生命科学、随机模拟、模糊数学、灰色理论以及分形理论等）交叉共生、相互耦合嫁接，以获得更广阔的发展空间。 从广义上来说，岩、土体的室内试验和原位试验也是一种模拟手段，本文称之为物理模拟。之所以如此称谓，是因为它们也是为较真实地近似再现岩、土体在其所赋存的环境中所处的受力状态所采用的一种手段。从这个意义上来说，它与数值模拟的基本原理是相同的，因此，可以将数值模拟称为虚拟实验室模拟。所不同的是，数值模拟除可以进行常规尺寸模型的模拟外，还可以进行宏观和细观两个层面尺寸模型的模拟，而其输入的参数则需通过物理模拟来提供。因此，数值模拟是与物理模拟并行发展、相互补充和相互验证的试验系统。 相较于其它方法，数值模拟具有可重复和操作性强，费用低廉，不受模型尺寸控制，可视化程度高的优点，能有效延伸和扩展分析人员的认知范围，为分析人员洞悉岩、土体内部的破坏机理提供了强有力的可视化手段。 当然作为一种分析方法，它也有自身的缺点，主要是易受制于岩、土体结构的描述和模型概化的准确性及合理性；受制于岩、土体物理试验模拟结果的准确性；受制于岩、土体本构关系与实际岩、土体力学响应特性拟合程度的高低。第二部分：数值模拟技术FLAC3D上机报告 FLAC3D数值模拟上机题 计算模型分别如图1、2、3所示，边坡倾角分别为30°、45°、60°，岩土体参数为： 密度ρ＝2500 kg/m3, 弹性模量E＝1×108 Pa，泊松比μ＝0.3， 抗拉强度σt＝0.8×106 Pa，内聚力C＝4.2×104 Pa，摩擦角φ＝17° 试用FLAC3D软件建立单位厚度的计算模型，并进行网格剖分，参数赋值，设定合理的边界条件，利用FLAC3D软件分别计算不同坡角情况下边坡的稳定性，并进行结果分析。 附 换算公式： 1 kN/m3= 100 kg/m3 剪切弹性模量：=3.846e7 体积弹性模量：=8.33e7 图1 倾角为30°的边坡(单位：m) 计算命令流如下： new gen zone brick p0 0 0 0 p1 100 0 0 p2 0 1 0 p3 0 0 40 size 50 1 10 gen zone brick p0 40 0 40 p1 100 0 40 p2 40 1 40 p3 74.64 0 60 p4 100 1 40 p5 74.64 1 60 p6 100 0 60 p7 100 1 60 size 30 1 10 fix x range x -0.1 0.1 fix x range x 99.9 100.1 fix y fix z range z 0.1 -0.1 model elas prop density 2500 bulk 3e9 shear 1e9 set gravity 0 0 -10 solve ini xdisp 0 ydisp 0 zdisp 0 ini xvel 0 yvel 0 zvel 0 model mohr prop density 2500 bulk 8.33e7