ANSYS岩土计算例子

ANSYS岩土计算例子ANSYS岩土计算例子ANSYS岩土计算例子ANSYS土工结构计算案例ANSYS－CHINA广州做事处2019年11月23日目录计算题目及计算要求说明···························1题目一·······································4一、计算说明···································4二、计算所用ANSYS邓肯－张的E－B模型说明···············5三、计算有限元模型及计算结果························6题目二·······································7一、用三维有限元模型计算···························7二、用三维有限元模型计算···························8题目三······································10一、计算说明··································10二、计算有限元模型及界线条件·······················10三、强夯地基固结计算·····························10题目四······································17一、计算说明··································17二、计算几何模型和有限元模型·······················17三、计算结果··································181、计算界线条件·······························182、计算结果··································193、结论···································201计算题目及计算要求说明题目一：高桩码头桩基与岸坡相互作用的线性有限元和非线性有线元解析题目二：大圆筒结构、波浪与地基的相互作用解析（大圆筒作为重力式码头结构，波浪为水平动荷载，门吊为竖向动荷载，地基为三层以上地基包括自抛碎石垫层、粘土层、粉细砂层和岩层，粉细砂层可能在波浪动荷载作用下液化造成圆筒推翻）题目三：（冲击荷载下）强夯地基固结有限元解析（供应固结方程或固结方程办理方案，孔隙水压力消逝计算方案、沉降计算方案及其他一些办理技巧）题目四：在降雨情况下土工格栅加筋土挡墙边坡上公路牢固解析（由上至下为公路面层，垫层，挡墙，挡墙面板采用预制混凝土块0.60.60.6m3，混凝土后方为钩挂式土工格栅，边坡比较陡，边坡有必然排水特点）。详尽办理方案包括：1、供应计算输入界面2、计算模型或采用本构情况3、前办理方案及网格划分技巧4、特别资料或模型嵌入技术5、计算技巧及解决方案6、后办理供应内容2详尽报价方案包括：进行简单报价，涉及以上题目的各模块的综合报价（包括前后办理及解决以上问题的结构与岩土问题的模块报价.能用通用模块计算尽量用通用模块,必定用CivilFEM模块计算的请注明.供应解决方案时间：2005年5月18日以前。3题目一高桩码头桩基与岸坡相互作用的线性有限元和非线性有线元解析一、计算说明高桩码头桩基与岸坡相互作用的线性有限元和非线性有线元解析，可以用空间有限元模型或平面有限元模型数值求解。空间有限元模型如图1-1所示。考虑到岸坡主若是平面问题，为了节约单元，可以只取四榀排架建立空间模型，并以中间两个排架的受力情况作为整个码头段的桩的受力标准。横梁可用梁单元或体单元，桩可用梁单元或壳单元或体单元，地基土可用体单元，桩土之间可加面－面接触单元Targe170和contal174以模拟桩土之间的相互作用。梁和桩作为线弹性资料，土可作为D－P资料按弹塑性解析或可用E－B资料以作非线性弹性解析。图1-1高桩码头桩基与岸坡相互作用空间有限元模型平面有限元模型如图1-2所示。计算模型宽度在码头前沿以前和后方堆场可各取码头宽的2.5倍，计算模型深度可从码头面到桩端以下5～10米。横梁和桩可用梁单元或平面单元，其模量和尺寸可按相似原理进行换算。无论空间还是平面模型，在计算模型范围内的原状土，自重引起的应力场应按初始应力场考虑，ANSYS供应了初应力计算的功能，解决初应力问题特别方便。4二、计算所用ANSYS邓肯－张的E－B模型说明本算例是一个高桩码头桩基与岸坡相互作用的线性有限元和非线性有线元解析，是一个平面计算的例子。本例的土体用邓肯－张的E－B模型模拟。E－B模型为非线性弹性模型，土的应力应变关系为双曲线，如图1-3所示。切线变形模量为Et。nRf(13)(123sin)EtKpapa1（1－1）2ccos23sin切线泊松比μ按下式确定：BEt（1－2）2)3(13m（1－3）BKbpa(pa)式中，K、n—试验常数；pa—大气压；B—体变模量；m—应力指数；1、3—大小主应力；Kb—体积模量系数；Rf－破坏比；c、－凝聚力和内摩擦角。图1-3土的应力应变关系上述非线性弹性问题的有限元求解方法，可用增量法或迭代法。本例采用5迭代法求解。ANSYS供应参数化设计语言APDL，可以方便用于按参数建模，也可以用来扩展有限元的解析能力，创办各种控制方案。别的，ANSYS的命令流也可以有效控制和实现很多求解过程和求解方案。本例结合ANSYS参数化设计语言及ANSYS邓肯－张的E－B模型来实现。三、计算有限元模型及计算结果图1-4为计算有限元模型，全用采用ANSYS自动四边形网格功能来划分，从而保证了较高的计算精度。图1-4高桩码头桩基与岸坡相互作用计算元模型图1-5高桩码头桩基与岸坡相互作用变形图6题目二大圆筒结构、波浪与地基的相互作用解析一、用三维有限元模型计算大圆筒在外荷载（包括动荷载）作用下和地基相互作用问题，可按空间问题求解。现以大圆筒在振动下沉时的动力解析为例阐述以下。大圆筒求解有限元模型如图2-1所示。用六面体等参元模拟土基础，用曲面壳体单元模拟钢圆筒，用梁单元模拟钢圆筒柱壳上的纵横向加强肋，在大圆筒和基土之间设置接触面单元以模拟筒体和地基的相互作用。对瞬时动荷载采用时程积分的方法，将一个周期的荷载分为足够多个时间步（比方300步）进行求解。图2-2所示为大圆筒下沉到十六米时在T＝11.5秒时的竖向应力云图。图2-1大圆筒有限元模型图2-2大圆筒下沉到十六米时在T＝11.5秒时的竖向应力云图大圆筒在侧向荷载（如波浪力）作用下的倾复问题也可用上述空间模型求7解。在侧向荷载作用下，达到必然程度时，圆筒和土体之间的接触单元会反响出筒体和土体间的分别现象。二、用三维有限元模型计算参照日本规范，大圆筒的位移和倾复问题，也可以用平面方法按以下模型求解。纵向取单宽按平面问题计算，横向宽度B可按纵向单元长度内面积相等的原则确定。墙后土压力按假想平面和库仑极限主动土压力计算，侧向摩阻力由主动土压力和土与墙背间摩擦角计算确定。墙前水平土抗力经过设置非线性弹簧元经非线性迭代计算确定。筒底地基反力经过设置筒底竖向弹簧元经计算确定；筒底基土的剪切反力，经过设置筒底水平向弹簧元计算确定，依照上述假定，在ANSYS上建立有限元模型如图2-3。图2-3经非线性迭代计算，获取土抗力和地基反力如图2-4所示。8图2-4位移图见图2-5，最大位移DMX＝10.7688cm。各要点点的位移如表一。表一要点点的位移关键点竖向位移（cm）水平位移（cm）水平赞同最大位移H＊1.5%（cm）筒顶2.88（向下）7.41（向外）＜29.25（满足要求）海底面2.84（向下）5.04（向外）筒底前趾2.42（向下）0.179（向内）筒底后趾1.73（向下）0.566（向外）图2-5位移云图9题目三（冲击荷载下）强夯地基固结有限元解析（供应固结方程或固结方程办理方案，孔隙水压力消逝计算方案、沉降计算方案及其他一些办理技巧）一、计算说明本题目对强夯地基固结进行有限元解析，可采用ANSYS－STRUCTURE解析模块来完成。二、计算有限元模型及界线条件图3－1有限元模型及模型尺寸图3－2冲击荷载加载图三、强夯地基固结计算图3－3冲击荷载作用下的变形图10图3－4冲击荷载作用下的等效应力云图依照Skempton土力学计算公式，每一个荷载增量以致的孔隙水压力的增量可用式(3-1)表示：uB[3A(13)],B1，A1/3(3-1)依照上式(3-1)，采用ANSYS可计算获取动载作用下孔隙水压力增量。图3-5动荷载末时的孔隙水压力增量u云图11图3－6动荷载末时的孔隙水压力增量u云图由式：uu0ughu（3－2）可计算获取孔隙水压力的分布，以以下图所示。图3－7动荷载末时的孔隙水压力u云图孔隙水压力与水头的关系可用式（3－3）表示，图3－8为冲击荷载作用末时水头分布云图。huy（3－3）w12图3－8冲击载荷末时水头h图3－101天后各点水头分布图3－114天后各点水头分布图3－125天后各点水头分布图3－1310天后各点水头分布13图3－1420天后各点水头分布图3－1530天后各点水头分布图3－1650天后各点水头分布图3－1760天后各点水头分布图3－1865天后各点水头分布图3－19100天后各点水头分布为了获取在冲击载荷作用下孔隙水压力的消逝过程，如图3－21～图3－23所示分别画出了图3－20所示3个点的孔隙水压力u随时间的变化曲线。14图3－20孔隙水压力u消逝过程三个监控点地址4.60E+044.50E+04)4.40E+04P(力压4.30E+04水隙孔4.20E+044.10E+044.00E+04-2.00E+060.00E+002.00E+064.00E+066.00E+068.00E+061.00E+07时间(秒))a

