ANSYS岩土计算例子

ANSYSReportANSYSt工结构计算案例ANSY&CHINA办事处2020年5月23日目录TOC\o"1-5"\h\z计算题目及计算要求说明 1\o"CurrentDocument"题目一 4\o"CurrentDocument"一、 计算说明 4\o"CurrentDocument"二、 计算所用ANSYSB肯一的E—B模型说明 5\o"CurrentDocument"三、 计算有限元模型及计算结果 6\o"CurrentDocument"题目二 7\o"CurrentDocument"一、 用三维有限元模型计算 7\o"CurrentDocument"二、 用三维有限元模型计算 8题目三 10\o"CurrentDocument"一、 计算说明 10\o"CurrentDocument"二、 计算有限元模型及边界条件 10\o"CurrentDocument"三、 强夯地基固结计算 10\o"CurrentDocument"题目四 17\o"CurrentDocument"一、 计算说明 17\o"CurrentDocument"二、 计算几何模型和有限元模型 17\o"CurrentDocument"三、 计算结果 181、 计算边界条件 182、 计算结果 193、 结论 20计算题目及计算要求说明题目一：高桩码头桩基与岸坡相互作用的线性有限元和非线性有线元分析题目二：大圆筒结构、波浪与地基的相互作用分析（大圆筒作为重力式码头结构，波浪为水平动荷载，门吊为竖向动荷载，地基为三层以上地基包括自抛碎石垫层、粘土层、粉细砂层和岩层，粉细砂层可能在波浪动荷载作用下液化造成圆筒倾覆）题目三：（冲击荷载下）强夯地基固结有限元分析（提供固结方程或固结方程处理方案，孔隙水压力消散计算方案、沉降计算方案及其他一些处理技巧）题目四：在降雨情况下土工格栅加筋土挡墙边坡上公路稳定分析 （由上至下为公路面层，垫层，挡墙，挡墙面板采用预制混凝土块0.60.60.6m3，混凝土后方为钩挂式土工格栅，边坡比较陡，边坡有一定排水特性）。具体处理方案包括：1、 提供计算输入界面2、 计算模型或采用本构情况3、 前处理方案及网格划分技巧4、 特殊材料或模型嵌入技术5、 计算技巧及解决方案6、 后处理提供容具体报价方案包括：进行简单报价，涉及以上题目的各模块的综合报价（包括前后处理及解决以上问题的结构与岩土问题的模块报价.能用通用模块计算尽量用通用模块，必须用CivilFEM模块计算的请注明.提供解决方案时间：2005年5月18日之前。题目一高桩码头桩基与岸坡相互作用的线性有限元和非线性有线元分析一、计算说明高桩码头桩基与岸坡相互作用的线性有限元和非线性有线元分析， 可以用空间有限元模型或平面有限元模型数值求解。空间有限元模型如图1-1所示。考虑到岸坡主要是平面问题，为了节省单元，可以只取四榀排架建立空间模型，并以中间两个排架的受力情况作为整个码头段的桩的受力标准。横梁可用梁单元或体单元，桩可用梁单元或壳单元或体单元，地基土可用体单元，桩土之间可加面—面接触单元 Targe170和contal174以模拟桩土之间的相互作用。梁和桩作为线弹性材料，土可作为 D-P材料按弹塑性分析或可用E—B材料以作非线性弹性分析。ANTTPS嘶＞；20G514:06:42rWTNO,TTPS图1-1 高桩码头桩基与岸坡相互作用空间有限元模型平面有限元模型如图1-2所示。计算模型宽度在码头前沿以前和后方堆场可各取码头宽的2.5倍，计算模型深度可从码头面到桩端以下5〜10米。横梁和桩可用梁单元或平面单元，其模量和尺寸可按相似原理进行换算。不管空间还是平面模型，在计算模型围的原状土，自重引起的应力场应按初始应力场考虑，ANSY提供了初应力计算的功能，解决初应力问题非常方便。、计算所用ANSY&B肯一的E—B模型说明

