Next Article in Journal
Fabrication of Porous Poly(3-hydroxybutyrate-co-3-hydroxyhexanoate) Monoliths via Thermally Induced Phase Separation
Next Article in Special Issue
Decorating Nanoparticle Surface for Targeted Drug Delivery: Opportunities and Challenges
Previous Article in Journal
Announcement of the 2016 Polymers Young Investigator Award
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Experimentation and Modeling of the Tension Behavior of Polycarbonate at High Strain Rates

1
Engineering Simulation and Aerospace Computing (ESAC), Northwestern Polytechnical University, Xi’an 710072, China
2
Mechanics Unit (UME), ENSTA-ParisTech, Paris, 91761 Palaiseau, France
*
Authors to whom correspondence should be addressed.
Polymers 2016, 8(3), 63; https://doi.org/10.3390/polym8030063
Submission received: 28 January 2016 / Revised: 23 February 2016 / Accepted: 24 February 2016 / Published: 29 February 2016
(This article belongs to the Special Issue Computational Modeling and Simulation in Polymer)

Abstract

:
A comprehensive understanding of the mechanical behavior of polycarbonate (PC) under high-rate loadings is essential for better design of PC products. In this work, the mechanical behavior of PC is studied during tensile loading at high strain rates, using a split Hopkinson tension bar (SHTB). A modified experimental technique based on the SHTB is proposed to perform the tension testing on PC at rates exceeding 1000 s−1. The effect of strain rates on the tension stress–strain law of PC is investigated over a wide range of strain rates (0.0005–4500 s−1). Based on the experiments, a physically based constitutive model is developed to describe the strain rate dependent tensile stress–strain law. The high rate tensile deformation mechanics of PC are further studied via finite element simulations using the LSDYNA code together with the developed constitutive model.

Graphical Abstract

1. Introduction

