EFFECTIVE SPECIFIC HEATS OF SPHERICAL PARTICULATE COMPOSITES WITH INHOMOGENEOUS INTERPHASES

基于细观力学复合球模型研究了含非均匀界面相粒子填充复合材料的有效热弹性性质,重点讨论了界面相性质的径向分布对有效比热的影响. 首先,将非均匀界面相沿径向离散为多个同心球壳,每个球壳内的材料性质假设是均匀的. 基于上述离散模型,利用含界面相的复合球模型,推导了复合材料的有效体积模量、有效热膨胀系数及有效比热的数值求解表达式;进一步,假设界面相的性质沿径向连续变化,建立了一组微分方程,上述有效性质依赖于该微分方程组的解. 特别地,当界面相杨氏模量为幂次分布时,通过求解该微分方程组得到了有效比热等热弹性性质的解析解. 算例结果表明,应用此方法预测的有效热膨胀系数与实验结果吻合良好;界面相热膨胀系数的径向分布对有效比热和有效热膨胀系数均有显著的影响,而界面相弹性模量的径向分布对有效比热有显著的影响,对有效热膨胀系数的影响相对较小.     

含非均匀界面相粒子填充复合材料的有效比热

基于细观力学复合球模型研究了含非均匀界面相粒子填充复合材料的有效热弹性性质,重点讨论了界面相性质的径向分布对有效比热的影响. 首先,将非均匀界面相沿径向离散为多个同心球壳,每个球壳内的材料性质假设是均匀的. 基于上述离散模型,利用含界面相的复合球模型,推导了复合材料的有效体积模量、有效热膨胀系数及有效比热的数值求解表达式;进一步,假设界面相的性质沿径向连续变化,建立了一组微分方程,上述有效性质依赖于该微分方程组的解. 特别地,当界面相杨氏模量为幂次分布时,通过求解该微分方程组得到了有效比热等热弹性性质的解析解. 算例结果表明,应用此方法预测的有效热膨胀系数与实验结果吻合良好;界面相热膨胀系数的径向分布对有效比热和有效热膨胀系数均有显著的影响,而界面相弹性模量的径向分布对有效比热有显著的影响,对有效热膨胀系数的影响相对较小.      

Estimation of the effective thermoelastic moduli of fibrous nanocomposites with cylindrically anisotropic phases

Recent developments in nanotechnology make it possible to fabricate nanofibers and identify their mechanical fibers. In particular, nanofibers are used as reinforcement in composites. The present work concerns unidirectional nanofibrous composites with cylindrically anisotropic phases and aims to analytically estimate their effective thermoelastic moduli. This objective is achieved by extending the classical generalized self-consistent model to the setting of thermoelasticity, to the case of cylindrically anisotropic phases, and to the incorporation of interface stress effect. Analytical closed-form estimations are derived for all the effective thermoelastic moduli, showing that these moduli depend on the fiber cross-section size. Numerical examples are provided to illustrate this size-dependent effect.     

Effective properties of elastic periodic composite media with fibers

A series solution to obtain the effective properties of some elastic composites media having periodically located heterogeneities is described. The method uses the classical expansion along Neuman series of the solution of the periodic elasticity problem in Fourier space, based on the Green's tensor, and exact expressions of factors depending on the shape of the inclusions. Some properties of convergence of the solution are presented, more specifically concerning the elasticity tensor of the reference medium, showing that the convergence occurs even for empty fibers. The solution is extended for rigid inclusions. A comparison is made with previous exact solutions for a fiber composite made of cylindrical fibers with circular cross-sections and with previous estimates. Different examples are presented for new situations concerning the study of fiber composites: composites with elliptic cross-sections and multi-phase fibrous composites.     

Inhomogeneity of the first and second statistical moments of stresses inside the heterogeneities of random structure matrix composites

