Mechanical Behavior Analysis of Composite Shell with Fracture Defect
-
摘要: 利用湿法缠绕制备了玻璃纤维增强环氧树脂基复合材料壳体,通过拉伸、双悬臂梁和三点端部开口弯曲试验对该复合材料层合板进行基本力学性能评估,得到强度和刚度参数并用于有限元模拟中。同时,在Abaqus中建立了含有不同深度断裂缺陷复合材料壳体的三维渐进损伤有限元模型,预测内压作用下壳体的力学响应。试验和模拟结果表明:复合材料主向拉伸强度为(222.7 ± 18)MPa,弹性模量为39.39 GPa,Ⅰ型和Ⅱ型断裂韧性层间强度分别为(4.67 ± 0.24)和(4.98 ± 0.26)kJ/m2。随着内压增大,Mises应力也不断增大;当断裂缺陷在最深层(接近内压的第一层,深度为18 mm)时,Mises应力最大;当内压为0.3 MPa时,Mises应力高达28.8 MPa,且周向应变小于纵向应变。Abstract: The glass fiber reinforced epoxy resin matrix composite shell was prepared by wet winding, the basic mechanical properties of the composite laminates were evaluated by tensile tests, double cantilever beams and three point end opening bending tests. And then the obtained strength and stiffness parameters were used in the finite element simulation. Meanwhile, a 3D progressive damage finite element model of composite shell containing fracture defects of different depths was established in Abaqus to predict the mechanical response of the shell under internal pressure. The results show that the tensile strength of the composite is (222.7 ± 18) MPa and the main elastic modulus is 39.39 GPa. The fracture toughness of typeⅠand type Ⅱ of interlayer strength are (4.67 ± 0.24)kJ/m2 and (4.98 ± 0.26)kJ/m2, respectively. As the internal pressure increases, Mises stress is increasing. The Mises stress is the largest when the fracture defect is in the deepest layer (the first layer close to the internal pressure, the depth is 18 mm). And when the internal pressure is 0.3 MPa, maximum Mises stress is up to 28.8 MPa, and the circumferential strain is less than the longitudinal strain.
-
Key words:
- composite material /
- shell /
- fracture defect /
- stress /
- strain
-
复合材料是一种性能优良、可设计性强的新型材料,具有比强度、比模量高等优点,被广泛应用于航空航天、汽车制造、风力发电叶片以及压力容器等领域[1-3]。其中玻璃纤维增强环氧树脂基壳体的应用率较高,也是近年来的研究重点[4-5]。用于管道和压力容器等的复合材料壳体因承受较大的力学载荷,壳体内部会发生基体开裂、纤维与基体脱粘、层合板分层和纤维基体同时断裂等损伤,损伤部位产生应力集中,当应力大于材料的许用应力后,会对整个系统的安全构成潜在的威胁[6-8]。因此,通过研究应力、应变等重要参数从而理解复合材料结构壳体的力学响应行为,对提高含有断裂缺陷的复合材料的安全性能具有重要意义[9-11]。
目前,科研人员对复合材料壳体的损伤已进行了大量的有限元模拟研究,通过建立渐进损伤有限元模型预测了复合材料壳体的损伤失效过程。Leh等[12]建立了连续和渐进损伤两种有限元模型,对纤维缠绕压力容器的爆破压力进行预测,与试验结果相比,渐进损伤模型能够更好地模拟复合材料储氢容器的爆破过程。Wang等[13]将材料递减规则和基于表面行为单元结合建立渐进损伤模型,预测了铝内胆碳纤维缠绕复合材料储氢气瓶的极限承载能力和复杂的失效行为。Liu等[14]对3种不同尺寸的碳纤维缠绕压力容器进行失效分析,采用Hashin失效准则并基于能量损伤演变规律,得到与实验结果相符的失效性质参数和爆破强度。杨斌等[15]为了分析混杂复合材料层合板在冲击载荷下的损伤演变过程,利用有限元分析软件Abaqus建立了预测复合材料层合板在低速冲击作用下损伤的三维有限元模型。
还有很多学者采用应变片或光纤埋入复合材料结构以实时监测的实验方法,研究了复合材料层合板和壳体在外压作用下的应变响应过程。Okabe等[16]采用光纤布拉格光栅传感器识别碳纤维增强塑料层合板的裂纹位置,将光纤光栅传感器嵌入层合板交叉层中,并在拉伸载荷作用下测量反射光谱,从测量的光谱判断裂纹位置,证明了小直径光纤光栅传感器也可用于裂纹位置的识别。肖飚等[17]在玻璃纤维缠绕金属内胆复合材料压力容器的制备过程中,将应变传感器埋在金属内胆与玻璃纤维/环氧树脂复合材料层之间,得到了具有原位监测功能的纤维缠绕压力容器。Kanerva等[18]原位测试了碳纤维增强复合材料的层间变化,结果表明埋入应变片可以精准监测复合材料的层间状态。应变参数可以有效反映复合材料层合板及壳体的受力状态,而将有限元模拟和应力和应变响应结合考虑复合材料壳体的损伤行为从而判断整个压力容器健康情况的相关研究还很少。开展含有断裂缺陷的复合材料壳体的力学行为分析,能够为复合材料壳体检测及容器的安全检测提供一定的数据基础。
本研究将在玻璃纤维缠绕环氧树脂基压力容器的圆柱体上截取试样,通过拉伸、双悬臂梁(DCB)和三点端部开口弯曲(3ENF)试验对复合材料层合板进行基本力学性能评估,编写进用于判断材料失效的VUMAT用户子程序中,并采用有限元软件Abaqus分析不同深度断裂缺陷对复合材料压力容器应力、应变的影响。
1. 试样制备
将圆柱形模具固定在旋转操作箱上,使旋转操作箱不停转动,将含有环氧树脂预浸料的玻璃纤维缠绕在模具表面,缠绕到一定厚度后脱模,形成完整的压力容器,如图1(a)所示。该容器的内径、高度和厚度分别为1600、1620和18 mm,容积为3 m3,设计压力为0.375 MPa。玻璃纤维部分包括玻璃纤维缠绕线层(Y)、布层(C)和毡层(M),3种纤维层的铺层顺序为M/Y/C/Y/C/Y/C/M/Y/C/Y/C/Y/C/M。按照ASTM D3039/D3039M-08标准,采用水切割方法在压力容器筒体段的0°、90°和45°方向截取试样,试样尺寸分别为220 mm × 25 mm × 18 mm、220 mm × 25 mm × 18 mm和160 mm × 25 mm × 18 mm,每种试样各取5份。
根据裂纹产生原因不同,复合材料可以划分为3种分层模式:Ⅰ型断裂又称张开型裂纹,由与层合板面垂直的力产生,分层扩展方向与外载荷方向垂直,Ⅱ型和Ⅲ型分别称滑移型断裂和撕开型断裂,与面内剪应力有关,分别由平行和垂直于裂纹扩展方向的面内力产生。考虑到复合材料压力容器是由内压而产生的力学行为特点,主要考虑Ⅰ型和Ⅱ型两种基本断裂模式,测量这两种类型分层的断裂韧性,从而确定该材料的层间性能。根据ISO 15024标准、ASTM D5528-13标准和航空工业HB7402-1996标准,采用DCB试验方法测量I型复合材料的层间断裂韧性,如图1(b)所示。试样采用单向层板,厚度(h)为15 mm,长度(l)为220 mm,预制分层长度(a0)为25 mm,试样宽度(b)为25 mm。根据ASTM D7905标准、日本JSAK7086标准和中国航空工业HB7403-1996标准,采用3ENF试验方法测定Ⅱ型复合材料的分层断裂韧性,如图1(c)所示。试样采用单向层板,厚度(h)为15 mm,长度(l)为100 mm,预制分层长度(a0)为25 mm,试样宽度(b)为10 mm。
2. 基本力学性能测试
力学性能测试在INSTRON-8032电液伺服动态万能材料疲劳试验机上进行,试验机的最大载荷为100 kN,测试速度为2 mm/min。通过拉伸试件测试获得试样的刚度参数:横向拉伸模量Ex、纵向拉伸模量Ey、泊松比vxy和剪切模量Gxy。由于所取试样的纤维方向与其主方向的夹角θ为12°,通过偏轴公式式(1)~式(4)得到最终的刚度参数:纤维主方向拉伸模量E1、与纤维垂直方向拉伸模量E2、泊松比v12和剪切模量G12。
1Ex=1E1cos4θ+(1G12−2v12E1)sin2θcos2θ+1E2sin4θ (1) 1Ey=1E1sin4θ+(1G12−2v12E1)sin2θcos2θ+1E2cos4θ (2) −vxyEx=−v12E1(sin4θ+cos4θ)+(1E1+1E2−1G12)sin2θcos2θ (3) 1Gxy=1G12(sin4θ+cos4θ)+4(1+2v21E1+1E2−12G12)sin2θcos2θ (4) 根据修正梁理论,DCB试验中Ⅰ型分层试样的层间应变能释放率
RⅠ\!\!c=3Fpcd2b(a+Δ) (5) 其中
F=1−310(da)2−23(dl1a2) (6) 式中:pc为试样端部施加的载荷,d为端部位移即分层开口位移,a为分层长度,F为考虑大位移和加载块影响时的修正因子,Δ为分层长度的修正量。对于Ⅱ型分层,应变能释放率
RⅡ\!c=9a20pcd2b(38l30+3a30) (7) 式中:l0为弯曲加载跨距的一半,d为端部位移即弯曲试样中心挠度。试样采用预制分层的单向层板,通过加载平台施加载荷。
3. 有限元模拟方法
如图2所示,利用ABAQUS建立含有不同深度断裂缺陷的复合材料壳体的有限元分析模型,预测其在内压载荷作用下的应力和应变响应。同时,编写用于定义材料参数及判断材料失效的VUMAT用户子程序,采用Hashin 3D失效准则,并考虑材料的刚度递减规律,参数的选取参考文献[19]。复合材料壳体的厚度为18 mm,根据纤维的缠绕工艺特点,将复合材料壳体分为6层。断裂缺陷的长度、宽度和厚度分别为50、2和3 mm。第1层断裂缺陷深度为18 mm,第2层为15 mm,第3层为12 mm,第4层为9 mm,第5层为6 mm,第6层为3 mm。内部压力分别设定为0.1、0.2和0.3 MPa。
4. 结果与讨论
玻璃纤维增强环氧树脂基复合材料的刚度和强度参数列于表1和表2中,表中数据为5个试样的平均值和方差。表中:E1、E2和E3分别为纤维方向、垂直纤维方向和试样厚度方向的拉伸模量,G12、G13和G23为3个面内的剪切模量,v12、v13和v23为3个面内的泊松比,XT和XC分别为纤维方向的拉伸和压缩强度,YT和YC分别为垂直纤维方向的拉伸和压缩强度,ZT和ZC分别为试样厚度方向的拉伸和压缩强度,S12、S13和S23表示3个面内的剪切强度,其中厚度方向的相关参数取自文献[20-21]。可以看到,复合材料主向的拉伸强度为(222.7 ± 18)MPa,弹性模量为39.39 GPa,横向的拉伸强度为(136 ± 22)MPa,弹性模量为18.1 GPa,层间强度Ⅰ型和Ⅱ型断裂韧性分别为(4.67 ± 0.24)kJ/m2和(4.98 ± 0.26)kJ/m2。
图3为不同深度断裂损伤的应力分布结果,H为断裂深度,pH为内压。如图3(a)所示,当断裂损伤深度为18 mm,内压载荷为0.1 MPa时,应力范围为3.621~13.140 MPa;内压载荷为0.2 MPa时,应力范围为5.002~21.840 MPa;内压载荷为0.3 MPa时,应力范围为5.116~28.880 MPa。最大应力分布在筒体的底端,主要是由于此处不仅受到内压载荷的作用,还受到位移的限制,并承受着筒体的整体重力。含有不同断裂深度壳体的应力结果均有此特点。如图3(b)所示,当断裂损伤深度为15 mm,内压载荷为0.1、0.2和0.3 MPa时,筒体底端的应力高达8.614、13.710和17.020 MPa,而断裂位置处的应力相对较小,主要是由于断裂缺陷位置与其相邻的壳体无接触,内部受到的载荷并没有传递到最外层壳体上,导致最外层壳体断裂缺陷位置处的应力较小。如图3(c)所示,当断裂损伤深度为12 mm,内压载荷为0.1、0.2和0.3 MPa时,缺陷位置处的应力分别为6.664、11.190和12.860 MPa,均小于最大应力。如图3(d)和图3(e)所示,当断裂损伤的深度分别为9和6 mm时,应力的分布较为相似,依然具有以上特点,不同的是随着深度减小,最大应力基本保持不变,内压载荷为0.1 MPa时,最大应力分别为8.066和8.130 MPa。这主要是因为断裂缺陷距离内压载荷面越来越远,对整个壳体应力分布的影响也越来越小,直到断裂缺陷出现在最外层壳体上,如图3(f)所示,当内压分别为0.1、0.2和0.3 MPa时,最大应力分别为8.122、11.410和15.150 MPa,断裂缺陷对整个壳体的应力分布影响最小。
由图3可以看出,随着内压载荷增大,最大应力也不断增大,且最大应力出现在壳体底端附近,缺陷位置处的应力相对较小。统计上述不同内压载荷下的最大应力值,如图4所示。可以明显看出,随着内压增加,应力也不断增大。当断裂缺陷发生在第1层时,复合材料压力容器的应力最大。内压为0.1 MPa时,最大Mises应力高达13.14 MPa;内压为0.2 MPa时,最大Mises应力高达21.84 MPa;内压为0.3 MPa时,最大Mises应力高达28.8 MPa。当断裂缺陷逐渐靠近外壁,即与内压作用面较远时,最大应力趋于不变,保持平稳状态,断裂缺陷在第1层产生的最大应力是在第6层时的2倍。
为了探究壳体外壁应变的变化情况,选择穿过断裂缺陷中心位置处的圆周路径,得到含有不同深度断裂缺陷的压力容器壳体的应变曲线,如图5所示。在图5中,
ε CS为周向应变,ε LS为纵向应变,横坐标零处为缺陷中心位置,横轴为与缺陷中心的周向距离。如图5(a)所示,当内压为0.1 MPa时,周向应变在3.0 × 10−4~1.50 × 10−3之间波动,曲线的形状与反正弦曲线相似,因此可进一步拟合,得到周向应变与距断裂缺陷距离的关系函数,而纵向应变在6.3 × 10−4~8.3 × 10−4之间波动。随着距离增加,纵向应变比较稳定,尤其是当断裂缺陷位于第1层和第6层时,几乎稳定在6.5 × 10−4~8.0 × 10−4范围内,这是由研究路径决定的,因研究路径位于壳体的最外层,当断裂缺陷位于第1层时,缺陷对研究路径上的应变影响较小,而缺陷位于第6层时,则是内压对最外层的影响小而导致应变分布较为均匀。如图5(b)和图5(c)所示,当内压为0.2和0.3 MPa时,周向和纵向应变的范围增大,且最大值也相应增加。当内压为0.2 MPa时,周向应变在5.0 × 10−4~2.50 × 10−3之间波动,纵向应变在1.35 × 10−3~1.70 × 10−3之间波动;当内压为0.3 MPa时,周向应变在5.0 × 10−4~3.50 × 10−3之间波动,纵向应变在2.20 × 10−3~2.60 × 10−3之间波动。综上所述,当内压为0.3 MPa时,周向应变的变化范围最大,可根据应变值确定断裂缺陷的位置和深度。因此,对该内压下的周向应变与距离曲线进行拟合,拟合采用waveform的sine公式,得到不同深度断裂缺陷的周向应变(
ε CS)与缺陷周向位置(x)的关系式(即式(8)~式(13)),可以较准确地确定不同深度断裂缺陷的周向位置和周向应变的关系。在实际工程实践中,采用应变检测方法确定壳体的应变分布,再通过关系式即可确定缺陷的周向位置和深度,为复合材料压力容器和管道检测提供一定的数据基础。18mm(1st):εCS=1989.3+1272.4sin(πx−2414.02546.6) (8) 15mm(2nd):εCS=2086.3+1390.2sin(πx−2421.22546.5) (9) 12mm(3rd):εCS=2026.8+1304.3sin(πx−2548.02543.0) (10) 9mm(4th):εCS=2027.8+1345.4sin(πx−2421.32552.7) (11) 6mm(5th):εCS=2035.4+1340.3sin(πx−2509.82564.6) (12) 3mm(6th):εCS=2045.2+1333.6sin(πx−2509.82564.6) (13) 5. 结 论
对玻璃纤维增强环氧树脂基复合材料层合板进行了基本力学性能试验,并利用Abaqus有限元软件建立了含有不同深度断裂缺陷的复合材料壳体的三维有限元模型,将所得的基本力学性能参数用于有限元模拟中,通过对模拟结果与实验结果进行对比分析,得到以下结论:
(1)复合材料主向的拉伸强度和弹性模量均大于横向,层间强度Ⅰ型和Ⅱ型断裂韧性相差不大;
(2)应力和应变参数可以反映出复合材料壳体断裂缺陷的位置,通过有限元模拟,能够很方便地得到含有断裂缺陷壳体的应力及应变分布情况,复合材料壳体的最大应力在筒体端部位置,缺陷位置处的应力相对较小;
(3)随着内压增加,复合材料壳体的Mises应力也不断增加,当断裂缺陷在第1层时,Mises应力最大,应变随着与缺陷位置距离的增加,呈现规律性变化,通过拟合公式,可以初步判断断裂损伤的周向位置和深度,为复合材料壳体的损伤检测提供一定的数据基础。
-
-
[1] 陈建良, 童永光. 复合材料在压力容器中的应用 [J]. 压力容器, 2001, 18(6): 47–50. doi: 10.3969/j.issn.1001-4837.2001.06.013CHEN J L, TONG Y G. The application of composite materials to pressure vessels [J]. Pressure Vessel Technology, 2001, 18(6): 47–50. doi: 10.3969/j.issn.1001-4837.2001.06.013 [2] HER S C, LIANG Y C. The finite element analysis of composite laminates and shell structures subjected to low velocity impact [J]. Composite Structures, 2004, 66(1/2/3/4): 277–285. doi: 10.1016/j.compstruct.2004.04.049 [3] 郑津洋, 傅强, 开方明, 等. 轻质高压贮氢容器的现状及发展趋势 [J]. 太阳能学报, 2004, 25(5): 576–581. doi: 10.3321/j.issn:0254-0096.2004.05.003ZHENG J Y, FU Q, KAI F M, et al. The state of art of lightweight high-pressure hydrogen storage tanks [J]. Acta Energiae Solaris Sinica, 2004, 25(5): 576–581. doi: 10.3321/j.issn:0254-0096.2004.05.003 [4] GRIMSLEY B W, CANO R J, JOHNSTON N J, et al. Hybrid composites for LH2 fuel tank structure [C]//33rd International SAMPE Technical Conference. Seattle, WA, 2001: 1224−1235. [5] MARTINS A T, ABOURA Z, HARIZI W, et al. Structural health monitoring for GFRP composite by the piezoresistive response in the tufted reinforcements [J]. Composite Structures, 2019, 209: 103–111. doi: 10.1016/j.compstruct.2018.10.091 [6] HAN M G, CHANG S H. Failure analysis of a Type Ⅲ hydrogen pressure vessel under impact loading induced by free fall [J]. Composite Structures, 2015, 127: 288–297. doi: 10.1016/j.compstruct.2015.03.027 [7] HALE K F. An optical-fibre fatigue crack-detection and monitoring system [J]. Smart Materials and Structures, 1992, 1: 156–161. doi: 10.1088/0964-1726/1/2/009 [8] GHIMIRE M, WANG C J, 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 [9] LU Z J, BLAHA F A. A fiber optic strain and impact sensor system for composite materials [C]//Proceedings Volume 1170, Fiber Optic Smart Structures and Skins II. Boston, United States: SPIE, 1989. [10] 赵海涛. 基于光纤传感技术的复合材料结构全寿命健康监测研究[D]. 哈尔滨: 哈尔滨工业大学, 2008.ZHAO H T. Research on whole life health monitoring of composites structure based on fiber optic sensing technology [D]. Harbin: Harbin Institute of Technology, 2008. [11] TOUGHIRY M. Examination of the nondestructive evaluation of composite gas cylinders [R]. Austin, TX: The Nondestructive Testing Information Analysis Center, 2002. [12] LEH D, SAFFRÉ P, FRANCESCATO P, et al. A progressive failure analysis of a 700-bar type IV hydrogen composite pressure vessel [J]. International Journal of Hydrogen Energy, 2015, 40(38): 13206–13214. doi: 10.1016/j.ijhydene.2015.05.061 [13] WANG L, ZHENG C X, LUO H Y, et al. Continuum damage modeling and progressive failure analysis of carbon fiber/epoxy composite pressure vessel [J]. Composite Structures, 2015, 134: 475–482. doi: 10.1016/j.compstruct.2015.08.107 [14] LIU P F, ZHENG J Y. Progressive failure analysis of carbon fiber/epoxy composite laminates using continuum damage mechanics [J]. Materials Science and Engineering: A, 2008, 485(1/2): 711–717. doi: 10.1016/j.msea.2008.02.023 [15] 杨斌, 章继峰, 周利民. 玻璃纤维-碳纤维混杂增强PCBT复合材料层合板的制备及低速冲击性能 [J]. 复合材料学报, 2015, 32(2): 435–443. doi: 10.13801/j.cnki.fhclxb.201502.004YANG 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 Materiae Compositae Sinica, 2015, 32(2): 435–443. doi: 10.13801/j.cnki.fhclxb.201502.004 [16] OKABE Y, TSUJI R, TAKEDA N. Application of chirped fiber Bragg grating sensors for identification of crack locations in composites [J]. Composites Part A: Applied Science and Manufacturing, 2004, 35(1): 59–65. doi: 10.1016/j.compositesa.2003.09.004 [17] 肖飚, 杨斌, 胡超杰, 等. 基于埋入式应变片的纤维缠绕压力容器的健康监测 [J]. 高压物理学报, 2019, 33(4): 043401. doi: 10.11858/gywlxb.20190726XIAO B, YANG B, HU C J, et al. Structural health monitoring of filament wound pressure vessel by embedded strain gauges [J]. Chinese Journal of High Pressure Physics, 2019, 33(4): 043401. doi: 10.11858/gywlxb.20190726 [18] KANERVA M, ANTUNES P, SARLIN E, et al. Direct measurement of residual strains in CFRP-tungsten hybrids using embedded strain gauges [J]. Materials and Design, 2017, 127: 352–363. doi: 10.1016/j.matdes.2017.04.008 [19] BAI H, YANG B, HUI H, et al. Experimental and numerical investigation of the strain response of the filament wound pressure vessels subjected to pressurization test [J]. Polymer Composites, 2019, 40(11): 4427–4441. doi: 10.1002/pc.25304 [20] ALSHAHRANI R F, MERAH N, KHAN S M A, et al. On the impact-induced damage in glass fiber reinforced epoxy pipes [J]. International Journal of Impact Engineering, 2016, 97: 57–65. doi: 10.1016/j.ijimpeng.2016.06.002 [21] YANG B, HE L, GAO Y. Simulation on impact response of FMLs: effect of fiber stacking sequence, thickness, and incident angle [J]. Science and Engineering of Composite Materials, 2018, 25(3): 621–631. doi: 10.1515/secm-2016-0226 -