Polycarbonate (PC) is a thermoplastic polymeric material with high transparency, high ductility, impact resistance and is comparatively lightweight. It has been widely used in impact protection applications including aircraft canopies, face shields, goggles, windshields and windows and blast shields [1,2]. Thus, impact response of the PC products is a subject of critical interest. High strain rates are commonly encountered during an impact event such as projectile striking or blast loading. Therefore, an accurate understanding and modeling of the mechanical behavior of PC at high strain rates is of great importance.
The mechanical behavior of materials at high strain rates has been extensively studied using various experimental procedures ranging from impact (gas gun impact, Taylor impact, and Izod impact) tests [3] to split Hopkinson pressure bar (SHPB) test [4]. Amongst these procedures, SHPB test has been instrumental in obtaining stress–strain behavior of materials at high strain rates up to 103 s−1.
Recently, experimental studies with SHPB have been performed to investigate the compressive behavior of PC at high strain rates [5,6,7,8,9,10,11]. Compared with metals, high-rate experiment of PC in compression with the SHPB is more complicated due to its low density, low modulus and low yield stress. The low densities and low elastic wave velocities ultimately result in low impedances, which lead to low amplitude of transmitted pulses and consequently increase the signal-to-noise ratio. To address these problems, some modifications such as the pulse shaping technique and the use of low-impendence aluminum, titanium or polymeric bars are proposed to perform the SHPB test [12]. Experimental results indicate that PC exhibits an obvious elastic-plastic stress–strain behavior, including linear elasticity, nonlinear elasticity, yielding, post-yield strain softening and strain hardening when subjected to high strain rates compressive loading. Moreover, it can be observed that the yield stress increases dramatically as the strain rate increases. Mulliken and Boyce [8] investigated the relationship between the yield stress and logarithm strain rate in a wide strain-rate region and indicated that there exists a transition threshold of rate sensitivity for yield behavior of PC.
Owing to the fact that there exist significant differences in tensile and compressive deformation behavior of polymers under high strain rate loadings [13], it is therefore necessary to study the dynamic tension behavior of PC. Moreover, tension tests provide more understandings on damage, failure and fracture behavior of materials. However, due to the experimental difficulties in tensile test with SHPB, very few studies have been reported on the tension behavior for polymers [14,15,16,17,18], especially for PC [17,18] at high strain rates. Design change of the SHPB into a split Hopkinson tension bar (SHTB) enabling tensile loading conditions is thus required for the high strain rate tensile test. Sarva and Boyce [17] used a split-collar type SHTB to investigate the behavior of PC during tensile loading at high strain rates. In their study, the threaded PC sample was attached between two bars and a split collar surrounding the sample was sandwiched between two bars. By using the split collar, the initial compressive pulse was transmitted into the second bar without loading the sample. The subsequent reflected tensile pulse (from the free rear end of the second bar) then loaded the sample. By using this split-collar type SHTB, the effects of varying strain rate, overall imposed strain magnitude and specimen geometry on the mechanical response were examined in details. Cao et al. [18] implemented high strain rate tension tests of PC using a modified SHTB. The PC sample was bonded to incident and transmitted bars by means of adhesion. The incident stress pulse was initiated by the impact of the hammer fixed on the high-speed rotating disk on the impact block, which caused the prefixed metal bar connected to the block and the incident bar to deform. The deformation of the prefixed metal bar generated the tensile loading stress pulse.
A constitutive model that accurately represents the material response of PC can provide convenient and useful guidelines on product design and effectively decrease the experimental cost. Significant advances have been achieved in developing phenomenological models and physically based models to describe the large strain, temperature, and rate-dependent deformation behavior of PC and other polymeric materials. Mulliken and Boyce [8] proposed a physically based high-rate constitutive model for the three-dimensional deformation of PC and Poly (methyl methacrylate) (PMMA) at strain rates ranging from 10−4 to 104 s−1. Uniaxial compression tests were used to characterize the rate-dependent yield and post-yield behavior of PC and PMMA. The proposed constitutive model was shown to correctly predict yield stress values, as well as the strain rate regime of the transition in the yield behavior. Sarva and Boyce [17] then extended this constitutive model to encompass high-rate tensile behavior of PC. Anand et al. [19] conducted large-strain compression experiments on PC and PMMA and presented a thermo-mechanically coupled elasto-viscoplasticity model to describe the strain rate and temperature dependent large-deformation response. Based on the high-rate tensile experiment investigation, Cao et al. [18] developed a physically based three-dimensional elastic-plastic constitutive model to characterize the rate-temperature dependent yield and post-yield behavior of PC when subjected to tension loading. Bouvard et al. [20] presented a constitutive model for amorphous polymers using a thermodynamic approach with physically motivated internal state variables. Model parameters were determined for PC, PMMA and polystyrene using high-rate compressive test results. Based on an elastic-viscoplastic rheological approach, Richeton et al. [21] developed a three-dimensional constitutive model to describe the material response of PC and PMMA for a large range of strain rates. The model parameters are fitted for each polymer on the experimental data of high-rate uniaxial compression tests. The aforementioned models have shown their abilities to capture the mechanical response of PC within a large strain rate range. However, these models are primarily calibrated using the results of high strain rate compressive data. Thus, most of the models can be only used to predict the high strain rate compressive behavior of PC.
From the literature review it can be concluded that characterizing the tensile behavior of PC at high strain rates is still insufficient at present, in both experimentation and modeling. Most of existed studies of the high rate deformation of PC concern the compressive behaviors. Due to the experimental difficulties in tensile test with SHPB, very few studies have been reported on the high rate (exceeding 1000 s−1) tensile behavior for PC. This is why we conduct a detailed investigation of the tensile behavior of PC at high strain rates in the present paper. A modified experimental technique based on the SHTB is proposed to perform the tension testing on PC at rates exceeding 1000 s−1. The chronological progression of the dynamic deformation is captured with a high-speed CCD camera. Quasi-static tensile tests are conducted as well to investigate the effects of varying strain rate on the mechanical response of PC. Based on the experimental investigation, a rate dependent constitutive model is developed to describe the strain-rate dependent tension stress–strain law. The high rate tensile deformation mechanics of PC are further studied via finite element simulations using the LSDYNA code together with the developed constitutive model.

2. Experiment

2.1. Material and Specimen

The material used in the present study is PC Lexan 141R. The mass density of PC is 1.2 × 103 kg/m3. The glass transition temperature (Tg) of PC determined by the Dynamic Mechanical Thermal Analysis (DMTA) is 148 °C. Geometry and dimension of the specimen used for both SHTB and quasi-static tensile tests are shown in Figure 1. Notice that tight-handed and left-handed threads are provided on either side of these specimens for easy mounting. All the specimens are carefully manufactured from 12 mm thick plates of PC and kept at room temperature for more than three days prior to testing.

2.2. High Strain-Rate Uniaxial Tension Test