In this paper linearly thermoelastic composite media are treated, which consist of a homogeneous matrix containing a statistically homogeneous random set of heterogeneities. Effective properties (such as compliance, thermal expansion, stored energy) as well as the first statistical moments of stresses in the phases are estimated for the general case of nonhomogeneity of the thermoelastic inclusion properties. The micromechanical approach is based on the generalization of the “multiparticle effective field” method (MEFM, see for references Buryachenko, Appl. Mech. Rev. (2001), 54, 1–47), previously proposed for the estimation of stress field averages in the phases. The method exploits as a background the new general integral equation proposed by the author before and makes it possible to abandon the use of the central concept of classical micromechanics such as effective field hypothesis as well as their satellite hypothesis of “ellipsoidal symmetry”. The implicit recursion representations of the effective thermoelastic properties and stress concentration factor are expressed through some building blocks described by numerical solutions for both the one and two inclusions inside the infinite medium subjected to the inhomogeneous effective fields evaluated from subsequent self-consistent estimations. One also estimates the inhomogeneous statistical moments of local stress fields which are extremely useful for understanding the evolution of nonlinear phenomena such as plasticity, creep, and damage. Just at some additional assumptions (such as an effective field hypothesis) the involved tensors can be expressed through the Green function, Eshelby tensor and external Eshelby tensor. These estimated inhomogeneities of effective fields lead to the detection of fundamentally new effects for the local stresses inside the heterogeneities.     

Acceleration waves and ellipticity in thermoelastic micropolar media

Acceleration waves in nonlinear thermoelastic micropolar media are considered. We establish the kinematic and dynamic compatibility relations for a singular surface of order 2 in the media. An analogy to the Fresnel–Hadamard–Duhem theorem and an expression for the acoustic tensor are derived. The condition for acceleration wave’s propagation is formulated as an algebraic spectral problem. It is shown that the condition coincides with the strong ellipticity of equilibrium equations. As an example, a quadratic form for the specific free energy is considered and the solutions of the corresponding spectral problem are presented.     

Effective mechanical properties of unidirectional composites in the presence of imperfect interface

In this paper, the equivalent inclusion method is implemented to estimate the effective mechanical properties of unidirectional composites in the presence of an imperfect interface. For this purpose, a representative volume element containing three constituents, a matrix, and interface layer, and a fiber component, is considered. A periodic eigenstrain defined in terms of Fourier series is then employed to homogenize non-dilute multi-phase composites. In order to take into account the interphase imperfection effects on mechanical properties of composites, a stiffness parameter in terms of a matrix and interphase elastic modulus is introduced. Consistency conditions are also modified accordingly in such a way that only the part of the fiber lateral stiffness is to be effective in estimating the equivalent composite mechanical properties. Employing the modified consistency equations together with the energy equivalence relation leads to a set of linear equations that are consequently used to estimate the average values of eigenstrain in non-homogeneous phases. It is shown that for composites with both soft and hard reinforcements, largest stiffness parameter that indicates complete fiber–matrix interfacial debonding causes the same equivalent lateral properties.     

Micromechanical analysis of the fully coupled finite thermoelastic response of rubber-like matrix composites

A micromechanical analysis for the prediction of the coupled thermoelastic response of multiphase composites that include rubber-like phases is presented. Rubber-like solids are highly nonlinear thermoelastic materials that exhibit anomalous behavior referred to as the thermoelastic inversion effect. Results are presented which show that the derived micromechanical model is capable of predicting this effect in nylon/rubber composites subjected to appropriate thermal loadings assuming one-way coupling. For full thermomechanical coupling, the nonlinear response and induced temperatures under several types of mechanical loading are investigated.     

Derivation of energy-momentum tensors in theories of micropolar hyperbolic thermoelasticity

In the framework of the classical field theory and using the theory of action variational symmetries, we consider the construction of canonical energy-momentum tensors for a coupled micropolar thermoelastic field taking account of the nonlocality of the Lagrangian density, which is typical of continuum micromechanics. We use the algorithms of group analysis to calculate the Noether currents and the energy-momentum tensors in three cases where the Lagrangian depends on the gradients of field variables of orders not exceeding 1, 2, and 3. In each of these cases, we present explicit formulas for the components of the canonical energy-momentum tensor. We construct the energy-momentum tensor for micropolar thermoelastic bodies in which the heat conduction process is characterized by a generalized heat equation of hyperbolic analytical type. In the equations of micropolar thermoelastic field, all possible restrictions on the microrotations are taken into account.     

