混凝土温度损伤的一种新扩展有限元模拟

1、混凝土细观温度损伤的一种新型扩展有限元模拟* 周储伟1，孙玉林，兴辰（机械结构强度与控制国家重点实验室，南京航空航天大学结构强度研究所，江苏南京210016）摘要：提出了一种单元子划分和子结构相结合的新型扩展有限元方法，该方法允许裂纹进入或穿过一个单元，因此裂纹的萌生和扩展可以不受有限元初始网格的限制，被裂纹进入或贯穿的单元依据裂纹几何进行子划分，覆盖一条裂纹的所有子划分单元组成一个子结构。建立了混凝土骨料、水泥砂浆、界面三相材料的细观力学模型，用上述扩展有限元方法模拟了混凝土温度荷载下微裂缝的萌生和发展，以及后继加载下的破坏情况，发现混凝土中温度损伤从弱界面开始，温度微裂缝是后继加载破坏的损

2、伤源。关键词： 混凝土；温度应力；细观力学模型；扩展有限元SIMULATION OF CONCRETE MICRO-SCALE THERMAL DAMAGE BY ANEW EXTENDED FINITE ELEMENT METHODZhou Chuwei, Sun Yulin, Xing Chen(State Key Laboratory of Mechanics and Control of Mechanical Structures ， Nanjing University of Aeronautics and Astronautics, Nanjing 210016)Abstract:

3、 A new extended finite element method (XFEM) has been developed which combining methods of element sub-partition and substructure. In this approach, crack can enter or split an element. Therefore, crack can nucleate and propagate free from the original FE mesh. The elements be notched or split by a

4、crack are sub-partitioned in respect to crack geometry and they group a substructure. Two micromechanical models of concrete with different volume fraction of coarse aggregate are employed to investigate initiation and propagation of micro scale crack under thermal load and consequent uni-axial comp

5、ressive load. The simulation demonstrates that week aggregate/cement interphase is vulnerable of micro crack initiation under thermal load, and this thermal damage will propagate further under consequent loads.Keyword: concrete, thermal stress, micromechanical model, extended finite element method1

6、一种新的扩展有限元方法细观有限元模型可以细致地考虑各种因素对混凝土破损的影响，模拟混凝土逐渐破坏的全过程1,2 ，因此为很多研究者采用。扩展有限元( XFEM ) 3,4是目前模拟裂纹扩展等动态不连续问题最有效的数值方法，但因需要构造裂尖应力奇异附加项，目前尚不能很好地模拟弹塑性等非线性问题。本文基于单元子划分和子结构的思想，提出一种新型的XFEM 方法 5,6，可以方便地推广到各种非线性问题。本文的方法将结构中的所有单元分为三类（图1）：裂纹穿过的单元（浅灰色）、包含裂尖的单元（深灰色）和普通单元（白色），前两种单元依照裂纹几何形态进行子划分（图1中虚线），第二种单元中在裂尖周围子划分一圈应

7、力奇异单元。包含一根裂纹的前两类单元组成一个子结构，子结构的内部结点（图1中实心的小圆点）的自由度，可以凝聚到子结构的边界结点（图1中加粗圆圈）的自由度，后者正是结构的初始网格点。子结构的平衡方程为：Ksst usstKout Kout-inuoutfout(1)Kin-out KinuinfinKsst 、 usst 为子结构的刚度矩阵和总结点位移，其中Kout 、 Kin 、 Kout-inKinT-out 分别为刚度阵中与边界结点位移项、内部结点位移项、边界与内部结点位移的交叉项相对应的部分，uout 、 uin 分别为边界结点基金项目：国家自然科学基金（11272147，1077207

8、8），江苏省高校优势学科建设工程项目资助作者简介： *周储伟（ 1964），男，江苏人，教授，博士，主要从事复合材料力学、结构强度、计算力学、细观力学等研究孙玉林（ 1990），女，江苏人，硕士研究生，从事计算力学研究；兴辰（1990），男，河南人，硕士研究生，从事计算力学研究。1和内部结点位移，fout 、 fin 分别是边界结点和内部结点荷载矢量。当裂纹表面不受荷载时有fin 0 ，由式 (1) 可以用边界结点自由度表达内部结点自由度：uin1(2)Kin Kin-out uout将内部结点自由度凝聚到边界结点自由度之后，子结构的刚度矩阵为：KoutKin1 Kin-outuout uou

