数值模拟报告(DOC)

《数值模拟报告(DOC)》由会员分享，可在线阅读，更多相关《数值模拟报告(DOC)（23页珍藏版）》请在装配图网上搜索。

1、数值模拟报告(DOC)第一部分：数值模拟技术研究文献综述浅析数值模拟技术1. 引言近年来，随着我国大规模地进行西部大开发”和尚水北调”等巨型工程，越来越多的岩土工程难题摆在我们面前，单纯依靠经验、解析法显然已不能有效指导工程问题的解决，迫切需要更强有力的分析手段来进行这些问题的研究和分析。自R.W.Clough上世纪60年代末首次将有限元引入某土石坝的稳定性分析以来，数值模拟技术在岩土工程领域取得了巨大的进步，并成功解决了许多重大工程问题。特别是个人电脑的普及及计算性能的不断提高，使得分析人员在室内进行岩土工程数值模拟成为可能。在这样的背景下，数值模拟特别是三维数值模拟技术逐渐成为当前中国岩土

2、工程研究和设计的主流方法之一，也使得岩土工程数值模拟技术成为当今高校和科研院所岩土工程专业学生学习的一个热点。采用大型通用软件对岩土工程进行数值模拟计算，在目前已成为项目科研、工程设计、风险评估等岩土类项目的必须，学习和掌握Ansys、FLAC3D、UDEC等数值计算软件已成为学校、科研院所对工程从业人员的基本要求。数值模拟方法主要有限元法、边界元法、加权余量法、半解析元法、刚体元法、非连续变形分析法、离散元法、无界元法和流形元法等，各种方法都有其对应的软件。2. 数值模拟的发展趋势可以说，继理论分析和科学试验之后，数值模拟已成为科学技术发展的主要手段之一。随着软件技术和计算机技术的发展，目前

3、国际上数值模拟软件发展呈现出以下一些趋势：(1) .由二维扩展为三维。早期计算机的能力十分有限，受计算费用和计算机储存能力的限制，数值模拟程序大多是一维或二维的，只能计算垂直碰撞或球形爆炸等特定问题。随着第三代、第四代计算机的出现，才开始研制和发展更多的三维计算程序。现在，计算程序一般都由二维扩展到了三维，如LS-DYNA2D和LS-DYNA3D、AUTODYN2D和AUTO-DYN3D。(2) .从单纯的结构力学计算发展到求解许多物理场问题。数值模拟分析方法最早是从结构化矩阵分析发展而来，逐步推广到板、壳和实体等连续体固体力学分析，实践证明这是一种非常有效的数值模拟方法。近年来数值模拟方法已

4、发展到流体力学、温度场、电传导、磁场、渗流等求解计算，最近乂发展到求解几个交义学科的问题。例如内爆炸时，空气冲击波使墙、板、柱产生变形，而墙、板、柱的变形乂反过来影响到空气冲击波的传播，这就需要用固体力学和流体动力学的数值模拟结果交义迭代求解。.由求解线性问题进展到分析非线性问题。随着科学技术的发展，线性理论已经远远不能满足设计的要求。诸如岩石、土壤、混凝土等，仅靠线性计算理论就不足以解决遇到的问题，只有采用非线性数值算法才能解决。众所周知，非线性的数值模拟是很复杂的，它涉及到很多专门的数学问题和运算技巧，很难为一般工程技术人员所掌握。为此，近年来国外一些公司花费了大量的人力和资金，开发了诸如

5、LS-DYNA3D、ABAQUS和AU-TODYN等专长求解非线性问题的有限元分析软件，并广泛应用丁工程实践。这些软件的共同特点是具有高效的非线性求解器以及丰富和实用的非线性材料库。3. 数值模拟的基本原理一般而言，岩、土体处丁三向受力状态，其破坏模式往往表现为压-剪破坏和拉伸破坏。要分析和预测岩、土体在外力作用下的变形、破坏，就需要对其变形、破坏情况进行较为直观地再现。岩土工程数值模拟正是从岩、土体的受力状态出发，来分析和预测岩、土体破坏情况的一种手段。其基本原理是以典型试样的物理试验（室内试验或现场试验）获得的强度来表征整个地质体的岩、土体强度，以边界条件替代地质体周围所受的约束条件，借由