High strain-rate tension tests at 103 s−1 strain rate are carried out using the SHTB shown in Figure 2. The SHTB setup consists of an air gun system, a striker tube, an incident bar with a transfer flange at one end, a transmission bar, a momentum trap bar and the data acquisition system. The striker tube and all the bars are made of 18Ni steel and behave elastically during the test. The striker tube is 380 mm in length, with an outer diameter of 26 mm and inner diameter of 19 mm. The incident and transmission tension bars are 2800 and 1200 mm in length, respectively; both are 19 mm in diameter. The momentum trap bar is 430 mm in length and 19 mm in diameter. Both the incident bar and transmission bar contain threaded holes to connect the specimen. The air gun system is used to propel the striker tube. A detailed illustration of the working principle of the air gun system can be found in [22].
The striker tube accelerated by the air gun slides along the incident bar and impacts the flange to generate a tensile pulse in the incident bar. The duration of the incident pulse can be controlled by adjusting the length of the striker tube. As shown in Figure 2, a gap is precisely preset to separate the momentum trap bar and the transfer flange of the incident bar. This gap is set such that the end of the momentum trap bar and the face of the transfer flange are brought in contact, once the tensile pulse generated by the striker tube is completely transferred into the incident bar through the transfer flange. The generated tensile pulse travels down the incident bar and then propagates into the specimen, where it is partly transmitted into the transmission bar, and is partly reflected as a compressive pulse back into the incident bar. This reflected compressive pulse is then transmitted into the momentum trap bar, and reflects off the free end of this bar as a tensile pulse. Since the contact interface with the transfer flange cannot support tension, this tensile pulse is trapped in the momentum trap bar. Therefore, the reloading of the specimen by the reflected pulse is effectively prevented.
The theoretical value of the separation between the momentum trap bar and the transfer flange, us, can be estimated by the following expression:
u s = 2 C 0 0 t ε c ( t ) d t   
where C0 is the longitudinal elastic wave velocity in the incident bar. ε c ( t ) is the compression strain in the flange. In practice, the calculated us will be further optimized by a few trials.
The incident strain ε i ( t ) , reflected strain ε r ( t ) and transmitted strain ε t ( t ) are recorded as functions of time t using strain gages G1 and G2 attached to the bars at two locations. The axial forces and displacements on the front and rear ends of the specimen can be calculated from ε i ( t ) , ε r ( t ) and ε t ( t ) through the following expressions:
F 1 ( t ) = E A ( ε i ( t ) + ε r ( t ) )        F 2 ( t ) = E A ε t ( t )
u 1 ( t ) = C 0 0 t ( ε i ( t ) ε r ( t ) ) d t        u 2 ( t ) = C 0 0 t ε t ( t ) d t   
where F1(t) and F2(t) are the forces and u1(t) and u2(t) are the displacements on the front end (right end of the incident bar) and rear end (left end of the transmitted bar) of the specimen, respectively. E and A are Young’s modulus and cross-sectional area of the bars. C0 is the longitudinal wave velocity in the bars.
The engineering stress σ s ( t ) , strain ε s ( t ) and strain rate ε ˙ s ( t ) in the specimen can further be obtained as:
σ s ( t ) = F 1 ( t ) + F 2 ( t ) 2 A s = E A 2 A s ( ε i ( t ) + ε r ( t ) + ε t ( t ) )   
ε s ( t ) = u 1 ( t ) u 2 ( t ) l s = C 0 l s 0 t ( ε i ( t ) ε r ( t ) ε t ( t ) ) d t
ε ˙ s ( t ) = C 0 l s ( ε i ( t ) ε r ( t ) ε t ( t ) )
where As and ls are the cross-sectional area and gage length of the specimen, respectively.
The principle of the SHTB test assumes that there exists a state of stress equilibrium and uniform deformation in the specimen during the process of impact loading. Therefore, as stress equilibrium is achieved during the test. The forces on the two ends of specimen are equal, i.e., F1(t) = F2(t). The relation ε i ( t ) + ε r ( t ) = ε t ( t ) can be consequently deduced and the engineering stress σ s ( t ) , strain ε s ( t ) and strain rate ε ˙ s ( t ) in the specimen can further be given as:
σ s ( t )    = E A A s ε t ( t )
ε s ( t ) = 2 C 0 l s 0 t ε r ( t ) d t
ε ˙ s ( t ) = 2 C 0 l s ε r ( t )
Then, the true stress and true strain in the specimen can be easily calculated from the engineering stress and engineering strain using the volume constancy law.

2.3. High-Speed Photography

A Photron SAX2 high-speed CCD camera, capable of acquiring images at a frame rate of one million frames per second, is used to photograph the dynamic deformation of the tensile specimens. For better illumination, high performance strobes are placed behind the specimen for silhouette lighting of the specimen. The camera is programmed to record a sequence of 30 separate images at prescribed time intervals, and images are acquired from a point of view normal to the sample.

2.4. Quasi-Static Uniaxial Tension Tests

Quasi-static tension tests at strain rates of 0.5 × 10−3, 1.0 × 10−3, 1.0 × 10−2 and 1.0 × 10−1 s−1 are performed on a DNS-100 electronic universal testing system. The specimen geometry used in quasi-static tension tests is exactly the same as that used in the high-rate tests. A fixture device including two threaded grip made of heat-treated high strength steel is used to connect the specimen, as shown in Figure 3.

3. Experimental Results

