基于细观结构混凝土损伤力学研究.docVIP

基于细观结构混凝土损伤力学研究 　　摘要：本文对混凝土的微观结构进行数值模拟，利用有限元软件ansys分析， 用构型力和守恒积分值来反映混凝土试件在加载过程中的损伤程度，并通过混凝土加载试件的CT实验结果，在细观层次上对比缺陷扩展、裂纹萌生等损伤现象，通过实验手段对这一混凝土新损伤理论进行验证。 　　关键词：混凝土，构型力，损伤力学 　　引言：作为主导建筑材料的混凝土，限于其强非均质性（骨料、水泥砂浆及其交界面组成的三相复合材料），如何从细观层次出发，正确预测和描述内部微裂纹、夹渣、气泡、空穴等多缺陷的损伤演化规律，仍存在着一些关键的科学和技术问题亟需解决。立足于其细观特性，本研究将发展一种新的适用于混凝土的细观损伤力学模型。借助于混凝土的微观结构，将细观力学和损伤力学相结合，并首次通过能量观点引入混凝土材料的构型力概念，为研究混凝土材料缺陷构型变化及其破坏损伤的微观机理，提供了一个极好的内变量；选取围绕缺陷的闭合曲线作为积分路径，对构型力进行数值积分，计算其守恒积分，作为衡量混凝土损伤水平的外变量。本研究将有助于完善当前混凝土细观损伤力学理论，研究成果将为混凝土材料损伤理论研究和工程应用提供重要的理论依据，对推动基于细观结构的混凝土损伤力学发展具有重要的意义。 　　1 常用混凝土细观数值模拟模型 　　格构模型 　　格构模型是典型的细观数值模型，已被使用了50多年，最初它被用来求解经典的弹性力学问题，但仅仅只停留在理论部分。格构模型从80年代后期开始，人们才用它来模拟非均匀材料的破坏过程[1]。此模型的原理是：在细观尺寸上，将连续介质离散成由梁单元或弹性杆件组合而成的格构系统，使用的网格是规则的三角形单元，单元由梁或杆组成，杆单元仅仅传递轴力，梁单元可以传递剪切力、轴力、弯矩，所以梁单元可以反应比较复杂的受力状态Van Mier[1]应用三维格构模型分别模拟了单轴拉伸、压缩试验以及联合拉剪试验。杨强教授以及他的团队[2]采用格构模型模拟了岩石类材料的开裂、破坏过程。Schlangen等人[3]应用格构模型（采用梁单元）模拟了一些非均质材料包括混凝土材料的破坏原则和裂纹贯通的过程。 　　2 有限元网格的模型与加载 　　已经建立的混凝土数值模型中，骨料和砂浆之间的粘结面单元划分采用多次网格划分技术，达到最大限度地降低计算规模，但是单元数量之多，计算量之巨大，给数值模拟带来了很大的困难。本文在骨料和砂浆单元之间加入面-面接触单元，模拟骨料和砂浆之间的粘结界面，以避免在计算中模拟的混凝土破坏区域出现本不存在的空洞，同时也避免事先设置接触面单元的盲目性。 　　2.1 Ansys的接触能力 　　（1）点-点接触单元。 　　点-点接触单元主要被用于模拟点和点之间的接触问题。必须要先知道接触的位置，点-点接触只适用于接触面之间相对滑动很小的情况（也包括几何非线性）。 　　（2）点-面接触单元 　　点-面接触单元是用来模拟点-面的接触行为。如果接触面是通过一组节点来生成多个单元，则可以用点-面的接触单元来给面-面的接触行为建模，接触面可以是柔体或者刚体。 　　（3）面-面接触单元 　　面-面接触单元主要给面和面的接触建模。分为：刚-柔和柔-柔接触。在ANSYS中，目标单元和接触单元组成一个接触对，程序通过一个共享的实常数号来识别接触对，为了建立接触对需要给目标单元和接触单元指定相同的实常数号。 　　2.2 有限元网格的模型生成步骤 　　由于细观混凝土三相材料之间性质的各不相同，对砂浆和骨料单独进行网格划分，在骨料和砂浆之间加入接触单元模拟粘结面。本文采用的是Matlab7.1，AutoCAD2007，ANSYS12.1软件来实现的，具体步骤如下：首先使用Matlab实现骨料随机几何模型，即编程实现骨料混凝土试件的数据、骨料轮廓线的数据。生成扩展名为.txt的文件，即AutoCAD可以识别的格式。然后，用AutoCAD进行前处理即：作面域，通过差集、分解等处理方法，将骨料面和砂浆面分别做成2个图形，都输出为扩展名为.sat的文件；最后把扩展名为.sat的2个图形文件，导入进ANSYS的界面中，分别导入两个文件，需要注意的是导入文件时需要区别开区域编号。此时，ANSYS里的图形存在2种面，一种是颗粒区域面，一种砂浆区域面 （只有一个），然后，按编号分别选元素种类、弹性模量、泊松比和划分的尺寸等，进行自由网格划分。再用接触单元来模拟砂浆与骨料之间的粘结面。 　　计算中，将均布荷载施加在混凝土试件的上边缘，对试件进行4步加载，每步的荷载值如下：5000N，第二步：5500N，第三步：6000N，第四步：6500N。 3 引用构型力概念对混凝土的损伤演化模型进行描述 　　3.1 构型力与守恒积分理论及数值计算方法