ANSYS岩土计算例子

1、实用文档ANSYS土工结构计算案例ANSYSCHINA办事处2020年 5 月 23日实用文档目录计算题目及计算要求说明 1题目一 4一、计算说明 4二、计算所用 ANSYS邓肯的 EB 模型说明 5三、计算有限元模型及计算结果 6题目二 7一、用三维有限元模型计算 7二、用三维有限元模型计算 8题目三 10一、计算说明 10二、计算有限元模型及边界条件 10三、强夯地基固结计算 10题目四 17一、计算说明 17二、计算几何模型和有限元模型 17三、计算结果 181、计算边界条件182、计算结果193、结论20实用文档计算题目及计算要求说明题目一：高桩码头桩基与岸坡相互作用的线性有限元和非线

2、性有线元 分析题目二：大圆筒结构、波浪与地基的相互作用分析（大圆筒作为重力 式码头结构，波浪为水平动荷载，门吊为竖向动荷载，地基为三层以 上地基包括自抛碎石垫层、粘土层、粉细砂层和岩层，粉细砂层可能 在波浪动荷载作用下液化造成圆筒倾覆） 题目三：（冲击荷载下）强夯地基固结有限元分析（提供固结方程或 固结方程处理方案， 孔隙水压力消散计算方案、 沉降计算方案及其他 一些处理技巧） 题目四：在降雨情况下土工格栅加筋土挡墙边坡上公路稳定分析 （由 上至下为公路面层，垫层，挡墙，挡墙面板采用预制混凝土块 0.6 0.6 0.6m3，混凝土后方为钩挂式土工格栅，边坡比较陡，边坡 有一定排水特性）。具体处

3、理方案包括：1、提供计算输入界面2、计算模型或采用本构情况3、前处理方案及网格划分技巧4、特殊材料或模型嵌入技术5、计算技巧及解决方案6、后处理提供容实用文档具体报价方案包括：进行简单报价， 涉及以上题目的各模块的综合报价 （包括前后处理及 解决以上问题的结构与岩土问题的模块报价 . 能用通用模块计算尽 量用通用模块 , 必须用 CivilFEM 模块计算的请注明 .提供解决方案时间： 2005年5 月18日之前。实用文档题目一高桩码头桩基与岸坡相互作用的线性有限元和非线性有线元分析一、计算说明高桩码头桩基与岸坡相互作用的线性有限元和非线性有线元分析， 可以用 空间有限元模型或平面有限元模型数

4、值求解。空间有限元模型如图 1-1 所示。考虑到岸坡主要是平面问题，为了节省单 元，可以只取四榀排架建立空间模型， 并以中间两个排架的受力情况作为整个码 头段的桩的受力标准。 横梁可用梁单元或体单元， 桩可用梁单元或壳单元或体单 元，地基土可用体单元，桩土之间可加面面接触单元 Targe170 和 contal174 以模拟桩土之间的相互作用。梁和桩作为线弹性材料，土可作为DP 材料按弹塑性分析或可用 EB 材料以作非线性弹性分析。图 1-1 高桩码头桩基与岸坡相互作用空间有限元模型平面有限元模型如图 1-2 所示。计算模型宽度在码头前沿以前和后方堆场 可各取码头宽的 2.5 倍，计算模型深度

5、可从码头面到桩端以下 510 米。横梁和 桩可用梁单元或平面单元，其模量和尺寸可按相似原理进行换算。不管空间还是平面模型，在计算模型围的原状土，自重引起的应力场应按 初始应力场考虑， ANSYS提供了初应力计算的功能，解决初应力问题非常方便。、计算所用 ANSYS邓肯的 EB 模型说明实用文档本算例是一个高桩码头桩基与岸坡相互作用的线性有限元和非线性有线元分析，是一个平面计算的例子。本例的土体用邓肯的 EB模型模拟。 EB模型为非线性弹性模型， 土的应力应变关系为双曲线， 如图 1-3 所示。切线变形模量为 EtEtKpapa2Rf ( 13)(1 sin )2ccos 2 3 sin11)1