6、本构关系表达岩、土体在外力作用下的应力-应变特性，最终了解、预测岩、土体变形破坏情况。它具有鲜明的时代特征，以计算机为实现平台，是信息化时代的产物。通过与其它方法（如人工智能、人工生命科学、随机模拟、模糊数学、灰色理论以及分形理论等）交义共生、相互耦合嫁接，以获得更广阔的发展空间。从广义上来说，岩、土体的室内试验和原位试验也是一种模拟手段，本文称之为物理模拟。之所以如此称谓，是因为它们也是为较真实地近似再现岩、土体在其所赋存的环境中所处的受力状态所采用的一种手段。从这个意义上来说，它与数值模拟的基本原理是相同的，因此，可以将数值模拟称为虚拟实验室模拟。所不同的是，数值模拟除可以进行常规尺寸模型

7、的模拟外，还可以进行宏观和细观两个层面尺寸模型的模拟，而其输入的参数则需通过物理模拟来提供。因此，数值模拟是与物理模拟并行发展、相互补充和相互验证的试验系统。相较丁其它方法，数值模拟具有可重复和操作性强，费用低廉，不受模型尺寸控制，可视化程度高的优点，能有效延伸和扩展分析人员的认知范围，为分析人员洞悉岩、土体内部的破坏机理提供了强有力的可视化手段。当然作为一种分析方法，它也有自身的缺点，主要是易受制丁岩、土体结构的描述和模型概化的准确性及合理性；受制丁岩、土体物理试验模拟结果的准确性；受制丁岩、土体本构关系与实际岩、土体力学响应特性拟合程度的高低。4. 数值分析方法中存在的问题到目前为止，研究

8、计算工程的文章很多，但真正用丁实际工程的数值分析方法（例如有限元法等）却较少。部分原因在丁有较多不成功应用的实例。为什么会有这种情况，原因是多方面的，下面列出几条仅供参考：（1）对岩土工程数值分析方法缺乏系统的知识和深入的理解，出现问题时不知道在什么情况下届丁理论问题或数学模型问题；在什么情况下是届丁计算方法问题或本构模型问题；在什么情况下是参数的确定问题或计算本身的问题等。（2）各种本构模型固有的局限性。具有多相性土的物理力学性质太复杂，难以准确地用数学模型和本构模型描述。例如邓肯一张模型不能反映剪胀性，不能反映压缩与剪切的交义影响；模型只能考虑硬化，不能反映软化；模型不能反映各向异性。剑桥

9、模型也仅能考虑硬化而不能反映软化，不能反映土的剪切膨胀和各向异性，不能用丁超固结土等。(3) 现有的试验手段和设备不能提供适当、合理和精确的参数。靠少数样本点所获得的参数难以准确地描述整个空间场地的物理力学性能；土的参数因土样扰动难以高质量的获取，其精度很差。有些模型要求较多的参数，但这些参数用常规的试验手段和设备难以获取等。岩土工程中如何应用精确的数学模型和本构模型是一个值得注意的问题。在一般结构分析中，因材料的力学性质简单、均匀，不确定性较小，一般采用较精确的数学模型会得到较精确的分析结果。但就土这种材料而言，因其不确定性非常大，其情况发生了很大的变化。众所周知，场地土性及其参数勘察结果的

10、精度和准确性是很差的，由此导致既使采用了很精确的数学模型，但因输入参数的精度不能与之相匹配，其计算结果同样会很差。采用精确的数学模型还会给人造成一种错觉，让人觉得其计算结果也一定会更好、更可靠。这样可能使人们忽略了精确的数学公式也照样会有出错的可能性。只有当输入参数的质量和精度很高，并能与数学模型的精度相匹配时，才有可能得到较为准确的计算结果。5. 结语20世纪60年代以后，由丁电子计算机的飞速发展使岩土工程数值分析方法得到不断发展和完善，并用丁岩土工程实践。虽然在工程实际使用中数值分析方法存在一些问题，但只要认活问题的实质，并采取措施去解决它，相信随着岩土工程数值分析方法的不断发展及其工程经

