第一章 计算机工程分析前处理和后处理.docVIP

第一章 计算机工程分析的前处理和后处理 本章教学任务要点：通过本章的系统学习，要求学生能够掌握有限单元法计算的前、后处理基本知识，基本原理和基本方法。 §1. 工程对象及选择 工程对象：进行西安城区地裂缝模拟： 研究目标：未来100年，地裂缝活动特点、以及对市政规划的影响。 选择对象： 工程对象：城市地铁工程隧道模拟 研究目标：未来100年，地裂缝活动特点、以及对市政规划的影响。 选择对象： A位置：隧道埋深2m，地面为繁华大街； B位置：隧道埋深10m，地面为高层建筑。 工程对象：煤层开采模拟 研究目标：采动后岩层运动规律、离层破坏高度 选择对象： §2. 计算模型的选取和建立 有限元计算模型的原则：计算模型要能全面反映工程对象的主要特点，又必须具备能适应计算模拟的功能。 工程对象：进行西安城区地裂缝模拟 城市地铁隧道 煤层开采模拟 有限元计算模型的建立主要指，1、将待分析的连续体，如结构物、固体等对象，用假想的点、线、面将连续体分割成有限多个、有限大小的子区域，这些子区域只有在特定点相互连接，从而使得连续体离散化为结构体。其中，这些子区域称为单元，单元与单元之间的连接点称为结点。 单元类型有： 直线单元 曲线元 A 一维单元 三角形 矩形 四边形 曲边形 B 二维单元 六面体单元 柱体单元 C 三维单元 结点类型：结点与结点之间的连接可以有铰接、固接。如果结构体的一个结点位移或一个方向被限制，则结点上可以安装铰支座、或杆支座等。 上机 上机操作熟练 §3. 网格化技术： 3.1 网格化并描述： 单元， 单元号， 单元节点， 节点坐标， 生成数据文件： 节点坐标信息，单元号。 例图巨型划分成小巨型格，也就是将巨型（1000X1000）划分成小巨型（5格 X 5格），其中如：单元n、节点为A、B、C、D。该单元其坐标为： A（400，400），B（600，400），C（600，600），D（400，600）。 用传统解析方法定量地处理岩石力学问题，由于无法考虑复杂的岩石性态及某些明显的地质构造的影响，而存在着几乎是不可克服的困难，正是在这方面有限单元法具有突出的效能。 有限单元法使连续体离散化；可以反映复杂实际对象。 单元应力状态可以指定，因而可以反映岩体实际状态。 处理格式一致。 3.2 网格单元形式 实践证明，采用常应变三角形单元或具有线性应变的任意四边形单元均能达到工程所要求的精度。 在计算精度方面，高阶应变特性的单元对提高计算精度是有限度的，而三角形单元精度略次于四边形单元，且处理岩石某些性态及进行非线性分析时程序处理上不如四边形方便。 3.3 离散化 正确确定边界及其位置条件，网格的细度，对有限元分析结果有很大影响。在保证要求精度的情况下，确定一个最低限度的网格要求是很必要的。例如根据大量计算经验，对地下硐室围岩分析，网格的适度细化，按节点数来说，应不少300节点。对具有一个对称轴时，节点总数不少于150~200个；对具有两个对称轴时，节点数不少于600个。 对网格疏密度，在应力变化处，可密；否则疏。例如对地下硐室围岩1~1.5倍直径区域加密。 对研究范围的大小及边界位置：如圆孔应力集中问题中，计算精度与精确解的误差在10%以内，必要的边界条件，一般应在距坑硐中心不少于坑硐直径3.5 ~ 4倍。又如，对弹模很低（1×105以下）的岩石，则所取范围尚应适当增大。 应指明：采用“翻转应力法”模拟开挖，边界均以固定点来考虑。 3.4 网格剖分的处理方法 地下硐室围岩应力分布问题，单元划分情况可分两种类型。即取全部或部分对称部分。主要适用于方形、圆形及马蹄形断面。坐标轴原点取在巷道中心点，以具有相等夹角的射线及同心圆环将所研究的范围划分成四边形网格，然后再以较短的一条对角线将其分为两个三角形。 两个对称轴、 有一个对称轴时， （对，具有一个对称轴而取时，将座标轴转90°） 无对称轴，Kn=7，Kn， 注：四边形割成两个三角形，其共合边中点处的值，即为四边形的值。 设划分单元的射线总数为Km，同心环包括巷道周边及外边界总数为Kn，则所得节点总数为： 单元总数目为： 对于性状不规则的断面，可以先按规则的轮廓划分，不规则部分作为附加单元处理。附加单元的编