Dynamics of Granular Crystals with Elastic-Plastic Contacts

We study the behavior of granular crystals subjected to impact loading that creates plastic deformation at the contacts between constituent particles. Granular crystals are highly periodic arrangements of spherical particles, arranged into densely packed structures resembling crystals. This special class of granular materials has been shown to have unique dynamics with suggested applications in impact protection. However, previous work has focused on very low amplitude impacts where every contact point can be described using the Hertzian contact law, valid only for purely elastic deformation. In this thesis, we extend previous investigation of the dynamics of granular crystals to significantly higher impact energies more suitable for the majority of applications. Additionally, we demonstrate new properties specific to elastic-plastic granular crystals and discuss their potential applications as well. We first develop a new contact law to describe the interaction between particles for large amplitude compression of elastic-plastic spherical particles including a formulation for strain-rate dependent plasticity. We numerically and experimentally demonstrate the applicability of this contact law to a variety of materials typically used in granular crystals. We then extend our investigation to one-dimensional chains of elastic-plastic particles, including chains of alternating dissimilar materials. We show that, using the new elastic-plastic contact law, we can predict the speed at which impact waves with plastic dissipation propagate based on the material properties of the constituent particles. Finally, we experimentally and numerically investigate the dynamics of two-dimensional and three-dimensional granular crystals with elastic-plastic contacts. We first show that the predicted wave speeds for 1D granular crystals can be extended to 2D and 3D materials. We then investigate the behavior of waves propagating across oblique interfaces of dissimilar particles. We show that the character of the refracted wave can be predicted using an analog to Snell's law for elastic-plastic granular crystals and ultimately show how it can be used to design impact guiding "lenses" for mitigation applications.

Investigation of the effects of extracellular osmotic pressure on morphology and mechanical 1 properties of individual chondrocyte

It has been demonstrated that most cells of the body respond to osmotic pressure in a systematic manner. The disruption of the collagen network in the early stages of osteoarthritis causes an increase in water content of cartilage which leads to a reduction of pericellular osmolality in chondrocytes distributed within the extracellular environment. It is therefore arguable that an insight into the mechanical properties of chondrocytes under varying osmotic pressure would provide a better understanding of chondrocyte mechanotransduction and potentially contribute to knowledge on cartilage degeneration. In this present study, the chondrocyte cells were exposed to solutions with different osmolality. Changes in their dimensions and mechanical properties were measured over time. Atomic Force Microscopy (AFM) was used to apply load at various strain-rates and the force-time curves were logged. The thin-layer elastic model was used to extract the elastic stiffness of chondrocytes at different strain-rates and at different solution osmolality. In addition, the porohyperelastic (PHE) model was used to investigate the strain-rate dependent responses under the loading and osmotic pressure conditions. The results revealed that the hypo-osmotic external environment increased chondrocyte dimensions and reduced Young’s modulus of the cells at all strain-rates tested. In contrast, the hyper-osmotic external environment reduced dimensions and increased Young’s modulus. Moreover, by using the PHE model coupled with inverse FEA simulation, we established that the hydraulic permeability of chondrocytes increased with decreasing extracellular osmolality which is consistent with previous work in the literature. This could be due to a higher intracellular fluid volume fraction with lower osmolality.

CAE-based prediction of projectile residual velocity for impact on single and multi-layered metallic armour plates

The present, paper deals with the CAE-based study Of impact of jacketed projectiles on single- and multi-layered metal armour plates using LS-DYNA. The validation of finite element modelling procedure is mainly based on the mesh convergence study using both shell and solid elements for representing single-layered mild steel target plates. It, is shown that the proper choice of mesh density and the strain rate-dependent material properties are essential for all accurate prediction of projectile residual velocity. The modelling requirements are initially arrived at by correlating against test residual velocities for single-layered mild steel plates of different depths at impact velocities in the ran.-c of approximately 800-870 m/s. The efficacy of correlation is adjudged, in terms of a 'correlation index', defined in the paper: for which values close to unity are desirable. The experience gained for single-layered plates is next; used in simulating projectile impacts on multi-layered mild steel target plates and once again a high degree of correlation with experimental residual velocities is observed. The study is repeated for single- and multi-layered aluminium target plates with a similar level of success in test residual velocity prediction. TO the authors' best knowledge, the present comprehensive study shows in particular for the first time that, with a. proper modelling approach, LS-DYNA can be used with a great degree of confidence in designing perforation-resistant single and multi-layered metallic armour plates.

Semi-empirical procedures for estimation of residual velocity and ballistic limit for impact on mild steel plates by projectiles

This paper deals with the development of simplified semi-empirical relations for the prediction of residual velocities of small calibre projectiles impacting on mild steel target plates, normally or at an angle, and the ballistic limits for such plates. It has been shown, for several impact cases for which test results on perforation of mild steel plates are available, that most of the existing semi-empirical relations which are applicable only to normal projectile impact do not yield satisfactory estimations of residual velocity. Furthermore, it is difficult to quantify some of the empirical parameters present in these relations for a given problem. With an eye towards simplicity and ease of use, two new regression-based relations employing standard material parameters have been discussed here for predicting residual velocity and ballistic limit for both normal and oblique impact. The latter expressions differ in terms of usage of quasi-static or strain rate-dependent average plate material strength. Residual velocities yielded by the present semi-empirical models compare well with the experimental results. Additionally, ballistic limits from these relations show close correlation with the corresponding finite element-based predictions.

Robust prediction of residual velocities and ballistic limits of projectiles for impact on thin aluminium plates

This paper deals with the simulation-driven study of the impact of hardened steel projectiles on thin aluminium target plates using explicit finite element analysis as implemented in LS-DYNA. The evaluation of finite element modelling includes a comprehensive mesh convergence study using shell elements for representing target plates and the solid element-based representation of ogivalnosed projectiles. A user-friendly automatic contact detection algorithm is used for capturing interaction between the projectile and the target plate. It is shown that the proper choice of mesh density and strain rate-dependent material properties is crucial as these parameters significantly affect the computed residual velocity. The efficacy of correlation with experimental data is adjudged in terms of a 'correlation index' defined in the present study for which values close to unity are desirable.By simulating laboratory impact tests on thin aluminium plates carried out by earlier investigators, extremely good prediction of experimental ballistic limits has been observed with correlation indices approaching unity. Additional simulation-based parametric studies have been carried out and results consistent with test data have been obtained. The simulation procedures followed in the present study can be applied with confidence in designing thin aluminium armour plates for protection against low calibre projectiles.

Tailoring energy absorption capacity of CNT forests through application of electric field

This study examines the effect of electric field on energy absorption capacity of carbon nanotube forests (CNTFs), comprising of vertically aligned multiwalled carbon nanotubes, under both quasistatic (strain rate, (epsilon) over dot = 10(-3) s(-1)) and dynamic ((epsilon) over dot = similar to 10(3) s(-1)) loading conditions. Under quasistatic condition, the CNTFs were cyclically loaded and unloaded while electric field was applied along the length of carbon nanotube (CNT) either throughout the loading cycle or explicitly during either the loading or the unloading segment. The energy absorbed per cycle by CNTF increased monotonically with electric field when the field was applied only during the loading segment: A 7 fold increase in the energy absorption capacity was registered at an electric field of 1 kV/m whereas no significant change in it was noted for other schemes of electro-mechanical loading. The energy absorption capacity of CNTF under dynamic loading condition also increased monotonically with electric field; however, relative to the quasistatic condition, less pronounced effect was observed. This intriguing strain rate dependent effect of electric field on energy absorption capacity of CNTF is explained in terms of electric field induced strengthening of CNTF, originating from the time dependent electric field induced polarization of CNT. (C) 2015 Elsevier Ltd. All rights reserved.

Shear Localization In Dynamic Deformation: Microstructural Evolution

Investigations made by the authors and collaborators into the microstructural aspects of adiabatic shear localization are critically reviewed. The materials analyzed are low-carbon steels, 304 stainless steel, monocrystalline Fe-Ni-Cr, Ti and its alloys, Al-Li alloys, Zircaloy, copper, and Al/SiCp composites. The principal findings are the following: (a) there is a strain-rate-dependent critical strain for the development of shear bands; (b) deformed bands and white-etching bands correspond to different stages of deformation; (c) different slip activities occur in different stages of band development; (d) grain refinement and amorphization occur in shear bands; (e) loss of stress-carrying capability is more closely associated with microdefects rather than with localization of strain; (f) both crystalline rotation and slip play important roles; and (g) band development and band structures are material dependent. Additionally, avenues for new research directions are suggested.

