Equivalent inclusion method for steady state heat conduction in composites

doi:10.1016/0020-7225(86)90011-X

International Journal of Engineering Science

Volume 24, Issue 7, 1986, Pages 1159-1172

https://doi.org/10.1016/0020-7225(86)90011-X Get rights and content

Abstract

The equivalent inclusion method is proposed to solve the steady state heat conduction problems in composites. This method is analogous to Eshelby's equivalent inclusion method in elasticity. Thus, the solution procedure by this method is very simple, involving only algebraic operations. The equivalent inclusion method developed here is also applicable to the case of high concentrations of ellipsoidal inhomogeneities. As a demonstration, several examples are worked out, including the temperature field in continuous fiber composite and the thermal conductivity of two-and three-phase composites.

References (21)

T. Mori et al.
Acta Met.
(1973)
K. Tanaka et al.
Acta Met.
(1970)
J.D. Eshelby
J.D. Eshelby
Y. Takao et al.
J. Appl. Mech.
(1982)
Y. Takao et al.
J. Appl. Mech.
(1984)
T. Mura
Micromechanics of Defects in Solids
(1982)
H.S. Cawslow et al.
Conduction of Heat in Solids
(1959)
G.E. Schneider et al.
ASME 79-WA/HT-42
(1979)
L.S. Han et al.
J. Heat Transfer
(1981)

There are more references available in the full text version of this article.

Cited by (295)

