基坑开挖数值模拟

《基坑开挖数值模拟》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、7数值模拟7.1数值模拟方法简介数值模拟技术作为一种研究手段，已经被广泛的应用于各行各业领域的研究中。目前，数值分析方法主要分为二大类：一类是以有限差分法为代表，其特点是直接求解基本方程和相应的定解条件的近似解；另一类数值分析方法是首先建立和原问题基本方程及相应定解条件等效的积分方法，然后据之建立近似解法。LS-DYNA作为世界上最著名的通用显示动力分析程序，能够模拟真实世界的各种复杂问题，特别适合求解各种二维三维非线性结构的高速碰撞，爆炸和金属成型等非线性动力冲击问题，同时可以求解传热，流体及流固耦合问题，在工程应用如汽车安全设计，武器系统设计，金属成型，跌落仿真等领域被

2、广泛应用。本次采用ANSYS/LS-DYNA，进行混凝土支撑梁结构爆破拆除数值模拟研究。在ANSYS/LS-DYNA环境下，数值模拟的实现总体上分为两个过程：在ANSYS中建立结构实体模型，完成有限元网格的划分，输出有限元模型信息即输出关键字文件；编辑关键字文件，在DYNA环境下完成对结构倒塌过程的数值模拟计算。对结构有限元模型的建立过程，数值模拟中采用的钢筋和混凝土材料模型、接触方式等各种计算控制项进行了阐述。LS-DYNA程序中主要提供如下几种计算方法：（1）Lagrange算法坐标固定在物质上或者说随物质一起运动和变形，处理自由面和物质界面非常直观，由于网格始终对应物

3、质，因此能够精确的跟踪材料边界和描述物质之间的界面，这是Lagrange算法的主要优点。但是，由于网格随材料流动而变形，一旦网格变形严重，就会引起数值计算的不稳定，甚至使得计算无法继续进行（如发生负体积或复杂声速等问题）。因此，Lagrange算法在处理大变形大位移问题时，有其无法克服的弊端。（2）Euler算法网格被固定在空间，是不变形的。物质通过网格边界流进流出，物质的大变形不直接影响时间步长的计算。因此，欧拉算法在处理大变形问题方面具有优势。欧拉方法通过输运项计算体积、质量、动量和能量的流动。欧拉计算可以直接通过在离散化格式中包括迁移导数项进行，或通过二步操作完成。二

4、步法操作的第一步主要是拉格朗日计算，第二步输运阶段是重分计算网格相当于回到它的原来状态。LS-DYNA程序采用后一种方法。欧拉算法的缺点是网格中物质边界不清晰，难以捕捉各物质界面。（3）ALE方法吸取了欧拉法和拉格朗日法两种方法的优点。ALE算法能够进行自动重分网格操作。它包括拉格朗日时间步，然后是一个输运步。输运步可以采用三种方法：1.发生合理的网格变形时空间网格不再重分（拉格朗日）；2.发生严重的网格变形时重分成原始形状（欧拉）；3.发生严重的网格变形时重分为合理的形状，因此允许网格拓扑（拉格朗日和欧拉）。混凝土是土木工程结构中应用极为广泛的材料，其最本质的特点是材料组

5、成的不均匀性，并且存在初始微裂缝。从混凝土受单轴压力时的应力应变关系来看，混凝土卸载时有残余变形，不符合弹性关系；如果对其应用弹塑性本构关系，又很难精确定义屈服条件。此外，混凝土在到达应力顶峰后，其应力-应变关系曲线有一下降段，即存在应变软化现象，所有这些都给建立混凝土的本构关系带来困难。多年以来，众多学者进行了大量的试验和理论研究，提出了各种各样的混凝土本构模型。第一：混凝土本构模型可以分为下面几种：（1）线弹性类本构模型。线弹性类本构模型是以弹性力学为基础的模型，当混凝土无裂缝时，将混凝土看成线弹性匀质材料而采用线弹性本构模型。虽然混凝土的变形特性是非线性的，但在一些特

6、定的情况下(比如描述混凝土受拉时的工作性能)，采用线弹性类本构模型进行分析还是有足够的精度的，其线弹性本构关系可用广义虎克定律来表示：（1）式中，为材料弹性常数，为四阶张量，共有81个常数。按照材料假设的不同，又可分为各向异性本构模型、正交各向异性本构模型、各向同性本构模型等，其中根据材料的不同而变换。（2）塑性理论类本构模型：塑性理论类本构模型是以塑性流动理论为基础，代表性的模型主要有：Mises条件的模型、理想弹塑性脆性断裂模型、应变或工作硬化塑性理论模型等。模型中考虑了混凝土加载路径和混凝土的硬化，在混凝土的应力-应变全曲线中，有上升段和下降段。自从Drucke公设和

7、Пyushin公设出现之后，经典塑性力学得到飞速发展，混凝土塑性力学模型也是基于这些公设建立的。以塑性理论为基础的混凝土本构模型，在对其加载面，包括初始屈服面，后续加载面和破坏包络面等特征面的研究中，这些特征面若以应力空间来表示时，当应力达到屈服后，材料发生应力松弛；若以应变空间表示时，当应变达到松弛面后，材料发生应变松弛。基于应力状态屈服面或破坏包络面的塑性理论类型的本构模型有弹性-全塑性模型、线弹性-硬化塑性-断裂模型等；基于松弛面的塑性理论类型的本构模型有塑性模型、塑性断裂模型、硬化断裂模型等。所有这些模型所