6、2)13)切线泊松比按下式确定：Et3(1 2 )B Kbpa(pa3 )m式中，K 、n 试验常数； pa 大气压； B 体变模量； m 应力指数； 1 、 3 大小主应力； K b 体积模量系数；Rf 破坏比； c、 凝聚力和摩擦角。上述非线性弹性问题的有限元求解方法，可用增量法或迭代法。本例选用 迭代法求解。 ANSYS提供参数化设计语言 APDL，可以方便用于按参数建模， 也可 以用来扩展有限元的分析能力，创建各种控制方案。另外， ANSYS的命令流也可 以有效控制和实现很多求解过程和求解方案。 本例结合 ANSYS参数化设计语言及 ANSYS邓肯的 EB 模型来实现。实用文档三、计算

7、有限元模型及计算结果图 1-4 为计算有限元模型，全用采用 ANSYS自动四边形网格功能来划分， 从而保证了较高的计算精度。图 1-4 高桩码头桩基与岸坡相互作用计算元模型图 1-5 高桩码头桩基与岸坡相互作用变形图题目二大圆筒结构、波浪与地基的相互作用分析实用文档一、用三维有限元模型计算大圆筒在外荷载（包括动荷载）作用下和地基相互作用问题，可按空间问 题求解。现以大圆筒在振动下沉时的动力分析为例阐述如下。大圆筒求解有限元模型如图 2-1 所示。用六面体等参元模拟土基础，用曲 面壳体单元模拟钢圆筒， 用梁单元模拟钢圆筒柱壳上的纵横向加强肋， 在大圆筒 和基土之间设置接触面单元以模拟筒体和地基的

8、相互作用。对瞬时动荷载采用时程积分的方法，将一个周期的荷载分为足够多个时间 步（比如 300 步）进行求解。图 2-2 所示为大圆筒下沉到十六米时在 T11.5 秒时的竖向应力云图。图 2-1 大圆筒有限元模型图 2-2 大圆筒下沉到十六米时在 T 11.5 秒时的竖向应力云图大圆筒在侧向荷载（如波浪力）作用下的倾复问题也可用上述空间模型求 解。在侧向荷载作用下， 达到一定程度时， 圆筒和土体之间的接触单元会反映出 筒体和土体间的分离现象。二、用三维有限元模型计算参考日本规，大圆筒的位移和倾复问题，也可以用平面方法按如下模型求实用文档解。纵向取单宽按平面问题计算，横向宽度 B 可按纵向单元长度

9、面积相等的原 则确定。墙后土压力按假想平面和库仑极限主动土压力计算，侧向摩阻力由主动土 压力和土与墙背间摩擦角 计算确定。墙前水平土抗力通过设置非线性弹簧元经非线性迭代计算确定。 筒底地基反力通过设置筒底竖向弹簧元经计算确定； 筒底基土的剪切反力， 通过设置筒底水平向弹簧元计算确定， 根据上述假定，在 ANSYS上建立有限元模 型如图 2-3 。经非线性迭代计算，得到土抗力和地基反力如图 2-4 所示实用文档图 2-4位移图见图 2-5 ，最大位移 DMX10.7688cm。各关键点的位移如表图 2-5 位移云图表一 关键点的位移关键点竖向位移（ cm）水平位移（ cm）水平允许最大位移H1.

10、5%（cm）筒顶2.88 （向下）7.41 （向外）29.25 （满足要求）海底面2.84 （向下）5.04 （向外）筒底前趾2.42 （向下）0.179 （向）筒底后趾1.73 （向下）0.566 （向外）实用文档题目三（冲击荷载下）强夯地基固结有限元分析（提供固结方程或固结方程处理方案，孔隙水压力消散计算方案、沉降计算方案及其他一些处理技巧）一、计算说明本题目对强夯地基固结进行有限元分析， 可采用 ANSYS STRUCTUR分E析模 块来完成。二、计算有限元模型及边界条件图 3 1 有限元模型及模型尺寸图 3 3 冲击荷载作用下的变形图图 3 2 冲击荷载加载图三、强夯地基固结计算实用文