The typical strain gage signals recorded on the incident and transmitted bars are displayed in Figure 4. Obviously, there exists a nearly flat plateau part in the reflected pulse, which means a steady strain rate condition. In addition, the velocity profiles of the incident bar–specimen interface and transmission bar–specimen interface during a test are shown in Figure 5.
Dynamic tensile tests are performed at strain rates of 1400, 2000, 3500 and 4500 s−1. Figure 6 show different strain rates of specimen during tests, we can see the strain rates are stable enough to guarantee the accuracy of the experiments.
Figure 7 presents the tensile true stress–true strain responses of PC for various strain rates ranging from quasi-static to dynamic loadings. It can be observed for all strain rates that a stress drop occurs following a plastic flow platform when the stress reaches the peak point. The response of PC is strongly dependent on the strain rate. The dynamic response of PC is significantly distinctive from the quasi-static one. In the case of quasi-static loading, the PC specimens are deformed in a ductile manner up to a strain of 0.60. The response includes linear elastic and nonlinear elastic homogeneous deformations and the strain softening and hardening. Furthermore, the yield stress and flow stress increase slightly with the increase of strain rate. In the case of high strain rate loading, the values of yield stress and strain at yield at high strain rates increase apparently than those under quasi-static loading, but they change slightly at different high strain rates investigated in the present paper. Compared with Figure 7a, there are more fluctuations in the dynamic tensile stress–strain curves, which is induced by the imperceptible gaps exist in the threaded connection of bars and specimen. This phenomenon is in agreement with the work of Boyce et al. [17]. Here, it should be mentioned that the hardening behaviors are not captured by the loading cases at strain rates of 1400 and 2000 s−1 (as shown in Figure 7). This is due to the limit of the split Hopkinson bar measurement in our experimental set-up. In our experiment, the duration of the incident stress pulse is restricted by the length of the prefixed steel bar. Since the duration of the incident stress pulse is constant in all the high rate tests, the total measurable uniaxial strain is decreased with the reduction of strain rate. Therefore, for strain rates of 1400 and 2000 s−1, only the softening behaviors are recorded. The hardening behavior after softening cannot be captured.
For a more comprehensive understanding of the influence of strain rate on the tension properties of PC, the true yield stress plotted against the logarithm of strain rate is shown in Figure 8. It is found that the yield stress of PC is quite sensitive to the strain rate. The yield stress increases linearly with the logarithm of strain rate at both quasi-static regime and high strain rate regime. However, the values of yield stress at dynamic loading conditions are much higher compared with quasi-static results. The experimental results indicate a bilinear relationship between the tensile yield stress and the logarithm of strain rate: at the low strain rates, PC displays weak strain rate dependence, while at high strain rates, the yield stress increases dramatically with increasing strain rate. This phenomenon is similar with the strain rate dependence of compressive yield stress of PC reported by Richetona et al. [9]. In addition, Richetona et al. [9] indicated that the increase of the compressive yield stress is correlated to secondary molecular processes. An increasing strain rate would decrease the molecular mobility of the polymer chains by making the chains stiffer. This can be also used to explain the strain rate dependence of tensile yield stress of PC.
Note that since the tension tests at strain rate range of 0.1–1000 s−1 (moderate rate) are not performed in this study, there is a gap in data about strain rate passing from a logarithm of strain rate of −1 to 3 (see Figure 8). A moderate strain-rate testing set-up is further needed for these tests. At present, we only carried out the quasi-static (below 0.1 s−1) and high-rate (exceeding 1000 s−1) tension tests. However, the bilinear relationship between the yield stress and the strain rate indicated by our experiment is similar with the strain rate dependence of compressive yield stress of PC reported in the literature. Of course, more experimental data could generate a more precise description of the strain rate dependence of tensile yield stress. In our future study, the moderate strain-rate tests will be employed to investigate the tension behavior of PC in detail.
Figure 9 shows the strain softening behavior of PC from static loading to dynamic loading. Strain softening is an intrinsic mechanical characteristic of amorphous polymers, which demonstrates a dramatic stress drop after yielding. This phenomenon is considered as a degradation of material caused by the rupture of molecular chains network. The upper yield point and lower yield point represent the yield strength and the minimal value after yield, respectively. The yield drop is defined as the difference between upper yield stress and lower yield stress. Although both upper yield stress and lower yield stress increase in bilinear relationship (as shown in Figure 9) with the logarithmic strain rate, the yield drop remains constant, indicating that strain softening is independent on strain rate. It should be noticed that due to the limitation of the measurable uniaxial strains for the tests at strain rates from 1400 to 2000 s−1, the lower yield points are not given in Figure 9 for these loading cases.

4. Constitutive Model

