Thermal buckling of circular cylindrical and conical thin-walled shells

"Aeronautical Research Laboratory. Contract AF 33(616)-3885, Project no. 7(8-1367), Task no. 70524."

On the buckling of circular cylindrical shells under external pressure,

"Project no. 9782. Contract AF 49(638)-430... CU-23-62-AF-430-CE."

Dynamic behavior of cylindrical shells under initial stress,

"Project no. 9782. Contract AF 49 (638)-430 ... CU-11-60-AF-430-CE."

Nonlinear theories for thin shells,

"Contract Nonr 1866(02)."

On the possibility of elastic strain localisation in a fault

The phenomenon of strain localisation is often observed in shear deformation of particulate materials, e.g., fault gouge. This phenomenon is usually attributed to special types of plastic behaviour of the material (e.g., strain softening or mismatch between dilatancy and pressure sensitivity or both). Observations of strain localisation in situ or in experiments are usually based on displacement measurements and subsequent computation of the displacement gradient. While in conventional continua the symmetric part of the displacement gradient is equal to the strain, it is no longer the case in the more realistic descriptions within the framework of generalised continua. In such models the rotations of the gouge particles are considered as independent degrees of freedom the values of which usually differ from the rotation of an infinitesimal volume element of the continuum, the latter being described for infinitesimal deformations by the non-symmetric part of the displacement gradient. As a model for gouge material we propose a continuum description for an assembly of spherical particles of equal radius in which the particle rotation is treated as an independent degree of freedom. Based on this model we consider simple shear deformations of the fault gouge. We show that there exist values of the model parameters for which the displacement gradient exhibits a pronounced localisation at the mid-layers of the fault, even in the absence of inelasticity. Inelastic effects are neglected in order to highlight the role of the independent rotations and the associated additional parameters. The localisation-like behaviour occurs if (a) the particle rotations on the boundary of the shear layer are constrained (this type of boundary condition does not exist in a standard continuum) and (b) the contact moment-or bending stiffness is much smaller than the product of the effective shear modulus of the granulate and the square of the width of the gouge layer. It should be noted however that the virtual work functional is positive definite over the range of physically meaningful parameters (here: contact stiffnesses, solid volume fraction and coordination number) so that strictly speaking we are not dealing with a material instability.

Influence of imperfect interfaces on bending and vibration of laminated composite shells

This paper is devoted to modeling elastic behavior of laminated composite shells, with special emphasis on incorporating interfacial imperfection. The conditions of imposing traction continuity and displacement jump across each interface are used to model imperfect interfaces. Vanishing transverse shear stresses on two free surfaces of a shell eliminate the need for shear correction factors. A linear theory underlying elastostatics and kinetics of laminated composite shells in a general configuration is presented from Hamilton's principle. In the special case of vanishing interfacial parameters, this theory reduces to the conventional third-order zigzag theory for perfectly bonded laminated shells. Numerical results for bending and vibration problems of laminated circular cylindrical panels are tabulated and plotted to indicate the influence of the interfacial imperfection. (C) 2000 Elsevier Science Ltd. All rights reserved.

Bistable attractors in a model of convection-driven spherical dynamos

The range of existence and the properties of two essentially different chaotic attractors found in a model of nonlinear convection-driven dynamos in rotating spherical shells are investigated. A hysteretic transition between these attractors is established as a function of the rotation parameter t. The width of the basins of attraction is also estimated. © 2012 The Royal Swedish Academy of Sciences.

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.

Phase-field analysis of finite-strain plates and shells including element subdivision

With the theme of fracture of finite-strain plates and shells based on a phase-field model of crack regularization, we introduce a new staggered algorithm for elastic and elasto-plastic materials. To account for correct fracture behavior in bending, two independent phase-fields are used, corresponding to the lower and upper faces of the shell. This is shown to provide a realistic behavior in bending-dominated problems, here illustrated in classical beam and plate problems. Finite strain behavior for both elastic and elasto-plastic constitutive laws is made compatible with the phase-field model by use of a consistent updated-Lagrangian algorithm. To guarantee sufficient resolution in the definition of the crack paths, a local remeshing algorithm based on the phase- field values at the lower and upper shell faces is introduced. In this local remeshing algorithm, two stages are used: edge-based element subdivision and node repositioning. Five representative numerical examples are shown, consisting of a bi-clamped beam, two versions of a square plate, the Keesecker pressurized cylinder problem, the Hexcan problem and the Muscat-Fenech and Atkins plate. All problems were successfully solved and the proposed solution was found to be robust and efficient.

Aplicação de indentação instrumentada na caracterização mecânica de poliuretana derivada de óleo de mamona

Neste trabalho são investigadas as propriedades mecânicas de poliuretana derivada do óleo de mamona, utilizando a técnica de indentação instrumentada com penetradores de geometrias piramidal e esférica. Foi analisada a influência da forma do penetrador utilizado nos ensaios de indentação instrumentada para se obter valores das propriedades mecânicas de polímero derivado de óleo de mamona. Os penetradores utilizados são de pontas piramidais dos tipos Berkovich e canto de cubo e esférico de raio igual a 150 μm em um Nanoindenter XP TM com cargas aplicadas entre 1 e 200 mN. As penetrações variam de acordo com o formato do penetrador, sendo maiores para pontas agudas. A dureza e o módulo de elasticidade foram determinados, utilizando o método de Oliver e Pharr. Verificou-se que os valores medidos para a dureza são maiores para penetradores mais agudos. Os valores obtidos com a ponta piramidal Berkovich foram de 0,14 GPa para pequenas penetrações e 0,12 GPa para maiores penetrações. Já os valores obtidos com ponta canto de cubo foram 25 a 30% maiores. Isso está relacionado com os volumes das regiões que apresentam deformações plásticas elevadas, no caso de penetradores agudos comparados com os volumes das regiões que sofrem deformações viscoelásticas. A viscosidade aparente determinada, utilizando penetrador esférico em testes de força aplicada constante, é igual a (22 ± 2) × 10(12) Pa.s.

