Micromechanical definition of an entropy for quasi-static deformation of granular materials

A micromechanical theory is formulated for quasi-static deformation of granular materials, which is based on information theory. A reasoning is presented that leads to the definition of an information entropy that is appropriate for quasi-static deformation of granular materials. This definition is based on the hypothesis that relative displacements at contacts with similar orientations are independent realisations of a random variable. This hypothesis is made plausible based on the results of Discrete Element simulations. The developed theory is then used to predict the elastic behaviour of granular materials in terms of micromechanical quantities. The case considered is that of two-dimensional assemblies consisting of non-rotating particles with an elastic contact constitutive relation. Applications of this case are the initial elastic (small-strain) deformation of granular materials. Theoretical results for the elastic moduli, relative displacements, energy distribution and probability density functions are compared with results obtained from the Discrete Element simulations for isotropic assemblies with various average numbers of contacts per particle and various ratios of tangential to normal contact stiffness. This comparison shows that the developed information theory is valid for loose systems, while a theory based on the uniform-strain assumption is appropriate for dense systems.     

The elastic response of a cohesive aggregate—a discrete element model with coupled particle interaction

A model is presented for the deformation of a cohesive aggregate of elastic particles that incorporates two important effects of large-sized inter-particle junctions. A finite element model is used to derive a particle response rule, for both normal and tangential relative deformations between pairs of particles. This model agrees with the Hertzian contact theory for small junctions, and is valid for junctions as large as half the nominal particle size. Further, the aggregate model uses elastic superposition to account for the coupled force–displacement response due to the simultaneous displacement of all of the neighbors of each particle in the aggregate. A particle stiffness matrix is developed, relating the forces at each junction to the three displacement degrees of freedom at all of the neighboring-particle junctions. The particle response satisfies force and moment equilibrium, so that the model is properly posed to allow for rigid rotation of the particle without introducing rotational degrees of freedom. A computer-simulated sintering algorithm is used to generate a random particle packing, and the stiffness matrix is derived for each particle. The effective elastic response is then estimated using a mean field or affine displacement calculation, and is also found exactly by a discrete element model, solving for the equilibrium response of the aggregate to uniform-strain boundary conditions. Both the estimate and the exact solution compare favorably with experimental data for the bulk modulus of sintered alumina, whereas Hertzian contact-based models underestimate the modulus significantly. Poisson's ratio is, however, accurately determined only by the full equilibrium discrete element solution, and shown to depend significantly on whether or not rigid particle rotation is permitted in the model. Moreover, this discrete element model is sufficiently robust, so it can be applied to problems involving non-homogeneous deformations in such cohesive aggregates.     

On the elastic moduli of three-dimensional assemblies of spheres: Characterization and modeling of fluctuations in the particle displacement and rotation

The elastic moduli of four numerical random isotropic packings of Hertzian spheres are studied. The four samples are assembled with different preparation procedures, two of which aim to reproduce experimental compaction by vibration and lubrication. The mechanical properties of the samples are found to change with the preparation history, and to depend much more on coordination number than on density.Secondly, the fluctuations in the particle displacements from the average strain are analyzed, and the way they affect the macroscopic behavior analyzed. It is found that only the average over equally oriented contacts of the relative displacement these fluctuations induce is relevant at the macroscopic scale. This average depends on coordination number, average geometry of the contact network and average contact stiffness. As far as the separate contributions from particle displacements and rotations are concerned, the former is found to counteract the average strain along the contact normal, while the latter do in the tangential plane. Conversely, the tangential components of the center displacements mainly arise to enforce local equilibrium, and have a small and generally stiffening effect at the macro-scale.Finally, the fluctuations and the shear modulus that result from two approaches available in the literature are estimated numerically. These approaches are both based on the equilibrium of a small-sized representative assembly. The improvement of these estimate with respect to the average strain assumption indicates that the fluctuations relevant to the macroscopic behavior occur with short correlation length.     

Effective elastic properties of periodic composite medium

