Multi-cycle deformation of semicrystalline polymers: Observations and constitutive modeling

Experimental data are reported on isotactic polypropylene in multi-cycle uniaxial tensile tests where a specimen is stretched up to some maximum strain and retracted down to the zero minimum stress, while maximum strains increase with number of cycles. Fading memory of deformation history is observed: when two samples are subjected to loading programs that differ along the first n − 1 cycles only, their stress–strain diagrams coincide starting from the nth cycle. Constitutive equations are developed in cyclic viscoelasticity and viscoplasticity of semicrystalline polymers, and adjustable parameters in the stress–strain relations are found by fitting the experimental data. Results of numerical simulation demonstrate that the model predicts the fading memory effect quantitatively. To confirm that the observed phenomenon is typical of semicrystalline polymers, experimental data are presented in tensile cyclic tests with large strains on low density polyethylene and compressive cyclic tests on poly(oxymethylene).     

Thermo-viscoelastic and viscoplastic behavior of high-density polyethylene

Observations are reported on high-density polyethylene in uniaxial tensile tests with constant strain rates and relaxation tests at various temperatures ranging from 25 to 90 °C. A constitutive model is derived for the nonlinear viscoelastic and viscoplastic behavior of semi-crystalline polymers at three-dimensional deformations. Adjustable parameters in the stress–strain relations are found by fitting the experimental data. It is demonstrated that (i) the model correctly approximates the observations and (ii) material parameters are independent of strain rate and change consistently with temperature.     

Plastic anisotropy and the role of non-basal slip in magnesium alloy AZ31B

Mechanistic explanations for the plastic behavior of a wrought magnesium alloy are developed using a combination of experimental and simulation techniques. Parameters affecting the practical sheet formability, such as strain hardening rate, strain rate sensitivity, the degree of anisotropy, and the stresses and strains at fracture, are examined systematically by conducting tensile tests of variously oriented samples at a range of temperatures (room temperature to 250 °C) and strain rates (10⁻⁵–0.1 s⁻¹). Polycrystal plasticity simulations are used to model the observed anisotropy and texture evolution. Strong in-plane anisotropy observed at low temperatures is attributed to the initial texture and the greater than anticipated non-basal cross-slip of dislocations with 〈a〉 type Burgers vectors. The agreement between the measured and simulated anisotropy and texture is further validated by direct observations of the dislocation microstructures using transmission electron microscopy. The increase in the ductility with temperature is accompanied by a decrease in the flow stress, an increase in the strain rate sensitivity, and a decrease in the normal anisotropy. Polycrystal simulations indicate that an increased activity of non-basal, 〈c + a〉, dislocations provides a self-consistent explanation for the observed changes in the anisotropy with increasing temperature.     

The characteristics of plastic flow and a physically-based model for 3003 Al–Mn alloy upon a wide range of strain rates and temperatures

The uniaxial compressive responses of 3003 Al–Mn alloy upon strain rates ranging from 0.001/s to about 10⁴/s with initial temperatures from 77 K to 800 K were investigated. Instron servohydraulic testing machine and enhanced split Hopkinson bar facilities have been employed in such uniaxial compressive loading tests. The maximum true strain up to 80% has been achieved. The following observations have been obtained from the experimental results: 1) 3003 Al–Mn alloy presents remarkable ductility and plasticity at low temperatures and high strain rates; 2) its plastic flow stress strongly depends on the applied temperatures and strain rates; 3) the temperature history during deformation strongly affects the microstructure evolution within the material. Finally, paralleled with the systematic experimental investigations, a physically-based model was developed based on the mechanism of dislocation kinetics. The model predictions are compared with the experimental results, and a good agreement has been observed.     

Integral-based constitutive equation for rubber at high strain rates

