留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

纤维增强复合板中声弹Lamb波的波结构分析

刘宏业 刘申 吕炎 何存富 周丽红

刘宏业, 刘申, 吕炎, 何存富, 周丽红. 纤维增强复合板中声弹Lamb波的波结构分析[J]. 工程力学. doi: 10.6052/j.issn.1000-4750.2019.09.0565
引用本文: 刘宏业, 刘申, 吕炎, 何存富, 周丽红. 纤维增强复合板中声弹Lamb波的波结构分析[J]. 工程力学. doi: 10.6052/j.issn.1000-4750.2019.09.0565
Hongye LIU, Shen LIU, Yan LYU, Cunfu He, Lihong ZHOU. WAVE STRUCTURE ANALYSIS OF ACOUSTOELASTIC LAMB WAVES IN FIBER REINFORCED COMPOSITE LAMINA[J]. Engineering Mechanics. doi: 10.6052/j.issn.1000-4750.2019.09.0565
Citation: Hongye LIU, Shen LIU, Yan LYU, Cunfu He, Lihong ZHOU. WAVE STRUCTURE ANALYSIS OF ACOUSTOELASTIC LAMB WAVES IN FIBER REINFORCED COMPOSITE LAMINA[J]. Engineering Mechanics. doi: 10.6052/j.issn.1000-4750.2019.09.0565

纤维增强复合板中声弹Lamb波的波结构分析

doi: 10.6052/j.issn.1000-4750.2019.09.0565
基金项目: 国家自然科学基金项目(51705325,11872082,11527801)
详细信息
    作者简介:

    刘 申(1996−),男,安徽人,硕士生,主要从事超声无损检测研究(E-mail: liushen223@foxmail.com)

    吕 炎(1985−),男,河南人,副教授,博士,硕导,主要从事智能传感器与结构健康监测研究(E-mail: lvyan@bjut.edu.cn)

    何存富(1958−),男,山西人,教授,博士,博导,主要从事超声无损检测新技术研究(E-mail: hecunfu@bjut.edu.cn)

    通讯作者: 刘宏业(1986−),男,山东人,讲师,博士,硕导,主要从事超声无损检测与结构健康监测研究(E-mail: liuhongye@usst.edu.cn)
  • 中图分类号: TB332