Determining the effective thermoelastic characteristics of regularly laminated composite in the asymmetric theory of elasticity

The asymmetric theory of elasticity is used to model a hybrid laminated composite of regular structure with all phases isotropic. The effective thermoelastic characteristics of the composite are determined. It is shown that the equations derived can be used to determine stress–strain state in all the phases of the composite using the average components of the tensors of force stresses, couple stresses, strains, and wryness in a layered material, which is of fundamental importance for the design of composites based on structural theories of failure     

Investigation on the Thermoelastic Constant of Carbon Fibre Composites

Orthotropic materials show different thermo-elastic constants depending on their fibre orientation. While most materials undergo a positive elongation with increasing temperature, carbon fibres present a heat-shrink behaviour, which in carbon fibre composites has an important consequence on thermoelastic constants. A decrease in thermoelastic constant with frequency has already been observed in glass fibre composites. Experiments made on uniaxial carbon fibre composites showed that the longitudinal thermoelastic constant increases with the frequency, while the transverse constant decreases. Furthermore, due to the opposite signs of the thermoelastic carbon fibre constant and of the surrounding resin matrix, the absolute values of the longitudinal thermoelastic constant were ten times lower in CFRPs than in GFRPs. An analytical model could successfully reproduce the frequency dependence of the longitudinal thermoelastic constant, thus helping to explain the reason for the observed behaviour. Two calibration samples were used to obtain the thermoelastic constants in the longitudinal and transverse directions. The values of the thermoelastic constants were then applied to a test sample with fibres forming an angle of 10° to the load direction. The expected theoretical results were compared with the results experimentally obtained, showing good agreement. A preliminary calibration of the longitudinal and transverse thermoelastic constants proved to be a useful approach to obtaining the correct value of the thermoelastic constant in a generic direction.     

Micromechanical analysis of the finite elastic–viscoplastic response of multiphase composites

Nonlinear thermoelastic–viscoplastic constitutive equations for large deformations with isotropic and directional hardening, are incorporated into a micromechanical finite strain analysis. As a result of this analysis, which is based on the homogenization technique for periodic microstructures, a global thermoinelastic constitutive law is established that governs the overall response of multiphase materials under finite deformations. This constitutive law is expressed in terms of the instantaneous effective mechanical and thermal stress tangent tensors together with the instantaneous global inelastic stress tensor that represents the viscoplastic effects. Results for a thermoinelastic matrix reinforced by a hyperelastic compressible material are given that illustrate the response of fibrous and particulate composites to various types of loading.     

Size-dependent effective thermoelastic properties of nanocomposites with spherically anisotropic phases

Composites made of semi-crystalline polymers and nanoparticles have a spherulitic microstructure which can be reasonably represented by a spherically anisotropic volume element. Due to the high surface-to-volume ratio of a nanoparticle, the particle-matrix interface stress, usually neglected in determining the effective elastic moduli of particle-reinforced composites, may have a non-negligible effect. To account for the latter in estimating the effective thermoelastic properties of a composite consisting of nanoparticles embedded in a semi-crystalline polymeric matrix, this work adopts a coherent interface model for the nanoparticle-matrix interface and proposes an extended version of the classical generalized-self consistent method. In particular, Eshelby's formulae widely used to calculate the elastic energy change of a homogeneous medium due to the introduction of an inhomogeneity are extended to the thermoelastic case. The nanoparticle size effect on the effective thermoelastic moduli of the composite are theoretically shown and numerically illustrated.     

Exactly solvable spherically anisotropic thermoelastic microstructures

Microstructures possessing local spherical anisotropy are considered in this paper. An example is a spherulitic polymer which can be modelled by an assemblage of spheres of all sizes in which a radial direction in every sphere is an axis of local transverse isotropy. Our purpose is to construct effectively isotropic microstructures, with spherically anisotropic and thermoelastic constituents, whose effective bulk modulus, thermal stress and specific heat can be exactly determined. The basic microstructure for which this is achieved in the present paper is the nested composite sphere assemblage of Milgrom and Shtrikman (J. Appl. Phys. 66 (1989) 3429) which was originally formulated for isotropic constituents, in the settings of conductivity and coupled fields with scalar potentials. Here, we allow the phases of this microstructure to be spherically thermoelastic with a symmetry class which can be trigonal, tetragonal, transversely isotropic, cubic or isotropic with respect to a local spherical coordinate system. A rich class of new exact results for two-phase composites and polycrystals are obtained.     

