岩土工程数值分析课程论文周森.doc
《岩土工程数值分析课程论文周森.doc》由会员分享,可在线阅读,更多相关《岩土工程数值分析课程论文周森.doc(41页珍藏版)》请在咨信网上搜索。
1、研 究 生 课 程 论 文(2009-2010学年第二学期)岩土工程数值分析课程论文研究生:周森提交日期: 2010年7月1日 研究生签名:学 号200920105032学 院土木与交通学院课程编号S0814007课程名称岩土工程数值分析学位类别硕士任课教师刘庭金教师评语: 成绩评定: 分 任课教师签名: 年 月 日目 录例案1基于Midas-GTS的条形、方形基础下附加应力及位移分析 (3)例案2基于强度折减法的Midas-GTS二维边坡稳定性分析 (13)例案3基于Midas-GTS的二维隧道衬砌模拟分析 (20)例案4基于Midas-GTS的三维基坑施工阶段模拟 (27)例案5 岩土工程
2、实例介绍分析广州陈家祠广场及周边环境综合整治地下空间工程基坑支护结构设计 (31)基于Midas-GTS的条形、方形基础下附加应力及位移分析周 森(华南理工大学土木与交通学院 09岩土工程)摘 要:运用了Midas-GTS分别建立了条形基础和方形基础沉降及基底附加应力计算模型。条形基础考虑平面应变问题,采用二维坐标建模;方形基础考虑荷载的分布情况,则采用三维实体建模。通过程序分析计算,得到了两种基础下位移、附加应力的分布规律,得到了位移分布等值云图和“应力泡”分布,与土力学经典理论分析结果进行了对比,发现数值模拟分析结果与理论分析结果具有一致性,从而也验证了理论分析结果的正确性。关键词:条形基
3、础;方形基础;位移;附加应力;规律1 引 言地基内的附加应力是由基底附加压力引起的,建筑物荷载在地基中产生的应力为附加应力,其大小直接影响着地基的变形和稳定性,因此地基附加应力的计算也显得尤为重要。在地基附加应力的计算方面,目前采用弹性理论推导出的公式,并引入地基是半无限空间弹性体和地基土是均匀、连续、各向同性两方面的假设1。但是地基土的成层分布及应力历史等因素,使的地基土性参数呈现空间变异性2。基于问题的复杂性和不确定性,对于地基附加应力的计算并没有简单成熟的方法。大量试验表明,当地基土作用的荷载不是很大、土中塑性变形区很小时,荷载和变形之间近似成线性关系,基于弹性理论的计算结果与实测值相差
4、不大。在地基沉降位移的计算方面,也存在很大程度的不精确性,目前确定性的计算方法主要有分层总和法3和建筑地基基础设计规范法4。随着有限元分析技术的日益成熟,目前也有一批学者运用有限元方法对地基沉降进行了研究。本文通过建立条形基础和方形基础下地基附加应力和位移计算模型,分析了各自附加应力分布规律及位移变形特点,与经典理论进行了对比,以更加生动、形象的方式验证了理论的正确性。2 模型建立2.1 条形基础模型运用Midas-GTS二维建模,建立了条形基础下地基附加应力和位移计算的模型,边界条件采用平面应变。条形基础宽度B为2m,而边界距离取宽度的23倍,地基土深度取影响范围之内约6倍宽度5;不考虑地基
5、土的空间变异性,取平均重度=20KN/m3,泊松比取为0.25,弹性模量E取5e3kpa,土体内摩擦角为15,粘聚强度取25kpa;条形基础施加给地基的线荷载为200KN/m。条形基础下地基附加应力和位移计算模型、施加支撑及荷载后网格分别如图1、2和3所示。2.2 方形基础模型运用Midas-GTS建立了方形基础下地基附加应力和位移计算的三维实体模型,边界支撑分别在实体底部限制Z方向位移,左右两侧限制X方向位移,前后两侧限制Y方向位移。方形基础宽度B为4m,边界距离取宽度的2倍即8m,地基土深度取为20m;不考虑地基土的空间变异性,取平均重度=20KN/m3,泊松比取为0.25,弹性模量E取5
6、e3kpa,土体内摩擦角为15,粘聚强度取25kpa;方形基础施加给地基的均布面荷载为200KN/m2。方形基础下地基附加应力和位移计算模型、施加支撑及载荷后网格分别如图4、5和6所示。 图1 条形基础计算模型 图2 施加支撑后条形基础模型 图3 施加载荷后条形基础计算模型 图4 方形基础计算模型 图5 方形基础计算模型 图6 方形基础计算模型(正立面)3 分析结果查看3.1 条形基础结果运行Midas-GTS对条形基础下地基土的附加应力及位移进行分析计算,在结果表单中查看位移、内力和应力的各项分布,总结出一些特点和规律。位移分布、内力分布和应力分布分别如图712所示。图7 条形基础下地基X方
7、向位移图8 条形基础下地基Y方向位移从位移图可以看出,地基土在条形基础作用下,X方向产生向侧下方的挤压变形,呈蝶形分布类型;Y方向产生竖向位移,呈椭圆环形分布,在条形基础中心下较小土体范围内,沉降最大。图9 条形基础下地基FX反力图10 条形基础下地基FY反力从内力分布可以看出,条形基础地基FX反力最大发生在地基土侧壁离地面约2/3位置,FY反力最大发生在正对条形基础下的地基土地面位置,大致呈现直线分布特征。图11 条形基础下地基附加应力分布图12 条形基础下地基附加应力分布从应力分布图可以看到,条形基础下地基附加应力分布呈蝶形分布,应力分布等值线表明,最大应力分布在条形基础下侧一较小范围内,
8、应力影响较明显区域大致分布在1.5倍基础宽度范围;条形基础下地基附加应力分布呈椭圆形分布,应力分布等值线表明,最大应力分布在条形基础下侧一较小范围内,应力从内到外按照从大到小环形扩散分布,应力影响较明显的区域大致分布在4倍基础宽度范围,与理论分析有一定出入,原因在于建模时没有考虑土体自重的影响。3.2方形基础结果运行Midas-GTS对方形基础下地基土的附加应力及位移进行分析计算,在结果表单中查看位移、内力和应力的各项分布,总结出一些特点和规律。位移分布、内力分布和应力分布分别如图1320所示。图13 方形基础下地基X方向位移图14 方形基础下沿垂直Y剖切所得X方向位移图14 方形基础下地基Y
9、方向位移图15 方形基础下沿垂直X剖切所得Y方向位移图16 方形基础下地基Z方向位移图17 方形基础下沿垂直X剖切所得Z方向位移从X、Y和Z方向位移分布云图得知,X、Y方向位移分布基本相同,产生侧下方挤压变形,呈现蝶形分布特征;Z方向位移在方形基础下侧某一局部区域最大,而后按照环形递减的顺序向外扩散。图18 方形基础下地基土分布图19 方形基础下地基土分布图20 方形基础下地基土分布从X、Y和Z应力分布等值线得知,方形基础下地基附加应力在X和Y方向分布大致相同,大致呈曲线分布,蝶形特征逐渐消失;Z方向应力分布呈大致环形递减,但影响范围较浅。4 结 论运用Midas-GTS建立了条形基础和方形基
10、础的地基附加应力和位移计算模型,条形基础采用二维建模,方形基础采用三维实体建模。通过程序,分析了两种基础下地基土附加应力和位移分布特征,得到了如下结论:1)条形基础作用下地基土在X方向产生向侧向下的挤压变形,呈蝶形分布特征;Y方向产生竖向位移,呈椭圆环形分布,在条形基础中心下较小土体范围内,沉降最大。2)条形基础下地基附加应力呈蝶形分布,最大应力分布在条形基础下侧一较小范围内,应力影响较明显区域大致分布在1.5倍基础宽度范围;附加应力分布呈椭圆形分布,最大应力分布在条形基础下侧一较小范围内,应力从内到外按照从大到小环形扩散分布,应力影响较明显的区域大致分布在4倍基础宽度范围。3)方形基础下地基
11、在X、Y方向位移分布基本相同,产生侧下方挤压变形,呈现蝶形分布特征;Z方向位移在方形基础下侧某一局部区域最大,而后按照环形递减的顺序向外扩散。4)方形基础下地基在X和Y方向应力分布大致相同,大致呈曲线分布,蝶形特征逐渐消失;Z方向应力分布呈大致环形递减,但影响范围较浅。参考文献15 肖绍然.土力学M.郑州:郑州大学出版社,2007,1.2 包承钢,高大钊,张庆华.地基可靠度分析的一般理论M.武汉:武汉测绘科技大学出版社,1997,1-18.基于强度折减法的Midas-GTS二维边坡稳定性分析周 森(华南理工大学土木与交通学院 09岩土工程)摘 要:介绍了有限元强度折减法的基本原理,阐述了该方法
12、建立模型的原则、分析边坡的优点、边坡整体失稳的依据。结合建立的二维边坡模型,基于有限元强度折减法,运用有限元软件 Midas/GTS对该边坡进行了数值分析,得到了边坡在自重作用下的水平位移等值云图、应力与应变云图及安全系数指标,依此确定了最危险滑面位置和形状以及整体失稳的安全系数。关键词:有限元强度折减法;边坡;滑动面;安全系数1 强度折减法的基本原理极限平衡法和塑性极限分析法是土坡稳定分析中的传统方法,基于强度折减的有限元法用于边坡稳定分析是较新的方法。有限元强度折减法 (以下简称有限元法 )的基本原理是将坡体强度参数粘聚力和内摩擦角值同时除以一个折减系数F1-2,得到一组新的和值,即经过折
13、减后的抗剪强度指标为 (1) (2)然后将和作为新的计算参数输入,再进行试算;当计算不收敛时,对应的F被称为坡体的最小稳定安全系数,此时坡体达到极限状态,发生剪切破坏,同时可以得到坡体的破坏滑动面。应用有限元法需满足以下条件3:1)要有一个成熟的有限元程序;2)选择可供实用的弹塑性模型和强度屈服准则;3)计算范围、边界条件、网格划分要满足计算要求。强度折减法是基于有限元计算理论之上的边坡整体稳定分析方法,因此它具有有限元法的一切优点。与传统的极限平衡法相比,边坡稳定分析的有限元法的优点总结如下4。1)破坏面的形状或位置不需要事先假定,破坏自然地发生在土的抗剪强度不能抵抗剪应力的地带(或者称为岩
14、土体剪切带),在求解安全系数时,不需要假定滑动面的形状和位置,也无须进行条分,而是由程序自动求出滑动面,从而避免了人工划分滑动面过程中存在误差;2)由于有限元法全面满足了静力许可、应变相容和应力应变之间的本构关系,因此不必引入假定条件,保持了严密的理论体系;3)采用数值分析方法,可以不受边坡几何形状的不规则和材料的不均匀性的限制,因此是一种比较理想的分析边坡应力、 变形和稳定性态的手段,有限元解提供了应力变形的全部信息。Midas-GTS在运用于岩土工程数值分析以来,其精度可以满足大多数工程分析需要,且提供了众多可供选择的强度准则。本文采用Mohr-coulomb强度准则。2 Mohr-cou
15、lomb强度准则Mohr-coulomb破坏准则表示如下 (3)式中,为土体极限剪应力,土体为粘聚强度,土体为内摩擦角。其该强度具有简单而且准确的优点,至今被广泛应用于土质材料的分析中。该强度准则用应力不变量I1、应力偏量不变量J2和罗德角表示为5 (4)图1 Mohr-coulomb强度准则在主应力空间中的表示莫尔库仑破坏准则在土质材料的分析上有两个主要缺点。第一是没有考虑第二主应力对破的影响,这与试验结果不相符;第二是莫尔图形的子午线和破坏包络线是直线,内摩擦角不随约束压力(或者静水压力)变化。所以当约束压力在限制范围内时,该准则结果准确;但是当约束压力在限制范围以外时,准确性将会降低5。
16、正因为该准则在实用的约束压力范围内具有较高的准确性,并且使用简单,所以在岩土分析中被广泛应用。3 二维边坡模型建立的二维边坡模型分两个界面。第一个界面为粘土层,弹性模量E取1.0e5 kpa,粘聚强度取50kpa,内摩擦角取36,土体重度=20KN/m3,泊松比取0.3;第二个界面为软弱夹层,弹性模量E取1.0e4 kpa,粘聚强度取30kpa,内摩擦角取15,土体重度=20KN/m3,泊松比取0.3。边坡总高度为20m,边坡模型示意如图2所示。图2 边坡模型示意图采用有限元软件进行边坡力学分析时,如果按照真实的边坡体模型进行分析,是一个三维力学问题,分析起来耗费时间和计算机资源,其分析结果未
17、必很理想。根据弹性力学理论,将这种在纵向比较长的的边坡体结构简化为平面应变问题来代替三维模型,结果偏于安全6。位移边界条件的简化视实际约束的强弱可以简化为固定约束或者铰支约束等。边界条件为:下部边界固定,左右边界水平约束,采用四边形单元划分网格,可在坡肩与坡趾区域局部以及粘土层和软弱夹层交界面处加密网格尺寸,有利于提高有限元计算结果的精确度。Midas所建边坡模型及模型网格划分、归并分别见图3、4和5。图3 Midas边坡模型示意图图4 二维边坡网格划分图5 二维边坡网格归并在定义边坡支撑条件时,为方便期间,直接选择地面支撑,这样就在坡底施加了水平及竖向约束,在坡体两侧施加了水平方向约束。定义
18、坡体土体自重方向沿Y轴负方向,即重力方向竖直向下。施加边界支撑及自重的网格划分如图6所示。图6 施加边界支撑及自重的网格划分4 分析结果查看 定义施工工况后,运行Midas-GTS分析按钮,则程序自动按照强度折减进行分析。分析历时约660.37s,合约12mins。首先可以从结果表单中,查得稳定安全系数Fs为1.9375,如图7所示。图7 边坡安全系数Fs4.1 坡体位移由结果表单提取边坡在水平方向和在竖直方向上的位移。结果表明,边坡在水平方向上的最大位移发生在边坡坡趾位置,最大位移为0.061913m,方向水平向右;边坡在竖直方向上的最大位移发生在坡肩位置,最大位移为0.073278m,方向
19、竖直向下。边坡水平方向上的位移分布及位移向量云图如图8、9所示,竖直方向上位移分布及位移向量云图如图10、11所示。图8 边坡水平方向的位移 图9 边坡水平位移向量云图 图10 边坡竖直方向上的位移 图11 边坡竖直位移向量云图4.2 坡体应力由结果表单提取边坡在水平方向和在竖直方向上的应力。结果表明,边坡在水平方向和竖直方向上的最大应力均发生在边坡坡底脚部,出现应力集中现象。绝对值最大应力分别为1.665102 kpa和3.897102 kp。边坡水平和竖直方向上的应力如图10、11所示。图10 边坡水平方向应力分布 图11 边坡竖直方向应力分布4.3 坡体应变边坡坡体水平、竖直方向和剪应变
20、分布分别如图12、13和14所示,其中水平最大应变为=6.5810-2,竖直方向最大应变为=5.7710-2,最大剪应变为=0.151。图12 边坡水平方向应变图13 边坡竖直方向应变 图14 边坡剪应变分布5 结 论运用Midas-GTS对所建立的二维含软弱夹层边坡进行了强度折减分析,根据分析结果,得到如下结论。1)运用Midas-GTS进行强度折减时,初始安全系数设为1.0,然后根据在边坡稳定性分析(SRM)设置选项中,定义折减步数,可运行程序自动得到安全系数指标,比Ansys等有限元软件分析过程智能化。本文所建模型,分析所得安全系数指标Fs为1.9375。2)边坡水平最大位移发生在边坡坡
21、趾位置,最大位移为0.061913m,方向水平向右;在竖直方向上的最大位移发生在坡肩位置,最大位移为0.073278m,方向竖直向下。3)边坡水平和竖向最大应力均发生在边坡坡体脚部处,且出现应力集中现象。参考文献1 时卫民,郑颖人,张鲁渝.岩石高边坡的有限元分析及其简化分析方法 J .地下空间,2001,21 (5) .2 张彩双.有限元强度折减法的边坡稳定分析 J .中国农村水利水电,2006 (5) :99.3 郑颖人,陈祖煜.边坡与滑坡工程治理 M .北京:人民交通出版社,2007,194 - 204.4 郑颖人,赵尚毅.有限元强度折减法在土坡与岩坡中的应用 J .岩土力学与工程学报,2
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 岩土 工程 数值 分析 课程 论文
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【a199****6536】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【a199****6536】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。