Dynamic and Static Tensile Mechanical Properties of Glass Fiber Reinforced Plastics
-
摘要: 为系统地研究环氧树脂玻璃钢在静、动态拉伸载荷作用下的力学性能,采用材料测试系统和分离式霍普金森拉杆对材料进行拉伸试验,获得0.001~0.1 s−1及1128~1840 s−1应变率下的应力-应变曲线和相应的力学参数。结果表明,动态加载下环氧树脂玻璃钢的应变率增强效应较为明显。为此,引入动态增强因子描述环氧树脂玻璃钢在高应变率下力学性能的增强。采用扫描电镜对损伤断面进行观测,发现动态加载下纤维束平整断裂,而非静态加载下纤维拔出失效。相较于静态加载,动态拉伸载荷作用下玻璃钢的基体-纤维界面断裂韧度更高。基于环氧树脂玻璃钢在动态拉伸下的力学响应,引入宏观损伤累积量,建立一种考虑损伤的非线性拉伸本构模型。拟合结果表明,该模型整体上可以反映环氧树脂玻璃钢在动态拉伸载荷作用下的力学响应。Abstract: To systematically study the mechanical properties of glass fiber reinforced plastics, the stress-strain curves and corresponding mechanical parameters of glass fiber reinforced plastic with strain rates of 0.001–0.1 s−1 and 1128–1840 s−1 were obtained by material testing system (MTS) and split Hopkinson tensile bar (SHTB). The results show that the glass fiber reinforced plastic exhibits significant strain rate strengthening effect under dynamic loading. Therefore, the dynamic increasing factor was introduced to describe the strengthening effect of strain rate on mechanical properties of glass fiber reinforced plastics. Scanning electron microscopy (SEM) was used to observe and analyze the damaged section, and it was found that the fiber bundle was flatly fracture under dynamic loading, rather than failed by pulling out under static loading. In addition, the results show that the fracture toughness of the fiber-matrix interface under dynamic tensile loading is higher than that under static loading. Finally, based on the dynamic tensile mechanical response of glass fiber reinforced plastic, a nonlinear tensile constitutive model considering damage was established by introducing the macroscopic damage accumulation. Compared with the experimental results, the model as a whole can be used to characterize the dynamic tensile mechanical response of glass fiber reinforced plastic.
-
It is very important for mining and civil construction to predict the morphology distribution of cracks induced by blasting. Hence, many researchers have paid their attention to dynamic fracture behavior of rocks due to drilling and blasting operations[1-3]. A number of experiments and numerical simulations have been conducted to investigate the blasting-induced fractures in the near borehole zone as well as in the far field[4-5]. In order to gain high fidelity in simulating the complex responses of rocks subjected to blasting loading, a realistic constitutive model is required. In the last 20 years, various macro-scale material models have been proposed, from relatively simple ones to more sophisticated, and their capabilities in describing actual nonlinear behavior of material under different loading conditions have been evaluated[6].
During blasting operation, chemical reactions of explosive in borehole occur rapidly, and instantaneously a shock/stress wave applies to borehole wall. Initially, a crushed zone around the borehole is developed by the shock/stress wave. Then, a radial shock/stress wave propagates away from borehole, and its magnitude decreases. Once the radial shock/stress drops below the local dynamic compressive strength, no shear damage occurs. At the same time, a tensile tangential stress with enough strength can be developed behind the radial compressive stress wave, which results in an extension of the existing flaws or a creation of new radial cracks. If there is a nearby free boundary, the incident compressive stress wave changes to a tensile stress wave upon reflection, and reflects back into the rock. In this case, if the dynamic tensile strength of rock is exceeded, spall cracks appear close to the free boundary.
The purpose of this paper is to conduct a numerical study on borehole blasting-induced fractures in rocks. First, a dynamic constitutive model for rocks based on the previous work of concrete[7] is briefly described, and the values of various parameters in the model for granite are estimated. The model is then employed to simulate the borehole blasting-induced fractures in granitic rocks. Comparisons between the numerical results and the experimental observations are made, and a discussion is given.
1. Dynamic Constitutive Model for Rocks
A number of models for concrete-like materials, such as TCK model[8], HJC model[9], RHT model[10], K&C model[11], have been developed. The sophisticated numerical models are increasingly used as they are capable of describing the material behavior under high strain rate loading. However, these models have been found to have some serious flaws, and cannot predict the experimentally observed crack patterns or exhibit improper behavior under certain loading conditions[7, 12-15].
In the following, a dynamic constitutive model for rocks is briefly described according to equation of state (EOS) and strength model, based on the previous work on concrete[7].
1.1 EOS
A typical form of EOS is the so-called p-
α relation, which is proved to be capable of representing brittle material’s response behavior at high pressures, and it allows for a reasonably detailed description of the compaction behavior at low pressure ranges as well, as shown schematically in Fig.1.pcrush corresponds to the pore collapse pressure beyond which plastic compaction occurs, andplock is the pressure when porosityα reaches 1,ftt is the tensile strength,ρ0 is the initial density,ρs0 refers to the density of the initial solid.Figure 1. Schematic diagram of EOS[7]The EOS for compression (p≥0) is given by
p=K1ˉμ+K2ˉμ2+K3ˉμ3 (1) where p denotes pressure, K1, K2, K3 are constants, and
ˉμ is defined byˉμ=ραρ0α0−1=αα0(1+μ)−1 (2) where
ρ is the current density,μ=ρ/ρ0−1 specifies the volumetric strain,α0 =ρ s0/ρ 0 andα=ρs/ρ represent the initial porosity and the current porosity, respectively.ρs refers to the density of fully compacted solid. Physically,α is a function of the hydro-static pressure p, and is expressed asα=1+(α0−1)(plock−pplock−pcrush)n (3) where n is the compaction exponent.
When material withstands hydro-static tension, the EOS for tension (p<0) is given by
p=K1ˉμ (4) α=α0(1+p/K1)/(1+μ) (5) 1.2 Strength Model
The strength model takes into account various effects, such as pressure hardening, damage softening, third stress invariant (Lode angle) and strain rate. The strength surface Y, shown schematically in Fig.2, can be written as[7]
Y={3(p+ftt)R(θ,e) p<0[3ftt + 3p(fcc−3ftt)/fcc]R(θ,e) 0⩽p⩽fcc/3{fcc+Bfc′[p/fc′−fcc/(3fc′)]N}R(θ,e) p>fcc/3 (6) where p is the hydro-static pressure, parameters B and N are constants, R(θ,e) is a function of the Lode angle θ and the tensile-to-compressive meridian ratio e,
fc′ is the static uni-axial compressive strength, the compressive strengthfcc and the tensile strengthftt are defined byfcc=fc′Dm_tηc (7) ftt=ftDtηt (8) where
ft is the static uni-axial tensile strength.Dm_t is the compression dynamic increase factor due to strain rate effect only, and can be expressed as[7, 16]Dm_t=(Dt−1)ft/fc′+1 (9) where
Dt is the tension dynamic increase factor determined byDt={tanh[(lg˙ε˙ε0−Wx)S][FmWy−1]+1}Wy (10) where
Fm ,Wx ,Wy andS are experimental constants,˙ε is the strain rate, and˙ε0 is the reference strain rate, usually taken˙ε0=1.0 s−1 .ηc is the damage function for compression, which can be expressed asηc={l+(1−l)η(λ)λ⩽λmr+(1−r)η(λ)λ>λm (11) where l and r are constants[7],
λm is the value of shear damage (λ ) when strength reaches its maximum value under compression.η(λ) is defined asη(λ)=aλ(λ−1)exp(−bλ) (12) in which a and b can be determined by setting
η(λ) = 1 and∂η∂λ=0 whenλ =λ m.ηt is the damage function for tension which can be written asηt=[1+(c1εtεfrac)3]exp(−c2εtεfrac)−εtεfrac(1+c31)exp(−c2) (13) where c1 and c2 are constants[7],
εt denotes the tensile strain andεfrac is the fracture strain.The residual strength (
rfc′ ) surface for rocks, shown schematically in Fig.2, can be obtained from Eq.(6) by settingftt=0 andfcc=rfc′ , vizY={3pR(θ,e) 0<p⩽rfc′/3{rfc′+Bfc′[p/fc′−rfc′/(3fc′)]N}R(θ,e) p>rfc′/3 (14) 2. Numerical Simulations
Granite is selected for investigating the dynamic fractures which result from borehole blast loading.
2.1 Evaluation of Various Parameters in the Model
Table 1 lists the values of the various parameters used in the dynamic constitutive model for granite. As to how to determine the values of the various parameters in the model, more details are presented in [7, 15-17].
p-α relation ρ0/(kg∙m−3) pcrush/MPa plock/GPa n K1/GPa K2/TPa K3/TPa 2660 50.5 3 3 25.7 −3 150 Strength surface Strain rate effect fc′/MPa ft/MPa B N G/GPa Fm Wx 161.5 7.3 2.59 0.66 21.9 10 1.6 Strain rate effect Shear damage Wy S ˙ε0/s−1 λs λm l r 5.5 0.8 1.0 4.6 0.3 0.45 0.3 Lode effect Tensile damage e1 e2 e3 c1 c2 εfrac 0.65 0.01 5 3 6.93 0.007 Fig.3 shows the comparison of the strength surface between Eq.(6) (with B=2.59, N=0.66) and the triaxial test data for granite[17]. It can be seen from Fig.3 that a good agreement is obtained. Similarly, Fig.4 shows the tensile strengths/dynamic increase factor obtained by Eq.(10) and the test results of various rocks at different strain rates[18-23]. It is clear from Fig.4 that a good agreement is achieved.
Figure 3. Comparison of the strength surface between Eq.(6) (with B=2.59, N=0.66) and the triaxial test data for granite[17]2.2 Numerical Results
In the following, numerical simulations are carried out for the response of the granite targets subjected to borehole blasting loading. The dynamic fracture behavior of two kinds of granite samples are studied, namely, cylindrical sample as reported in the literature and square sample as examined in our own laboratory.
2.2.1 Cylindrical Rock Sample
In consideration of the sizes of the cylindrical granite samples prepared for laboratory-scale blasting experiments by Dehghan Banadaki and Mohanty[17] (with a diameter of 144 mm, a height of 150 mm and a borehole diameter of 6.45 mm), a circular plane strain model with an outer diameter of 144 mm is made in our simulation, as shown in Fig.5, being a scaled close-up view of the borehole region. Multi-material Euler solver is used for modeling PETN explosive, polyethylene and air. Lagrangian descriptions are used for modeling the copper tube and granite.
The material model and the properties of PETN explosive, polyethylene, air and copper tube used in the simulation are given in Ref.[17]. The values of various parameters in the constitutive model for granite are listed in Table 1.
Fig.6 shows the comparison of the peak pressures between our simulation results of the present model, the numerical results[17], and the experimental results by Dehghan Banadaki and Mohanty[24]. It can be seen from Fig.6 that a good agreement is obtained.
In order to characterize the damping behavior of stress in granite, the peak pressure p in granite is expressed in an exponential form as
pp0=(dd0)−γ (15) where p0 is the peak pressure on the borehole wall, d0 is the initial radius of the borehole, d is the distance from the center point of the borehole,
γ is an index. It is evident from Fig.6 that Eq.(15) withγ =1.6 correlates well with the experimental results.Fig.7 shows the comparison between the crack patterns predicted numerically based on the present model and the one observed experimentally in the cylindrical granite sample[17]. It is clear that a relatively good agreement on the crack pattern is obtained. It is also clear that the stress waves produce three distinct crack regions in the cylindrical granite sample: densely populated smaller cracks around the borehole, a few large radial cracks propagating towards the outer boundary, and circumferential cracks close to the sample boundary which are due to the reflected tension stress.
Figure 7. Comparison of the crack patterns between the numerical prediction and the experiment of the cylindrical granite sample[17]In order to make an assessment of the contributions of both the compression/shear stress and the tensile stress to the crack patterns, Fig.8 shows the numerically predicted crack pattern which results from the tension stress only. Fom Fig.8 and Fig.7(a), it is apparent that there are virtually few changes in crack patterns, both having the large radial and circumferential cracks caused by the same tensile stress, however, the smaller cracks around the borehole in Fig.8 are much less than those in Fig.7(a). In another word, crack patterns are mainly caused by tensile stress, and smaller cracks around borehole are created largely by compression/shear stress.
2.2.2 Square Granite Sample
The laboratory-scale single-hole blasting tests are also carried out in order to validate further the accuracy and the reliability of the present model. Two square granite samples with a side length of 400 mm and a height of 100 mm are employed in the experiments. The borehole diameter is 4 mm and a series of concentric rings is drawn on the top surface of the samples (see Fig.9), so that the damage regions induced by detonation can be assessed visually.
A cylindrical RDX explosive enclosed by an aluminum sheath (see Fig.10) is tightly installed in the borehole of the No.1 sample, while an unwrapped RDX explosive is inserted into the borehole of the No.2 sample. The density of the RDX is 1700 kg/m3, and the material model and properties of the RDX explosive and the aluminum sheath used in the experiments are given in ref.[25]. The values of the various parameters in the constitutive model for granite are listed in Table 1.
Fig.11 and Fig.12 show the comparisons of the crack patterns between the numerical predictions from the present model and the ones observed experimentally in the square granite samples. It can be seen from Fig.11 and Fig.12 that good agreements are obtained. It should be mentioned here that No.1 sample receives less damage due to less RDX explosive used in the test, and that No.2 sample is broken up into four major pieces due to more RDX explosive employed in the experiment. Severe damages and small cracks are induced in the vicinity of the boreholes of both samples, as can be seen clearly from Fig.11(b) and Fig.12(b).
3. Conclusions
A numerical study on the borehole blasting-induced fractures in rocks is conducted in this paper, using a dynamic constitutive model developed previously for concrete. Two kinds of granite rocks are simulated numerically, one in the cylindrical form and the other in the square form. The numerical results are compared with the corresponding experiments. Main conclusions can be drawn as follows.
(1) The crack patterns predicted numerically from the present model are found to be in good agreement with the experimental observations, both in cylindrical and square granite samples subjected to borehole blasting loading.
(2) The peak pressures predicted numerically based on the present model are found to be in good agreement with the test data.
(3) Crack pattern observed experimentally in the rock sample is mainly caused by the tensile stress, while the smaller cracks in the vicinity of the borehole are created largely by compression/shear stress.
(4) The consistency between the numerical results and the experimental observations demonstrates the accuracy and reliability of the present model. Thus the model can be used in the numerical simulations of the response and the failure of rocks under blasting loading.
-
表 1 环氧树脂玻璃钢的准静态拉伸力学参数
Table 1. Quasi-static tensile mechanical properties of glass fiber reinforced plastics
Strain rate/
s−1Young’s modulus/
GPaTensile strength/
MPaFailure strain/
%0.001 16.66 526.74 4.57 0.01 16.70 565.75 4.90 0.1 16.79 586.69 4.98 表 2 不同应变率下环氧树脂玻璃钢的拉伸力学参数
Table 2. Tensile mechanical properties of glass fiber reinforced plastics at various strain rates
Strain rate/
s−1Tensile strength/
MPaFailure strain/
%0.001 526.74 4.57 1128 527.15 3.58 1366 528.04 3.47 1547 533.99 3.51 1735 549.33 3.54 1840 568.09 4.12 -
[1] DEAN J, S-FALLAH A, BROWN P M, et al. Energy absorption during projectile perforation of lightweight sandwich panels with metallic fibre cores [J]. Composite Structures, 2011, 93(3): 1089–1095. doi: 10.1016/j.compstruct.2010.09.019 [2] LIU Y J, JIANG Z, WEN H M. Predicting impact induced delamination of FRP laminates [J]. International Journal of Impact Engineering, 2020, 137: 103436. doi: 10.1016/j.ijimpeng.2019.103436 [3] 田宏伟, 郭伟国. 平纹机织玻璃纤维增强复合材料面内压缩力学行为及破坏机制 [J]. 复合材料学报, 2010, 27(2): 133–139.TIAN H W, GUO W G. In-plane compressive mechanics behavior and failure mechanism for SW200/LWR-2 glass-woven composite [J]. Acta Materiae Compositae Sinica, 2010, 27(2): 133–139. [4] 朱文墨, 李刚, 杨小平, 等. 连续纤维增强树脂复合材料纵向压缩强度预测模型的发展及其影响因素 [J]. 复合材料学报, 2020, 37(1): 1–15.ZHU W M, LI G, YANG X P, et al. Development of prediction model and influencing factors of longitudinal compressive strength for continuous fiber reinforced polymer composites [J]. Acta Materiae Compositae Sinica, 2020, 37(1): 1–15. [5] 夏源明, 王兴, 杨报昌. 单向玻璃纤维增强环氧树脂在冲击拉伸时的一维本构方程 [J]. 复合材料学报, 1994, 11(4): 110–116.XIA Y M, WANG X, YANG B C. One dimensional constitutive equation of unidirectional glassfiber reinforced epoxy under tensile impact [J]. Acta Materiae Compositae Sinica, 1994, 11(4): 110–116. [6] SHOKRIEH M M, OMIDI M J. Tension behavior of unidirectional glass/epoxy composites under different strain rates [J]. Composite Structures, 2009, 88(4): 595–601. doi: 10.1016/j.compstruct.2008.06.012 [7] 吴健, 王纬波, 李泓运, 等. 中等应变率下单向玻璃纤维增强环氧树脂基复合材料的剪切本构关系 [J]. 复合材料学报, 2018, 35(2): 304–310.WU J, WANG W B, LI H Y, et al. Shear constitutive relationship of unidirectional glass fiber reinforced epoxy composites under intermediate strain rate [J]. Acta Materiae Compositae Sinica, 2018, 35(2): 304–310. [8] 刘子尚, 杨喆, 魏延鹏, 等. 单向增强玻璃钢复合材料静/动态拉伸实验研究 [J]. 爆炸与冲击, 2019, 39(9): 093101.LIU Z S, YANG Z, WEI Y P, et al. Static/dynamic tensiletest of unidirectional reinforced GFRP composites [J]. Explosion and Shock Waves, 2019, 39(9): 093101. [9] FERESHTEH-SANIEE F, MAJZOOBI G H, BAHRAMI M. An experimental study on the behavior of glass-epoxy composite at low strain rates [J]. Journal of Materials Processing Technology, 2005, 162: 39–45. [10] 张磊, 孙清, 王虎长, 等. E玻璃纤维增强环氧树脂基复合材料力学性能试验研究 [J]. 电力建设, 2010, 31(9): 118–121. doi: 10.3969/j.issn.1000-7229.2010.09.030ZHANG L, SUN Q, WANG H C, et al. Experimental study on the mechanical properties of E-glass fiber/epoxy composite material [J]. Electric Power Construction, 2010, 31(9): 118–121. doi: 10.3969/j.issn.1000-7229.2010.09.030 [11] 张硕, 姚宁, 吴继平, 等. 玻璃纤维增强环氧树脂复合材料的力学性能 [J]. 电工材料, 2016(1): 11–14, 19.ZHANG S, YAO N, WU J P, et al. Mechanical properties of glass fiber reinforced epoxy resin composite material [J]. Electrical Engineering Materials, 2016(1): 11–14, 19. [12] CHEN W S, MENG Q F, HAO H, et al. Quasi-static and dynamic tensile properties of fiberglass/epoxy laminate sheet [J]. Construction and Building Materials, 2017, 143: 247–258. doi: 10.1016/j.conbuildmat.2017.03.074 [13] NAIK N K, YERNAMMA P, THORAM N M, et al. High strain rate tensile behavior of woven fabric E-glass/epoxy composite [J]. Polymer Testing, 2010, 29(1): 14–22. doi: 10.1016/j.polymertesting.2009.08.010 [14] STAAB G H, GILAT A. High strain rate response of angle-ply glass/epoxy laminates [J]. Journal of Composite Materials, 1995, 29(10): 1308–1320. doi: 10.1177/002199839502901003 [15] GAO G F, LI Y C. Dynamic behavior of a woven glass-fiber-reinforced polymer composite at high strain rates and its dynamic constitutive relationship [J]. Mechanics of Advanced Materials and Structures, 2017, 24(13): 1086–1093. doi: 10.1080/15376494.2016.1227493 [16] 许沭华, 王肖钧, 张刚明, 等. Kevlar纤维增强复合材料动态压缩力学性能实验研究 [J]. 实验力学, 2001, 16(1): 26–33.XU S H, WANG X J, ZHANG G M, et al. Experimental investigation on the dynamic compression properties of kevlar fiber-rainforced composite laminates [J]. Journal of Experimental Mechanics, 2001, 16(1): 26–33. [17] HU J X, YIN S, YU T X, et al. Dynamic compressive behavior of woven flax-epoxy-laminated composites [J]. International Journal of Impact Engineering, 2018, 117: 63–74. doi: 10.1016/j.ijimpeng.2018.03.004 [18] POURNOORI N, SOARES G C, ORELL O, et al. Adiabatic heating and damage onset in a pultruded glass fiber reinforced composite under compressive loading at different strain rates [J]. International Journal of Impact Engineering, 2021, 147: 103728. [19] 哈维尔, 刘问, 孟鑫淼, 等. 基于断裂理论的竹纤维复合材料拉伸强度的尺寸效应研究 [J]. 复合材料科学与工程, 2021(2): 11–18.NCHAMA J, LIU W, MENG X M, et al. Size effect study of tensile strength of bamboo fiber composites on basis of fracture theory [J]. Composites Science and Engineering, 2021(2): 11–18. [20] 阮班超, 史同亚, 王永刚. E玻璃纤维增强环氧树脂基复合材料轴向拉伸力学性能的应变率效应 [J]. 复合材料学报, 2018, 35(10): 2715–2722.RUAN B C, SHI T Y, WANG Y G. Influence of strain rate on tensile mechanical behavior of E glass fiber reinforced epoxy resin composites [J]. Acta Materiae Compositae Sinica, 2018, 35(10): 2715–2722. [21] CHEN C Y, ZHANG C, LIU C L, et al. Rate-dependent tensile failure behavior of short fiber reinforced PEEK [J]. Composites Part B: Engineering, 2018, 136: 187–196. doi: 10.1016/j.compositesb.2017.10.031 [22] NANDLALL D, WILLIAMS K, VAZIRI R. Numerical simulation of the ballistic response of GRP plates [J]. Composites Science and Technology, 1998, 58(9): 1463–1469. doi: 10.1016/S0266-3538(98)00030-X [23] 段茜. 复合应力状态下聚酰胺的粘弹-塑性本构研究 [D]. 太原: 太原理工大学, 2018: 14−18.DUAN Q. Viscoelastic-plastic constitutive model of polyamide under combined stress state [D]. Taiyuan: Taiyuan University of Technology, 2018: 14−18. -