基于数值模拟混凝土系列有效断裂判据_邓爱民.pdf

基于数值模拟的混凝土系列有效断裂判据* 1 1 邓爱民 , 徐道远 （1. 河海大学力学与材料学院， 南京 210098 ） 摘 要：基于ANSYS平台，对某系列混凝土断裂破坏试验的加载过程进行有限元数值模拟，获得 了荷载-加力点位移曲线。利用获得的曲线，在线性渐进叠加假定的基础上，确定了试件在失稳断 裂时的临界有效裂缝长度。将模拟峰值荷载和临界有效裂缝长度代入由LEFM推导出的公式计算了 考虑断裂过程区影响的有效断裂韧度。最后，通过对系列试验有效断裂韧度的分析，建立了基于 数值模拟的混凝土系列有效断裂判据。 关键词：数值模拟；混凝土断裂；有效断裂韧度；断裂判据 中图分类号：TU37-02 1 引 言 在混凝土断裂方面的研究中发现，混凝土与金属材料的断裂特性不同，混凝土结构受力后， 在宏观裂缝的尖端部位会出现一个微裂纹区，称之为断裂过程区，结构失稳断裂时，由于断裂过 程区的发育，宏观裂缝已经有了一定的扩展，根据峰值荷载与初始裂缝长度计算断裂韧度显然是 不合适的，因此在混凝土断裂判据建立过程中，断裂韧度的计算必须考虑断裂过程区的影响[1，2] 。 随着计算机技术和数值模拟方法的进步，对混凝土材料的断裂破坏过程进行有限元数值模拟已成 为可能。在数值模拟结构与实际结构具有足够近似物理条件下，数值模拟比真实试验省时、省力、 无仪器限制，可以重复、大批量进行，能够模拟许多加载设备难以实现的复杂加载路径，而且能 够不受试验中各种因素对于结果的影响。对混凝土结构而言，根据数值模拟实验获得的荷载-位移 曲线也可以反映微裂纹的出现、发展等行为即断裂过程区的发育情况，因此由该曲线可确定对应 于峰值荷载的临界有效裂缝长度，并根据模拟峰值荷载与临界有效裂缝长度，计算考虑断裂过程 区影响的有效断裂韧度，从而建立考虑断裂过程区影响的有效断裂判据[3,4,5] 。 近年来，ANSYS 作为大型通用有限元软件，被大量应用于各类型结构的数值分析计算。考虑 到混凝土材料的非线性特性，ANSYS 中引入了 SOLID65 单元，该单元允许发生塑性变形及三个方 向的开裂及压碎，因此常用来模拟混凝土结构的断裂破坏过程。该单元的单元属性本身包括了受 拉区裂缝的模拟和用来说明受压区混凝土的压碎概率的塑性算法。每个单元有 8 个积分点，在这 些点处完成裂缝和压碎的检查。在没有达到混凝土抗拉强度或者抗压强度之前，单元表现为线性。 一旦单元主应力之一在积分点超过了混凝土抗拉强度或者抗压强度，单元裂缝或者压碎开始出现。 随着应力在局部的重分布，在垂直于相应主应力方向形成裂缝区或者压碎区。由于这样的单元是 非线性的，软件使用了迭代求解器[6,7,8] 。 2 系列混凝土断裂破坏试验的数值模拟 本文基于 ANSYS 软件平台，采用上述介绍的 SOLID65 混凝土单元，对某系列混凝土断裂破 坏试验的加载过程进行了有限元数值模拟。数值模型参照文献[9]，分为两大类，即梁式试件和柱 式试件。试件初始裂缝长度 a0 均为 50mm，采用自底而上的建模方式，先创建关键点，由关键点 连成面，在面上先剖分成平面等参 4 节点单元（PLANE42 ），再拉伸成 8 节点空间等参单元 SOLID65。 梁式试件（图 1.a）用于张拉Ⅰ型、Ⅱ型、Ⅲ型、Ⅱ-Ⅲ复合型及各种含张拉Ⅰ型成分的复合 型裂缝的模拟。柱式试件（图 1.b）用于压缩Ⅰ型和各种含压缩Ⅰ型成分的复合型裂缝的模拟。模 国家自然科学基金资助项目） 作者简介：邓爱民（1977-），男，汉族，湖南永州人，讲师、博士，主要从事工程结构的断裂与损伤方面研究． 拟峰值荷载参考选取文献[9]中结果，并分成 5-6个荷载步施加，每个荷载步内由程序自动分为若 干子步。根据数值模拟结果，绘出的典型荷载-加力点位移曲线如图 2 所示。 由模拟得到的荷载-加力点位移曲线可以看出，由于混凝土材料在加载过程中的非线性特性， 曲线在上升段的后阶段斜率已经有了一定的降低，说明在混凝土试件失稳断裂前，试件的柔度逐 渐降低，也就是说，在失稳断裂时试件的初始裂缝已经有了一定的扩展。若仍根据峰值荷载和初 始裂缝长度计算其应力强度因子显然是不合适的，因此需确定对应最大荷载时的裂缝长度即临界 有效裂缝长度