图3－21第1个监控点孔隙水压力u随时间的消逝过程8.35E+048.30E+048.25E+048.20E+04(力压8.15E+04水隙8.10E+04孔8.05E+048.00E+047.95E+04-2.00E+060.00E+002.00E+064.00E+066.00E+068.00E+061.00E+07时间(秒)图3－22第2个监控点孔隙水压力u随时间的消逝过程151.24E+051.23E+05)1.23E+05aP1.22E+05(力压1.22E+05水隙1.21E+05孔1.21E+051.20E+051.20E+050.00E+1.00E+2.00E+3.00E+4.00E+5.00E+6.00E+7.00E+8.00E+9.00E+00060606060606060606时间(秒)图3－23第3个监控点孔隙水压力u随时间的消逝过程16题目四：在降雨情况下土工格栅加筋土挡墙边坡上公路牢固解析（由上至下为公路面层，垫层，挡墙，挡墙面板采用预制混凝土块0.60.60.6m3，混凝土后方为钩挂式土工格栅，边坡比较陡，边坡有必然排水特点）。一、计算说明题目四是一渗流－结构耦合计算问题。而ANSYS拥有功能富强的渗流－结构耦合计算功能。下边是针对题目四，采用ANSYS所完成的：在降雨情况下土工格栅加筋土挡墙边坡上公路牢固解析。二、计算几何模型和有限元模型以下各图分别是计算模型及有关尺寸、资料的说明图4-1计算几何模型图4-2计算模型尺寸及有关资料17图4-3渗流－结构耦合计算有限元模型该坝体结构共有六种资料，上图4-3中分别用不同样颜色表示结构不同样的材料。三、计算结果、计算界线条件图4-4渗流－结构耦合计算有限元模型18以上是边坡路基在降雨情况下牢固解析界线条件说明。除了考虑降雨，同时路