Homogenization of acoustic waves in strongly heterogeneous porous structures

We consider acoustic waves in fluid-saturated periodic media with dual porosity. At the mesoscopic level, the fluid motion is governed by the Darcy flow model extended by inertia terms and by the mass conservation equation. In this study, assuming the porous skeleton is rigid, the aim is to distinguish the effects of the strong heterogeneity in the permeability coefficients. Using the asymptotic homogenization method we derive macroscopic equations and obtain the dispersion relationship for harmonic waves. The double porosity gives rise to an extra homogenized coefficient of dynamic compressibility which is not obtained in the upscaled single porosity model. Both the single and double porosity models are compared using an example illustrating wave propagation in layered media.     

均匀化方法在粘弹性多层复合材料中的应用

主要研究了由各向同性线弹性加强体和各向同性线粘弹性基体组成的多层复合材料的问题,在已有的线弹性多层材料的均匀化方法的基础上,应用弹性一粘弹性对应原理,在Carson域中求解粘弹性多层材料的问题。通过Burgers模型表示线粘弹性基体材料,反演得到了多层材料的有效松弛模量和有效泊松比在时间域中的表达式,并且与实验结果和其他结果进行了比较。     

Second-order homogenization of periodic materials based on asymptotic approximation of the strain energy: formulation and validity limits

In this paper a second-order homogenization approach for periodic material is derived from an appropriate representation of the down-scaling that correlates the micro-displacement field to the macro-displacement field and the macro-strain tensors involving unknown perturbation functions. These functions take into account of the effects of the heterogeneities and are obtained by the solution of properly defined recursive cell problems. Moreover, the perturbation functions and therefore the micro-displacement fields result to be sufficiently regular to guarantee the anti-periodicity of the traction on the periodic unit cell. A generalization of the macro-homogeneity condition is obtained through an asymptotic expansion of the mean strain energy at the micro-scale in terms of the microstructural characteristic size ?; the obtained overall elastic moduli result to be not affected by the choice of periodic cell. The coupling between the macro- and micro-stress tensor in the periodic cell is deduced from an application of the generalised macro-homogeneity condition applied to a representative portion of the heterogeneous material (cluster of periodic cell). The correlation between the proposed asymptotic homogenization approach and the computational second-order homogenization methods (which are based on the so called quadratic ansätze) is obtained through an approximation of the macro-displacement field based on a second-order Taylor expansion. The form of the overall elastic moduli obtained through the two homogenization approaches, here proposed, is analyzed and the differences are highlighted. An evaluation of the developed method in comparison with other recently proposed in literature is carried out in the example where a three-phase orthotropic material is considered. The characteristic lengths of the second-order equivalent continuum are obtained by both the asymptotic and the computational procedures here analyzed. The reliability of the proposed approach is evaluated for the case of shear and extensional deformation of the considered two-dimensional infinite elastic medium subjected to periodic body forces; the results from the second-order model are compared with those of the heterogeneous continuum.     

复合材料中的渐近均匀化方法

本文将非均质弹性体的渐近均匀化方法应用于复合材料的宏观与细观分析之中。该方法基于平均化的思想，将复合材料视作由周期性的细观结构所构成，其场变量依赖于宏观和细观两个尺度的坐标变量而变化。通过建立位移和应力的渐近表达式，推导出关于周期性基元的细观平衡方程和细观本构关系，并与有限元数值方法相结合，得到材料的宏观等效性能和细观应力分布。对典型算例的分析，反映出该方法的有效性及准确性。     

Variational asymptotic micromechanics modeling of heterogeneous porous materials