High-speed experiments were conducted to characterize the deformation and failure of Styrene Butadiene Rubber at impact rates. Dynamic tensile stress–strain curves of uniaxial strip specimens and force–extension curves of thin sheets were obtained from a Charpy tensile impact apparatus. Results from the uniaxial tension tests indicated that although the rubber became stiffer with increasing strain rates, the stress–strain curves remained virtually the same above 280 s⁻¹. Above this critical strain rate, strength, fracture strain and toughness decreased with increasing strain rates. When strain rates were below 180 s⁻¹, the initial modulus, tensile strength and breaking extension increased as the strain rate increased. Between strain rates of 180 and 280 s⁻¹, the initial modulus and tensile strength increased with increasing strain rates but the extension at break decreased with increasing strain rates. A hyper-viscoelastic constitutive relation of integral form was used to describe the rate-dependent material behavior of the rubber. Two characteristic relaxation times, 5 ms and 0.25 ms, were needed to fit the proposed constitutive equation to the data. The proposed constitutive equation was implemented in ABAQUS Explicit via a user-defined subroutine and used to predict the dynamic response of the rubber sheets in the experiments. Numerical predictions for the transient deformation and failure of the rubber sheet were within 10% of experimental results.     

Thermo-mechanical response of nylon 101 under uniaxial and multi-axial loadings: Part I,Experimental results over wide ranges of temperatures and strain rates

A comprehensive study of the thermo-mechanical response of a thermoplastic polymer, nylon 101 is presented. Quasi-static and dynamic compression uniaxial and multi-axial experiments (stress states) were performed at a wide range of strain rates (10⁻⁵ to 5000 s⁻¹) and temperatures (−60 to 177 °C or −76 to 350 °F). The material is found to be non-linearly dependent on strain rate and temperature. The change in volume after plastic deformation is investigated and is found to be negligibly small. The relaxation and creep responses at room temperature are found to be dependent on strain rate and the stress–strain level at which these phenomena are initiated. Total deformation is decomposed into visco-elastic and visco-plastic components; these components have been determined at different levels of deformation. Results from non-proportional uniaxial to biaxial compression, and torsion experiments, are also reported for three different strain rates at room temperature. It is shown that nylon 101 has a response dependent on the hydrostatic pressure.     

Influence of temperature and strain rate on the mechanical behavior of three amorphous polymers: Characterization and modeling of the compressive yield stress

Uniaxial compression stress–strain tests were carried out on three commercial amorphous polymers: polycarbonate (PC), polymethylmethacrylate (PMMA), and polyamideimide (PAI). The experiments were conducted under a wide range of temperatures (−40 °C to 180 °C) and strain rates (0.0001 s⁻¹ up to 5000 s⁻¹). A modified split-Hopkinson pressure bar was used for high strain rate tests. Temperature and strain rate greatly influence the mechanical response of the three polymers. In particular, the yield stress is found to increase with decreasing temperature and with increasing strain rate. The experimental data for the compressive yield stress were modeled for a wide range of strain rates and temperatures according to a new formulation of the cooperative model based on a strain rate/temperature superposition principle. The modeling results of the cooperative model provide evidence on the secondary transition by linking the yield behavior to the energy associated to the β mechanical loss peak. The effect of hydrostatic pressure is also addressed from a modeling perspective.     

Residual stresses in thin film systems: Effects of lattice mismatch,thermal mismatch and interface dislocations

This paper explores the mechanisms of the residual stress generation in thin film systems with large lattice mismatch strain, aiming to underpin the key mechanism for the observed variation of residual stress with the film thickness. Thermal mismatch, lattice mismatch and interface misfit dislocations caused by the disparity of the material layers were investigated in detail. The study revealed that the thickness-dependence of the residual stresses found in experiments cannot be elucidated by thermal mismatch, lattice mismatch, or their coupled effect. Instead, the interface misfit dislocations play the key role, leading to the variation of residual stresses in the films of thickness ranging from 100 nm to 500 nm. The agreement between the theoretical analysis and experimental results indicates that the effect of misfit dislocation is far-reaching and that the elastic analysis of dislocation, resolved by the finite element method, is sensible in predicting the residual stress distribution. It was quantitatively confirmed that dislocation density has a significant effect on the overall film stresses, but dislocation distribution has a negligible influence. Since the lattice mismatch strain varies with temperature, it was finally confirmed that the critical dislocation density that leads to the measured residual stress variation with film thickness should be determined from the lattice mismatch strain at the deposition temperature.     

