Experimental research on the mechanical properties of fiber-reinforced autoclaved aerated concrete under cyclic loading
-
摘要:
为研究纤维增强蒸压加气混凝土(FAAC)循环受压力学行为，共设计11组棱柱体试件进行单调及循环受压试验，分析纤维种类(玄武岩纤维(BF)、碳纤维(CF))和纤维掺量对FAAC破坏形态、应力-应变全曲线特征、塑性应变、刚度退化率、应力退化率等力学性能指标的影响规律。研究结果表明：循环荷载作用下FAAC的破坏模式主要为剪切破坏和竖向劈裂破坏，随纤维掺量增加，试件破坏模式由剪切破坏转向竖向劈裂破坏；纤维掺量为0.4wt%时，FAAC的峰值应力增幅最大，BF/蒸压加气混凝土(AAC)的单调加载曲线和循环加载曲线峰值应力分别增加了24.29%、29.16%，CF/AAC的单调加载曲线和循环加载曲线峰值应力则分别增加了31.45%、37.81%；纤维掺量为0.5wt%时，FAAC的峰值应变增幅最大，BF/AAC的单调加载曲线和循环加载曲线峰值应变分别增加了28.12%、28.77%，CF/AAC的单调加载曲线和循环加载曲线峰值应变则分别增加了37.17%、41.50%；两种纤维均小幅度增加了AAC的累积塑性应变，但纤维掺量与卸载刚度及应力退化率之间未表现出明显的规律。基于试验结果，采用幂函数对FAAC标准化塑性应变与卸载点之间的关系进行拟合；提出应力退化率及加、卸载曲线双折线简化模型；最后，建立了循环荷载作用下FAAC的应力-应变曲线计算方程。
Abstract:In order to investigate the compressive mechanical properties of fiber-reinforced autoclaved aerated concrete (FAAC) under cyclic loading, a total of 11 sets of prismatic specimens were designed for uniaxial monotonic and cyclic compression tests. The effects of fiber types (basalt fiber (BF), carbon fiber (CF)) and fiber content on the mechanical performance indicators of FAAC such as failure mode, stress-strain curve characteristics, plastic strain, stiffness degradation rate, and stress degradation rate, were analyzed. The research results indicate that the failure mode of FAAC under cyclic loading contains oblique shear failure and vertical splitting failure, and with the fiber content increasing, the failure mode of FAAC turns from shear failure to splitting failure. When the fiber content is 0.4%, the peak stresses of FAAC reach their maximum values, the peak stresses of BF/autoclaved aerated concrete (AAC) under monotonic and cyclic loading increase 24.29% and 29.16%, respectively, while that of CF/AAC increase 31.45% and 37.81%, respectively. When the fiber content is 0.5%, the peak strains of FAAC reach their maximum values, the peak strains of BF/AAC under monotonic and cyclic loading increase 28.12% and 28.77%, respectively, while that of CF/AAC increase 37.17% and 41.50%, respectively. Both BF and CF slightly improve the cumulative plastic strain of AAC, but there is no significant relationship between fiber content, unloading stiffness, and stress degradation rate. Based on the experimental results, a power function is used to fit the relationship between the standardized plastic stain and unloading strain. Simplified double line models for stress degradation rate, unloading and reloading curves are proposed. Finally, the stress-strain curve calculation equation for FAAC under cyclic loading is established.
-
Keywords:
- autoclaved aerated concrete /
- fiber /
- cyclic loading /
- stress-strain curve /
- plastic strain /
- calculation equation
-
柔性压力传感器因其轻便性、柔韧性和生物兼容性等特点在航空航天、智能机器人、汽车、生物医学等领域应用广泛[1-4]。现有的柔性压力传感器按传感机制的不同可以分为压阻式、压电式、压容式和摩擦电式，其中压容式传感器因具有结构简单、稳定性强、响应迅速等优点而被广泛关注。灵敏度[5]作为压力传感单元性能中最为关键的指标之一，反映了传感器整体的力学、电学性能的优异程度。然而，制造具备出色灵敏度的柔性压力传感器也受到传感系统复杂性的严重制约，因此柔性压力传感器的灵敏度性能优化具有重要的研究价值。
触觉传感单元又称传感器触觉敏感单元，作为柔性压力传感器最为核心的组成部件，主要由聚合物纳米复合材料制备而成。聚合物纳米复合材料是以聚合物为基体连续相，以纳米级填充物为分散相的一种复合材料，其多样的基体、填料组合形式及不同的制备方式，赋予了聚合物纳米复合材料与众不同的力学[6]、电学[7]等特性，目前也广泛应用于应力/应变传感器[8]、超级电容[9]、生物医学[10]等工程领域。石墨烯和碳纳米管[11]等碳基材料因其优异的导电性和柔顺性，常作为聚合物填料应用于柔性压力传感单元的制备中。药芳萍等[12]通过设计双层砂纸与石墨烯的微组合结构制备了一种电容式柔性压力传感器，在0~25 kPa的压力载荷下灵敏度最高可达到0.451 kPa−1。Xu等[13]将碳纳米管与MXene结合到多孔微结构的聚二甲基硅氧烷海绵上研制出一种高灵敏度、宽传感范围、快响应速度的应变传感器。吴其皓等[14]以碳纳米管(CNT)与聚二甲基硅氧烷(PDMS)作为填充材料，并引入不同结构的砂纸作为模板对其电极结构进行优化，制备的轻薄型柔性传感器的灵敏度在0~1 kPa下可达到0.602 kPa−1。Ma等[15]将聚氨酯泡沫(PU)制备成“树枝状”骨架并用CNT和BaTiO3进行覆盖，在0~100 kPa的宽测量范围内灵敏度最高可达到2.51 kPa−1。Ke等[16]在热塑性聚氨酯(TPU)基体中以CNT和石墨烯片(MLG)为导电填料，并且通过优化填料比例使传感器在0~1.2 MPa压力下灵敏度最高可以达到2.05 MPa−1。Bi等[17]通过调节镍粉颗粒、碳纳米管、石墨粉混合比例的方式，优化了柔性压力传感器的灵敏度与蠕变性能，使其平均灵敏度在0~16 kPa的压力下达到
0.0589 kPa−1。王苏等[18]则采用商用聚氨酯海绵(PU)为基底，通过多次浸渍涂覆的方法将MXene/多壁碳纳米管(MWCNT)复合材料负载在海绵基底上，制备出一种柔性压力传感器，其在0~50 kPa的宽压力范围内灵敏度达到0.16 kPa−1。根据上述研究可知，设计具有不同微观结构的复合材料和掺杂具有高介电性能的填料颗粒是提高传感器灵敏度的重要手段，如引文[12-15]所述，但其材料结构设计过程复杂，耗时较长，且部分高介电颗粒的成本昂贵，故想要批量化生产此类传感单元十分困难。而基于共混方式制备传感单元的方法操作简单，耗时较少，且适合大批量生产，如引文[16-18]所述，但对于灵敏度性能的提升较微结构设计相对较低。针对以上问题，本文旨在采用操作简易的机械共混方法，以成本较低的碳纳米管(CNT)和石墨烯(MLG)作为导电填料，结合下述的人工智能优化算法，实现无微结构的柔性电容式传感单元的制备与性能优化。在确定所需的填料类型之后，由于填充颗粒之间的协同效应和反应机制的复杂性，这些填料的添加比例仍然难以确定，且在多种因素的交互影响下，不依靠数学模型对制造条件进行估计也是一件耗时漫长、成本高昂的工作。在材料合成和工艺优化领域，传统实验合成方法的分析和优化一直是重要的研究课题。最常用的统计分析模型响应面法(RSM)[19]，由Box和Wilson在1951年开发而来，其主要特点是通过少量试验创建预测模型来验证单个因素或多个因素的相互作用，用于评估和预测多参数实验条件，目前广泛应用于食品和化学领域[20]。但对于复杂的非线性问题，RSM二次多项式有时也无法达到所需的预测精度。支持向量机(SVR)作为机器学习中非线性回归模型，对非线性和非连通数据具有良好的预测能力。通过变换目标函数、改善方程的条件、降低计算的复杂度等操作[21]，可以提高其运算速度，提升预测精度与泛化能力，能在一定程度上克服RSM的缺点。虽然RSM和SVR模型已被用于各类工程领域的过程预测，但鲜有基于此类研究方法对传感器性能进行预测。蜣螂优化算法(DBO)[22]是一种新型群体智能优化算法，相较于传统的遗传算法[23] 和粒子群算法[24]，该算法的求解精度更高、收敛速度更快、全局搜索能力更强。目前该算法在三维路径规划[25]、电力信号预测[26]等问题上已有了初步应用。
综上所述，本文研究了电容式柔性压力传感单元导电机制，建立了相应的解析模型，确定了影响灵敏度性能的主要因素。采用实验设计法(DOE)中的中心组合设计法(CCD)[27]研究了CNT含量、MLG含量、搅拌时间和成型温度4个因素对灵敏度的影响。分别采用RSM和SVR机器学习算法对实验输入条件和输出性能建立预测模型，最后采用改进的DBO算法(IDBO)对该目标模型进行全局迭代寻优，确定传感器的最佳制备条件，最终实现无微结构的柔性电容式传感单元的性能优化。
1. 实验材料及方法
1.1 原材料及流程
本文所需原料包括硅橡胶基体(深圳红叶杰科技有限公司，双组分室温硫化硅橡胶RTV-2)、填充颗粒为多层石墨烯(MLG，苏州碳丰层数5~10层，厚度3.4~8 nm)、碳纳米管(CNT，苏州碳峰，纯度大于95wt%，长度大于3 μm，外径为8~15 nm)，二甲基硅油、偶联剂(KH550)与硅橡胶配套固化剂，皆购于深圳红叶杰科技有限公司。图1为压力传感器的制备流程。
![]()
图 1 触觉敏感单元制备流程图CNT—Carbon nanotubes; MLG—Graphene; KH550—Coupling agentFigure 1. Flow chart for preparation of pressure sensor首先，在室温下(25~30℃)，将一定比例的CNT和MLG加入研磨器研磨搅拌5~10 min，再加入到体积比为3∶1的乙醇水溶液中，超声分散15 min，把预处理后的填料颗粒在高温100℃左右进行干燥处理1 h。随后加入一定比例的偶联剂(KH550)进行表面改性处理，以增强分布宽度，减弱团聚效应[28]。并置于100℃高温干燥箱(DZF-6050B，上海合恒仪器设备有限公司)中进行烘干(1~3 h)，干燥后利用研钵器轻研，即得到改性后的导电填料颗粒。随后将导电填料与45 g RTV基体进行混合，同时加入质量比为5%~8%二甲基硅油进行稀释，并使用搅拌器在
1800 r/min的转速下搅拌一定的时间。搅拌完成后按混合比例1.5wt%加入硅胶固化剂(深圳红叶科技，主要成分为聚硅氧烷)再次搅拌1~2 min，并将搅拌混合后的聚合材料注入模具并放置于RS-1型真空箱(河南沃林仪器设备有限公司)中进行排气处理15 min，然后将模具放入设定的高温干燥箱中加热3 h，最后取出放置于室温下待自然冷却成形。图2为厚度4 mm、直径20 mm的圆形片状体，图3为触觉敏感单元在透射电镜(FEI Tecnai F20，日本电子) 100 nm放大倍数下改性前后的微观结构，从图3(a)可以明显看出，改性前CNT/MLG团聚现象较明显，且分布相对随机。而在改性后，如图3(b)所示，团聚明显减少，导电颗粒联系紧密且在基体中的分布相对更加均匀，保证了其稳定的传感性能。![]()
图 3 触觉敏感单元改性前(a)和改性后(b)对比Figure 3. Comparison between before (a) and after (b) modification of the tactile sensitivity unit1.2 试验方法
1.2.1 灵敏度测量及模型推理
灵敏度是衡量电容式传感器性能的重要指标，它代表了压力传感器电容量与载荷的相对变化关系[29]，式(1)为灵敏度的函数表达式。如图4所示，测量过程中使用微机控制电子万能试验机(FBS200N，弗布斯检测设备有限公司)与LCR数字电桥(UC2836B，优策电子科技有限公司)构成压力传感器电容测试回路，将样品放置在试验机中，使用试验机对样品施加0~30 kPa的单轴压力使样品产生塑性变形，同时通过LCR数字电桥和采集软件监测电容变化数据。最后得到不同压力情况下的压力传感器电容变化：
S=δ(ΔC/C0)δP=(εA/d−εA/d0)/(εA/d0)EΔd/d0 (1) 式中：S表示灵敏度(kPa−1)；C0为不施加压力时的初始电容值(pF)；ΔC=C−C0表示电容变化量，其中C为单轴压力条件下的电容值；P表示施加的压力(kPa)；E为弹性模量；d0是复合材料的初始厚度，Δd=d−d0是压缩过程中的材料厚度变化；复合压敏材料介电常数为ε；电极正对面积为A。
根据渗流理论[30-31]，在以CNT、MLG为填料的渗滤复合材料中，介电常数不再是固定的常数。渗滤复合材料作为介电层的柔性电容式压力传感器的工作原理如图5所示，在压缩过程中，填料浓度局部增加。介电常数随着浓度的增加而增加，当浓度接近渗流阈值时，介电常数增加速率变得更快，最终形成导电通路。
![]()
图 5 渗滤复合材料作为介电层的柔性电容式压力传感器的工作原理图d—Initial thickness; ε—Initial dielectric constant; P—Applied pressure; d–Δd—Material thickness after compression; ε+Δε—Dielectric constant of the composite pressure-sensitive material after compressionFigure 5. Working principle diagram of flexible capacitive pressure sensor with percolation composite as dielectric layer基于如下幂律方程所示：
ε=εd(fC−f)−t/(f<fC) (2) 其中：f和fC分别是导电填料的浓度和渗流阈值；ε和εd分别为复合材料和基体的介电常数；t表示比例常数。
故传感器的电容变化率可表示为下式，其中压敏单元厚度压缩比λ=d/d0：
ΔCC0=(εA/d−εA/d0)εA/d0=ελε0−1=[(fC−f0)/(fC−f)]tλ−1 (3) 综合上述表达式得到最终的传感器灵敏度解析模型如下：
S=[(fC−f0)/(fC−f)]t/λ−1E(λ−1) (4) 式中，f0表示初始状态的导电填料浓度。
由式(4)可知，电容式压力触觉敏感单元的灵敏度与聚合物纳米复合材料的介电性能和弹性模量密切相关。而合适的填料含量及填料在基体中的分布是否均匀是保证传感器优异介电性能的重要因素，同时根据Yang等[32]的研究可知，弹性模量E的屈服强度与温度密切相关。故本文将以填料含量、搅拌时间、成型温度作为主要影响因素，对触觉单元的灵敏度性能进行优化。
1.2.2 单因素实验设计
单因素实验设计是分别研究不同填料配比、搅拌时间、成型温度对传感器灵敏度性能的影响，为后续多因素实验设计提供一个合理的数值范围，如表1所示。试验时，首先固定MLG的含量为0.6 g，搅拌时间15 min，成型温度为 80℃，改变CNT的含量，分别是0.4、0.5、0.6、0.7、0.8、0.9、1.0、1.1 g；同理，固定 CNT的含量为0.7 g，保持搅拌时间与成型温度不变，改变MLG的含量，分别是0.4、0.5、0.6、0.7、0.8、0.9、1.0、1.1 g；而后保持CNT、MLG的含量与成型温度不变，分别为0.7 g、0.6 g、80℃，改变压力传感器制备过程中的搅拌时间分别为10、15、20、25、30 min；最后在保持填料含量与搅拌时间不变的情况下改变固化温度分别为40、60、80、100、120℃。实验重复 3 次，每个指标在相同条件下测定 3 次，最终得到同一条件下制备的传感器的平均灵敏度。
表 1 单因素实验方案Table 1. Single-factor protocolIndex Content of CNT/g Content of MLG/g Mixing time/min Molding temperature/℃ 1 0.4, 0.5, 0.6, 0.7, 0.8, 0.9,
1.0, 1.10.6 15 80 2 0.7 0.4, 0.5, 0.6, 0.7, 0.8, 0.9,
1.0, 1.115 80 3 0.7 0.6 10, 15, 20, 25, 30 80 4 0.7 0.6 15 40, 60, 80, 100, 120 1.3 预测与优化方法
1.3.1 响应面法(RSM)
响应面法(RSM)是一种用以评估一个或多个输出和输入参数之间关系建模方法，其计算量低，采用随机无偏设计，可以减少许多不必要的实验消耗，节约实验的时间与经济成本[33]。RSM中输入参数与输出变量之间的函数关系如下：
y=α0+N∑i=1αixi+N∑i=1αiix2i+∑i<j∑αijxixj (5) 其中：y是响应变量；xi是输入参数；α0表示偏差；αi表示线性效应；αii表示因素的平方效应；αij表示影响因素之间的交互效应。
RSM的实验设计法(DEO)主要分为两类：中心复合设计法(CCD)和Box-Behnken设计法(BBD)。为尽可能多的获取最佳实验条件，本文采用CCD法，进行多因素实验设计，其中每个输入参数有5个编码水平，分别为−2、−1、0、+1、+2。
图6(a)~6(d)分别为在得到每个单因素的实验结果后，根据其实验结果，以影响效果最大的点为中心点，确定了每个影响因素的实验范围，然后在实验范围内应用RSM确定多因素的实验方案，每个输入的具体水平因子显示在表2中，具体的CCD多因素实验方案如表3所示。
![]()
图 6 输入参数的单因素影响：(a) CNT含量；(b) MLG含量；(c) 搅拌时间；(d) 成型温度Figure 6. Single factor effect of input parameter: (a) CNT content; (b) MLG content; (c) Mixing time; (d) Molding temperature表 2 4个因素的水平和编码Table 2. Level and code of four factorsCode Levels factor CNT/g MLG/g Mixing
time/minMolding
temperature/℃−2 0.7 0.6 10 40 −1 0.8 0.7 15 60 0 0.9 0.8 20 80 1 1.0 0.9 25 100 2 1.1 1.0 30 120 表 3 多因素中心组合实验设计(CCD)与结果Table 3. Multifactorial central combinatorial experimental design (CCD) and resultsRun Factor Sensitivity/
(kPa−1)CNT MLG Mixing
timeMolding
temperature1 1 −1 −1 1 0.227 2 0 0 2 0 0.130 3* 0 0 0 0 0.118 4 −1 −1 −1 1 0.058 5 1 1 1 1 0.191 6 1 −1 −1 −1 0.208 7 −1 1 1 1 0.025 8* 0 0 0 0 0.145 9 −1 1 −1 1 0.050 10 −1 −1 1 −1 0.125 11 −2 0 0 0 0.050 12 0 −2 0 0 0.185 13 0 0 −2 0 0.193 14 −1 1 1 −1 0.164 15 1 1 −1 −1 0.375 16 1 −1 1 −1 0.180 17* 0 0 0 0 0.135 18 1 1 1 −1 0.310 19 0 1 0 0 0.198 20 2 0 0 0 0.478 21 −1 1 −1 −1 0.178 22 0 0 0 0 0.154 23 1 1 −1 1 0.326 24 0 0 0 −2 0.124 25* 0 0 0 0 0.117 26* 0 0 0 0 0.112 27 −1 −1 −1 −1 0.078 28 1 −1 1 1 0.248 29 −1 −1 1 1 0.090 30 0 0 0 2 0.073 Note: * indicates the central repeated trials. 1.3.2 支持向量回归(SVR)
支持向量机(SVM)是一种机器学习方法，如图7所示，它通过在高维空间中构建一个或一组超平面来实现分类或回归，其中用于回归的方法称为支持向量回归(SVR)。在SVR模型中核函数用以解决线性不可分的问题，其基本的数学表达式如下：
![]()
图 7 支持向量机回归(SVR)模型示意图Figure 7. Support vector machine regression (SVR) model diagramxn represents the independent variable; wij represents the weight coefficient; (βi, βi*) represents the Lagrange multiplier used for dual space transformation; K(x, xi) represents the kernel function used to solve the quadratic equation; b represents the deviation; y represents the dependent variabley=n∑i=1(βi−β∗i)K(x,xi)+b (6) 其中：(βi, β∗i)表示用于对偶空间变换的拉格朗日乘数；K(x,xi)表示用于求解二次方程的核函数；b表示偏差。
由于核函数是影响SVR模型精度的重要因素之一，因此为SVR模型选择合适的核函数尤为重要。目前，支持向量回归模型的核函数主要分为4种类型，即线性、多项式、径向基函数(RBF)和sigmoid。此外，核函数选取后，算法的超参数如g、C和ε也会影响SVR模型的性能，因此超参数的选取也是SVR建模的关键步骤。在本文中，RBF被用作模型计算的核函数，超参数优化则采用K折交叉验证的网格搜索方法(图8)。
1.3.3 蜣螂优化算法及其改进(IDBO)
蜣螂优化算法(DBO)作为一种新型启发式群智能优化算法，通过模拟蜣螂的滚球、跳舞、觅食、偷窃和繁殖等行为实现种群全局范围内最优值的选择。其适应度值作为衡量优化结果优劣的标准，优化的具体流程如图9所示。
![]()
图 9 蜣螂优化算法(DBO)流程图Figure 9. Dung beetle optimization (DBO) algorithm flowchartt—Number of iterations at the present stage; T—Number of iterations set by the algorithm同其他群智能优化算法类似，DBO通过随机生成位置的方式实现种群的初始化，但这容易出现种群分布不均匀，导致种群多样性减少，种群质量不高，从而影响算法的收敛速度及搜索性能。为了能够让个体在算法开始时有较高的全局搜索能力，需要让种群的位置均匀分布在整个问题的解空间内，混沌作为一种非线性的自然现象，因其混沌序列具有遍历性、随机性等优点，被广泛用于优化搜索问题。利用混沌变量搜索显然比无序随机搜索具有更大的优越性，因此本文采用改进的Tent混沌映射[34]替换原算法的初始化。原算法初始化与改进的Tent混沌初始化表达式如下：
xi,j=lbj+ri,j(ubj−lbj) (7) 式中：ri,j表示随机数列；lbj、ubj分别表示算法种群维度的上下界。
{2xi+rand(0,1)1N, 0⩽ (8) 式中， N 是序列内粒子的个数。引入随机变量 \mathrm{r}\mathrm{a}\mathrm{n}\mathrm{d}\left(0,1\right)/N 不仅仍然保持了Tent混沌映射的随机性、遍历性、规律性，而且能够有效避免迭代落入小周期点和不稳定周期点内。
在对种群初始化进行优化后，为了提升算法在迭代的不同时期全局/局部搜索能力。本文参考了文献[35]采用自适应 t 分布扰动策略，使DBO算法获得更优异的搜索性能。改进后的全局搜索位置更新公式为
\left\{\begin{array}{l}\mathrm{if\ rand} < p \\ x_i^t=x_i+x_it\left(i_{\mathrm{t}\mathrm{e}\mathrm{r}}\right)\end{array}\right. (9) 式中： p=0.5-0.1\left(T-t\right)/T 表示扰动概率； {x}_{{{i}}}^{t} 为突变后的蜣螂位置； {x}_{{{i}}} 为第 i 只蜣螂个体位置； t\left({i}_{\mathrm{t}\mathrm{e}\mathrm{r}}\right) 是以迭代次数 {i}_{\mathrm{t}\mathrm{e}\mathrm{r}} 为自由度的 t 分布函数；T为算法设定的迭代次数。
通过自适应 t 分布突变扰动策略虽然能增强算法全局搜索和跳出局部最优的能力，但是无法确定扰动之后得到的新位置一定比原位置的适应度值要好，因此在进行变异扰动更新后，加入贪婪策略[36]，通过比较新旧两个位置的适应度值，确定是否要更新位置。贪婪策略表达式如下：
{x}_{{{i}}}^{\mathrm{n}\mathrm{e}\mathrm{w}}\left(t\right)=\left\{\begin{array}{l}{x}_{{{i}}}^{}\left(t\right),\left(f\right({x}_{{{i}}}^{}\left(t\right)) < f({x}_{{\mathrm{g}}}\left(t\right)))\\ {x}_{{\mathrm{g}}}\left(t\right),\left(f\right({x}_{{{i}}}^{}\left(t\right)) > f({x}_{{\mathrm{g}}}\left(t\right)))\end{array}\right. (10) 其中： {x}_{{{i}}}^{\mathrm{n}\mathrm{e}\mathrm{w}}\left(t\right) 为经过贪婪规则更新后的个体位置； {x}_{{{i}}}\left(t\right) 为扰动后的个体； {x}_{{\mathrm{g}}}\left(t\right) 为扰动前的个体； f 代表蜣螂适应度函数。
IDBO算法的具体流程如图10所示。
![]()
图 10 改进蜣螂优化算法(IDBO)流程图P— Disturbance probabilityFigure 10. Improved dung beetle optimization (IDBO) algorithm flowchart为了验证所提蜣螂优化算法改进策略的有效性，本文将改进的IDBO与4种被广泛研究的优化算法(粒子群算法(PSO)、遗传算法(GA)、灰狼优化算法(GWO)、DBO)进行了比较。选取多种标准测试函数[37]评价其性能表现。其中，F1~F3为高维单峰测试函数，F4~F6为高维多峰测试函数。为保证实验公平性，所有算法的迭代次数与种群数都设置为100，结果如图11所示。
![]()
图 11 算法性能对比收敛曲线图GA—Genetic algorithm; PSO—Particle swarm optimization; GWO—Grey wolf optimizerFigure 11. Algorithm performance comparison convergence plot由图11可以看出，改进后的蜣螂算法在搜索精度及收敛速度上皆优于其他算法，其性能得到了显著提高。
1.3.4 预测模型精度的评价指标
在本文中，由于RSM、SVR对于传感单元灵敏度的性能预测本质上是一种回归问题，因此采用决定系数(R2)、均方根误差(Rmse)和平均绝对误差(Mae)用于评估二者所建立的预测模型的准确性，具体表达式如下：
{R}^{2}=1-\frac{\displaystyle\sum _{i=1}^{n}{\left({Y}_{{{i}}}-{{Y}_{\mathrm{p}\mathrm{r}\mathrm{e}\mathrm{d}}}_{{{i}}}\right)}^{2}}{\displaystyle\sum _{i=1}^{n}{\left({Y}_{{{i}}}-\overline{Y}\right)}^{2}} (11) {R}_{\mathrm{m}\mathrm{s}\mathrm{e}}=\sqrt{\frac{1}{N}\sum _{i=1}^{N}{\left({Y}_{{{i}}}-{{Y}_{\mathrm{p}\mathrm{r}\mathrm{e}\mathrm{d}}}_{{{i}}}\right)}^{2}} (12) {M}_{\mathrm{a}\mathrm{e}}=\frac{1}{N}\sum _{i=1}^{N}\left|{Y}_{{{i}}}-{{Y}_{\mathrm{p}\mathrm{r}\mathrm{e}\mathrm{d}}}_{{{i}}}\right| (13) 其中： {Y}_{{{i}}} 和 {{Y}_{\mathrm{p}\mathrm{r}\mathrm{e}\mathrm{d}}}_{{{i}}} 分别是通过实验测量的实际灵敏度和通过应用估算模型估算的灵敏度； \overline{Y} 是通过实验所测得灵敏度的平均值；N为本次实验测量的灵敏度样品的总个数。
2. 结果与讨论
2.1 响应面(RSM)模型的建立与分析
本实验采用Design-Expert软件对数据进行处理和分析，并通过方差分析(ANOVA)以检查RSM模型对数据的拟合情况，从而确定自变量对灵敏度的影响程度。表4为灵敏度方差分析结果。可以看出，该回归模型的F-value为32.91，P-value小于
0.0001 ，而失拟项的P-value为0.1422， 大于0.05，表明整体回归模型高度显著。RSM模型的R2为0.9685 ，证明了该模型的拟合程度较高。该回归模型的方差系数(CV)为14.86%，低于15%，CV值较低，表明其置信度较高，模型能很好地反映真实的实验数据。基于F-value的显著性排序为CNT>MLG>成型温度>搅拌时间。根据CCD的ANOVA回归系数，灵敏度与填料配比、搅拌时间、成型温度的响应关系可表示为表 4 灵敏度响应的方差分析Table 4. Analysis of variance for sensitivity responseSource Sum of squares df Mean square F-value P-value Status Model 0.2877 14 0.0205 32.91 < 0.0001 Significant {x}_{1} -CNT 0.1931 1 0.1931 309.36 < 0.0001 {x}_{2} -MLG 0.0117 1 0.0117 18.789 0.0006 {x}_{3} -Mixing time 0.0036 1 0.0036 5.73 0.0302 {x}_{4} -Temperature 0.0106 1 0.0106 17.02 0.0009 {x}_{1}{x}_{2} 0.0047 1 0.0047 7.46 0.0155 {x}_{1}{x}_{3} 0.0038 1 0.0038 6.11 0.0259 {x}_{1}{x}_{4} 0.0036 1 0.0036 5.81 0.0292 {x}_{2}{x}_{3} 0.0060 1 0.0060 9.68 0.0071 {x}_{2}{x}_{4} 0.0136 1 0.0136 21.83 0.0003 {x}_{3}{x}_{4} 0.0001 1 0.0001 0.2211 0.6449 {{x}_{1}}^{2} 0.0267 1 0.0267 42.79 < 0.0001 {{x}_{2}}^{2} 0.0085 1 0.0085 13.65 0.0022 {{x}_{3}}^{2} 0.0008 1 0.0008 1.26 0.2795 {{x}_{4}}^{2} 0.0030 1 0.0030 4.88 0.0432 Residual 0.0094 15 0.0006 Lack of fit 0.0079 10 0.0008 2.70 0.1422 Not significant Pure error 0.0015 5 0.0003 Cor total 0.2970 29 {R}^{2} 0.9685 {R}_{\mathrm{A}\mathrm{d}\mathrm{j}}^{2} 0.9390 {R}_{\mathrm{P}\mathrm{r}\mathrm{e}\mathrm{d}}^{2} 0.8311 CV/% 14.86 Adequate precision 23.8126 Notes: Cor total—Total number of correlation coefficients; {R}_{\mathrm{A}\mathrm{d}\mathrm{j}}^{2} —Adjusted coefficient of determination; {R}_{\mathrm{P}\mathrm{r}\mathrm{e}\mathrm{d}}^{2} —Predicted coefficient of determination; Adequate precision is an important statistical metric for evaluating the predictive ability of response surface models; CV—Coefficient of variance; df—Number of independent variables. \begin{split} Y=&0.1324+0.0897{x}_{1}+0.0250{x}_{2}-0.0122{x}_{3}-\\ &0.0210{x}_{4}+0.0171{x}_{1}{x}_{2}-0.0154{x}_{1}{x}_{3}+0.0151{x}_{1}{x}_{4}- \\ &0.0194{x}_{2}{x}_{3}-0.0292{x}_{2}{x}_{4}-0.0029{x}_{3}{x}_{4}+0.0309{{x}_{1}}^{2}+\\ &0.0228{{x}_{2}}^{2}+0.0053{{x}_{3}}^{2}-0.0104{{x}_{4}}^{2}\\[-1pt]\end{split} (14) 式中： Y 表示响应值，代表预测灵敏度； {x}_{1} 、 {x}_{2} 分别表示CNT和MLG的含量； {x}_{3} 代表搅拌时间； {x}_{4} 表示成型温度。
图12展示了各因素之间交互作用对灵敏度的影响，当研究两个因素之间交互时，其余两个因素处于中心水平状态。
![]()
图 12 灵敏度三维响应面图：(a) CNT/MLG相互作用；(b) CNT/搅拌时间相互作用；(c) CNT/成型温度相互作用；(d) MLG/搅拌时间相互作用；(e) MLG/成型温度相互作用；(f)搅拌/成型温度相互作用Figure 12. Three-dimensional response surface plots of sensitivity: (a) Interaction between CNT/MLG content; (b) Interaction between CNT/mixing time; (c) Interaction between CNT/molding temperature; (d) Interaction between MLG/mixing time; (e) Interaction between MLG/molding temperature; (f) Interaction between mixing time/molding temperature由图12(a)可以看出，随着CNT和MLG含量的提升，灵敏度的提升十分显著，这是由于这类碳基材料具有优异的导电性和柔顺性，两者结合能够实现基于二维石墨烯蜂窝结构与一维碳纳米管的三维导电网络，从而进一步增强其分散性[38]，增加填料与基体间的接触面积[39]，从而提升传感器的介电性能，这也与式(4)的推论相互印证。根据图12(b)、图12(d)可知，在CNT或MLG含量一定时，传感器的灵敏度性能随着搅拌时间的增加改变并不明显，而随着搅拌时间的继续增加，灵敏度有下降的趋势。这是由于长时间的混合剪切破坏了填料内部的三维导电网络结构[40]，这也表明适合的搅拌时间对压力传感器灵敏度的重要性。从图12(c)、图12(e)可知，在固化过程中随着温度的升高灵敏度逐渐下降。这是由于在高温固化的过程中基体持续软化使CNT和MLG之间的平均间距减小，这样便容易发生导电颗粒的团聚现象，从而对压力传感器的力学性能、电学性能产生了负面影响[41]。而从图12(f)可以了解到，在没有导电填料的情况下，仅对硅橡胶(RTV)基体进行搅拌和高温固化并不会对传感器的性能造成实质性的影响。
2.2 支持向量机回归(SVR)模型的建立与优化
基于上述CCD的数据样本，以CNT含量、MLG含量、搅拌时间和成型温度为输入参数，灵敏度为预测的输出变量建立预测模型。如图13所示，利用网格搜索法优化SVR的超参数，当C为
68.5935 ，g为0.7071 时，在5折交叉验证下的均方误差(MSE)达到最小值，平均R2值最大，即模型基于数据样本的预测性能达到最优。其中，基于SVR灵敏度预测模型的R2、Rmse和Mae值如表5所示。表 5 SVR模型的5折交叉验证预测性能评价指标Table 5. Performance evaluation index of SVR model by 5-fold cross-validationIndex {R}^{2} Rmse Mae 1 0.9893 0.0253 0.0193 2 0.9841 0.0328 0.0301 3 0.9735 0.0354 0.0231 4 0.9769 0.0389 0.0321 5 0.9882 0.0188 0.0151 Average 0.9824 0.0320 0.0239 Notes: Rmse—Root mean square error; Mae—Mean error rate. 从表5可以看出，在灵敏度预测性能评估中，平均R2、平均Rmse和平均Mae分别为
0.9824 、0.0320 和0.0239 ，这表明了预测值和实验值之间的高度相关性。此外，R2、Rmse和Mae的范围是0.9735 ~0.9893 、0.0188 ~0.0389 和0.0151 ~0.0321 ，这表明模型预测的波动范围极小，证明了此预测模型的高鲁棒性、高置信度。2.3 RSM与SVR预测模型对比
如表6所示，为比较两种模型的预测性能，基于式(10) RSM回归模型对相同数据进行了5折交叉验证。
表 6 响应面法(RSM)模型的5折交叉验证预测性能评价指标Table 6. Performance evaluation index of response surface methodology (RSM) model by 5-fold cross-validationIndex R2 Rmse Mae 1 0.9563 0.04002 0.03623 2 0.9642 0.02701 0.02443 3 0.9778 0.03103 0.02671 4 0.9669 0.04731 0.03113 5 0.9782 0.03312 0.02608 Average 0.96868 0.035698 0.028916 根据表5、表6的平均比较结果表明，SVR模型基于灵敏度性能的预测平均R2大于RSM模型，表明其适用性更好。SVR模型基于灵敏度预测的平均Rmse和Mae均小于RSM模型，更加表明超参优化后的SVR模型的预测精度优于RSM模型。
图14为SVR和RSM模型的5折交叉验证的雷达图。雷达图中可更加清楚地了解5折交叉验证下两预测模型的性能表现，图14(a)中SVR模型灵敏度的R2面积大于RSM模型，图14(b)、图14(c)中Rmse和Mae的面积均小于RSM模型，充分说明SVR模型在灵敏度的预测上具有更好的准确性。因此，通过此次比较分析表明，RSM和SVR模型在灵敏度性能预测方面都表现出良好的适应性和鲁棒性，而相较于RSM模型，SVR模型具有更好的预测效果。
![]()
图 14 SVR和RSM模型的5折交叉验证的雷达图Figure 14. Radar diagrams of SVR and RSM models with 5 fold cross validation2.4 基于改进蜣螂算法(IDBO)的模型优化
此优化过程是为了获得与灵敏度相适应的最佳材料配比及制作成型的工艺参数。通过上述的研究分析可知，由于通过超参优化后的SVR模型预测性能相优于RSM模型，故选择SVR模型用作优化的目标函数模型。通过IDBO算法不断迭代寻找最优的工艺参数，其中种群数量设为50，变量维度为4，迭代次数100次，迭代过程中的群体最优适应度值如图15所示。可知，随着迭代次数的增加，种群最优的适应度值也在不断上升，直至找到SVR模型的最优解，经过10次迭代后达到最佳适应度，即灵敏度性能达到最优，随后基本保持不变。最优的材料配比及制作工艺为CNT占比2.3wt%、MLG占比1.9wt%、搅拌时间15 min、成型温度78℃，仿真灵敏度为
0.5512 kPa−1。2.5 实验验证与对比
为了验证上述实验分析的合理性，根据上述优化结果的数据进行传感单元样品的制作，图16为所制作传感单元的灵敏度性能图，拟合直线斜率表示压力传感器在0~30 kPa时灵敏度大小。可知，采用优化后数据所制备的传感单元灵敏度S相较未优化时有大幅提升，达到
0.5371 kPa−1，与仿真灵敏度的相对误差为2.625%，也证明了基于SVR-IDBO优化预测模型的可靠性。![]()
图 16 灵敏度测试电容输出响应曲线S—Sensitivity; ΔC/C0—Relative change of the capacitanceFigure 16. Sensitivity test capacitor output response curves表7为本文所制备的柔性电容传感单元与近年来相关柔性压力传感器的灵敏度性能比较情况[12, 42-45]。结果表明，在相同压力载荷下，本文的传感器单元灵敏度性能处于较高的水平。其中，在同等压力下与张鹏等[45]的研究相比，灵敏度性能提升了5倍，与Yan等[44]相比提升了1.4倍，与药芳萍等[12]、赵珂[42]相比较，灵敏度性能分别提升了1.7倍与2.5倍。且与原始实验相比灵敏度最高提升了20倍。进一步表明结合SVR-IDBO方法进行触觉敏感单元性能优化的有效性与优越性。
表 7 传感单元灵敏度性能对比Table 7. Comparison of sensitivity performance of sensing unitsMethod Sensitivity/ \mathrm{kP}\mathrm{a}^{-1} Pressure range/ \mathrm{k}\mathrm{P}\mathrm{a} Comparison/ \mathrm{kP}\mathrm{a}^{-1} Ref. Microstructure design 0.451 0-25 0.7665 [12] Microstructure design 0.69 0-10 1.7531 [42] Microstructure design 1.01 0-25 0.7665 [43] Mechanical blending 0.3850 0-30 0.5371 [44] Microstructure design 0.341 0-5 1.9698 [45] Mechanical blending 0.5371 0-30 – This study 3. 相关性能测试
为满足传感器在实际场景的应用，本文的柔性触觉敏感单元除了在宽传感范围内具备较高的灵敏度之外，还具备其他较理想的特性，包括快速响应、高重复稳定性、良好的耐用性等。
如图17(a)所示，使用万能试验机对触觉敏感单元进行加载和卸载，在10 mm/min的较慢载荷进行速度下，其响应时间和恢复时间均能达到1 s，且每一循环峰值都极为接近，表明其在快速响应的同时也保证了电容信号的稳定传输。
如图17(b)、图17(c)所示，在15 kPa、30 kPa、45 kPa的载荷范围内进行10次循环测试，均能表现出良好的电容输出稳定性。且通过施加0~30 kPa的300次循环载荷也体现出触觉单元在受到载荷的中后期具有良好的重复稳定性。
4. 结 论
采用机器学习与智能优化算法相结合的方式改善了基于机械共混的柔性压力传感单元的制备条件。主要工作如下：
(1)通过推导压容式传感器灵敏度的机制，确定了影响其性能的主要因素有碳纳米管(CNT)含量、石墨烯(MLG)含量、搅拌时间和成型温度。采用中心复合实验设计法(CCD)动态检测了各种因素对输出性能的影响，结果表明，所有的因素均对输出性能有显著影响，并得到了响应面法(RSM)回归模型；
(2)基于CCD的实验数据作为输入，利用超参优化后的支持向量机(SVR)算法对灵敏度进行预测建模，并与RSM预测模型进行了比较，结果表明两种方法对灵敏度性能均表现出极高的预测精度，特别是超参优化后的SVR模型的预测结果与实际结果十分接近，这也表明所提出的方法在预测多因素影响和传感器性能之间的复杂关系方面具有巨大的应用潜力；
(3)利用改进的蜣螂优化算法(IDBO)对SVR模型进行迭代优化，逆向解析出了柔性传感器单元最佳的制备条件，实验结果表明，灵敏度性能大幅提升，且与仿真结果的相对误差较低。
综上所述，本文所采用的研究方法为传感单元性能优化提供了参考，同时为智能优化算法在柔性压力传感器方面的应用提供了理论依据，并且对实现柔性传感单元的大规模、工业化生产具有一定的指导意义。
-
![]()
图 3 循环荷载作用下纤维增强蒸压加气混凝土(FAAC)破坏过程
Figure 3. Failure process of fiber-reinforced autoclaved aerated concrete (FAAC) under cyclic loading
{\varepsilon }_{\mathrm{c}} , {\sigma }_{\mathrm{c}} , and {\varepsilon }_{\mathrm{p}} are the peak strain, peak stress, and plastic strain of the curve, respectively; {\varepsilon }_{\mathrm{u}\mathrm{n}-1} and {\sigma }_{\mathrm{u}\mathrm{n}-1} are the unloading strain and unloading stress of the previous level unloading curve, respectively; {\varepsilon }_{\mathrm{u}\mathrm{n}} and {\sigma }_{\mathrm{u}\mathrm{n}} are the unloading strain and unloading stress of the next level unloading curve
![]()
图 5 循环荷载作用下FAAC最终的破坏现象
Figure 5. Final failure phenomenon of FAAC under cyclic loading
![]()
图 6 单调及循环荷载作用下FAAC应力(σ)-应变(ε)曲线
Figure 6. Stress (σ)-strain (ε) curves of FAAC under monotonic and cyclic loading
![]()
图 7 FAAC塑性应变与卸载点应变关系
R2—Fitting degree
Figure 7. Relationship between plain strain and unloading strain of FAAC
![]()
图 10 单调荷载作用下FAAC的归一化应力-应变曲线
Figure 10. Normalized stress-strain curves of FAAC under monotonic loading
![]()
图 11 纤维掺量与单调加载曲线拟合参数之间的关系
c, d—Fitting parameters; wb, wc—Mass content of basalt fibers and carbon fibers, respectively
Figure 11. Relationship between fiber content and fitting parameters of stress-strain curves under monotonic loading
![]()
图 12 卸载、再加载曲线简化模型
Figure 12. Simplified model for unloading and reloading curves
{\sigma }_{\mathrm{u}\mathrm{n},\mathrm{r}} —Stress of the point on reloading curve whose strain is equal to that of unloading strain of previous unloading curve
![]()
图 13 FAAC试验曲线与简化计算方程对比
Figure 13. Comparison between test curve and calculation equation of FAAC
表 1 纤维基本物理力学性能
Table 1 Basic physical and mechanical properties of fiber
Fiber Density/(g·cm−3) Length/mm Diameter/μm Tensile strength/MPa Melt point/℃ Elastic modulus/GPa CF 1.75 3 7 4900 800-900 230 BF 2.63-2.65 3 7-15 3000 -4800 1050 91-110 Notes: CF—Carbon fiber; BF—Basalt fiber. 
下载: 导出CSV
表 2 蒸压加气混凝土(AAC)试件配合比
Table 2 Mix proportion of autoclaved aerated concrete (AAC)
Tailing sand/wt% Lime/wt% Cement/wt% Gypsum/wt% Aluminum powder/wt% Water/Solid materials 58.3 12.9 25.8 3 0.08 0.50
下载: 导出CSV
表 3 AAC试件参数
Table 3 Specimen design parameters of AAC
No. Fiber Fiber content/wt% No. Fiber Fiber content/wt% AAC – – AAC – – BF/AAC-0.1 BF 0.1 CF/AAC-0.1 CF 0.1 BF/AAC-0.2 0.2 CF/AAC-0.2 0.2 BF/AAC-0.3 0.3 CF/AAC-0.3 0.3 BF/AAC-0.4 0.4 CF/AAC-0.4 0.4 BF/AAC-0.5 0.5 CF/AAC-0.5 0.5 Notes: BF/AAC-n—Basalt fiber reinforced autoclaved aerated concrete with fiber content of nwt%, n varies from 0.1 to 0.5; CF/AAC-n—Carbon fiber reinforced autoclaved aerated concrete with fiber content of nwt%, n varies from 0.1 to 0.5.
下载: 导出CSV
表 4 单调及循环荷载作用下FAAC试件的峰值应力及峰值应变
Table 4 Peak stress and peak strain of FAAC under monotonic and cyclic loading
No. Monotonic loading Cyclic loading No. Monotonic loading Cyclic loading {\varepsilon }_{\mathrm{c}} /10−3 {\sigma }_{\mathrm{c}} /MPa {\varepsilon }_{\mathrm{c}} /10−3 {\sigma }_{\mathrm{c}} /MPa {\varepsilon }_{\mathrm{c}} /10−3 {\sigma }_{\mathrm{c}} /MPa {\varepsilon }_{\mathrm{c}} /10−3 {\sigma }_{\mathrm{c}} /MPa AAC 3.958 1.968 3.945 1.883 AAC 3.958 1.968 3.945 1.883 BF/AAC-0.1 4.268 2.056 4.293 2.015 CF/AAC-0.1 4.327 2.194 4.418 2.152 BF/AAC-0.2 4.385 2.178 4.477 2.170 CF/AAC-0.2 4.545 2.341 4.490 2.246 BF/AAC-0.3 4.773 2.310 4.840 2.347 CF/AAC-0.3 5.089 2.518 5.020 2.454 BF/AAC-0.4 4.805 2.446 4.831 2.432 CF/AAC-0.4 5.343 2.587 5.496 2.595 BF/AAC-0.5 5.071 2.357 5.080 2.282 CF/AAC-0.5 5.429 2.444 5.582 2.399
下载: 导出CSV
-
[1] WANG C L, NI W, ZHANG S Q, et al. Preparation and properties of autoclaved aerated concrete using coal gangue and iron ore tailings[J]. Construction and Building Materials, 2016, 104: 109-115. DOI: 10.1016/j.conbuildmat.2015.12.041
[2] 陈潇, 张浩宇, 薛鑫, 等. 固体废弃物在蒸压加气混凝土中的应用现状综述[J]. 硅酸盐通报, 2023, 42(2): 541-553. DOI: 10.3969/j.issn.1001-1625.2023.2.gsytb202302018 CHEN Xiao, ZHANG Haoyu, XUE Xin, et al. Review on application of solid wastes in autoclaved aerated concrete[J]. Bulletin of the Chinese Ceramic Society, 2023, 42(2): 541-553(in Chinese). DOI: 10.3969/j.issn.1001-1625.2023.2.gsytb202302018
[3] PEHLIVANLI Z O, IBRAHIM U. Effect of polypropylene fiber length on mechanical and thermal properties of autoclaved aerated concrete[J]. Construction and Building Materials, 2022, 322: 126506. DOI: 10.1016/j.conbuildmat.2022.126506
[4] 黄炜, 张敏, 江永涛, 等. 装配式混凝土墙抗震性能试验研究[J]. 建筑结构学报, 2015, 36(10): 88-95. HUANG Wei, ZHANG Min, JIANG Yongtao, et al. Experimental study on seismic behavior of precast concrete walls[J]. Journal of Building Structures, 2015, 36(10): 88-95(in Chinese).
[5] MIAO X W, HUANG W, FAN Z H, et al. Mechanical property test and numerical analysis of a novel precast shear wall[J]. Engineering Structures, 2024, 300: 117236. DOI: 10.1016/j.engstruct.2023.117236
[6] HUANG W, AN Y J N, MIAO X W, et al. Research on seismic behavior and shear strength of precast composite walls with different reinforced concrete braces[J]. Structures, 2024, 61: 106067. DOI: 10.1016/j.istruc.2024.106067
[7] 陈国新, 黄炜, 张荫. 内填不同材料填充砌块生态复合墙体抗震性能对比[J]. 中南大学学报(自然科学版), 2012, 43(11): 4491-4500. CHEN Guoxin, HUANG Wei, ZHANG Yin. Comparison on seismic behavior of ecological composite walls filled with different materials block[J]. Journal of Central South University (Science and Technology), 2012, 43(11): 4491-4500(in Chinese).
[8] BONAKDAR A, BABBITT F, MOBASHER B. Physical and mechanical characterization of fiber-reinforced aerated concrete (FRAC)[J]. Cement and Concrete Composites, 2013, 38: 82-91.
[9] LAUKAITIS A, KERIENE J, MIKULSKIS D, et al. Influence of fibrous additives on properties of aerated autoclaved concrete forming mixtures and strength characteristics of products[J]. Construction & Building Materials, 2009, 23(9): 3034-3042.
[10] QUAN W L, HUANG W, AN Y J N, et al. The effect of natural bamboo fiber and basalt fiber on the properties of autoclaved aerated concrete[J]. Construction and Building Materials, 2023, 377: 131153. DOI: 10.1016/j.conbuildmat.2023.131153
[11] XU R S, HE T S, DA Y Q, et al. Utilizing wood fiber produced with wood waste to reinforce autoclaved aerated concrete[J]. Construction and Building Materials, 2019, 208: 242-249. DOI: 10.1016/j.conbuildmat.2019.03.030
[12] 张杰, 黄斐, 刘文地, 等. 改性竹纤维加气混凝土的制备与界面特性[J]. 建筑材料学报, 2022, 25(7): 686-692. DOI: 10.3969/j.issn.1007-9629.2022.07.005 ZHANG Jie, HUANG Fei, LIU Wendi, et al. Preparation and interfacial characteristics of modified bamboo fibers reinforced autoclaved aerated concrete[J]. Journal of Building Materials, 2022, 25(7): 686-692(in Chinese). DOI: 10.3969/j.issn.1007-9629.2022.07.005
[13] HUANG F, ZHANG J, ZHENG X Y, et al. Preparation and performance of autoclaved aerated concrete reinforced by dopamine-modified polyethylene terephthalate waste fibers[J]. Construction and Building Materials, 2022, 348: 128649. DOI: 10.1016/j.conbuildmat.2022.128649
[14] 彭军芝, 彭小芹, 桂苗苗, 等. 蒸压加气混凝土孔结构表征的图像分析方法[J]. 材料导报, 2011, 25(2): 125-129. PENG Junzhi, PENG Xiaoqin, GUI Miaomiao, et al. Pore structure characterization of autoclaved aerated concrete using image analysis method[J]. Materials Reports, 2011, 25(2): 125-129(in Chinese).
[15] CHEN G L, LI F L, GENG J Y, et al. Identification, generation of autoclaved aerated concrete pore structure and simulation of its influence on thermal conductivity[J]. Construction and Building Materials, 2021, 294: 123572. DOI: 10.1016/j.conbuildmat.2021.123572
[16] 孟宏睿. 生态轻质水泥基墙体材料性能及密肋复合墙体弹塑性分析模型研究[D]. 西安: 西安建筑科技大学, 2007. MENG Hongrui. Research on performance of ecological lightweight cement-based wall materials & elasto-plastic analysis model of multi-ribbed wall[D]. Xi'an: Xi'an University of Architecture & Technology, 2007(in Chinese).
[17] 熊耀清, 姚谦峰. 轻质多孔混凝土受压应力-应变全曲线试验研究[J]. 四川建筑科学研究, 2010, 36(2): 228-232. DOI: 10.3969/j.issn.1008-1933.2010.02.059 XIONG Yaoqing, YAO Qianfeng. Experimental study on the total stress-strain curve of porous lightweight concrete[J]. Sichuan Building Science, 2010, 36(2): 228-232(in Chinese). DOI: 10.3969/j.issn.1008-1933.2010.02.059
[18] 陈国新. 内填不同材料生态复合墙体基于统一强度理论的非线性损伤分析[D]. 西安: 西安建筑科技大学, 2010. CHEN Guoxin. Analysis on nonlinear damage of ecological composite walls filled with different materials based on the twin shear unified strength theory[D]. Xi'an: Xi'an University of Architecture & Technology, 2010(in Chinese).
[19] 中国国家标准化管理委员会. 蒸压加气混凝土性能试验方法: GB/T 11969—2020[S]. 北京: 中国标准出版社, 2020. Standardization Administration of the People's Republic of China. Test method for performance of autoclaved aerated concrete: GB/T 11969—2020[S]. Beijing: Standards Press of China, 2020(in Chinese).
[20] 陈宇良, 李浩, 叶培欢, 等. 循环荷载作用下钢纤维再生混凝土力学性能试验[J]. 复合材料学报, 2022, 39(11): 5574-5585. CHEN Yuliang, LI Hao, YE Peihuan, et al. Experimental study on mechanical behavior of steel fiber recycled concrete under cyclic compression[J]. Acta Materiae Compositae Sinica, 2022, 39(11): 5574-5585(in Chinese).
[21] 徐礼华, 黄彪, 李彪, 等. 循环荷载作用下聚丙烯纤维混凝土受压应力-应变关系研究[J]. 土木工程学报, 2019, 52(4): 1-12. XU Lihua, HUANG Biao, LI Biao, et al. Study on the stress-strain relation of polypropylene fiber reinforced concrete under cyclic compression[J]. China Civil Engineering Journal, 2019, 52(4): 1-12(in Chinese).
[22] 徐礼华, 李长宁, 李彪, 等. 循环受压状态下钢纤维混凝土一维弹塑性损伤本构模型研究[J]. 土木工程学报, 2018, 51(11): 77-87. XU Lihua, LI Changning, LI Biao, et al. Investigation on 1D elasto-plastic constitutive model of steel fiber reinforced concrete under uniaxial cyclic compression[J]. China Civil Engineering Journal, 2018, 51(11): 77-87(in Chinese).
[23] 陈宗平, 覃钦泉, 梁莹, 等. 聚丙烯纤维珊瑚海水混凝土循环受压试验及应力-应变本构关系[J]. 复合材料学报, 2024, 41(8): 4259-4271. CHEN Zongping, QIN Qinquan, LIANG Ying, et al. Cyclic compression test and stress-strain constitutive relationship of polypropylene fiber coral seawater concrete[J]. Acta Materiae Compositae Sinica, 2024, 41(8): 4259-4271(in Chinese).
[24] 过镇海, 张秀琴, 张达成, 等. 混凝土应力-应变全曲线的试验研究[J]. 建筑结构学报, 1982, 3(1): 1-12. GUO Zhenhai, ZHANG Xiuqin, ZHANG Dacheng, et al. Experimental investigation of the complete stress-strain curve of concrete[J]. Journal of Building Structures, 1982, 3(1): 1-12(in Chinese).
[25] 徐子豪, 胡晓斌, 张文良, 等. 重复加载下再生混凝土单轴受压应力应变关系[J]. 工业建筑, 2018, 48(7): 115-1211. XU Zihao, HU Xiaobin, ZHANG Wenliang, et al. The uniaxial compressive stress-strain relation of recycled concrete under cyclic loading[J]. Industrial Construction, 2018, 48(7): 115-121(in Chinese).
-
其他相关附件
-
目的
近年来，随着装配式建筑的发展，蒸压加气混凝土(AAC)被广泛应用于建筑内隔墙、屋面板等非结构部位及装配式复合结构体系中的内填块材等结构部位。然而，实际工程结构常常遭受地震等循环荷载作用，为改善AAC的力学性能并使其更好地应用于实际工程结构中，本文在AAC基体中掺入纤维增强材料，通过循环受压试验探究纤维增强蒸压加气混凝土(FAAC)的循环受压力学行为。
方法以纤维种类(玄武岩纤维(BF)；碳纤维(CF))和纤维掺量(0.1wt%, 0.2wt%, 0.3wt%, 0.4wt%, 0.5wt%)为变量，共设计11组棱柱体(100mm×100mm×300mm)试件进行单轴单调受压试验和单轴循环受压试验，分析FAAC破坏形态、应力-应变全曲线特征、塑性应变、刚度退化率、应力退化率等力学性能指标，进而建立循环荷载作用下FAAC的应力-应变曲线计算方程。
结果(1)试验现象：单调及循环荷载作用下各试件破坏过程相似；随着纤维掺量的增加，FAAC最终破坏模式由斜向剪切破坏向竖向劈裂破坏转变。(2)应力-应变曲线特征：①单调及循环荷载作用下各试件应力-应变曲线的外包络线相近，均包含上升段、下降段和稳定阶段。②随着纤维掺量增加，玄武岩增强蒸压加气混凝土(BF/AAC)和碳纤维增强蒸压加气混凝土(CF/AAC)的峰值应力均呈先增加后减小的趋势，峰值应变则呈波动上升趋势；纤维掺量为0.4%时，BF/AAC和CF/AAC试件的峰值应力均达到最大值，BF/AAC试件单调及循环应力-应变曲线的峰值应力分别增加了24.29%和29.16%，CF/AAC试件则分别增加了31.45%和37.81%；纤维掺量为0.5%时，BF/AAC和CF/AAC试件的峰值应变均达到最大值，BF/AAC试件单调及循环应力-应变曲线的峰值应变分别增加了28.12%和28.77%，CF/AAC试件则分别增加了37.17%和41.50%。(3)力学性能指标：①随着纤维掺量越大，BF/AAC和CF/AAC试件的塑性应变均呈小幅度增加趋势。②FAAC的刚度退化过程可分为上升阶段、下降阶段和稳定阶段。③随荷载增加，各组试件的应力退化率先缓慢下降，随后上升，最后进入稳定阶段。(4)应力-应变曲线计算方程：①采用幂函数对FAAC标准化塑性应变与卸载点之间的关系进行拟合。②提出了应力退化率及加、卸载曲线双折线简化模型。③采用过镇海本构模型建立了FAAC单调受压本构模型；基于①和②，建立了循环荷载作用下FAAC的应力-应变曲线计算方程。
结论随纤维掺量增加，FAAC的破坏模式由斜向剪切破坏转向竖向劈裂破坏；随纤维掺量增加，BF/AAC和CF/AC的峰值应力均呈先增加后减小的趋势，峰值应变则呈波动上升趋势，在纤维掺量为0.4%时，BF/AAC和CF/AAC试件的峰值应力均达到最大值。相同纤维掺量条件下，CF对AAC的增强效果优于BF。纤维的掺入会小幅度增加AAC的塑性应变，但纤维掺量对卸载刚度及应力退化率的影响规律不明显。基于试验结果，提出了应力退化率及加、卸载曲线双折线简化模型，建立的FAAC单轴循环受压应力-应变曲线计算方程具有较好的精度。
-
蒸压加气混凝土(AAC)是以硅砂、水泥、石灰为主要原材料并经静停、搅拌、浇筑、蒸养等工艺形成的多孔硅酸盐材料。AAC具有轻质、保温、节能、耐火等优点，但材料自身强度低，脆性大，延性及韧性不足。因此，目前AAC主要被应用于建筑非结构部位；对其研究也多见于配合比、蒸压养护条件等工艺参数优化以及材料干密度、孔结构、导热系数等物理性质指标方面。
前期，课题组提出了以AAC等轻质材料为内填砌块的装配式复合墙体结构，拓宽了AAC的应用途径。课题组前期研究表明，内填砌块材料的物理力学性能对装配式复合墙体的抗震力学性能指标有较大影响。但目前缺乏循环荷载作用下AAC力学行为及其本构模型的研究，这亦是装配式复合墙体非线性分析的前提。
本文在基体材料中掺入玄武岩纤维(BF)、碳纤维(CF)增强材料以改善AAC的承载力、延性及韧性，通过单调加载及循环加载试验得到了纤维增强蒸压加气混凝土(FAAC)的应力-应变曲线，并获取了其峰值应力、峰值应变、累积塑性应变、刚度退化、应力退化率等力学性能指标；基于试验分析结果建立了循环荷载作用下FAAC的应力-应变曲线计算方程。
循环荷载作用下BF/AAC部分力学性能指标 (a)BF-AAC塑性应变 (b)计算模型验证
计量
- 文章访问数: 197
- HTML全文浏览量: 101
- PDF下载量: 18
