On the M-Eigenvalues of Elasticity Tensor and the Strong Ellipticity

Total Page:16

File Type:pdf, Size:1020Kb

On the M-Eigenvalues of Elasticity Tensor and the Strong Ellipticity On the M-eigenvalues of elasticity tensor and the strong ellipticity condition Hua Xianga∗ Liqun Qib† Yimin Weic‡ a School of Mathematics and Statistics, Wuhan University, Wuhan, 430072, P.R. China b Department of Applied Mathematics, The Hong Kong Polytechnic University, Hong Kong c School of Mathematical Sciences and Shanghai Key Laboratory of Contemporary Applied Mathematics, Fudan University, Shanghai, 200433, P. R. of China January 23, 2018 Abstract Strong ellipticity is an important property in the elasticity theory. In 2009, M-eigenvalues were introduced for the elasticity tensor. It was shown that M-eigenvalues are invariant under coordinate system choices, and the strong ellipticity condition holds if and only if all the M-eigenvalues of the elasticity tensor are positive. Thus, M-eigenvalues are some in- trinsic parameters of the elasticity tensor. In this paper, we show that the M-eigenvalues of the elasticity tensor are closely related with some elastic moduli, such as the bulk modulus, the shear modulus, Lam´e’s first parameter, the P-wave modulus, etc, and the positiveness of the M- eigenvalues are corresponding to some existing conditions for strong ellip- ticity in some special cases, such as the isotropic case, the cubic case, the polar anisotropic case and the tetragonal case. We also present new suf- ficient conditions for the strong ellipticity of the orthotropic case. These, in a certain sense, further reveal the physical meanings of M-eigenvalues. Keywords. M-eigenvalue, Z-eigenvalue, Elasticity tensor, Strong ellipticity, Positive definiteness, Elastic modulus arXiv:1708.04876v2 [math.RA] 22 Jan 2018 1 Introduction The elasticity tensor is a fourth-rank tensor, namely a four-linear map : 3 3 3 3 C C R R R R R. With respect to the basis ∂i in a local coordinate × × × → { } system, the components of is expressed by cijkl = (∂i, ∂j , ∂k, ∂l), where C C ∗E-mail: [email protected]. H. Xiang is supported by the National Natural Science Foundation of China under grants 11571265, 11471253 and NSFC-RGC No.11661161017. †E-mail: [email protected]. L. Qi is supported by the Hong Kong Research Grant Council (Grant No. PolyU 15302114, 15300715, 15301716 and 15300717). ‡E-mail: [email protected]. Y. Wei is supported by the National Natural Science Foun- dation of China under grant 11771099 and International Cooperation Project of Shanghai Municipal Science and Technology Commission under grant 16510711200. 1 i,j,k,l 1, 2, 3 . With respect to the dual basis dxi , the elasticity tensor ∈ { } i j k l { } can be written as = cijkldx dx dx dx . Using these components, we can reformulate the tensorC as C i j k l (w,x,y,z)= wxyz = cijklw x y z , C C i i i i where w = w ∂i, x = x ∂i, y = y ∂i, z = z ∂i, and we adopt the Einstein convention of summing over the repeated indices. When working on the or- thonormal Cartesian coordinate system with a Euclidean metric given by the Kronecker delta, the covariant and contravariant types do not differ, and the indices can be upper or lower, and hence we have the option to work with only subscripts in the remainder of this article. For the fourth-order elasticity or stiffness tensor , the strong ellipticity is defined by the following function to be positive. C 3 f(x, y)= xyxy = Cijklxiyjxkyl > 0, x, y R . C ∀ ∈ Such strong ellipticity condition ensures that the governing differential equations for elastostatics problems be completely elliptic. It is an important property in the elasticity theory associated with uniqueness, instability, wave propagation, etc, and has been studied extensively [2, 3, 8, 9, 10, 11, 16, 17, 21, 22, 23, 25, 26, 27, 28, 31, 32, 33]. In 2009, Qi, Dai and Han [24] introduced the M-eigenvalues for the elas- ticity tensor. Here, the letter “M” means mechanics. It was shown that the strong ellipticity condition holds if and only if the optimal value of a global polynomial optimization problem is positive, which further naturally leads to the M-eigenvalue problem from the optimality condition. The M-eigenvalues θ of the elasticity tensor are defined as follows. C yxy = θx, (1) C· xyx = θy, (2) C · T T under the constraints x x = 1 and y y = 1. Here, ( yxy)i = Cijklyjxkyl and C· ( xyx )l = Cijklxiyjxk in the Einstein summation notation. It was shown that M-eigenvaluesC · are invariant under coordinate system choices, and the strong ellipticity condition holds if and only if all the M-eigenvalue of the elasticity tensor are positive [24]. Thus, M-eigenvalues are some intrinsic parameters of the elasticityC tensor. The M-eigenvalue problem can be reduced to the Z-eigenvalue problem. x3 = ηx, (3) C where xT x = 1, and η is the Z-eigenvalue. As the M-eigenvalue is associated with the strong ellipticity, the Z-eigenvalue is related to the rank-one positive definiteness of . All the Z-eigenvalues are M-eigenvalues of , but not vice versa. Obviously, theC positiveness of Z-eigenvalues provides the necessaC ry conditions for the strong ellipticity of material. We note that there exists the major symmetry Cijkl = Cklij , and the mi- nor symmetry Cijkl = Cjikl and Cijkl = Cijlk . The requirement of the sym- metry reduces the number of different elements to 21. Checking the invari- ance of under a specific transformation Q can further reduces this num- ber for someC special cases. If the material response is to be the same, i.e., 2 Cijkl = QimQjnQkpQlqCmnpq, then it yields the material symmetry relation, which will provide a system of equations that allows reduction in the number of independent elastic moduli [20, 30]. For example, as to an orthotropic material, it has three mutually perpendicular symmetry planes, and this type of material has nine independent material constants. Common examples of such materials include wood and fiber-reinforced composites. Omitting the zero elements, the M-eigenvalue problem (1)-(2) for this case reads C1111y1x1y1 + (C1122 + C1221)y1x2y2 + (C1133 + C1331)y1 x3y3 + C1313 y3x1y3 + C1212 y2x1y2 = θx1, (C2211 + C2112)y2x1y1 + C2222 y2x2y2 + (C2233 + C2332)y2 x3y3 + C2323 y3x2y3 + C2121 y1x2y1 = θx2, (C3311 + C3113)y3x1y1 + (C3322 + C3223)y3 x2y2 + C3333 y3x3y3 + C3232 y2x3y2 + C3131 y1x3y1 = θx3, C1111 x1y1x1 + (C2211 + C1221 )x2y2x1 + (C3311 + C1331 )x3y3x1 + C3131 x3y1x3 + C2121x2y1x2 = θy1 , (C1122 + C2112)x1y1x2 + C2222x2y2x2 + (C3322 + C2332 )x3y3x2 + C3232 x3y2x3 + C1212x1y2x1 = θy2 , (C1133 + C3113)x1y1x3 + (C2233 + C3223)x2y2x3 + C3333x3y3x3 + C2323 x2y3x2 + C1313x1y3x1 = θy3 . To keep accordance with the standard shorthand notation in the theory of elasticity, we need the Voigt notation or the Kelvin notation. Since the minor symmetries are valid, the fourth-order stiffness tensor Cijkl may be written as a6 6 stiffness matrix (cαβ) (a tensor of second order). Due to the major sym- metry,× this matrix is also symmetric. In the Voigt notation, each pair of indices is compressed into one index, and the standard mapping for tensor indices reads as follows. ij = 11 22 33 23, 32 13, 31 12, 21 α⇓ = ⇓⇓⇓⇓1234 ⇓5 ⇓6 For example, in a properly chosen coordinate system, the Voigt matrix of the orthotropic material has the following form. c11 c12 c13 0 0 0 c12 c22 c23 0 0 0 c13 c23 c33 0 0 0 (cαβ )= . 0 0 0 c44 0 0 0 0 0 0 c 0 55 0 0 0 0 0 c66 The elasticity tensor is described by nine parameters: c11 = C1111, c22 = C2222, c33 = C3333, c23 = C2233, c13 = C1133, c12 = C1122, c44 = C2323 = C3232, c55 = C1313 = C3131, and c66 = C1212 = C2121. Considering the symmetries in the elasticity tensor , we naturally introduce the following eigenproblem. C Cijklxkl = ζ xij , (4) where the eigentensor xij is symmetric. Such eigenproblem is closely related to the positive definiteness of , a less general hypothesis than the strong ellipticity, which also guarantees theC uniqueness of solutions in problems of elasticity [13, P.159]. The elasticity stiffness tensor must be positive definite [30, P.48][29, P.34]. That is, The strain energy densityC or elastic potential satisfies 1 C ǫ ǫ > 0, 2 ijkl ij kl 3 where ǫij is any symmetric strain tensor. It physically means that energy is needed to deform an elastic body from its unloaded equilibrium position [14]. The condition for the positiveness of strain energy can be replaced by the less restrictive condition Cijklγij γkl > 0, where the γij need not be symmetric [29, P.34]. Considering the symmetry of C, we rewrite (4) in its equivalent matrix form. C1111 C1122 C1133 √2C1123 √2C1113 √2C1112 x11 x11 C2211 C2222 C2233 √2C2223 √2C2213 √2C2212 x22 x22 C C C √2C √2C √2C x33 x33 3311 3322 3333 3323 3313 3312 = ζ . √2C √2C √2C 2C 2C 2C √2x √2x 2311 2322 2333 2323 2313 2312 23 23 √2C √2C √2C 2C 2C 2C √2x √2x 1311 1322 1333 1323 1313 1312 13 13 √ √ √2C1211 √2C1222 √2C1233 2C1223 2C1213 2C1212 2x12 2x12 The positiveness condition corresponds to impose that all the eigenvalues ζ are positive [30]. What’s more, the positive definiteness of implies strong ellipticity, while the converse statement is not true [13, P.159].C In other words, the positive definiteness conditions are sufficient for the strong ellipticity. In the following, we consider five further simplified cases of the orthotropic material: the isotropic case, the cubic case, the polar anisotropic case, the tetragonal case and the orthotropic case, which have two, three, five, six and nine independent elastic parameters respectively. In this paper, we will show that the M-eigenvalues are closely related with some elastic moduli, such as the bulk modulus, the shear modulus, Lam´e’s first parameter, the P-wave modulus, etc. In fact, all the elastic constants in the diagonal of stiffness matrix (cαβ) are M-eigenvalues.
Recommended publications
  • An Approximate Nodal Is Developed to Calculate the Change of •Laatio Constants Induced by Point Defect* in Hep Metals
    1 - INTBOPUCTIOB the elastic conatants.aa well aa othernmechanieal pro partita of* IC/79/lW irradiated materiale (are vary sensitive to tne oonoantration of i- INTERNAL REPORT (Limited distribution) rradiation produced point defeota.One of the firat eatimatee of thla effect waa done by Dienea [l"J who aiaply averaged over the whole la- International Atomic Energy Agency ttice the locally changed interatoalo bonds due to the pxesense of" and the defect.With tola nodal ha predioted an inoreaee of the alaatle United Nations Educational Scientific and Cultural Organization constant* of about 10* par atonic f of interatitlala in Ou and a da. oreaee of l]t par at. % of vacanolea. INTERNATIONAL CENTRE FOE THEORETICAL PHYSICS Later on,experlaental etudiea by Konlg at al.[2]and wanal lilgar* very large decrease a of about 5O}( per at.)t of Vrenlcal dafaota.fha theory waa than iaproved in order to relate the change of a la* tic con a tan t a to the defect lnduoed change of force oonatanta and tno equivalent aethoda were devalopedi the energy-»athod of ludwig [4] CHANGE OF ELASTIC CONSTANTS and the t-aatrix method of Slllot et al.C ?]. INDUCED BY FOIMT DEFECTS IN hep CRYSTALS * Iheoretioal eetlnates for oublo oryetala have been oarrie* out by Ludwig[4] for the caae of vacancleafby Piatorlueld for intarati- Carlos Tome •• tiala and by Sederloaa et al.t?] for duabell interotitiala. International Centre for Theoretical Physics, Trieste, Italy. Re thaoretloal work baa been done ao far for hexagonal cryatmlaj and the experimental neaeurentanta (available only for Kg) ax* eona- ABSTRACT what crude t8,9i,ayan though in the last few.
    [Show full text]
  • Gravitation in the Surface Tension Model of Spacetime
    IARD 2018 IOP Publishing IOP Conf. Series: Journal of Physics: Conf. Series 1239 (2019) 012010 doi:10.1088/1742-6596/1239/1/012010 Gravitation in the surface tension model of spacetime H A Perko1 1Office 14, 140 E. 4th Street, Loveland, CO, USA 80537 E-mail: [email protected] Abstract. A mechanical model of spacetime was introduced at a prior conference for describing perturbations of stress, strain, and displacement within a spacetime exhibiting surface tension. In the prior work, equations governing spacetime dynamics described by the model show some similarities to fundamental equations of quantum mechanics. Similarities were identified between the model and equations of Klein-Gordon, Schrödinger, Heisenberg, and Weyl. The introduction did not explain how gravitation arises within the model. In this talk, the model will be summarized, corrected, and extended for comparison with general relativity. An anisotropic elastic tensor is proposed as a constitutive relation between stress energy and curvature instead of the traditional Einstein constant. Such a relation permits spatial geometric terms in the mechanical model to resemble quantum mechanics while temporal terms and the overall structure of tensor equations remain consistent with general relativity. This work is in its infancy; next steps are to show how the anisotropic tensor affects cosmological predictions and to further explore if geometry and quantum mechanics can be related in more than just appearance. 1. Introduction The focus of this research is to find a mechanism by which spacetime might curl, warp, or re-configure at small scales to provide a geometrical explanation for quantum mechanics while remaining consistent with gravity and general relativity.
    [Show full text]
  • Arxiv:1704.01012V1 [Cond-Mat.Mtrl-Sci] 1 Apr 2017 Is the Second Most Important Material After Silicon
    Symmetry and Piezoelectricity: Evaluation of α-Quartz coefficients C. Tannous Laboratoire des Sciences et Techniques de l'Information, de la Communication et de la Connaissance, UMR-6285 CNRS, Brest Cedex3, FRANCE Piezoelectric coefficients of α-Quartz are derived from symmetry arguments based on Neumann's Principle with three different methods: Fumi, Landau-Lifshitz and Royer-Dieulesaint. While Fumi method is tedious and Landau-Lifshitz requires additional physical principles to evaluate the piezo- electric coefficients, Royer-Dieulesaint is the most elegant and most efficient of the three techniques. PACS numbers: 77.65.-j, 77.65.Bn, 77.84.-s Keywords: Piezoelectricity, piezoelectric constants, piezoelectric materials I. INTRODUCTION AND MOTIVATION Physics students are exposed to various types of symmetry [1] and conservation laws in Graduate/Undergraduate Mechanics and Electromagnetism with Lorentz transformation and Gauge symmetries, in Graduate/Undergraduate Quantum Mechanics during the study of Atoms and Molecules. In undergraduate courses such as Special Relativity, Lorentz Transformation is used to unify symmetries between Mechanics and Electromagnetism. In Graduate High Energy Physics, the CPT theorem where C denotes charge conjugation (Q ! −Q), P is parity (r ! −r) and T is time reversal (t ! −t) as well as Gauge symmetry (Ai ! Ai + @iχ) provide an important insight into the role of symmetry in the building blocks of matter and unification of fundamental forces and interaction between particles. Graduate/undergraduate Solid State Physics provide a direct illustration of how Crystal Symmetry plays a fun- damental role in the determination of physical constants and transport coefficients as well as conservation and sim- plification of physical laws. The relation between symmetry and dispersion relations through Kramers theorem (T symmetry) is another example of the power of symmetry in Solid State physics.
    [Show full text]
  • A Manifold Learning Approach to Data-Driven Computational Mechanics
    A Manifold Learning Approach to Data-Driven Computational Mechanics Ruben Ibanez, Emmanuelle Abisset-Chavanne, Jose Vicente Aguado, David Gonzalez, Elías Cueto, Francisco Chinesta To cite this version: Ruben Ibanez, Emmanuelle Abisset-Chavanne, Jose Vicente Aguado, David Gonzalez, Elías Cueto, et al.. A Manifold Learning Approach to Data-Driven Computational Mechanics. 13e colloque national en calcul des structures, Université Paris-Saclay, May 2017, Giens, Var, France. hal-01926477 HAL Id: hal-01926477 https://hal.archives-ouvertes.fr/hal-01926477 Submitted on 19 Nov 2018 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. CSMA 2017 13ème Colloque National en Calcul des Structures 15-19 Mai 2017, Presqu’île de Giens (Var) A Manifold Learning Approach to Data-Driven Computational Me- chanics R. Ibañez1, E. Abisset-Chavanne1, J.V. Aguado1, D. Gonzalez2, E. Cueto2, F. Chinesta1 1 ICI Institute, Ecole Centrale Nantes, {Ruben.Ibanez-Pinillo,Emmanuelle.Abisset-Chavanne,Jose.Aguado-Lopez,Francisco.Chinesta}@ec-nantes.fr 2 Aragon Institute of Engineering Research, Universidad de Zaragoza, Spain, {gonzal;ecueto}@unizar.es Résumé — Standard simulation in classical mechanics is based on the use of two very different types of equations.
    [Show full text]
  • Irreducible Matrix Resolution of the Elasticity Tensor for Symmetry Systems
    Irreducible matrix resolution of the elasticity tensor for symmetry systems Yakov Itin Inst. Mathematics, Hebrew Univ. of Jerusalem, Givat Ram, Jerusalem 91904, Israel, and Jerusalem College of Technology, Jerusalem 91160, Israel, email: [email protected] December 11, 2018, file ElastMatrix1.tex Abstract In linear elasticity, a fourth order elasticity (stiffness) tensor of 21 independent com- ponents completely describes deformation properties of a material. Due to Voigt, this tensor is conventionally represented by a 6 × 6 symmetric matrix. This classical matrix representation does not conform with the irreducible decomposition of the elasticity tensor. In this paper, we construct two alternative matrix representations. The 3 × 7 matrix representation is in a correspondence with the permutation transformations of indices and with the general linear transformation of the basis. An additional repre- sentation of the elasticity tensor by three 3 × 3 matrices is suitable for description the irreducible decomposition under the rotation transformations. We present the elasticity tensor of all crystal systems in these compact matrix forms and construct the hierarchy diagrams based on this representation. Key index words: anisotropic elasticity tensor, irreducible decomposition, matrix repre- sentation. Mathematics Subject Classification: 74B05, 15A69, 15A72. Contents arXiv:1812.03367v1 [physics.class-ph] 8 Dec 2018 1 Introduction2 2 Matrix representations3 2.1 Voigt’s representation . .4 2.2 An alternative representation . .4 1 Y. Itin Irreducible matrix resolution of the elasticity tensor 2 3 Irreducible decomposition5 3.1 GL(3)-decomposition . .5 3.2 SO(3; R)-decomposition . .8 3.2.1 Cauchy part . .8 3.2.2 Non-Cauchy part . 10 3.3 Irreducible decomposition .
    [Show full text]
  • APPLIED ELASTICITY, 2Nd Edition Matrix and Tensor Analysis of Elastic Continua
    APPLIED ELASTICITY, 2nd Edition Matrix and Tensor Analysis of Elastic Continua Talking of education, "People have now a-days" (said he) "got a strange opinion that every thing should be taught by lectures. Now, I cannot see that lectures can do so much good as reading the books from which the lectures are taken. I know nothing that can be best taught by lectures, expect where experiments are to be shewn. You may teach chymestry by lectures. — You might teach making of shoes by lectures.' " James Boswell: Lifeof Samuel Johnson, 1766 [1709-1784] ABOUT THE AUTHOR In 1947 John D Renton was admitted to a Reserved Place (entitling him to free tuition) at King Edward's School in Edgbaston, Birmingham which was then a Grammar School. After six years there, followed by two doing National Service in the RAF, he became an undergraduate in Civil Engineering at Birmingham University, and obtained First Class Honours in 1958. He then became a research student of Dr A H Chilver (now Lord Chilver) working on the stability of space frames at Fitzwilliam House, Cambridge. Part of the research involved writing the first computer program for analysing three-dimensional structures, which was used by the consultants Ove Arup in their design project for the roof of the Sydney Opera House. He won a Research Fellowship at St John's College Cambridge in 1961, from where he moved to Oxford University to take up a teaching post at the Department of Engineering Science in 1963. This was followed by a Tutorial Fellowship to St Catherine's College in 1966.
    [Show full text]
  • Cambridge University Press 978-0-521-64250-7 — Group Theory with Applications in Chemical Physics Patrick Jacobs Index More Information
    Cambridge University Press 978-0-521-64250-7 — Group Theory with Applications in Chemical Physics Patrick Jacobs Index More Information Index Abelian group 2 central extension 336, 337, 367 acoustic mode 393 centralizer 14, 19, 434 active representation 23 centre 19 adiabatic potential 173 character 74, 99 adjoint 54 character system 74 of an operator 102 character tables 76–78, 80, 447 adjoint matrix 418 character vector 259 algebra of turns 228 characteristic equation 420, 441 ambivalent class 435 charge overlap 107 angular momentum 184, 189 charge transfer 178 anharmonicity 160 chemical bond 106 constant of 160 class 5, 19 antibonding orbitals 106, 125 class algebra 439 antiferromagnetic crystal 265 class constants 436 antilinear operator 252 class property 440 antipole 223 Clebsch–Gordan series 209, 277, 385 antisymmetrical direct product 100 closed shell 172 À2 antisymmetrizing operator 141 closo Bn Hn 51 antiunitary operator 252, 267, 405 closure 1, 393 associated Legendre functions 194 co-factor 413 associative 2, 220 coincidence 162 axial groups 82 column matrix 415 axial tensor 283 combination bands 160 axial vector 82 commutation relations (CRs) 131, 187 compatibility relations 362 basic domain 331 complementary IR 303 basis 53, 96 complementary minor 413 of a lattice 308 complementary operators 265 basis functions, construction of 97 complex conjugate 218 bcc see body-centred cubic complex conjugation operator 253 benzene 104, 109, 174 complex number 218 Bethe 80, 150, 151 complex plane 219 bilateral binary (BB) 232 complex
    [Show full text]
  • Continuum Mechanics Using Mathematica® Fundamentals, Methods, and Applications
    Modeling and Simulation in Science, Engineering and Technology Antonio Romano Addolorata Marasco Continuum Mechanics using Mathematica® Fundamentals, Methods, and Applications Second Edition Modeling and Simulation in Science, Engineering and Technology Series Editor Nicola Bellomo Politecnico di Torino Tor ino, Italy Editorial Advisory Board K.J. Bathe P. Koumoutsakos Department of Mechanical Engineering Computational Science & Engineering Massachusetts Institute of Technology Laboratory Cambridge, MA, USA ETH Zürich Zürich, Switzerland M. Chaplain Division of Mathematics H.G. Othmer University of Dundee Department of Mathematics Dundee, Scotland, UK University of Minnesota Minneapolis, MN, USA P. Degond Department of Mathematics, K.R. Rajagopal Imperial College London, Department of Mechanical Engineering London, United Kingdom Texas A&M University College Station, TX, USA A. Deutsch Center for Information Services T.E. Tezduyar and High-Performance Computing Department of Mechanical Engineering & Technische Universität Dresden Materials Science Dresden, Germany Rice University Houston, TX, USA M.A. Herrero Departamento de Matematica Aplicada A. Tosin Universidad Complutense de Madrid Istituto per le Applicazioni del Calcolo Madrid, Spain “M. Picone” Consiglio Nazionale delle Ricerche Roma, Italy More information about this series at http://www.springer.com/series/4960 Antonio Romano • Addolorata Marasco Continuum Mechanics using MathematicaR Fundamentals, Methods, and Applications Second Edition Antonio Romano Addolorata Marasco Department
    [Show full text]
  • Valley Degeneracy in Biaxially Strained Aluminum Arsenide Quantum Wells
    PHYSICAL REVIEW B 84, 125319 (2011) Valley degeneracy in biaxially strained aluminum arsenide quantum wells S. Prabhu-Gaunkar,1 S. Birner,2 S. Dasgupta,2 C. Knaak,2 and M. Grayson1 1Electrical Engineering and Computer Science, Northwestern University, Evanston, Illinois 60208, USA 2Walter Schottky Institut, Technische Universitat¨ Munchen,¨ Garching DE-85748, Germany (Received 13 September 2010; revised manuscript received 19 May 2011; published 16 September 2011) This paper describes a complete analytical formalism for calculating electron subband energy and degeneracy in strained multivalley quantum wells grown along any orientation with explicit results for AlAs quantum wells (QWs). In analogy to the spin index, the valley degree of freedom is justified as a pseudospin index due to the vanishing intervalley exchange integral. A standardized coordinate transformation matrix is defined to transform between the conventional-cubic-cell basis and the QW transport basis whereby effective mass tensors, valley vectors, strain matrices, anisotropic strain ratios, piezoelectric fields, and scattering vectors are all defined in their respective bases. The specific cases of (001)-, (110)-, and (111)-oriented aluminum arsenide (AlAs) QWs are examined, as is the unconventional (411) facet, which is of particular importance in AlAs literature. Calculations of electron confinement and strain for the (001), (110), and (411) facets determine the critical well width for crossover from double- to single-valley degeneracy in each system. The biaxial Poisson ratio is calculated for the high-symmetry lower Miller index (001)-, (110)-, and (111)-oriented QWs. An additional shear-strain component arises in the higher Miller index (411)-oriented QWs and we define and solve for a shear-to-biaxial strain ratio.
    [Show full text]
  • Graphical Representation of the Generalized Hooke's Law Izei®Ei H
    TECHNISCHE MECHANIK, Band 21, Heft 2, (2001), 145-158 Manuskripteingang: 30. Januar 2001 Graphical Representation of the Generalized Hooke’s Law T. Böhlke, C. Brüggemann The anisotropic linear elastic behavior of single crystals can be described equivalently by a 4th-0rder elasticity tensor or two functions E(d) and K These functions represent Young’s modulus and a generalized bulk modulus as functions of the tensile direction d in a tension test. In the present paper three— and two—dimensional graphical representations of Young’s modulus and the generalized bulk modulus are given for single crystals belonging to one of the following symmetry groups: monoclinic, rhombic, trigonal, tetragonal, hemagonal, and cubic symmetry. 1 Introduction The generalized Hooke’s law is the geometrical and physical linear relation between stress and strain of anisotropic elastic solids. It is one of the oldest and best known constitutive relations in continuum mechanics. One dimensional formulations have been discussed first in the 17th century. Over a long time the scientific community agreed to differ about the most general form in linear elasticity. Finally Voigt (1850-1919) answered the open question how many constants have to be determined in the isotropic and anisotropic case. He proved by experiments Green’s (1793-1841) hypothesis that in the isotropic case two and in the general anisotropic hyperelastic case 21 constants have to be determined. Based on Hooke’s law structures are analyzed in all branches of mechanical and civil engineering. In most cases the isotropic version of Hooke’s law is used. But in the last decades more and more anisotropies have been taken into account, which are especially important if the component parts are laminated or made from single crystals.
    [Show full text]
  • Modeling and Characterization of the Elastic Behavior Of
    MODELING AND CHARACTERIZATION OF THE ELASTIC BEHAVIOR OF INTERFACES IN NANOSTRUCTURED MATERIALS: FROM AN ATOMISTIC DESCRIPTION TO A CONTINUUM APPROACH A Dissertation Presented to The Academic Faculty by Rémi Dingreville In Partial Fulfillment of the Requirements for the Degree Doctor of Philosophy in the School of Mechanical Engineering Georgia Institute of Technology December 2007 MODELING AND CHARACTERIZATION OF THE ELASTIC BEHAVIOR OF INTERFACES IN NANOSTRUCTURED MATERIALS: FROM AN ATOMISTIC DESCRIPTION TO A CONTINUUM APPROACH Approved by: Dr. Jianmin Qu, Advisor Dr. Mo Li GWW School of Mechanical Engineering School of Materials Science & Engineering Georgia Institute of Technology Georgia Institute of Technology Dr. David L. McDowell Dr. Elisa Riedo GWW School of Mechanical Engineering School of Physics Georgia Institute of Technology Georgia Institute of Technology Dr. Min Zhou GWW School of Mechanical Engineering Georgia Institute of Technology Date Approved: July 24 th , 2007 To my wife Stephani, my family and my friends. ACKNOWLEDGEMENTS This work is the accomplishment of both an individual and collective effort to find a path in the world of “nano ” and the meanders of interface theory. If one had to find the roots of this journey, I would probably cite two particular events. The first one would be the reading of an essay by Gleiter on nanocrystalline materials in Progress in Materials Science during my Masters study in France. I discovered a whole new world made of nanotubes and nanoribbons, a world ruled by abnormal behaviors. It seemed to me there were many interesting things to discover and issues to tackle. The second event occurred a few months earlier during an internship at Georgia Tech in the summer of 2000.
    [Show full text]
  • V “Type” Notation and Voigt Or Mandel Matrices
    T June 6, 2015 9:08 pm A F Tensors R a n n o n D c a B r R e b e c 1 n 2 Vm “type” notation and Voigt or Mandel matrices 3 4 This section defines a special “Vn ” tensor classification notation. Any tensor of type 5 m 6 n n n 7 Vm will be seen to have m components. Furthermore, any tensor of type Vm is also of 8 1 n 9 type VN where N= m . This latter identification helps us to convert indicial formulas 10 11 involving 3u 3 tensor matrices into matrix formulas involving 9u 1 arrays. Similarly, 12 13 indicial formulas involving 3 3 3 third-order tensor matrices may be alternatively cast 14 u u 15 in terms of 3 9 , 9 3 or even 27 1 arrays, depending on what is most convenient for 16 u u u 17 18 actual computations using linear algebra libraries in modern computing environments. 19 20 th st 21 Scalars are often called 0 -order tensors. Vectors are sometimes called 1 -order ten- 22 th n * 23 sors. In general, an n order engineering tensor has 3 components, and we say that 24 n 2 25 these tensors are of type V3 . For example, stress is a second-order tensor, so it has 3 26 27 (nine) components. Stress is symmetric, so these components are not all independent, but 28 29 that doesn’t change the fact that there are still nine components. 30 31 32 When solving a problem for which all tensors have symmetry with respect to some 2D 33 34 35 plane embedded in 3D space, it is conventional to set up the basis so that the third basis 36 37 vector points perpendicular to that plane.
    [Show full text]