Micromechanics with Mathematica (eBook)

Seiichi Nomura (Autor)

eBook Download: EPUB

2016
John Wiley & Sons (Verlag)
978-1-118-38570-8 (ISBN)

Infos

Demonstrates the simplicity and effectiveness of Mathematica as the solution to practical problems in composite materials

Designed for those who need to learn how micromechanical approaches can help understand the behaviour of bodies with voids, inclusions, defects, this book is perfect for readers without a programming background. Thoroughly introducing the concept of micromechanics, it helps readers assess the deformation of solids at a localized level and analyse a body with microstructures. The author approaches this analysis using the computer algebra system Mathematica, which facilitates complex index manipulations and mathematical expressions accurately.

The book begins by covering the general topics of continuum mechanics such as coordinate transformations, kinematics, stress, constitutive relationship and material symmetry. Mathematica programming is also introduced with accompanying examples. In the second half of the book, an analysis of heterogeneous materials with emphasis on composites is covered.

• Takes a practical approach by using Mathematica, one of the most popular programmes for symbolic computation
• Introduces the concept of micromechanics with worked-out examples using Mathematica code for ease of understanding
• Logically begins with the essentials of the topic, such as kinematics and stress, before moving to more advanced areas
• Applications covered include isotropic materials, plates and shells and thermal stress
• Features a problem and solution section on the book's companion website, useful for students new to the programme

Autor

Seiichi Nomura, The University of Texas at Arlington, USA Professor Nomura studied his BSc and MSc at the University of Tokyo, Japan before completing his PhD at the University of Delaware. He is now Professor in the Department of Mechanical and Aerospace Engineering at the University of Texas at Arlington, USA. Professor Nomura has undertaken research on the analysis of mechanical and thermal properties of heterogeneous materials including composite materials and functionally graded materials (FGMs) and has published extensively in this area.

Inhalt

Preface

1 Coordinate transformation and tensors 1

1.1 Index notation 1

1.2 Coordinate Transformations (Cartesian tensors) 10

1.3 Definition of tensors 12

References 20

2 Field Equations 21

2.1 Concept of Stress 21

2.1.1 Properties of stress 23

2.1.2 (Stress) Boundary conditions 25

2.1.3 Principal stresses 26

2.1.4 Stress deviator 30

2.1.5 Mohr's Circle 33

2.2 Strain 35

2.2.1 Shear deformation 42

2.3 Compatibility condition 43

2.4 Constitutive Relation, Isotropy, Anisotropy 45

2.4.1 Isotropy 46

2.4.2 Elastic modulus 48

2.4.3 Orthotropy 50

2.5 Constitutive relation for fluids 52

2.6 Derivation of field equations 53

2.6.1 Divergence theorem (Gauss theorem) 53

2.6.2 Material derivative 54

2.6.3 Equation of continuity 56

2.6.4 Equation of motion 56

2.6.5 Equation of energy 57

2.6.6 Isotropic solids 59

2.6.7 Isotropic fluids 59

2.6.8 Thermal effects 60

2.7 General coordinate system 61

2.7.1 Introduction to tensor analysis 61

2.7.2 Definition of tensors in curvilinear systems 62

References 71

3 Inclusions in infinite media 73

3.1 Eshelby's solution for an ellipsoidal inclusion problem 74

3.1.1 Eigenstrain problem 77

3.1.2 Eshelby tensors for an ellipsoidal inclusion 78

3.1.3 Inhomogeneity (inclusion) problem 87

3.2 Multi-layered inclusions 96

3.2.1 Background 96

3.2.2 Implementation of index manipulation in Mathematica 97

3.2.3 General formulation 100

3.2.4 Exact solution for two-phase materials 106

3.2.5 Exact solution for three-phase materials 113

3.2.6 Exact solution for four-phase materials 125

3.2.7 Exact solution for 2-D multi-phase materials 129

3.3 Thermal stress 130

3.3.1 Thermal stress due to heat source 131

3.3.2 Thermal stress due to heat flow 139

3.4 Airy's stress function approach 154

3.4.1 Airy's stress function 154

3.4.2 Mathematica programming of complex numbers 157

3.4.3 Multi-phase inclusion problems using Airy's stree function 161

3.5 Effective properties 173

3.5.1 Upper and lower bounds of effective properties 174