A constitutive model is essential for the modeling and prediction of the deformation and failure of structures and materials. In this study, the experimental tension stress–strain responses at various strain rates show a nonlinear deformation behavior including elasticity, yield, strain softening and strain hardening. In addition, the experimental investigations demonstrate bilinear relationship between yield stress and the logarithm of strain rate. The similar deformation behavior can be also found for PC under compression loading conditions. In our previous study [23], within the framework of the thermodynamics of irreversible processes, a physically based constitutive model for the high rate behavior of PC was developed based on the uniaxial compression tests at high strain rates. The high rate compressive behavior of PC, including elasticity, yield, strain softening, strain hardening and yield-rate behavior, has been successfully captured by the constitutive model. In this study, this model is extended to describe the strain rate dependent tension response of PC.

4.1. Strain Rate Dependent Yield

To model the strain rate dependent yield behavior over the whole strain rate range, a power exponent function is employed here to approximate the bilinear relationship:
σ y = K ( 1 + C r ε ¯ ˙ e q m )
where K is the yield stress under uniaxial static loading (minimum strain rate in this study). Cr and m are strain rate relative constants. ε ¯ ˙ e q is the equivalent strain rate and defined as:
ε ¯ ˙ e q = 2 3 d ε i j d ε i j
where d ε i j is the deviator strain increment.

4.2. Strain Softening and Hardening

The plastic flow of PC consists of strain softening and strain hardening. These two phenomena are caused by different micro-mechanisms. Strain softening is mainly due to the fracture of molecular main chains and the disentanglement of branches, which can be regarded as the damage of molecular network [24]. Meanwhile, strain hardening results from the orientation of main chains along the applied stress direction. In fact, the orientation of main chains and the disentanglement of branches evolve simultaneously. Hence, it is believed that this competitive mechanism between hardening and softening determines the plastic flow of amorphous polymers. In this study, a hardening function combined with the damage evolution is proposed to describe the plastic flow rule of PC.
The hardening function of the virgin material (without damage) can be approximated by a concise function:
σ e f f = C h ( ε ¯ p ) γ
where Ch and γ are hardening constants identified by experimental results. σ e f f is the effective stress of the virgin material. ε ¯ p is accumulated plastic strain given by:
ε ¯ p = 0 t 2 3 d ε i j p d ε i j p d t
where d ε i j p is the plastic strain increment.
For the purpose of measuring damage macroscopically, Lemaitre [25] proposed an important hypothesis indicating that the deformation of damaged material can be modeled in the form of the constitutive relation of virgin material. In a damaged material, the equivalent stress σ ¯ e q can be expressed as a function of damage variable D and effective stress σ e f f of the virgin material:
σ ¯ e q = ( 1 D ) σ e f f
The damage variable D evolves according to the following equation [25]:
D = C d ( ε ¯ p ) χ
where Cd and χ are damage evolution parameters. According to [25], the evolution of damage can be reflected by the decrease of elastic modulus at macroscopic scale, which is easily obtained by cyclic loading/unloading uniaxial test.
Combining Equations (12)–(15), the equivalent stress σ ¯ of a damaged material can be expressed as:
σ ¯ e q = ( 1 C d ( 0 t 2 3 d ε i j p d ε i j p d t ) χ ) ( C h ( 0 t 2 3 d ε i j p d ε i j p d t ) γ )
In this study, Equation (16) is used to describe the plastic flow rule of PC. The first term of this equation describes the strain softening behavior of PC. The second term reflects the strain hardening.

5. Finite Element Simulation of Tension Responses

5.1. Finite Element Models

The constitutive model is programmed into the explicit finite element program LS-DYNA by employing a user-defined material subroutine (UMAT). Numerical simulations are conducted to study the deformation behavior of PC specimens during quasi-static and dynamic tensile loading.
The material constants and model parameters are determined by experimental measurement and curve fitting. They are listed as follows.
  • Material constants
    • Elastic modulus: E = 2271 MPa
    • Poisson’s ratio: 0.4
    • Density: ρ = 1190 kg/m3
  • Model parameters
The model parameters used in the finite element simulation can be divided into three categories: yield parameters, hardening parameters and damage parameters. The identified values of these parameters are
  • Yield parameters: K = 62.04 MPa, Cr = 0.307, m = 0.108
  • Hardening parameters: Ch = 288.1 MPa, γ = 1.994
  • Damage parameters: Cd = 0.365, χ = 0.399
The specimen geometry used in the numerical simulations is identical to that in the tensile tests. The specimen is meshed by eight-node brick element, as shown in Figure 10. Since the deformation is mostly concentrated in the gauge section, a finer mesh is chosen for this region. For simplicity, the Hopkinson bars are not included in simulation.
The simulations are performed at constant engineering strain rates corresponding to the actual tension tests. One end of the connection part of the specimen is fixed and the velocity boundary conditions are applied to the other end of the specimen connection part. To keep the consistence with the actual experimental measurement, the engineering stress is defined as the tensile reaction forces divided by the original cross-sectional area in the middle of the specimen. The engineering strain is obtained from the relative displacements between the two ends of the specimen gage section. Finally, the true stress and true strain are calculated from the engineering stress and engineering strain through the volume constancy law.

