Loading [MathJax]/jax/output/SVG/jax.js
ZHAO Fuqi, XU Peibao, WEN Heming. Influence of Specimen Size in SHPB Tests on Concrete[J]. Chinese Journal of High Pressure Physics, 2018, 32(1): 014101. doi: 10.11858/gywlxb.20170532
Citation: WANG Bi, AN Erfeng, CHEN Pengwan, ZHOU Qiang, GAO Xin. Fabrication of W-Al Energetic Structural Materials by Explosive Consolidation and Investigation of Its Quasi-Static Compression Properties[J]. Chinese Journal of High Pressure Physics, 2019, 33(6): 063401. doi: 10.11858/gywlxb.20190753

Fabrication of W-Al Energetic Structural Materials by Explosive Consolidation and Investigation of Its Quasi-Static Compression Properties

doi: 10.11858/gywlxb.20190753
  • Received Date: 01 Apr 2019
  • Rev Recd Date: 13 Apr 2019
  • Issue Publish Date: 25 Sep 2019
  • Nearly fully dense W-Al energetic structural material (ESM) was successfully prepared by explosive sintering with W and Al powder with different particle sizes. It was found that shock wave pressure is the dominant factor for powder densification and the particle size of powder has a significant influence on the final density and microstructure of the compacted ESM. The smaller the W particle size is, the more severely the W particle agglomerate, which hinders the densification, leading to the formation of continuous W phase in the compacted ESM. The maximum compressive strength and failure strain of the sample reach 288 MPa and 20%, respectively. The mechanical properties and fracture mode of the consolidated material depend on the continuous phase, the ESM with continuous Al phase presents low compressive strength and good ductility with anaxial split failure, while the one with continuous W phase shows brittleness and high compressive strength with a shear failure, which is consistent with the properties of Al and W, respectively.

     

  • Dynamic compressive strengths of concrete-like materials are usually obtained by conducting laboratory tests such as split Hopkinson pressure bar (SHPB) tests[1-5].Many empirical formulae for dynamic increase factor (DIF) based on laboratory test data were reported in the related literature Ref.[6-9].

    It has to be mentioned here that the SHPB test data available for concrete materials in the literature were very scattered[7-10] due to the combined effects of strain rate, inertia and specimen size.More comprehensive investigation and discussion of the possible influencing factors can be also obtained in Ref.[1, 2, 7-11].The inertia (confinement) effect has been widely investigated experimentally[12], theoretically[13-14], numerically[5, 15] and an empirical equation for the dynamic increase factor due to inertia (confinement) effect was suggested in Ref.[16], which took no account of the effect of concrete specimen size.As strain rate effects on the compressive strengths of concrete-like materials play an important part in the construction of material constitutive models which, in turn, exert a great deal of influence on the numerical simulations of concrete structures subjected to intense dynamic loadings, it is, therefore necessary to obtain the pure strain rate effect data by eliminating the data due to the effects of inertia (confinement) and specimen size from SHPB tests.

    In the present study, numerical simulations with a rate-independent material model are carried out on the influence of specimen size in SHPB tests on concrete and a new empirical equation for the dynamic increase factors due to inertia (confinement) effect is proposed which takes account of specimen size effect through its volume.Comparisons are made between the results from the numerical simulation and those from the new empirical formula and discussed.

    The computational constitutive model for concrete developed in Ref.[1] is used to simulate the SHPB tests on concrete in the present study.This material model consists of equation of state (EOS), strength model including Lode effect, damage criteria and strain rate effects, etc.

    The porous equation of state is used in the present study and can be expressed as[17]

    ˉμ=ραρ0α01=αα0(1+μ)1 (1)

    where μ=ρ/ρ0-1 is the volumetric strain, in which ρ and ρ0 are respectively the current and initial densities; α=ρs/ρ and α0=ρs0/ρ0 are the current and initial porosities in which ρs and ρs0 are the current and initial densities of solid (fully-compacted) material.

    For μ > 0, concrete material is under compression condition

    {α=max{1,min[α0,1+(α01)(plockpplockpcrush)n]}p=K1ˉμ+K2ˉμ2+K3ˉμ3 (2)

    where K1, K2 and K3 are the bulk moduli for fully compacted concrete material, pcrush is the pressure at which pore collapse occurs, plock is the pressure beyond which concrete material is fully compacted, n is the compaction exponent.

    For μ < 0, concrete material is under tension condition.Hence, the pressure is

    p=K1ˉμ (3)

    The strength surface of concrete can be written in the following form[1]

    Y={3(p+ftt)r(θ,e)p0[3ftt+3p(fcc3ftt)fcc]r(θ,e)0<p<fcc3[fcc+Bfc(pfcfcc3fc)N]r(θ,e)pfcc3 (4)

    where fcc=fc'δm_tηc and ftt=ftδtηt in which fc' and ft are the static compressive and tensile strength, δm_t and δt[2] are the dynamic increase factors due to strain rate effects only in compression and tension, ηc and ηt[18-20] are shape functions which represent shear damage and tensile softening of concrete respectively; B and N are empirical constants; r(θ, e)[21] is the Lode effect with θ and e being the Lode angle, the ratio of the tensile meridian to the compression meridian, respectively.

    By setting ftt=0 and fcc=fc'r (residual strength), the residual strength surface for concrete can be obtained from Eq.(4)

    Y={3pr(θ,e)0<p<fcr3[fcr+Bfc(pfcr3)N]r(θ,e)pfcr3 (5)

    Other information about shear, tensile damage, Lode effect and strain rate effect can be found in Ref.[1-2].

    In this section, the commercial software LS-DYNA3D with user-defined subroutines is used to carry out the numerical study on the influence of specimen size in SHPB tests on concrete using a recently developed constitutive model for concrete[1] by setting the material DIF (δ) equal to 1.The inertia effect dynamic enhancement factor (Rc) can be obtained by dividing the strength results of simulation by quasi-static strength.The values of various parameters used in the material model are listed in Table 1[1].In this case, the numerically obtained strength increment of the studied specimen is attributed to the inertia (confinement) effect only.

    Table  1.  Material parameters of concrete[1]
    Parameters of EOS Parameters of constitutive model
    ρ0/(kg·m-3) ρs0/(kg·m-3) pcrush/MPa plock/GPa n K1/GPa K2/GPa K3/GPa Strength surface Shear damage Tensile damage Lode effect
    fc'/MPa ft/MPa B N G/GPa λs l r λm c1 c2 εfrac e1 e2 e3
    2 400 2 680 15.2 3 3 13.9 30 10 45.6 3.8 1.7 0.7 10.5 4.6 0.45 0.3 0.3 3 6.93 0.007 0.65 0.01 5
     | Show Table
    DownLoad: CSV

    An SHPB system contains incident and transmitter pressure bars with a short specimen between them, as shown in Fig. 1.A stress pulse of trapezium shape is applied to the incident pressure bar as shown schematically in Fig. 2 in which the incident stress starts from 0, quickly rises to the peak value of ppeak at t1, keeps as a constant for t2, then drops back to 0 after t3 and are summarized in Table 2.

    Figure  1.  Finite element model of SHPB system
    Figure  2.  Loading function used in finite element analysis
    Table  2.  Load function parameters for direct compression analyses
    t1/μs t2/μs t3/μs ppeak/MPa
    25 200 25 Varies
     | Show Table
    DownLoad: CSV

    Friction is an important factor which needs to be considered.However, in SHPB tests on concrete the effect of friction is negligibly small[22] in terms of its contribution to the total dynamic increase factor and is usually ignored in numerical simulations[23].Moreover, in SHPB tests measures are usually taken (i.e. by applying lubricant to both ends of concrete samples) to further reduce the effect of friction.Hence, in the present study the effect of friction is also ignored in the numerical simulations.

    In this study, the dynamic strength increase factor is employed to reflect the effect of specimen size in SHPB tests on concrete.If there is no size effect there will be no difference for concrete samples with different sizes or volumes in terms of the dynamic increase factors at the same strain rate.Otherwise there will be size effect in SHPB tests on concrete.

    Fig. 3 shows the numerical results for the dynamic increase factor due to inertia (confinement) effect Rc for concrete specimens with the same volume of different length/diameter ratios (i.e., Ø51 mm×51 mm, Ø64 mm×32 mm, Ø74 mm×24 mm) whilst Fig. 4 shows the numerical results for concrete samples with different volume of the same length/diameter ratios (i.e., Ø64 mm×32 mm, Ø80 mm×40 mm).It can be seen from Fig. 3 that the inertia effect is mainly related to the specimen volume and is insensitive to the length/diameter ratio when the volume and the material parameters of the specimens are kept the same.It also can be seen from Fig. 4 that the dynamic increase factor due to inertia (confinement) effect increases with the increase of the concrete specimen volume, which indicates that size effect does exist in SHPB tests on concrete as described above.

    Figure  3.  Comparison of Eq.(6) with numerical results fordynamic increase factor due to inertia (confinement)effect for concrete specimens with the same volumeof different length/diameter ratios
    Figure  4.  Comparison of Eq.(6) with numerical results fordynamic increase factor due to inertia (confinement)effect for concrete specimens with different volumeof the same length/diameter ratio

    On the basis of the numerical results discussed in the previous section, a new empirical equation which takes consideration of the specimen size effect in SHPB tests on concrete is suggested here to describe the dynamic increase factor due to inertia (confinement) effect Rc by using the monotonically increasing and continuous properties of exponential functions in a simple and easy to use form, namely

    Rc=Slg(˙ε˙ε0)+1βlg(VV0)W+1 (6)

    where S, W, β are the constants to be determined numerically, ˙ε is the strain rate, ˙ε0 is the reference strain rate for a reference specimen usually taken to be ˙ε0=1.0 s-1, V is the volume of a concrete specimen under investigation, V0 is the volume of a concrete sample with a reference size, say, Ø51 mm×51 mm.Set V=V0 in Eq.(6), one obtains the dynamic increase factor due to inertia effect for the reference concrete specimen Rc0.

    In the following, first, one determines the values of S and W in Eq.(6) using the numerical results for the concrete specimen with the same volume (i.e. V=V0) as presented in Fig. 3 and, then, the value of β in Eq.(6) using the numerical results for the concrete samples with different sizes/volumes as given in Fig. 4.

    Fig. 3 shows the comparison between the numerical results and Eq.(6) with S=6, W=2.8 and V=V0 as indicated by the dashed line.It can be seen from Fig. 3 that reasonable agreement is obtained.Further examination reveals that the value of S in Eq.(6) should not be a constant but a function of strain rate.The slope of the curve is too small at the low strain rate, and too large at the high strain rate.By using the centrosymmetric properties of the hyperbolic tangent and the x-axis paralleled, S is taken as the following form to increase the slope of the curve at the low strain rate and decrease the slope at the high strain rate, namely

    S=Fitanh[lg(˙ε˙ε0Wi)Si]+Gi (7)

    where Fi, Wi, Si, Gi are the constants determined using the numerical results for concrete samples with constant volume as presented in Fig. 3, namely, Fi=6.0, Wi=2.8, Si=0.8, Gi=8.5.The solid line in Fig. 3 is predicted from Eq.(6) with V=V0 together with Eq.(7).It is clear from Fig. 3 that good agreement is obtained.

    Fig. 4 shows the comparison between Eq.(6) with β=2.7 and the numerical results obtained for the concrete samples with different sizes/volumes.The values of all the other parameters in Eq.(6) are listed in Table 3.It is evident from Fig. 4 that the present model agrees well with the numerical simulations.

    Table  3.  Values of various parameters in Eq.(6) and Eq.(7)
    ˙ε0/s-1 Fi Si Wi Gi W β
    1.0 6.0 0.8 2.8 8.5 2.8 2.7
     | Show Table
    DownLoad: CSV

    In order to verify the validity of Eq.(6) more numerical simulations are performed on concrete specimens with different sizes. The numerical results are also presented in Fig. 4 and comparisons are also made between Eq.(6) and the numerical results.It can be seen from Fig. 4 that good agreement is obtained.

    Fig. 5 shows the normalization of all the numerical results for concrete specimens with different sizes as given in Fig. 3 and Fig. 4 with respect to those for the reference concrete sample, namely, Ø51 mm×51 mm.It is clear from Fig. 5 that all the numerical results collapse into one line.It lends further support to the validity of Eq.(6) for the dynamic increase factor due to inertia (confinement) effect which takes into consideration the influence of specimen size.

    Figure  5.  Variation of normalized numericallyobtained dynamic increase factor due to inertia(confinement) effect with strain rate

    The influence of specimen size in SHPB tests on concrete has been investigated numerically using a rate-independent material model.A new empirical equation for the dynamic increase factor due to inertia (confinement) effect has also been proposed which takes into account the specimen size effect through its volume.It is demonstrated that the new empirical formula agrees well with the numerical results for SHPB tests on concrete with different specimen sizes.

  • [1]
    张先锋, 赵晓宁. 多功能含能结构材料研究进展 [J]. 含能材料, 2009, 17(6): 731–739. doi: 10.3969/j.issn.1006-9941.2009.06.021

    ZHANG X F, ZHAO X N. Research progress of multifunctional energetic structural materials [J]. Chinese Journal of Energetic Materials, 2009, 17(6): 731–739. doi: 10.3969/j.issn.1006-9941.2009.06.021
    [2]
    HUGH E. Reactive fragment: US3961576 [P]. 1973–06–25.
    [3]
    刘晓俊, 任会兰. 一种反应材制备及准态力学特性研究 [J]. 北京理工大学学报, 2016, 36(4): 365–369.

    LIU X J, REN H L. Preparation of a reactive material and its quasi-state mechanical properties [J]. Journal of Beijing Institute of Technology, 2016, 36(4): 365–369.
    [4]
    WANG H X, LI Y C, FENG B. Compressive properties of PTFE/Al/Ni composite under uniaxial loading [J]. Journal of Materials Engineering and Performance, 2017, 26(5): 2331–2336. doi: 10.1007/s11665-017-2666-y
    [5]
    XU F Y, LIU S B, ZHENG Y F. Quasi-static compression properties and failure of PTFE/Al/W reactive materials [J]. Advanced Engineering Materials, 2016, 19(1): 1600350.
    [6]
    XU S, YANG S, ZHANG W. The mechanical behaviors of polytetrafluorethylene/Al/W energetic composites [J]. Journal of Physics Condensed Matter: AnInstitute of Physics Journal, 2009, 21(28): 285401. doi: 10.1088/0953-8984/21/28/285401
    [7]
    GE C, MAIMAITITUERSUN W, DONG Y. A study on the mechanical properties and impact-induced initiation characteristics of brittle PTFE/Al/W reactive materials [J]. Materials, 2017, 10(5): 452. doi: 10.3390/ma10050452
    [8]
    刘晓俊, 任会兰, 宁建国. Zr-W多功能含能结构材料的制备及动态压缩特性 [J]. 复合材料学报, 2016, 33(10): 2297–2302.

    LIU X J, REN H L, NING J G. Preparation and dynamic compression properties of Zr-W multifunctional energetic structural materials [J]. Journal of Composite Materials, 2016, 33(10): 2297–2302.
    [9]
    王占磊. 爆炸压实 W-Cu 纳米合金及其聚能破甲应用研究 [D]. 大连: 大连理工大学, 2012.

    WANG Z L. Explosive compaction W-Cu nano-alloy and its application of energy-absorbing armor [D]. Dalian: Dalian University of Technology, 2012.
    [10]
    MAMALIS A G, VOTTEA I N, MANOLAKOS D E. On the modelling of the compaction mechanism of shock compacted powders [J]. Journal of Materials Processing Technology, 2001, 108(2): 165–178. doi: 10.1016/S0924-0136(00)00748-2
    [11]
    MORRIS D G. Bonding processes during the dynamic compaction of metallic powders [J]. Materials Science & Engineering, 1983, 57(2): 187–195.
    [12]
    FARINHA A R, MENDES R, BARANDA J, et al. Behavior of explosive compacted/consolidated of nanometric copper powders [J]. Journal of Alloys & Compounds, 2009, 483(1–2): 235–238.
    [13]
    EAKINS D, THADHANI N N. Shock-induced reaction in a flake nickel + spherical aluminum powder mixture [J]. Journal of Applied Physics, 2006, 100(11): 113521. doi: 10.1063/1.2396797
    [14]
    EAKINS D, THADHANI N N. Discrete particle simulation of shock wave propagation in a binary Ni+Al powder mixture [J]. Journal of Applied Physics, 2007, 101(4): 3635.
    [15]
    EAKINS D E, THADHANI N N. The shock-densification behavior of three distinct Ni+Al powder mixtures [J]. Applied Physics Letters, 2008, 92(11): 111903. doi: 10.1063/1.2896653
    [16]
    EAKINS D E, THADHANI N N. Mesoscale simulation of the configuration-dependent shock-compression response of Ni+Al powder mixtures [J]. Acta Materialia, 2008, 56(7): 1496–1510.
    [17]
    YANG R Y, YU A B, CHOI S K. Agglomeration of fine particles subjected to centripetal compaction [J]. Powder Technology, 2008, 184(1): 122–129. doi: 10.1016/j.powtec.2007.08.010
    [18]
    金属材料室温压缩试验方法: GB/T 7314–2005 [S]. 北京: 国家质量监督检验检疫总局, 2005.
    [19]
    张庆明, 刘彦, 黄风雷.材料的动力学行为 [M]. 北京: 国防工业出版社, 2006: 103-182.

    ZHANG Q M, LIU Y, HUANG F L.Dynamic behavior of materials [M]. Beijing: National Defence Industry Press, 2006: 103–182.
    [20]
    WEI C T, VITALI E, JIANG F. Quasi-static and dynamic response of explosively consolidated metal-aluminum powder mixtures [J]. Acta Materialia, 2012, 60(3): 1418–1432.
  • Relative Articles

    [1]GAO Linyu, DU Shiyu, CHANG Hui, ZHANG Tuanwei, WANG Zhihua. Strain-rate and temperature dependent compressive deformation behavior of CrCoNiSi0.3 medium-entropy alloy[J]. Chinese Journal of High Pressure Physics. doi: 10.11858/gywlxb.20251047
    [2]SUN Yan-Yun, LIU Fu-Sheng, ZHANG Ming-Jian, XU Li-Hua. Three-Body Interactions and Shock Compression Properties of Condensed Nitrogen[J]. Chinese Journal of High Pressure Physics, 2009, 23(2): 137-142 . doi: 10.11858/gywlxb.2009.02.010
    [3]SUN Yu-Chao. Research on Thermodynamic Mechanism of Diamond Surface Metallization[J]. Chinese Journal of High Pressure Physics, 1992, 6(1): 37-47 . doi: 10.11858/gywlxb.1992.01.006
    [4]DONG Lian-Ke, WANG Xiao-Wei, WANG Ke-Gang, Lü Guo-Hao. Uncertain Problem in Application of Fractal to Fracture Toughness[J]. Chinese Journal of High Pressure Physics, 1990, 4(2): 118-122 . doi: 10.11858/gywlxb.1990.02.007
    [5]ZHOU Kai-Yong, YU Xin-Lu. A Material Test Technique for Super-High Pressure Gasket Material[J]. Chinese Journal of High Pressure Physics, 1990, 4(1): 7-16 . doi: 10.11858/gywlxb.1990.01.002
    [6]YAO Bin, ZHANG Qiang, SU Wen-Hui, XU Da-Peng. A Study on the Formation and Stability of Quasicrystal Al4Mn under High Static Pressure[J]. Chinese Journal of High Pressure Physics, 1990, 4(1): 50-56 . doi: 10.11858/gywlxb.1990.01.008
    [7]LONG Qi-Wei. Perimeter-Area Relation of Fractal Island[J]. Chinese Journal of High Pressure Physics, 1990, 4(4): 259-262 . doi: 10.11858/gywlxb.1990.04.004
    [8]ZHAO Yong-Nian, ZOU Guang-Tian. High-Pressure Raman Spectra and Structure Phase Transition for Molecule Crystals[J]. Chinese Journal of High Pressure Physics, 1989, 3(4): 269-278 . doi: 10.11858/gywlxb.1989.04.002
    [9]LI Qiang-Min, SU Wen-Hui, LONG Xiang, WU Dai-Ming, GAO Zhong-Min. High-Pressure and Temperature Syntheses and Structural Stability Studies of Some Rare-Earth Oxides[J]. Chinese Journal of High Pressure Physics, 1989, 3(1): 42-50 . doi: 10.11858/gywlxb.1989.01.006
    [10]ZHANG Chun-Bin. Molecular Dynamics Study on Equation of State[J]. Chinese Journal of High Pressure Physics, 1989, 3(1): 51-57 . doi: 10.11858/gywlxb.1989.01.007
    [11]YU Wan-Rui. Molecular Dynamics Studies of Material Behavior at High Strain Rates[J]. Chinese Journal of High Pressure Physics, 1989, 3(2): 143-147 . doi: 10.11858/gywlxb.1989.02.006
    [12]HU Xiao-Mian. Crystal structure Stability Study by Molecular Dynamics Method with Variable Cell[J]. Chinese Journal of High Pressure Physics, 1989, 3(2): 132-142 . doi: 10.11858/gywlxb.1989.02.005
    [13]WANG Yi-Feng, SU Wen-Hui, QIAN Zheng-Nan, MA Xian-Feng, YAN Xue-Wei. A Study on the High Temperature-High Pressure Stability and the X-Ray Phase Analysis for Compounds R2Fe4/3W2/3O7 with Pyrochlore Structure[J]. Chinese Journal of High Pressure Physics, 1988, 2(4): 296-304 . doi: 10.11858/gywlxb.1988.04.002
    [14]ZHANG Qiang, SU Wen-Hui. A Study on the Formation and the Stability of Al6Mn Quasicrystal under High Static Pressure[J]. Chinese Journal of High Pressure Physics, 1988, 2(1): 58-66 . doi: 10.11858/gywlxb.1988.01.008
    [15]YU Wan-Rui, LIU Ge-San. Molecular Dynamic Investigation of Shock Waves in the Solid[J]. Chinese Journal of High Pressure Physics, 1988, 2(1): 73-78 . doi: 10.11858/gywlxb.1988.01.010
    [16]YU Wan-Rui, ZHANG Ruo-Qi, LIU Ge-San. Monte-Carlo Investigation of Interaction of Relativistic Electron-Beam with Solid and Calculation of Its Mechanical Effects[J]. Chinese Journal of High Pressure Physics, 1987, 1(1): 58-70 . doi: 10.11858/gywlxb.1987.01.008
    [17]CUI Shuo-Jing, XU Da-Peng, SU Wen-Hui, WANG Fu-Song, WANG Sheng-Long. An Investigation on Electric Properties of the Polyaniline Complexes under High Pressure[J]. Chinese Journal of High Pressure Physics, 1987, 1(1): 71-75 . doi: 10.11858/gywlxb.1987.01.009
    [18]SHEN Zhong-Yi, CHEN Gui-Yu, ZHANG Yun, YIN Xiu-Jun, SHEN De-Jiu, WU Hao-Quan. The Thermal Stability and Variation of Precipitated Phase in Amorphous Zr70Cu30 Alloy at High Pressure[J]. Chinese Journal of High Pressure Physics, 1987, 1(2): 165-172 . doi: 10.11858/gywlxb.1987.02.010
    [19]ZHANG Chun-Bin, ZHANG Shi-Ze, LI Shao-Meng, LIN Shao-Ming, ZHANG Wan-Xiang. EOS and Phase Structure Stability of Lithium[J]. Chinese Journal of High Pressure Physics, 1987, 1(2): 150-154 . doi: 10.11858/gywlxb.1987.02.008
    [20]CHEN Dong-Quan, XIE Guo-Qiang. Molecular Dynamics Simulation of Polymorphous Transitions[J]. Chinese Journal of High Pressure Physics, 1987, 1(1): 50-57 . doi: 10.11858/gywlxb.1987.01.007
  • 加载中

Catalog

    通讯作者: 陈斌, bchen63@163.com
    • 1. 

      沈阳化工大学材料科学与工程学院 沈阳 110142

    1. 本站搜索
    2. 百度学术搜索
    3. 万方数据库搜索
    4. CNKI搜索

    Figures(11)  / Tables(1)

    Article Metrics

    Article views(8979) PDF downloads(47) Cited by()
    Proportional views
    Related

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return