iMechanica - inclusionshttps://imechanica.org/taxonomy/term/351
enOn the stress singularities generated by anisotropic eigenstrains and the hydrostatic stress due to annular inhomogeneitieshttps://imechanica.org/node/17614
<div class="field field-name-taxonomy-vocabulary-6 field-type-taxonomy-term-reference field-label-hidden"><div class="field-items"><div class="field-item even"><a href="/taxonomy/term/76">research</a></div></div></div><div class="field field-name-taxonomy-vocabulary-8 field-type-taxonomy-term-reference field-label-hidden"><div class="field-items"><div class="field-item even"><a href="/taxonomy/term/4575">Geometric elasticity</a></div><div class="field-item odd"><a href="/taxonomy/term/351">inclusions</a></div><div class="field-item even"><a href="/taxonomy/term/10248">Anisotropic eigenstrain</a></div><div class="field-item odd"><a href="/taxonomy/term/2455">residual stresses</a></div><div class="field-item even"><a href="/taxonomy/term/10249">Stress singularity</a></div></div></div><div class="field field-name-body field-type-text-with-summary field-label-hidden"><div class="field-items"><div class="field-item even"><p>The problems of singularity formation and hydrostatic stress created by an inhomogeneity with eigenstrain in an incompressible isotropic hyperelastic material are considered. For both a spherical ball and a cylindrical bar with a radially-symmetric distribution of finite possibly anisotropic eigenstrains, we show that the anisotropy of these eigenstrains at the center (the center of the sphere or the axis of the cylinder) controls the stress singularity. If they are equal at the center no stress singularity develops but if they are not equal then stress always develops a logarithmic singularity. In both cases, the energy density and strains are everywhere finite. As a related problem, we consider annular inclusions for which the eigenstrains vanish in a core around the center. We show that even for an anisotropic distribution of eigenstrains, the stress inside the core is always hydrostatic. We show how these general results are connected to recent claims on similar problems in the limit of small eigenstrains.</p>
</div></div></div><div class="field field-name-upload field-type-file field-label-hidden"><div class="field-items"><div class="field-item even"><table class="sticky-enabled">
<thead><tr><th>Attachment</th><th>Size</th> </tr></thead>
<tbody>
<tr class="odd"><td><span class="file"><img class="file-icon" alt="PDF icon" title="application/pdf" src="/modules/file/icons/application-pdf.png" /> <a href="https://imechanica.org/files/AnnularInclusion_YaGo14.pdf" type="application/pdf; length=383316">AnnularInclusion_YaGo14.pdf</a></span></td><td>374.33 KB</td> </tr>
</tbody>
</table>
</div></div></div>Sun, 07 Dec 2014 22:14:59 +0000Arash_Yavari17614 at https://imechanica.orghttps://imechanica.org/node/17614#commentshttps://imechanica.org/crss/node/17614Nonlinear elastic inclusions in isotropic solidshttps://imechanica.org/node/15307
<div class="field field-name-taxonomy-vocabulary-6 field-type-taxonomy-term-reference field-label-hidden"><div class="field-items"><div class="field-item even"><a href="/taxonomy/term/76">research</a></div></div></div><div class="field field-name-taxonomy-vocabulary-8 field-type-taxonomy-term-reference field-label-hidden"><div class="field-items"><div class="field-item even"><a href="/taxonomy/term/351">inclusions</a></div><div class="field-item odd"><a href="/taxonomy/term/2455">residual stresses</a></div><div class="field-item even"><a href="/taxonomy/term/4575">Geometric elasticity</a></div></div></div><div class="field field-name-body field-type-text-with-summary field-label-hidden"><div class="field-items"><div class="field-item even"><p>We introduce a geometric framework to calculate the residual stress fields and deformations of nonlinear solids with inclusions and eigenstrains. Inclusions are regions in a body with different reference configurations from the body itself and can be described by distributed eigenstrains. Geometrically, the eigenstrains define a Riemannian 3-manifold in which the body is stress-free by construction. The problem of residual stress calculation is then reduced to finding a mapping from the Riemannian material manifold to the ambient Euclidean space. Using this construction, we find the residual stress fields of three model systems with spherical and cylindrical symmetries in both incompressible and compressible isotropic elastic solids. In particular, we consider a finite spherical ball with a spherical inclusion with uniform dilatational eigenstrain and we show that the stress in the inclusion is uniform and hydrostatic. We also show how singularities in the stress distribution emerge as a consequence of a mismatch between radial and circumferential eigenstrains at the centre of a sphere or the axis of a cylinder.</p>
</div></div></div><div class="field field-name-upload field-type-file field-label-hidden"><div class="field-items"><div class="field-item even"><table class="sticky-enabled">
<thead><tr><th>Attachment</th><th>Size</th> </tr></thead>
<tbody>
<tr class="odd"><td><span class="file"><img class="file-icon" alt="PDF icon" title="application/pdf" src="/modules/file/icons/application-pdf.png" /> <a href="https://imechanica.org/files/Inclusions_YavGor2013_0.pdf" type="application/pdf; length=489942" title="Inclusions_YavGor2013.pdf">Inclusions_YavGor2013.pdf</a></span></td><td>478.46 KB</td> </tr>
</tbody>
</table>
</div></div></div>Fri, 13 Sep 2013 15:07:24 +0000Arash_Yavari15307 at https://imechanica.orghttps://imechanica.org/node/15307#commentshttps://imechanica.org/crss/node/15307ISDMM11 - Upcoming Abstract Deadline January 30th 2011 https://imechanica.org/node/9686
<div class="field field-name-taxonomy-vocabulary-6 field-type-taxonomy-term-reference field-label-hidden"><div class="field-items"><div class="field-item even"><a href="/taxonomy/term/74">conference</a></div></div></div><div class="field field-name-taxonomy-vocabulary-8 field-type-taxonomy-term-reference field-label-hidden"><div class="field-items"><div class="field-item even"><a href="/taxonomy/term/351">inclusions</a></div><div class="field-item odd"><a href="/taxonomy/term/420">cracks</a></div><div class="field-item even"><a href="/taxonomy/term/499">dislocations</a></div><div class="field-item odd"><a href="/taxonomy/term/3087">phase boundaries</a></div><div class="field-item even"><a href="/taxonomy/term/3088">shape optimization</a></div><div class="field-item odd"><a href="/taxonomy/term/5791">precipitates</a></div></div></div><div class="field field-name-body field-type-text-with-summary field-label-hidden"><div class="field-items"><div class="field-item even"><p class="title">
5th International Symposium on Defect and Material Mechanics, Seville (Spain), <span>from June 27 to July 1, 2011. </span>
</p>
<p>
<span>Scientific Committee:</span><span><span> </span></span>
</p>
<p>
<span>P. Ariza (University of Sevilla, Spain) <strong>Chair</strong><br />
D. Bigoni (University of Trento, Italy) <br />
R. Kienzler (University of Bremen, Germany)<br />
X. Markenscoff (University of California-San Diego, USA)<br />
A. Needleman (University of North Texas, USA)<br />
M. Ortiz (California Institute of Technology, USA)<br />
P. Steinmann (University of Erlangen-Nuremberg, Germany)<br />
C. Stolz (Ecole Polytechnique, France)<br />
V. Tvergaard (Technical University of Denmark, Denmark)</span><span> </span>
</p>
<p>
<span>Plenary Speakers (confirmed plenay speakers to date)</span>
</p>
<p>
<span>Gilles Francfort</span> (Universite Paris Nord)<br /><span>Ben Freund</span> (Brown University)<br /><span>Richard James</span> (University of Minnesota)<br /><span>Rodolfo Miranda</span> (LASUAM Madrid)<br /><span>Nick Schryvers</span> (University of Antwerp)
</p>
<p>
</p>
<p>
<span>For more information visit </span><a href="http://congreso.us.es/isdmm2011/"><span>http://congreso.us.es/isdmm2011/</span></a>
</p>
<p>
</p>
<p>
</p>
</div></div></div>Mon, 24 Jan 2011 17:21:53 +0000jsanz9686 at https://imechanica.orghttps://imechanica.org/node/9686#commentshttps://imechanica.org/crss/node/96865th International Symposium on Defect and Material Mechanics - Abstract Deadline January 30th 2011 https://imechanica.org/node/9344
<div class="field field-name-taxonomy-vocabulary-6 field-type-taxonomy-term-reference field-label-hidden"><div class="field-items"><div class="field-item even"><a href="/taxonomy/term/74">conference</a></div></div></div><div class="field field-name-taxonomy-vocabulary-8 field-type-taxonomy-term-reference field-label-hidden"><div class="field-items"><div class="field-item even"><a href="/taxonomy/term/351">inclusions</a></div><div class="field-item odd"><a href="/taxonomy/term/420">cracks</a></div><div class="field-item even"><a href="/taxonomy/term/499">dislocations</a></div><div class="field-item odd"><a href="/taxonomy/term/3087">phase boundaries</a></div><div class="field-item even"><a href="/taxonomy/term/3088">shape optimization</a></div><div class="field-item odd"><a href="/taxonomy/term/5791">precipitates</a></div></div></div><div class="field field-name-body field-type-text-with-summary field-label-hidden"><div class="field-items"><div class="field-item even"><p>
<span>We are delighted to invite you to submit your contribution to ISDMM11. The conference will take place in Seville, Spain, June 27 – July 1st 2011. </span>
</p>
<p>
<span>ISDMM11 is the fifth international meeting devoted to Mechanics of Material Forces, following the workshops held at Kaiserslautern (2003), Symi (2005), Aussois (2007) and Trento (2009), and will be held in the city of Seville, Spain, from June 27 to July 1, 2011. The Symposium is intended to bring together researchers in the areas of the mechanics of defects, in a broad sense: cracks, dislocations, inclusions, precipitates, phase boundaries, and shape optimization.</span>
</p>
<p>
<span>Scientific Committee:</span><span><span> </span></span>
</p>
<p>
<span></span><span>P. Ariza (University of Sevilla, Spain) <strong>Chair</strong><br />
D. Bigoni (University of Trento, Italy) <br />
R. Kienzler (University of Bremen, Germany)<br />
X. Markenscoff (University of California-San Diego, USA)<br />
A. Needleman (University of North Texas, USA)<br />
M. Ortiz (California Institute of Technology, USA)<br />
P. Steinmann (University of Erlangen-Nuremberg, Germany)<br />
C. Stolz (Ecole Polytechnique, France)<br />
V. Tvergaard (Technical University of Denmark, Denmark)</span><span> </span>
</p>
<p>
<span>Plenary Speakers (confirmed plenay speakers to date)</span></p>
<p><span>Gilles Francfort</span> (Universite Paris Nord)<br /><span>Ben Freund</span> (Brown University)<br /><span>Richard James</span> (University of Minnesota)<br /><span>Rodolfo Miranda</span> (LASUAM Madrid)<br /><span>Nick Schryvers</span> (University of Antwerp)
</p>
<p>
<span>For more information visit </span><a href="http://congreso.us.es/isdmm2011/"><span>http://congreso.us.es/isdmm2011/</span></a>
</p>
</div></div></div>Fri, 19 Nov 2010 18:58:09 +0000jsanz9344 at https://imechanica.orghttps://imechanica.org/node/9344#commentshttps://imechanica.org/crss/node/9344Modelling inclusion (second phase particles) by Ansyshttps://imechanica.org/node/7182
<div class="field field-name-taxonomy-vocabulary-6 field-type-taxonomy-term-reference field-label-hidden"><div class="field-items"><div class="field-item even"><a href="/taxonomy/term/76">research</a></div></div></div><div class="field field-name-taxonomy-vocabulary-8 field-type-taxonomy-term-reference field-label-hidden"><div class="field-items"><div class="field-item even"><a href="/taxonomy/term/351">inclusions</a></div></div></div><div class="field field-name-body field-type-text-with-summary field-label-hidden"><div class="field-items"><div class="field-item even"><p>
Hi all
</p>
<p>
</p>
<p>
I am phd student and working on aluminium alloy 2214. i wanted to know if anyone here has experience of FE modelling of inclusions present in alloys by Ansys? i have characterized my material n i know inclusion distribution n size. can anyone help me in this regard?
</p>
<p>
</p>
<p>
</p>
<p>
Merci en avance
</p>
</div></div></div><div class="field field-name-upload field-type-file field-label-hidden"><div class="field-items"><div class="field-item even"><table class="sticky-enabled">
<thead><tr><th>Attachment</th><th>Size</th> </tr></thead>
<tbody>
<tr class="odd"><td><span class="file"><img class="file-icon" alt="Image icon" title="image/jpeg" src="/modules/file/icons/image-x-generic.png" /> <a href="https://imechanica.org/files/1.JPG" type="image/jpeg; length=13210" title="1.JPG">1.JPG</a></span></td><td>12.9 KB</td> </tr>
</tbody>
</table>
</div></div></div>Tue, 01 Dec 2009 17:22:09 +0000msa4u7182 at https://imechanica.orghttps://imechanica.org/node/7182#commentshttps://imechanica.org/crss/node/7182On Eshelby's two classicshttps://imechanica.org/node/3529
<div class="field field-name-taxonomy-vocabulary-8 field-type-taxonomy-term-reference field-label-hidden"><div class="field-items"><div class="field-item even"><a href="/taxonomy/term/173">Eshelby</a></div><div class="field-item odd"><a href="/taxonomy/term/351">inclusions</a></div><div class="field-item even"><a href="/taxonomy/term/352">inhomogeneities</a></div><div class="field-item odd"><a href="/taxonomy/term/2638">elastic stress</a></div><div class="field-item even"><a href="/taxonomy/term/2639">misfit stresses</a></div><div class="field-item odd"><a href="/taxonomy/term/2640">eigenstrain</a></div></div></div><div class="field field-name-body field-type-text-with-summary field-label-hidden"><div class="field-items"><div class="field-item even"><p>
Recently, <a href="http://ontheshouldersofgiants.wordpress.com/about/">a new carnival called The Giant's Shoulders</a> has been started and <a href="http://scienceblogs.com/clock/2008/07/the_giants_shoulders_1.php">the first edition of the same is out at A blog around the clock</a>. A post of mine on the <a href="http://mogadalai.wordpress.com/2008/07/15/elastic-stresses-due-to-inclusions-and-inhomogeneities/">elastic stresses due to inclusions and inhomogeneities</a> made it to the carnival. I am cross-posting the piece here since it might also be of interest to the readers of iMechanica (though I did post a short note earlier <a href="node/456">here</a> which forms the core of this long post too).
</p>
<p><strong>Crystallanity</strong></p>
<p>
Almost all the metals and alloys that are used in practical applications are crystalline -- that is, the atoms or molecules that make up the metal or alloy are arranged periodically in space. For the sake of simplicity (and, without loss of generality), in this post, I am assuming that this periodic arrangement can be built up of cubes -- the corners and face centers of which are populated by atoms. This specific crystal structure is known as <a href="http://en.wikipedia.org/wiki/Face-centred_cubic">face centered cubic</a> (fcc).
</p>
<p>
The periodic arrangement of atoms/molecules in a crystal results in many important and interesting properties. One of them is the lattice parameter, which, in our case, is the size of the cube, or the distance between the centers of two atoms that are occupying the cube corners.
</p>
<p>
The elasticity also follows from the crystalline structure rather naturally. In a crystalline solid, the atoms act as if they are connected by springs, and the crystalline structure that a particular metal or alloy chooses is dictated, at some level, also by these “springs” and their strength. Hence, in a crystal, if you try to move any atom from its equilibrium position (determined by the lattice parameter), it tries to go back; if it can't, the springs that attach it to the other atoms are either stretched or compressed; these stretchings and compressions are what make the crystal elastic. As soon as the forces on these atoms (which moved them away from their equilibrium positions) are removed, the atoms go back to their original positions. The elastic constants of a material tells us, for a given force, how much these atoms can be strained.
</p>
<p>
The crystal structure, by definition, makes the solid anisotropic (that is, if you sit on an atom and look at different directions, its properties are different); hence, the elastic constants naturally inherit the anisotropy of the underlying crystalline structure. In the case of cubic crystals, the elastic constants are obviously cubic anisotropic. What this means in practical terms is that, if you look in the directions of cube edges and the cube diagonals, the elastic properties are different; more specifically, either the cube edge direction or the diagonal direction is elastically softer as compared to the other; that is, for a given force, the atoms in the softer direction are relatively more pliable.
</p>
<p><strong>Dual phases</strong></p>
<p>
Many of the metallic materials used in practical applications are not only alloys (that is, they consist of more than one type of elements) and polycrystalline (that is, each material consists of several crystallites), but also consist of more than one phase (that is, consists of solid material that has different physical properties). The different phases and the different combination of crystallites give rise to a wide variety of interesting microstructural features (features initially noticed at the micrometre scale, and hence the name), which, in turn, give rise to several interesting (and some times important) properties to the material. Thus, it is no wonder that a large fraction of materials scientists and engineers are interested in studying the microstructural features, their effects on properties and ways of tuning both.
</p>
<p>
For this particular post, I am going to consider a specific model alloy which consists of two elements -- nickel and aluminium (in practice, several other elements too -- but for our purposes, it is sufficient to deal with the alloy as if it consists of only these two elements). It also consists of two phases; one of them, the nickel rich phase, has fcc crystal structure; the other, which consists of the specific fraction of three nickel atoms for each aluminium atom, crystallises in a structure that is very close to fcc called $latex L1_2$; in fact, it is the fcc structure, except that the aluminium atoms prefer to occupy the cube corners while nickel atoms occupy the face centers. Naturally, these two phases have different lattice parameters and elastic constants; however, since both these phases are cubic crystalline structure based, the elastic anisotropy is the same for both phases. The typical microstructure in this material consists of cuboids of the $latex L1_2$ phase (precipitates) distributed in the fcc material (matrix), making it look like a miniature version of bricks ($latex L1_2$) and mortar (fcc) in masonry.
</p>
<p><strong>Misfit and elastic inhomogeneity</strong></p>
<p>
In the case of dual phase materials like the one described above, the complexity of the microstructure of the materials also leads to several other important materials properties and parameters, of which, two are of specific interest to us. One is known as the misfit, which gives the difference in lattice parameter between the matrix and precipitate phases (normalised by the matrix phase). The second property is called inhomogeneity -- that is, at different parts of the material, the properties (specifically, the elastic constants) are different (since the phases are different).
</p>
<p>
Both the <a href="http://en.wikipedia.org/wiki/Precipitation_strengthening">misfit and elastic inhomogeneity play a key role in strengthening of alloys</a>; in fact, this is the key process that leads to the superior mechanical properties of <a href="http://en.wikipedia.org/wiki/Superalloy">superalloys -- alloys which are used in aerospace industry and in the making of gas turbines</a>; and, the $latex L1_2 Ni_3Al$ precipitates in nickel rich fcc is in fact <em>the</em> most important ingredients of nickel-base superalloys.
</p>
<p><strong>The Problem</strong></p>
<p>
There are several interesting questions that one can ask about the microstructure and its evolution in dual phase (for simplicity, single crystalline) alloy materials of the type described above. In this post, we will ask one such question, namely, what are the elastic stress and strain fields associated with the microstructure in such materials? This question of the elastic stresses and strains is of interest both from (a) the point of view of understanding these materials and their properties, and, (b) from the point of view of using these materials in practical applications; thus, it is no wonder that the two papers that J D Eshelby wrote outlining a process for obtaining these fields (albeit for the case of some special geometries) have become classics in the field.
</p>
<p>
Both these were published in the <em>Proceedings of the Royal Society of London: Series A. Mathematical and Physical Sciences</em>. The first, published in 1957, is titled <strong><em>The determination of the elastic field of an ellipsoidal inclusion, and related problems</em></strong> [1]. The second, published two years after the first in 1959, is titled <strong><em>The elastic field outside an ellipsoidal inclusion</em></strong> [2]. These two papers are very accessible and are a pleasure to read - and are a must read for anybody who is interested in theoretical materials science.
</p>
<p>
Of course, Eshelby is not the first scientist to look at the problem; in his first paper, he does list several authors who have discussed many special cases of the problem; however, as Eshelby himself notes, not only is his method rather simple and straight-forward, it is also the most general; this generality allows for easier generalization in applying Eshelby's methodology for shapes that are not simple and obtaining elastic solutions numerically in such cases. Further, the simplicity of the methodology, and the exact solutions given by him for the special cases also serve another important purpose, namely, that of validating codes written for evaluating stress and strain fields associated with more complex shapes. In fact, Eshelby's paper also notes several mistakes in the usual expressions for these fields given for certain special cases prior to this work -- some by Eshelby himself, and some by masters in the field such as Landau and Lifshitz in their classic textbook. Thus, these two classics not only clarified several important issues, but, in some sense, actually gave the road map for future work in the area of micromechanics of defects in solids for the next several decades; even today, not only the attempts to push Eshelby's methodology to more and more general (and, more and more complex) situations, but also the efforts to obtain similar, exact analytical expressions for shapes other than those discussed by Eshelby, namely, ellipsoids of revolution, are being pursued.
</p>
<p>
As with all classics, the appreciation of the work, in the minds of its readers, is always enhanced with a little bit of knowledge about the authors. So, before I proceed with the post, I would like to draw the attention of the readers to the memoir on Eshelby, published by B A Bilby in the <em>Biographical memoirs of Fellows of Royal Society</em> titled <strong><em>John Douglas Eshelby. 21 December 1916 – 10 December 1981</em></strong> [3]. There is also a short, one page appreciation (more about the man than about his work) of Eshelby, published by Alex D King in the <strong><em>Posterminaries</em></strong> column of the MRS Bulletin in July 1999 [4]; this piece, among other things, tells as to why Eshelby chose to work on theoretical problems instead of doing experiments, and <a href="http://gururajan.mp.googlepages.com/eshelby">how he came to be awarded FRS</a>; both of these stories are extremely interesting and funny (if a bit apocryphal). Finally, I should also mention the book <strong><em>The coming of materials science</em></strong> by Robert W Cahn [5], which puts this and other classics that Eshelby wrote (Yes; though he wrote only 56 papers or so in his entire career, several of them are considered classics and continue to be studied with great interest – and, at times, with awe and reverence) in the general perspective of growth of materials science as a discipline.
</p>
<p><strong>Elastic stresses during phase transformations</strong></p>
<p>
The easiest way to imagine as to why there might be elastic stresses and strains in such dual phase materials, and what these stresses or strains might look like, consider a material that is initially in the fcc crystal structure; now, let us assume that one small part of this material transforms into $latex L1_2$; since this part now has a different lattice parameter and hence a misfit as compared to the original fcc phase, the transformed region either wants to shrink or expand (depending on the sign of the misfit).
</p>
<p>
If suppose the fcc phase is infinitely rigid, all the strain associated with the expansion or shrinkage will be accommodated by the transformed region; and hence, it will be in a distorted geometry as compared to what it would in free state —that is, if it was not surrounded by this rigid fcc phase.
</p>
<p>
On the other hand, it is also possible to imagine that the fcc phase is infinitely compliant as compared to the transformed region; in this case, the transformed region would look exactly like what it would look without the surrounding material; however, all the strain associated with the transformation is now accommodated by the straining of the surrounding fcc phase.
</p>
<p>
In reality, neither of the phases are infinitely rigid; in fact, a parameter called inhomogeneity ratio $latex \delta$, can be defined which tells how relatively rigid the transformed region is as compared to the original phase; and, corresponding to $latex \delta$, the strain due to misfit is distributed in both the phases; thus, the problem is to find the exact value of the stresses and strains at various points given the misfit and the shape of the transforming region.
</p>
<p>
In fact, even in the case where we assume that the different phases have the same elastic constants and differ only in the lattice parameter, the strain fields are not easily evaluated; and, depending on the geometry of the transformed region, the resulting strain fields might have many subtle and aesthetically pleasing properties (as was pointed out by Eshelby for the first time): for example, for isotropic, circular inclusions, if the eigenstrain is dilatational (no shear components), the shear stresses outside the inclusions are zero; the principal stresses are equal and opposite in sign (and, one of them is discontinuous while the other continuous at the inclusion-matrix boundary).
</p>
<p>
As the titles of the papers indicate, Eshelby only dealt with a particular class of shapes, namely, ellipsoids of revolution. However, his methodology forms the basis for a numerical evaluation of the strain fields for arbitrary shapes (which is one the reasons why the paper is so essential and influential). Having said that, note that ellipsoids of revolution, contain within them, as special cases, several important shapes that one is interested in general for evaluating many material properties – like needles, spheres, plates, etc. So, the paper and some of the results presented in them are of great use by themselves.
</p>
<p><strong>Elastic solutions (using the elegant Eshelbian cuts, strains, and weldings)</strong></p>
<p>
The exact problem that Eshelby solved ("with the help of a simple set of imaginary cutting, straining and welding operations"), in his own words, is the following:
</p>
<ul><li>The transformation problem</li>
</ul><p>
A region (the "inclusion") in an infinite homogeneous isotropic elastic medium undergoes a change of shape and size which, but for the constraint imposed by its surroundings (the "matrix"), would be an arbitrary homogeneous strain. What is the elastic state of inclusion and matrix?
</p>
<p>
The homogeneous strain is known as "eigenstrain" or "transformational strain". In the same paper, Eshelby also introduced the concept of "equivalent inclusion" for solving the transformation problem when the matrix and the region of eigenstrain (the "inhomogeneity") have different elastic constants.
</p>
<p>
The operations that Eshelby used to solve the transformation problem are the following:
</p>
<ul><li>Remove the region of interest from the matrix.</li>
<li>Allow it to take the eigenstrain.</li>
<li>Restore the region to its original shape and size by applying suitable surface tractions and put it back into the matrix and rejoin.</li>
<li>Remove the body force on the interface between the inclusion and matrix by applying an equal and opposite layer of body force.</li>
</ul><p>
In step (3), the stress is zero in the matrix and is a known constant in the inclusion. The additional stress introduced in step (4) is found by the integration from the expression for the elastic field of a point force.
</p>
<p>
Eshelby (with his cuts, strains and weldings) has shown that the transformation problem is equivalent to solving for the equations of elastic equilibrium of a homogeneous body with a known body force distribution; he also gave exact expressions for the stress and strain fields provided the region of interest is an ellipsoid of revolution; his methodology has been generalised to arbitrary geometries and multiple inclusions/inhomogeneities (with the resultant equations being solved numerically, of course) by Khachaturyan and his co-workers [6].
</p>
<p>
For homogeneous bodies with known body force distribution, the equations of elastic equilibrium are solved using the elastic Green function. Rob Philips and Mura describe the Green function approach in great detail [7,8].
</p>
<p>
In case you are interested in looking at some solutions of inclusion problems, Rob Philips describes the radial displacements associated with a spherical inclusion of radius "a" with dilatational eigenstrain obtained using Green functions (See the figure 10.14 on page 524) and indicates that the elastic energy of a spherical inclusion with dilatational misfit scales as the volume of the inclusion. Solutions for elliptic <a href="http://gururajan.mp.googlepages.com/inclusions">inclusions</a> and <a href="http://gururajan.mp.googlepages.com/inhomogeneities">inhomogeneinities</a> are also available at <a href="http://gururajan.mp.googlepages.com/inhomclusions">my googlepage</a>; and, sometime soon, I will upload the code used for obtaining these fields, and give a link to the page here.
</p>
<p>
Mura's classic is a singular testimony to the power of a combination of Green function and the eigenstrain approach of Eshelby. The Green function approach finally results in the evaluation of elliptic integrals for obtaining the displacements. This is not surprising since we are integrating the body forces over an elliptic geometry (Remember, the inclusions were ellipsoids/ellipses). The gradients of displacement give us the strain - That means we need to differentiate the Green functions. Thus, it is rather cumbersome to calculate the elastic stress, strain, or displacement fields using the Eshelby-Green approach.
</p>
<p><strong>A digression</strong></p>
<p>
By the way, Green is another fascinating figure in the annals of science/mathematics, and I understand that <a href="http://gururajan.mp.googlepages.com/green">the paper in which Green introduced the idea of the functions that bear his name is a classic by itself</a>.
</p>
<p><strong>'Equivalent inclusion' for inhomogeneities</strong></p>
<p>
The idea of "equivalent inclusion" is simple - It is one of those nice mathematical tricks where we solve a problem by reducing it to another which has already been solved.
</p>
<p>
Let the inhomogeneity have an elastic constant that is different from that of the matrix. The idea is to replace the inhomogeneity with an inclusion - The eigenstrain in the hypothetical inclusion is such that it exerts the same actions on the matrix as the original inhomogeneity. Mathematically, finding out the eigenstrain in the equivalent inclusion amounts to solving for a set of three equations in three unknowns (in 2D) or six equations in six unknowns (in 3D).
</p>
<p><strong>Going beyond Eshelby (at least in 2D): the complex variable formalism</strong></p>
<p>
If we are interested in ellipses and not ellipsoids (that is, 2D problems), it is possible to avoid the cumbersome integrals mentioned earlier. In 1960, in a paper published in the <em>Proceedings of Cambridge Philosophical Society</em>, Jaswon and Bhargava showed how to avoid the elliptic integrals using a complex variable formulation [9].
</p>
<p>
Jaswon and Bhargava motivated their complex variable formulation by making the following observation:
</p>
<ul><li>Although Eshelby has proved some general theorems of great interest, using elegant methods, his solutions involve analytically intractable integrals of a formidable nature.</li>
</ul><p>
Eshelby himself felt the same way, since in his 1959 paper [2] he says:
</p>
<ul><li>It has to be admitted that, except in the simplest cases, a calculation of the external field is laborious.</li>
</ul><p>
Jaswon and Bhargava built their complex variable formalism on the "ingeneous attack" on the transformation problem by Eshelby "utilizing the point-force concept". In view of the "novelty and importance" of the approach, they also give a brief description of Eshelby's arguements; and, their description is by far one of the best that I have seen in the literature.
</p>
<p>
The solution of Jaswon and Bhargava is based on the following ideas/results:
</p>
<ul><li>Eshelby's method involves integrals of point forces on the matrix-inclusion boundary.</li>
<li>The expression for the displacement at any point $latex x$<em> </em>due to a point force $latex F$<em> </em>acting at the point $latex y$ being known, the Eshelby problem now reduces to an integration of a continuous distribution of the forces over the inclusion surface.</li>
<li>Green and Zerna [10] and Mushkelishvili [11] give the expressions for writing down the contour integrals!</li>
</ul><p><strong>A classic, for all times!</strong></p>
<p>
Bilby, writing in 1990 [3], notes that the first of Eshelby's paper discussed here was in the second group of 100 most highly cited papers in all fields of science covered by the Science Citation Index, 1955--1986. However, citations tell only one part of the story. The ideas and results due to Eshelby have become text book material and are being used continually; thus, sometimes references are not made to his paper but to some text book that describes the same methodology.
</p>
<p>
Bilby, in the same memoir, also notes that
</p>
<p>
This work on inclusions and inhomogneities has been applied by others, not only to calculate the stress fields and interactions of inclusions, inhomogeneities, precipitates, twins, martensite plates, cavities and cracks, but also to find the bulk elastic properties of bodies containing distributions of inhomogeneities and cavities and to discuss the properties of polycrystals and composite materials. Eshelby sketched many of these applications and noted also that the method could be applied to find the perturbation caused in a slow viscous flow by the presence of a rigid or deforming ellipsoid. (...) He was particularly interested in the viscous problem because the fact that an ellipsoid remains an ellipsoid under homogeneous deformationmeans that its finite change of shape can be studied without the need to find the details of the complicated flow outside it. This application, which has seen considerable further development, is relevant to the theory of the homogenization of glass, the determination of strain in rocks, the deformation of voids and the flow of suspensions containing rigid or defromable particles.
</p>
<p>
After nearly 18 years, today, if, had Bilby been writing an appreciation of these papers, he would have included several more fields where his results cotninue to play a crucial role -- one of them being the study of microstructural evolution in elastically inhomogeneous solids with defects (to which, Khachaturyan and his co-workers have contributed immensely [6]).
</p>
<p>
However, as Eshelby himself seems to have noted, the most important reason why these two papers are classics are not for their results, but for the methodology that was developed in them, which continue to be of use fifty years after he published these papers [3]:
</p>
<p>
However, he liked to regard himself as a humble 'supplier of tools for the trade' and often left their detailed use to others.
</p>
<p>
In the process, <em>this supplier of tools</em>, also made the tools themselves more respectable (as is noted in the context of a prior paper of Eshelby which also uses the imaginary cutting, straining and welding operations) [3]:
</p>
<p>
Notable in this paper of 1951 is his derivation of results by the use of imaginary cutting, straining and welding operations, a technique that he used frequently with great effect. The method was not regarded as quite respectable by some with a more formal mathematical training.
</p>
<p>
Finally, as with all classics, these two papers of Eshelby should be read not only because of their relevance and use, but also for their beauty and elegance, which gives the reader so much of pleasure!
</p>
<p><strong>References</strong></p>
<p>
[1] J D Eshelby, <em>The determination of the elastic field of an ellipsoidal inclusion, and related problems</em>, Proceedings of the Royal Society of London: Series A. Mathematical and Physical Sciences, <strong>241</strong>, p. 376, 1957.
</p>
<p>
[2] J D Eshelby, <em>The elastic field outside an ellipsoidal inclusion</em>, Proceedings of the Royal Society of London: Series A. Mathematical and Physical Sciences, <strong>252</strong>, p. 561, 1959.
</p>
<p>
[3] B A Bilby, <em>John Douglas Eshelby. 21 December 1916 – 10 December 1981</em>, Biographical memoirs of Fellows of Royal Society, <strong>36</strong>, p. 126, 1990.
</p>
<p>
[4] A H King, <em>Posternimnaries: Lessons from J D Eshelby</em>, M R S Bulletin, p. 80, July 1999.
</p>
<p>
[5] R W Cahn, <em>The coming of materials science</em>, Pergamon materials series, Elsevier Science Publishers, 2003.
</p>
<p>
[6] A G Khachaturyan, <em>Theory of structural transformations in solids</em>, John Wiley &amp; Sons (p. 198), 1983; A G Khachaturyan, S Semenovskaya and T Tsakalakos, <em>Elastic strain energy of inhomogeneous solids</em>, Physical Review B, <strong>52</strong>, p. 15909, 1992.
</p>
<p>
[7] R Philips, <em>Crystals, defects and microstructures: Modeling across scales </em>, Cambridge University Press (p. 520), 2001.
</p>
<p>
[8] T Mura, <em>Micromechanics of defects in solids</em>, Kluwer academic publishers (p. 74), 1987.
</p>
<p>
[9] M A Jaswon and R D Bhargava, <em>Two-dimensional elastic inclusion problems</em>, Proceedings of Cambridge Philosophical Society, <strong>57</strong>, p. 669, 1960.
</p>
<p>
[10] A E Green and W Zerna, <a href="http://books.google.com/books?id=7q-Bi1A1fXMC&amp;pg=PT1&amp;dq=Green+and+Zerna+Theory+of+elasticity&amp;ei=xyl7SIS3H4eusgOqsfy5Aw&amp;sig=ACfU3U0V5u0n5HcoHfRhAT81_DYa-e6JWg"><em>Theoretical elasticity</em></a>, Oxford University Press, 1968.
</p>
<p>
[11] N I Muskhelishvili, <a href="http://books.google.com/books?hl=en&amp;id=WsOLds7xSnwC&amp;dq=Muskhelishvili+some+basic+problems&amp;printsec=frontcover&amp;source=web&amp;ots=JY-PtzSbm2&amp;sig=cIudFlQAThtBsK4ZAcCDLc_KyEU&amp;sa=X&amp;oi=book_result&amp;resnum=1&amp;ct=result#PPP1,M1"><em>Some basic problems of the mathematical theory of elasticity</em></a>, Springer, 1975.
</p>
<p>
</p>
</div></div></div>Fri, 18 Jul 2008 02:31:17 +0000Mogadalai Gururajan3529 at https://imechanica.orghttps://imechanica.org/node/3529#commentshttps://imechanica.org/crss/node/3529Eshelby and his two classics (and some more on the side)https://imechanica.org/node/456
<div class="field field-name-taxonomy-vocabulary-8 field-type-taxonomy-term-reference field-label-hidden"><div class="field-items"><div class="field-item even"><a href="/taxonomy/term/347">elasticity</a></div><div class="field-item odd"><a href="/taxonomy/term/351">inclusions</a></div><div class="field-item even"><a href="/taxonomy/term/352">inhomogeneities</a></div><div class="field-item odd"><a href="/taxonomy/term/353">Green function</a></div></div></div><div class="field field-name-body field-type-text-with-summary field-label-hidden"><div class="field-items"><div class="field-item even"><p>Eshelby and the inclusion/inhomogeneity problems
</p><p>Any materials scientist interested in mechanical behaviour would be aware of the contributions of J.D. Eshelby. With 56 papers, Eshelby revolutionised our understanding of the theory of materials. The problem that I wish to discuss in this page is the elastic stress and strain fields due to an ellipsoidal inclusion/inhomogeneity - a problem that was solved by Eshelby using an elegant thought experiment.</p>
<p id="exid"> In two papers published in the Proceedings of Royal Society (A) in 1957 and 1959 (Volume 241, p. 376 and Volume 252, p. 561) Eshelby solved the following problem ("with the help of a simple set of imaginary cutting, straining and welding operations"): In his own words, </p>
<ul><li> <strong> The transformation problem </strong> <br /> A region (the "inclusion") in an infinite homogeneous isotropic elastic medium undergoes a change of shape and size which, but for the constraint imposed by its surroundings (the "matrix"), would be an arbitrary homogeneous strain. What is the elastic state of inclusion and matrix? </li>
</ul><p id="exid"> The homogeneous strain is known as "eigenstrain" or "transformational strain". In the same paper, Eshelby also introduced the concept of "equivalent inclusion" for solving the transformation problem when the matrix and the region of eigenstrain (the "inhomogeneity") have different elastic constants. These two papers are very accessible and are a pleasure to read - and in case I forgot to mention, these two papers are a must read for anybody who is interested in theoretical materials science. </p>
<p> Thought experiment of Eshelby<br /></p><p id="exid">The operations that Eshelby used to solve the transformation problem are the following:</p>
<ul><li>Remove the region of interest from the matrix. </li>
<li> Allow it to take the eigenstrain.</li>
<li>Restore the region to its original shape and size by applying suitable surface tractions and put it back into the matrix and rejoin. </li>
<li>Remove the body force on the between the inclusion and matrix by applying an equal and opposite layer of body force.</li>
</ul><p>In step (3), the stress is zero in the matrix and is a known constant in the inclusion. The additional stress introduced in step (4) is found by the integration from the expression for the elastic field of a point force. </p>
<p id="exid"> The details of these calculations are found in the following (advanced) texts: <br /> (1) Rob Philips, <em> Crystals, defects and microstructures: Modeling across scales </em>, Cambridge University Press (2001), p. 520. <br /> (2) Toshio Mura, <em> Micromechanics of defects in solids</em>, Kluwer academic publishers (1987), p. 74. <br /> (3) A.G. Khachaturyan, <em> Theory of structural transformations in solids</em>, John Wiley &amp; Sons (1983), p. 198. <br /> While Rob philips and Mura invoke the elastic Green function to solve the transformation problem for a single inclusion, Khachaturyan describes the generalisation of the inclusion problem to multiple inclusions/eigenstrains.</p>
<p id="exid"> Eshelby (with his cuts, strains and weldings) has shown that the transformation problem is equivalent to solving for the equations of elastic equilibrium of a homogeneous body with a known body force distribution.</p>
<p> The Green function approach<br /></p><p id="exid"> For homogeneous bodies with known body force distribution, the equations of elastic equilibrium are solved using the elastic Green function. Rob Philips and Mura describe the Green function approach in great detail. </p>
<p id="exid"> In case you are interested in looking at some solutions of inclusion problems, Rob Philips describes the radial displacements associated with a spherical inclusion of radius "a" with dilatational eigenstrain obtained using Green functions (See the figure 10.14 on page 524) and indicates that the elastic energy of a spherical inclusion with dilatational misfit scales as the volume of the inclusion. Of course, Mura's classic is a singular testimony to the power of a combination of Green function and the eigenstrain approach of Eshelby. </p>
<p id="exid"> The Green function approach finally results in the evaluation of elliptic integrals for obtaining the displacements. This is not surprising since we are integrating the body forces over an elliptic geometry (Remember, the inclusions were ellipsoids/ellipses). The gradients of displacement give us the strain - That means we need to differentiate the Green functions. Thus, it is rather cumbersome to calculate the elastic stress, strain, or displacement fields using the Eshelby-Green approach.</p>
<p> The complex variable formalism<br /></p><p id="exid"> However, if we are interested in ellipses and not ellipsoids (that is, 2D problems), it is possible to avoid the cumbersome integrals. In 1960, in a paper published in the Proceedings of Cambridge Philosophical Society, (Volume 57, p. 669), Jaswon and Bhargava showed how to avoid the elliptic integrals using a complex variable formulation! The solution of Jaswon and Bhargava is based on the following ideas/results: </p>
<ul><li> Eshelby's method involves integrals of point forces on the matrix-inclusion boundary. </li>
<li>The expression for the displacement at any point <em> x </em> due to a point force <em> F </em> acting at the point <em> y </em> being known, the Eshelby problem now reduces to an integration of a continuous distribution of the forces over the inclusion surface. </li>
<li>Green and Zerna <em> (Theoretical elasticity)</em> and Mushkelishvili <em> (Some basic problems of the mathematical theory of elasticity) </em> give the expressions for the required mathematical arsenal for writing down the contour integrals!</li>
</ul><p id="exid"> Jaswon and Bhargava motivated their complex variable formulation by making the following observation: </p>
<ul><li> Although Eshelby has proved some general theorems of great interest, using elegant methods, his solutions involve analytically intractable integrals of a formidable nature. </li>
</ul><p> Eshelby himslef felt the same way, since in his 1959 paper he says: </p>
<ul><li> It has to be admitted that, except in the simplest cases, a calculation of the external field is laborious. </li>
</ul><p id="exid"> Jaswon and Bhargava built their complex variable formalism on the "ingeneous attack" on the transformation problem by Eshelby "utilizing the point-force concept". In view of the "novelty and importance" of the approach, they also give a brief description of Eshelby's arguements - And their description is by far one of the best that I have seen in the literature. Add this paper to your reading list, if you have not done it already! </p>
<p> 'Equivalent inclusion' for inhomogeneities<br /></p><p id="exid">The idea of "equivalent inclusion" is simple - It is one of those nice mathematical tricks where we solve a problem by reducing it to another which has already been solved. </p>
<p id="exid"> Let the inhomogeneity have an elastic constant that is different from that of the matrix. The idea is to replace the inhomogeneity with an inclusion - The eigenstrain in the hypothetical inclusion is such that it exerts the same actions on the matrix as the original inhomogeneity. Mathematically, finding out the eigenstrain in the equivalent inclusion amounts to solving for a set of three equations in three unknowns (in 2D) or six equations in six unknowns (in 3D). </p>
<p id="exid">There is some more information <a href="http://gururajan.mp.googlepages.com/inhomclusions">on the inclusion/inhomogeneity problem in my homepage here</a> (as you might have noticed, this post is a slightly modified version of the introduction page of that entry). Have fun! </p>
</div></div></div>Sun, 19 Nov 2006 12:45:52 +0000Mogadalai Gururajan456 at https://imechanica.orghttps://imechanica.org/node/456#commentshttps://imechanica.org/crss/node/456