WAVE STRUCTURE ANALYSIS OF ACOUSTOELASTIC LAMB WAVES IN FIBER REINFORCED COMPOSITE LAMINA

  • 摘要: 基于线性三维弹性理论和“增量变形力学”理论,采用勒让德正交多项式展开法,推导了在水平和垂直方向施加初应力时,沿纤维增强复合板的非主对称轴方向传播的声弹Lamb波的波动方程,并对波动方程进行数值求解。为了验证方法的准确性,首先将本文方法求解的各向同性材料的相速度频散曲线与Disperse®软件的计算结果进行比较,两者吻合良好。然后以单层单向纤维增强复合材料板为例,计算了无初应力状态下的波结构应力曲线,并与应力自由边界初始条件的一致性进行了比较。最后,研究了水平和垂直方向初应力效应对Lamb波频散曲线的影响。针对声弹效应较为敏感的Lamb波A0模态,着重分析了初应力效应对波结构位移分布曲线的影响。
  • 图  1  初应力作用下纤维增强复合板空间坐标示意图

    Figure  1.  Schematic diagram of space coordinates for fiber reinforced composite plate under initial stress.

    图  2  本文方法所计算的铝板相速度频散曲线与Disperse®软件计算结果对比图

    Figure  2.  Comparison of phase velocity dispersion curves for aluminum plate by the authors’ program with available data calculated from Disperse®.

    图  3  300 kHz频率下Lamb波A0模态的应力波结构曲线

    Figure  3.  Stress distributions of Lamb wave A0 mode at 300 kHz.

    图  4  截止项αΓ由3变化到7时各模态频散曲线的收敛

    Figure  4.  Convergence of dispersion curves of all the modes with the change of cut-off term αΓ from 3 to 7.

    图  5  S0模态频散曲线收敛过程的放大图

    Figure  5.  Zoom-in graph of convergence process of the S0 mode.

    图  6  ±1 GPa初应力作用下,单层单向纤维增强复合材料板中沿0°方向传播的Lamb波相速度频散曲线

    Figure  6.  Phase velocity dispersion curves of Lamb waves propagating along 0° direction in fiber reinforced composite lamina under the initial stresses of ±1 GPa.

    图  7  传播方向为15°时,频率300 kHz的A0模态在不同水平初应力作用下的位移分布

    Figure  7.  Displacement distributions of A0 mode at 300 kHz along 15° propagation direction under different horizontal initial stresses

    图  8  传播方向为45°时,频率300 kHz的A0模态在不同水平初应力作用下的位移分布

    Figure  8.  Displacement distributions of A0 mode at 300 kHz along 45° propagation direction under different horizontal initial stresses

    图  9  传播方向为75°时,频率300 kHz的A0模态在不同水平初应力作用下的位移分布

    Figure  9.  Displacement distributions of A0 mode at 300 kHz along 75° propagation direction under different horizontal initial stresses

    图  10  ±1 GPa垂直初应力作用下,单层单向纤维增强复合材料板中沿0°方向传播的Lamb波相速度频散曲线

    Figure  10.  Phase velocity dispersion curves of Lamb waves propagating along 0° direction in fiber reinforced composite lamina under the vertical initial stresses of ±1 GPa.

    图  11  传播方向为15°时,频率300 kHz的A0模态在不同垂直初应力作用下的位移分布

    Figure  11.  Displacement distributions of A0 mode at 300 kHz along 15° propagation direction under different vertical initial stresses

    图  12  传播方向为45°时,频率300 kHz的A0模态在不同垂直初应力作用下的位移分布

    Figure  12.  Displacement distributions of A0 mode at 300 kHz along 45° propagation direction under different vertical initial stresses

    图  13  传播方向为75°时,频率300 kHz的A0模态在不同垂直初应力作用下的位移分布

    Figure  13.  Displacement distributions of A0 mode at 300 kHz along 75° propagation direction under different vertical initial stresses

    表  1  碳纤维增强复合材料T300/914的弹性常数

    Table  1.   Elastic constants of carbon-fiber reinforced composites T300/914.

    单位:cij(GPa)
    C11C12C13C22C23C33C44C55C66
    143.86.26.213.36.513.33.65.75.7
    下载: 导出CSV
  • [1] 邢丽英, 包建文, 礼嵩明, 等. 先进树脂基复合材料发展现状和面临的挑战[J]. 复合材料学报, 2016, 33(7): 1327 − 1338.

    Xing Liying, Bao Jianwen, Li Songming, et al. Development status and facing challenge of advanced polymer matrix composites [J]. Acta Materiae Compositae Sinica, 2016, 33(7): 1327 − 1338. (in Chinese)
    [2] Deng M, Yang J. Characterization of elastic anisotropy of a solid plate using nonlinear Lamb wave approach [J]. Journal of Sound and Vibration, 2007, 308(1): 201 − 211.
    [3] Cao X, Shi J, Jin F. Cut-off frequencies of Lamb waves in various functionally graded thin films [J]. Applied Physics Letters, 2011, 99(12): 121907(1) − 121907(3).
    [4] 郑阳, 何存富, 周进节, 等. 超声Lamb波在缺陷处的二维散射特性研究[J]. 工程力学, 2013, 30(8): 236 − 243.

    Zheng Yang, He Cunfu, Zhou Jinjie, et al. 2Dscattering characteristics of Lamb wave at defects [J]. Engineering Mechanics, 2013, 30(8): 236 − 243. (in Chinese)
    [5] 刘增华, 徐营赞, 何存富, 等. 板状结构中基于Lamb波单模态的缺陷成像试验研究[J]. 工程力学, 2014, 31(4): 232 − 238.

    Liu Zenghua, Xu Yingzan, He Cunfu, et al. Experimental study on defect imaging based on single LAMB wave mode in plate-like structures [J]. Engineering Mechanics, 2014, 31(4): 232 − 238. (in Chinese)
    [6] Pei N, Bond L J. Higher order acoustoelastic Lamb wave propagation in stressed plates [J]. Journal of the Acoustical Society of America, 2016, 140(5): 3834 − 3843. doi:  10.1121/1.4967756
    [7] 周正干, 冯占英, 高翌飞, 等. 超声导波在大型薄铝板缺陷检测中的应用[J]. 航空学报, 2008, 29(4): 1044 − 1048. doi:  10.3321/j.issn:1000-6893.2008.04.045

    Zhou Zhenggan, Feng Zhanying, Gao Yifei, et al. Application of Ultrasonic Guided Waves to Defect Inspection of Large Thin Aluminum Plate [J]. Acta Aeronautica Et Astronautica Sinica, 2008, 29(4): 1044 − 1048. (in Chinese) doi:  10.3321/j.issn:1000-6893.2008.04.045
    [8] 张海燕, 刘凡杰, 范国鹏, 等. 各向同性板中盲孔缺陷的兰姆波拓扑成像[J]. 声学学报, 2018, 43(6): 968 − 976.

    Zhang Haiyan, Liu Fanjie, Fan Guopeng, et al. Lamb wave topological imaging for blind hole defects in isotropic plates [J]. Acta Sinica, 2018, 43(6): 968 − 976. (in Chinese)
    [9] 曹小彬, 金峰, 王子昆. 功能梯度材料板中Lamb波传播特性研究[J]. 固体力学学报, 2009, 30(1): 35 − 41.

    Cao Xiaoshan, Jin Feng, Wang Zikun. Propagation properties of Lamb waves in a functionally graded material plate [J]. Chinese Journal of Solid Mechanics, 2009, 30(1): 35 − 41. (in Chinese)
    [10] 王国锋, 李富才, 刘志强, 等. 超声导波在圆管结构损伤定位中的应用[J]. 振动、测试与诊断, 2017, 37(3): 440 − 448.

    Wang Guofeng, Li Fucai, Liu Zhiqiang, et al. Application of ultrasonic guided wave in damage localization of circular tube structure [J]. Journal of Vibration, Measurement & Diagnosis, 2017, 37(3): 440 − 448. (in Chinese)
    [11] NG C T, Veidt M. Scattering analysis of fundamental anti-symmetric Lamb wave at delaminations in composite laminates [J]. Australian Journal of Mechanical Engineering, 2011(129): 1288 − 1296.
    [12] 张海燕, 刘镇清, 马小松. 界面层对层状各向异性复合结构中Lamb波的影响[J]. 物理学报, 2003(10): 2492 − 2499. doi:  10.3321/j.issn:1000-3290.2003.10.026

    Zhang Haiyan, Liu Zhenqing, Ma Xiaosong. The influences of interface layer characteristics on Lamb waves in layered anisotropic media [J]. Acta Physica Sinica, 2003(10): 2492 − 2499. (in Chinese) doi:  10.3321/j.issn:1000-3290.2003.10.026
    [13] 徐春铃, 王鑫伟. Lamb波在复合材料板中传播的谱有限元模拟[J]. 力学季刊, 2011, 32(1): 10 − 18.

    Xu Chunling, Wang Xinwei. Finite element simulation of Lamb wave propagation in composite plate [J]. Chinese Quarterly of Mechanics, 2011, 32(1): 10 − 18. (in Chinese)
    [14] 何存富, 刘宏业, 刘增华, 等. 正交各向异性板中非主轴方向的Lamb波[J]. 固体力学学报, 2013, 34(1): 55 − 62. doi:  10.3969/j.issn.0254-7805.2013.01.008

    He Cunfu, Liu Hongye, Liu Zenghua, et al. Lamb waves propagation along non-principal directions in orthotropic plate [J]. Chinese Journal of Solid Mechanics, 2013, 34(1): 55 − 62. (in Chinese) doi:  10.3969/j.issn.0254-7805.2013.01.008
    [15] He C F, Liu H Y, Liu Z H, et al. The propagation of coupled Lamb waves in multilayered arbitrary anisotropic composite laminates [J]. Journal of Sound and Vibration, 2013, 332: 7243 − 7256. doi:  10.1016/j.jsv.2013.08.035
    [16] Du J K, Jin X Y, Wang J. Love wave propagation in layered magneto-electro-elastic structures with initial stress [J]. Acta Mechanica, 2007, 192: 169 − 189. doi:  10.1007/s00707-006-0435-3
    [17] Navneet G, Jennifer E M, Sang J L. Aoustoelastic Lamb wave propagation in biaxially stressed plates [J]. Journal of the Acoustical Society of America, 2012, 132(3): 1284 − 1293. doi:  10.1121/1.4740491
    [18] Yu J G, Zhang Ch. Influences of initial stresses on guided waves in functionally graded hollow cylinders [J]. Acta Mechanica, 2013, 224(4): 745 − 757. doi:  10.1007/s00707-012-0748-3
    [19] Yu J G, Li S L. Dispersion of guided wave in initially stressed layered plates [J]. Mechanics of materials and structures, 2013, 8: 2 − 4.
    [20] Lefebvre J. E., Zhang V., Gazalet J., et al. Acoustic wave propagation in continuous functionally graded plates: an extension of the Legendre polynomial approach [J]. IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, 2001, 48: 1332 − 1339. doi:  10.1109/58.949742
  • [1] 张春晓, 何翔, 刘国权, 王世合, 李磊, 李砚召.  泡沫陶瓷球壳与高黏弹沥青热压复合板材研制及其抗爆性能试验研究 . 工程力学, doi: 10.6052/j.issn.1000-4750.2016.03.S029
    [2] 黄天朗, 宫凤强, 伍婷玉.  小样本岩土参数概率分布的正交多项式推断方法 . 工程力学, doi: 10.6052/j.issn.1000-4750.2016.03.0169
    [3] 张宇, 闫云聚, 余龙, 袭著有.  CSM在基于Lamb波的结构损伤检测中的应用 . 工程力学, doi: 10.6052/j.issn.1000-4750.2013.12.1159
    [4] 许斌, 辛璐璐, 贺佳.  基于切比雪夫多项式模型的多自由度结构非线性恢复力时域识别 . 工程力学, doi: 10.6052/j.issn.1000-4750.2013.05.0428
    [5] 胡剑虹, 唐志峰, 吕福在, 姜晓勇.  基于SAFE的层状结构导波传播特性研究 . 工程力学, doi: 10.6052/j.issn.1000-4750.2013.05.0457
    [6] 何志全, 周丽, 孙虎.  基于传递阻抗能量的无基准Lamb波裂纹检测 . 工程力学, doi: 10.6052/j.issn.1000-4750.2013.01.0113
    [7] 刘增华, 徐营赞, 何存富, 吴斌.  板状结构中基于Lamb波单模态的缺陷成像试验研究 . 工程力学, doi: 10.6052/j.issn.1000-4750.2012.11.0842
    [8] 孔艳平, 岳鹏阳, 刘金喜.  正交各向异性弹性层/压磁半空间结构中Love波的传播 . 工程力学, doi: 10.6052/j.issn.1000-4750.2012.08.0608
    [9] 张小明, 王裕清, 禹建功.  粘弹正交各向异性空心圆柱体周向SH波 . 工程力学, doi: 10.6052/j.issn.1000-4750.2011.11.0736
    [10] 郑阳, 何存富, 周进节, 张也弛.  超声Lamb波在缺陷处的二维散射特性研究 . 工程力学, doi: 10.6052/j.issn.1000-4750.2012.03.0181
    [11] 骞朋波 尹晓春 沈煜年 孔德平.  梁撞击弹塑性波传播的动态子结构方法的研究 . 工程力学, doi: 10.6052/j.issn.1000-4750.2011.03.0165
    [12] 孙虎, 周丽.  基于谱元法的裂Lamb纹梁波传播特性研究 . 工程力学, doi: 10.6052/j.issn.1000-4750.2010.12.0921
    [13] 林晓平, 刘增华, 雷振坤, 武湛君.  水下带粘弹性层输油管道中纵向导波的传播特性 . 工程力学,
    [14] 乔红威, 吕震宙, 关爱锐, 刘旭华.  平稳随机激励下随机结构动力可靠度分析的多项式逼近法 . 工程力学,
    [15] 刘亚男, 肖亚飞, 刘应华, 岑章志.  基于DB小波无网格法:二维弹塑性问题 . 工程力学,
    [16] 于国友, 罗晓健, 苗 青.  复杂结构损伤检测的分段小波变换法 . 工程力学,
    [17] 宫凤强, 李夕兵.  基于Legendre正交多项式逼近法的结构可靠性分析 . 工程力学,
    [18] 吴 斌, 禹建功, 何存富.  正交各向异性空心圆柱体中的周向超声导波 . 工程力学,
    [19] 邱颖, 任青文, 朱建华.  基于小波奇异性的梁结构损伤诊断 . 工程力学,
    [20] 鞠彦忠, 阎贵平, 陈建斌, 陈景彦, 陈建华.  用小波神经网络检测结构损伤 . 工程力学,
  • 加载中