3.5.2 Self-consistent approximation 176

References 180

4 Inclusions in finite matrix 181

4.1 General Approaches for Numerically Solving Boundary Value Problems 182

4.1.1 Method of Weighted Residuals 182

4.1.2 Rayleigh-Ritz Method 194

4.1.3 Sturm-Liouville System 196

4.2 Steady-State Heat Conduction Equations 199

4.2.1 Derivation of permissible functions 200

4.2.2 Finding temperature field using permissible functions 213

4.3 Elastic Fields with Bounded Boundaries 218

4.4 Numerical Examples 224

4.4.1 Homogeneous medium 225

4.4.2 Single inclusion 226

References 235

Appendix

Textauszug

Chapter 1
Coordinate Transformation and Tensors

To describe the state of the deformation for a deformable body, the coordinate transformation plays an important rule, and the most appropriate way to represent the coordinate transformation is to use tensors. In this chapter, the concept of coordinate transformations and the introduction to tensor algebra in the Cartesian coordinate system are presented along with their implementation in Mathematica. As this book is not meant to be a textbook on continuum mechanics, the readers are referred to some good reference books including Romano et al. (2006) and Fung (1965), among others. Manipulation involving indices requires a considerable amount of algebra work when the expressions become lengthy and complicated. It is not practical to properly handle and evaluate quantities that involve tensor manipulations by conventional scientific/engineering software such as FORTRAN, C, and MATLAB. Software packages capable of handling symbolic manipulations include Mathematica (Wolfram 1999), Maple (Garvan 2001), and others. In this book, Mathematica is exclusively used for implementation and evaluation of derived formulas. A brief introduction to the basic commands in Mathematica is found in the appendix, which should be appropriate to understand and execute the Mathematica code used in this book.

1.1 Index Notation

If one wants to properly express the deformation state of deformable bodies regardless of whether they are solids or fluids, the use of tensor equations is essential. There are several different ways to denote notations of tensors, one of which uses indices and others without using indices at all. In this book, the index notation is exclusively used throughout to avert unnecessary abstraction at the expense of mathematical sophistication.

The following are the main compelling reasons to mandate the use of tensor notations in order to describe the deformation state of bodies correctly.

1. 1. The principle of physics stipulates that a physically meaningfulobject must be described independent of the frame of references.1 If the equation for a physically meaningful object changes depending on the coordinate system used, that equation is no longer a correct equation.
2. 2. Tensor equations can be shown to be invariant under the coordinate transformation. Tensor equations are thus defined as those equations that are unchanged from one coordinate system to another.

Hence, by combining the two aforementioned statements, it can be concluded that only tensor equations can describe the physical objects properly. In other words, if an equation is not in tensorial format, the equation does not represent the object physically.

The index notation, also known as the Einstein notation (Einstein et al. 1916)2 or the summation convention, is the most widely used notation to represent tensor quantities, which will be used in this book. The index notation in the Cartesian coordinate system is summarized as follows:

1. 1. For mathematical symbols that are referred to quantities in the , , and directions, use subscripts, 1, 2, 3, as in or , instead of or . The subscripted numbers 1, 2, and 3, refer to the , , and directions, respectively. Obviously, the upper limit of the number is 2 for 2-D and 3 for 3-D.
2. 2. If there are twice repeated indices in a term of products such as , the summation with respect to that index () is always assumed. For example,

   There is no exception to this rule. An expression such as is not allowed as the number of repetitions is 3 instead of 2.

A repeated index is called the dummy index as it does not matter what letter is used, and an unrepeated index is called the free index.3 For example,

all of which represent a summation . An unrepeated index such as (or or ) stands for one of , , or .

It should be noted that the notations and conventions introduced are valid for the Cartesian coordinate system only. In a curvilinear coordinate system such as the spherical coordinate system, the length of base vectors is not necessarily unity, and this mandates the aforementioned index notation to be modified to reflect the difference between the contravariant components and the covariant components, which will be discussed in Chapter 2.

1.1.1 Some Examples of Index Notation in 3-D

1. 1.

   As the index is repeated, the summation symbol, , must be added in front, i.e.,

   Note that is different from . While represents a single expression with three terms, represents one of the three expressions (, , or ).

