应力空间下混凝土动态劈拉损伤演化方程及损伤界点的确定.pdf

第 3 3 卷 第 8期 2 0 1 6年 8月 长江科学院 院报 J o u r n a l o f Ya n g t z e R i v e r S c i e n t i fi c Re s e a r c h I n s t i t u t e V0 1 33 No 8 A u g 2 0 1 6 d o i ： 1 0 1 1 9 8 8 c k y y b 2 0 1 5 0 4 6 2 2 0 1 6 ， 3 3 ( O 8 ) ： 1 2 5 1 2 9 应力空问下混凝土动态劈拉损伤演化方程 及损伤界点的确定 吴彬 ， 胡伟华 ， 唐克静 ( 1 湖北工业职业技术学院 建筑工程系， 湖北 十堰4 4 2 0 0 0 ； 2 十堰市建筑设计研究院， 湖北 十堰442 0 0 0 ； 3 三峡地区地质灾害与生态环境湖北省协同创新中心 ， 湖北 宜昌443 0 0 2 ) 摘要： 为研究混凝土动态劈拉特性 ， 进行了不同应变速率下的混凝土动态劈拉试验， 建立劈拉应力水平与声发射 能量数的数学关系式， 明确了混凝土劈拉损伤演化方程， 提出了损伤界点的概念。基于医学 R O C曲线确定最佳临 界点理论 ， 分析出了混凝土劈拉破坏损伤界点所对应的应力水平。研究结果表明： 高应变速率下， 损伤第 1阶段声 发射信号开始产生时的应力滞后， 当应力水平超过一定数值后 ， 损伤第 2段的声发射现象与低应变速率时的相似 ； 损伤界点随着应变速率的提高向应力水平较高处推移， 其变化范围为 0 7 1 ， 0 7 7 ， 且推移的规律服从线性分布， 应 力空间下混凝土劈拉损伤界点为0 7 5 。 关键词 ： 混凝土劈拉试验； 声发射； 损伤演化 ； R O C曲线； 损伤界点 中图分 类号 ： T U 5 0 2 6 文献标志码 ： A 文章编号 ： 1 0 0 1 5 4 8 5 ( 2 0 1 6 ) 0 8 0 1 2 5 0 5 1 问题 的提出 自从 D o u j i l l l 1 将损伤力学应用到混凝土材料研 究领域 ， 国内外学者对混凝土材料的损伤研究一直 没有间断过 。关于混凝土损伤的研究 ， 遵循 了从单 轴应力损伤到多轴应力损伤、 从弹性损伤到弹塑性 损伤、 从静力损伤到动力损伤 、 从各 向同性损伤到各 向异性损伤 的思路 。D u b e等 从与应变速率无 关的损伤模 型 中推导 出与加载 速率相关 的损伤模 型 ， 并建立了考虑 P e r z y n黏塑性影 响的损伤本构方 程。E i b l 等 在试验结果基础上 ， 对混凝土损伤特 性进行了分析 ， 认为应变速率改变了混凝土在加载 过程 中的损 伤路径。 肖诗云等 构建 以切线模 量 为参量的损伤变量 ， 分析 了混凝 土动 态损伤 门槛值 随应变率的变化规律和混凝土 的动态损伤演 化规 律。李庆斌等 从 混凝土动力 损伤与静力损 伤之 间关系的角度 ， 提 出了考虑初始 弹性模量变化的混 凝土动力损 伤本构模型 的建 立方法 。李杰等 采 用受拉损伤变量和受剪损伤变量建立了基于能量的 弹塑性损伤本构模 型， 并进行 了数值计算和试验验 证。王春来等 根据 We i b u l l 统计分 布理论和等效 应变假定原理 ， 建立 了单轴受压状态下钢纤维混 凝 土损伤本构模型。 目前对于混凝土材料损伤的研究主要集中在损 伤特l生与损伤模型方面， 但是对于混凝土损伤演化规 律鲜有报道 ， 而对于损伤阶段的划分和损伤界点的确 定几乎无相关的研究成果。已有研究成果阐明混凝 土单轴压缩损伤过程损伤存在着明显的 3 个阶段， 即 损伤萌芽阶段、 损伤稳定发展 阶段、 损伤不稳定发展 阶段 。损伤阶段所 占整个损伤过程 的比例对于研 究混凝土损伤演化过程的规律具有重大意义。 因此， 无论是混凝土单轴压缩还是劈拉试验， 以 何种方式来 明确划分损伤阶段显得至关重要。鉴于 此， 进行了 4种应变速率( 1 0 ， 1 O 一， 1 0 一， 1 0 S ) 下 混凝土动态劈拉试验 ， 对劈拉过程中的混凝土声发射 信号进行实时采集， 通过确定劈拉应力水平与声发射 能量数之间的关系 ， 分析了混凝土劈拉损伤演化规 律 ， 并提出了损伤界点这一新概念 ， 即损伤界点为前 一 损伤阶段的终点或者后一阶段的起点。 2试验概况 2 1 试验设备 试验采用的设备 为 1 0 MN微机控制 电液伺 服 大型多功能动静力三轴仪。该设备可进行多种试件 收稿 日期 ： 2 0 1 5 0 6 0 2 ； 修回 日期 2 0 1 5 0 8 1 2 作者简介 ： 吴彬 ( 1 9 8 9 一 ) ， 女 ， 湖北 十堰人 ， 硕士研究生 ， 主要从事混凝土材料研究 ， ( 电话 ) 1 8 7 7 2 8 4 6 3 7 0( 电子信箱)1 6 6 0 6 2 1 6 2 8 q q C O IX I 。 通讯作者： 胡伟华 ( 1 9 8 8 - ) ， 男， 湖南益阳人， 硕士研究生， 主要从事结构抗震设计和混凝土材料研究 ， ( 电话) 1 8 7 7 2 8 0 3 7 4 9( 电子信箱) 4 64 9 1 3 98 8 q q t o m。 长江科学院院报 2 0 1 6年 尺寸大小 的混凝土动静力加载试验 ， 具有多功能性。 控制器最大采样频率为0 0 0 0 1 S ， 系统最大数据存 储量为 3 O万组样本 ， 可实时 自动采集数据及记录当 前情况并反馈相关信息， 实现存储、 处理试验数据同 步进行 。声发射信号采集设备为 S A E U 2 S系统 。 2 2 试件的制作与加工 2 2 1 试 件材 料 的选 用 试验采用4 2 5 普通硅酸 盐水泥 ， 搅拌混 凝土所 用的水为饮用 自来水。 2 2 2 配合比的确定及试件 的养护与加工 试件采用 1 5 0 m mX 1 5 0 mm圆柱体试件 。配合 比根据 普通混凝土配合 比设计规程 ( J G J 5 5 2 0 1 1 ) 进行 理 论计 算 与 试 配 ， 最 终 确 定 配合 比为 1： 0 5 6： 2 5 5：3 8 3( 水 泥 ：水 ：粗骨 料 ：细骨 料) 。混凝土拌合料坍落度值在 2 0 4 0 m m之间试 件成型后 4 8 h 脱模， 标准养护。 2 3 试验步骤 ( 1 )装样 。确保试件的中心线和传力柱严格对 中。 ( 2 )声发射探头安装 。通过耦合剂将声发射探 头贴在试件平整的部位。 ( 3 )声发射 参数设 置及 检查。参 数设 置完毕 后， 进行声发射信号的检查 ， 检查各个通道是否正常 接收信号。针对本试验， 声发射参数设置内容见 表 1 。 表 1 声发射仪 器的参数设置 T a b l e 1 P a r a me t e r s e t t i n g o f AE( a c o u s ti c e mi s s i o n) i ns t r u m e nt ( 4 ) 变形计的安装与检查 。 ( 5 ) 试验加载。针对 4种不同工况预先设好竖 向加载的速率 ， 在点击发送指令 的同时开始采集声 发射信号， 加载到试件破坏为止。 3基 于声 发 射 能 量 的劈 拉 损 伤 演 化 分析 3 1 试验结果分析 试验得到的声发射撞击数与应力水平的对应关 系如 图 1 。其中 。为峰值应力 ， 田 为能量数最大 值 ， 为能量数 瞬时值 。所谓应力水平是 指应 力与 峰值应力 的比值 ， 即 V=c r c r 。 。 l 0 。 。 0 0 ( a 】 应变速率为l ff S s 1 0 。 害 。 O O 1 O 2 0以3 O 4 0 5 O ( b 】 应变速率为1 0 4 s O l 2 t s 3 4 u u l I S u z u j ( c ) 应变速率为1 0 s( d 】 应变速率为1 0 z s _-能量水平 应力水平 图 1 不 同应 变速率下的应力水平 与声 发射 能量水平的比较 Fi g 1 Re l a t i o n s h i p b e t we e n s t r e s s l e v e l a n d e n e r g y n u m be r o f AE u n d e r d i ffe r e n t s t r a i n r a t e s 由图 1 可知 ， 声发射数与混凝土应力水平 存在 着一种更为直接的关系， 现对混凝土劈拉试验的声发 射能量水平 n与相对应力水平 的关系进行进一步分 析。声发射能量水平 n是指从初始状态到某一应力 水平下声发射能量总数与应力水平为 1时声发射能 量总数之 比， 即 n = N 眦。不同应变速率下的声发 射能量水平的具体数值详见表 2 ， 声发射相对能量水 平与相对应力水平的关系如图 2所示 。 3 2 试验结果拟合 从 图 2中曲线的走势可 以明显地看出 ， 应力水 表 2不 同应变速 率下声 发射能量水平统计 Ta b l e 2 Re l a t i v e e n e r g y n u mbe r s o f AE u n d e r d i ffe r e n t s t r a i n r a t es 第8期 吴彬 等 应力空间下混凝土动态劈拉损伤演化方程及损伤界点的确定 1 1 0 O 0 0 2 0 4 ， 0 6 0 8 1 0 an ( a ) 应变速率为1 0 s 1 2 1 0 Q 0 8 o 6 0 4 0 2 O o |d 6 | 。 a ( b ) 3 1 0 N s ( a ) 应变速率为l 0 。 V s ( b ) 应变速率为l 0 4 s a。 d。 q ( c ) 应变速率为l 0 。 V s ( d ) 应变速率为1 0 。 s 图 2不同应变速率下能量水平 与应 力水 平的关 系 F i g 2 Re l a t i o n s h i p b e t we e n e n e r g y l e v e l a n d s t r e s s l e v e l u n d e r d i ff e r e n t s t r a i n r a t e s 平 0 8 时， 声发射相对能量总数或者 损伤值急剧上升。混凝土劈拉过程中混凝土损伤演 化规律存在着明显 的 2阶段 ， 即损伤发展 阶段 ( 损 伤第 1阶段) 、 损伤不稳定阶段( 损伤第 2阶段 ) 。 文献 9 假定应力水平从 到 V + d V的过程 中 材料所产生的声发射数 的概率 密度 函数为 ) ， 则 厂 ( 满足如下关系 ： d V=d N Nm =凡。 ( 1 ) 文献 1 0 进一步假定密度函数 厂 ( 具有如下 形式 ： 厂 ( =a V+b 。 ( 2 ) 根据式( 1 ) 和式 ( 2 ) 可得到声发射水平 n与应 力水平 之 间的关系为 乃=c We x p ( b v ) 。 ( 3 ) 式中 a ， b ， C 为试验参数。 文献 1 1 阐明混凝土损伤可通过声发射能量数 来表征并定义损伤变量为 D= N Nm ， 结合本文给 出的能量水平定义发现 D= n， 故损伤的表达式为 D( = =c We x p ( 6 。 ( 4 ) 现利用式 ( 4 ) 对劈拉试验声发射损伤与相对应 力水平进行拟合。拟合参数详见表 3 ， 拟合效 果如 图 3所 示 。 表 3 拟合参数值 Ta b l e 3 Fi t t i n g p a r a m e t e r s i n s p l i t tin g t e s t 座变逵壅 ： ! 壁 1 0 00 3 7 4 1 l 1 3 2 51 0 9 9 1 4 1 00 01 9 1 0 5 5 3 9 6 8 0 9 9 9 1 1 00 0 5 6 3 6 4 8 2 8 8 8 0 9 8 4 9 1 0 0 0 3 8 4 1 2 9 3 2 7 7 0 9 9 1 4 1 2 1 0 0 4 O 2 0 1 2 1 0 Q 0 8 o 6 0 4 O 2 O 口。q d | On ( c ) 应变速率为1 0 。 s ( d ) 应变速率为1 0 2 s -试验值 拟合 曲线 图3 不同应变速率下声发射损伤与应力水平拟合关系 F i g 3 F i t tin g r e l a ti o n s h i p b e t we e n AE d ama g e an d s t r e s s l e v e l u n d e r d i f f e r e n t s t r a i n r a tes 从拟合 的效果来看 ， 式 ( 4 ) 能较好地反 映混凝 土在劈拉过程 中的损伤演化规律 。因此 ， 式 ( 4 ) 可 用来进一步确定劈拉损伤界点 。 4 基于 R O C曲线的劈拉损伤界点确定 R O C( r e c e i v e r o p e r a t i n g c h a r a c t e r i s t i c ) 曲线是医 学上将诊断试验结果 划分为若干临界点 ， 以每个临 界点对应 的灵敏度为纵坐标 ， 特异度为横坐标 ， 作图 得到的曲线 ， 是一种全面 、 准确评价诊断试验的有效 工具H 。R O C曲线的另一个作用是确定检测的最 佳阈值 。 R O C曲线法确定临界点多数情况下， 选 择 曲线上尽 量靠近左上方 的点来 确定临界点 为最 佳 。 本文将医学确定 临界值 的 R O C曲线方法运用 到损伤界点的确定过程中来 ， 规定应力水平对应灵 敏度 ， 损伤对应特异度 ， 即对劈拉损伤界点的确定转 化为损伤演化曲线上的点到右下方点 ( 1 ， O ) 的距离 最短的点。 损 伤 演 化 曲线 上 任 意 一 点 ( ，D( )到 ( 1 ， O ) 的距离的平方值为 F ( =( 一1 ) +D ( 。 ( 5 ) 对式( 5 ) 求一阶导数为 F ( =2 一1+ c We x p ( b y )】 ( a V+b ) 。 ( 6 ) F ( 的最小值 点必 然在 F ( =0处产 生 ， 故 由 F ( =0得到 一 1+ c V a e x p ( b Y )】 ( a V+b )=0 。 ( 7 ) 长 江科 学院 院报 2 0 1 6年 上述方程为超越 方程 ， 用 目前的求根方法无法 获得其精确解 ， 只能对其进行数值分析运用改进的 区间搜索法即二分法对其进行近似求解 。 本文规定精度为 =0 0 0 1 ， 由式 ( 7 ) 可知 区间 搜索 的次数 k=7 ， 相关数据 如表 4所示 ， 损伤界 点 应力水平与应变速率的变化关系及拟合曲线如图 4 所 示 。 表 4 区间搜索相关参数 T a b l e 4 Re l a t e d p a r a me ter s f o r i n t e r V a l s e a r c h 图 4损伤 界点应力水平与应变速 率的变化关 系 及拟合 曲线 F i g 4 Re l a t i o n s h i p b e t we e n s t r e s s l e v e l a t c r i t i c a l d a ma g e p o i n t a n d s t r a i n r a t e 由表 4及图 4 ( a ) 易知 ， 随着应变 速率 的提高 ， 混凝土劈拉损伤界点 向着应力水平较 高处推移 ， 界 点的变化引起 了损伤 2个 阶段 应力水 平 区间的变 化。损伤界点 的变化范 围为0 7 1 0 7 7 。由不 同应 变速率下损伤界点应力水平变化图的走势可知， 损 伤界点应力水平与应变速率呈现较为明显 的线性关 系 ， 故可对表 4中的数据进行线性 回归分析 ， 发现损 伤界点应力水平与应变速率的关系均满足关系式： y = A x + B， 其中 = l g ， Y= t r t r 0 ， 相应的拟合效果 图 与拟合度详见 图 4 ( b ) 。损伤界点变化 幅度范 围为 2 3 8 一 8 6 8 ， 界点极差 为8 7 8 ， 小 于 3 0 ， 因此 应力空间下劈拉损伤界点可取0 7 5 。 5 结 论 本文对劈拉过程混凝土材料所产生的声发射信 号进行分析 ， 建立 了劈拉损伤值与应力水平数学关 系 ， 确定 了混凝土劈拉损伤演化规律演化方程 ， 提出 了应力空间条件下混凝土劈拉损伤界点 ， 得 出的主 要结论如下 ： ( 1 )在高应变速率下 ， 声发射发展 阶段 明显滞 后 ， 当应力水平超过一定值后 的声发射现象与低应 变速率时的规律相似。 ( 2 )试件劈 拉破 坏过程 中， 当相对应 力水 平 O 8 时， 损伤值急剧上升 ，说明损伤界点 的应力水平在 0 8 附近波动。 ( 3 )损伤界点随应变速率 的变化范 围为0 7 1 0 7 7， 变化的幅度范 围为2 3 8 8 6 8 ， 应力空间条 件下的劈拉损伤界点可取0 7 5 。 参考文献 ： 1 D O U G I L L J W O n S t a b l e P r o g r e s s i v e l y F r a c t u ri n g S o l i d s J J o u r n a l o f A p p l i e d Ma t h e ma t i c s a n d P h y s i c s ，1 9 7 6 ， 2 7 ( 4 ) ： 4 2 3 - 4 3 7 2 封伯昊， 张立翔， 李桂青 混凝土损伤研究综述 J 昆 明理工大学学报(自然科学版) ， 2 0 0 1 ，2 6( 3 ) ： 2 1 3 0 3 D U B E J F ，P I J A U D I E R C A B O T G，B O R D E R I E C L Ra t e De p e n d e n t Da ma g e Mo d e l f o r Co n c r e t e i n Dy n a mi c s J J o u rna l o f E n g i n e e r i n g Me c h a n i c s ， 1 9 9 6 ， 1 2 2 ( 1 0 ) ： 9 39-9 47 4 E I B L J ，S C H MI D T H U R T I E N N E B S t r a i n r a t e s e n s i t i v e C o n s ti t u ti v e L a w f o r C o n p c r e t e J J o u r n a l f o r E n g i n e e r - i n g M e c h ani c s ， 1 9 9 9 ， 1 2 5 ( 1 2 ) ： 1 4 1 1 1 4 2 0 5 肖诗云， 张剑 不同应变率下混凝土受压损伤试验研 究 J 土木工程学报 ， 2 0 1 0 ， 4 3 ( 3 ) ： 4 0 - 4 5 6 李庆斌， 邓宗才， 张立翔 考虑初始弹模变化的混凝土 动力损伤本构模型 J 清华大学学报( 自然科学版) ， 2 0 0 3 ， 4 3 ( 8 ) ： 1 0 8 8 1 0 9 1 7 李杰， 吴建营 混凝土弹塑性损伤本构模型研究 I： 基本公式 J 土木工程学报， 2 0 0 5 ，3 8 ( 9 ) ：1 5 - 2 0 8 王春来 ， 徐必根 ， 李庶林 ， 等 单轴受压状态下钢纤 维混凝土损伤本构模型研究 J 岩土力学 ，2 0 0 6 ， 2 7 ( 1 ) ：1 5 1 - 1 5 4 9 纪洪广， 张天森， 蔡美峰， 等 混凝土材料损伤的声发射 动态检测试验研究 J 岩石力学与工程学报 ， 2 0 0 0 ， l 9 ( 2 ) ： 1 6 5 - 6 8 1 O 1 MA T S U Y A M A K， I S H I B A S H I A， O H T S O M R a t e P r o c e s s An a l y s i s o f AE Ac ti v i t y i n Co r e t e s t o f De t e ri o r a t e d C o n - c r e t e C J a p a n e s e S o c i e t y f o r N o n d e s t r u c t i v e I n s p e c - t i o n Proc e e d i n g s o f t h e l l th I n t e r n a t i o n al Ac o u s t i c E mi s s i o n S y mp o s i u m， F u k u o k a ， J a p an ， O c t o b e r 2 6 2 9， 1 9 92： 1 1 2-11 5 1 1 张明， 李仲奎 ， 杨强， 等 准脆性材料声发射的损伤 模型及统计分析 J 岩石 力学与工程 学报 ， 2 0 0 6 ， 2 5 ( 1 2 ) ： 2 4 9 3 2 5 0 1 1 2 Z WE I G M H， C A M P B E L L G R e c e i v e r o p e r a t i n g C h a r a c 第 8期 吴彬 等应力空间下混凝土动态劈拉损伤演化方程及损伤界点的确定 1 2 9 t e r i s t i c( R O C)P l o t s ： A F u n d a m e n t a l E v al u a t i o n T o o l i n C l i n i c al M e d i c i n e J C l i n i c al C h e m i s t r y ， 1 9 9 3 ， 3 9 ( 4 ) ： 5 6 1 -5 6 7 1 3 王骏， 吴虹桥， 宋兆琪 受试者作业特征曲线 ( R O C) 及其 在影 像学 中的应用 J 放射 学 实践，2 0 0 1 ， 1 5( 3 ) ： 2 0 1 2 0 2 1 4 】陈卫中， 倪宗瓒， 潘晓平， 等 用 R O C曲线确定最佳临 界点和可疑值范 围 J 现代预 防医学， 2 0 0 5 ， ( 7) ： 7 2 9 7 3 1 De t e r m i n a t i o n o f Cr i tic a l Da m a g e Po i n t a n d Ev o l Co nc r e t e s Dy n a m i c S p l i t ting Da m a g e u n d e r ut i o n Eq S t r e s s S ( 编辑 ： 姜小兰 ) u a t i o n o f p a c e wu B i n ，H U We i h u a ，T A N G K e - j i n g ( 1 D e p a r t m e n t o f A r c h i t e c t u r e ， H u b e i I n d u s t ri a l P o l y t e c h n i c ， S h i y a n 4 4 2 0 0 0， C h i n a ； 2 S h i y a n A r c h i t e c t u r a l D e s i g n I n s t i t u t e，S h i y a n 4 4 2 0 0 0，C h i n a ；3 C o l l a b o r a t i v e I n n o v a t i o n C e n t e r f o r G e o h a z a r d s a n d E c o E n v i r o n me n t i n T h r e e G o r g e s A r e a o f Hu b e i P rov i n c e ， Y i c h a n g 4 4 3 0 0 2 ， C h i n a ) Ab s t r a c t ： I n o r d e r t o e x p l o r e t h e f e a t u r e s o f c o n c r e t e d y n a mi c s p l i t t i n g， we c a r ri e d o u t s p l i t t i n g t e s t u n d e r d i f f e r e n t s t r a i n r a t e s a n d e s t a b l i s h e d t h e r e l a t i o n s h i p b e t w e e n s t r e s s l e v e l a n d A E ( a c o u s t i c e m i s s i o n )e n e r g y n u mb e r T h e n ， w e p u t for w a r d t h e c o n c e p t o f c rit i c a l d a ma g e p o i n t a n d d e t e r mi n e d t h e e v o l u t i o n e q u a t i o n o f d y n a mi c s p l i t t i n g d am - a g e o f c o n c r e t e O n t h e b a s i s o f d e t e rmi n i n g o p t i ma l c ri t i c a l p o i n t b y u s i n g m e d i c a l R O C ( r e c e i v e r o p e r a t i n g c h ara c t e r i s t i c )c u r v e ， w e g a v e t h e s t r e s s l e v e l o f c o n c r e t e a t c ri t i c a l d am a g e p o i n t i n s p l i t t i n g t e s t R e s u l t s s h o w e d t h a t 1 ) u n d e r h i g h s tr a i n r a t e ，s t r e s s l a g s w h e n e m i s s i o n s i g n al o c c u r s a t t h e fi r s t s t a g e o f d a m a g e ； 2 )A E p h e n o me n a a t t h e s e c o n d s t a g e w e r e s i m i l a r t o t h o s e w i t h l o w s t r a i n r a t e a f t e r s tr e s s l e v e l e x c e e d e d a g i v e n v al u e ； 3 )w i t h t h e i n c r e a s e o f s t r a i n r a t e ， t h e s tr e s s l e v e l a t c ri t i c a l d a m a g e p o i n t i n c r e a s e d ， a n d i t s r a n g e w a s 0 7 1 ， 0 7 7 ； 4 )t h e r e i s fl l i n e a r r e l a t i o n s h i p b e t w e e n s t r e s s l e v e l a t c ri t i c al d am a g e p o i n t a n d l o g a r i t h mi c f u n c t i o n o f s t r a i n r a t e ，an d t h e s tra i n r a t e o f c o n c r e t e a t c rit i c al d a ma g e p o i n t i n s p l i t t i n g t e s t u n d e r s t r e s s s p a c e wa s 0 7 5 K e y wo r ds： s p l i t t i ng t e s t ；a c o u s t i c e mi s s i o n； da ma g e e v o l u t i o n；ROC c u r v e；c rit i c al da ma g e p o i n t ( 上接第 1 2 4页 ) s t r e n g t h d e t e ri o r a t i o n o f c o r r o d e d c o n c r e t e t h r o u g h t e s t s o f c o r r o s i o n( b y a m mo n i u m s u l f a t e s o l u t i o n )o n c o n c r e t e s mi x e d wi t h fl y a s h，s l a g，a n d b o t h fl y a s h a n d s l a g ，r e s p e c t i v e l y F u r t h e rm o r e， we e mp l o y e d the t h e o r y o f g r a y c o r r e l a t i o n t o a n aly z e t h e c o r r e l a t i o n b e t w e e n s tre n g t h a n d mi n e r al a d mi x t u r e s o f c o n c r e t e Re s u l t s s h o w t h a t t h e c o n c r e t e mi x e d wi t h fl y a s h a n d s l a g r e p l a c i n g c e me n t i n e q u al a mo u n t c o ul d s i gn i fic a n tly i n c r e a s e t h e r e s i s t a n c e o f c o n c r e t e t o am mo n i u m s u l f a t e c o r r o s i o n；t h e s t r e n gth o f c o n c r e t e wi t h h i g h c o n t e n t o f mi n e r al a d mi x t u r e s g r o w s s l o wl y，b u t t h e r e s i s t a nc e t o a mmo n i u m s u lfa t e c o r r o s i o n i s be t t e r ；mi n e r al a dmi x t u r e s h a v e g o o d e f f e c t i n d e l a y i n g the d e c l i n e o f c o mp r e s s i v e s t r e n gth o f c o rr o d e d c o n c r e t e b u t s h o ws wo r s e e f f e c t f 0 r fl e x u r al s tr e n gth T h e r e s u l t s o f gre y a n aly s i s r e v e al t h a t c o n c r e t e s C 1( 1 0 d o s a g e o f f l y a s h ) ，C 4( 4 0 d o s a g e o f s l a g ) ， a n d C 5( 5 0 d o s a g e o f s l a g )h a v e t h e b e s t p e r f o rm anc e s i n r e s i s t i n g am m o n i u m s u l f a t e c o r r o s i o n K e y wo r ds ： c o nc r e t e；a mmo n i u m s u lfa t e；fl y a s h；s l a g；d o u b l e mi x i n g；c o r r e l a t i o n