This paper presents a numerical technique to predict the effective elastic properties of heterogeneous fluid-filled porous media where the heterogeneity may result from dissimilar solid and fluid phase properties or due to mismatch in porous microstructure. The technique is based on the variational asymptotic method of homogenization where finite element method is employed for discretization. Biot’s theory of poroelasticity is used to describe porous media where both solid and fluid phase motions (u ? U formulation) are considered with associated strain measures. The method estimates the poroelastic constitutive law in single analysis which makes it very efficient compared to other finite element based homogenization techniques. The method is also general enough to compute all 28 elements of an anisotropic constitutive matrix. Other than estimating the effective properties the micro-stress/strain distribution is also obtained at no additional cost.The method is successfully applied for homogenization of porous media, fluid-filled cavity and finally for effective property estimation of bone lamella. In absence of any other direct method of porous media homogenization, the present technique is compared with classical homogenization methods with fluid approximated as solid of very high Poisson’s ratio. The suitability of this approximation and various other alternatives are also discussed. It is shown that the present homogenization method can be an efficient tool for bone property estimation where fluid-filled porous hierarchical micro-/nanostructure must be respected at all steps.     

Derivation of the Forchheimer Law via Homogenization

In this paper we derive the Forchheimer law via the theory of homogenization. In particular, we study the nonlinear correction to Darcy's law due to inertial effects on the flow of a Newtonian fluid in rigid porous media. A general formula for this correction term is derived directly from the Navier–Stokes equation via homogenization. Unlike other studies based on the same approach that concluded for the nonlinear correction to be cubic in velocity for isotropic media, the present work shows that the nonlinear correction is quadratic. An example is constructed to illustrate our theory. In this example, the analytic solution to the Navier–Stokes equation is obtained and is utilized to show the validity of the quadratic correction. Both incompressible and compressible fluids are considered.     

Novel numerical implementation of asymptotic homogenization method for periodic plate structures

The present paper develops and implements finite element formulation for the asymptotic homogenization theory for periodic composite plate and shell structures, earlier developed in  and , and thus adopts this analytical method for the analysis of periodic inhomogeneous plates and shells with more complicated periodic microstructures. It provides a benchmark test platform for evaluating various methods such as representative volume approaches to calculate effective properties. Furthermore, the new numerical implementation (Cheng et al., 2013) of asymptotic homogenization method of 2D and 3D materials with periodic microstructure is shown to be directly applicable to predict effective properties of periodic plates without any complicated mathematical derivation. The new numerical implementation is based on the rigorous mathematical foundation of the asymptotic homogenization method, and also simplicity similar to the representative volume method. It can be applied easily using commercial software as a black box. Different kinds of elements and modeling techniques available in commercial software can be used to discretize the unit cell. Several numerical examples are given to demonstrate the validity of the proposed methods.     

Asymptotic analysis of laminated plates and shallow shells

It was noted long ago [1] that the material strength theory develops both by improving computational methods and by widening the physical foundations. In the present paper, we develop a computational technique based on asymptotic methods, first of all, on the homogenization method [2, 3]. A modification of the homogenization method for plates periodic in the horizontal projection was proposed in [4], where the bending of a homogeneous plate with periodically repeating inhomogeneities on its surface was studied. A more detailed asymptotic analysis of elastic plates periodic in the horizontal projection can be found, e.g., in [5, 6]. In [6], three asymptotic approximations were considered, local problems on the periodicity cell were obtained for them, and the solvability of these problems was proved. In [7], it was shown that the techniques developed for plates periodic in the horizontal projection can also be used for laminated plates. In [7], this was illustrated by an example of asymptotic analysis of an isotropic plate symmetric with respect to the midplane.     

Variational asymptotic method for unit cell homogenization of periodically heterogeneous materials