5.2. Numerical Results and Discussions

Firstly, the comparisons of the numerical results of true stress–true strain responses with the experimental data for the uniaxial tension at different strain rates are presented in Figure 11. It can be found that the presented model is capable of capturing the typical behaviors of the tension deformation of PC, including the elastic stage, the yield peak, the post-yield softening and hardening. In addition, the model has proved its ability to describe the rate dependent behavior of PC under tension loading within a wide range of strain rates.
Figure 12 shows a comparison of the experimentally observed true yield stress values of PC and the corresponding model predictions. The predicted results coincide well with the experimental data. The yield stress increases with the logarithm of strain rate and shows increased rate sensitivity at high rates.
To clearly understand the mechanics of tension deformation of PC at high strain rate, the numerical results and the high-speed photographs of the gauge region of specimen at various stages of deformation are listed in Table 1. The contours of axial stress at different instants during the deformation at strain rate of 4500 s−1 are displayed. The time intervals are carefully chosen to enable a direct comparison with the photographs. It can be observed that to some extent the simulation replicates the deformation profile observed in the photographs. As shown at 20 μs, the displacement initiates from the incident end, resulting in stress fields emanating from that end. After the initial stress wave reverberations, the stress is found to be relatively uniform in the gauge region as shown at 55 μs, which corresponds to elastic deformation. As observed in the photographs from 20 to 55 μs, the gauge uniformly elongates. Yielding then occurs at a stress level of nearly 115 MPa at 65 μs. A neck is found to initiate at the middle region of gauge in the photograph of 65 μs. The stress levels are then higher in the neck region due to the reduced cross-sectional area.

6. Conclusions

The tension behaviors of PC at high strain rates up to 4500 s1 are investigated in this paper. A comprehensive experimental and finite element study of PC is conducted to investigate the mechanics of deformations within a wide range of strain rates of high-rate tension. The high strain rate tension tests are carried out using a modified SHTB. The experimental tension stress–strain responses at various strain rates show a nonlinear deformation behavior including elasticity, yield, strain softening and strain hardening. In addition, the experimental results indicate that the tension behavior of PC is sensitive to strain rate and the values of yield stress increase dramatically with the increase of strain rate. Based on the experimental observations, a physically based constitutive model is developed and combined with finite element simulations to describe the tension response of PC. The numerical results show that the developed model is capable of capturing the typical behaviors of the tension deformation of PC and describing the rate dependent behavior of PC under tension loading within a wide range of strain rates. Moreover, the numerical results of deformation evolutions at strain rate of 4500 s1 are compared with the high-speed photographs. It is shown that the simulation replicates the deformation profile observed in the photographs. The computed axial stresses and the high-speed photographs are used to discuss the stress and deformation evolution of PC during a high rate tension.

Acknowledgments

This work is supported by 973 Program (2012CB025904), the Fundamental Research Funds for the Central Universities (3102014JCS05003) and National Natural Science Foundation of China (11302174, 11432011).

Author Contributions

