Structural Health Monitoring of Filament Wound Pressure Vessel by Embedded Strain Gauges
-
摘要: 在玻璃纤维缠绕金属内胆复合材料压力容器的制备过程中,将应变传感器埋在金属内胆与玻璃纤维/环氧树脂复合材料层之间,得到了具有原位监测功能的纤维缠绕压力容器。对该纤维缠绕压力容器开展水压疲劳和爆破实验。疲劳压力的最大值和最小值分别为25 MPa和2 MPa,最大疲劳周次为5700;打压爆破压力为零到爆破压力,打压速率为2 MPa/s。实验过程中,利用埋入式应变传感器原位监测了压力容器的应变变化,建立了不同载荷作用下纤维缠绕压力容器的应变与受载情况之间的关联。结果表明:采用埋入式应变传感器监测纤维缠绕压力容器的健康状况具有可行性;该方法在保护应变传感器不受外载荷破坏的前提下,原位监测了压力容器在疲劳和爆破实验中的应变变化趋势。Abstract: During the manufacturing process of a filament wound pressure vessel, we embed the strain gauges between the metal tank and glass fiber reinforced epoxy composite layer to obtain the capability of in-situ monitoring . Experiments with a full-scale composite pressure vessel during hydraulic fatigue cycling and pressurization are performed. The maximum and minimum pressures in the fatigue test are set as 25 and 2 MPa, and the maximum cycle number is set as 5700 cycles, respectively. The pressurization speed is set as 2 MPa/s from 0 MPa to busting pressure. The strain of the pressure vessel in the two loading tests is monitored by the embedded strain gauge. The relationship between the stain and the loading conditions of the pressure vessel was thus built. Results show that, by embedding the strain gauges during the processing, it is possible to monitor the health status of the vessel under hydraulic fatigue cycling and pressurization load without hurting the sensors by the external load.
-
盐类矿物是指以矿石或岩石形式存在于自然界的碱金属/碱土金属的卤化物、硫酸盐、碳酸盐、硝酸盐等,被广泛应用于工业和农业领域[1]。该类矿物在高温高压下的结构相变在物理与地球科学领域备受关注,例如:2013 年理论与实验研究发现,典型离子晶体 NaCl 在高压下转变为 Na3Cl和NaCl3[2];Santos等[3]发现,作为地球深部碳循环重要载体的碳酸盐在高温高压下具有10余种结构。硝酸根离子(
NO−3 )和碳酸根离子(CO2−3 )具有等电子结构,形成的盐类矿物均为离子晶体。离子晶体中的离子键易受外部温度和压强调控而发生结构相变,鉴于NO−3 和CO2−3 在结构和电子性质上具有相似之处,高温高压下碱金属硝酸盐是否具有丰富的结构相变是值得深入探讨的问题。碱金属硝酸盐被广泛用于多个领域,如作为传热、储热介质等运用于工业炸药、冶金淬火硬化、太阳能和核能领域[4–6]。碱金属硝酸盐在较低压强下已经展现出丰富的相变行为[7],引起了许多关注。研究发现,在常压下硝酸铷(RbNO3)从熔融态退火至常温的过程中历经了4个不同的结构,即Ⅰ、Ⅱ、Ⅲ、Ⅳ相。Ⅰ相在熔点至564 K的区间内稳定,为NaCl型的立方无序结构(空间群Fm
¯3 m)[8];Ⅱ相在564~492 K区间稳定,具有菱方结构(空间群R¯3 m)[7];Ⅲ相在492~437 K区间稳定,一般认为是CsCl型立方结构(空间群Pm¯3 m)[8];Ⅳ相在常温下稳定,为菱方结构(空间群P31)[7]。2001年,Liu等[9]利用分子动力学模拟方法从理论上发现了RbNO3在Ⅳ→Ⅲ相变过程中NO−3 基团从定向有序变为无序。在高压实验研究方面,1984年Dean等[10]通过原位X射线衍射研究发现了RbNO3的Ⅳ→Ⅴ结构相变,其相变压强约为1.2 GPa,高压Ⅴ相为正交结构(空间群Pmmn)。与RbNO3的高温研究相比,RbNO3的高压研究并不丰富。本研究将基于第一性原理计算,结合CALYPSO晶体结构预测软件,系统地探索零温下RbNO3在0~100 GPa区间的结构与物理性质。
1. 计算方法
在0 K和0、20、40、60、80和100 GPa的压强条件下,使用CALYPSO软件[11–12]进行RbNO3的结构搜索。结构预测过程中演化代数为20,模拟采用2倍胞和4倍胞,每个分子式的体积在32~70 Å3区间,体积随压强的升高而减小。基于第一性原理计算,交换关联泛函采用广义梯度近似(generalized gradient approximation,GGA)[13]的Perdew-Burke-Ernzerhof(PBE)泛函[13]、Original vdW-DF泛函[14]、Revised PBE for solids(PBEsol)泛函[15]以及Revised PBE Hammer et al.(RPBE)泛函[16],数值模拟在VASP[17]软件包中实现,采用投影缀加波(projected augmented wave,PAW)[18]方法,Rb、N、O的价电子分别选为4s24p65s1、2s22p3和2s22p4,平面波截断能取520 eV,对于R3m、Pnma和Pmmn相,布里渊区采样的网格尺寸分别为10×10×10、6×8×6和10×8×10,总能量收敛至小于1×10−5 eV/atom。采用有限位移法结合PHONOPY软件[19]计算声子谱,扩胞数分别为3×3×2、2×3×2和3×3×3。弹性模量计算采用Voigt-Reuss-Hill(VRH)平均法[20–21],晶体结构和声子模式的原子位移可视化通过VESTA软件[22]和V_Sim软件[23]实现。
2. 结果与讨论
RbNO3的常温常压结构为RbNO3-Ⅳ相,空间群为P31。采用4种不同交换关联泛函优化0 GPa下RbNO3-Ⅳ相的平衡体积与晶胞参数,结果列于表1。对比实验结果可以发现,4种交换关联泛函的准确性从低到高依次为RPBE、vdW-DF、 PBE、 PBEsol。因此,后续研究将主要基于PBE和PBEsol泛函展开。
表 1 RbNO3-Ⅳ相的晶格常数的理论值与实验值的比较Table 1. Comparison between experimental and calculated lattice parameters of RbNO3-Ⅳ phaseMethod Condition a/Å c/Å Volume/Å3 Simulation PBE 10.68(+1.62%) 7.60(+1.88%) 867.56(+5.29%) vdW-DF 10.72(+2.00%) 7.62(+2.14%) 876.29(+6.35%) PBEsol 10.39(−1.14%) 7.39(−0.94%) 797.69(−3.19%) RPBE 11.14(+5.99%) 7.91(+6.03%) 982.63(+19.26%) Experiment T=296 K[24] 10.47 7.44 815.58 T=298 K[25] 10.55 7.47 831.43 Average 10.51 7.46 823.96 采用PBE和PBEsol泛函在0~100 GPa压强范围内对能量较低的备选结构进行精细优化后发现,R3m、Pnma、P31、Pmmn相在不同的压强区间内具有热力学稳定性,其焓差随压强变化曲线如图1所示。在0 GPa下,2种泛函的预测结果均表明R3m相的能量最稳定;随压强的升高,R3m相均转变为Pnma相。继续加压,PBE泛函与PBEsol泛函预测的相变序列不同:PBE泛函预测结果表明,Pnma相先相变至Ⅳ相(即P31结构),再相变至Ⅴ相(即Pmmn结构);PBEsol泛函预测结果表明,Pnma相直接转变为Ⅴ相。压强-体积关系如图2所示,其中,ΔV/V为体积坍塌率。除PBE泛函预测的Pnma→Ⅳ相变中不能看到明显的体积坍缩以外,PBE泛函所预测的其余2个相变和PBEsol预测的2个相变均存在明显的体积坍塌,这种体积不连续表明它们均为一级相变。
考虑到PBEsol泛函对实验Ⅳ相的平衡体积与晶胞参数的描述更为准确,后续研究将主要采用PBEsol泛函预测零温下RbNO3的相变序列(R3m→Pnma→Pmmn)。该序列所对应的2个相变点的压强分别为1.7和8.2 GPa,相应的体积坍塌率分别为3.73%和2.54%。由于碳酸钙的结构与碱金属硝酸盐的结构相似,两者具有等电子结构,因此,对比了它们的高压相变序列。碳酸钙的相变序列为:R
¯3 c→Pnma→Pmmn→P21/c,相变压强分别为1、40和67 GPa[26]。可以发现,RbNO3的Pnma→Pmmn的相变压强远低于CaCO3发生相同相变时的压强,其可能的原因如下:一方面,Rb+与NO−3 之间的离子相互作用比Ca2+与CO2−3 之间的离子作用更弱;另一方面,Rb+的半径大于Ca2+的半径,Rb+产生的“化学预压缩”更强。上述零温理论预测出了能量比常温常压下的Ⅳ相更小的2个低压新相(R3m相和Pnma相),表明低压下RbNO3可能存在不同于实验给出的常温常压下Ⅳ相的低温新相。预测获得的0 GPa下的R3m相、4 GPa下的Pnma相和12 GPa下的Pmmn相的平衡态晶格常数、原子位置列于表2,相应的晶体结构见图3。3个结构存在明显差异,三方相的单胞包含3个RbNO3分子式,而正交相的单胞则分别包含4和2个RbNO3分子式。3个相在0、4、12 GPa压强下的弹性常数列于表3,可以发现,三方相和正交相分别有7和9个独立的弹性常数。在R3m相中,C44随压强的升高先减小后增大,其余的弹性常数均与压强呈正相关。3个相的弹性模量,即体积模量B、剪切模量G以及杨氏模量E也列于表3中,分别表示抵抗脆性变形、抵抗塑性变形和抵抗张力的能力。同一压强下Pmmn相的弹性模量最大,表明Pmmn相的抗变形能力比另2个相更强。此外,根据“伯恩-黄昆”弹性稳定性判据[27–28]分析3个相的机械稳定性,判别条件见式(1)和式(2)。力稳定性分析显示:12 GPa时,Pnma相的C55 < 0,不满足弹性稳定性判据,因而Pnma相不是弹性稳定结构;除此之外,各相在相应压强下的弹性常数均满足弹性稳定判据。
表 2 高压下R3m、Pnma和Pmmn相的晶格参数和原子位置Table 2. Lattice parameters and atomic coordinates of R3m, Pnma and Pmmn phases at different pressuresPhase Pressure/GPa Lattice parameters Wyckoff positions R3m 0 a=b=5.64 Å, c=9.48 Å,
α=β=90°, γ=120°Rb1: 3a(0.3333, 0.6667, 0.1699)
N1: 3a(0.3333, 0.6667, 0.7275)
O1: 9b(0.5373, 0.0746, 0.3953)Pnma 4 a=7.18 Å, b=5.61 Å, c=6.77 Å,
α=β=γ=90°Rb1: 4c(0.483, 0.250, 0.686)
N1: 4c(0.350, 0.250, 0.122)
O1: 4c(0.448, 0.250, 0.275)
O5: 8d(0.302, 0.444, 0.043)Pmmn 12 a=4.67 Å, b=5.34 Å, c=4.63 Å,
α=β=γ=90°Rb1: 2b(0.500, 0, 0.400)
N1: 2a(0, 0, 0.992)
O1: 2a(0, 0, 0.721)
O3: 4e(0, 0.204, 0.128)表 3 高压下R3m、Pnma和Pmmn相的弹性常数和弹性模量Table 3. Elastic constant and elastic moduli of R3m, Pnma and Pmmn phase at high pressureGPa Phase Pressure C11 C12 C13 C14 C22 C23 C33 C44 C55 C66 B G E R3m 0 37 14 11 −4 13 7 12 15 7 19 4 75 28 21 −8 39 6 23 35 11 29 12 124 46 42 −15 94 12 39 66 19 52 Pnma 0 22 7 9 41 9 32 7 10 9 16 10 24 4 50 19 30 73 19 69 12 17 19 36 17 44 12 96 43 59 129 35 119 19 −5 34 69 −57 −236 Pmmn 0 26 9 9 50 16 54 15 5 5 21 10 26 4 60 16 25 84 32 83 33 20 9 40 21 53 12 113 25 53 143 59 129 60 37 16 72 35 91 C11>|C12|,C44>0,C213<12C33(C11+C12),C214<12C44(C11−C12) (1) C11>0,C11C22>C212,C11C22C33+2C12C13C23−C11C223−C22C213−C33C212>0C44>0,C55>0,C66>0 (2) 为了确定RbNO3的3个相的动力学稳定性,采用有限位移法计算了它们的声子谱。若声子谱具有负频率,则意味着动力学不稳定[29]。在简谐近似下,分别计算了R3m、Pnma和Pmmn相在不同压强下的声子色散曲线以及声子分波态密度(phonon density of states, PHDOS),结果如图4所示。可以看出,声子谱均未出现负频率,表明3个相在相应压强下均具有动力学稳定性。
声子态密度分析研究表明,RbNO3的3个相的声子态密度分布趋势皆与常压下的CaCO3[30]十分相似,例如:当RbNO3(CaCO3)的振动频率低于10 THz(12 THz)时,声子态密度的贡献主要来自于O原子与金属阳离子Rb+(Ca2+);在振动频率32 THz附近,声子态密度的贡献主要来自于O原子;在振动频率大于40 THz的高频区,声子态密度的贡献来自于N―O键或者C―O键。采用V_Sim软件对Γ点处声子的振动模式进行可视化分析,结果表明:当频率小于10 THz时,N原子与O原子之间的距离不变,Rb和
NO−3 协同振动;在32 THz附近,R3m和Pnma相的N―O振动和C―O振动均为对称伸缩振动,与方解石和文石的特性类似;在21 THz附近,N―O振动和C―O振动均为平面内弯曲振动;在RbNO3(CaCO3)的振动频率约为24 THz(25 THz)时,N―O振动(C―O振动)均为平面外弯曲振动;在大于40 THz的范围内,N―O振动和C―O振动均为反对称弯曲振动。相对于0 GPa的R3m相和4 GPa的Pnma相而言,12 GPa的Pmmn相的声子振动频率整体升高,但不同模式的频率分布范围大致相同。为了研究R3m、Pnma和Pmmn相中原子间的成键行为,计算了0 GPa的R3m相、4 GPa的Pnma相和12 GPa的Pmmn相的电子局域函数(electron localization functionm, ELF),计算结果如图5所示。三维ELF图中等值面的取值为0.8,二维ELF图的晶格切面分别为(1.000 14.005 31.363)、(20.145 2.048 −1.000)和(4.542 1.000 −1.516)。可以看出,N原子与O原子间的ELE值趋近于1,表明N―O之间存在强共价键。
此外,计算了0 GPa的R3m相、4 GPa的Pnma相和12 GPa的Pmmn相的晶体轨道哈密顿布居(crystal orbital Hamiton population, COHP),计算结果如图6所示。COHP为负表示成键态。相对于O―O键和N―O键而言,Rb―O键对COHP的贡献非常小,几乎可以忽略不计,符合Rb+与
NO−3 之间为离子相互作用的事实。在−1.5~0 eV范围内,O―O键的COHP为正值,表明其主要来源为反键轨道,包含NO−3 从Rb+处获取的1个价电子的贡献。在0 eV以下,N―O键的COHP几乎全部为负值,表明NO−3 中的N―O键是极为稳定的强共价键,与ELF的分析结果一致。在0 GPa的R3m相、4 GPa的Pnma相和12 GPa的Pmmn相中,N原子与O原子的键长分别为1.26、1.25和1.26 Å。计算Bader电荷转移以确定Rb原子、N原子和O原子之间的电子转移情况,计算结果如表4所示。这3个结构中,硝酸根离子分别从Rb+获得了0.90e、0.87e和0.86e的电荷,表明压强抑制Rb+的电荷向NO−3 离子转移。表 4 高压下R3m、Pnma和Pmmn相的元胞内的Bader电荷转移Table 4. Calculated Bader charges within R3m, Pnma and Pmmn phases primitive cells at different pressuresPhase Pressure/GPa Atoms Number Charge/e Charge transfer/e R3m 0 Rb 1 8.10 0.90 N 1 4.16 0.84 O 3 6.58 −0.58 Pnma 4 Rb 4 8.13 0.87 N 4 4.13 0.87 O1 4 6.57 −0.57 O2 4 6.58 −0.58 O3 4 6.59 −0.59 Pmmn 12 Rb 2 8.14 0.86 N 2 4.07 0.93 O1 2 6.58 −0.58 O2 2 6.59 −0.59 O3 2 6.62 −0.62 图7给出了0 GPa的R3m相、4 GPa的Pnma相和12 GPa的Pmmn相的电子能带结构和电子态密度( density of states, DOS)。电子态密度结果表明,价带主要来自O的p轨道,导带主要来自N-p与O-p的杂化。3个结构在相应压强下的带隙分别为3.06、3.40和3.10 eV。带隙的存在说明它们均为半导体材料,然而3个结构的带隙类型不同,其中,R3m结构具有直接带隙,Pnma和Pmmn结构具有间接带隙。此外,还可以发现,在结构相变过程中,带隙随压强的增加呈现非单调变化的规律,即第1次相变带隙增大,而第2次相变带隙减小。
3. 结 论
采用第一性原理计算与CALYPSO结构预测相结合的方法系统地研究了RbNO3在0~100 GPa压强范围内的高压结构与物理性质。以RbNO3常温常压结构的实验数据为参照,测试了PBE、vdW-DF、RPBE、PBEsol泛函对结构描述的准确性,发现PBEsol泛函给出的结果最准确。基于PBEsol泛函预测出零温下RbNO3的高压相变序列为R3m→Pnma→Pmmn,相变压强分别约为1.7 和8.2 GPa;2次相变均为一级结构相变,相应的体积坍塌率分别为3.73%和2.54%。预测零温下RbNO3存在2个低压新相,即R3m相和Pnma相,它们在能量上均优于早期实验在常温常压下给出的P31结构,表明降温可能会诱导RbNO3发生结构相变。0 GPa的R3m相、4 GPa的Pnma相和12 GPa的Pmmn相均具有力学稳定性和动力学稳定性,其带隙分别为3.06、3.40和3.10 eV,对比分析发现,压强会抑制Rb+的电荷向
NO−3 转移。本研究获得的硝酸铷结构与电子性质随压强的变化规律可以为其他碱金属硝酸盐的理论与实验研究提供参考。 -
图 1 具有原位监测功能的复合材料压力容器制备过程:(a)纤维缠绕工艺,(b)粘贴在金属内胆的应变片,(c)带有埋入式应变片的纤维缠绕压力容器
Figure 1. Preparation of the filament wound pressure vessel with the embedded strain gauges: (a) the fiber winding processing; (b) the strain gauges attached on the inner steel tank; (c) the pressure vessel with the embedded stain gauges
表 1 纤维缠绕压力容器中金属内胆参数
Table 1. Parameters of the inner metal tank in the filament wound pressure vessel
Material External
diameter/mmWall thickness
/mmHeight
/mmWeight
/kgVolume
/LWork pressure
/MPaTop pressure
/MPa30CrMo steel 325 5 750 45.3 55 20 34 表 2 纤维缠绕压力容器打压疲劳试验参数
Table 2. Parameters of the hydraulic fatigue cycling
Media Cycling rate Minimum pressure/MPa Rising pressure time/s Minimum pressure holding time/s Water 7.8 2 3.1 2 Temperature/℃ Cycle times Maximum pressure/MPa Pressure drop time/s Maximum pressure holding time/s 11.2 5700 25 0.6 2 表 3 不同纤维缠绕压力容器在不同压力下的最大应变
Table 3. Maximum strain of the pressure vessel at various pressures
Pressure/MPa Maximum strain Without fatigue and groove Without fatigue, with groove Fatigue, without groove Fatigue, with groove 5 –120 –100 –400 –1000 10 –160 –150 –1200 –1200 15 –300 –300 –1200 –1200 20 –500 –500 –1500 –1500 25 –600 –650 -
[1] 缑林虎, 郑锡涛, 程勇. 平面缠绕炭纤维压力容器大变形有限元分析 [J]. 固体火箭技术, 2010, 33(2): 205–208. doi: 10.3969/j.issn.1006-2793.2010.02.019GOU L H, ZHENG X T, CHENG Y. Large deformation finite element analysis of planar carbon fiber wound composite pressure vessel [J]. Journal of Solid Rocket Technology, 2010, 33(2): 205–208. doi: 10.3969/j.issn.1006-2793.2010.02.019 [2] 郑津洋, 开方明, 刘仲强, 等. 轻质高压储氢容器 [J]. 化工学报, 2004, 55(Suppl 1): 130–133.ZHENG J Y, KAI F M, LIU Z Q, et al. Lightweight high-pressure hydrogen tank [J]. Journal of Industry and Engineering, 2004, 55(Suppl 1): 130–133. [3] 杨斌, 章继峰, 梁文彦, 等. 玻璃纤维表面纳米SiO改性对GF/PCBT复合材料力学性能的影响 [J]. 复合材料学报, 2015, 32(3): 691–698.YANG B, ZHANG J F, LIANG W Y, et al. Effects of glass fiber surface modified by nano-SiO2 on mechanical properties of GF/PCBT composites [J]. Acta Metallurgica Sinica, 2015, 32(3): 691–698. [4] 杨斌, 章继峰, 周利民. 玻璃纤维-碳纤维混杂增强PCBT复合材料层合板的制备及低速冲击性能 [J]. 复合材料学报, 2015, 32(2): 435–443.YANG B, ZHANG J F, ZHOU L M. Preparation and low-velocity impact properties of glass fiber-carbon fiber hybrid reinforced PCBT composite laminate [J]. Acta Metallurgica Sinica, 2015, 32(2): 435–443. [5] 路智敏, 李强, 李卓. 基于爆破试验的CFRP固体火箭发动机壳体的可靠性设计 [J]. 复合材料学报, 2009, 26(2): 176–180. doi: 10.3321/j.issn:1000-3851.2009.02.031LU Z M, LI Q, LI Z. Reliability design of CFRP solid rocket motor vessel based on the burst experiment [J]. Acta Metallurgica Sinica, 2009, 26(2): 176–180. doi: 10.3321/j.issn:1000-3851.2009.02.031 [6] OZEVIN D, HARDING J. Novel leak localization in pressurized pipeline networks using acoustic emission and geometric connectivity [J]. International Journal of Pressure Vessels and Piping, 2012, 92: 63–69. doi: 10.1016/j.ijpvp.2012.01.001 [7] CHOU H Y, MOURITZ A P, BANNISTER M K, et al. Acoustic emission analysis of composite pressure vessels under constant and cyclic pressure [J]. Composites Part A: Applied Science and Manufacturing, 2015, 70: 111–120. doi: 10.1016/j.compositesa.2014.11.027 [8] KHAN A, KO D K, LIM S C, et al. Structural vibration-based classification and prediction of delamination in smart composite laminates using deep learning neural network [J]. Composites Part B: Engineering, 2019, 161: 586–594. doi: 10.1016/j.compositesb.2018.12.118 [9] 王晓勇, 熊建平, 高义广. X射线切线照相检测技术在纤维缠绕压力容器检测中的应用 [J]. 航天制造技术, 2011(6): 65–68.WANG X Y, XIONG J P, GAO Y G. Application of X-ray inspection technique in detection of filament-wound pressure vessel [J]. Aerospace Manufacturing Technology, 2011(6): 65–68. [10] 杜善义, 冷劲松, 顾震隆. 用应力波技术对配橡胶内衬的复合材料板壳进行无损检测 [J]. 复合材料学报, 1993, 10(1): 65–69.DU S Y, LENG J S, GU Z L. Non-destructive testing for composite plate and shell with rubber liner using stress wave technique [J]. Acta Metallurgica Sinica, 1993, 10(1): 65–69. [11] 乔业程, 王福强. 压力容器氢损伤的监测与检测方法 [J]. 橡塑技术与装备, 2018, 44(20): 54–56.QIAO Y C, WANG F Q. Monitoring and detection of hydrogen damage in pressure vessels [J]. China Rubber/Plastics Technology and Equipment, 2018, 44(20): 54–56. [12] 赵海涛, 张博明, 武湛君, 等. 纤维缠绕复合材料压力容器健康监测研究进展 [J]. 压力容器, 2007, 24(3): 48–61. doi: 10.3969/j.issn.1001-4837.2007.03.012ZHAO H T, ZHANG B M, WU Z J, et al. Development of health monitoring for filament wound composite pressure vessels [J]. Pressure Vessel Technology, 2007, 24(3): 48–61. doi: 10.3969/j.issn.1001-4837.2007.03.012 [13] BELLAN F, BULLETTI A, CAPINERI L, et al. A new design and manufacturing process for embedded Lamb waves interdigital transducers based on piezopolymer film [J]. Sensors and Actuators A, 2005, 123: 379–387. [14] AI D, ZHU H, LUO H. Sensitivity of embedded active PZT sensor for concrete structural impact damage detection [J]. Construction and Building Materials, 2016, 111: 348–357. doi: 10.1016/j.conbuildmat.2016.02.094 [15] ANNAMDAS V G M, SOH C K. Embedded piezoelectric ceramic transducers in sandwiched beams [J]. Smart Materials and Structures, 2006, 15(2): 538–549. doi: 10.1088/0964-1726/15/2/037 [16] DZIENDZIKOWSKI M, KURNYTA A, DRAGAN K, et al. In situ barely visible impact damage detection and localization for composite structures using surface mounted and embedded PZT transducers: a comparative study [J]. Mechanical Systems and Signal Processing, 2016, 78: 91–106. doi: 10.1016/j.ymssp.2015.09.021 [17] GHIMIRE M, WANG C, DIXON K, et al. In situ monitoring of prestressed concrete using embedded fiber loop ringdown strain sensor [J]. Measurement, 2018, 124: 224–232. doi: 10.1016/j.measurement.2018.04.017 [18] WANG Y, WANG Y, HAN B, et al. Strain monitoring of concrete components using embedded carbon nanofibers/epoxy sensors [J]. Construction and Building Materials, 2018, 186: 367–378. doi: 10.1016/j.conbuildmat.2018.07.147 [19] CHOWDHURY N T, JOOSTEN M W, PEARCE G M K. An embedded meshing technique (SET) for analysing local strain distributions in textile composites [J]. Composite Structures, 2019, 210: 294–309. doi: 10.1016/j.compstruct.2018.11.026 [20] KANERVA M, ANTUNES P, SARLIN E, et al. Direct measurement of residual strains in CFRP-tungsten hybrids using embedded strain gauges [J]. Materials & Design, 2017, 127: 352–363. -