计量
  • 文章访问数:  55
  • HTML全文浏览量:  60
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-09-25
  • 修回日期:  2020-02-19
  • 网络出版日期:  2020-05-21

纤维增强复合板中声弹Lamb波的波结构分析

doi: 10.6052/j.issn.1000-4750.2019.09.0565
    基金项目:  国家自然科学基金项目(51705325,11872082,11527801)
    作者简介:

    刘 申(1996−),男,安徽人,硕士生,主要从事超声无损检测研究(E-mail: liushen223@foxmail.com)

    吕 炎(1985−),男,河南人,副教授,博士,硕导,主要从事智能传感器与结构健康监测研究(E-mail: lvyan@bjut.edu.cn)

    何存富(1958−),男,山西人,教授,博士,博导,主要从事超声无损检测新技术研究(E-mail: hecunfu@bjut.edu.cn)

    通讯作者: 刘宏业(1986−),男,山东人,讲师,博士,硕导,主要从事超声无损检测与结构健康监测研究(E-mail: liuhongye@usst.edu.cn)
  • 中图分类号: TB332

摘要: 基于线性三维弹性理论和“增量变形力学”理论,采用勒让德正交多项式展开法,推导了在水平和垂直方向施加初应力时,沿纤维增强复合板的非主对称轴方向传播的声弹Lamb波的波动方程,并对波动方程进行数值求解。为了验证方法的准确性,首先将本文方法求解的各向同性材料的相速度频散曲线与Disperse®软件的计算结果进行比较,两者吻合良好。然后以单层单向纤维增强复合材料板为例,计算了无初应力状态下的波结构应力曲线,并与应力自由边界初始条件的一致性进行了比较。最后,研究了水平和垂直方向初应力效应对Lamb波频散曲线的影响。针对声弹效应较为敏感的Lamb波A0模态,着重分析了初应力效应对波结构位移分布曲线的影响。

English Abstract

