Effect of Stiffener Shape on Low-Velocity Impact Behavior of Composite Stiffened Panels
-
摘要: 针对复合材料加筋壁板的抗冲击问题,着重考虑筋条形状对其冲击动态响应的影响。通过局部细化的方法分别建立T型、工字型和帽型筋条的复合材料加筋壁板的有限元分析模型,采用黏结接触模拟层与层之间的界面,分析了3种筋条形状的复合材料加筋壁板在低速冲击下的动态响应和损伤情况。根据数值计算结果发现:结构的动态响应与冲击位置的刚度密切相关,T型和工字型筋条复合材料加筋壁板在冲击位置以外的区域发生翘曲;此时T型和工字型筋条的基体拉伸损伤不仅形状相似,而且损伤面积相近;蒙皮与筋条界面相邻两层的铺层角影响分层损伤的面积。Abstract: In this research, we focused on the effect of the stiffener shape on the dynamic response of composite stiffened panels.We first established the finite element analysis model for composited stiffened structures with 3 kinds of stiffeners using the local refinement with the interface between the layers simulated by bonding interface, and then analyzed the dynamic response and the delamination damage of these composite stiffened panels.According to the results of the numerical simulation, it is found that the dynamic response is closely related to the structure stiffness, the composite stiffened panels with T-shaped and I-shaped stiffeners are warped in the region beyond the impacted part of the panels, and their damage distribution and damage area are similiar.The results also show that the layer angle of the two layers adjacent to the surface of the skin and the stiffener affects the area of the delamination damage.
-
Key words:
- composite stiffened panel /
- low-velocity impact /
- delamination damage /
- dynamic response
-
随着海洋资源的开发与水下爆破技术的发展, 水下爆破安全防护问题逐渐受到人们的重视。在陆地岩石中爆破, 地震波能量一般不足爆破总能量的10%;而在水下介质中爆破时, 地震波能量可达20%。水下爆破时结构的振动效应较陆地岩石中更复杂, 一般是低频波与高频波共同作用的结果, 其主要特征是振动频带宽、高频丰富、振幅大[1-3]。在浅水域近距离处开展水下钻孔爆破工程时, 通过采用水下预裂爆破、气泡帷幕防护、空气层防护、多段延时起爆等措施削减地震波和水中冲击波[4]。由于空气的波阻抗比水小很多, 当水下冲击波传播到空气与水的交界面时, 水下冲击波将在介质交界面处发生反射, 其透射系数较小, 少部分能量透射到空气介质中, 可有效缓冲水下冲击波和衰减地震能量[5-6]。樊自建等[7]、贾虎等[8-9]、寇晓枫等[10]开展实验研究了空气隔层对水中冲击波的衰减效果, 得出结论:空气隔层能有效衰减冲击波峰值压力、冲击波冲量及总能量, 对保护对象有较好的防护效果。姚熊亮等[11-12]、刘欣等[13]建立了多介质水下爆炸数值模型进行模拟计算, 其研究结果表明, 空气层结构不仅可以防护接触爆炸产生的冲击波, 而且对非接触爆炸也有很好的防护效果。
本课题组前期针对不同段别雷管的水下爆炸振动信号进行了Hilbert-Huang变换(Hilbert-Huang Transform, HHT)分析[14], 但未进行有机玻璃-空气层结构对水下爆炸振动影响的研究, 且国内外相关文献在此方面报道甚少。基于此考虑, 本工作在小水池中进行有机玻璃-空气层结构约束下的水下爆炸振动测试, 研究空气层厚度对水下爆炸最大振动速度的影响, 通过HHT分析其地基振动信号, 探究空气层厚度对水下爆炸振动全局频率的影响, 为爆炸水池减振、隔振以及减弱爆炸振动对建筑物的破坏作用提供理论支撑。
1. 振动信号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)=n∑i=1Ci(t)+rn(t) (1) 式中: Ci(t)为经EMD分解得到的第i阶IMF分量, rn(t)为原始信号减去n个IMF分量之和的残差。
对得到的IMF分量做Hilbert变换
H[Ci(t)]=1πPV∫∞−∞Ci(t′)t−t′dt′ (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)=tan−1H[Ci(t)]Ci(t) (5) 原始信号分解后可表示为
X(t)=Ren∑i=1ai(t)ejφi(t)=Ren∑i=1ai(t)e∫ωi(t)dt (6) Hilbert谱可表示为
H(ω,t)=Ren∑i=1ai(t)e∫ωi(t)dt (7) 式中:ω表示角频率。
根据以上理论, 采用Matlab编写程序代码, 加载相应模块即可作出相关Hilbert谱。
2. 有机玻璃-空气层结构条件下小水池水下爆炸振动试验
2.1 试验装置与仪器
设计了4种类型的有机玻璃-空气层结构, 空气层厚度(直径)分别为40、80、120和160 mm。中心装药管为壁厚1 mm的有机玻璃管, 外层为壁厚2 mm的有机玻璃壳体, 装药中心与复合结构中心在同一高度, 如图 1所示。
采用NUBOX-6016型爆破振动测试仪记录振动信号, 其传感器使用TP-3V-10144型号, 采用单次采样方式, 采样频率和时间分别为5 kHz、2 s, 延时点数-256, 触发方式为通道出窗触发, 触发条件为0.05 cm/s, 耦合方式为直流耦合。
2.2 试验系统
在直径为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.3 试验结果与分析
根据《爆破安全规程》(GB 6722-2014)规定, 选取x、y、z 3个方向中振动速度最大的z方向(与地面垂直方向)振动信号进行分析。测试得到不同距离、不同厚度空气层条件下的水下爆炸振动参数, 如表 1所示。为了探究有机玻璃壳体对水下爆炸振动速度的影响, 进行填充水和空气层的对比试验。表 2为不同距离、不同厚度水隔层情况下的水下爆炸振动参数。
表 1 不同距离、不同厚度空气层对水下爆炸振动的影响Table 1. Influence of distances and thickness of air interlayer on underwater explosion vibrationDistance/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 表 2 不同距离、不同厚度水隔层对水下爆炸振动的影响Table 2. Influence of distances and thickness of water interlayer on underwater explosion vibrationDistance/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 由表 1和表 2可看出:相对于水隔层, 不同厚度空气层条件下各测点处主振频率皆有所降低; 部分情况下出现主振频率值不变的现象, 可能由于爆炸当量较小及测振仪灵敏度较小等原因, 因此需要进一步对信号进行频谱分析。图 3、图 4分别为不同厚度空气层、水隔层结构对各测点振动速度的影响。
从图 3、图 4可看出, 对比相同厚度的空气层和水隔层结构, 空气层结构对水下爆炸振动的衰减更明显。相同厚度水隔层结构条件下:3 m处的振动速度衰减程度最大, 振动速度值最小; 2 m处速度值次之; 1 m处速度值最大。而相同厚度空气层结构条件下:1 m处的振动速度衰减程度最大, 振动速度值最小; 3 m处速度值次之; 2 m处速度值最大。可见空气层结构对近区1 m处的振动速度影响较大。
3. 振动信号的HHT分析
目前国内使用的振动信号采集仪主要采集最大振幅所对应的频率, 即主振频率, 无法对振动信号整个过程的振动频率进行监测; 但爆破工程中需要保护的对象有其自身的固有频率, 为了避免共振造成的破坏, 需要测试振动信号的全局频率, 以观测危险频段的振幅值是否需要进一步防护。本试验爆炸当量较小, 振动信号相对较弱, 同时可能由于测振仪灵敏度较小, 出现部分情况主振频率值不变的现象。因此, 有必要对振动信号进行Hilbert谱分析。
3.1 EMD分解
通过Matlab软件编写相应程序, 将振动测试仪采集到的振动信号读取到矩阵中, 并根据试验设置的采集频率、采集时长等设置运行参数。对事先装载的EMD工具箱, 只需调用命令函数即可使用EMD分解算法。以空气层厚度为40 mm、距离1 m处的垂直方向振动信号为例(见图 5), 利用EMD分解该爆炸振动信号得到各IMF分量, 如图 6所示。
从图 6可看出, 采用EMD分解得到12个IMF分量C1~C12和一个余量r12。其中:C1和C2分量频率最高, 所含能量较低, 代表振动测试过程中周围环境的高频噪声; 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 decompositionIMF 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 3.2 振动信号Hilbert谱分析
对垂向爆炸振动信号经EMD分解得到的IMF分量进行Hilbert变换, 得到频率-时间-振幅三维Hilbert谱。以距离1 m、水隔层厚80 mm的情况为例, 对其振动信号的IMF分量进行Hilbert变换得到Hilbert谱, 如图 7所示。
从图 7可看出:幅值较大(即颜色较深)的样点主要分布于7 000时间采样点(即1.4 s)之前, 且800时间采样点(即0.16 s)前振动幅值最大, 相对应的频率为60~65 Hz, 与测试仪采集到的主振频率接近; 频率为5~15 Hz处对应的幅值较大, 且作用时间较长, 该频率段与房屋等建筑物的固有频率较为接近, 因此可能对建筑物造成更大的破坏, 在实际工程爆破中需要加以防护。为此, 采用不同厚度的空气层结构进行防护, 并对采集到的振动信号进行Hilbert谱分析, 如图 8所示。
从图 8可看出, 在有机玻璃-空气层结构下, 最大振幅所对应的主振频率为45~55 Hz, 与测振仪采集到的主振频率接近。相比于水隔层, 有机玻璃-空气层结构下5~15 Hz区间段的幅值减弱较明显, 且作用时间得到缩短, 但随着空气层厚度的增加, 该频带所对应的幅值有所增加。可见, 针对本试验条件, 为了减小水下爆炸引起的地基振动, 并非空气层越厚越好, 当其厚度为120 mm时防护效果较好。
4. 结论
(1) 对于有机玻璃-空气层结构, 随着空气层厚度的增加, 水下爆炸最大振动速度呈先减小后增大趋势, 且当空气层厚度为120 mm时, 振动速度值最小, 对地基振动的衰减效果较好。
(2) HHT法能有效解决振动测试仪对小药量振动频率监测不敏感等问题, 同时可以得到全局频率所对应的幅值, 直观得出所需要防护的振动频率。
(3) 通过HHT法分析得出, 有机玻璃-空气层结构作用下, 能有效衰减5~15 Hz低频区间段的振动速度, 且能够缩短作用时间, 有效防止与建/构筑物之间产生共振。
-
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 -
[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=QK201002862699ZHU 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.009WANG 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.001ZHAO 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/hkxb201307013SUN 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 -