A new method is presented for calculating the bulk effective elastic stiffness tensor of a two-component composite with a periodic microstructure. The basic features of this method are similar to the one introduced by Bergman and Dunn (1992) for the dielectric problem. It is based on a Fourier representation of an integro-differential equation for the displacement field, which is used to produce a continued-fraction expansion for the elastic moduli. The method enabled us to include a much larger number of Fourier components than some previously proposed Fourier methods. Consequently our method provides the possibility of performing reliable calculations of the effective elastic tensor of periodic composites that are neither dilute nor low contrast, and are not restricted to arrays of nonoverlapping inclusions. We present results for a cubic array of nonoverlapping spheres, intended to serve as a test of quality, as well as results for a cubic array of overlapping spheres and a two dimensional hexagonal array of circles (a model for a fiber reinforced material) for comparison with previous work.     

Effective elastic modulus of film-on-substrate systems under normal and tangential contact

Load and depth sensing indentation methods have been widely used to characterize the mechanical properties of the thin film-substrate systems. The measurement accuracy critically depends on our knowledge of the effective elastic modulus of this heterogeneous system. In this work, based on the exact solution of the Green's function in Fourier space, we have derived an analytical relationship between the surface tractions and displacements, which depends on the ratio of the film thickness to contact size and the generalized Dundurs parameters that describe the modulus mismatch between the film and substrate materials. The use of the cumulative superposition method shows that the contact stiffness of any axisymmetric contact is the same as that of a flat-ended punch contact. Therefore, assuming a surface traction of the form of [1−(r/a)²]^−1/2 with radial coordinate r and contact size a, we can obtain an approximate representation of the effective elastic moduli, which agree extremely well with the finite element simulations for both normal and tangential contacts. Motivated by a recently developed multidimensional nanocontact system, we also explore the dependence of the ratio of tangential to normal contact stiffness on the ratio of film thickness to contact radius and the Dundurs parameters. The analytical representations of the correction factors in the relationship between the contact stiffness and effective modulus are derived at infinite friction conditions.     

Displacement of surrounding rock in a deep circular hole considering double moduli and strength-stiffness degradation

The problem of cavity stability widely exists in deep underground engineering and energy exploitation. First, the stress field of the surrounding rock under the uniform stress field is deduced based on a post-peak strength drop model considering the rock's characteristics of constant modulus and double moduli. Then, the orthogonal non-associative flow rule is used to establish the displacement of the surrounding rock under constant modulus and double moduli, respectively, considering the stiffness degradation and dilatancy effects in the plastic region and assuming that the elastic strain in the plastic region satisfies the elastic constitutive relationship. Finally, the evolution of the displacement in the surrounding rock is analyzed under the effects of the double modulus characteristics, the strength drop, the stiffness degradation, and the dilatancy. The results show that the displacement solutions of the surrounding rock under constant modulus and double moduli have a unified expression. The coefficients of the expression are related to the stress field of the original rock, the elastic constant of the surrounding rock, the strength parameters, and the dilatancy angle. The strength drop, the stiffness degradation, and the dilatancy effects all have effects on the displacement. The effects can be characterized by quantitative relationships.     

Bounds on the Elastic Moduli of Completely Random Two-Dimensional Polycrystals

Explicit bounds on the elastic moduli of completely random planar polycrystals, the shape and crystalline orientations of the constituent grains of which are uncorrelated, are derived and calculated for a number of crystals of general two-dimensional anisotropy. The bounds on the elastic two-dimensional bulk modulus happen to coincide with the simple third order (in anisotropy contrast) bounds for the subclass of idealistic circular cell polycrystals. The bounds on the shear modulus are close to the much simpler bounds for circular cell polycrystals, which approximate aggregates of equiaxed grains.     

Tailored heterogeneity increases overall stability regime of composites having a negative-stiffness inclusion