2. 2.

   Note that the index is a dummy index (repeated twice) while the index is a free index (no repeat). Therefore,

3. 3.

   This is not a valid tensor expression as the number of repeated indices must be 2.

1.1.2 Mathematica Implementation

As Mathematica itself does not support tensor manipulation natively, it is necessary to devise a way to handle index notation and tensor manipulation. In this book, a list or a list of lists (a nested list) is used to represent tensor quantities. Using a nested list to define a tensor of any rank is straightforward but at the same time limited to the Cartesian tensors. For tensors defined in a curvilinear coordinate system, a slightly different approach is needed.

When running Mathematica first time, a default directory should be selected so that all the notebook files can be saved and accessed in this directory. By default, Mathematica looks for all the files stored in c:users<user> where <user> is the user's home directory.4The SetDirectory command can change this location. For example, if you want to change the default directory to c:tmp, the SetDirectory command can specify the default directory as

It is noted that the directory delimiter needs to be entered as “/” (forward slash) even though the Windows delimiter character is “” (backslash).

To enter a three-dimensional vector, , the following Mathematica command can be entered to create a list with braces (curly brackets) as

An individual component of can be referenced using double square brackets ([[]]) as

The partial derivative of with respect to can be entered as

You can also differentiate an individual component as

Implementation of the coordinate component, , into Mathematica can be done by using the Table function. To define a position vector, , whose components are , enter

The given Table function generates a list of elements. For example, the following command generates a sequence of for .

In the definition of the position vector, , it is noted that the coordinate components, , are entered as x[1], x[2], x[3] instead of x[[1]], x[[2]], x[[3]]. It is important to distinguish a single square bracket ([]) and a double square bracket ([[]]). The single square bracket ([]) is for a parameter used in a function. The quantity, x[1], stands for a function, x, with the argument of 1. By using a single square bracket, the quantities such as x[1], x[2] can stand for themselves, meaning that initial values do not have to be preassigned. On the other hand, if x[[1]] were used instead of x[1], 0 would be returned unless a list x is previously defined.

To define a function in Mathematica, use the following syntax:

In the aforementioned example code, a user-defined function, f[x], that returns is defined. The syntax is such that variables of the function must be presented to the left of the equal sign with the underscore and the definition of the function is given to the right of a colon and an equal sign (). It is important to note that a function in Mathematica returns itself if no prior definition is given.

In the aforementioned example, f[x, y] is defined as a function that takes two variables returning . However, when f is called with only one variable, 1, it returns itself, i.e., f[1], as f with only one variable has not been defined. When g[a] is entered, it returns itself without evaluation as there is no prior definition of g[x] given. It is this property of a function in Mathematica that enables manipulating index notations.

As an example of using x[i] as the coordinate components, here is how to implement the summation convention. As Mathematica does not have support for the summation convention built in, if x[i]x[i]meant as is entered as

Mathematica does not automatically expand and reduce the result to . Hence, it is necessary to explicitly use the Sum[] command as

for .

However, it is possible to implement the summation convention in Mathematica with the following procedure.

The aforementioned three-line code adds a new rule to automatically replace x[i]x[i] by through pattern matching. When Mathematica evaluates the product of two quantities, it calls its internal function, Times, which tries to simplify...

Produkt-Details

• Erscheint lt. Verlag: 2.3.2016
• Sprache: englisch
• Themenwelt: Mathematik / Informatik ► Mathematik ► Computerprogramme / Computeralgebra
• Themenwelt: Naturwissenschaften ► Physik / Astronomie ► Mechanik
• Themenwelt: Technik ► Maschinenbau
• Schlagworte: applied mechnanics • composite materials</p> • Composites • Computational / Numerical Methods • Continuum Mechanics • Deformation • Festkörpermechanik • Festkörpermechanik • isotropic materials • Kinematics • <p>micromechanics • Maschinenbau • Materials Science • Materialwissenschaften • Mathematica • mechanical engineering • Mikromechanik • Plates • Rechnergestützte / Numerische Verfahren im Maschinenbau • Rechnergestützte / Numerische Verfahren im Maschinenbau • shells thermal stress • solid mechanics • solids • Verbundwerkstoffe
• ISBN-10: 1-118-38570-5 / 1118385705
• ISBN-13: 978-1-118-38570-8 / 9781118385708