筋条形状对复合材料加筋壁板低速冲击动态响应的影响

丁铁 张晓晴 姚小虎

任会兰, 宁建国, 许香照. 不同炸药量在工事中爆炸的三维数值模拟[J]. 高压物理学报, 2013, 27(2): 216-222. doi: 10.11858/gywlxb.2013.02.008
引用本文: 丁铁, 张晓晴, 姚小虎. 筋条形状对复合材料加筋壁板低速冲击动态响应的影响[J]. 高压物理学报, 2017, 31(6): 769-777. doi: 10.11858/gywlxb.2017.06.012
REN Hui-Lan, NING Jian-Guo, XU Xiang-Zhao. The 3-D Numerical Simulation for Different Explosive Charges in the Fortifications[J]. Chinese Journal of High Pressure Physics, 2013, 27(2): 216-222. doi: 10.11858/gywlxb.2013.02.008
Citation: DING Tie, ZHANG Xiao-Qing, YAO Xiao-Hu. Effect of Stiffener Shape on Low-Velocity Impact Behavior of Composite Stiffened Panels[J]. Chinese Journal of High Pressure Physics, 2017, 31(6): 769-777. doi: 10.11858/gywlxb.2017.06.012

筋条形状对复合材料加筋壁板低速冲击动态响应的影响

doi: 10.11858/gywlxb.2017.06.012
基金项目: 

国家自然科学基金 11472110

国家自然科学基金 11372113

国家自然科学基金 11472100

华南理工大学中央高校专项资金 2015ZP042

广东省普通高校特色创新项目 2014ktscx015

爆炸科学与技术国家重点实验室开放基金 KFJJ15-20M

详细信息
    作者简介:

    丁铁(1991—), 男,硕士研究生,主要从事复合材料与冲击动力学研究.E-mail:ctd.tie@mail.scut.edu.cn

    通讯作者:

    姚小虎(1974—), 男,教授,博士生导师,主要从事复合材料与冲击动力学研究.E-mail:yaoxh@scut.edu.cn

  • 中图分类号: O347.3;V214.8