Yingjie Xu and Weihong Zhang proposed the concept of the paper and developed the model. Tenglong Gao performed experiments and finite element computations. Jun Wang and Yingjie Xu derived the constitutive model. Yingjie Xu and Weihong Zhang analyzed the results and wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Xu, Y.J.; Zhang, Q.W.; Zhang, W.H.; Zhang, P. Optimization of injection molding process parameters to improve the mechanical performance of polymer product against impact. Int. J. Adv. Manuf. Technol. 2015, 76, 2199–2208. [Google Scholar] [CrossRef]
  2. Wang, J.; Xu, Y.J.; Zhang, W.H. Finite element simulation of PMMA aircraft windshield against bird strike by using a rate and temperature dependent viscoelastic constitutive model. Compos. Struct. 2014, 108, 21–30. [Google Scholar] [CrossRef]
  3. Xu, Y.J.; Lu, H.; Gao, T.L.; Zhang, W.H. Predicting the low-velocity impact behavior of polycarbonate: Influence of thermal history during injection molding. Int. J. Impact. Eng. 2015, 86, 265–273. [Google Scholar] [CrossRef]
  4. Gray, G.T. Classic split-Hopkinson pressure bar testing. In Mechanical Testing and Evaluation, Metals Handbook; American Society for Metals: Materials Park, OH, USA, 2008; Volume 8. [Google Scholar]
  5. Rietsch, F.; Bouette, B. The compression yield behaviour of polycarbonate over a wide range of strain rates and temperatures. Eur. Polym. J. 1990, 26, 1071–1075. [Google Scholar] [CrossRef]
  6. Duan, Y.; Saigal, A.; Greif, R.; Zimmerman, M.A. A uniform phenomenological constitutive model for glassy and semicrystalline polymers. Polym. Eng. Sci. 2001, 41, 1322–1328. [Google Scholar] [CrossRef]
  7. Siviour, C.R.; Walley, S.M.; Proud, W.G.; Field, J.E. The high strain rate compressive behaviour of polycarbonate and polyvinylidene difluoride. Polymer 2005, 46, 12546–12555. [Google Scholar] [CrossRef]
  8. Mulliken, A.D.; Boyce, M.C. Mechanics of the rate-dependent elastic–plastic deformation of glassy polymers from low to high strain rates. Int. J. Solids Struct. 2006, 43, 1331–1356. [Google Scholar] [CrossRef]
  9. Richeton, J.; Ahzi, S.; Vecchio, K.S.; Jiang, F.C.; Adharapurapu, R.R. Influence of temperature and strain rate on the mechanical behaviour of three amorphous polymers: Characterization and modeling of the compressive yield stress. Int. J. Solids Struct. 2006, 43, 2318–2335. [Google Scholar] [CrossRef]
  10. Omar, M.F.; Akil, H.M.; Ahmad, Z.A. Measurement and prediction of compressive properties of polymers at high strain rate loading. Mater. Des. 2011, 32, 4207–4215. [Google Scholar] [CrossRef]
  11. Dar, U.A.; Zhang, W.H.; Xu, Y.J.; Wang, J. Thermal and strain rate sensitive compressive behavior of polycarbonate polymer-experimental and constitutive analysis. J. Polym. Res. 2014, 21, 519–530. [Google Scholar] [CrossRef]
  12. Khan, A.S.; Lopez-Pamies, O.; Kazmi, R. Thermo-mechanical large deformation response and constitutive modeling of viscoelastic polymers over a wide range of strain rates and temperatures. Int. J. Plast. 2006, 22, 581–601. [Google Scholar] [CrossRef]
  13. Chen, W.; Lu, F.; Cheng, M. Tension and compression tests of two polymers under quasi-static and dynamic loading. Polym. Test. 2002, 21, 113–121. [Google Scholar] [CrossRef]
  14. Cheng, M.; Chen, W. Experimental investigation of the stress-stretch behavior of EPDM rubber with loading rate effects. Int. J. Solids Struct. 2003, 40, 4749–4768. [Google Scholar] [CrossRef]
  15. Rae, P.; Brown, E. The properties of poly(tetrafluoroethylene) (PTFE) in tension. Polymer 2005, 46, 8128–8140. [Google Scholar] [CrossRef]
  16. Yin, Z.N.; Wang, T.J. Deformation response and constitutive modeling of PC, ABS and PC/ABS alloys under impact tensile loading. Mater. Sci. Eng. A 2010, 527, 1461–1468. [Google Scholar] [CrossRef]
  17. Sarva, S.S.; Boyce, M.C. Mechanics of polycarbonate during high-rate tension. J. Mech. Mater. Struct. 2007, 2, 1853–1880. [Google Scholar] [CrossRef]
  18. Cao, K.; Wang, Y.; Wang, Y. Experimental investigation and modeling of the tension behavior of polycarbonate with temperature effects from low to high strain rates. Int. J. Solids Struct. 2014, 51, 2539–2548. [Google Scholar] [CrossRef]
  19. Anand, L.; Ames, N.M.; Srivastava, V.; Chester, S.A. A thermo-mechanically coupled theory for large deformations of amorphous polymers, part I: Formulation. Int. J. Plast. 2009, 25, 1474–1494. [Google Scholar] [CrossRef]
  20. Bouvard, J.L.; Ward, D.K.; Hossain, D.; Marin, E.B.; Bammann, D.J.; Horstemeyer, M.F. A general inelastic internal state variable model for amorphous glassy polymers. Acta Mech. 2010, 213, 71–96. [Google Scholar] [CrossRef]
  21. Richeton, J.; Ahzi, S.; Vecchio, K.S.; Jiang, F.C.; Makradi, A. Modeling and validation of the large deformation inelastic response of amorphous polymers over a wide range of temperatures and strain rates. Int. J. Solids Struct. 2007, 44, 7938–7954. [Google Scholar] [CrossRef]
  22. Xu, Z.; Li, Y.; Huang, F. Application of split Hopkinson tension bar technique to study of dynamic fracture properties of materials. Acta Mech. Sin. 2012, 28, 424–431. [Google Scholar] [CrossRef]
  23. Wang, J.; Xu, Y.J.; Moumni, Z.; Zhang, W.H. A new damage-based model for the non-linear behavior of polycarbonate polymers. Mater. Des. 2016. under review. [Google Scholar]
  24. De Souza Neto, E.A.; Peric, D.; Owen, D.R.J. Computational Methods for Plasticity; Wiley: New York, NY, USA, 2008. [Google Scholar]
  25. Lemaitre, J.; Lippmann, H. A Course on Damage Mechanics; Springer: Berlin, Germany, 1996. [Google Scholar]