本算例是一个高桩码头桩基与岸坡相互作用的线性有限元和非线性有线元

分析，是一个平面计算的例子。本例的土体用邓肯—的 E—B模型模拟。E—B模型为非线性弹性模型，土的应力应变关系为双曲线，如图1-3所示。切线变形模量为EtEtKPaPEtKPaPa2Rf(1 3)(1Sin)2ccos23sin(1—1)(1—(1—2)(1—3)切线泊松比卩按下式确定:Et3(1 2)BKbpa(h式中，K、n—试验常数；pa—大气压；B—体变模量；m—应力指数；1、 3—大小主应力；心一体积模量系数；Rf—破坏比；c、—凝聚力和摩擦角。上述非线性弹性问题的有限元求解方法，可用增量法或迭代法。本例选用迭代法求解。ANSYS!供参数化设计语言APDL可以方便用于按参数建模，也可以用来扩展有限元的分析能力，创建各种控制方案。另外， ANSYS勺命令流也可以有效控制和实现很多求解过程和求解方案。本例结合ANSYS^数化设计语言及ANSYSP肯—的E—B模型来实现。三、计算有限元模型及计算结果题目二大圆筒结构、波浪与地基的相互作用分析一、用三维有限元模型计算大圆筒在外荷载（包括动荷载）作用下和地基相互作用问题，可按空间问题求解。现以大圆筒在振动下沉时的动力分析为例阐述如下。大圆筒求解有限元模型如图2-1所示。用六面体等参元模拟土基础，用曲

面壳体单元模拟钢圆筒，用梁单元模拟钢圆筒柱壳上的纵横向加强肋， 在大圆筒和基土之间设置接触面单元以模拟筒体和地基的相互作用。对瞬时动荷载采用时程积分的方法，将一个周期的荷载分为足够多个时间步（比如300步）进行求解。图2-2所示为大圆筒下沉到十六米时在T=11.5秒时的竖向应力云图。大圆筒在侧向荷载（如波浪力）作用下的倾复问题也可用上述空间模型求解。在侧向荷载作用下，达到一定程度时，圆筒和土体之间的接触单元会反映出筒体和土体间的分离现象。二、用三维有限元模型计算参考日本规，大圆筒的位移和倾复问题，也可以用平面方法按如下模型求解。纵向取单宽按平面问题计算，横向宽度B可按纵向单元长度面积相等的原则确定。墙后土压力按假想平面和库仑极限主动土压力计算，侧向摩阻力由主动土压力和土与墙背间摩擦角 计算确定。墙前水平土抗力通过设置非线性弹簧元经非线性迭代计算确定。筒底地基反力通过设置筒底竖向弹簧元经计算确定；筒底基土的剪切反力，通过设置筒底水平向弹簧元计算确定，根据上述假定，在ANSY上建立有限元模型如图2-3。经非线性迭代计算，得到土抗力和地基反力如图 2-4所示

siao/731^位移图见图2-5，最大位移DMXf10.7688cm。各关键点的位移如表表一关键点的位移关 键 点竖向位移（cm）水平位移（cm）水平允许最大位移H*1.5%（cm）筒 顶2.88（向下）7.41（向外）V29.25（满足要求）海 底 面2.84（向下）5.04（向外）筒底前趾2.42（向下）0.179（向）筒底后趾1.73（向下）0.566（向外）CT.：：F：-¥FT!FT77STBF^lSUB=1TIbgDt£<=OLimftBDn

题目二（冲击荷载下）强夯地基固结有限元分析（提供固结方程或固结方程处理方案，孔隙水压力消散计算方案、沉降计算方案及其他一些处理技巧）一、计算说明本题目对强夯地基固结进行有限元分析， 可采用ANSYSSTRUCTURE析模块来完成。二、计算有限元模型及边界条件图3—2冲击荷载加载图三、强夯地基固结计算ANar?-idTIHEh,1

NQDMGEONSTIP=1NQDMGEONSTIP=1SUB=10TIKE-.LSEQVIMflG)Dffi：-.7«E-03£MH=21^7=114591XAF102(05L2i&4i4O图3—4冲击荷载作用下的等效应力云图根据Skempton土力学计算公式，每一个荷载增量导致的孔隙水压力的增量可用式(3-1)表示：uB[3A(i 3)],B1,A1/3 (3-1)根据上式(3-1)，采用ANSY创计算得到动载作用下孔隙水压力增量。图3-5动荷载末时的孔隙水压力增量 U云图图3-6动荷载末时的孔隙水压力增量u云图由式：uu0ughu (3-2)可计算得到孔隙水压力的分布，如下图所示。图3—7动荷载末时的孔隙水压力 U云图孔隙水压力与水头的关系可用式(3—3)表示，图3—8为冲击何载作用末时水头分布云图。(3—3)Ef：r*JLLSOLWIOI；ANxnrw±00^UlltlQlSTJBTIMK-1TEXP■ar■曹丄u空斗=10TiH710.Q1J L9・1$6 1叽彳 丄0』士电£ 10.&a-?图3—8冲击载荷末时水头h图3-101天后各点水头分布图3-114天后各点水头分布图3-125天后各点水头分布图3—1310天后各点水头分布图3—1420天后各点水头分布图3—1530天后各点水头分布图3-1650天后各点水头分布图3—1760天后各点水头分布图3-1865天后各点水头分布 图3-19100天后各点水头分布为了得到在冲击载荷作用下孔隙水压力的消散过程， 如图3-21〜图3-23所示分别画出了图3-20所示3个点的孔隙水压力u随时间的变化曲线。10m10m20m图3-20孔隙水压力U消散过程三个监控点位置

4.60E+04►4.50E+0444.40E+04 |4.30E+044.20E+044 *4.10E+04—♦ I400E+04-2.00E+06 0.00E+002.00E+064.00E+066.00E+068.00E+061.00E+07时间（秒）图3-21第1个监控点孔隙水压力 U随时间的消散过程图3-22第2个监控点孔隙水压力U随时间的消散过程0.00E+ 1.00E+ 2.00E+ 3.00E+4.00E+5.00E+6.00E+7.00E+ 8.00E+ 9.00E+00060606060606060606时间（秒）图3-23第3个监控点孔隙水压力U随时间的消散过程题目四:在降雨情况下土工格栅加筋土挡墙边坡上公路稳定分析（由上至下为公路面层，垫层，挡墙，挡墙面板采用预制混凝土块0.60.60.6m3,混凝土后方为钩挂式土工格栅，边坡比较陡，边坡有一定排水特性）。一、计算说明题目四是一渗流—结构耦合计算问题。而ANSYS具有功能强大的渗流—结构耦合计算功能。下面是针对题目四，采用ANSY新完成的：在降雨情况下土工格栅加筋土挡墙边坡上公路稳定分析。二、计算几何模型和有限元模型以下各图分别是计算模型及有关尺寸、材料的说明

MftrHUMMftrHUM图4-3渗流-结构耦合计算有限元模型该坝体结构共有六种材料，上图4-3中分别用不同颜色表示结构不同的材料。三、计算结果1、计算边界条件

以上是边坡路基在降雨情况下稳定分析边界条件说明。除了考虑降雨，同时路基两侧考虑了10m常水头的作用。2、计算结果STEPwlSUB7TIME-iSTEPwlSUB7TIME-i卜雨情况下压力水头等值云图HODAL8OLUTICNffTEP-1卜雨情况下压力水头等值云图HODAL8OLUTICNffTEP-1SUB=1TIME-lSK?V 刖SWI-2€0.4«L吕HK=l^OBOdF雨情况下孔隙水压力WvDAuuJ-OLUTLl#WvDAuuJ-OLUTLl#旨3VEIwlTIME->2盘网 iATO)ere? 酉SMH-7&3.-1：gm>M0L2：4DSKK£i-4eiE4«IS4169 5212&9 315219 4.01240F雨情况下结构等效应力云图同时，为了对下雨和无雨情况下边坡路基的稳定性进行对比，下面是无雨情