On thermoelastic dielectrics

Constitutive equations for a linear thermoelastic dielectric are derived from the energy balance equation assuming dependence of the stored energy function on the strain tensor, the polarization vector, the polarization gradient tensor and entropy. A method is indicated for constructing a hierarchy of constitutive equations for materials with arbitrary symmetry by introducing various thermodynamic potentials. Maxwell's relations are constructed for the thermodynamic potential W^L. The entropy inequality is used to obtain stability conditions for an elastic dielectric in equilibrium under prescribed boundary constraints. Frequencies are explicitly determined for a plane wave propagating along the x₁-axis in an infinite centro-symmetric isotropic thermoelastic dielectric.     

On the effective behavior of nonlinear inelastic composites: I. Incremental variational principles

A new method for determining the overall behavior of composite materials comprised of nonlinear inelastic constituents is presented. Upon use of an implicit time-discretization scheme, the evolution equations describing the constitutive behavior of the phases can be reduced to the minimization of an incremental energy function. This minimization problem is rigorously equivalent to a nonlinear thermoelastic problem with a transformation strain which is a nonuniform field (not even uniform within the phases). In this first part of the study the variational technique of Ponte Castañeda is used to approximate the nonuniform eigenstrains by piecewise uniform eigenstrains and to linearize the nonlinear thermoelastic problem. The resulting problem is amenable to simpler calculations and analytical results for appropriate microstructures can be obtained. The accuracy of the proposed scheme is assessed by comparison of the method with exact results.     

Inspection of free energy functions in gradient crystal plasticity

The dislocation density tensor computed as the cud of plastic distortion is regarded as a new constitutive variable in crystal plasticity. The dependence of the free energy function on the dislocation density tensor is explored starting from a quadratic ansatz. Rank one and logarithmic dependencies are then envisaged based on considerations from the statistical theory of dislocations. The rele- vance of the presented free energy potentials is evaluated from the corresponding analytical solutions of the periodic two-phase laminate problem under shear where one layer is a single crystal material undergoing single slip and the second one remains purely elastic.     

Axisymmetric free vibrations of infinite micropolar thermoelastic plate

The propagation of axisymmetric free vibrations in an infinite homogeneous isotropic micropolar thermoelastic plate without energy dissipation subjected to stress free and rigidly fixed boundary conditions is investigated. The secular equations for homogeneous isotropic micropolar thermoelastic plate without energy dissipation in closed form for symmetric and skew symmetric wave modes of propagation are derived. The different regions of secular equations are obtained. At short wavelength limits, the secular equations for symmetric and skew symmetric modes of wave propagation in a stress free insulated and isothermal plate reduce to Rayleigh surface wave frequency equation. The results for thermoelastic, micropolar elastic and elastic materials are obtained as particular cases from the derived secular equations. The amplitudes of displacement components, microrotation and temperature distribution are also computed during the symmetric and skew symmetric motion of the plate. The dispersion curves for symmetric and skew symmetric modes and amplitudes of displacement components, microrotation and temperature distribution in case of fundamental symmetric and skew symmetric modes are presented graphically. The analytical and numerical results are found to be in close agreement.     

Thermodynamic restrictions and wave features of a non-linear Maxwell model

A non-linear rate-type constitutive equation, established by Rajagopal, provides a generalization of the Maxwell fluid. This note embodies such a constitutive equation within the scheme of materials with internal variables thus allowing also for solids with both dissipative and thermoelastic mechanisms. The compatibility with the second law of thermodynamics, expressed by the Clausius–Duhem inequality, is examined and the restrictions on the evolution equations are determined. Next the propagation condition of discontinuity waves is derived, for shock waves and acceleration waves, by regarding the body as a definite conductor. Infinitesimal shock waves and acceleration waves show similar effects. The effective acoustic tensor proves to be the sum of a thermoelastic tensor and a tensor arising from the rate-type equation.