Figure 1. Geometry and dimension of the specimen for SHTB test.
Figure 1. Geometry and dimension of the specimen for SHTB test.
Polymers 08 00063 g001
Figure 2. The SHTB test system: (a) schematic diagram; and (b) photograph.
Figure 2. The SHTB test system: (a) schematic diagram; and (b) photograph.
Polymers 08 00063 g002
Figure 3. The fixture device connecting the specimen in quasi-static tension tests.
Figure 3. The fixture device connecting the specimen in quasi-static tension tests.
Polymers 08 00063 g003
Figure 4. Typical signals measured from the gages on the bars.
Figure 4. Typical signals measured from the gages on the bars.
Polymers 08 00063 g004
Figure 5. Velocity profiles of the incident bar–specimen interface and transmission bar–specimen interface.
Figure 5. Velocity profiles of the incident bar–specimen interface and transmission bar–specimen interface.
Polymers 08 00063 g005
Figure 6. The fixture device connecting the specimen in quasi-static tension tests.
Figure 6. The fixture device connecting the specimen in quasi-static tension tests.
Polymers 08 00063 g006
Figure 7. True tensile stress–true strain laws of PC: (a) quasi-static; and (b) dynamic.
Figure 7. True tensile stress–true strain laws of PC: (a) quasi-static; and (b) dynamic.
Polymers 08 00063 g007
Figure 8. Yield stress of PC over a wide range of strain rates.
Figure 8. Yield stress of PC over a wide range of strain rates.
Polymers 08 00063 g008
Figure 9. Strain softening behavior of PC for different strain rates.
Figure 9. Strain softening behavior of PC for different strain rates.
Polymers 08 00063 g009
Figure 10. Finite element model of specimen.
Figure 10. Finite element model of specimen.
Polymers 08 00063 g010
Figure 11. Numerical predictions and experimental results of stress–strain responses of PC at different strain rates: (a) 0.01 s−1; (b) 0.1 s−1; (c) 3500 s−1; and (d) 4500 s−1.
Figure 11. Numerical predictions and experimental results of stress–strain responses of PC at different strain rates: (a) 0.01 s−1; (b) 0.1 s−1; (c) 3500 s−1; and (d) 4500 s−1.
Polymers 08 00063 g011
Figure 12. Numerical predictions and experimental results of yield stresses of PC.
Figure 12. Numerical predictions and experimental results of yield stresses of PC.
Polymers 08 00063 g012
Table 1. Numerical results and the high-speed photographs during the test at strain rate of 4500 s1.
Table 1. Numerical results and the high-speed photographs during the test at strain rate of 4500 s1.
Time (μs)The contours of axial stress (MPa)Photograph of the gauge region
20 Polymers 08 00063 i001 Polymers 08 00063 i002 Polymers 08 00063 i003
30 Polymers 08 00063 i004 Polymers 08 00063 i005 Polymers 08 00063 i006
55 Polymers 08 00063 i007 Polymers 08 00063 i008 Polymers 08 00063 i009
65 Polymers 08 00063 i010 Polymers 08 00063 i011 Polymers 08 00063 i012
75 Polymers 08 00063 i013 Polymers 08 00063 i014 Polymers 08 00063 i015
100 Polymers 08 00063 i016 Polymers 08 00063 i017 Polymers 08 00063 i018
120 Polymers 08 00063 i019 Polymers 08 00063 i020 Polymers 08 00063 i021
140 Polymers 08 00063 i022 Polymers 08 00063 i023 Polymers 08 00063 i024

Share and Cite

MDPI and ACS Style

Xu, Y.; Gao, T.; Wang, J.; Zhang, W. Experimentation and Modeling of the Tension Behavior of Polycarbonate at High Strain Rates. Polymers 2016, 8, 63. https://doi.org/10.3390/polym8030063

AMA Style

Xu Y, Gao T, Wang J, Zhang W. Experimentation and Modeling of the Tension Behavior of Polycarbonate at High Strain Rates. Polymers. 2016; 8(3):63. https://doi.org/10.3390/polym8030063

Chicago/Turabian Style

Xu, Yingjie, Tenglong Gao, Jun Wang, and Weihong Zhang. 2016. "Experimentation and Modeling of the Tension Behavior of Polycarbonate at High Strain Rates" Polymers 8, no. 3: 63. https://doi.org/10.3390/polym8030063

APA Style

Xu, Y., Gao, T., Wang, J., & Zhang, W. (2016). Experimentation and Modeling of the Tension Behavior of Polycarbonate at High Strain Rates. Polymers, 8(3), 63. https://doi.org/10.3390/polym8030063

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop