The high strain mechanical response of the wet Kelvin model for open-cell foams

Surface Evolver software was used to create the three-dimensional geometry of a Kelvin open-cell foam, to simulate that of polyurethane flexible foams. Finite Element Analysis (FEA) with 3D elements was used to model large compressive deformation in the [0 0 1] and [1 1 1] directions, using cyclic boundary conditions when necessary, treating the polyurethane as an elastic or elastic–plastic material. The predicted foam Young’s moduli in the [0 0 1] direction are double those of foams with uniform Plateau border cross-section edges, for the same foam density and material properties. For compression in the [1 1 1] direction, the normalized Young’s modulus increases from 0.9 to 1.1 with foam relative density, and the predicted stress–strain relationship can have a plateau, even for a linearly-elastic polymer. As the foam density increases, the predicted effects of material plasticity become larger. For foam of relative density 0.028, edge-to-edge contact is predicted to occur at a 66% strain for [1 1 1] direction compression. The foam is predicted to contract laterally when the [1 1 1] direction compressive strain exceeds 25%.     

Analytical and numerical techniques to predict carbon nanotubes properties

In this paper, two different approaches for modeling the behaviour of carbon nanotubes are presented. The first method models carbon nanotubes as an inhomogeneous cylindrical network shell using the asymptotic homogenization method. Explicit formulae are derived representing Young’s and shear moduli of single-walled nanotubes in terms of pertinent material and geometric parameters. As an example, assuming certain values for these parameters, the Young’s modulus was found to be 1.71 TPa, while the shear modulus was 0.32 TPa. The second method is based on finite element models. The inter-atomic interactions due to covalent and non-covalent bonds are replaced by beam and spring elements, respectively, in the structural model. Correlations between classical molecular mechanics and structural mechanics are used to effectively model the physics governing the nanotubes. Finite element models are developed for single-, double- and multi-walled carbon nanotubes. The deformations from the finite element simulations are subsequently used to predict the elastic and shear moduli of the nanotubes. The variation of mechanical properties with tube diameter is investigated for both zig-zag and armchair configurations. Furthermore, the dependence of mechanical properties on the number of nanotubules in multi-walled structures is also examined. Based on the finite element model, the value for the elastic modulus varied from 0.9 to 1.05 TPa for single and 1.32 to 1.58 TPa for double/multi-walled nanotubes. The shear modulus was found to vary from 0.14 to 0.47 TPa for single-walled nanotubes and 0.37 to 0.62 for double/multi-walled nanotubes.     

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.     

Buckling and progressive crushing of laterally loaded honeycomb

This paper presents a comprehensive study of the lateral compressive response of hexagonal honeycomb panels from the initial elastic regime to a fully crushed state. Expanded aluminum alloy honeycomb panels with a cell size of 9.53 mm, a relative density of 0.026, and a height of 15.9 mm are laterally compressed quasi statically between rigid platens under displacement control. The cells buckle elastically and collapse at a higher stress due to inelastic action. Deformation then first localizes at mid-height and the cells crush by progressive formation of folds; associated with each fold family is a stress undulation. The response densifies when the whole panel height is consumed by folds. The buckling and crushing events are simulated numerically using finite element models involving periodic domains of a single or several characteristic cells. The models idealize the microstructure as hexagonal, with double walls in one direction. The nonlinear behavior is initiated by elastic buckling while inelastic collapse that leads to the localization observed in the experiments occurs at a significantly higher load. The collapse stress is found to be mildly sensitive to various problem imperfections. The subsequent folding can be reproduced numerically using periodic domains but requires a fine mesh capable of capturing the complexity of the folds. The calculated crushing response is shown to better resemble measured ones when a 4 × 4 cell domain is used. However, the average crushing stress can be captured with engineering accuracy even from a single cell domain.     

Discrete element modeling of a Mars Exploration Rover wheel in granular material

Three-dimensional discrete element method (DEM) simulations were developed for the Mars Exploration Rover (MER) mission to investigate: (1) rover wheel interactions with martian regolith; and (2) regolith deformation in a geotechnical triaxial strength cell (GTSC). These DEM models were developed to improve interpretations of laboratory and in situ rover data, and can simulate complicated regolith conditions. A DEM simulation was created of a laboratory experiment that involved a MER wheel digging into lunar regolith simulant. Sinkage and torques measured in the experiment were compared with those predicted numerically using simulated particles of increasing shape complexity (spheres, ellipsoids, and poly-ellipsoids). GTSC simulations, using the same model regolith used in the MER simulations, indicate a peak friction angle of approximately 37–38° compared to internal friction angles of 36.5–37.7° determined from the wheel digging experiments. Density of the DEM regolith was 1820 kg/m³ compared to 1660 kg/m³ for the lunar simulant used in the wheel digging experiment indicating that the number of grain contacts and grain contact resistance determined bulk strength in the DEM simulations, not density. An improved correspondence of DEM and actual test regolith densities is needed to simulate the evolution of regolith properties as density changes.     

Cracked infinite cylinder with two rigid inclusions under axisymmetric tension

This paper considers the problem of an axisymmetric infinite cylinder with a ring shaped crack at z = 0 and two ring-shaped rigid inclusions with negligible thickness at z = ±L. The cylinder is under the action of uniformly distributed axial tension applied at infinity and its lateral surface is free of traction. It is assumed that the material of the cylinder is linearly elastic and isotropic. Crack surfaces are free and the constant displacements are continuous along the rigid inclusions while the stresses have jumps. Formulation of the mixed boundary value problem under consideration is reduced to three singular integral equations in terms of the derivative of the crack surface displacement and the stress jumps on the rigid inclusions. These equations, together with the single-valuedness condition for the displacements around the crack and the equilibrium equations along the inclusions, are converted to a system of linear algebraic equations, which is solved numerically. Stress intensity factors are calculated and presented in graphical form.     

Dynamic crushing of sandwich panels with prismatic lattice cores

The dynamic out-of-plane compressive response of stainless steel corrugated and Y-frame sandwich cores have been investigated for impact velocities ranging from quasi-static to 200 ms⁻¹. Laboratory-scale sandwich cores of relative density 2.5% were manufactured and the stresses on the front and rear faces of the dynamically compressed sandwich cores were measured using a direct impact Kolsky bar. Direct observational evidence is provided for micro-inertial stabilisation of both topologies against elastic buckling at impact velocities below 30 ms⁻¹. At higher impact velocities, plastic waves within the core members result in the front face stresses increasing with increasing velocity while the rear face stresses remain approximately constant. While the finite element calculations predict the rear face stresses and dynamic deformation modes to reasonable accuracy, the relatively slow response time of the measurement apparatus results in poor agreement between the measured and predicted front face stresses. The finite element calculations also demonstrate that material strain-rate effects have a negligible effect upon the dynamic compressive response of laboratory-scale and full-scale sandwich cores.     

Dynamic coupling of fluid and structural mechanics for simulating particle motion and interaction in high speed compressible gas particle laden flow

In this work, structural finite element analyses of particles moving and interacting within high speed compressible flow are directly coupled to computational fluid dynamics and heat transfer analyses to provide more detailed and improved simulations of particle laden flow under these operating conditions. For a given solid material model, stresses and displacements throughout the solid body are determined with the particle–particle contact following an element to element local spring force model and local fluid induced forces directly calculated from the finite volume flow solution. Plasticity and particle deformation common in such a flow regime can be incorporated in a more rigorous manner than typical discrete element models where structural conditions are not directly modeled. Using the developed techniques, simulations of normal collisions between two 1 mm radius particles with initial particle velocities of 50–150 m/s are conducted with different levels of pressure driven gas flow moving normal to the initial particle motion for elastic and elastic–plastic with strain hardening based solid material models. In this manner, the relationships between the collision velocity, the material behavior models, and the fluid flow and the particle motion and deformation can be investigated. The elastic–plastic material behavior results in post collision velocities 16–50% of their pre-collision values while the elastic-based particle collisions nearly regained their initial velocity upon rebound. The elastic–plastic material models produce contact forces less than half of those for elastic collisions, longer contact times, and greater particle deformation. Fluid flow forces affect the particle motion even at high collision speeds regardless of the solid material behavior model. With the elastic models, the collision force varied little with the strength of the gas flow driver. For the elastic–plastic models, the larger particle deformation and the resulting increasingly asymmetric loading lead to growing differences in the collision force magnitudes and directions as the gas flow strength increased. The coupled finite volume flow and finite element structural analyses provide a capability to capture the interdependencies between the interaction of the particles, the particle deformation, the fluid flow and the particle motion.     

Calculating the dissipation in fluid dampers with non-newtonian fluid models

The present paper gives a comparison of the Maxwell, Upperconvected Maxwell and the Oldroyd-B model for the calculation of dissipation in high shear-rate cases. Usage of viscodampers in the automotive industry is the most common. There is a good scope of the computing this power in the case of Newtonian fluids. When a polymeric liquid is considered that part of energy that is irreversible cannot be calculated as P_diss. = τ : d. For fluids where the separation into a solvent and a polymer part is not available but the deformation gradient tensor must be separated into two parts. One part consists of only the elastic deformation while the other is the non-elastic. This paper shows this separation using the Maxwell and the UCM models. A simple problem is shown, solving both analytically and numerically. The steady state temperature distribution of a damper then is validated with measurement.     

Buckling behavior of Kelvin open-cell foams under [0 0 1], [0 1 1] and [1 1 1] compressive loads

This paper describes buckling modes and stresses of elastic Kelvin open-cell foams subjected to [0 0 1], [0 1 1] and [1 1 1] uniaxial compressions. Cubic unit cells and cell aggregates in model foams are analyzed using a homogenization theory of the updated Lagrangian type. The analysis is performed on the assumption that the struts in foams have a non-uniform distribution of cross-sectional areas as observed experimentally. The relative density is changed to range from 0.005 to 0.05. It is thus found that long wavelength buckling and macroscopic instability primarily occur under [0 0 1] and [0 1 1] compressions, with only short wavelength buckling under [1 1 1] compression. The primary buckling stresses under the three compressions are fairly close to one another and almost satisfy the Gibson–Ashby relation established to fit experiments. By also performing the analysis based on the uniformity of strut cross-sectional areas, it is shown that the non-uniformity of cross-sectional areas is an important factor for the buckling behavior of open-cell foams.     

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.     

Stress relaxation of PBI based membrane electrode assemblies

Stress relaxation in the membrane electrode assemblies (MEA) in PEM fuel cells subjected to compressive loads is analyzed. This behavior is important because nonzero contact stress is required to maintain low electric resistivity in the fuel cell stack. Experimental results are used to guide the choice of the viscoelastic properties of the constituents of the MEA, the membrane and the gas diffusion layer (GDL), needed for the model. These properties are incorporated into the model that treats the membrane as a porous-viscoelastic solid, and the gas diffusion layer as a nonlinear elastic solid. Using numerical simulations (finite element method), the stress relaxation curves for the MEA are obtained for different fluid flow boundary conditions, variations in the material properties of the membrane and the GDL. The results are compared to experimental stress relaxation curves. Most of the experimental data were obtained at a temperature of 180 °C, corresponding to operating conditions, so in the model the temperature was considered fixed and equal to this value.     

An analytical and computational study of crack initiation under transient creep conditions

An analytical method has been developed to predict creep crack initiation (CCI), based on the accumulation of a critical level of damage at a critical distance. The method accounts for the re-distribution of stress from the elastic or elastic–plastic field, experienced on initial loading, to a steady state creep stress distribution, via a transient creep region. The method has been applied to predict CCI times in a fracture specimen of type 316H stainless steel at 550 °C. The failure model has been also been implemented into a finite element (FE) framework. Reasonable and conservative predictions of CCI time can be obtained from the analytical solution relative to FE solutions. Conservative predictions of experimental CCI times are obtained when stress redistribution is taking into account. However, CCI times predicted from a steady state creep model are found to be non-conservative.     

Elastic and piezoelectric fields due to polyhedral inclusions

We derive explicit, closed-form expressions describing elastic and piezoelectric deformations due to polyhedral inclusions in uniform half-space and bi-materials. Our analysis is based on the linear elasticity theory and Green’s function method. The method involves evaluation of volume and surface integrals of harmonic and bi-harmonic potentials. In case of polyhedra, such integrals are expressed through algebraic functions. Our results generalize numerous studies on this subject, and they allow to obtain fully analytical solutions for a number of physical and engineering problems. In the limiting case of an infinite space, our relations have an essentially more compact form, than relations obtained by other authors. We present solutions to classical Mindlin and Cherruti problems. We describe the elastic relaxation of a misfitting polygonal quantum dot in bi-materials assuming isotropic and vertically isotropic properties. It is explained how to analyze non-hydrostatic and non-uniform inclusions. We also study piezoelectric fields induced by inclusions in materials with cubic and hexagonal lattices. Among other results, we have found that a cubic inclusion in an isotropic material reproduces fields of quantum dots in GaAs (0, 0, 1) and GaAs (1, 1, 1) depending on the orientation of the cube. This suggests that one can qualitatively model crystals with different lattices by choosing an appropriate inclusion shape.     

Elastic modulus measurements of vertically aligned multi walled carbon nanotubes carpets by using the nanoindentation technique

Electrical, thermal and mechanical properties of Vertically Aligned Multi Walled Carbon NanoTubes (VA-MWCNT) make them an ideal candidate to replace some of conventional materials in micro and nano-electronic components. Integrating this material in micro components requires a good knowledge of their properties. As the electrical and thermal properties, the MWCNT mechanical properties are difficult to assess. Several techniques have been developed to estimate the CNT Young's modulus and the obtained results cover a large range of scale. In this study, we propose an indirect technique for MWCNT carpet Young's modulus measurements by using the nanoindentation technique. Nanoindentation tests are performed on a metallic film deposited on MWCNT. The measured equivalent reduced modulus takes into account the elastic properties of the metallic thin film and those of the MWCNT substrate. Bec et al. model, introduced in 2006, is used to separate elastic properties, and thus determine the MWCNT reduced Young’s modulus which is estimated between 329 and 352 GPa. Knowing the indenter mechanical properties, we estimate the Young’s modulus in the 461 to507 GPa range.     

Computational homogenization of periodic beam-like structures

This paper is concerned with the computation of the effective elastic properties of periodic beam-like structures. The homogenization theory is used and leads to an equivalent anisotropic Navier–Bernoulli–Saint-Venant beam. The overall behavior is obtained from the solution of basic cell problems posed on the three-dimensional period of the structure and solved using three-dimensional finite element implementation.This procedure is first applied to two corrugated zigzag and sinus beams subjected to in-plane loading. Next, the axial elastic properties of a stranded ‘6 + 1’ wire-cable are computed. The effective properties values obtained appear to be very close to analytical reference results showing the efficiency of the approach.     

Irreversible thermodynamics of triple junctions during the intergranular void motion under the electromigration forces

A rigorous reformulation of internal entropy production and the rate of entropy flow is developed for multi-component systems consisting of heterophases, interfaces and/or surfaces. The result is a well-posed moving boundary value problem describing the dynamics of curved interfaces and surfaces associated with voids and/or cracks that are intersected by grain boundaries. Extensive computer simulations are performed for void configuration evolution during intergranular motion. In particular we simulate evolution resulting from the action of capillary and electromigration forces in thin film metallic interconnects having a “bamboo” structure, characterized by grain boundaries aligned perpendicular to the free surface of the metallic film interconnects. Analysis of experimental data utilizing previously derived mean time to failure formulas gives consistent values for interface diffusion coefficients and enthalpies of voids. 3.0 × 10⁻⁶ exp(−0.62 eV/kT) m² s⁻¹ is the value obtained for voids that form in the interior of the aluminum interconnects without surface contamination. 6.5 × 10⁻⁶ exp(−0.84 eV/kT) m² s⁻¹ is obtained for those voids that nucleate either at triple junctions or at the grain boundary-technical surface intersections, where the chemical impurities may act as trap centers for hopping vacancies.     

Superiority of PANS compared to LES in predicting a rudimentary landing gear flow with affordable meshes

Results of simulations of the flow around a rudimentary landing gear are presented in the paper. A newly proposed improved Partially-Averaged Navier–Stokes (PANS) method using k − ε − ζ − f turbulence model is used for prediction of the flow. The results are compared with the experimental data but also with the results of two LES simulations performed using the PANS computational grids. PANS simulations predicted the flow in good agreement with the experimental data. LES predicted a non-physical creation of separation over the front wheels that does not exist in the PANS prediction and was not observed in the experimental oil film. PANS simulations showed low sensitivity to the grid refinement. They show clear advantage compared with the LES simulations when the computational grid is inadequate for resolution of the near-wall flow structures.     

Flow of mono-dispersed particles through horizontal bend

Pipeline slurry flow of mono-dispersed particles through horizontal bend is numerically simulated by implementing Eulerian two-phase model in FLUENT software. A hexagonal shape and Cooper type non-uniform three-dimensional grid is chosen to discretize the entire computational domain, and a control volume finite difference method is used to solve the governing equations. The modeling results are compared with the experimental data collected in 53.0 mm diameter horizontal bend with radius of 148.4 mm for concentration profiles and pressure drops. Experiments are performed on narrow-sized silica sand with mean diameter of 450 μm and for flow velocity up to 3.56 m/s (namely, 1.78, 2.67 and 3.56 m/s) and four efflux concentrations up to 16.28% (namely, 0%, 3.94%, 8.82% and 16.28%) by volume for each velocity. Eulerian model gives fairly accurate predictions for both the pressure drop and concentration profiles at all efflux concentrations and flow velocities.     

Analytical modeling of synthetic fiber ropes subjected to axial loads. Part I: A new continuum model for multilayered fibrous structures

Synthetic fiber ropes are characterized by a very complex architecture and a hierarchical structure. Considering the fiber rope architecture, to pass from fiber to rope structure behavior, two scale transition models are necessary, used in sequence: one is devoted to an assembly of a large number of twisted components (multilayered), whereas the second is suitable for a structure with a central straight core and six helical wires (1 + 6). The part I of this paper first describes the development of a model for the static behavior of a fibrous structure with a large number of twisted components. Tests were then performed on two different structures subjected to axial loads. Using the model presented here the axial stiffness of the structures has been predicted and good agreement with measured values is obtained. A companion paper (Ghoreishi, S.R. et al., in press. Analytical modeling of synthetic fiber ropes, part II: A linear elastic model for 1 + 6 fibrous structures, International Journal of Solids and Structures, doi:10.1016/j.ijsolstr.2006.08.032) presents the second model to predict the mechanical behavior of a 1 + 6 fibrous structure.