Temperature dependent bulge test for elastomers

The bulge test is a particularly convenient testing method for characterizing elastomers under biaxial loading. In addition, it is convenient to utilize this test for validating material models in simulation due to the heterogeneous strain field induced during inflation. During the bulge test the strain field for elastomers covers uniaxial tension at the border to pure shear and equibiaxial tension at the pole. Elastomeric materials exhibit a hyperelastic material behavior, with a dependency on temperature and loading rate. The temperature effect on the mechanical behavior during biaxial loading is considered in the present study. A bulge test setup combined with a temperature chamber is developed in order to characterize this effect, and an exemplary temperature dependent characterization of a poly(norbornene) elastomer is performed with this setup. The equibiaxial stress–strain curves measured at 60 °C, 20 °C and −20 °C are presented.     

Lattice orientation effects on void growth and coalescence in fcc single crystals

Void growth and coalescence in fcc single crystals were studied using crystal plasticity under uniaxial and biaxial loading conditions and various orientations of the crystalline lattice. A 2D plane strain unit cell with one and two cylindrical voids was employed using three-dimensional 12 potentially active slip systems. The results were compared to five representative orientations of the tensile axis on the stereographic triangle. For uniaxial tension conditions, the void volume fraction increase under the applied load is strongly dependent on the crystallographic orientation with respect to the tensile axis. For some orientations of the tensile axis, such as [1 0 0] or [1 1 0], the voids exhibited a growth rate twice as fast compared with other orientations ([1 0 0], [2 1 1]). Void growth and coalescence simulations under uniaxial loading indicated that during deformation along some orientations with asymmetry of the slip systems, the voids experienced rotation and shape distortion, due mainly to lattice reorientation. Coalescence effects are shown to diminish the influence of lattice orientation on the void volume fraction increase, but noteworthy differences are still present. Under biaxial loading conditions, practically all differences in the void volume fraction for different orientations of the tensile axes during void growth vanish. These results lead to the conclusion that at microstructural length scales in regions under intense biaxiality/triaxiality conditions, such as crack tip or notched regions, the plastic anisotropy due to the initial lattice orientation has only a minor role in influencing the void growth rate. In such situations, void growth and coalescence are mainly determined by the stress triaxiality, the magnitude of accumulated strain, and the spatial localization of such plastic strains.     

Wavelet decomposition method decoupled boiling/evaporation oscillation mechanisms over two to three timescales: A study for a microchannel with pin fin structure

Boiling/evaporation heat transfer in a microchannel with pin fin structure was performed with water as the working fluid. Simultaneous measurements of various parameters were performed. The chip wall temperatures were measured by a high spatial-time resolution IR image system, having a sensitivity of 0.02 °C. The flow pattern variations synchronously changed wall temperatures due to ultra-small Bi number. The wavelet decomposition method successfully identified the noise signal and decoupled various temperature oscillations with different amplitudes and frequencies. Three types of temperature oscillations were identified according to heat flux q and mass flux G. The first type of oscillation occurred at q/G < 0.62 kJ/kg. The approximation coefficient of wavelet decomposition decided the dominant cycle period which was ∼3 times of the fluid residence time in the microchannel, behaving the density wave oscillation characteristic. The detail coefficients of wavelet decomposition decided the dominant cycle period, which matched the flow pattern transition determined value well, representing the flow pattern transition induced oscillation. For the second type of oscillation, the wavelet decomposition decoupled the three oscillation mechanisms. The pressure drop oscillation caused the temperature oscillation amplitudes of 5–10 °C and cycle periods of 10–15 s. The density wave oscillation and flow pattern transition induced oscillation are embedded with both the pressure rise and decrease stages of the pressure drop oscillation. The third type of oscillation happened at q/G > 1.13 kJ/kg, having the density wave oscillation coupled with the varied liquid film evaporation induced oscillation. The liquid island, retention bubble induced nucleation sites and cone-shape two-phase developing region are unique features of microchannel boiling with pin fin structure. This study illustrated that pressure drop oscillation and density wave oscillation, usually happened in large size channels, also take place in microchannels. The flow pattern transition and varied liquid film evaporation induced oscillations are specific to microchannel boiling/evaporation flow.     