Recent theoretical and experimental results have shown the possibility of enormous increases in composite material overall elastic stiffness, damping, thermal expansion, piezoelectricity, etc., when the composite contains a tuned non-positive-definite (i.e., negative stiffness) constituent. For such composite materials to have practical utility, they must be stable. Recent research has shown they can be, for a limited range of constituent negative stiffness. This research has treated linear elastic composite materials with homogeneous phases, via the energy method and full dynamic stability analyses.In the present work, we first show how to analyze the composites previously treated by the comprehensive but simpler static stability approach, obtaining closed-form results. We then employ this approach to show that permitting heterogeneity of the positive-definite phase can substantially increase the range of constituent negative stiffness while maintaining overall composite stability. We first treat the positive-definite phase heterogeneity as piecewise homogeneous, and then treat it as continuously-varying. In the continuously-varying heterogeneity case, we seek the radially optimal distribution of the elastic moduli in the coatings, under constant coating average moduli constraint, to permit the most negative possible inclusion stiffness while maintaining overall composite stability. This is accomplished for three coating cases: constant bulk modulus but arbitrarily radially-varying shear modulus; constant shear modulus but arbitrarily radially-varying bulk modulus; and both moduli arbitrarily radially varying. We find the optimal coatings to be: a heterogeneous one with shear modulus being a specific continuously decreasing function of radius for the first case; a homogeneous one for the second case; and a heterogeneous one with both moduli being either Dirac-delta or Heaviside-step decreasing functions of radius for the last case (if the coating moduli are unrestricted in magnitude or have upper limits, respectively). The results show a substantial increase in the permissible inclusion negative stiffness range is provided by coating heterogeneity, while maintaining overall composite stability. Such an increased range of constituent negative stiffness provides an enlarged tuning parameter range for the development of novel, high-performance composite materials.     

Thin interphase/imperfect interface in elasticity with application to coated fiber composites

The imperfect interface conditions which are equivalent to the effect of a thin elastic interphase are derived by a Taylor expansion method in terms of interface displacement and traction jumps. Plane and cylindrical interfaces are analyzed as special cases. The effective elastic moduli of a unidirectional coated fiber composite are obtained on the basis of the derived imperfect interface conditions. High accuracy of the method is demonstrated by comparison of solutions of several problems in terms of the imperfect interface conditions or explicit presence of interphase as a third phase. The problems considered are transverse shear of a coated infinite fiber in infinite matrix and effective transverse bulk and shear moduli and effective axial shear modulus of a coated fiber composite. Unlike previous elastic imperfect interface conditions in the literature, the present ones are valid for the entire range of interphase stiffness, from very small to very large.     

ELASTIC BEHAVIOR ANALYSIS OF 3D ANGLE-INTERLOCK WOVEN CERAMIC COMPOSITES

A micromechanical model for elastic behavior analysis of angle-interlock woven ceramic composites is proposed in this paper. This model takes into account the actual fabric structure by considering the fiber undulation and continuity in space, the cavities between adjacent yarns and the actual cross-section geometry of the yarn. Based on the laminate theory, the elastic properties of 3D angle-interlock woven ceramic composites are predicted. Different numbers of interlaced wefts have almost the same elastic moduli. The thickness of ceramic matrix has little effect on elastic moduli. When the undulation ratio increases longitudinal modulus decreases and the other Young's moduli increase. Good agreement between theoretical predictions and experimental results demonstrates the feasibility of the proposed model in analyzing the elastic properties of 3D angle-interlock woven ceramic composites. The results of this paper verify the fact that the method of analyzing polyester matrix composites is suitable for woven ceramic composites.     

Micromechanical study of elastic moduli of three-dimensional granular assemblies

In micromechanics of granular materials, relationships are investigated between micro-scale characteristics of particles and contacts and macro-scale continuum characteristics. For three-dimensional isotropic assemblies the macro-scale elastic characteristics are described by the bulk and the shear modulus, which depend on the micro-scale characteristics of the coordination number (i.e. the average number of contacts per particle) and the interparticle contact stiffnesses in directions normal and tangential to the contact.     

基于敏度分析的拉压不同模量桁架问题的数值分析

利用光滑函数技术,提出光滑化的拉压不同弹性模量问题的应力应变关系,与有限元方法相结合,建立了拉压不同模量一维连续体与桁架结构的数值求解模型,推导了敏度计算公式,采用Newton-Raphson算法进行求解.数值结果表明,本文算法具有较高的计算精度和收敛速度.     

Analytical test functions for free vibration analysis of rotating non-homogeneous Timoshenko beams

In this paper, the governing equations for free vibration of a non-homogeneous rotating Timoshenko beam, having uniform cross-section, is studied using an inverse problem approach, for both cantilever and pinned-free boundary conditions. The bending displacement and the rotation due to bending are assumed to be simple polynomials which satisfy all four boundary conditions. It is found that for certain polynomial variations of the material mass density, elastic modulus and shear modulus, along the length of the beam, the assumed polynomials serve as simple closed form solutions to the coupled second order governing differential equations with variable coefficients. It is found that there are an infinite number of analytical polynomial functions possible for material mass density, shear modulus and elastic modulus distributions, which share the same frequency and mode shape for a particular mode. The derived results are intended to serve as benchmark solutions for testing approximate or numerical methods used for the vibration analysis of rotating non-homogeneous Timoshenko beams.     