Dynamic buckling of an inclined strut

The dynamic compressive response of a sandwich plate with a metallic corrugated core is predicted. The back face of the sandwich plate is held fixed whereas the front face is subjected to a uniform velocity, thereby compressing the core. Finite element analysis is performed to investigate the role of material inertia, strain hardening and strain rate hardening upon the dynamic collapse of the corrugated core. Three classes of collapse mode are identified as a function of impact velocity: (i) a three-hinge plastic buckling mode of wavelength equal to the strut length, similar to the quasi-static mode, (ii) a 'buckle-wave' regime involving inertia-mediated plastic buckling of wavelength less than that of the strut length, and (iii) a 'stubbing' regime, with shortening of the struts by local fattening at the front face. The presence of strain hardening reduces the regime of dominance of the stubbing mode. The influence of material strain rate sensitivity is evaluated by introducing strain rate dependent material properties representative of type 304 stainless steel. For this choice of material, strain rate sensitivity has a more minor influence than strain hardening, and consequently the dynamic collapse strength of a corrugated core is almost independent of structural dimension. © 2012 Elsevier Ltd. All rights reserved.

A viscous froth model for dry foams in the Surface Evolver

Cox, S.J. (2005) A Viscous Froth Model for Dry Foams in the Surface Evolver. Coll. Surf. A 263:81-89. Sponsorship: Partial support from both the Ulysses France?Ireland Exchange Scheme and the European Space Agency, Contract 14308/00/NL/SH (AO-99-031) CCN 002 MAP Project AO-99-075.

Characterization of viscoplastic yielding of asphalt concrete

A temperature and strain rate dependent yield surface model was proposed to characterize the viscoplastic yielding of asphalt concrete. Laboratory tests were conducted on specimens that have two binders, two air void contents, and three aging periods. Strain decomposition was performed to obtain viscoplastic strain and stress-pseudostrain curves were constructed to determine the model parameters accurately and efficiently. Results indicate that a stiffer asphalt concrete has greater cohesion and strain hardening amplitude, both of which decline as temperature increases or strain rate decreases. The temperature and strain rate factors of the yield surface can be accurately determined solely by the peak stress of the strength tests. © 2013 Elsevier Ltd. All rights reserved.

Temperature dependent compressive behavior of graphene mediated three-dimensional cellular assembly

Mechanical behavior of three-dimensional cellular assembly of graphene foam (GF) presented temperature dependent characteristics evaluated at both low temperature and room temperature conditions. Cellular structure of GF comprised of polydimethyl siloxane polymer as a flexible supporting material demonstrated 94% enhancement in the storage modulus as compared to polymer foam alone. Evaluation of frequency dependence revealed an increase in both storage modulus and tan delta with the increase in frequency. Moreover, strain rate independent highly reversible behavior is measured up to several compression cycles at larger strains. It is elucidated that the interaction between graphene and polymer plays a crucial role in thermo-mechanical stability of the cellular structure. (C) 2015 Elsevier Ltd. All rights reserved.

Increased time-dependent room temperature plasticity in metallic glass nanopillars and its size-dependency

Room temperature, uniaxial compression creep experiments were performed on micro-/nano-sized pillars (having diameters in the range of 250-2000 nm) of a Zr-based bulk metallic glass (BMG) to investigate the influence of sample size on the time-dependent plastic deformation behavior in amorphous alloys. Experimental results reveal that plastic deformation indeed occurs at ambient temperature and at stresses that are well below the nominal quasi-static yield stress. At a given stress, higher total strains accrue in the smaller specimens. In all cases, plastic deformation was found to be devoid of shear bands, i.e., it occurs in homogeneous manner. The stress exponent obtained from the slope of the linear relation between strain rate and applied stress also shows a strong size effect, which is rationalized in terms of the amount of free volume created during deformation and the surface-to-volume ratio of the pillar. (C) 2012 Elsevier Ltd. All rights reserved.

Compressive flow behavior of Al-Si based alloy: Role of heat treatment

The flow characteristics of a near eutectic Al-Si based cast alloy have been examined in compression at strain rates varying from 3 x 10(-4) to 10(2) s(-1) and at three different temperatures, i.e., room temperature (RT), 100 degrees C and 200 degrees C. The dependence of the flow behavior on heat treatment is studied by testing the alloy in non-heat treated (NHT) and heat treated (HT) conditions. The heat treatment has strong influence on strain rate sensitivity (SRS), strength and work hardening behavior of the alloy. It is observed that the strength of the alloy increases with increase in strain rate and it increases more rapidly above the strain rate of 10(-1) s(-1) in HT condition at all the temperatures, and at 100 degrees C and 200 degrees C in NHT condition. The thermally dependent process taking place in the HT matrix is responsible for the observed greater SRS in HT condition. The alloy in HT condition exhibits a larger work hardening rate than in NHT condition during initial stages of straining. However, the hardening rate decreases more sharply at higher strains in HT condition due to precipitate shearing and higher rate of Si particle fracture. Thermal hardening is observed at 200 degrees C in NHT condition due to precipitate formation, which results in increased SRS at higher temperatures. Thermal softening is observed in HT condition at 200 C due to precipitate coarsening, which leads to a decrease in SRS at higher temperatures. Stress simulations by a finite element method support the experimentally observed particle and matrix fracture behavior. A negative SRS and serrated flow are observed in the lower strain rate regime (3 x 10(-4)-10(-2) s(-1)) at RT and 100 degrees C, in both NHT and HT conditions. The observations show that both dynamic strain aging (DSA) and precipitate shearing play a role in serrated flow. (C) 2015 Elsevier B.V. All rights reserved.

Role of Si modification on the compressive flow behavior of Al-Si based alloy

The flow characteristics of a near-eutectic heat-treated Al-Si based cast alloy have been examined in compression at strain rates varying from 3 x 10(-4) to 10(2) s(-1) and at three different temperatures, i.e., room temperature (RT), 100 degrees C and 200 degrees C. The dependence of flow behavior on modification is examined by testing the alloy in both the unmodified and modified conditions. Modification has strong influence on strain rate sensitivity (SRS), strength and work hardening behavior of the alloy. The strength of the alloy is found to increase with increase in strain rate for both the conditions. The increase is more rapid above the strain rate of 10(-1) s(-1) for the unmodified alloy at all the temperatures. This rapid increase is observed at 1 s(-1) at RT and 100 degrees C, and at 10(-2) s(-1) at 200 degrees C for the modified alloy. The thermally dependent process of the Al matrix is rate controlling in the unmodified alloy. On the other hand, the thermally dependent process of both Al matrix and Si particles are rate controlling, which is responsible for the higher strain rate sensitivity (SRS) in the modified alloy. The unmodified alloy exhibits a larger work hardening rate than the modified alloy during the initial stages of straining due to fiber loading of unmodified Si particles. However, the hardening rate decreases sharply at higher strains for the unmodified alloy due to a higher rate of Si particle fracture. Thermal softening is observed for both alloys at 200 degrees C due to precipitate coarsening, which leads to a decrease in SRS at higher temperatures. Stress simulations by microstructure based finite element method support the experimentally observed particle and matrix fracture behavior. Negative SRS and serrated flow are observed at lower strain rate regime (3 x 10(-4) to 10(-2) s(-1)) at RT and 100 degrees C, in both alloys. The critical onset strain is found to be lower and the magnitude of serration is found to be higher for the modified alloy, which suggests that, in addition to dynamic strain aging, Si particle size and morphology also play a role in serrated flow. (C) 2015 Elsevier Inc All rights reserved.

Strain Gradient Effects on Deformation Strengthening Behavior of Particle Reinforced Metal Matrix Composites

Dislocation models with considering the mismatch of elastic modulus between matrix and reinforcing particles are used to determine the effective strain gradient \ita for particle reinforced metal matrix composites (MMCp) in the present research. Based on Taylor relation and the kinetics of dislocation multiplication, glide and annihilation, a strain gradient dependent constitutive equation is developed. By using this strain gradient-dependent constitutive equation, size-dependent deformation strengthening behavior is characterized. The results demonstrate that the smaller the particle size, the more excellent in the reinforcing effect. Some comparisons with the available experimental results demonstrate that the present approach is satisfactory.