11、验的不断积累，在工程实践中将会得到越来越多的应用，它必将成为岩土工程分析中的有力工具。1 参考文献张森，言志信，段建.边坡开挖数值模拟及其稳定性评价研究J.西部探矿工程.(3).汪军，刘海波.边坡稳定性的有限元数值模拟建模J.华北科技学院学报(0).陈印东，刘叔灼.基丁强度折减法的边坡稳定性分析J.科学技术与工程(0).2 王浩.类土质路堑高边坡典型失稳机制与加固工程对策的数值模拟研究D.铁道部科学研究院，2004.3 张超，杨春和.有限差分强度折减法求解边坡稳定性J.土木工程与管理学报.4 郑颖人，赵尚毅，宋雅坤.有限元强度折减法研究进展J.后勤工程学院学报.(0).5 邹宝祥，李明，唐伟华

12、.某大桥边坡稳定性FLAC3D数值模拟分析J.山西建筑.郭辉.山西晋城土质垂直高边坡稳定性计算及数值模拟研究D.西安科技大学,2011.6 郭志柳，陈建东，吴鹏.填土物理力学性质对路堤边坡稳定性影响的数值模拟J.江西理工大学学报.(9).第二部分：数值模拟技术FLAC3D上机报告flaCd数值模拟上机题计算模型分别如图1、2、3所示，边坡倾角分别为30、45、60,岩土体参数为：密度p=2500kg/m3,弹性模量E=1x108Pa,泊松比1=0.3,抗拉强度bt=0.8X106Pa,内聚力C=4.2x104Pa,摩擦角4=17试用FLAC3D软件建立单位厚度的计算模型，并进行网格剖分，参数赋

13、值，设定合理的边界条件，利用FLAC3D软件分别计算不同坡角情况下边坡的稳定性，并进行结果分析。附换算公式：1kN/m3=100kg/m3剪切弹性模量：g由土=3.846e72(1)体积弹性模量：k=8.33e73(12)图1倾角为30的边坡（单位：m）计算命令流如下：newgenzonebrickp0000p110000p2010p30040size50110genzonebrickp040040p1100040p240140p374.64060p4100140&p574.64160p6100060p7100160size30110fixxrangex-0.10.1fixxrangex99.

14、9100.1fixyfixzrangez0.1-0.1modelelaspropdensity2500bulk3e9shear1e9setgravity00-10solveinixdisp0ydisp0zdisp0inixvel0yvel0zvel0modelmohrpropdensity2500bulk8.33e7shear3.846e7c42000fric17ten800000solvefosfileslope30.savassociated计算结果如下：图1-a,网格剖分图图1-b,速度矢量图图1-c,速度等值线图图1-d,位移等值线图最终计算边坡的稳定性系数为：Fs=1.49分析：30

15、边坡稳定性系数采用的是FLACd内置的强度折减法求解，稳定性系数1.491,从稳定性系数系数可以判断该边坡处于安全状态。坡面最大速度为0.25mm/s,随着深度的增加，竖向应力逐渐增大。坡肩处出现下沉，最大值达到3.54m。图2倾角为45的边坡（单位：m）计算命令流如下：newgenzonebrip0000p110000p2020p30040size50110genzonebrip040040p1100040p240240p360060p4100240p560260p6100060&p7100160size30110modelmohrpropdensity2500.0bulk8.3e7shea