A new micromechanics model, namely, the variational asymptotic method for unit cell homogenization (VAMUCH), is developed to predict the effective properties of periodically heterogeneous materials and recover the local fields. Considering the periodicity as a small parameter, we can formulate a variational statement of the unit cell through an asymptotic expansion of the energy functional. It is shown that the governing differential equations and periodic boundary conditions of mathematical homogenization theories (MHT) can be reproduced from this variational statement. In comparison to other approaches, VAMUCH does not rely on ad hoc assumptions, has the same rigor as MHT, has a straightforward numerical implementation, and can calculate the complete set of properties simultaneously without using multiple loadings. This theory is implemented using the finite element method and an engineering program, VAMUCH, is developed for micromechanical analysis of unit cells. Many examples of binary composites, fiber reinforced composites, and particle reinforced composites are used to demonstrate the application, power, and accuracy of the theory and the code of VAMUCH.     

Recent advances on the fast multipole accelerated boundary element method for 3D time-harmonic elastodynamics

This article is mainly devoted to a review on fast BEMs for elastodynamics, with particular attention on time-harmonic fast multipole methods (FMMs). It also includes original results that complete a very recent study on the FMM for elastodynamic problems in semi-infinite media. The main concepts underlying fast elastodynamic BEMs and the kernel-dependent elastodynamic FM-BEM based on the diagonal-form kernel decomposition are reviewed. An elastodynamic FM-BEM based on the half-space Green’s tensor suitable for semi-infinite media, and in particular on the fast evaluation of the corresponding governing double-layer integral operator involved in the BIE formulation of wave scattering by underground cavities, is then presented. Results on numerical tests for the multipole evaluation of the half-space traction Green’s tensor and the FMM treatment of a sample 3D problem involving wave scattering by an underground cavity demonstrate the accuracy of the proposed approach. The article concludes with a discussion of several topics open to further investigation, with relevant published work surveyed in the process.     

Homogenization of aligned “fuzzy fiber” composites

The aim of this work is to study composites in which carbon fibers coated with radially aligned carbon nanotubes are embedded in a matrix. The effective properties of these composites are identified using the asymptotic expansion homogenization method in two steps. Homogenization is performed in different coordinate systems, the cylindrical and the Cartesian, and a numerical example are presented.     

The elastodynamic Green's tensor for a half-space with an embedded anisotropic layer

The two-dimensional elastodynamic Green's tensor is found in the form of a generalized ray expansion, for an isotropic half-space in which is embedded an anisotropic layer. Particular attention is paid to the displacement of the free surface when the source transmits waves through the layer. The first motion approximation, in which individual terms are replaced by their asymptotic forms close to their arrival times, is shown to provide a fair representation without the massive computation that the full solution requires. An example which lacks symmetry shows that the layer can transmit significant SH waves from a source of P-SV type; this phenomenon is relevant to studies of the earth's upper mantle.     

Fundamental solutions in antiplane elastodynamic problem for anisotropic medium under moving oscillating source

In this paper we study the antiplane problem of concentrated point force moving with constant velocity and oscillating with constant frequency in unbounded homogeneous anisotropic elastic medium.The explicit representation of the elastodynamic Green's function is obtained by using Fourier integral transform techniques for all rates of source motion as a sum of the integrals over the finite interval. The dynamic and quasistatic components of the Green's function are extracted. The stationary phase method is applied to derive an asymptotic approximation at the far wave field. The simple formulae for Poynting energy flux vectors for moving and fixed observers are presented too.It is shown that the motion brings some differences in the far field properties, such as, for example, fast and slow waves appearance under superseismic motion and modification of the wave propagation zones and their numbers.The case of isotropic medium is considered separately. For isotropic material all main formulae are obtained in explicit forms.     

Effective Diffusion Coefficient: From Homogenization to Experiment

In this paper, an example of the application of the homogenization approach (asymptotic expansion technique) to predict the effective diffusion coefficient for an equivalent continuum, together with the experimental verification of the theoretical results is presented. The experimental setup was constructed for the measurements of diffusion in a model periodic porous medium made of Plexiglas. The computer program using the FEM was elaborated to solve the local boundary value problem for a period and to calculate the effective diffusion coefficient. The comparison between the theory and the experiment indicates good agreement between the numerical and experimental values of the effective diffusion coefficient. Interpretation of the test data from the point of view of the homogenization theory is also incorporated.     

