2020年新版基坑开挖数值模拟.docxVIP

PAGE PAGE # / 33 7 数值模拟 数值模拟方法简介 数值模拟技术作为一种研究手段， 已经被广泛的应用于各行各业 领域的研究中。目前，数值分析方法主要分为二大类：一类是以有限 差分法为代表，其特点是直接求解基本方程和相应的定解条件的近似 解；另一类数值分析方法是首先建立和原问题基本方程及相应定解条 件等效的积分方法，然后据之建立近似解法。 LS-D YNA乍为世界上最著名的通用显示动力分析程序，能够模拟 真实世界的各种复杂问题， 特别适合求解各种二维三维非线性结构的 高速碰撞， 爆炸和金属成型等非线性动力冲击问题， 同时可以求解传 热，流体及流固耦合问题， 在工程应用如汽车安全设计，武器系统设 计，金属成型，跌落仿真等领域被广泛应用。 本次采用 ANSYS/LS-DYN，A 进行混凝土支撑梁结构爆破拆除数值模拟研究。在ANSYS/LS-DYN环 境下，数值模拟的实现总体上分为两个过程： 在ANS丫芽建立结构实 体模型，完成有限元网格的划分， 输出有限元模型信息即输出关键字 文件；编辑关键字文件，在 DYNA环境下完成对结构倒塌过程的数值 模拟计算。 对结构有限元模型的建立过程， 数值模拟中采用的钢筋和混凝土 材料模型、接触方式等各种计算控制项进行了阐述。 LS-D YNA程序中主要提供如下几种计算方法： （ 1） Lagrange 算法 坐标固定在物质上或者说随物质一起运动和变形， 处理自由面和 物质界面非常直观， 由于网格始终对应物质， 因此能够精确的跟踪材 料边界和描述物质之间的界面，这是 Lagrange 算法的主要优点。但 是，由于网格随材料流动而变形，一旦网格变形严重，就会引起数值 计算的不稳定， 甚至使得计算无法继续进行 （如发生负体积或复杂声 速等问题）。因此， Lagrange 算法在处理大变形大位移问题时，有其 无法克服的弊端。 （ 2）Euler 算法 网格被固定在空间，是不变形的。物质通过网格边界流进流出， 物质的大变形不直接影响时间步长的计算。 因此，欧拉算法在处理大 变形问题方面具有优势。欧拉方法通过输运项计算体积、质量、动量 和能量的流动。 欧拉计算可以直接通过在离散化格式中包括迁移导数 项进行， 或通过二步操作完成。 二步法操作的第一步主要是拉格朗日 计算，第二步输运阶段是重分计算网格相当于回到它的原来状态。 LS-D YNA程序采用后一种方法。欧拉算法的缺点是网格中物质边界不 清晰，难以捕捉各物质界面。 （3）ALE方法 吸取了欧拉法和拉格朗日法两种方法的优点。 ALE算法能够进行 自动重分网格操作。它包括拉格朗日时间步，然后是一个输运步。输 运步可以采用三种方法： 1. 发生合理的网格变形时空间网格不再重 分（拉格朗日）；2. 发生严重的网格变形时重分成原始形状 （欧拉）； 3. 发生严重的网格变形时重分为合理的形状， 因此允许网格拓扑（拉 格朗日和欧拉）。 混凝土是土木工程结构中应用极为广泛的材料， 其最本质的特点 是材料组成的不均匀性， 并且存在初始微裂缝。 从混凝土受单轴压力 时的应力应变关系来看， 混凝土卸载时有残余变形， 不符合弹性关系； 如果对其应用弹塑性本构关系，又很难精确定义屈服条件。此外，混 凝土在到达应力顶峰后，其应力 - 应变关系曲线有一下降段，即存在 应变软化现象， 所有这些都给建立混凝土的本构关系带来困难。 多年 以来，众多学者进行了大量的试验和理论研究， 提出了各种各样的混 凝土本构模型。 第一：混凝土本构模型可以分为下面几种： （1）线弹性类本构模型。线弹性类本构模型是以弹性力学为基 础的模型， 当混凝土无裂缝时， 将混凝土看成线弹性匀质材料而采用 线弹性本构模型。 虽然混凝土的变形特性是非线性的， 但在一些特定 的情况下 （比如描述混凝土受拉时的工作性能 ） ，采用线弹性类本构模 型进行分析还是有足够的精度的， 其线弹性本构关系可用广义虎克定 律来表示： ij Cijkl kl （ 1） 式中，Cjki为材料弹性常数，为四阶张量，共有 81个常数。按照 材料假设的不同， 又可分为各向异性本构模型、 正交各向异性本构模 型、各向同性本构模型等，其中 Cijkl 根据材料的不同而变换。 （2）塑性理论类本构模型：塑性理论类本构模型是以塑性流动 理论为基础，代表性的模型主要有： Mises 条件的模型、理想弹塑性 脆性断裂模型、 应变或工作硬化塑性理论模型等。 模型中考虑了混凝 土加载路径和混凝土的硬化，在混凝土的应力 - 应变全曲线中，有上 升段和下降段。自从Drucke公设和nyushin公设出现之后，经典塑 性力学得到飞速发展， 混凝土塑性力学模型也是基于这些公设建立的。 以塑性理论为基础的混凝土本构模型， 在对其加载面， 包括初始屈服 面