Determination of the Infinitesimal Strain Tensor from Shears and Elongations

It is assumed that at a point P in a body, the longitudinal strains (elongations) along three non-coplanar directions are known from observation and that the shears of the three pairs of infinitesimal material line elements along the three non-coplanar directions are also known. With this information the strain tensor e at P is determined explicitly. The strain tensor e takes a simpler form in the special case when the three shears are zero. This simpler form is precisely the form obtained by Boulanger and Hayes in their study (Boulanger and Hayes, Proc R Ir Acad 103A:113–141, 2003) of the consequences of writing the displacement gradient at P as the sum of a skew symmetric tensor and a tensor with three real eigenvalues. The special case when the three elongations are zero is also considered. Dedicated to Franz Ziegler on the occasion of his seventieth birthday.     

On the stiffness of materials containing a disordered array of microscopic holes or hard inclusions

We study the macroscopic mechanical behavior of materials with microscopic holes or hard inclusions. Specifically, we deal with the effective elastic moduli of composites whose microgeometry consists of either soft or hard isolated inclusions surrounded by an elastic matrix. We approach this problem by taking the stiffness of the inclusion phase to be a complex variable, which we eventually evaluate at the soft or hard limits. Our main result states that there is a certain class of non-physical, negative-definite values of the elastic moduli of the inclusion phase for which the effective tensor does not have infinities or become otherwise singular.We present applications of this result to the estimation of effective moduli and to homogenization theorems. The first application involves using complexanalytic methods to obtain rigorous and accurate bounds on the effective moduli of the high-contrast composites under consideration. We also discuss the variational estimates of Rubenfeld & Keller, which yield a complementary set of bounds on these moduli. The best bounds are given by a combination of the analytical and variational results. As a second application, we show that certain known theorems of homogenization for materials with holes are simple consequences of our main result, and in this connection we establish corresponding new theorems for materials with hard inclusions. While our rederivation of the homogenization theorems for materials with holes can be closely related to other known constructions, it appears that certain elements provided by our main result are essential in the proof of homogenization for the hard-inclusion case.     

Strong-contrast expansion correlation approximations for the effective elastic moduli of multiphase composites

Following the previous approach of Pham and Torquato (J Appl Phys 94:6591–6602, 2003) and Torquato (J Mech Phys Solids 45:1421–1448, 1997; Random heterogeneous media, Springer, Berlin, 2002), we derive the strong-contrast expansions for the effective elastic moduli K _e,G _e of d-dimensional multiphase composites. The series consists of a principal reference part and a fluctuation part (perturbation about a homogeneous reference or comparison material), which contains multi-point correlation functions that characterize the microstructure of the composite. We propose a three-point correlation approximation for the fluctuation part with an objective choice of the reference phase moduli, such that the fluctuation terms vanish. That results in the approximations for the effective elastic moduli of isotropic composites, which coincide with the well-known self-consistent and Maxwell approximations for two-phase composites having respective microstructures. Applications to some two-phase materials are given.     

A classical problem revisited: rheology of nematic polymer monodomains in small amplitude oscillatory shear

We revisit the classical problem of the viscoelastic response of nematic (liquid crystal) polymers to small amplitude oscillatory shear. A multiple time scale perturbation analysis is applied to the Doi–Hess mesoscopic orientation tensor model to describe key features observed of longtime experiments, both physical (Moldenaers and Mewis, J Rheol, 30:567–584, 1986; Larson and Mead, J Rheol, 33:1251–1281, 1989b) and numerical (herein). First, there is a very slow time scale drift in the envelope of oscillations of the major director; we characterize the mean director angle and the envelope of oscillation. Second, there are bistable asymptotic orientational states, distinguished in that they are precisely the zero-stress orientational distributions noted in Larson and Mead (J Rheol, 33:185–206, 1989a). Third, the drift dynamics and asymptotic mean director angle are determined by the initial orientation of the director, not by material properties; we characterize the domain of attraction of each bistable state. Finally, the director drift leads to a predicted longtime decrease in the storage and loss moduli, consistent with experimental observations.

---

Identification of Damping Characteristic of a Structural Adhesive by Extended Modal Based Direct Model Updating Method

In order to identify damping characteristics (loss factors in bending and shear modes versus frequency) of an adhesive material, modal based direct model updating method used in previous paper (Jahani and Nobari, Exp Mech, 48:2008) is modified here. Two approaches adopted in the identification process. In the first one, the pre identified moduli from undamped model were used in identification of loss factors and in the second approach loss factors and moduli are identified simultaneously. The loss factors of a typical adhesive were identified using both approaches leading to acceptable results in both cases. The results show that loss factors and Young’s and shear moduli of adhesive are frequency dependent. All experiments and subsequent identifications are conducted both in bending and shear modes. The repeatability and consistency of method is proved by repeating the identification process for several times.     