The mechanics of pyramids

The elastostatic fields and free-vibration characteristics are studied for a class of continuous solids in the shape of homogeneous or layered pentahedral pyramids. The pyramids possess an arbitrary number of linear elastic layers containing dissimilar elastic constants. The layers are assumed to be perfectly bonded in the out-of-plane direction of the pyramid. The change in elastic stiffness across each interface requires a model that allows for a jump in displacement gradient at these location. A discrete-layer representation is used that combines the Ritz method with polynomial in-plane approximations with one-dimensional Lagrangian interpolation polynomials in the thickness direction. The free-vibration characteristics are examined for a variety of isotropic and anisotropic materials and representative bulk stiffness estimates are given for homogeneous pyramids under static deformation.     

Mechanical properties of single-walled carbon nanotube bundles as bulk materials

Single-walled carbon nanotubes (SWNTs) in crystalline bundles may exhibit a transition in which the cross-sections of tubes turn from perfectly circular to hexagonal, depending upon the tube diameter and externally applied pressure, and this structural instability leads to an abrupt change in the bulk elastic properties of SWNT bundles. This paper presents a hybrid atom/continuum model to study the bulk elastic properties of SWNT bundles, and the predicted characteristics of this structural instability agree well with the experimental observations available in the literature. Linearized bulk elastic properties of SWNT bundles with respect to a stable configuration are transversely isotropic and hence can be characterized by five independent elastic moduli. A complete set of these five moduli is predicted for the first time. It is found that the deformability of tube cross-sections play a dominant role in characterizing the transverse moduli.     

Spring-network model of red blood cell: From membrane mechanics to validation

Red blood cell membrane is highly elastic and proper modeling of this elasticity is essential for biomedical applications that involve computational experiments with blood flow. Inseparable and often some of the most difficult parts of modeling process are verification and validation. In this work, we present a revised model, which uses a spring network to represent the cell membrane immersed in a fluid and has been successfully used in blood flow simulations. We demonstrate the validation steps by first deriving the theoretical relations between the bulk properties of elastic membranes—shear modulus and area compressibility modulus—and parameters of the model that enter the nonlinear stretching and local area conservation computational moduli. We verify the theoretically derived relations using computer simulations of deformable triangular mesh. We calibrate the model by performing a computational version of the optical tweezers experiment. And finally, we validate the modeled cell behavior by investigating the cell rotation frequency when it is subjected to shear flow and cell deformation in narrow channels. The supplementary material contains an extensive dataset that can be used for setting different elastic properties for each cell in simulations of dense suspensions, while still conforming to the biological data. This work contains a complete model development process: From modelling of basic mechanical concepts (the spring network) and advanced biomechanical concepts (such as elasticity of the membrane), through calibration process towards the final stage of model validation.     

复合材料模量细观分析的一般性模型

基于细观力学理论,提出了横观各向异性复合材料弹性常数的分析模型,并给出了纵向(纤维长度方向) 与横向弹性模量、横截面内剪切模量和泊松比等材料性能参数的理论推导和计算公式. 此方法简洁准确地反映了材料横观各向异性性能,改进了常规细观力学模型,也为工程分析提供了理论依据;同时在复合材料力学课程教学中,有助于深化学生对复合材料细观分析的理解,具有一定的启发性.     

A new micromechanical approach of micropolar continuum modeling for 2-D periodic cellular material

In this paper, we present a new united approach to formulate the equivalent micropolar constitutive relation of two-dimensional (2-D) periodic cellular material to capture its non-local properties and to explain the size effects in its structural analysis. The new united approach takes both the displacement compatibility and the equilibrium of forces and moments into consideration, where Taylor series expansion of the displacement and rotation fields and the extended aver-aging procedure with an explicit enforcement of equilibrium are adopted in the micromechanical analysis of a unit cell. In numerical examples, the effective micropolar constants obtained in this paper and others derived in the literature are used for the equivalent micropolar continuum simulation of cellular solids. The solutions from the equivalent analysis are compared with the discrete simulation solutions of the cellu-lar solids. It is found that the micropolar constants developed in this paper give satisfying results of equivalent analysis for the periodic cellular material.