Asymptotic method of homogenization of fissured elastic plates

This paper deals with the homogenization of thin elastic plates weakened by periodically distributed fissures. The classical Kirchhoff theory of bending plates admits five different types of unilateral fissures. To derive effective properties of fissured plates we employ the asymptotic method. As a result we obtain five effective hyperelastic plates. An illustrative example concerns the homogenization of the plate with unidirectional fissures parallel to a straight line. The constitutive equation describing the effective plate is found provided that fissures are of the flexural type resembling that observed in reinforced concrete plates in bending.     

Homogenization for wave propagation in periodic fibre-reinforced media with complex microstructure. I—Theory

The effective elastic properties of periodic fibre-reinforced media with complex microstructure are determined by the method of asymptotic homogenization via a novel solution to the cell problem. The solution scheme is ideally suited to materials with many fibres in the periodic cell. In this first part of the paper we discuss the theory for the most general situation—N arbitrarily anisotropic fibres within the periodic cell. For ease of exposition we then restrict attention to isotropic phases which results in a monoclinic composite material with 13 effective moduli and expressions for each of these are determined. In the second part of this paper we shall discuss results for a variety of specific microstructures.     

Reflection and transmission of elastic waves by the spatially periodic interface between two solids (Theory of the integral-equation method)

The linear theory of two-dimensional reflection and transmission of time-harmonic, elastic waves by the spatially periodic interface between two perfectly elastic media is developed. A given phase progression of the incident wave in the direction of periodicity induces a modal structure in the elastodynamic field and leads to the introduction of the so-called spectral orders. The main tools in the analysis are the elastodynamic Green-type integral relations. They follow from the two-dimensional form of the elastodynamic field reciprocity theorem, where in the latter a Green state adjusted to the periodicity of the structure at hand is used. One of these relations is a vectorial integral equation from which the elastodynamic field quantities can be determined.

The consequences of field reciprocity in the structure and of the conservation of energy are developed in view of their serving as a check on numercal results to be obtained from the relevant integral equations.

The formalism thus developed applies to profiles, if periodic, of arbitrary shape and size and can without too serious difficulties be implemented on a computer. The major difficulty in this respect is the relevant Green function, the series representation of it being slowly convergent. Its evaluation becomes tractable after an appropriate technique for accelerating the convergence. The only practical limitations are then put by the speed of the computer and its storage capacity.     

各向异性轴对称问题的弹性动力学解

呈现一种简便的解析方法，求解了各向异性轴对称有界结构的弹性动力学问题在实例中，计算了各向异性带孔圆板的动应力响应历程和分布规律     

Analysis of Composite Molding Using the Homogenization Method

In this paper, we present an application of the homogenization method to the analysis of Resin Transfer Molding (RTM) and Structural Reaction Injection Molding (SRIM). RTM and SRIM are relatively new molding processes for manufacturing continuous fiber reinforced polymer composites. First, the mold flow is analyzed. In the molding process, the resin experiences significant temperature changes as it fills the mold and forms a free boundary with air as it pushes out the air. In addition, the flow domain is the mold cavity packed with fiber perform, which is a porous medium. Here, the homogenization method is used to model the non-isothermal flow through porous media with free boundaries. A computer program is developed which is capable of simulating a three-dimensional mold flow using the finite element approximation. An example is provided for a three-dimensional part. Then, an analysis of the residual stress developed in the curing stage is given. The curing stage starts when the mold is completely filled and it involves chemical reaction and large temperature variation. In curing, the resin part undergoes larger volume shrinkage than the fiber part, and the residual stresses are developed due to this volume mismatch. In some cases, these stresses are large enough to cause micro-cracking and to exhaust the strength of the material. Here, a brief discussion of the application of the homogenization method to a residual stress analysis is given and one example is provided.