《岩土工程数值分析学习笔记》由会员分享，可在线阅读，更多相关《岩土工程数值分析学习笔记（20页珍藏版）》请在金锄头文库上搜索。
1、岩土工程数值分析读书笔记学号：13102081401003 姓名：贾红晶岩土工程数值分析读书笔记摘要：阅读笔记分为两部分：理论学习和plaxis模拟相关问题。理论部分0岩土工程数值分析简介岩土工程问题解析分析是以弹塑性力学理论和结构力学作为理论依据，适用于解决连续介质、各向同性材料、未知量少、边界条件简单的工程问题，存在很大的局限性。岩土工程问题数值分析是借助于计算机的计算能力，适用于解决材料复杂、边界条件复杂、任意荷载、任意几何形状，适用范围广。岩土工程数值分析发展过程：20世纪40年代，使用差分法解决了土工中的渗流及固结问题，如土坝渗流及浸润线的求法、土坝及地基的固结等。20世纪60年代，
2、使用有限元法成解决了土石坝的静力问题的求解。20世纪70年代，使用有限元法解决了土石坝及高楼(包括地基)的抗震分析。20世纪80年代，边界元法异军突起，解决了半无限域的边界问题；地基的静力及动力问题都使用边界元法得到了有效地解决。岩土工程数值分析的方法有两类，一类方法是将土视为连续介质，随后又将其离散化，如有限单元法、有限差分法、边界单元法、有限元线法、无单元法以及各种方法的耦合。另一类计算方法是考虑岩土材料本身的不连续性，如裂缝及不同材料间界面的界面模型和界面单元的使用，离散元法，不连续变形分析，流形元法，颗粒流等数值计算方法。1数值分析过程中存在的问题及解决措施问题：(1)对岩土工程数值分
3、析方法缺乏系统的知识和深入的理解，出现问题时不知道在什么情况下属于理论问题或数学模型问题；在什么情况下是属于计算方法问题或本构模型问题；在什么情况下是参数的确定问题或计算本身的问题等。(2)各种本构模型固有的局限性。具有多相性土的物理力学性质太复杂，难以准确地用数学模型和本构模型描述。例如邓肯一张模型不能反映剪胀性，不能反映压缩与剪切的交叉影响；(3)现有的试验手段和设备不能提供适当、合理和精确的参数。靠少数样本点所获得的参数难以准确地描述整个空间场地的物理力学性能；土的参数因土样扰动难以高质量的获取，其精度很差。(4)数学模型还会给人造成一种错觉，让人觉得其计算结果也一定会更好、更可靠。这样
4、可能使人们忽略了精确的数学公式也照样会有出错的可能性。只有当输入参数的质量和精度很高，并能与数学模型的精度相匹配时，才有可能得到较为准确的计算结果。措施：(1)加强对土的本构模型的教学与培训，了解和掌握各种土的本构模型的优点和局限性以及模型参数的离散性。(2)在使用数值分析方法的同时，不断地积累使用经验，包括他人的经验。(3)在丰富的使用经验基础上，建立相应的使用规范。2 有限元法软件plaxis理论学习（1）模型的选用和分析：Mohr-Coulomb模型（MC）：这种模型被推荐用于问题的初步分析，对于每个土层，可以估计出一个平均刚度常数，由于这个刚度是常数，计算往往会相对较快，可以得到变形的
5、一个初步印象。五个输入参数，即：表示土体弹性的E和，表示土体塑性的和c，以及剪胀角。节理岩石模型（JR）：节理模型是一种各向异性的弹塑性模型，特别适用于模拟包括层理尤其是断层方向在内的岩层行为等。两个参数和c。Hardening-Soil模型（HS）：是一种改进了的模拟岩土行为的模型，适合所有类型的土；不能用来模拟滞后或者反复循环加载情形。三个参数三轴加载刚度E50、三轴卸载刚度Eur和固结仪加载刚度Eoed。我们一般取Eur=3E50和Eoed=E50作为不同土体类型的平均值。软土蠕变模型（SSC）：软土蠕变模型是一个新近开发的应用于地基和路基等的沉陷问题的模型，通常会过高地预计弹性岩土的行
6、为范围。软土模型（SS）：软土模型是一种Cam-Clay类型的模型，特别适用于接近正常固结的粘性土的主压缩，可以被Hardening-Soil模型取代。改进的Cam-Clay模型（MCC）：主要用于模拟接近正常固结的粘性土，在实际应用中是不被推荐的。如果要对所考虑的问题进行一个简单迅速的初步分析，我们建议使用Mohr-Coulomb模型。当缺乏好的土工数据时，进一步的高级分析是没有用的。在许多情况下，当你拥有主导土层的好的数据时，可以利用Hardening-Soil模型来进行一个额外的分析。毫无疑问，同时拥有三轴试验和固结仪试验结果的可能性是很小的。但是,原位实验数据的修正值对高质量实验数据来
7、说是一个有益的补充最后，软土蠕变模型可以用于分析蠕变（即：极软土的次压缩）。用不同的土工模型来分析同一个岩土问题显得代价过高，但是它们往往是值得的。首先，用Mohr-Coulomb模型来分析是相对较快而且简单的；其次，这一过程通常会减小计算结果的误差。（2）各模型的参数：MOHR-COULOMB模型的基本参数：E：杨氏模量；v：泊桑比；：内摩擦角；c：内聚力；：剪胀角，这些参数可以从土样的基本试验得到。杨氏模量E:在土力学中，初始斜率用E0表示，50%强度处的割线模量由E50表示。对于具有大范围线弹性行为的材料来说，使E0是符合实际的，但是对于土体加载问题一般使用E50。如果考虑隧道和开挖问题
8、中的卸载问题要用Eur替换E50.Tu 图2 强度参数设置Tu 图1 项目设置泊桑比V:在许多情况下得到的v值是介于0.3和0.4之间的。一般地说，除了一维压缩，这个范围的值还可以用在加载条件下，在卸载条件下，使用0.15和0.25之间的值更为普遍。内摩擦角:PLAXIS可以处理无粘性砂土c=0，但是有一些选项的执行不太好。为了避免复杂性，我们建议不熟练的用户至少输入一个较小值（使用c0.2kPa）。剪胀角:以度的方式指定的。除了严重的超固结土层以外，粘性土通常没有什么剪胀性（=0）。砂土的剪胀性依赖于密度和摩擦角。对于石英砂土来说，=30，的值比的值小30度，然而，剪胀角在多数情况下为零。的
9、小的负值仅仅对极松的砂土是实际的。HARDENING-SOIL模型的参数：Tu 图2 HS模型强度参数设置软土模型参数：Tu 图3 软土模型强度参数设置（3）单位和符号规定符号规定：在平面应变分析里，zz指向平面外。轴对称分析里，x代表径向坐标，y代表轴向坐标，z代表切向。此时，xx表示径向应力，zz表示环向应力。在所有输出数据里，压应力（包括孔隙压力）和压力设为负值，而拉应力和拉力设为正值。家 建模过程：绘制几何轮廓线 设定边界条件添加荷载添加材料特性划分网格设置初始条件计算工序分层计算 警 警告：如果修改已经建好的模型，要意识到必须重新生成网格，有时还包括重新生成初始条件，以便和修改后的模
10、型一致。Tu 图4 基本界面（4）荷载问题:荷载子菜单包含的选项，用于在几何模型里引入分布荷载、线荷载或集中荷载以及指定位移。荷载和指定位移既可用于模型边界，也可用于模型内部。指定位移是为了控制某些点的位移而强加于模型的特殊条件。大小为零的指定位移即为约束。注意：如果在一条几何线上同时施加指定位移和荷载，那么在计算过程中指定位移优先于荷载，在全约束线上施加指定位移，约束则优先于指定位移。因而，在全约束线段上没有必要施加指定位移、集中荷载、分布荷载。如果只有一个位移方向是预加的，而其他方向是自由的，那么就可以在自由位移的方向上施加集中荷载、分布荷载。（5）材料性状种类材料类型排水性状：该设置不产
11、生超静水压。很明显，它适用于干土，以及由于高渗透性（砂土）和/或缓慢加载的完全排水条件。该选项也可以用来模拟长期岩土性质，其中不需要模拟不排水加载和固结的精确历史。不排水性状：该设置用于研究超静水压的完全发展过程。由于低透水性（粘土）和/或快速加载，孔隙水渗流有时可以忽略。注意，要输入有效的模型参数。非多孔性状：使用该设置的类组，不论是初始孔压还是超静水压，都不予考虑。有关应用可以在模拟混凝土或结构性状的例子里找到。非多孔性状通常和线弹性模型组合应用。输入的饱和容重和渗透性参数对非多孔材料是无效的。（6）网格生成网格的基本单元类型为15节点的三角形单元或6节点的三角形单元。全局疏密度水平：很粗
12、疏、粗疏、中等、细密和很细密。默认情况下，全局疏密度水平为粗疏。局部疏密度：在可能出现强烈的应力集中或大变形梯度的区域，最好是用更精确（细密）的有限元元网格来模拟，这通常是在几何模型里出现边缘、角或结构对象的时候用到。这时候的局部单元尺寸因子就应当等于0.5。关于网格生成的几点建议：为了优化有限元计算，可以先用相对粗疏的网格作初始分析。这样可以检查几何模型的大小是否合适，还可以预计发生应力集中和大变形梯度的位置。所得成果主要用来加密有限元网格，生成加密的有限元模型。为了能够顺利地生成详细的有限元网格，首先应当从网格子菜单里选择需要的全局疏密度。此外，如果还要进行局部加密，那么应当首先加密类组，
13、然后加密几何线，最后再加密几何点。需要时，可以直接给几何点设局部单元尺寸因子。（7）初始条件在几何模型和有限元网格都建好之后，必须明确初始应力状态和初始构造。初始条件由两个不同的模式组成：一个模式是生成初始水压（水力条件模式），另一个模式是定义初始几何构造和生成初始有效应力场（几何构造模式）。注意：从初始条件返回到几何图形创建模式是允许的，但是一般不这么做，因为这样做可能会丢失一些有关初始条件的信息。（8）计算计算分类：分三种基本计算类型，塑性计算、固结分析和Phi-c折减（安全分析）。塑性计算：在弹-塑性变形分析当中，如果不需要考虑超静水压随时间的减小，那么应当选塑性计算。一般性塑性计算根据
14、未变形的几何图形得出刚度矩阵。这类计算适用于大多数岩土工程的实际应用情况。固结分析：要分析饱和粘性土的超静水压随时间的变化与消散，应该选择固结计算。一般情况下，执行完不排水塑性计算之后，再执行不添加另外荷载的固结分析。也可以在固结分析过程中施加荷载。但是，在接近破坏时要引起注意，因为此时迭代过程可能会不收敛。Phi-c折减（安全分析）：可以用减小土的强度参数的方法来执行安全分析。这个过程就叫做Phi-c折减，是一个单独的计算类型。如果要计算全局安全系数，就应当选择Phi-c折减。建议在计算工序列表的末尾再定义所有安全分析，并且用起始工序号参数注明是在给哪个计算工序求算安全系数。plaxis例题模拟算例一砂土层上圆形基础的沉降(选用柔性基础)图1表示放置在4m厚砂土层上半径为1m的一个圆形基础。砂土层下是深厚的坚硬岩石层。计算土体在上部荷载作用下产生的位移和应力。选用轴对称模型，竖直方向取4m(到岩石层顶部)，水平方向取5m(5倍的基础半径)。计算过程如下： 图5 新建项目参数设置 20 / 20 图6 几何模型 图7 基础模型参数设置 图8 地基模型的初始应力场 图9 计算参数设置放大 图10 定义计算位置图11 变形网格图12 变形总位移场