Finite element micromechanics model of impact compression of closed-cell polymer foams

Finite element analysis, of regular Kelvin foam models with all the material in uniform-thickness faces, was used to predict the compressive impact response of low-density closed-cell polyethylene and polystyrene foams. Cell air compression was analysed, treating cells as surface-based fluid cavities. For a typical 1 mm cell size and 50 s^?1 impact strain rate, the elastic buckling of cell faces, and pop-in shape inversion of some buckled square faces, caused a non-linear stress strain response before yield. Pairs of plastic hinges formed across hexagonal faces, then yield occurred when trios of faces concertinaed. The predicted compressive yield stresses were close to experimental data, for a range of foam densities. Air compression was the hardening mechanism for engineering strains <0.6, with face-to-face contact also contributing for strains >0.7. Predictions of lateral expansion and residual strains after impact were reasonable. There were no significant changes in the predicted behavior at a compressive strain rate of 500 s^?1.     

An experimental and numerical study of the localization behavior of tantalum and stainless steel

In general, the shear localization process involves initiation and growth. Initiation is expected to be a stochastic process in material space where anisotropy in the elastic–plastic behavior of single crystals and inter-crystalline interactions serve to form natural perturbations to the material’s local stability. A hat-shaped sample geometry was used to study shear localization growth. It is an axi-symmetric sample with an upper “hat” portion and a lower “brim” portion with the shear zone located between the hat and brim. The shear zone length is 870–890 μm with deformation imposed through a split-Hopkinson pressure bar system at maximum top-to-bottom velocity in the range of 8–25 m/s. We present experimental results of the deformation response of tantalum and 316L stainless steel samples. The tantalum samples did not form shear bands but the stainless steel sample formed a late stage shear band. We have also modeled these experiments using both conductive and adiabatic continuum models. An anisotropic elasto-viscoplastic constitutive model with damage evolution was used within the finite element code EPIC. A Mie-Gruneisen equation of state and the rate and temperature sensitive MTS flow stress model together with a Gurson flow surface were employed. The models performed well in predicting the experimental data. The numerical results for tantalum suggested a maximum equivalent strain rate on the order of 7 × 10⁴ s⁻¹ in the gage section for an imposed top surface displacement rate of 17.5 m/s. The models also suggested that for an initial temperature of 298 K a temperature in the neighborhood of 900 K was reached within the shear section. The numerical results for stainless steel suggest that melting temperature was reached throughout the shear band shortly after peak load. Due to sample geometry, the stress state in the shear zone was not pure shear; a significant normal stress relative to the shear zone basis line was developed.     

Three dimensional analytical solution for finite circular cylinders subjected to indirect tensile test

This paper derives a new three-dimensional (3-D) analytical solution for the indirect tensile tests standardized by ISRM (International Society for Rock Mechanics) for testing rocks, and by ASTM (American Society for Testing and Materials) for testing concretes. The present solution for solid circular cylinders of finite length can be considered as a 3-D counterpart of the classical two dimensional (2-D) solutions by Hertz in 1883 and by Hondros in 1959. The contacts between the two steel diametral loading platens and the curved surfaces of a cylindrical specimen of length H and diameter D are modeled as circular-to-circular Hertz contact and straight-to-circular Hertz contact for ISRM and ASTM standards respectively. The equilibrium equations of the linear elastic circular cylinder of finite length are first uncoupled by using displacement functions, which are then expressed in infinite series of some combinations of Bessel functions, hyperbolic functions, and trigonometric functions. The applied tractions are expanded in Fourier–Bessel series and boundary conditions are used to yield a system of simultaneous equations. For typical rock cylinders of 54 mm diameter subjected to ISRM indirect tensile tests, the contact width is in the order of 2 mm (or a contact angle of 4°) whereas for typical asphalt cylinders of 101.6 mm diameter subjected to ASTM indirect tensile tests the contact width is about 10 mm (or a contact angle of 12°). For such contact conditions, 50 terms in both Fourier and Fourier–Bessel series expansions are found sufficient in yielding converged solutions. The maximum hoop stress is always observed within the central portion on a circular section close to the flat end surfaces. The difference in the maximum hoop stress between the 2-D Hondros solution and the present 3-D solution increases with the aspect ratio H/D as well as Poisson’s ratio ν. When contact friction is neglected, the effect of loading platen stiffness on tensile stress in cylinders is found negligible. For the aspect ratio of H/D = 0.5 recommended by ISRM and ASTM, the error in tensile strength may be up to 15% for both typical rocks and asphalts, whereas for longer cylinders with H/D up to 2 the error ranges from 15% for highly compressible materials, and to 60% for nearly incompressible materials. The difference in compressive radial stress between the 2-D Hertz solution or 2-D Hondros solution and the present 3-D solution also increases with Poisson’s ratio and aspect ratio H/D. In summary, the 2-D solution, in general, underestimates the maximum tensile stress and cannot predict the location of the maximum hoop stress which typically locates close to the end surfaces of the cylinder.     

