《隧道工程章隧道有限元.pptx》由会员分享，可在线阅读，更多相关《隧道工程章隧道有限元.pptx（51页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、中国矿业大学建筑工程学院-李元海1主要内容7.1 有限元是什么？7.2 有限元法用到哪些基础知识？7.3 隧道围岩弹塑性有限元分析？如何用7.4 实战解析（新大瑶山隧道）第1页/共51页中国矿业大学建筑工程学院-李元海2前 言1965年“有限元”名词首次出现，到今天在工程上广泛应用，经历40多年发展史，理论和算法都日趋完善。有限元核心思想是结构离散化，即将实际结构假想地离散为有限数目的规则单元组合体，实际结构物理性能可以通过对离散体进行分析，得出满足工程精度的近似结果来替代对实际结构的分析，这样可以解决很多实际工程需要解决而理论分析又无法解决的复杂问题。第2页/共51页中国矿业大学建筑工程学院
2、-李元海3思考题1.什么是有限元分析？2.简述有限元的分析步骤？3.岩土材料本构关系包含哪四个部分？4.隧道中最常用有限元分析软件有哪些？第3页/共51页中国矿业大学建筑工程学院-李元海47.1 概述7.1.1 数值分析方法简介1.隧道工程问题解析解求解困难。2.数值分析方法是一种相对于弹性力学精确解析解的近似解求解方法。3.数值分析方法包含有限元法、有限差分法、边界元法等。4.有限元法是将弹性理论、计算数学和计算机软件有机结合。第4页/共51页中国矿业大学建筑工程学院-李元海57.1 概述有限元分析是利用数学近似方法对真实物理系统（几何和载荷工况）进行模拟。还利用简单而又相互作用的元素，即单
3、元，用有限数量的未知量去逼近无限未知量的真实系统。第5页/共51页中国矿业大学建筑工程学院-李元海6有限元模型 是真实系统理想化的数学抽象。真实物理系统有限元模型第6页/共51页中国矿业大学建筑工程学院-李元海7节点和单元，有限元模型节点:空间中的坐标位置，具有一定自由度和 存在相互物理作用。单元:一组节点自由度间相互作用的数值、矩阵 描述（称为刚度或系数矩阵)。单元有线、面或实体以及二维或三维的单元等种类。有限元模型由一些简单形状的单元组成，单元之间通过节点连接，并承受一定载荷。载荷载荷第7页/共51页中国矿业大学建筑工程学院-李元海8单元、边界、几何模型单元节点边界条件真实物理系统有限元模
4、型第8页/共51页中国矿业大学建筑工程学院-李元海97.1 概述7.1.2 有限元的发展概况（P221）7.1.3 有限元法软件简介1.大型通用软件：ANSYS，MARC2.专用软件：FLAC，MADAS3.自编软件：第9页/共51页中国矿业大学建筑工程学院-李元海107.1 概述v国际早在20世纪50年代末、60年代初就投入大量人力和物力开发有强大功能的有限元分析程序。最为著名的是由美国国家宇航局（NASA）在1965年委托美国计算科学公司和贝尔航空公司开发的NASTRAN有限元分析系统。该系统发展至今已有几十个版本，如MARC,是目前世界上规模最大、功能最强的有限元分析系统。此外，还有美国
5、的ANSYS，ABAQUS、ADINA等公司产品。第10页/共51页中国矿业大学建筑工程学院-李元海117.1.3 有限元法软件简介-ANSYSv美国ANSYS公司成立于1970年，创始人是John Swanson博士，总部位于美国宾西法尼亚州匹兹堡。vANSYS是个通用有限元仿真分析软件，早期只提供热分析和线性结构分析功能，只能运行在大型计算机上。v20世纪70年代后，增加非线性计算功能。随PC机和交互OS的出现，ANSYS建立了交互式菜单环境，程序除求解器，同时提供前后处理器。vANSYS12.0，能同时进行结构、热、流体、电磁、声学高级多物理场耦合分析，多场耦合分析技术在现今世界首届一指
6、。第11页/共51页中国矿业大学建筑工程学院-李元海12第12页/共51页中国矿业大学建筑工程学院-李元海137.1.3 有限元法软件简介-FLAC-岩土工程专用FLAC3D简介：FLAC3D由美国Itasca公司开发的。二维计算程序V3.0以前为DOS版本，V2.5版本仅能够使用计算机的基本内存（64K），程序求解的最大结点数仅限于2000个以内。1995年，FLAC2D已升级为V3.3的版本，程序能够使用扩展内存。FLAC3D是一个三维有限差分程序，目前已发展到V4.0版本。FLAC3D是能够进行土 质、岩石和其它材料三维结构受力特性模拟和塑性流动分析。能进行大变形分析。第13页/共51页
7、中国矿业大学建筑工程学院-李元海14第14页/共51页中国矿业大学建筑工程学院-李元海157.2 有限元基础7.2.1 有限元分析步骤（原理上）1.受力连续体“离散化”，仅节点联系，节点传力；2.单元外力转化为等效节点力3.建立节点力的平衡方程式4.加入位移边界条件求解方程组5.得到全部未知位移后，进而求得的单元应变与应力第15页/共51页中国矿业大学建筑工程学院-李元海167.2.2 平面问题的有限元分析以弹性力学平面问题有限元分析为例，介绍有限元基本思想、原理和分析步骤。核心是建立节点平衡方程组第1步：结构离散 （网格划分）第16页/共51页中国矿业大学建筑工程学院-李元海177.2.2
8、平面问题的有限元分析第2步：单元分析（P223）-荷载与位移关系（1）节点位移与节点力的表示形式。（2）节点位移与单元任一点位移关系。（3）节点位移与应变关系。（4）节点位移与应力关系。（5）节点位移与节点力的关系。第17页/共51页中国矿业大学建筑工程学院-李元海187.2.2 平面问题的有限元分析第3步：整体分析（1）单元贡献矩阵：一个单元（2）整体刚度矩阵集成：多个单元第4步：荷载移置将不在节点上的外荷载按虚功等效原则将荷载移到节点上。第18页/共51页中国矿业大学建筑工程学院-李元海197.2.2 平面问题的有限元分析第5步：引入支承或边界条件第6步：解方程组求出节点位移第7步：求出单
9、元应变与应力7.2.3 较精密的平面单元第19页/共51页中国矿业大学建筑工程学院-李元海207.3 隧道弹塑性有限元分析隧道岩土和结构材料具有弹塑性性质材料应力与应变关系非线性特征明显关键建立岩土材料的弹塑性本构关系求解非线性方程组第20页/共51页中国矿业大学建筑工程学院-李元海217.3.2 非线性问题的求解方法非线性方程组求解方法：直接迭代、切线刚度1.线性问题 ：K中元素为常量2.非线性问题：K中元素为变量（1）材料非线性（2）几何非线性第21页/共51页中国矿业大学建筑工程学院-李元海227.3.2 非线性问题的求解方法非线性方程组求解方法：直接迭代、切线刚度2，1.99999，2
10、.00001第22页/共51页中国矿业大学建筑工程学院-李元海237.3.3 岩土材料的弹塑性本构关系岩土材料本构关系包括四个组成部分（1）屈服条件和破坏条件（2）硬化定律（3）流动法则（4）加载和卸载准则第23页/共51页中国矿业大学建筑工程学院-李元海241.几种常用的屈服准则（1）摩尔-库仑屈服准则固体内任一点发生剪切破坏时，破坏面上的剪应力()应等于或大于材料本身的抗切强度(C)和作用于该面上由法向应力引起的摩擦阻力(tg)之和。第24页/共51页中国矿业大学建筑工程学院-李元海251.几种常用的屈服准则按照库仑-纳维尔理论，岩石的强度包络线是一条斜直线，破坏面与最小主平面的夹角恒等于
11、45/2。库仑-纳维尔判据适用于坚硬、较坚硬的脆性岩石产生剪切破坏的情况，而不适用于拉破坏的情况。该判据没有考虑中主应力2的影响。第25页/共51页中国矿业大学建筑工程学院-李元海261.几种常用的屈服准则德鲁克-普拉格（rucker-Prager）判据适用于以延性破坏为主的岩石。优点是考虑了中间主应力的作用。第26页/共51页中国矿业大学建筑工程学院-李元海277.3.2 非线性问题的求解方法2.硬化法则：材料进入塑性变形后应力应变关系3.流动法则：塑性应变增量分量和应力分量以及应力增量分量之间关系。4.加载、卸载准则：与本构关系选择相关5.弹塑性应力与应变关系：本构关系第27页/共51页中
12、国矿业大学建筑工程学院-李元海287.3.4 隧道围岩弹塑性有限元分析方法与步骤1.分析步骤1）确定岩体结构材料参数和施工模拟方法；2）划分围岩和结构有限元网格；3）分阶段计算等效荷载；4）分阶段计算节点与单元的位移、应变与应力；5）叠加各阶段计算值；6）计算结果分析与围岩结构的稳定度及稳定性评价第28页/共51页中国矿业大学建筑工程学院-李元海297.3.4 隧道围岩弹塑性有限元分析方法与步骤2.施工阶段 A.A.开挖前（初始应力）开挖前（初始应力）B.B.上台阶开挖，初衬上台阶开挖，初衬C.C.中台阶开挖，初衬中台阶开挖，初衬D.D.下台阶开挖，初衬下台阶开挖，初衬E.E.二衬二衬F.F.
13、仰拱仰拱第29页/共51页中国矿业大学建筑工程学院-李元海307.4 工程实例分析实例一：分阶段开挖和衬砌的弹塑性分析1.施工过程分析2.网格划分3.参数选择、力学模型选择4.计算5.结果汇总、整理、分析6.结论实例二：锚杆与围岩相互作用的弹塑性分析第30页/共51页中国矿业大学建筑工程学院-李元海317.4 工程实例分析实例一：分阶段开挖和衬砌的弹塑性分析1.施工方法及施工阶段分析全断面施工；初支：喷砼；二次支护：施作混凝土；（1）开挖全断面至边墙底部，喷射10cm混凝土；（2）施作二衬（3）开挖底部（4）浇灌仰拱第31页/共51页中国矿业大学建筑工程学院-李元海32网格划分第32页/共51
14、页中国矿业大学建筑工程学院-李元海33实例一：分阶段开挖和衬砌的弹塑性分析输进垃圾，输出垃圾材料参数选择一定要正确第33页/共51页中国矿业大学建筑工程学院-李元海34实例一：分阶段开挖和衬砌的弹塑性分析衬砌轮廓变形第34页/共51页中国矿业大学建筑工程学院-李元海35实例一：分阶段开挖和衬砌的弹塑性分析围岩位移等值线第35页/共51页中国矿业大学建筑工程学院-李元海36实例一：分阶段开挖和衬砌的弹塑性分析衬砌关键部位最大最小应力第36页/共51页中国矿业大学建筑工程学院-李元海37实例二：锚杆与围岩相互作用的弹塑性分析第37页/共51页中国矿业大学建筑工程学院-李元海38实例二：锚杆与围岩相
15、互作用的弹塑性分析位移场应力场第38页/共51页中国矿业大学建筑工程学院-李元海39计算实例隧道开挖方法优化三维有限元分析工作步骤第一步 明确研究目的 第39页/共51页中国矿业大学建筑工程学院-李元海40第二步 了解工程基础资料 隧道隧道断面断面第40页/共51页中国矿业大学建筑工程学院-李元海41第三步 了解隧道施工方法 隧道开挖方法第41页/共51页中国矿业大学建筑工程学院-李元海42第四步第四步 数值几何模型建模范围选择数值几何模型建模范围选择第42页/共51页中国矿业大学建筑工程学院-李元海43第五步第五步 创建三维几何模型创建三维几何模型-综合建模方法综合建模方法 第43页/共51
16、页中国矿业大学建筑工程学院-李元海44第六步 建立三维“力学”模型 第44页/共51页中国矿业大学建筑工程学院-李元海45第七步第七步 锚杆与衬砌布置与模拟方法锚杆与衬砌布置与模拟方法第45页/共51页中国矿业大学建筑工程学院-李元海46第八步第八步 整理和分析位移计算结果整理和分析位移计算结果 第46页/共51页中国矿业大学建筑工程学院-李元海47第八步第八步 整理和分析位移计算结果整理和分析位移计算结果第47页/共51页中国矿业大学建筑工程学院-李元海48第八步第八步 整理和分析位移计算结果整理和分析位移计算结果第48页/共51页中国矿业大学建筑工程学院-李元海49第八步第八步 整理和分析位移计算结果整理和分析位移计算结果第49页/共51页中国矿业大学建筑工程学院-李元海50第九步第九步 编写计算分析报告编写计算分析报告第50页/共51页中国矿业大学建筑工程学院-李元海51感谢您的观看！第51页/共51页