A constitutive model for non-Newtonian materials based on the persistence-of-straining tensor

We present a constitutive equation for non-Newtonian materials which is capable of predicting, independently, steady state rheological material functions both in shear and in extension. The basic assumption is that the extra-stress tensor is a function of both the rate-of-strain tensor, D, and the persistence-of-straining tensor, -\boldsymbol{P}=\boldsymbol{D}\overline{\boldsymbol{W}}-\overline{\boldsymbol {W}}\boldsymbol{D}, introduced in Thompson and de Souza Mendes (Int. J. Eng. Sci. 43(1–2):79–105, 2005). The resulting equation falls within the category of constitutive equations of the form t=t(D,[`(W)])\boldsymbol{\tau}=\boldsymbol{\tau}(\boldsymbol {D},\overline{\boldsymbol{W}}), with the advantage of eliminating the undesirable stress jumps that may occur when [`(W)]\overline {\boldsymbol{W}} becomes locally undetermined. We also show that this formulation is not restricted to motions with constant relative principle stretch history (MWCRPSH), in contrast to what is suggested in the literature. The same basis of tensors that comes from representation theorems also arises from an elastic constitutive equation based on the difference between the Jauman and the Harnoy convected time derivatives, in the limit of small values of the Deborah number.     

Bounds on the non-local effective elastic properties of composites

The paper deals with the effective linear elastic behaviour of random media subjected to inhomogeneous mean fields. The effective constitutive laws are known to be non-local. Therefore, the effective elastic moduli show dispersion, i.e¹ they depend on the “wave vector” k of the mean field. In this paper the well-known Hashin-Shtrikman bounds (1962) for the Lamé parameters of isotropic multi-phase mixtures are generalized to inhomogeneous mean fields k ≠ 0. The bounds involve two-point correlations of random elastic moduli. In the limit k → ∞ the bounds converge to the exact result. The interest is focussed on composites with cell structures and on binary mixtures. To illustrate the results, numerical evaluations are carried out for a binary cell material composed of nearly spherical grains of equal size.     

A constrained theory for single crystal shape memory wires with application to restrained recovery

The theory of thin wires developed in Dret and Meunier (Comptes Rendus de l’Académie des Sciences. Série I. Mathématique 337:143–147, 2003) is adapted to phase-transforming materials with large elastic moduli in the sense discussed in James and Rizzoni (J Elast 59:399–436, 2000). The result is a one-dimensional constitutive model for shape memory wires, characterized by a small number of material constants. The model is used to analyze self-accommodated and detwinned microstructures and to study superelasticity. It also turns out that the model successfully reproduces the behavior of shape memory wires in experiments of restrained recovery (Tsoi et al. in Mater Sci Eng A 368:299–310, 2004; Tsoi in 50:3535–3544, 2002; S̆ittner et al. in Mater Sci Eng A 286:298–311, 2000; vokoun in Smart Mater Struct 12:680–685, 2003; Zheng and Cui in Intermetallics 12:1305–1309, 2004; Zheng et al. in J Mater Sci Technol 20(4):390–394, 2004). In particular, the model is able to predict the shift to higher transformation temperatures on heating. The model also captures the effect of prestraining on the evolution of the recovery stress and of the martensite volume fraction.     

The role of dissipation and defect energy in variational formulations of problems in strain-gradient plasticity. Part 2: single-crystal plasticity

Variational formulations are constructed for rate-independent problems in small-deformation single-crystal strain-gradient plasticity. The framework, based on that of Gurtin (J Mech Phys Solids 50: 5–32, 2002), makes use of the flow rule expressed in terms of the dissipation function. Provision is made for energetic and dissipative microstresses. Both recoverable and non-recoverable defect energies are incorporated into the variational framework. The recoverable energies include those that depend smoothly on the slip gradients, the Burgers tensor, or on the dislocation densities (Gurtin et al. J Mech Phys Solids 55:1853–1878, 2007), as well as an energy proposed by Ohno and Okumura (J Mech Phys Solids 55:1879–1898, 2007), which leads to excellent agreement with experimental results, and which is positively homogeneous and therefore not differentiable at zero slip gradient. Furthermore, the variational formulation accommodates a non-recoverable energy due to Ohno et al. (Int J Mod Phys B 22:5937–5942, 2008), which is also positively homogeneous, and a function of the accumulated dislocation density. Conditions for the existence and uniqueness of solutions are established for the various examples of defect energy, with or without the presence of hardening or slip resistance.     

Mean velocity and moments of turbulent velocity fluctuations in the wake of a model ship propulsor

Pod drives are modern outboard ship propulsion systems with a motor encapsulated in a watertight pod, whose shaft is connected directly to one or two propellers. The whole unit hangs from the stern of the ship and rotates azimuthally, thus providing thrust and steering without the need of a rudder. Force/momentum and phase-resolved laser Doppler anemometry (LDA) measurements were performed for in line co-rotating and contra-rotating propellers pod drive models. The measurements permitted to characterize these ship propulsion systems in terms of their hydrodynamic characteristics. The torque delivered to the propellers and the thrust of the system were measured for different operation conditions of the propellers. These measurements lead to the hydrodynamic optimization of the ship propulsion system. The parameters under focus revealed the influence of distance between propeller planes, propeller frequency of rotation ratio and type of propellers (co- or contra-rotating) on the overall efficiency of the system. Two of the ship propulsion systems under consideration were chosen, based on their hydrodynamic characteristics, for a detailed study of the swirling wake flow by means of laser Doppler anemometry. A two-component laser Doppler system was employed for the velocity measurements. A light barrier mounted on the axle of the rear propeller motor supplied a TTL signal to mark the beginning of each period, thus providing angle information for the LDA measurements. Measurements were conducted for four axial positions in the slipstream of the pod drive models. The results show that the wake of contra-rotating propeller is more homogeneous than when they co-rotate. In agreement with the results of the force/momentum measurements and with hypotheses put forward in the literature (see e.g. Poehls in Entwurfsgrundlagen für Schraubenpropeller, 1984; Schneekluth in Hydromechanik zum Schiffsentwurf, 1988; Breslin and Andersen in Hydrodynamics of ship propellers, 1996; Schneekluth and Bertram in Ship design for efficiency and economy, 1998), the co-rotating propellers model showed a much stronger swirl in the wake of the propulsor. The anisotropy of turbulence was analyzed using the anisotropy tensor introduced by Lumley and Newman (J Fluid Mech 82(1):161–178, 1977). The invariants of the anisotropy tensor of the wake flow were computed and were plotted in the Lumley–Newman-diagram. These measurements revealed that the anisotropy tensor in the wake of ship propellers is located near to the borders of the invariant map, showing a large degree of anisotropy. They will be presented and will be discussed with respect to applications of turbulence models to predict swirling flows.     

Bounds for effective elastic moduli of inhomogeneous solid bodies

In connection with the extensive use of various kinds of inhomogeneous materials (glass, carbon and boron reinforced plastics, cermets, concrete, reinforced materials, etc.) in technology, there arises a need to calculate the elastic properties of such systems. Here in each case it is necessary to work out specific methods for finding both elastic fields and effective moduli. Since, as a rule, such methods do not take into account the character of distribution of inhomogeneities in space, which is reflected on the form of the central moment functions [1], they can be referred to a single class and, consequently, can be obtained by a common method [2], In the given paper, by means of the method of solution of stochastic problems for microinhomogeneous solid bodies proposed in the work of the author [2], we find elastic fields and effective moduli in an arbitrary approximation. Depending on the choice of parameters, the latter form bounds within which there lie the exact values of the effective moduli. It is shown that the conditions used earlier for finding these parameters [3] are not the best ones. The effective elastic moduli of an inhomogeneous medium are calculated, and bounds, narrower than the bounds formed in [3], are found for them.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhniki, No. 5, pp. 144–150, September–October, 1973.     

Full 3D correlation tensor computed from double field stereoscopic PIV in a high Reynolds number turbulent boundary layer

The turbulence structure near a wall is a very active subject of research and a key to the understanding and modeling of this flow. Many researchers have worked on this subject since the fifties Hama et al. (J Appl Phys 28:388–394, 1957). One way to study this organization consists of computing the spatial two-point correlations. Stanislas et al. (C R Acad Sci Paris 327(2b):55–61, 1999) and Kahler (Exp Fluids 36:114–130, 2004) showed that double spatial correlations can be computed from stereoscopic particle image velocimetry (SPIV) fields and can lead to a better understanding of the turbulent flow organization. The limitation is that the correlation is only computed in the PIV plane. The idea of the present paper is to propose a new method based on a specific stereoscopic PIV experiment that allows the computation of the full 3D spatial correlation tensor. The results obtained are validated by comparison with 2D computation from SPIV. They are in very good agreement with the results of Ganapthisubramani et al. (J Fluid Mech 524:57–80, 2005a).     

Bounds for effective elastic moduli of disordered materials

Recently P.H. Dederichs and R. Zeller (1973) have developed a formal theory of the bounds of odd order n for the effective elastic moduli of linearly elastic, disordered materials. The bounds are established by use of statistical information given in terms of correlation functions up to order n (= 1, 3, 5,…). This theory is extended to include the bounds of even order n. It is indicated how these bounds can be made optimum under the given statistical information. The results for bounds of even and odd order are obtained in forms which resemble Neumann series, containing multiple integrals up to order (n?1). These integrals can be calculated for certain materials which are classified in terms of a gradual statistical homogeneity, isotropy and disorder. Materials which possess these properties up to the correlation functions of nth order are called overall grade n materials. The optimum bounds for overall grade 2 and grade 3 materials are given explicitly. Optimum bounds for materials which are of grade ∞ in homogeneity and isotropy (i.e. (statistically) perfectly homogeneous and isotropic) and, at the same time, disordered of grade 2 or 3 are also derived. Those for grade 2 in disorder are the Z. Hashin and S. Shtrikman's (1963) bounds. Those for grade 3 are the narrowest, explicit bounds so far derived for random elastic materials. They contain within themselves the so-called self-consistent elastic moduli.     

��Hypertractions and hyperstresses convey the same mechanical information Continuum Mech. Thermodyn. (2010) 22:163�C176�� by Prof. Podio Guidugli and Prof. Vianello and some related papers on higher gradient theories

In this commentary, we try to make clearer the state of the art concerning the relation between mechanical contact interactions and the different notions of stresses. We emphasize the importance of the concept of virtual displacements. Its role has been recognized in Mechanics and in Continuum Mechanics long ago (see e.g., Vailati in Il principio dei lavori virtuali da Aristotele a Erone d’Alessandria, 113–128, 1987; Russo in The forgotten revolution, Springer, Berlin, 2003, or Cosserat and Cosserat in Sur la Théorie des Corps Déformables, Herman, Paris, 1909; Cosserat and Cosserat in Note sur la théorie de l.action euclidienne, Gauthier-Villars, Paris, 1908), and it is central as well when starting with an expression of the power expended by internal stresses and deducing the form of contact interactions as when starting with some form of the contact interactions and developing a representation theorem for these contact interactions based on the Cauchy tetrahedron construction.     

Bounds on the elastic and transport properties of two-component composites

Weshow in detail how the McCoy bounds on the effective shear modulus of a statistically isotropic composite, can be simplified and expressed in terms of the volume fraction, f₁, and two geometric parameters. ζ₁ and η₁. We simplify Silnutzer's bounds on the effective elastic moduli of fibre-reinforced composites and find they can be expressed in terms f₁and two geometric parameters, ζ'₁ and η'₁. The parameter ζ'₁ also determines bounds on the transport and optical constants of such composites. Also, the Elsayed-McCoy bounds on the transport properties of fibre-reinforced, symmetric-cell materials are shown to depend on three geometric parameters.     

Comment on: I-Shih Liu: Constitutive theory of anisotropic rigid heat conductors

In I-Shih Liu’s paper (J Math Phys 50:083506, 2009), the compatibility of anisotropy and material frame indifference of a rigid heat conductor is investigated. For this purpose, the deformation gradient is introduced into the domain of the constitutive mapping. Because of the presupposed rigidity, the deformation gradient is here represented by an orthogonal tensor. The statement, that the usual procedure—not to introduce the deformation gradient into the state space of rigid heat conductors—causes isotropy because of material frame indifference, is misleading.     

On the instability of steady motion

This paper deals with the instability of steady motions of conservative mechanical systems with cyclic coordinates. The following are applied: Kozlov’s generalization of the first Lyapunov’s method, as well as Rout’s method of ignoration of cyclic coordinates. Having obtained through analysis the Maclaurin’s series for the coefficients of the metric tensor, a theorem on instability is formulated which, together with the theorem formulated in Furta (J. Appl. Math. Mech. 50(6):938–944, 1986), contributes to solving the problem of inversion of the Lagrange-Dirichlet theorem for steady motions. The cases in which truncated equations involve the gyroscopic forces are solved, too. The algebraic equations resulting from Kozlov’s generalizations of the first Lyapunov’s method are formulated in a form including one variable less than was the case in existing literature.     

Optimal anisotropic three-phase conducting composites: Plane problem

The paper establishes tight lower bound for effective conductivity tensor K₁ of two-dimensional three-phase conducting anisotropic composites and defines optimal microstructures. It is assumed that three materials are mixed with fixed volume fractions and that the conductivity of one of the materials is infinite. The bound expands the Hashin–Shtrikman and translation bounds to multiphase structures, it is derived using a combination of translation method and additional inequalities on the fields in the materials; similar technique was used by Nesi, 1995, Cherkaev, 2009 for isotropic multiphase composites. This paper expands the bounds to the anisotropic composites with effective conductivity tensor K₁. The lower bound of conductivity (G-closure) is a piece-wise analytic function of eigenvalues of K₁, that depends only on conductivities of components and their volume fractions. Also, we find optimal microstructures that realize the bounds, developing the technique suggested earlier by Albin et al., 2007a, Cherkaev, 2009. The optimal microstructures are laminates of some rank for all regions. The found structures match the bounds in all but one region of parameters; we discuss the reason for the gap and numerically estimate it.