16、r3.8e7c42000.0tens0.8E6&friction17dilation20fixxyzrangez-0.10.1fixxrangex99.9100.1fixxrangex-0.10.1fixysetgravity=10.0plotaddaxeredplotcondis;定义循环终止条件defcalfosait1=0.02k11=0.6k12=2.0loopwhile(k12-k11)ait1fs=(k12+k11)/2.0refric=(atan(tan(17*pi/180)/fs)*180/pirecoh=42000/fs;折减实现过程commandinisxx0.0syy0.

17、0szz0.0sxy0.0sxz0.0syz0.0inixvel0.0yvel0.0zvel0.0inixdis0.0ydis0.0zdis0.0profricrefriccohrecohsetmechratio1e-5solvestep5000printfsend_commandaa=mech_ratioifaa1,从稳定性系数系数可以判断该边坡处于安全状态。坡面最大速度为0.64mm/s,随着深度的增加，竖向应力逐渐增大。坡肩处出现下沉，最大值达到5,34m。100图3倾角为60的边坡（单位：m）计算命令流如下：newgenzonebrickp0000pl10000p2010p30（40s

18、ize50110genzonebrickp040040p1100040p240140p351.55060p4100140p551.55160&p6100060p7100160size30110attachfacerangez39.940.1fixxrangex-0.10.1fixxrangex99.9100.1fixyfixzrangez0.1-0.1modelelaspropdensity2500bulk3e9shear1e9setgravity00-10solveinixdisp0ydisp0zdisp0inixvel0yvel0zvel0modelmohrpropdensity2500b

19、ulk8.33e7shear3.846e7c42000fric17ten800000solvefosfileslope60.savassociated计算结果如下：图3-a,网格剖分图图3-b,速度矢量图FHlpVHlfliMfeWMUNOWXi；BBWhiKMJBMIla:UOb-BDFB|IBM心0皿图3-c,速度等值线图图3-d,位移等值线图最终计算边坡的稳定性系数为：Fs=0.94分析：60边坡稳定性系数采用的是FLACd内置的强度折减法求解，稳定性系数0.941，从稳定性系数系数可以判断该边坡处于安全状态。坡面最大速度为0.51mm/s,随着深度的增加，竖向应力逐渐增大。坡肩处出现下

20、沉，最大值达到3.24m。第三部分：数值模拟技术研究应用实例分析基于FLAC3D在不同土体参数条件下的边坡稳定性研究摘要：边坡稳定性的影响因素很多，最直接的是边坡岩土体的性质、地下水、边坡坡脚等。对丁人工填土的路堤边坡而言，主要是填土的物理参数。本文介绍了数值模拟法的工程应用的历史及现状，数值模拟在边坡中的应用。然后利用flaCd软件，采用控制变量法，在保持其他土体物理参数不变的情况下，分别依次改变弹性模量、泊松比、密度、凝聚力、内摩擦角，计算边坡的稳定性，观察边坡稳定性系数变化，以及各点位移量的变化，从而得出土体物理参数变化对路堤边坡稳定性的影响，为路堤边坡填土土类提供依据和建议。关键词：边

21、坡、物理参数、FLAC3D1研究背景21世纪以来，中国处丁快速发展的阶段，国内基础建设蓬勃发展。随着水利工程、公路等基础设施建设的大力开展，尤其是我国西部大开发战略的实施，大量公路建设深入西部山区，严重地破坏了局部区域内地质环境的平衡，大量的工程活动对地质环境的改变日益加剧，导致了大量地质灾害的发生，所以边坡的稳定性研究显得越来越重要。边坡稳定性评价一直是边坡工程的一项主要内容，也是边坡工程设计和施工的基础。边坡稳定性计算理论和判别方法可靠与否，关系到工程的安全问题，一旦边坡失稳，不仅会给国家带来巨大的经济损失，而且会危及人民生命财产安全。2国内外研究状态目前边坡稳定性分析方法研究主要集中在与

22、计算机技术、岩土力学、数学模型结合上，也产生和发展了一些新方法，其中三维稳定性分析是研究热点之一，以有限元法为代表的数值分析法以及各种不确定性分析方法发展迅速，而传统的极限平衡法主要以改进为主。根据分析认为边坡稳定性分析研究在以后需要解决以下主要问题：继续完善和发展现有理论和方法，扬长避短；建立具有普遍意义的边坡失稳机理和稳定性评价方法；统一评价标准，增加各方法之间的对比性；建立多因素的综合评价方法；建立反映边坡各个时段稳定状态的全过程评价方法；建立符合边坡稳定分析的理论体系或组合理论体系；重视人类活动动态，加强人类活动与边坡稳定的相互作用研究。3FLACD模拟计算3.1基本模型以下面的简单的

23、模型为基础，计算不同物理参数下的稳定性系数，来进行本次研究，模型及形体参数见图1。基本物理参数见表1。2040./路堤、/42地基、80/图1基本模型及形体参数表1基本模型的物理参数土体物理参数参数弹性模量E/MPa泊松比密度P/(g,cm3)凝聚力C/kPa内摩擦4/（。）角抗拉强度bt/MPa剪切弹性模量/MPaG2（1体积弹性模量/MPaE-K)3(12)100.32.542170.838.583.3在flaCd中建立模型，并计算边坡稳定性；创建几何模型：newgenzonebrickp0000pl8000p2020p30020size40110genzonebrick&p040020p