刘宏业, 刘申, 吕炎, 何存富, 周丽红. 纤维增强复合板中声弹Lamb波的波结构分析[J]. 工程力学. doi: 10.6052/j.issn.1000-4750.2019.09.0565
引用本文: 刘宏业, 刘申, 吕炎, 何存富, 周丽红. 纤维增强复合板中声弹Lamb波的波结构分析[J]. 工程力学. doi: 10.6052/j.issn.1000-4750.2019.09.0565
Hongye LIU, Shen LIU, Yan LYU, Cunfu He, Lihong ZHOU. WAVE STRUCTURE ANALYSIS OF ACOUSTOELASTIC LAMB WAVES IN FIBER REINFORCED COMPOSITE LAMINA[J]. Engineering Mechanics. doi: 10.6052/j.issn.1000-4750.2019.09.0565
Citation: Hongye LIU, Shen LIU, Yan LYU, Cunfu He, Lihong ZHOU. WAVE STRUCTURE ANALYSIS OF ACOUSTOELASTIC LAMB WAVES IN FIBER REINFORCED COMPOSITE LAMINA[J]. Engineering Mechanics. doi: 10.6052/j.issn.1000-4750.2019.09.0565
  • 纤维增强复合材料因其优异的力学性能,近些年被广泛应用于航空航天、军工、汽车、土木建筑等重要领域[1-2]。由于纤维增强复合材料不同组分之间固有的热学性能差异,复合材料部件自身固化成型工艺的限制和服役过程中环境、工况载荷的变化,使得该类材料结构中不可避免的会产生初应力。为实现纤维增强复合材料的出厂质量控制和结构健康监测,寻找一种可靠的无损检测方法至关重要。Lamb波因其传播距离长、检测快速、经济高效等优势被广泛认为是适用于纤维增强复合材料无损检测和结构健康监测的重要技术之一[3-6]。因此,研究初应力效应对纤维增强复合材料中潜在Lamb波传播模态的影响具有十分重要的意义。

    在利用Lamb波对各向异性复合材料进行无损检测和结构健康监测时,通常采用分布式传感器阵列通过多点多传播路径对材料中的缺陷进行定位和成像。对于Lamb波在各向同性材料中的传播问题,由于各Lamb波模态传播特征的各向一致性,往往不涉及“沿非主对称轴传播”的问题。周正干等[7]推导了Lamb波在铝板中传播的理论模型,通过对铝板中Lamb波波结构的分析实现了薄铝板中缺陷的实验检测。张海燕等[8]将拓扑成像用于各向同性板中盲孔缺陷的兰姆波检测,提出的成像方法对板类结构中非对称缺陷的Lamb波检测具有很强的适应能力。曹小杉等[9]采用幂级数法研究了材料性能参数沿厚度方向连续变化的横观各向同性板中Lamb波的传播问题。王国峰[10]等通过利用数值法对频散方程的求解,研究了圆管结构中纵向模态和周向Lamb波的频散曲线,对圆管中导波的传播机理及特征进行了分析。

    但对于各向异性的纤维增强复合材料而言,Lamb波多沿其非主对称轴方向传播。Ng等[11]采用三维有限元法对复合材料板中A0模态的频散特性进行了分析。张海燕等[12]通过建立各向异性界面层的弹簧模型,研究了不同界面条件对层状各向异性复合板中Lamb波频散特性以及位移分布的影响。徐春铃和王鑫伟[13]采用谱有限元法,对复合材料板中对称模式与反对称模式下的Lamb波传播速度进行了求解,数值结果表明较好地模拟出Lamb波在复合材料板结构中的传播特性。何存富等[14-15]关注沿非主对称轴方向传播的Lamb波,研究了单层与多层正交各向异性板中耦合Lamb波的传播特性。

    考虑初应力的影响,“声弹”描述了弹性体中声波波速随初应力效应而改变的现象。Du等[16]考虑在初应力作用下压电材料中Love波的传播问题,通过对Love波的相速度进行数值计算,分析了初应力对Love波相速度的影响,研究结果亦表明初应力作用对波的传播速度有重要影响。Navneet等[17]发展了声弹Lamb波的理论,研究了双轴应力场作用下初应力效应对各向同性板中声弹Lamb波传播速度的影响。Yu和Zhang[18]利用“增量变形力学”理论,研究在初应力作用下空心圆柱体中导波的传播特性,发现初应力对扭转波和纵向波的频散特性和波结构特性的影响不同。同时,Yu和Zhang[19]研究了由Si3N4和Zinc构成的功能梯度板中类Lamb波和SH波的传播问题,发现重力效应的影响可以忽略不计。

    本文以线性三维弹性理论和“增量变形力学”理论为基础,利用勒让德正交多项式展开法对单层单向纤维增强复合材料板中沿非主对称轴传播的声弹Lamb波展开研究。通过分别考虑水平和垂直方向初应力的作用,研究了纤维铺层方向的改变对Lamb频散特性的影响。最后,着重对声弹Lamb波的波结构特性进行了分析。

    • 基于线性三维弹性体理论,图1所示为水平方向无限长,上下表面均为自由表面的纤维增强复合材料板,对于任意纤维铺层方向的复合材料,其本构关系可以看作单斜晶系材料予以表示。其中,O点是纤维增强复合材料板上表面任意一点,以此作为坐标原点,(x1´,x2´,x3´)$ ({x}_{1}^{\prime}{\text{,}}{x}_{2}^{\prime}{\text{,}}{x}_{3}^{\prime}) $是以材料晶轴方向建立的局部坐标系,(x1x2x3)$ ({x}_{1}{\text{,}}{x}_{2}{\text{,}}{x}_{3}) $是把(x1´,x2´,x3´)$ ({{x}}_{1}^{\prime}{\text{,}}{{x}}_{2}^{\prime}{\text{,}}{{x}}_{3}^{\prime}) $坐标系绕x3´轴顺时针旋转$\varphi $角度得到的全局坐标系,复合板厚度为h0。Lamb波沿x1轴方向传播。假定施加的水平初应力表示为${S_{11}} = - P$,垂直初应力表示为${S_{33}} = - Q$

      图  1  初应力作用下纤维增强复合板空间坐标示意图

      Figure 1.  Schematic diagram of space coordinates for fiber reinforced composite plate under initial stress.

      在全局坐标系(x1x2x3)中,材料的弹性常数可以由局部坐标系(x1´,x2´,x3´)下的弹性常数根据式(1)导出,其中${\gamma _{ij}}$为两坐标系下坐标轴xi´和xj之间的方向余弦。

      $${C_{mnop}} = {\gamma _{mi}}{\gamma _{nj}}{\gamma _{ok}}{\gamma _{pl}}{C'_{ijkl}}$$ (1)

      O-x1x2x3坐标系下,单斜晶系材料的本构方程可表示为:

      $$\begin{split}& {T_{11}} = {C_{11}}{\varepsilon _{11}} + {C_{12}}{\varepsilon _{22}} + {C_{13}}\varepsilon _{33}^{} + 2{C_{16}}{\varepsilon _{12}} \\& {T_{22}} = {C_{12}}{\varepsilon _{11}} + {C_{22}}{\varepsilon _{22}} + {C_{23}}\varepsilon _{33}^{} + 2{C_{26}}{\varepsilon _{12}} \\& {T_{33}} = {C_{13}}{\varepsilon _{11}} + {C_{23}}{\varepsilon _{22}} + {C_{33}}\varepsilon _{33}^{} + 2{C_{36}}{\varepsilon _{12}} \\& {T_{23}} = 2{C_{44}}{\varepsilon _{23}} + 2{C_{45}}{\varepsilon _{13}} \\& {T_{13}} = 2{C_{45}}{\varepsilon _{23}} + 2{C_{55}}{\varepsilon _{13}} \\& {T_{12}} = {C_{16}}{\varepsilon _{11}} + {C_{26}}{\varepsilon _{22}} + {C_{36}}\varepsilon _{33}^{} + 2{C_{66}}{\varepsilon _{12}} \end{split} $$ (2)

      在小变形假设下,应变-位移关系为:

      $$\begin{split}& {\varepsilon _{{\rm{11}}}} = \frac{{\partial {{\rm{u}}_{\rm{1}}}}}{{\partial {x_1}}},\;{\varepsilon _{22}} = \frac{{\partial {u_2}}}{{\partial {x_2}}},\;{\varepsilon _{33}} = \frac{{\partial {u_3}}}{{\partial {x_3}}},\;\\& {\varepsilon _{23}} = \frac{1}{2}\left(\frac{{\partial {u_2}}}{{\partial {x_3}}} + \frac{{\partial {u_3}}}{{\partial {x_2}}}\right),\;\\& {\varepsilon _{13}} = \frac{1}{2}\left(\frac{{\partial {u_1}}}{{\partial {x_3}}} + \frac{{\partial {u_3}}}{{\partial {x_1}}}\right),\;\\& {\varepsilon _{12}} = \frac{1}{2}\left(\frac{{\partial {u_1}}}{{\partial {x_2}}} + \frac{{\partial {u_2}}}{{\partial {x_1}}}\right) \end{split}$$ (3)

      在无重力的影响下,施加初应力的纤维增强复合板中控制波传播的方程为:

      $$\begin{split}& \frac{{\partial {T_{{\rm{11}}}}}}{{\partial {x_1}}} + \frac{{\partial {T_{12}}}}{{\partial {x_2}}} + \frac{{\partial {T_{13}}}}{{\partial {x_3}}} - \frac{1}{2}\left(\frac{{\partial {u_1}}}{{\partial {x_3}}} + \frac{{\partial {u_3}}}{{\partial {x_1}}}\right)\frac{{\partial {S_{xx}}}}{{\partial {x_3}}} +\\& {S_{xx}}\frac{{\partial {\omega _3}}}{{\partial {x_2}}} + ({S_{zz}} - {S_{xx}})\frac{{\partial {\omega _2}}}{{\partial {x_3}}} = \rho \frac{{{\partial ^2}{u_1}}}{{\partial {t^2}}},\\& \frac{{\partial {T_{12}}}}{{\partial {x_1}}} + \frac{{\partial {T_{22}}}}{{\partial {x_2}}} + \frac{{\partial {T_{23}}}}{{\partial {x_3}}} - {S_{zz}}\frac{{\partial {\omega _1}}}{{\partial {x_3}}} + \\& {S_{xx}}\frac{{\partial {\omega _3}}}{{\partial {x_1}}} = \rho \frac{{{\partial ^2}{u_2}}}{{\partial {t^2}}},\\& \frac{{\partial {T_{13}}}}{{\partial {x_1}}} + \frac{{\partial {T_{23}}}}{{\partial {x_2}}} + \frac{{\partial {T_{33}}}}{{\partial {x_3}}} + ({S_{zz}} - {S_{xx}})\frac{{\partial {\omega _2}}}{{\partial {x_1}}} - \\& {S_{zz}}\frac{{\partial {\omega _1}}}{{\partial {x_2}}} = \rho \frac{{{\partial ^2}{u_3}}}{{\partial {t^2}}} \end{split}$$ (4)

      其中,

      $$\begin{split}& {\omega _{\rm{1}}} = \frac{{\rm{1}}}{{\rm{2}}}\left(\frac{{\partial {u_3}}}{{\partial {x_2}}} - \frac{{\partial {u_2}}}{{\partial {x_3}}}\right),\\& {\omega _2} = \frac{1}{2}\left(\frac{{\partial {u_1}}}{{\partial {x_3}}} - \frac{{\partial {u_3}}}{{\partial {x_1}}}\right),\\& {\omega _3} = \frac{1}{2}\left(\frac{{\partial {u_2}}}{{\partial {x_1}}} - \frac{{\partial {u_1}}}{{\partial {x_2}}}\right) \end{split}$$ (5)

      式(1)-式(5)中,TijƐij分别代表应力和应变;ui是位移量;CIJ是介质弹性常数;ρ代表物质密度。这里约定,除非另行说明,本文中所有小写下角标的取值为1,2,3,例如ij;而所有大写下角标的取值为1,2,…,6,例如IJ

      考虑复合材料板的应力自由边界条件(即在x3=0和x3=h0处,T33=T13=T23=0)这里引入矩形窗函π(x3):

      $$\pi ({x_3}) = \left\{ {\begin{array}{*{20}{c}} {0,}&{{\text{其它}}}\\ {1,}&{{\rm{0}} \le {x_{\rm{3}}} \le {h_0}} \end{array}} \right.$$

      因为材料特性参数受到位置的约束,故材料弹性常数CIJ、物质密度ρ、水平方向初应力S11和垂直方向初应力S33均可通过矩形窗函数表示为位置的函数,即:

      $$\begin{split}& {C_{IJ}} = {C_{IJ}}\pi \left( {{x_3}} \right) \\& \rho = \rho \pi \left( {{x_3}} \right) \\& {S_{{\rm{11}}}} = {\rm{ - }}P\pi ({x_3}) \\& {{\rm{S}}_{{\rm{33}}}} = {\rm{ - }}Q\pi ({x_3}) \end{split} $$ (6)

      对于纤维增强复合板,假定沿x2方向无限长,则满足平面应变问题的条件,质点位移的自由谐波解可以表示为:

      $$\begin{split}& {u_1} = U({x_3}) \cdot {e^{i(k{x_1} - \omega t)}} \\& {u_2} = V({x_3}) \cdot {e^{i(k{x_1} - \omega t)}} \\& {u_3} = W({x_3}) \cdot {e^{i(k{x_1} - \omega t)}} \end{split} $$ (7)

      其中,U(x3)、V(x3)、W(x3)代表质点位移在x1x2x3方向上的振动幅度,k是波数,ω是角频率。

      将式(2)、式(3)、式(5)-式(7)代入式(4)中,得到单层纤维增强复合板中Lamb波传播的波动方程为:

      $$\tag{7a} \begin{split}& {\{ [({C_{55}} + 0.5P - 0.5Q)U'' + 0.5P'U' + {{(ik)}^2}{C_{11}}U] + }\\& {[{C_{45}}V'' + {{(ik)}^2}{C_{16}}V] + }\\& {[(ik)({C_{13}} + 0.5Q - 0.5P + {C_{55}})W' + 0.5(ik)P'W]\} \pi ({x_3})}+ \\& { \{ {C_{55}}U' + {C_{45}}V' + (ik){C_{55}}W\} \cdot [\delta ({x_3} - a) - \delta ({x_3} - b)]} = \\& {- \rho {\omega ^2}U\pi ({x_3})} \end{split} $$
      $$\tag{7b}\begin{split}& {\{ [({C_{44}} - 0.5Q)V'' + {{(ik)}^2}({C_{66}} - 0.5P)V] + }\\& {[{C_{45}}U'' + {{(ik)}^2}{C_{16}}U] + [(ik)({C_{36}} + {C_{45}})W']\} \pi ({x_3}) + }\\& {\{ {C_{45}}U' + {C_{44}}V' + (ik){C_{45}}W\} \cdot [\delta ({x_3} - a) - \delta ({x_3} - b)]} = \\& {- \rho {\omega ^2}V\pi ({x_3})} \end{split} $$
      $$ \tag{7c}\begin{split}& {\{ (ik)({C_{55}} + {C_{13}} + 0.5P - 0.5Q)U' + (ik)({C_{45}} + {C_{36}})V'}\\& { + [{C_{33}}W'' + {{(ik)}^2}({C_{55}} + 0.5Q - 0.5P)W]\} \pi ({x_3}) + }\\& {\{ (ik){C_{13}}U + (ik){C_{36}}V + {C_{33}}W'\} \cdot [\delta ({x_3} - a) - \delta ({x_3} - b)]}= \\& { - \rho {\omega ^2}W\pi ({x_3})} \end{split} $$

      此时的Lamb波由准对称模态、准反对称模态和准SH模态三种模式耦合而成。为得到方程(7a)、(7b)和(7c)控制的耦合波动方程的解,把质点在x1x2x3方向上的振动幅度U(x3)、V(x3)和W(x3)用勒让德正交多项式展开[20],即:

      $$\begin{split}& U({x_3}) = \sum\limits_{\alpha = 0}^\infty {p_\alpha ^1{Q_\alpha }({x_3})} \\& V({x_3}) = \sum\limits_{\alpha = 0}^\infty {p_\alpha ^2{Q_\alpha }({x_3})} \\& W({x_3}) = \sum\limits_{\alpha = 0}^\infty {p_\alpha ^3{Q_\alpha }({x_3})} \end{split} $$ (8)

      这里,$p_\alpha ^i$(i=1,2,3)是展开系数。

      $${Q_\alpha }({x_3}) = \sqrt {\frac{{2\alpha + 1}}{h}} {P_\alpha }\left(\frac{{2{x_3}}}{h} - 1\right)$$

      其中,Pα代表第α阶勒让德正交多项式,则Qα形成了一组完备的正交归一的多项式组。就理论上而言,α取从0到无穷大值,但实际计算中α取一有限值αΓ,更高阶的项认为是高阶小量,可以忽略不计。

      在式(7a)、(7b)和(7c)两边同乘以$Q_\beta ^*({x_3})$,然后各式对x3从0到h积分,利用勒让德正交多项式的正交性,可以得到3(αΓ+1)个方程:

      $$\begin{split}& (A_{11}^{\beta ,\alpha } + {\omega ^2}M_\alpha ^\beta )P_\alpha ^1 + A_{12}^{\beta ,\alpha }P_\alpha ^2 + A_{13}^{\beta ,\alpha }P_\alpha ^3 = 0 \\& A_{21}^{\beta ,\alpha }P_\alpha ^1 + (A_{22}^{\beta ,\alpha } + {\omega ^2}M_\alpha ^\beta )P_\alpha ^2 + A_{23}^{\beta ,\alpha }P_\alpha ^3 = 0 \\& A_{31}^{\beta ,\alpha }P_\alpha ^1 + A_{32}^{\beta ,\alpha }P_\alpha ^2 + (A_{33}^{\beta ,\alpha } + {\omega ^2}M_\alpha ^\beta )P_\alpha ^2 = 0 \end{split} $$ (9)

      这里,假定公式中的哑标适用爱因斯坦求和约定。$A_{ij}^{\beta ,\alpha }$$M_\alpha ^\beta $是系数矩阵元素,αβ的取值均为从0到αΓ范围内的正整数。要想使方程(9)有非零解,则方程(9)系数矩阵的行列式必须为零,即:

      $$\left| {\begin{array}{*{20}{c}} {A_{11}^{\beta {\rm{,}}\alpha } + {\omega ^2}M_\alpha ^\beta }&{A_{12}^{\beta {\rm{,}}\alpha }}&{A_{13}^{\beta {\rm{,}}\alpha }} \\ {A_{21}^{\beta {\rm{,}}\alpha }}&{A_{22}^{\beta {\rm{,}}\alpha } + {\omega ^2}M_\alpha ^\beta }&{A_{23}^{\beta {\rm{,}}\alpha }} \\ {A_{31}^{\beta {\rm{,}}\alpha }}&{A_{32}^{\beta {\rm{,}}\alpha }}&{A_{33}^{\beta {\rm{,}}\alpha } + {\omega ^2}M_\alpha ^\beta } \end{array}} \right| = {\rm{0}}$$ (10)

      方程(10)的频散方程可以化归为一个特征值问题,即

      $$\left( {\begin{array}{*{20}{c}} {A_{11}^{\beta ,\alpha }}&{A_{12}^{\beta ,\alpha }}&{A_{13}^{\beta ,\alpha }} \\ {A_{21}^{\beta ,\alpha }}&{A_{22}^{\beta ,\alpha }}&{A_{23}^{\beta ,\alpha }} \\ {A_{31}^{\beta ,\alpha }}&{A_{32}^{\beta ,\alpha }}&{A_{33}^{\beta ,\alpha }} \end{array}} \right)\left( {\begin{array}{*{20}{c}} {p_\alpha ^1} \\ {p_\alpha ^2} \\ {p_\alpha ^3} \end{array}} \right) = {\omega ^2}M_\alpha ^\beta \left( {\begin{array}{*{20}{c}} {p_\alpha ^1} \\ {p_\alpha ^2} \\ {p_\alpha ^3} \end{array}} \right)$$ (11)

      上式即为声弹Lamb波沿纤维增强复合板任意方向传播时的频散方程。随着αΓ值的增加,其特征值亦会相应增加,该问题只取随着αΓ值的增加,数值上趋于收敛的解。

    • 基于上述理论公式,利用Mathematica®软件编写程序对纤维增强复合板中沿任意方向传播的声弹Lamb波的波动方程进行求解。为验证计算结果的正确性,首先将程序应用于无初应力的各向同性材料中,并与Disperse®软件(全局矩阵法)的计算结果进行比较。然后以碳纤维增强复合材料T300/914单层板为例,计算了无初应力各向异性板的波结构应力曲线,验证了应力分布曲线的连续性及边界应力与初始条件的一致性。最后,分析了碳纤维增强复合材料T300/914单层板中初应力效应和纤维铺层方向改变时,对Lamb波频散特性和波结构特性的影响。

    • 为了验证本文方法所编写程序的正确性,首先计算了无初应力作用时,各向同性铝板6061-T6中Lamb波的频散曲线。所采用铝板的厚度为1 mm,其它材料特性参数如下:密度ρ = 2700 kg/m3,弹性模量E=69 GPa,泊松比v = 0.3。图2为本文方法所计算的铝板相速度频散曲线与Disperse®软件计算结果对比图,其中实心圆点为本文计算结果,实线为Disperse®计算结果。从图2可以看出,两种方法计算结果吻合良好。

      图  2  本文方法所计算的铝板相速度频散曲线与Disperse®软件计算结果对比图

      Figure 2.  Comparison of phase velocity dispersion curves for aluminum plate by the authors’ program with available data calculated from Disperse®.

      进一步验证本文方法对单斜晶系材料的适用性,以单层单向碳纤维增强复合材料T300/914为例,图3计算了无初应力情况下纤维铺层方向0°时,Lamb波A0模态在300 kHz频率下的应力波结构曲线,算例中材料的密度${\rho _f}$ = 1560 kg/m3,厚度${h_f}$ = 1 mm,弹性常数如表1所示[15]。由图3可见,应力分布曲线呈现出较好的连续性,在上下表面处应力分量T33T13T23均与应力自由边界条件一致。应力分量以T11为主,说明Lamb波在传播过程中A0模态能流主要沿x1方向传播。

      图  3  300 kHz频率下Lamb波A0模态的应力波结构曲线

      Figure 3.  Stress distributions of Lamb wave A0 mode at 300 kHz.

      表 1  碳纤维增强复合材料T300/914的弹性常数

      Table 1.  Elastic constants of carbon-fiber reinforced composites T300/914.

      单位:cij(GPa)
      C11C12C13C22C23C33C44C55C66
      143.86.26.213.36.513.33.65.75.7
    • 考虑到截止项αΓ对本文算例精确性的影响,图4给出了碳纤维增强复合材料T300/914单层板在无初应力状态时,截止项αΓ由3变化到7时各模态频散曲线的收敛情况。从图4中可以看出,当截止项αΓ增大至5时,各模态的频散曲线漂移已基本趋于稳定。图5S0模态频散曲线收敛过程的放大图。由图5可见,随着截止项αΓ由3到7变化,S0模态频散曲线的收敛过程存在着“跳跃”现象,即当截止项αΓ由3变化到5时,计算结果迅速收敛;而当截止项αΓ由5变化到7时,计算结果的收敛速度则十分缓慢。这主要是因为渐进求解方法计算结果的精度受截止项舍断累计误差的影响,其收敛过程呈现非线性变化的特性,当截止项取值较小时,计算结果与理论真值之间误差较大。为保证本文所讨论频率范围内计算结果的精度,下文算例中均取截止项αΓ = 7进行计算。

      图  4  截止项αΓ由3变化到7时各模态频散曲线的收敛

      Figure 4.  Convergence of dispersion curves of all the modes with the change of cut-off term αΓ from 3 to 7.

      图  5  S0模态频散曲线收敛过程的放大图

      Figure 5.  Zoom-in graph of convergence process of the S0 mode.

    • 为研究水平方向上初应力对单层单向纤维增强复合板中沿任意(含非主对称轴)方向传播的耦合Lamb波的影响,图6计算了沿水平方向施加±1 GPa的初应力时,Lamb波沿0°方向传播的相速度频散曲线。注意,符号“+”表示在水平方向施加压应力,符号“-”表示在水平方向施加拉应力。这里,需要提及地是1 GPa的拉应力粗略地相当于该材料所能承受极限强度的30%。由图6中可以看出,对于低阶模态而言,低频段水平初应力对A0模态有着更为显著的影响,SH0模态其次,对S0模态的影响最小。高频段,高阶模态的截止频率在压应力的影响下呈现增加趋势(向右漂移),在拉应力的影响下呈现减少趋势(向左漂移),压应力与拉应力对高阶模态截止频率的影响效果相反。

      图  6  ±1 GPa初应力作用下,单层单向纤维增强复合材料板中沿0°方向传播的Lamb波相速度频散曲线

      Figure 6.  Phase velocity dispersion curves of Lamb waves propagating along 0° direction in fiber reinforced composite lamina under the initial stresses of ±1 GPa.

      上文分析已表明,A0模态在低频段受水平初应力的影响更为明显。鉴于A0模态的易激发性和工程适用的潜在可能性,下文着重对A0模态的波结构位移曲线在不同程度水平初应力作用时的变化情况进行分析。

      图7图8图9分别为纤维增强复合板在0 GPa、0.8 GPa和1.6 GPa水平初应力状态下,A0模态在激励频率300 kHz时沿15°、45°和75°方向(即局部坐标系x1´轴与全局坐标系x1轴的夹角分别为15°、45°和75°)传播的位移分布曲线。由图7图8图9可以看出,在单层纤维增强复合材料板中,A0模态的质点位移以平面的法向位移分量u3为主,且u3曲线关于板的中心平面呈对称特性,而u1u2曲线则呈现出反对称特性。不同于各向同性板中A0模态传播的情形,纤维增强复合材料板中A0模态的质点位移包含了u1、u2u3三个分量,暗示出此时的A0模态实际上为准Lamb波模态。

      图  7  传播方向为15°时,频率300 kHz的A0模态在不同水平初应力作用下的位移分布

      Figure 7.  Displacement distributions of A0 mode at 300 kHz along 15° propagation direction under different horizontal initial stresses

      图  8  传播方向为45°时,频率300 kHz的A0模态在不同水平初应力作用下的位移分布

      Figure 8.  Displacement distributions of A0 mode at 300 kHz along 45° propagation direction under different horizontal initial stresses

      图  9  传播方向为75°时,频率300 kHz的A0模态在不同水平初应力作用下的位移分布

      Figure 9.  Displacement distributions of A0 mode at 300 kHz along 75° propagation direction under different horizontal initial stresses

      进一步观察图7可以发现,水平方向初应力效应对u2曲线的影响最为明显,u1其次,对u3曲线的影响最弱。表明尽管总体而言,Lamb波A0模态对水平初应力的影响更为敏感,但由于Lamb波沿非主轴方向传播时各模态耦合效应的存在,水平初应力主要通过影响u2位移分量实现对A0模态的调控。如图8所示,当Lamb波A0模态的传播方向变为45°时,水平初应力的影响较15°时更加明显。对比图9图7图8则发现,随着传播方向角度超过45°逐渐增大到75°,位移波结构u3曲线发生了明显的向下弯曲。同时,u2位移曲线在靠近上下表面处的幅度逐渐增大,u1位移曲线靠近上下表面处的幅度逐渐减小,表明此时耦合Lamb波的能量分布随着传播方向的改变而发生变化。考虑水平初应力的影响,当传播方向超过45°角时,水平初应力效应对波结构位移曲线的整体影响趋弱,该结果表明对于实际的无损检测和结构健康检测而言,水平初应力对沿不同方向传播的Lamb波存在不同程度的影响。

    • 本小节继续研究了垂直方向上初应力效应对单层单向纤维增强复合板中沿非主对称轴方向传播的耦合Lamb波的影响。如图10所示,分别计算了纤维增强复合板当在垂直方向受到±1 GPa初应力作用时的相速度频散曲线。由图10的计算结果可见,对于低阶模态而言,低频段垂直初应力依然对A0模态有着最为明显的影响,对SH0模态和S0模态的影响较小。高频范围内对高阶模态的影响与图6的水平初应力效应相反。整体来看垂直初应力对高阶模态的影响较水平初应力更大,表现为高频段频散曲线的漂移程度更大。此外,压应力与拉应力的影响效果也完全相反。

      图  10  ±1 GPa垂直初应力作用下,单层单向纤维增强复合材料板中沿0°方向传播的Lamb波相速度频散曲线

      Figure 10.  Phase velocity dispersion curves of Lamb waves propagating along 0° direction in fiber reinforced composite lamina under the vertical initial stresses of ±1 GPa.

      同样地,针对垂直方向初应力效应,图11图12图13着重对A0模态在300 kHz频率下,沿不同方向传播时在不同程度垂直初应力状态下的波结构位移分布变化情况进行分析。由三幅图的结果可见,当Lamb波A0模态分别沿15°、45°和75°方向传播时,随着垂直初应力由0 GPa向0.8 GPa和1.6 GPa变化,垂直初应力效应对A0模态整体波结构位移分布曲线的影响有限,其变化量较水平初应力的影响要小。相较于垂直初应力的影响,纤维铺层方向的改变更能引起Lamb波A0模态位移波结构曲线的改变,即此时Lamb波A0模态波结构的变化主要由Lamb波传播方向的改变引起,而非垂直初应力的影响。相比较而言,当Lamb波沿45°方向传播时,垂直初应力效应对u1u2的影响最为明显,表明垂直初应力效应依然是通过耦合效应凭借影响u1u2位移分量的大小,实现对A0模态沿板厚方向上的能量调控。

      图  11  传播方向为15°时,频率300 kHz的A0模态在不同垂直初应力作用下的位移分布

      Figure 11.  Displacement distributions of A0 mode at 300 kHz along 15° propagation direction under different vertical initial stresses

      图  12  传播方向为45°时,频率300 kHz的A0模态在不同垂直初应力作用下的位移分布

      Figure 12.  Displacement distributions of A0 mode at 300 kHz along 45° propagation direction under different vertical initial stresses

      图  13  传播方向为75°时,频率300 kHz的A0模态在不同垂直初应力作用下的位移分布

      Figure 13.  Displacement distributions of A0 mode at 300 kHz along 75° propagation direction under different vertical initial stresses

    • 采用勒让德正交多项式展开法,研究了水平和垂直方向上施加初应力时,沿纤维增强复合板非主对称轴方向传播的声弹Lamb波。通过本文的研究发现,利用勒让德正交多项式展开法的正交特性,可以把声弹Lamb波耦合波动方程的求解问题方便地转化为特征值问题,合理的截止项取值使得对求解过程获得了很快的收敛解。对于单层单向纤维增强复合板来说,在水平方向施加压应力与拉应力的影响效果相反,但水平方向的拉应力与垂直方向压应力的影响效果基本相同。对于低阶的Lamb波A0模态而言,初应力效应主要通过耦合效应影响波结构位移分量的大小来实现Lamb波传播能量的调控。相较于纤维铺层方向的改变,即Lamb波传播方向的变化,Lamb波高阶模态的截止频率对初应力效应有着更为敏感的变化。对于纤维增强复合板来讲,初应力效应对沿不同方向传播的Lamb波模态亦有着不同程度的影响。因此,在无损检测和结构健康监测中,应妥善考虑Lamb波接收角度和初应力效应的影响,相关结论亦可进一步用于传感器监测网络的布局优化。

WeChat 关注分享

返回顶部

目录

    /

    返回文章
    返回