Numerical modeling and experimental validation of heat and mass transfer within decomposing carbon fibers/epoxy resin composite laminates
2024, International Journal of Thermal Sciences
CFRP composite laminates decompose by pyrolysis/oxidation when exposed to high heat fluxes. While gaseous and solid residues created within the composite material are ignored in most models, a multi-species approach is proposed to describe transient decompositions from the virgin to the final charred state of the material. The model is based on a porous medium description that includes gas creation and transport within the laminate, respectively driven by Arrhenius and Darcy’s laws. Orthotropic heat transfer and decomposition gas transport within the laminate up to the surface remain often uninvestigated while it can drastically act upon the dynamics of a fire event. An experiment is carried out to analyze multi-dimensional heat and mass transfer during the charring processes of a composite laminate exposed to a high power laser beam. Transient temperature measurements on the unexposed surface of the material are confronted to 3D unsteady numerical results to validate the model formulation.
An isoparametric inclusion model for determining the thermo-elastic fields produced by varying Eigen-temperature gradients
2024, Applied Mathematical Modelling
The existence of inclusions or inhomogeneities may disrupt heat flow, resulting in increased stresses and temperature fluctuations at the material interface. Analyzing the thermo-elastic fields surrounding such microstructures is crucial for unveiling the intricate failure mechanisms within materials. Based on the method of Green's function and contour integral, this work presents a complete set of explicit analytical solutions for thermal and elastic fields of an arbitrary polygonal inclusion subjected to linearly varying eigen-temperature gradients or thermal eigenstrains. In contrast to the previous studies on potentials showing analogy with anti-plane elasticity, a mathematical analogy of Green's functions between steady-state heat conduction and thermal inclusion for plane elasticity is accordingly established. By employing the linear interpolation of the nodal eigen-temperature gradients or eigenstrains, the work further extends the traditional isoparametric element idea to the thermos-elastic inclusions, and proposes a novel isoparametric inclusion model to solve arbitrarily shaped inclusions with any distributed eigen-temperature gradients. Unlike the finite element method, which requires meshing for both the matrix and inclusion, the present model exhibits efficiency, flexibility, and versatility with the mesh discretization only conducted inside the inclusion domain. The inclusion based isoparametric method may have potential applications in mechanical engineering and material science for analyzing the thermo-structural behavior of various microstructures.
Thermomechanical modeling of functionally graded materials based on bimaterial fundamental solutions
2024, International Journal of Engineering Science
The Green’s function technique has been used to directly calculate the local fields of a functionally graded material (FGM) under thermomechanical loading, thus predicting its effective material properties. For a bi-phase FGM continuously switching the particle and matrix phases, the particle size and material gradation play a complex role in its effective material behavior. Using Eshelby’s equivalent inclusion method, particles are simulated by a source of eigen-fields in a bounded bi-layered domain, while the boundary effects are evaluated by the boundary integrals of the fundamental solutions. Using the volume integral of Green’s functions, over 10,000 particles are used to simulate an FGM under thermal and mechanical loading, respectively. The dual equivalent inclusion method is used to solve for the temperature and stress fields coupled with temperature loading. The averaged thermomechanical field distribution in the gradation direction is evaluated under different loading conditions. The effective stiffness, thermal expansion coefficient, and heat conductivity significantly change with the loading condition, particle size, and material gradation. The homogenization methods, which approximate an FGM as a continuously graded material with thermoelastic properties depending on the volume fraction only, cannot capture these micromechanical features of FGMs, while the present cross-scale approach with the inclusion-based boundary element method (iBEM) directly evaluates local fields and predicts effective material behaviors with high fidelity and efficiency.
Multiscale simulation of polymer curing of composites combined mean-field homogenisation methods at large strains
2024, International Journal of Solids and Structures
This work introduces a coupled thermo-chemo-mechanical framework for simulating the polymerisation process of composites, accounting for large strains. The framework encompasses the microscale, mesoscale and, macroscale. The microscale represents a step-growth polymerisation mechanism, which is modelled as a system of multiple growing spheres. These spheres are further decomposed into concentric layers of solid, resin, and curing agent material. At the mesoscale, the composite is divided into inclusions and the polymer matrix. The macroscale is considered as a homogeneous material. In order to describe the behaviour of the heterogeneous scales, mean-field homogenisation methods are employed. These methods have been developed before for linear elastic materials. However, considering finite deformations, a special treatment is necessary to accurately account for the effects of such deformations. Effective properties at all scales are obtained using the Mori–Tanaka, Voigt, Reuss, and Voigt–Reuss–Hill average mean-field homogenisation method. The results demonstrate how the different homogenisation methods influence the curing process and, consequently, the macroscopic response of the composite material.
A novel method for depolarization tensor and average form of an arbitrarily shaped inclusion: Extension to different physical fields and their effective transport properties of composites
2023, Composites Science and Technology
Depolarization tensor (DT), exclusively controlled by the inclusion shape, is crucial in assessing composite transport behavior. However, the solution of DT for a non-ellipsoidal inclusion exhibits nonuniformity and lacks a closed-form expression. Average DT is then introduced to ensure uniformity and seamlessly integrate the evaluation of effective transport properties of composites. In this work, we devise a novel method to obtain DT and average DT of an arbitrarily shaped inclusion by leveraging mathematical similarity from Greenʼs function to degenerate the Eshelby tensor and average Eshelby tensor in an elastic field. The method circumvents singularity and complexity in direct integration and successfully extracts all tensor components. As natural extensions, we also derive the Eshelby conduction tensor, polarizability tensor, and their averages. The accuracy of the method is validated by comparing the average polarizabilities of Platonic solids and stem cells with existing literature data. Furthermore, we incorporate average DT into a developed homogenization model to predict effective transport properties of two-phase composites comprising randomly distributed congruent irregular inclusions. Comparisons with finite element simulations demonstrate accurate predictions of the effective transport properties like effective permeability, diffusivity, conductivity, and permittivity. The proposed framework offers valuable guidance for tailoring inclusion shapes to design particulate/fibrous/porous/biological composites.
Micromechanical analysis for effective thermal and electromagnetic properties of composite materials containing anisotropical and ellipsoidal fillers oriented randomly
2023, Heliyon
In this study, micromechanical modeling will be performed for composite materials containing fillers oriented randomly in the matrix. The purpose of this study is to derive more general and explicit solutions for the effective thermal and electromagnetic properties of such composite materials without restricting the properties and shapes of the fillers. For this purpose, it is assumed that the physical properties of the filler are the same anisotropic properties as orthorhombic materials, and the shape of the fillers is ellipsoidal. This model is analyzed by micromechanics combining the Eshelby's equivalent inclusion method with the self-consistent method or the Mori-Tanaka's theory. Solutions of the effective thermal and electromagnetic properties both for composite materials containing many kinds of fillers with different shapes and physical properties and for polycrystalline materials can be also derived. Using the obtained solutions, the effect of the shape, the anisotropy, and the volume fraction of the filler on the effective thermal conductivity is examined for the carbon filler / polyethylene and the two types of quartz particles (and voids) / polyethylene. As a result, for the carbon filler / polyethylene, it is found that the effective thermal conductivity of the material when the shape of the filler is flat is about 20% higher than that when the shape of the filler is fibrous. Furthermore, when the shape of the carbon filler is flat, the result when the carbon filler is assumed to be isotropic is significantly different from that when the filler is assumed to be anisotropic. From the above, when the filler is oriented randomly in the material, it is found that simultaneously considering not only the shape of the filler but also its anisotropic properties is important to accurately evaluate the effective physical properties of the composite material. For two types of quartz particles (and voids) / polyethylene materials, the experimental result agrees better with the result of the Mori-Tanaka's theory than that of the self-consistent method, even if the volume fraction of the filler is more than 50%. From the above results, it is found that the analytical solutions of this study can generally explain the experimental results and can be applied to actual materials.

View all citing articles on Scopus

^†: On leave from Mitsubishi Electric Corporation, Polymer Engineering Department, 1-1-57 Miyashimo, Sagamihara, Kanagawa 229, Japan.

View full text

Equivalent inclusion method for steady state heat conduction in composites

Abstract

Acta Met.

Acta Met.

J. Appl. Mech.

J. Appl. Mech.

Micromechanics of Defects in Solids

Conduction of Heat in Solids

ASME 79-WA/HT-42

J. Heat Transfer