24、l80020p240220p360040&p480220p560240p680040p780240&size20110;赋予材料模型届性modelelassetgravity0,0,-10propertybulk=8e9shear=5e9density=2500;施加边界约束fixxrangex-0.10.1fixxrangex79.980.1fixyfixzrangez-0.10.1；试算solve;设置重力场inixdis=0ydis=0zdis=0inixvel=0yvel=0zvel=0;赋予材料模型届性modelmohrpropertybulk=0.833e8shear=0.385e

25、8friction=17propertycohesion=0.42e5tension=0.8e6；fos计算solvefosfileslopefos.savassociatedrLAC33.0022：17:0$MwI4w263013Canted即烦mX4.000*4001X；0000，工Y：0.000Z：2(KHte+0O1Z：ODOOBist22228+002Msg-iAJig22500CoincurofShetrSirairRateMagac-0。网+30CGrmCalculathriJTlO：8e-OQB(gQOOOOhWO05DDO&*QIKtJSDDKfeW5000De-CClcI.

26、WKs-035iOODDe-OQSio1OCe-W515DDJ&41C5Ifi2.0000=352DDO3&-Ca5ic2.600*035i&QOOe-OOCn?9.的5弱-的5irtsat=5De-OKFoSFoSvalUBii109ItascaCcisultjngGroup,Inc.Mkin的叫liA.HNUSA图2基本模型的稳定性系数和剪应变增量图3.2填土物理力学性质对路堤稳定性的影响分析3.2.1填土弹性模量对路堤稳定性的影响在表1基本参数的情况下，改变路堤填土弹性模量E的大T，考虑弹性模量为10MPa、15MPa、20MPa、25MPa、30MPa等5种取值，采用FLAC3件分析路

27、堤边坡安全系数Fs随弹性模量E的变化情况及其弹性模量E对路堤沉降的影响情况，计算得到的填土不同弹性模量E下路堤边坡的值Fs见下表2。表2填土不同弹性模量E下路堤边坡的值Fs弹性模量E/MPa1015202530稳定性系数Fs1.091.091.091.091.09从表2可以看出安全系数为定值Fs=1.10,即路堤填土的弹性模量E对路堤的稳定性几乎没有影响.-ri5图3路堤表面在不同弹性模量E下的z方向位移情况图3为路堤表面在不同弹性模量E下的z方向位移(即路堤表面沉降)情况，从中可以看出：(1) 同一弹性模量E下，路堤表面各点的沉降量由路堤中心向边缘逐渐减小，并不成线性变化。(2) 路堤表面距

28、路堤边缘10m内沉降量随填土弹性模量E的增加而大幅增大。10m范围外，路堤表面沉降量缓慢增加，距路堤边缘位置13m位置沉降量最大，即此处为路堤滑动面贯通位置。图4路堤表面各点在不同弹性模量E下的z方向位移情况从图4当中可以得出结论：随着弹性模量的增大，各点沉降量逐渐增大，且增加速率也逐步增大。3.2.2填土泊松比对路堤稳定性的影响在表1基本参数的情况下，改变路堤填土泊松比H的大小，考虑泊松比为0.15、0.20、0.25、0.30、0.35种取值，采用FLAC3D软件分析路堤边坡随泊松比两勺变化情况，计算得到填土不同泊松比m下路堤边坡的Fs值见表3。表3填土不同泊松比H下路堤边坡的值Fs泊松比

29、0.150.200.250.300.35稳定性系数Fs1.091.091.091.091.09从表3可以看出安全系数为定值Fs=1.09,即路堤填土的泊松比V对路堤边坡的稳定性几乎没有影响.图5为出=0.3下的路堤边坡剪应变增量及其位移图。从中可以看出，边坡滑动位置贯穿丁坡底和左半幅路基表面1/4处的位置，边坡的最大水平位移在坡底位置。3.00N22$6MorNce2013RKathn:0000Y0000Z0q邮句1562Z50OCsnterX.4.W0心1Y:I.OOOetOOOLZDtMt+WI服四蜘ContourofShearStrainIncntmentMsghc=OOOOs-fQOO

30、Sa阐酎O.OOOOf4COO_W.0QQ5C2-rOOOOfrOQ；tfl4DOS-aO2S4C0CW02*EOC妇贮$000W0扫S.toK*-F仲B6J2MgcfeiFrapeS22:30.09MpnUov253013Center.RotationK:4OfMlX:0000Y:rDQOvtOODY：0000l浏(WXJ1IQW0DisiW由卜OWUgg1.56冲Z2KOCor.wcfSKsrSiiIncrcmenfMbbc-00001*000jradirrcaicuatfl&乾由比O.OOOGi+COCOOOOOetOOOts200e-(W22.0000002kl4.0ggjwo-otc

31、tutwwe-(26.WC0-00CKlECDO(B2iKODeKIfl1OKrt-OQII.QOODa-OOlki1.30006-0011,2QOte-OOiIfi1.：M3e-K1=i.Ot-OOiFaSraSva.j10Vclocir1Mulrnum：16454-004g秘lutciOoneuibraGroup.In.MinnsapO帆MNUSA图9c=42时路堤边坡的剪应变增量云图、安全系数和速度欠量图从图8和图9可以看出，随着粘聚力的增大，边坡的安全系数呈线性增大。分析原因主要是粘聚力的增大，边坡的抗剪强度也增大，使得边坡的抗滑能力也增加，从而使边坡趋丁稳定，即安全系数增大。图9还可

32、以明显看出塑性贯通区域，即路堤边坡的潜在滑动面，速度欠量图表明滑动面外侧区域各网格点的速度明显大丁其他区域，说明这一区域已经出现滑动，即该区域发生了破坏。3.2.4填土内摩擦角对路堤稳定性的影响在表1基本参数的情况下，改变路堤填土内摩擦角小的大小，考虑填土内摩擦角为13、15、17、19、21等5种取值，采用FLACd软件分析路堤边坡安全系数Fs随填土内摩擦角小的变化情况，绘制小-Fs关系曲线及其拟合公式如图10Fs内摩擦角内摩擦角图10巾-Fs关系曲线及其拟合公式从图10可以看出，随着粘聚力的增大，边坡的安全系数呈线性增大。对比图9和图10,尤其是其中的拟合公式，由图10的拟合公式中斜率大丁

33、图9中的拟合公式中的斜率，故知内摩擦角对稳定性系数的影响大丁凝聚力对稳定性系数的影响。4结论经过对影响路堤边坡稳定性的填土物理力学性质进行数值模拟分析，可以得出以下结论：(1) 路堤填土弹性模量和泊松比对路堤边坡的稳定性没有影响，安全系数为定值.但对路堤的沉降有一定程度的影响，沉降量随填土弹性模量、泊松比的增加而成非线性减小。路堤边坡安全系数与填土粘聚力、内摩擦角分别成线性关系，安全系数随着粘聚力、内摩擦角的增大而成线性增大，边坡趋丁稳定，通过分析c-Fs和小-Fs的关系，得出摩擦角比粘聚力对路堤边坡的稳定性影响较大。(2) 路堤边坡的安全系数随边坡的坡率、填土高度的增加而减小。考虑路堤填料的物理力学性质对路堤边坡稳定性的影响，对公路建设过程中路堤填料的选择具有一定的指导作用。