A positional FEM Formulation for geometrical non-linear analysis of shells

This work presents a fully non-linear finite element formulation for shell analysis comprising linear strain variation along the thickness of the shell and geometrically exact description for curved triangular elements. The developed formulation assumes positions and generalized unconstrained vectors as the variables of the problem, not displacements and finite rotations. The full 3D Saint-Venant-Kirchhoff constitutive relation is adopted and, to avoid locking, the rate of thickness variation enhancement is introduced. As a consequence, the second Piola-Kirchhoff stress tensor and the Green strain measure are employed to derive the specific strain energy potential. Curved triangular elements with cubic approximation are adopted using simple notation. Selected numerical simulations illustrate and confirm the objectivity, accuracy, path independence and applicability of the proposed technique.

Analysis of the tip roundness effects on the micro- and macroindentation response of elastic-plastic materials

In this work, the effects of indenter tip roundness oil the load-depth indentation curves were analyzed using finite element modeling. The tip roundness level was Studied based on the ratio between tip radius and maximum penetration depth (R/h(max)), which varied from 0.02 to 1. The proportional Curvature constant (C), the exponent of depth during loading (alpha), the initial unloading slope (S), the correction factor (beta), the level of piling-up or sinking-in (h(c)/h(max)), and the ratio h(max)/h(f) are shown to be strongly influenced by the ratio R/h(max). The hardness (H) was found to be independent of R/h(max) in the range studied. The Oliver and Pharr method was successful in following the variation of h(c)/h(max) with the ratio R/h(max) through the variation of S with the ratio R/h(max). However, this work confirmed the differences between the hardness values calculated using the Oliver-Pharr method and those obtained directly from finite element calculations; differences which derive from the error in area calculation that Occurs when given combinations of indented material properties are present. The ratio of plastic work to total work (W(p)/W(t)) was found to be independent of the ratio R/h(max), which demonstrates that the methods for the Calculation of mechanical properties based on the *indentation energy are potentially not Susceptible to errors caused by tip roundness.

Analysis of the effects of conical indentation variables on the indentation response of elastic-plastic materials through factorial design of experiment

In this work, the effects of conical indentation variables on the load-depth indentation curves were analyzed using finite element modeling and dimensional analysis. A factorial design 2(6) was used with the aim of quantifying the effects of the mechanical properties of the indented material and of the indenter geometry. Analysis was based on the input variables Y/E, R/h(max), n, theta, E, and h(max). The dimensional variables E and h(max) were used such that each value of dimensionless Y/E was obtained with two different values of E and each value of dimensionless R/h(max) was obtained with two different h(max) values. A set of dimensionless functions was defined to analyze the effect of the input variables: Pi(1) = P(1)/Eh(2), Pi(2) = h(c)/h, Pi(3) = H/Y, Pi(4) = S/Eh(max), Pi(6) = h(max)/h(f) and Pi(7) = W(P)/W(T). These six functions were found to depend only on the dimensionless variables studied (Y/E, R/h(max), n, theta). Another dimension less function, Pi(5) = beta, was not well defined for most of the dimensionless variables and the only variable that provided a significant effect on beta was theta. However, beta showed a strong dependence on the fraction of the data selected to fit the unloading curve, which means that beta is especially Susceptible to the error in the Calculation of the initial unloading slope.

Using the ultrasound and instrumented indentation techniques to measure the elastic modulus of engineering materials

Currently, the acoustic and nanoindentation techniques are two of the most used techniques for material elastic modulus measurement. In this article fundamental principles and limitations of both techniques are shown and discussed. Last advances in nanoindentation technique are also reviewed. An experimental study in ceramic, metallic, composite and single crystals was also done. Results shown that ultrasonic technique is capable to provide results in agreement with those reported in literature. However, ultrasonic technique does not allow measuring the elastic modulus of some small samples and single crystals. On the other hand, the nanoindentation technique estimates the elastic modulus values in reasonable agreement with those measured by acoustic methods, particularly in amorphous materials, while in some policristaline materials some deviation from expected values was obtained.

Quasiharmonic elastic constants corrected for deviatoric thermal stresses

The quasiharmonic approximation (QHA), in its simplest form also called the statically constrained (SC) QHA, has been shown to be a straightforward method to compute thermoelastic properties of crystals. Recently we showed that for noncubic solids SC-QHA calculations develop deviatoric thermal stresses at high temperatures. Relaxation of these stresses leads to a series of corrections to the free energy that may be taken to any desired order, up to self-consistency. Here we show how to correct the elastic constants obtained using the SC-QHA. We exemplify the procedure by correcting to first order the elastic constants of MgSiO(3) perovskite and MgSiO(3) postperovskite, the major phases of the Earth's lower mantle. We show that this first-order correction is quite satisfactory for obtaining the aggregated elastic averages of these minerals and their velocities in the lower mantle. This type of correction is also shown to be applicable to experimental measurements of elastic constants in situations where deviatoric stresses can develop, such as in diamond-anvil cells.