9、tfout(3)图1 包含裂纹结构中的三类单元起裂采用最大应力（应变）准则，一旦某普通单元中的平均主应力达到了材料强度，就假设萌生一条穿过单元形心且与主应力方向垂直的裂纹。扩展采用应力强度因子准则， 应力强度因子由交叉积分求出3 ：KI KIKII KIIEijuiijuiW1 jq(x1, x2 )dA(4)2 Axxx11jWij ijijij(5)其中K I 、K II 为待求的应力强度因子，、为裂尖附近的应力和位移，、ijui可由有限元计算得到， K I K IIij 、 ui 分别为辅助场的应力强度因子及裂尖附近的应力和位移。裂纹扩展采用简化的工程判据：f2 (K I , K II

10、)K IK IIK IC1(6)其中 K IC 为材料的断裂韧性。裂纹扩展方向采用最大周向拉应力准则5。2 混凝土温度损伤的细观力学模拟当混凝土处于变温环境时，除了承受外荷载的作用，还在很大程度上受到温度应力的影响。大体积混凝土内外的温度场梯度可能导致混凝土开裂。此外，混凝土从细观角度上看是非均质材料，其中的骨料、水泥砂浆的热胀（冷缩）系数不一致，这将导致变温时存在细观的热失配应力并可能造成微裂缝损伤。图2为混凝土二维细观力学模型（粗骨料体积含量分别为35%、 51%），模型将混凝土细观结构简化为三相材料：粗骨料（浅灰色）、水泥砂浆（白色）、两者之间的界面（深灰色），这里的水泥砂浆实际上合并了

11、水泥砂浆和其他级配的细骨料。三相材料的物理、力学性能见表1。对两种模型均施加两组荷载，第一组先降温 50o C ，然后施加竖向单向压缩荷载（工况1）；第二组先升温50o C ，再竖向压缩（工况2）。图 3、图 4分别为两种模型在不同工况、不同阶段的细观裂缝扩展情况。表 1 混凝土中各相材料的物理、力学参数材料弹模 / GPa泊松比抗拉强度 / MPa抗压强度 / MPa断裂能 / N·m热胀系数-5 o10/C骨料600.166.048.02800.7砂浆250.223.530.01701.0界面300.222.520.01201.32(a)粗骨料体积含量35%(b)粗骨料体积含量

12、51%图 2混凝土细观有限元模型(a) 降温 25o C(b) 降温 50o C ，y15MPa(c) 降温 50o C ，y30MPa(e) 升温 30o C(f) 升温 50o C ，y15MPa(g) 升温 50o C ，y30MPa图 3 粗骨料体积含量35% 的模型各工况、各阶段的细观裂缝扩展(a) 降温 30o C(b) 降温 50o C ，y15MPa(c) 降温 50o C ，y34MPa3(e) 升温 30o C(f) 升温 50o C ，y15MPa(g) 升温 50o C ，y33MPa图 4 粗骨料体积含量51% 的模型各工况、各阶段的细观裂缝扩展从图 4、5 中可以发

13、现，变温时微裂缝起始于界面相，对于降温的情况更明显，因为本文模型界面相的强度比水泥砂浆更弱，而降温时界面相受拉应力。从图中还可以发现，后继单轴压缩荷载下，除了会产生新的裂缝外，原来的温度微裂缝将继续发展，并且最终的贯穿主裂缝大都发源于温度裂缝，说明温度引起的微裂缝损伤对混凝土的最终强度有明显的影响。图 5 为两种模型后继单轴压缩荷载下的平均应力应变曲线，从中可以看出，因降温是界面相和水泥砂浆受拉应力为主，比升同样的温度情况下的微裂缝损伤程度更严重一些，所以后继因裂缝扩展造成的模量下降也更明显一些。图 5两种模型 ( 温度加载 ) 后继单轴压缩的平均应力- 应变曲线参考文献1马怀发 , 陈厚群

14、, 黎保琨 . 混凝土细观力学研究进展及评述J. 中国水利水电科学研究院学报 , 2004, 2(2): 124 130.2陈爱玖 , 章青 , 刘仲秋 , 夏晓舟 . 混凝土细观分析的研究现状J. 混凝土 , 2008, 6(224): 21 25.3 T. Belytschko, T. Black. Elastic crack growth in finite elements with minimal remeshing J. International Journal for Numerical Methods in Engineering. 1999, 45: 601 620.4 Nicolas, Dolbow, Belytschko. A finite element method for crack growth withou