Elastic–plastic behavior of steel sheets under in-plane cyclic tension–compression at large strain

Elastic–plastic behavior of two types of steel sheets for press-forming (an aluminum-killed mild steel and a dual-phase high strength steel of 590 MPa ultimate tensile strength) under in-plane cyclic tension–compression at large strain (up to 25% strain for mild steel and 13% for high strength steel) have been investigated. From the experiments, it was found that the cyclic hardening is strongly influenced by cyclic strain range and mean strain. Transient softening and workhardening stagnation due to the Bauschinger effect, as well as the decrease in Young's moduli with increasing prestrain, were also observed during stress reversals. Some important points in constitutive modeling for such large-strain cyclic elasto-plasticity are discussed by comparing the stress–strain responses calculated by typical constitutive models of mixed isotropic–kinematic hardening with the corresponding experimental observations.     

Experimental investigation of the transverse mechanical properties of a single Kevlar® KM2 fiber

A new experimental setup is developed to investigate the transverse mechanical properties of Kevlar^® KM2 fibers, which has been widely used in ballistic impact applications. Experimental results for large deformation reveal that the Kevlar^® KM2 fibers possess nonlinear, pseudo-elastic transverse mechanical properties. A phenomenon similar to the Mullins effect (stress softening) in rubbers exists for the Kevlar^® KM2 fibers. Large transverse deformation does not significantly reduce the longitudinal tensile load-bearing capacity of the fibers. In addition, longitudinal tensile loads stiffen the fibers' transverse nominal stress–strain behaviors at large transverse deformation. Loading rates have insignificant effects on their transverse mechanical properties even in the finite deformation range. An analytical relationship between transverse compressive force and displacement is derived at infinitesimal strain level. This relation is used to estimate the transverse elastic modulus of the Kevlar^® KM2 fibers, which is 1.34 ± 0.35 GPa.     

High strain-rate behavior of ice under uniaxial compression

In the present study, a modified split Hopkinson pressure bar (SHPB) is employed to investigate the dynamic response of ice under uniaxial compression in the range of strain rates from 60 to 1400 s^?1 and at initial test temperatures of ?10 and ?30 °C. The compressive strength of ice shows positive strain-rate sensitivity over the range of strain rates employed; a slight influence of ice microstructure is observed, but it is much less than that reported previously for ice deformation under quasi-static loading conditions [Schulson, E.M., IIiescu, D., Frott, A., 2005. Characterization of ice for return-to-flight of the space shuttle. Part 1 – Hard ice. NASA CR-2005-213643-Part 1]. Specimen thickness, within the range studied, was found to have little or no effect on the peak (failure) strength of ice, while lowering the test temperature from ?10 to ?30 °C had a considerable effect, with ice behaving stronger at the lower test temperature. Moreover, unlike in the case of uniaxial quasi-static compression of ice, the effect of specimen end-constraint during the high rate compression was found to be negligible. One important result of these experiments, which may have important implications in modeling ice impacts, involves the post “peak-stress” behavior of the ice in that the ice samples do not catastrophically lose their load carrying capacity even after the attainment of peak stress during dynamic compression. This residual (tail) strength of the damaged/fragmented ice is sizable, and in some cases is larger than the quasi-static compression strength reported for ice. Moreover, this residual strength is observed to be dependent on sample thickness and the strain rate, being higher for thinner samples and at higher strain-rates during dynamic compression.     