Effect of Stiffener Shape on Low-Velocity Impact Behavior of Composite Stiffened Panels

  • 摘要: 针对复合材料加筋壁板的抗冲击问题,着重考虑筋条形状对其冲击动态响应的影响。通过局部细化的方法分别建立T型、工字型和帽型筋条的复合材料加筋壁板的有限元分析模型,采用黏结接触模拟层与层之间的界面,分析了3种筋条形状的复合材料加筋壁板在低速冲击下的动态响应和损伤情况。根据数值计算结果发现:结构的动态响应与冲击位置的刚度密切相关,T型和工字型筋条复合材料加筋壁板在冲击位置以外的区域发生翘曲;此时T型和工字型筋条的基体拉伸损伤不仅形状相似,而且损伤面积相近;蒙皮与筋条界面相邻两层的铺层角影响分层损伤的面积。

     

  • 随着海洋资源的开发与水下爆破技术的发展, 水下爆破安全防护问题逐渐受到人们的重视。在陆地岩石中爆破, 地震波能量一般不足爆破总能量的10%;而在水下介质中爆破时, 地震波能量可达20%。水下爆破时结构的振动效应较陆地岩石中更复杂, 一般是低频波与高频波共同作用的结果, 其主要特征是振动频带宽、高频丰富、振幅大[1-3]。在浅水域近距离处开展水下钻孔爆破工程时, 通过采用水下预裂爆破、气泡帷幕防护、空气层防护、多段延时起爆等措施削减地震波和水中冲击波[4]。由于空气的波阻抗比水小很多, 当水下冲击波传播到空气与水的交界面时, 水下冲击波将在介质交界面处发生反射, 其透射系数较小, 少部分能量透射到空气介质中, 可有效缓冲水下冲击波和衰减地震能量[5-6]。樊自建等[7]、贾虎等[8-9]、寇晓枫等[10]开展实验研究了空气隔层对水中冲击波的衰减效果, 得出结论:空气隔层能有效衰减冲击波峰值压力、冲击波冲量及总能量, 对保护对象有较好的防护效果。姚熊亮等[11-12]、刘欣等[13]建立了多介质水下爆炸数值模型进行模拟计算, 其研究结果表明, 空气层结构不仅可以防护接触爆炸产生的冲击波, 而且对非接触爆炸也有很好的防护效果。

    本课题组前期针对不同段别雷管的水下爆炸振动信号进行了Hilbert-Huang变换(Hilbert-Huang Transform, HHT)分析[14], 但未进行有机玻璃-空气层结构对水下爆炸振动影响的研究, 且国内外相关文献在此方面报道甚少。基于此考虑, 本工作在小水池中进行有机玻璃-空气层结构约束下的水下爆炸振动测试, 研究空气层厚度对水下爆炸最大振动速度的影响, 通过HHT分析其地基振动信号, 探究空气层厚度对水下爆炸振动全局频率的影响, 为爆炸水池减振、隔振以及减弱爆炸振动对建筑物的破坏作用提供理论支撑。

    HHT法适用于处理非线性、非平稳信号, 主要由经验模态分解(Empirical Mode Decomposition, EMD)和Hilbert变换两部分组成[15-17]。相对于小波分析, EMD分解克服了小波基选择问题, 避免了由于基函数不同所带来的影响。另外, 原始信号经EMD分解后得到的固有模态函数(Intrinsic Mode Function, IMF)具有不同特征尺度, 该过程是自适应的、高效的。通过EMD分解得到的各IMF分量经Hilbert变换后, 得到各IMF分量的瞬时频谱, 综合所有IMF分量的瞬时频谱得到Hilbert谱。

    爆炸振动原始信号中包含3个方向的信号, 依据《爆破安全规程》(GB 6722-2014)中的13.2.2条规定, 选取振动速度值最大的垂向振动信号(记为X(t))进行分析。基于HHT法对爆炸振动信号进行EMD分解以及Hilbert变换, 过程如下。

    原始信号X(t)经n次EMD分解, 可表示为

    X(t)=ni=1Ci(t)+rn(t)
    (1)

    式中: Ci(t)为经EMD分解得到的第i阶IMF分量, rn(t)为原始信号减去n个IMF分量之和的残差。

    对得到的IMF分量做Hilbert变换

    H[Ci(t)]=1πPVCi(t)ttdt
    (2)

    式中:PV表示柯西主值。

    构造解析信号Z(t)

    Z(t)=Ci(t)+jH[Ci(t)]=ai(t)ejφi(t)
    (3)

    式中:ai(t)为幅值函数, φi(t)为相位函数。

    ai(t)=C2i(t)+H2[Ci(t)]
    (4)
    φi(t)=tan1H[Ci(t)]Ci(t)
    (5)

    原始信号分解后可表示为

    X(t)=Reni=1ai(t)ejφi(t)=Reni=1ai(t)eωi(t)dt
    (6)

    Hilbert谱可表示为

    H(ω,t)=Reni=1ai(t)eωi(t)dt
    (7)

    式中:ω表示角频率。

    根据以上理论, 采用Matlab编写程序代码, 加载相应模块即可作出相关Hilbert谱。

    设计了4种类型的有机玻璃-空气层结构, 空气层厚度(直径)分别为40、80、120和160 mm。中心装药管为壁厚1 mm的有机玻璃管, 外层为壁厚2 mm的有机玻璃壳体, 装药中心与复合结构中心在同一高度, 如图 1所示。

    图  1  有机玻璃-空气层结构
    Figure  1.  Picture of plexiglass-air interlayer structure

    采用NUBOX-6016型爆破振动测试仪记录振动信号, 其传感器使用TP-3V-10144型号, 采用单次采样方式, 采样频率和时间分别为5 kHz、2 s, 延时点数-256, 触发方式为通道出窗触发, 触发条件为0.05 cm/s, 耦合方式为直流耦合。

    在直径为5.5 m、高3.62 m的小型爆炸水池中进行水下爆炸试验。入水深度的选取是基于消除爆炸水池边界效应的考虑。一般装药置于水池中心总深度的2/3处时, 来自水面和池底的边界效应可互相抵消[18], 因此装药中心的入水深度选取为2.4 m。选用3台测振仪, 分别布设在距爆炸水池筒体外壁1、2、3 m处(地基振动速度在3 m处已衰减到较小数值, 而1、2、3 m处测点为等间距的典型观测点, 分别标记为1#、2#、3#)。试验系统如图 2所示。

    图  2  振动测试系统(1.爆炸水池; 2.固定架; 3.8号标准雷管; 4.空气隔层; 5.发爆器; 6.1#、2#、3#测振仪)
    Figure  2.  Vibration test system (1.Explosion pool; 2.Fixed mount; 3.Detonator; 4.Air interlayer; 5.Initiator; 6.1#, 2#, 3# vibration measurer)

    根据《爆破安全规程》(GB 6722-2014)规定, 选取xyz 3个方向中振动速度最大的z方向(与地面垂直方向)振动信号进行分析。测试得到不同距离、不同厚度空气层条件下的水下爆炸振动参数, 如表 1所示。为了探究有机玻璃壳体对水下爆炸振动速度的影响, 进行填充水和空气层的对比试验。表 2为不同距离、不同厚度水隔层情况下的水下爆炸振动参数。

    表  1  不同距离、不同厚度空气层对水下爆炸振动的影响
    Table  1.  Influence of distances and thickness of air interlayer on underwater explosion vibration
    Distance/m Direction Main vibration frequency/Hz Maximum vibration velocity/(cm·s-1)
    ∅40 mm ∅80 mm ∅120 mm ∅160 mm ∅40 mm ∅80 mm ∅120 mm ∅160 mm
    1 z 56.763 56.763 56.763 56.763 0.086 0.070 0.054 0.070
    2 z 56.763 56.763 56.763 44.556 0.118 0.097 0.082 0.101
    3 z 56.763 56.763 56.763 56.763 0.080 0.071 0.061 0.086
    下载: 导出CSV 
    | 显示表格
    表  2  不同距离、不同厚度水隔层对水下爆炸振动的影响
    Table  2.  Influence of distances and thickness of water interlayer on underwater explosion vibration
    Distance/m Direction Main vibration frequency/Hz Maximum vibration velocity/(cm·s-1)
    ∅40 mm ∅80 mm ∅120 mm ∅160 mm ∅40 mm ∅80 mm ∅120 mm ∅160 mm
    1 z 62.256 62.256 62.256 56.763 0.086 0.169 0.184 0.185
    2 z 62.256 62.256 62.256 62.256 0.155 0.145 0.160 0.161
    3 z 62.256 62.256 62.256 62.256 0.101 0.097 0.107 0.116
    下载: 导出CSV 
    | 显示表格

    表 1表 2可看出:相对于水隔层, 不同厚度空气层条件下各测点处主振频率皆有所降低; 部分情况下出现主振频率值不变的现象, 可能由于爆炸当量较小及测振仪灵敏度较小等原因, 因此需要进一步对信号进行频谱分析。图 3图 4分别为不同厚度空气层、水隔层结构对各测点振动速度的影响。

    图  3  不同厚度空气层对各测点振速的影响
    Figure  3.  Influence of different thicknesses of air interlayer on vibration velocity
    图  4  不同厚度水隔层对各测点振速的影响
    Figure  4.  Influence of different thicknesses of water interlayer on vibration velocity

    图 3图 4可看出, 对比相同厚度的空气层和水隔层结构, 空气层结构对水下爆炸振动的衰减更明显。相同厚度水隔层结构条件下:3 m处的振动速度衰减程度最大, 振动速度值最小; 2 m处速度值次之; 1 m处速度值最大。而相同厚度空气层结构条件下:1 m处的振动速度衰减程度最大, 振动速度值最小; 3 m处速度值次之; 2 m处速度值最大。可见空气层结构对近区1 m处的振动速度影响较大。

    目前国内使用的振动信号采集仪主要采集最大振幅所对应的频率, 即主振频率, 无法对振动信号整个过程的振动频率进行监测; 但爆破工程中需要保护的对象有其自身的固有频率, 为了避免共振造成的破坏, 需要测试振动信号的全局频率, 以观测危险频段的振幅值是否需要进一步防护。本试验爆炸当量较小, 振动信号相对较弱, 同时可能由于测振仪灵敏度较小, 出现部分情况主振频率值不变的现象。因此, 有必要对振动信号进行Hilbert谱分析。

    通过Matlab软件编写相应程序, 将振动测试仪采集到的振动信号读取到矩阵中, 并根据试验设置的采集频率、采集时长等设置运行参数。对事先装载的EMD工具箱, 只需调用命令函数即可使用EMD分解算法。以空气层厚度为40 mm、距离1 m处的垂直方向振动信号为例(见图 5), 利用EMD分解该爆炸振动信号得到各IMF分量, 如图 6所示。

    图  5  原始振动信号
    Figure  5.  Original vibration signal
    图  6  EMD分解后得到各IMF分量
    Figure  6.  IMF components by EMD decomposition

    图 6可看出, 采用EMD分解得到12个IMF分量C1~C12和一个余量r12。其中:C1C2分量频率最高, 所含能量较低, 代表振动测试过程中周围环境的高频噪声; C3分量代表信号的高频成分; C4~C9分量幅值最大, 所含的能量最多, 是信号的优势频率段, 包含信号的最显著特征; C10~C12为信号的低频成分; r12幅值很小, 可能是信号微弱的趋势或仪器的飘零。

    对水下爆炸振动信号经EMD分解后得到的各IMF分量与原始信号进行相关分析, 得到相关分析参数, 如表 3所示。从表 3可看出:各IMF分量的方差值几乎全为零, 且包含信号主要信息的IMF分量段的贡献率达到96.77%。可见, EMD分解过程具有准确性、高效性和自适应性, 保证了信号经分解后仍保留水下爆炸振动信号的非平稳特性, 表明HHT分析方法适用于水下爆炸振动信号处理。

    表  3  信号经EMD分解后相关参数
    Table  3.  Related parameters of the signal by EMD decomposition
    IMF Variance Proportion/%
    C1 0 1.646 3
    C2 0 0.862 9
    C3 0 0.626 8
    C4 0.000 1 70.307 9
    C5 0 20.200 7
    C6 0 2.208 9
    C7 0 1.254 8
    C8 0 1.081 4
    C9 0 1.719 6
    C10 0 0.015 0
    C11 0 0.060 6
    C12 0 0.015 0
    下载: 导出CSV 
    | 显示表格

    对垂向爆炸振动信号经EMD分解得到的IMF分量进行Hilbert变换, 得到频率-时间-振幅三维Hilbert谱。以距离1 m、水隔层厚80 mm的情况为例, 对其振动信号的IMF分量进行Hilbert变换得到Hilbert谱, 如图 7所示。

    图  7  1 m处充80 mm水隔层垂向振动Hilbert谱
    Figure  7.  Vertical vibration Hilbert spectrum of the 80 mm water-filled interlayer at 1 m

    图 7可看出:幅值较大(即颜色较深)的样点主要分布于7 000时间采样点(即1.4 s)之前, 且800时间采样点(即0.16 s)前振动幅值最大, 相对应的频率为60~65 Hz, 与测试仪采集到的主振频率接近; 频率为5~15 Hz处对应的幅值较大, 且作用时间较长, 该频率段与房屋等建筑物的固有频率较为接近, 因此可能对建筑物造成更大的破坏, 在实际工程爆破中需要加以防护。为此, 采用不同厚度的空气层结构进行防护, 并对采集到的振动信号进行Hilbert谱分析, 如图 8所示。

    图  8  不同空气层厚度条件下的垂向振动Hilbert谱
    Figure  8.  Hilbert spectrum of vertical vibration of different thicknesses of air interlayer

    图 8可看出, 在有机玻璃-空气层结构下, 最大振幅所对应的主振频率为45~55 Hz, 与测振仪采集到的主振频率接近。相比于水隔层, 有机玻璃-空气层结构下5~15 Hz区间段的幅值减弱较明显, 且作用时间得到缩短, 但随着空气层厚度的增加, 该频带所对应的幅值有所增加。可见, 针对本试验条件, 为了减小水下爆炸引起的地基振动, 并非空气层越厚越好, 当其厚度为120 mm时防护效果较好。

    (1) 对于有机玻璃-空气层结构, 随着空气层厚度的增加, 水下爆炸最大振动速度呈先减小后增大趋势, 且当空气层厚度为120 mm时, 振动速度值最小, 对地基振动的衰减效果较好。

    (2) HHT法能有效解决振动测试仪对小药量振动频率监测不敏感等问题, 同时可以得到全局频率所对应的幅值, 直观得出所需要防护的振动频率。

    (3) 通过HHT法分析得出, 有机玻璃-空气层结构作用下, 能有效衰减5~15 Hz低频区间段的振动速度, 且能够缩短作用时间, 有效防止与建/构筑物之间产生共振。

  • 图  层内损伤模型

    Figure  1.  Intralaminar damage model

    图  数值模拟所得冲击力时程曲线与实验数据的比较

    Figure  2.  comparison of impact force time-history between numerical and experimental results

    图  复合材料加筋壁板尺寸示意(单位:mm)

    Figure  3.  Configuration of composite stiffened panel (Unit:mm)

    图  复合材料加筋壁板模型细节

    Figure  4.  Details of composite stiffened panel model

    图  复合材料加筋壁板的冲击力时程

    Figure  5.  Impact force versus time of composite stiffened panels

    图  复合材料加筋壁板的位移时程

    Figure  6.  Displacement versus time of composite stiffened panels

    图  复合材料加筋壁板沿路径DI上点的位移

    Figure  7.  Displacements along path DI of composite stiffened panel

    图  工字型筋条复合材料加筋壁板的有限元模型

    Figure  8.  Finite element model of I-shaped stiffened panel

    图  复合材料加筋壁板蒙皮的基体拉伸损伤(1、15、17、19、20为铺层序号)

    Figure  9.  Skin's matrix tensile damage of composite stiffened panel (1, 15, 17, 19, 20 are ply numbers)

    图  10  复合材料加筋壁板筋条的基体拉伸损伤(1、2、3、5、9为铺层序号)

    Figure  10.  Stiffener's matrix tensile damage of composite stiffened panel (1, 2, 3, 5, 9 are ply numbers)

    图  11  复合材料加筋壁板蒙皮的基体拉伸损伤分布

    Figure  11.  Distribution of skin's matrix tensile damage of composite stiffened panels

    图  12  蒙皮与筋条间的分层损伤

    Figure  12.  Delamination damage between skin and stiffener

    图  13  不同复合材料加筋壁板筋条铺层之间的分层损伤对比

    Figure  13.  Comparison of delamination damage between stiffener's plies of different composite stiffened panels

    表  1  复合材料单层材料参数[12]

    Table  1.   Material properties of composite[12]

    Elastic modulus Shear modulus Poisson's ratio
    E11/(GPa) E22/(GPa) E33/(GPa) G12/(MPa) G13/(MPa) G23/(MPa) ν12 ν13 ν23
    145 10.3 12.1 5 300 5 275 3 950 0.301 0.5 0.495
    Strength Fracture energy
    XT/(MPa) XC/(MPa) YT/(MPa) YC/(MPa) S12/(MPa) G1Ct/(N/mm) G1Cc/(N/mm) G2Ct/(N/mm) G2Cc/(N/mm)
    2 000 1 600 64 290 98 12.5 12.5 1 1
    下载: 导出CSV
  • [1] ROESELER W G, SARH B, KISMARTON M U, et al.Composite structures-the first 100 years[C]//16th International Conference on Composite Materials, 2007.
    [2] CHOI H Y, CHANG F K.A model for predicting damage in graphite/epoxy laminated composites resulting from low-velocity impact[J].J Compos Mater, 1992, 26(14):2134-2169. doi: 10.1177/002199839202601408
    [3] 朱炜垚, 许希武.复合材料层合板低速冲击损伤的有限元模拟[J].复合材料学报, 2010, 27(6): 200-207. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK201002862699

    ZHU W Y, XU X W.Finite element simulation of low velocity impact damage on composite laminates[J].Acta Materiae Compositae Sinica, 2010, 27(6): 200-207. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK201002862699
    [4] 王一飞, 张晓晶, 汪海.复合材料层压板低速冲击响应与损伤参数关系研究[J].固体力学学报, 2013, 34(1):63-72. doi: 10.3969/j.issn.0254-7805.2013.01.009

    WANG Y F, ZHANG X J, WANG H.Energy and contact force study of composite laminate impact damage response[J].Chinese Journal of Solid Mechanics, 2013, 34(1):63-72. doi: 10.3969/j.issn.0254-7805.2013.01.009
    [5] 赵维涛, 杨其蛟, 刘炜华.复合材料加筋板低速冲击响应分析[J].沈阳航空航天大学学报, 2015, 32(3): 1-5. doi: 10.3969/j.issn.2095-1248.2015.03.001

    ZHAO W T, YANG Q J, LIU W H.Analysis on the low-velocity impact response of stiffened composite panels[J].Journal of Shenyang Aerospace University, 2015, 32(3): 1-5. doi: 10.3969/j.issn.2095-1248.2015.03.001
    [6] FOUND M S, HOWARD I C, PARAN A P.Impact perforation of thin stiffened CFRP panels[J].Compos Struct, 2000, 48(1/2/3):95-98. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=CC026031071
    [7] FOUND M S, HOWARD I C, PARAN A P.Impact behaviour of stiffened CFRP sections[J].Compos Struct, 1997, 39(3/4):229-235. http://www.sciencedirect.com/science/article/pii/S0263822397001177
    [8] WIGGENRAAD J F M, ZHANG X, DAVIES G A O.Impact damage prediction and failure analysis of heavily loaded, blade-stiffened composite wing panels[J].Compos Struct, 1999, 45(2):81-103. doi: 10.1016/S0263-8223(98)00132-9
    [9] HASHIN Z, ROTEM A.A fatigue failure criterion for fiber reinforced materials[J].J Compos Mater, 1973, 7(4):448-464. doi: 10.1177/002199837300700404
    [10] BENZEGGAGH M L, KENANE M.Measurement of mixed-mode delamination fracture toughness of unidirectional glass/epoxy composites with mixed-mode bending apparatus[J].Compos Sci Technol, 1996, 56(4):439-449. doi: 10.1016/0266-3538(96)00005-X
    [11] GREENHALGH E, MEEKS C, CLARKE A, et al.The effect of defects on the performance of post-buckled CFRP stringer-stiffened panels[J].Compos Part A:Appl Sci Manuf, 2003, 34(7):623-633. doi: 10.1016/S1359-835X(03)00098-8
    [12] 孙晶晶, 张晓晶, 宫占峰, 等.复合材料帽型筋条脱粘的失效机理分析[J].航空学报, 2013, 34(7):1616-1626. http://d.old.wanfangdata.com.cn/Periodical/hkxb201307013

    SUN J J, ZHANG X J, GONG Z F, et al.Failure mechanism study on omega stringer debonding[J].Acta Aeronautica et Astronautica Sinica, 2013, 34(7):1616-1626. http://d.old.wanfangdata.com.cn/Periodical/hkxb201307013
  • 加载中
图(13) / 表(1)
计量
  • 文章访问数:  7411
  • HTML全文浏览量:  3009
  • PDF下载量:  170
出版历程
  • 收稿日期:  2017-01-11
  • 修回日期:  2017-03-29

目录

/

返回文章
返回