11、档图 34 冲击荷载作用下的等效应力云图根据 Skempton土力学计算公式， 每一个荷载增量导致的孔隙水压力的增量 可用式 (3-1) 表示：u B 3 A( 13) , B 1，A 1/3 (3-1)根据上式 (3-1) ，采用 ANSYS可计算得到动载作用下孔隙水压力增量。图 3-5 动荷载末时的孔隙水压力增量u 云图实用文档图 3 6 动荷载末时的孔隙水压力增量u 云图由式： u u0 u gh u（32）可计算得到孔隙水压力的分布，如下图所示。图 37 动荷载末时的孔隙水压力 u 云图孔隙水压力与水头的关系可用式（ 3 3）表示，图 38 为冲击荷载作用末 时水头分布云图。33）实用

12、文档图 3 8 冲击载荷末时水头 h实用文档图 3 10 1天后各点水头分布图 3 11 4 天后各点水头分布实用文档图 3 12 5 天后各点水头分布图 3 13 10天后各点水头分布图 3 14 20 天后各点水头分布 图 3 15 30 天后各点水头分布实用文档图 3 16 50 天后各点水头分布 图 3 17 60 天后各点水头分布实用文档图 3 18 65 天后各点水头分布图 3 19 100 天后各点水头分布为了得到在冲击载荷作用下孔隙水压力的消散过程， 如图 321图 323 所示分别画出了图 320所示 3个点的孔隙水压力 u随时间的变化曲线。图 3 20 孔隙水压力 u 消散

13、过程三个监控点位置实用文档4.60E+044.50E+044.40E+044.30E+044.20E+044.10E+044.00E+04-2.00E+06 0.00E+002.00E+064.00E+066.00E+068.00E+061.00E+07时间（秒）图 321 第 1 个监控点孔隙水压力 u随时间的消散过程时间（秒）)aP(力压水隙孔图 322 第 2 个监控点孔隙水压力 u随时间的消散过程)aP(力压水隙0.00E+ 1.00E+ 2.00E+ 3.00E+ 4.00E+ 5.00E+ 6.00E+ 7.00E+ 8.00E+ 9.00E+00 06 06 06 06 06 0

14、6 06 06 06时间（秒）图 323 第 3 个监控点孔隙水压力 u随时间的消散过程实用文档题目四：在降雨情况下土工格栅加筋土挡墙边坡上公路稳定分析（由上至下为公路 面层，垫层，挡墙，挡墙面板采用预制混凝土块 0.6 0.6 0.6m3，混凝土后方为 钩挂式土工格栅，边坡比较陡，边坡有一定排水特性） 。一、计算说明题目四是一渗流结构耦合计算问题。而 ANSYS具有功能强大的渗流结 构耦合计算功能。 下面是针对题目四， 采用 ANSYS所完成的： 在降雨情况下土工 格栅加筋土挡墙边坡上公路稳定分析。二、计算几何模型和有限元模型以下各图分别是计算模型及有关尺寸、材料的说明图 4-1 计算几何模型图 4-2 计算模型尺寸及有关材料实用文档图 4-3 渗流结构耦合计算有限元模型该坝体结构共有六种材料，上图 4-3 中分别用不同颜色表示结构不同的材 料。三、计算结果1、计算边界条件图 4-4 渗流结构耦合计算有限元模型实用文档以上是边坡路基在降雨情况下稳定分析边界条件说明。除了考虑降雨，同 时路基两侧考虑了 10m常水头的作用。2、计算结果雨情况下压力水头等值