Spatially resolved wall temperature measurements during flow boiling in microchannels

Spatial and temporal variations of channel wall temperature during flow boiling microchannel flows using infrared thermography are presented and analyzed. In particular, the top channel wall temperature in a branching microchannel silicon heat sink is measured non-intrusively. Using this technique, time-averaged temperature measurements, with a spatial resolution of 10 μm, are presented over an 18 mm × 18 mm area of the heat sink. Also presented, within a specific sub-region of the heat sink, are intensity maps that are recorded at a rate of 120 frames per second. Time series data at selected point locations in this sub-region are analyzed for their frequency content, and dominant temperature fluctuations are extracted using proper orthogonal decomposition.Results at low-vapor-quality boiling condition indicate that temperatures can be determined from recorded radiation intensities with a temperature uncertainty varying from 0.9 °C at 25 °C to 1.0 °C at 125 °C. The time series data indicate periodic wall temperature fluctuations of approximately 2 °C that are attributed to the passage of vapor slugs. A dominant band of frequencies around 2–4 Hz is suggested by the frequency analysis. Proper orthogonal decomposition results indicate that first six orthogonal modes account for approximately 90% of the variance in temperature. The first mode reconstruction accounts for temporal variations in the dataset in the sub-region analyzed; however the magnitude of fluctuations and spatial variations in temperature are not accurately captured. A reconstruction using the first 25 modes is considered sufficient to capture both the temporal and spatial variations in the data.     

Superelastic shape memory alloy cables: Part II – Subcomponent isothermal responses

This paper constitutes the second part of our experimental study of the thermo-mechanical behavior of superelastic NiTi shape memory alloy cables. Part I introduced the fundamental, room temperature, tensile responses of two cable designs (7 × 7 right regular lay, and 1 × 27 alternating lay). In Part II, each cable behavior is studied further by breaking down the response into the contributions of its hierarchical subcomponents. Selected wire strands were extracted from the two cable constructions, and their quasi-static tension responses were measured using the same experimental setup of Part I. Consistent with the shallow wire helix angles in the 7 × 7 construction, the force–elongation responses of the core wire, 1 × 7 core strand and full 7 × 7 cable were similar on a normalized basis, with only a slight decrease in transformation force plateaus and slight increase in plateau strains in this specimen sequence. By contrast, each successive 1 × 27 component (1 × 6 core strand, 1 × 15 strand, and full cable) included an additional outer layer of wires with a larger number of wires, greater helix radius, and deeper helix angle, so the normalized axial load responses became significantly more compliant. Each specimen in the sequence also exhibited progressively larger strains at failure, reaching 40% strain in the full 1 × 27 cable.Stress-induced phase transformations involved localized strain/temperature and front propagation in all of the tested 7 × 7 components but none of the 1 × 27 components aside from the 1 × 27 core wire. Stereo digital image correlation measurements revealed finer features within a global transformation front of the 1 × 7 core strand than the 7 × 7 cable, consisting of an staggered pattern of individual wire fronts that moved in lock-step during elongation. Although the 1 × 27 multi-layer strands exhibited temperature/strain localizations in a distributed pattern during transformations, the localizations did not propagate and their cause was traced back to contact indentations (stress concentrations) arising from the cable’s fabrication. The normalized axial torque responses of the multi-layer 1 × 27 components during transformation were distinctly non-monotonic and complex, due to the alternating handedness of the layers. Force and torque contributions of individual wire layers were deduced by subtracting 1 × 27 component responses, which helped to clarify the transformation kinetics within each layer and explain the unusual force and torque undulations seen in the 1 × 27 cable response of Part I.