Features of the free and open source toolbox MTEX for texture analysis
|
|
- Basil Patrick
- 5 years ago
- Views:
Transcription
1 Features of the free and open source toolbox MTEX for texture analysis Ralf Hielscher 1 Florian Bachmann 2 Helmut Schaeben 3 David Mainprice 4 1 Applied Functional Analysis, Technische Universität Chemnitz, Germany 2 Geoscience Mathematics and Informatics, Technische Universität Bergakademie Freiberg, Germany 3 Géosciences Montpellier UMR 5242, Université Montpellier 2 & CNRS, France Internal MTEX Training Géosciences Montpellier, September 2012
2 Outline 1 Introduction MTEX project MATLAB environment 2 EBSD Data Analysis Importing EBSD Data 3 Visualization Grains Analysis EBSD - to - ODF Reconstruction
3 History of MTEX : A Texture Calculation Toolbox 1 Originally developed by Ralf Hielscher as part of this Ph.D. thesis The Radon Transform on the Rotation Group : Inversion and Application to Texture Analysis (2007) under the direction of Helmut Schaeben at Technische Universität Bergakademie Freiberg, Germany. 2 The toolbox has been extended to include the analysis of 2D and 3D EBSD data by Florian Bachmann and others. 3 The physical properties of 2 nd, 3 rd and 4 th rank tensors have been incorporated by David Mainprice and others.
4 Features of MTEX crystal geometry: all kind of symmetries, different Euler angle conventions, import from crystallographic information files (CIF); pole figure data analysis: 20 data formats, data correction; pole figure to ODF inversion: all symmetries, any sampling of pole figures including irregular grids, incomplete and non-normalized pole figures, ghost correction option, zero range option; EBSD data analysis: 5 data formats, data correction; grain detection: fabric analysis, misorientation analysis; ODF estimation from EBSD data: automized determination of kernel width, arbitrarily many individual orientations; ODF analysis: modal orientations, difference ODFs, volume portions, entropy, texture index, Fourier coefficients; ODF modeling: any composition of uniform, unimodal, fibre and Bingham ODFs, simulation of pole figure and individual orientation data; material property tensors: average tensors from EBSD data and ODFs; elasticity tensors: elastic stiffness tensor, elastic compliance tensor, Young s modulus, shear modulus, Poisson s ratio, linear compressibility, compressional and shear elastic wave velocities, velocity ratios.
5 MTEX : The Concept ODF Estimation EBSD Data Analysis Plotting MTEX Pole Figure Analysis Texture Characteristics ODF Modelling MATLAB toolbox Fourier and numerical methods on SO3 Pole figure inversion route to ODF Individual orientation route to ODF Spatially referenced data in 2D or 3D EBSD... > 1000 functions > 20 data formats imported Publication ready plots Batch processing FREE! for everyone
6 MTEX : Webpage an import source of information
7 MTEX : Not so frequently asked questions
8 MTEX : Issues - a contact with the development team
9 MTEX : Central role of Orientation Distribution Function (ODF)
10 MTEX : Analyzing and Visualizing ODFs
11 Overview of MATLAB 1 High-level (line by line interpreter) language for technical computing 2 A large number of mathematical functions integrated in the language for matrix algebra, statistics, Fourier analysis, D and 3-D graphics functions for visualizing data 4 Extensions using commercial or free open source toolboxes (e.g. MTEX ) 5 Functions for integrating MATLAB based algorithms with external languages, such as C, C++, FORTRAN...
12 Vectors and Matrices Define vectors: v = [ 1 ; 2 ; 3 ; 4 ] ; % a v e r t i c a l v e c t o r h = [ ] ; % a h o r i z o n t a l v e c t o r M = [ [ 1 ; 2 ] [ 3 ; 4 ] ] ; % a 2 by 2 Matrix s = 1 : 4 ; % the same as h i n l i n e 2 Calculate with vectors: u = 5 v % m u l t i p l y v e c t o r s by a s c a l a r w = u v % s u b t r a c t two v e c t o r s x = u. v % e l e m e n t w i s e m u l t i p l i c a t i o n u v, u v % i n n e r & o u t e r p roduct y = [ u v ] % c o n c a t e n a t i o n o f v e c t o r s Example: polar to Euclidean coordinates x = s i n ( t h e t a ). cos ( rho ) y = s i n ( t h e t a ). s i n ( rho )
13 Vectors Indexing by numbers: v ( 1 ) % the f i r s t element v ( 2 : end ) % the l a s t t h r e e e l e m e n t s M( 2, 1 ) % the bottom l e f t element M( :, 2 ) % the second column Indexing by conditions: v ( v > 1 ) % a l l e l e m e n t s t h a t a r e l a r g e r then 1 v ( u == v ) % a l l e l e m e n t s where u and v c o i n c i d e Changing elements v ( 1 : 3 ) = 1 % s e t the f i r s t t h r e e e l e m e n t s to 1 v ( v == 0 ) = [ ] % remove a l l z e r o s
14 MATLAB programming style - Orientation tensor 1 For a pole figure displaying one crystallographic direction h of number (N) individual orientation pixels or grains the orientation tensor is traditionally defined in the geological literature as: M ij = 1 N N i xi 2 N i x i y N i i x i z i N i x i y N i i yi 2 N i y i z i N i x i z N i i y i z N i i zi 2 where r = (x i y i z i ) are direction cosines of a crystallographic direction h in specimen coordinates r with Cartesian reference axes X, Y and Z.
15 MATLAB programming style - Orientation tensor 2 The same Orientation tensor can also be formulated as a symmetric matrix and the product of the multiplication of a vector containing all r by its transpose: M ij = M 11 M 12 M 13 M 12 M 22 M 23 M 13 M 23 M 33 = 1 N N r T i r i i where r = (x i y i z i ) are direction cosines of a crystallographic direction h in specimen coordinates r with Cartesian reference axes X, Y and Z.
16 MATLAB programming style - Orientation tensor 3 The two methods of calculation are illustrated by small MATLAB script. The method using the formation of covariance matrix from a column and row vector of r is 300 times faster. Efficient programming requires adapting to the MATLAB style.
17 MATLAB programming style - Orientation tensor 4 An example of c-axes pole figure of ice 1h at 1171m Talos Dome ice core, Antarctica. Eigen-vectors (red squares) : E1, E2, E3 Eigen-values (magnitudes) : EM1 > EM2 > EM3 G50 fabric analyser data courtesy of Maurine Montagnat LGGE UMR5183, CNRS / UJF Grenoble, France
18 Single Orientation Measurements in MTEX Parameters of a single EBSD phase in MTEX arbitrary as list: crystal symmetry specimen symmetry list of individual orientations phase example olivine spatial coordinates x y Band Contrast diffraction signal MAD Mean Angular Deviation... Olivine in torsion: EBSD map data courtesy of Sylvie Demouchy,Montpellier Colorcode for crystallographic directions parallel to shear direction
19 Importing EBSD Data with The MTEX Import Wizard Two ways to import EBSD data 1 From command window link Import EBSD data 2 Type in the command window import_wizard( EBSD ) Formats supported by MTEX : *.ang TSL-EDAX *.ctf HKL-Oxford *.csv Oxford *.sor LaboTEX *.txt (Generic ASCII files generated by user)
20 Importing EBSD Data - Using a script This script was generated by the EBSD import wizard: % C r y s t a l s y m m e t r i e s d e f i n e a, b, c, alpha, beta, gamma % and c r y s t a l d i r e c t i o n s used f o r E u l e r a n g l e C a r t e s i a n r e f e r e n c e frame CS = {... n o t I n d e x e d,... symmetry ( mmm, [ ], m i n e r a l, O l i v i n e ),... symmetry ( mmm, [ ], m i n e r a l, E n s t a t i t e ),... symmetry ( 2/m, [ ], [ , , ] degree,... X a, Y b, Z c, m i n e r a l, D i o p s i d e ),... symmetry ( m 3m, m i n e r a l, Magnetite ),... symmetry ( m 3m, m i n e r a l, Chromite ) } ; % Specimen symmetry SS = symmetry ( 1 ) ; % path to f i l e s pname = / MatLab Programs / % which f i l e s to be i m p o r t e d fname = {[ pname 920D 13. c t f ], } % c r e a t e an EBSD v a r i a b l e c o n t a i n i n g t h e data ebsd = loadebsd ( fname, CS, SS, i n t e r f a c e, c t f )
21 Visualize EBSD Data in MTEX : orientation space Scatter plots in orientation space (Rodrigues space, axis-angle space or Euler angle space) % Not v e r y u s e f u l as a x e s a r e not l a b e l l e d s c a t t e r ( ebsd ( O l i v i n e ), AXISANGLE ) s c a t t e r ( ebsd ( O l i v i n e ), RODRIGUES ) % More u s e f u l w i t h l a b e l l e d s e c t i o n s p l o t o d f ( ebsd ( O l i v i n e ), sigma, SECTIONS, 1 8, p o i n t s,320744)
22 Visualize EBSD Data in MTEX : specimen space Scatter plots of pole figures: crystal direction h in specimen space r % O l i v i n e p o l e f i g u r e s % p l o t t i n g c o n v e n t i o n p l o t x 2 e a s t % X=East Y=North Z=C e n t e r upper h e m i s p h e r e p l o t x 2 e a s t % s p e c i f y c r y s t a l d i r e c t i o n s h w i t h h k l f o r p o l e and uvw d i r e c t i o n s h = [ M i l l e r ( 1, 0, 0, uvw ), M i l l e r ( 0, 1, 0, uvw ), M i l l e r ( 0, 0, 1, uvw ) ] f i g u r e ( p o s i t i o n, [ ] ) p l o t p d f ( ebsd ( O l i v i n e ), h, a n t i p o d a l, p o i n t s,320744) % s p e c i f y specimen d i r e c t i o n s r % you can used predefined xvector, yvector or zvector % o r d e f i n e a r b i t a r y d i r e c t i o n s w i t h v e c t o r 3 d ( 1. 0, 0. 0, 0. 0 ) r = [ v e c t o r 3 d ( 1. 0, 0. 0, 0. 0 ), v e c t o r 3 d ( 0. 0, 1. 0, 0. 0 ), v e c t o r 3 d ( 0. 0, 0. 0, 1. 0 ) ] a n n o t a t e ( r, a l l, l a b e l,{ X, Y, Z }, BackgroundColor, w ) Option : Antipodal (non-polar) Option : Complete (polar)
23 Visualize EBSD Data in MTEX : crystal space Scatter plots of inverse pole figures: specimen direction r in crystal space h % O l i v i n e i n v e r s e p o l e f i g u r e s : specimen d i r e c t i o n r i n c r y s t a l s p a c e h % s p e c i f y specimen d i r e c t i o n s r r = [ v e c t o r 3 d ( 1. 0, 0. 0, 0. 0 ), v e c t o r 3 d ( 0. 0, 1. 0, 0. 0 ), v e c t o r 3 d ( 0. 0, 0. 0, 1. 0 ) ] f i g u r e ( p o s i t i o n, [ ] ) p l o t i p d f ( ebsd ( O l i v i n e ), r, a n t i p o d a l, p o i n t s,320744) % s p e c i f y c r y s t a l d i r e c t i o n s h h = [ M i l l e r ( 1, 0, 0, uvw ), M i l l e r ( 0, 1, 0, uvw ), M i l l e r ( 0, 0, 1, uvw ) ] a n n o t a t e ( h, a l l, l a b e l e d, BackgroundColor, w )
24 Visualize EBSD Data in MTEX : Spatial maps Spatial plots of EBSD data Example hkl color coded orientation map % i p f o r i e n t a t i o n map p l o t ( ebsd, c o l o r c o d i n g, i p d f ) % phase d i s t r i b u t i o n map p l o t ( ebsd, p r o p e r t y, phase ) % d i f f r a c t e d i n t e n s i t y map p l o t ( ebsd, p r o p e r t y, bc ) % mean a n g u l a r d e v i a t i o n e r r o r map p l o t ( ebsd, p r o p e r t y, mad ) Color Codings: bunge, angle, sigma, ihs, ipdf, rodriguesquat, rodriguesinverse, euler, bunge2, hkl (Olivine torsion data courtesy of Sylvie Demouchy, GM, UM2, Montpellier)
25 Grain Analysis with MTEX Grain reconstruction in MTEX is done via the command calcgrains. As an optional argument the desired threshold angle for misorientation defining a grains boundary can be specified. A grain is defined as a region, in which the misorientation of at least one neighbouring measurement site is smaller than a given threshold angle. The command grainsize returns the number of indexed pixels within a grain. % F i l l s a l l non i n d e x e d a r e a s. Angle i s t h e same f o r a l l p h a s e s. g r a i n s = c a l c G r a i n s ( ebsd, a n g l e,10 degree ) ; % Angle i s 0,10 & 15 f o r p h a s e s 1, 2, 3. F i l l s a l l non i n d e x e d a r e a s. g r a i n s = c a l c G r a i n s ( ebsd, a n g l e, [ ] degree ) ; % Conserves non indexed areas. Angle i s the same f o r a l l phases. g r a i n s = c a l c G r a i n s ( ebsd, a n g l e,10 degree, keepnotindexed ) ; EBSD Map Grain Map : calcgrains(ebsd( Olivine ), angle,10*degree)
26 Grain Properties EBSD properties p l ot ( g r a i n s ) % p l o t a l l g r a i n s p l ot ( g r a i n s ( O l i v i n e ) % p l o t o l i v i n e g r a i n s p l ot ( g r a i n s, p r o p e r t y, phase ) % p l o t the phase o f eac plotboundary ( g r a i n s ) % p l o t g r a i n b o u n d a r i e s Plot grain boundaries: p l o t b o u n d a r y ( g r a i n s,/+<o p t i o n s>+/) p l o t b o u n d a r y ( g r a i n s, p r o p e r t y, a n g l e ) p l o t b o u n d a r y ( g r a i n s, p r o p e r t y, r o t ) Access to properties of grains phase = get ( g r a i n s, phase ) g r a i n s q u a r t z = g r a i n s ( Quartz ) o = get ( g r a i n s q u a r t z, o r i e n t a t i o n )
27 Grain Properties EBSD properties p l ot ( g r a i n s,/+<o p t i o n s>+/) /+p l o t ( g r a i n s, p r o p e r t y, phase )+/ Plot grain boundaries: p l o t b o u n d a r y ( g r a i n s,/+<o p t i o n s>+/) p l o t b o u n d a r y ( g r a i n s, p r o p e r t y, a n g l e ) p l o t b o u n d a r y ( g r a i n s, p r o p e r t y, r o t ) Access to properties of grains phase = get ( g r a i n s, phase ) g r a i n s q u a r t z = g r a i n s ( Quartz ) o = get ( g r a i n s q u a r t z, o r i e n t a t i o n )
28 Grain Properties: Geometry Basic functions of grain geometry % commands that i n c l u d e non index points ( holes ) p e r i m e t e r ( g r a i n s ) ; % p e r i m e t e r o f a g r a i n p o l y g o n a r e a ( g r a i n s ) ; % a r e a o f a g r a i n p o l y g o n c e n t r o i d ( g r a i n s ) ; % b a r y c e n t e r o f t h e g r a i n p o l y g o n % e q u i v a l e n t r a d i u s and p e r i m e t e r e q u i v a l e n t r a d i u s ( g r a i n s ) ; % d e f i n e d as s q r t ( a r e a / p i ) e q u i v a l e n t p e r i m e t e r ( g r a i n s ) ; % d e f i n e d as 2 p i e q u i v a l e n t r a d i u s % g r a i n geometry d i a m e t e r ( g r a i n s ) ; % l o n g e s t d i s t a n c e between two v e r t i c e s o f g r a i n a s p e c t r a t i o ( g r a i n s ) ; % r a t i o o f two p r i n c i p a l components o f a g r a i n s h a p e f a c t o r ( g r a i n s ) ; % d e f i n e d as p e r i m e t e r / e q u i v a l e n t p e r i m e t e r g r a i n S i z e ( g r a i n s ) ; % t h e number o f i n d e x e d measurements p e r g r a i n Grain-size distribution A = a r e a ( g r a i n s ) ; % h i s t o g r a m bar ( h i s t (A, exp ( 1. 5 : 6. 5 ) ) ) Other properties: h a s h o l e s, h a s s u b f r a c t i o n
29 Working with EBSD Data Directly Access EBSD properties o r i e n t a t i o n s = get ( ebsd, o r i e n t a t i o n ) x = get ( ebsd, x ) Rotate EBSD data ebsd = rotate ( ebsd, r o t a t i o n (... a x i s, xvector, a n g l e,90 degree ) ) Omit one pixel grains % g r a i n s c o n t a i n i n g more than one i n d e x e d p i x e l g r a i n s = g r a i n s ( g r a i n s i z e ( g r a i n s ) > 1)
30 EBSD to ODF Reconstruction in MTEX MTEX uses kernel density estimation to compute an ODF from EBSD data. The sensitive parameter of this method is the kernel function. Syntax: odf = calcodf ( ebsd,<o p t i o n s>) Options: k e r n e l % the k e r n e l to be used h a l f w i d t h % h a l f width o f the k e r n e l r e s o l u t i o n % r e s o l u t i o n o f a p p r o x i m a t i o n g r i d F o u r i e r % ODF by i t s F o u r i e r c o e f f i c i e n t s e x a c t % e x a c t g r i d
Calculating anisotropic physical properties from texture data using the MTEX open source package
Calculating anisotropic physical properties from texture data using the MTEX open source package David Mainprice 1, Ralf Hielscher 2, Helmut Schaeben 3 February 10, 2011 1 Geosciences Montpellier UMR CNRS
More informationGeological Society, London, Special Publications
Geological Society, London, Special Publications Calculating anisotropic physical properties from texture data using the MTEX open-source package David Mainprice, Ralf Hielscher and Helmut Schaeben Geological
More informationL3: Review of linear algebra and MATLAB
L3: Review of linear algebra and MATLAB Vector and matrix notation Vectors Matrices Vector spaces Linear transformations Eigenvalues and eigenvectors MATLAB primer CSCE 666 Pattern Analysis Ricardo Gutierrez-Osuna
More informationA Generalised Spherical Harmonic Deconvolution to Obtain Cubic Texture from Ultrasonic Wave Speed
A Generalised Spherical Harmonic Deconvolution to Obtain Cubic Texture from Ultrasonic Wave Speed Bo Lan 1, Michael J.S. Lowe 2, and Fionn P.E. Dunne 2 1 Department of Engineering Science, University of
More informationTexture, Microstructure & Anisotropy A.D. (Tony) Rollett
1 Carnegie Mellon MRSEC 27-750 Texture, Microstructure & Anisotropy A.D. (Tony) Rollett Last revised: 5 th Sep. 2011 2 Show how to convert from a description of a crystal orientation based on Miller indices
More informationTensors and Anisotropic physical properties of rocks Part I : Tensor representation of the physical properties of single crystals
Tensors and Anisotropic physical properties of rocks Part I : Tensor representation of the physical properties of single crystals David Mainprice Lab. Tectonophysique, Université Montpellier II France
More informationExercise: concepts from chapter 8
Reading: Fundamentals of Structural Geology, Ch 8 1) The following exercises explore elementary concepts associated with a linear elastic material that is isotropic and homogeneous with respect to elastic
More informationMath 307 Learning Goals. March 23, 2010
Math 307 Learning Goals March 23, 2010 Course Description The course presents core concepts of linear algebra by focusing on applications in Science and Engineering. Examples of applications from recent
More informationChapter 5. Vibration Analysis. Workbench - Mechanical Introduction ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.
Workbench - Mechanical Introduction 12.0 Chapter 5 Vibration Analysis 5-1 Chapter Overview In this chapter, performing free vibration analyses in Simulation will be covered. In Simulation, performing a
More informationImage Compression Using the Haar Wavelet Transform
College of the Redwoods http://online.redwoods.cc.ca.us/instruct/darnold/laproj/fall2002/ames/ 1/33 Image Compression Using the Haar Wavelet Transform Greg Ames College of the Redwoods Math 45 Linear Algebra
More informationECE521 W17 Tutorial 1. Renjie Liao & Min Bai
ECE521 W17 Tutorial 1 Renjie Liao & Min Bai Schedule Linear Algebra Review Matrices, vectors Basic operations Introduction to TensorFlow NumPy Computational Graphs Basic Examples Linear Algebra Review
More informationElements of Continuum Elasticity. David M. Parks Mechanics and Materials II February 25, 2004
Elements of Continuum Elasticity David M. Parks Mechanics and Materials II 2.002 February 25, 2004 Solid Mechanics in 3 Dimensions: stress/equilibrium, strain/displacement, and intro to linear elastic
More informationQ. WANG, Q-K. XIA, S. Y. O REILLY, W. L. GRIFFIN, E. E. BEYER AND H. K. BRUECKNER
Pressure- and stress-induced fabric transition in olivine from peridotites in the Western Gneiss Region (Norway): implications for mantle seismic anisotropy Q. WANG, Q-K. XIA, S. Y. O REILLY, W. L. GRIFFIN,
More informationINTERSECTIONS OF PLANES
GG303 Lab 4 9/17/14 1 INTERSECTIONS OF PLANES Lab 4 (125 points total) Read each exercise completely before you start it so that you understand the problem-solving approach you are asked to execute This
More informationCopyright(c)JCPDS-International Centre for Diffraction Data 2000,Advances in X-ray Analysis,Vol ISSN
Copyright(c)JCPDS-International Centre for Diffraction Data 2000,Advances in X-ray Analysis,Vol.43 129 MATHEMATICAL OF DIFFRACTION PROPERTIES POLE FIGURES ABSTRACT Helmut Schaeben Mathematics and Computer
More informationINTERSECTIONS OF PLANES
GG303 Lab 4 10/9/13 1 INTERSECTIONS OF PLANES Lab 4 Read each exercise completely before you start it so that you understand the problem-solving approach you are asked to execute This will help keep the
More informationApplied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB Contents Getting Started Creating Arrays Mathematical Operations with Arrays Using Script Files and Managing Data Two-Dimensional Plots Programming in
More informationHigh-resolution finite volume methods for hyperbolic PDEs on manifolds
High-resolution finite volume methods for hyperbolic PDEs on manifolds Randall J. LeVeque Department of Applied Mathematics University of Washington Supported in part by NSF, DOE Overview High-resolution
More informationInferential statistics of electron backscatter diffraction data from within individual crystalline grains
Journal of Applied Crystallography ISSN 2-8898 Editor: Anke R. Kaysser-Pyzalla Inferential statistics of electron backscatter diffraction data from within individual crystalline grains Florian Bachmann,
More informationConstitutive Equations
Constitutive quations David Roylance Department of Materials Science and ngineering Massachusetts Institute of Technology Cambridge, MA 0239 October 4, 2000 Introduction The modules on kinematics (Module
More informationFrequently Asked Questions
Frequently Asked Questions Why do we have to make the assumption that plane sections plane? How about bars with non-axis symmetric cross section? The formulae derived look very similar to beam and axial
More informationChapter 1. Gaining Knowledge with Design of Experiments
Chapter 1 Gaining Knowledge with Design of Experiments 1.1 Introduction 2 1.2 The Process of Knowledge Acquisition 2 1.2.1 Choosing the Experimental Method 5 1.2.2 Analyzing the Results 5 1.2.3 Progressively
More information(Mathematical Operations with Arrays) Applied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB (Mathematical Operations with Arrays) Contents Getting Started Matrices Creating Arrays Linear equations Mathematical Operations with Arrays Using Script
More informationSpherical Wavelets with an Application in Preferred Crystallographic Orientation
Spherical Wavelets with an Application in Preferred Crystallographic Orientation Helmut Schaeben ( ), Daniel Potts ( ), and Jürgen Prestin ( ) ( ) Mathematics and Computer Sciences in Geology Freiberg
More informationLecture 4 Implementing material models: using usermat.f. Implementing User-Programmable Features (UPFs) in ANSYS ANSYS, Inc.
Lecture 4 Implementing material models: using usermat.f Implementing User-Programmable Features (UPFs) in ANSYS 1 Lecture overview What is usermat.f used for? Stress, strain and material Jacobian matrix
More informationFinite Element Simulations of Microbeam Bending Experiments
Finite Element Simulations of Microbeam Bending Experiments Master s thesis in Applied Mechanics JOHN WIKSTRÖM Department of Applied Mechanics CHALMERS UNIVERSITY OF TECHNOLOGY Göteborg, Sweden 2017 MASTER
More informationSuperMix2 features not available in HLM 7 Contents
SuperMix2 features not available in HLM 7 Contents Spreadsheet display of.ss3 files... 2 Continuous Outcome Variables: Additional Distributions... 3 Additional Estimation Methods... 5 Count variables including
More informationQuadratic and cubic monocrystalline and polycrystalline materials: their stability and mechanical properties
Journal of Physics: Conference Series Quadratic and cubic monocrystalline and polycrystalline materials: their stability and mechanical properties To cite this article: C Jasiukiewicz et al 010 J. Phys.:
More informationPhysics 584 Computational Methods
Physics 584 Computational Methods Introduction to Matlab and Numerical Solutions to Ordinary Differential Equations Ryan Ogliore April 18 th, 2016 Lecture Outline Introduction to Matlab Numerical Solutions
More informationRigid Geometric Transformations
Rigid Geometric Transformations Carlo Tomasi This note is a quick refresher of the geometry of rigid transformations in three-dimensional space, expressed in Cartesian coordinates. 1 Cartesian Coordinates
More informationIDENTIFICATION OF THE ELASTIC PROPERTIES ON COMPOSITE MATERIALS AS A FUNCTION OF TEMPERATURE
IDENTIFICATION OF THE ELASTIC PROPERTIES ON COMPOSITE MATERIALS AS A FUNCTION OF TEMPERATURE Hugo Sol, hugos@vub.ac.be Massimo Bottiglieri, Massimo.Bottiglieri@vub.ac.be Department Mechanics of Materials
More informationImage Compression. 1. Introduction. Greg Ames Dec 07, 2002
Image Compression Greg Ames Dec 07, 2002 Abstract Digital images require large amounts of memory to store and, when retrieved from the internet, can take a considerable amount of time to download. The
More informationMATERIAL MECHANICS, SE2126 COMPUTER LAB 4 MICRO MECHANICS. E E v E E E E E v E E + + = m f f. f f
MATRIAL MCHANICS, S226 COMPUTR LAB 4 MICRO MCHANICS 2 2 2 f m f f m T m f m f f m v v + + = + PART A SPHRICAL PARTICL INCLUSION Consider a solid granular material, a so called particle composite, shown
More informationCharacterizations of Aluminum Alloy Sheet Materials Numisheet 2005
Characterizations of Aluminum Alloy Sheet Materials Numisheet 25 John C. Brem, Frederic Barlat, Robert E. Dick, and Jeong-Whan Yoon Alcoa Technical Center, PA, 1569-1, USA Abstract. This report reproduces
More informationExperiment 1: Linear Regression
Experiment 1: Linear Regression August 27, 2018 1 Description This first exercise will give you practice with linear regression. These exercises have been extensively tested with Matlab, but they should
More informationModule 10: Free Vibration of an Undampened 1D Cantilever Beam
Module 10: Free Vibration of an Undampened 1D Cantilever Beam Table of Contents Page Number Problem Description Theory Geometry 4 Preprocessor 6 Element Type 6 Real Constants and Material Properties 7
More informationCAAM 335 Matrix Analysis Planar Trusses
CAAM 5 Matrix Analysis Planar Trusses September 1, 010 1 The Equations for the Truss We consider trusses with m bars and n nodes. Each node can be displaced in horizontal and vertical direction. If the
More informationIntroduction to Finite Element Analysis Using Pro/MECHANICA Wildfire 5.0
Introduction to Finite Element Analysis Using Pro/MECHANICA Wildfire 5.0 Randy H. Shih Oregon Institute of Technology SDC PUBLICATIONS Schroff Development Corporation www.schroff.com Better Textbooks.
More informationFundamental principles
Dynamics and control of mechanical systems Date Day 1 (03/05) - 05/05 Day (07/05) Day 3 (09/05) Day 4 (11/05) Day 5 (14/05) Day 6 (16/05) Content Review of the basics of mechanics. Kinematics of rigid
More informationPHYSICAL PROPERTIES OF CRYSTALS
PHYSICAL PROPERTIES OF CRYSTALS THEIR REPRESENTATION TENSORS AND MATRICES BY By J. F. NYE, F.R.S. CLARENDON PRESS OXFORD NOTATION INTRODUCTION xiii xv PART 1. GENERAL PRINCIPLES I. THE GROUNDWORK OF CRYSTAL
More information3D Stress Tensors. 3D Stress Tensors, Eigenvalues and Rotations
3D Stress Tensors 3D Stress Tensors, Eigenvalues and Rotations Recall that we can think of an n x n matrix Mij as a transformation matrix that transforms a vector xi to give a new vector yj (first index
More informationIntroduction. Finite and Spectral Element Methods Using MATLAB. Second Edition. C. Pozrikidis. University of Massachusetts Amherst, USA
Introduction to Finite and Spectral Element Methods Using MATLAB Second Edition C. Pozrikidis University of Massachusetts Amherst, USA (g) CRC Press Taylor & Francis Group Boca Raton London New York CRC
More informationCrystallographic Symmetry. Jeremy Karl Cockcroft
Crystallographic Symmetry Jeremy Karl Cockcroft Why bother? To describe crystal structures Simplifies the description, e.g. NaCl structure Requires coordinates for just 2 atoms + space group symmetry!
More informationPiezo materials. Actuators Sensors Generators Transducers. Piezoelectric materials may be used to produce e.g.: Piezo materials Ver1404
Noliac Group develops and manufactures piezoelectric materials based on modified lead zirconate titanate (PZT) of high quality and tailored for custom specifications. Piezoelectric materials may be used
More informationSingular Value Decomposition and Digital Image Compression
Singular Value Decomposition and Digital Image Compression Chris Bingham December 1, 016 Page 1 of Abstract The purpose of this document is to be a very basic introduction to the singular value decomposition
More informationLinear Inverse Problems. A MATLAB Tutorial Presented by Johnny Samuels
Linear Inverse Problems A MATLAB Tutorial Presented by Johnny Samuels What do we want to do? We want to develop a method to determine the best fit to a set of data: e.g. The Plan Review pertinent linear
More informationUNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES
UNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES A Thesis by WOORAM KIM Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the
More informationMultimedia & Computer Visualization. Exercise #5. JPEG compression
dr inż. Jacek Jarnicki, dr inż. Marek Woda Institute of Computer Engineering, Control and Robotics Wroclaw University of Technology {jacek.jarnicki, marek.woda}@pwr.wroc.pl Exercise #5 JPEG compression
More informationDAMOCO: MATLAB toolbox for multivariate data analysis, based on coupled oscillators approach Version 1.0
DAMOCO: MATLAB toolbox for multivariate data analysis, based on coupled oscillators approach Version 1.0 Björn Kralemann 1, Michael Rosenblum 2, Arkady Pikovsky 2 1 Institut für Pädagogik, Christian-Albrechts-Universität
More informationLindgren CRYSTAL SYMMETRY AND ELASTIC CONSTANTS MICHAEL WANDZILAK. S.B., Massachusetts Institute of Technology (196'7)
CRYSTAL SYMMETRY AND ELASTIC CONSTANTS by MICHAEL WANDZILAK S.B., Massachusetts Institute of Technology (196'7) Submitted in partial fulfillment of the requirements for the degree of Master of Science
More informationLinear Algebra Review. Fei-Fei Li
Linear Algebra Review Fei-Fei Li 1 / 37 Vectors Vectors and matrices are just collections of ordered numbers that represent something: movements in space, scaling factors, pixel brightnesses, etc. A vector
More informationConstitutive Relations
Constitutive Relations Dr. Andri Andriyana Centre de Mise en Forme des Matériaux, CEMEF UMR CNRS 7635 École des Mines de Paris, 06904 Sophia Antipolis, France Spring, 2008 Outline Outline 1 Review of field
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Lecture - 06
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Lecture - 06 In the last lecture, we have seen a boundary value problem, using the formal
More informationProbabilistic Latent Semantic Analysis
Probabilistic Latent Semantic Analysis Seungjin Choi Department of Computer Science and Engineering Pohang University of Science and Technology 77 Cheongam-ro, Nam-gu, Pohang 37673, Korea seungjin@postech.ac.kr
More informationMath 671: Tensor Train decomposition methods
Math 671: Eduardo Corona 1 1 University of Michigan at Ann Arbor December 8, 2016 Table of Contents 1 Preliminaries and goal 2 Unfolding matrices for tensorized arrays The Tensor Train decomposition 3
More informationFeb 21 and 25: Local weighted least squares: Quadratic loess smoother
Feb 1 and 5: Local weighted least squares: Quadratic loess smoother An example of weighted least squares fitting of data to a simple model for the purposes of simultaneous smoothing and interpolation is
More informationSpectral Graph Theory
Spectral Graph Theory Aaron Mishtal April 27, 2016 1 / 36 Outline Overview Linear Algebra Primer History Theory Applications Open Problems Homework Problems References 2 / 36 Outline Overview Linear Algebra
More informationCOSY INFINITY Version 7
COSY INFINITY Version 7 Kyoko Makino and Martin Berz Department of Physics and Astronomy and National Superconducting Cyclotron Laboratory Michigan State University, East Lansing, MI 48824 Abstract An
More informationLab 2 Worksheet. Problems. Problem 1: Geometry and Linear Equations
Lab 2 Worksheet Problems Problem : Geometry and Linear Equations Linear algebra is, first and foremost, the study of systems of linear equations. You are going to encounter linear systems frequently in
More informationComposite models 30 and 131: Ply types 0 and 8 calibration
Model calibration Composite Bi-Phase models 30 and 3 for elastic, damage and failure PAM-CRASH material model 30 is for solid and 3 for multi-layered shell elements. Within these models different ply types
More informationCONVERSION OF COORDINATES BETWEEN FRAMES
ECS 178 Course Notes CONVERSION OF COORDINATES BETWEEN FRAMES Kenneth I. Joy Institute for Data Analysis and Visualization Department of Computer Science University of California, Davis Overview Frames
More informationTechnische Universität Kaiserslautern WS 2016/17 Fachbereich Mathematik Prof. Dr. J. Franke 19. January 2017
Technische Universität Kaiserslautern WS 2016/17 Fachbereich Mathematik Prof. Dr. J. Franke 19. January 2017 Exercises for Computational Modelling with Statistics IV This exercise is about the spectral
More information(Linear equations) Applied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB (Linear equations) Contents Getting Started Creating Arrays Mathematical Operations with Arrays Using Script Files and Managing Data Two-Dimensional Plots
More informationMOOSE Workshop MOOSE-PF/MARMOT Training Mechanics Part Nashville, February 14th Daniel Schwen
www.inl.gov MOOSE Workshop MOOSE-PF/MARMOT Training Mechanics Part Nashville, February 14th Daniel Schwen www.inl.gov III. Mechanics Mechanics Outline 1. Introduction 2. Tensor-based System 1. Compute
More informationLecture contents. Stress and strain Deformation potential. NNSE 618 Lecture #23
1 Lecture contents Stress and strain Deformation potential Few concepts from linear elasticity theory : Stress and Strain 6 independent components 2 Stress = force/area ( 3x3 symmetric tensor! ) ij ji
More informationRigid body Rigid body approach
Rigid body Rigid body approach Model molecule is a fragment, which is placed arbitrarily and does not contribute to structure factors Model molecule is transformed to Actual positions by translation vector
More informationComputational Foundations of Cognitive Science
Computational Foundations of Cognitive Science Lecture 14: Inverses and Eigenvectors in Matlab; Plotting and Graphics Frank Keller School of Informatics University of Edinburgh keller@inf.ed.ac.uk February
More information! EN! EU! NE! EE.! ij! NN! NU! UE! UN! UU
A-1 Appendix A. Equations for Translating Between Stress Matrices, Fault Parameters, and P-T Axes Coordinate Systems and Rotations We use the same right-handed coordinate system as Andy Michael s program,
More informationPhysical properties in MTEX, current status and future developments
Physical properties in MTEX, current status and future developments Chemnitz MTEX Workshop 9-13 February 2015 Technische Universität Chemnitz, Germnay David Mainprice (Géosciences Montpellier, Université
More informationLinear Algebra Review. Fei-Fei Li
Linear Algebra Review Fei-Fei Li 1 / 51 Vectors Vectors and matrices are just collections of ordered numbers that represent something: movements in space, scaling factors, pixel brightnesses, etc. A vector
More informationarxiv: v1 [physics.ins-det] 4 Jun 2018
Noname manuscript No. (will be inserted by the editor) Simulation of force-insensitive optical cavities in cubic spacers Eugen Wiens Stephan Schiller arxiv:186.2176v1 [physics.ins-det] 4 Jun 218 Received:
More informationSelecting protein fuzzy contact maps through information and structure measures
Selecting protein fuzzy contact maps through information and structure measures Carlos Bousoño-Calzón Signal Processing and Communication Dpt. Univ. Carlos III de Madrid Avda. de la Universidad, 30 28911
More informationDeriving Principal Component Analysis (PCA)
-0 Mathematical Foundations for Machine Learning Machine Learning Department School of Computer Science Carnegie Mellon University Deriving Principal Component Analysis (PCA) Matt Gormley Lecture 11 Oct.
More informationA Brief Introduction To. GRTensor. On MAPLE Platform. A write-up for the presentation delivered on the same topic as a part of the course PHYS 601
A Brief Introduction To GRTensor On MAPLE Platform A write-up for the presentation delivered on the same topic as a part of the course PHYS 601 March 2012 BY: ARSHDEEP SINGH BHATIA arshdeepsb@gmail.com
More informationCHAPTER 6 MECHANICAL PROPERTIES OF METALS PROBLEM SOLUTIONS
CHAPTER 6 MECHANICAL PROPERTIES OF METALS PROBLEM SOLUTIONS Concepts of Stress and Strain 6.1 Using mechanics of materials principles (i.e., equations of mechanical equilibrium applied to a free-body diagram),
More informationCraFT user s guide. Hervé Moulinec
CraFT user s guide Hervé Moulinec December 12, 2014 1 Introduction In this chapter, one will describe how to run CraFT. That means what sort of input data CraFT needs, how to specify it to CraFT, what
More informationTextures in experimentally deformed olivine aggregates: the effects of added water and melt.
Textures in experimentally deformed olivine aggregates: the effects of added water and melt. F. Heidelbach 1, a, B. Holtzman 2, b, S. Hier-Majumder 2, c and D. Kohlstedt 2, d 1 Bayerisches Geoinstitut,
More informationSYMBOLIC AND NUMERICAL COMPUTING FOR CHEMICAL KINETIC REACTION SCHEMES
SYMBOLIC AND NUMERICAL COMPUTING FOR CHEMICAL KINETIC REACTION SCHEMES by Mark H. Holmes Yuklun Au J. W. Stayman Department of Mathematical Sciences Rensselaer Polytechnic Institute, Troy, NY, 12180 Abstract
More informationThe tensors useful for geophysics and estimation of anisotropic polycrystalline physical properties
The tensors useful for geophysics and estimation of anisotropic polycrystalline physical properties David Mainprice Laboratoire de Tectonophysique, ISTEEM, CNRS UMR 5558, Université Montpellier II, 34095
More informationGSI Geometric Science of Information, Paris, August Dimensionality reduction for classification of stochastic fibre radiographs
GSI2013 - Geometric Science of Information, Paris, 28-30 August 2013 Dimensionality reduction for classification of stochastic fibre radiographs C.T.J. Dodson 1 and W.W. Sampson 2 School of Mathematics
More informationTransducer Design and Modeling 42 nd Annual UIA Symposium Orlando Florida Jay Sheehan JFS Engineering. 4/23/2013 JFS Engineering
42 nd Annual UIA Symposium Orlando Florida 2013 Jay Sheehan JFS Engineering Introduction ANSYS Workbench Introduction The project format Setting up different analysis Static, Modal and Harmonic Connection
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Module - 01 Lecture - 13
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras (Refer Slide Time: 00:25) Module - 01 Lecture - 13 In the last class, we have seen how
More informationThree-Dimensional Electron Microscopy of Macromolecular Assemblies
Three-Dimensional Electron Microscopy of Macromolecular Assemblies Joachim Frank Wadsworth Center for Laboratories and Research State of New York Department of Health The Governor Nelson A. Rockefeller
More informationLearning goals: students learn to use the SVD to find good approximations to matrices and to compute the pseudoinverse.
Application of the SVD: Compression and Pseudoinverse Learning goals: students learn to use the SVD to find good approximations to matrices and to compute the pseudoinverse. Low rank approximation One
More informationMath 671: Tensor Train decomposition methods II
Math 671: Tensor Train decomposition methods II Eduardo Corona 1 1 University of Michigan at Ann Arbor December 13, 2016 Table of Contents 1 What we ve talked about so far: 2 The Tensor Train decomposition
More information3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is
More informationEE Experiment 11 The Laplace Transform and Control System Characteristics
EE216:11 1 EE 216 - Experiment 11 The Laplace Transform and Control System Characteristics Objectives: To illustrate computer usage in determining inverse Laplace transforms. Also to determine useful signal
More informationElaboration on the Hexagonal Grid and Spiral Method for Data Collection Via Pole Figures
Student Projects Elaboration on the Hexagonal Grid and Spiral Method for Data Collection Via Pole Figures Anthony C. Rizzie Anthony Rizzie is a junior majoring in Mathematics (Opt. 1) with a double minor
More informationIntroduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams.
Outline of Continuous Systems. Introduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams. Vibrations of Flexible Strings. Torsional Vibration of Rods. Bernoulli-Euler Beams.
More informationLinear Neural Networks
Chapter 10 Linear Neural Networks In this chapter, we introduce the concept of the linear neural network. 10.1 Introduction and Notation 1. The linear neural cell, or node has the schematic form as shown
More informationTEMScripts Real-time Crystallography Manual. TEMScripts LLC. Last updated: 11/4/2016
TEMScripts Real-time Crystallography Manual TEMScripts LLC. Last updated: 11/4/2016 Close Digital Micrograph Installation Copy following files to \\Gatan\DigitalMicrograph\PlugIns (normally under C:\Program
More informationExploring Piezoelectric Properties of Wood and Related Issues in Mathematical Description. Igor Dobovšek
Exploring Piezoelectric Properties of Wood and Related Issues in Mathematical Description Igor Dobovšek University of Ljubljana Faculty of Mathematics and Physics Institute of Mathematics Physics and Mechanics
More informationECE 661: Homework 10 Fall 2014
ECE 661: Homework 10 Fall 2014 This homework consists of the following two parts: (1) Face recognition with PCA and LDA for dimensionality reduction and the nearest-neighborhood rule for classification;
More informationImage Data Compression
Image Data Compression Image data compression is important for - image archiving e.g. satellite data - image transmission e.g. web data - multimedia applications e.g. desk-top editing Image data compression
More informationIntroduction to Matlab
History of Matlab Starting Matlab Matrix operation Introduction to Matlab Useful commands in linear algebra Scripts-M file Use Matlab to explore the notion of span and the geometry of eigenvalues and eigenvectors.
More informationJanuary 18, 2008 Steve Gu. Reference: Eta Kappa Nu, UCLA Iota Gamma Chapter, Introduction to MATLAB,
Introduction to MATLAB January 18, 2008 Steve Gu Reference: Eta Kappa Nu, UCLA Iota Gamma Chapter, Introduction to MATLAB, Part I: Basics MATLAB Environment Getting Help Variables Vectors, Matrices, and
More information03 - Basic Linear Algebra and 2D Transformations
03 - Basic Linear Algebra and 2D Transformations (invited lecture by Dr. Marcel Campen) Overview In this box, you will find references to Eigen We will briefly overview the basic linear algebra concepts
More informationRotational motion of a rigid body spinning around a rotational axis ˆn;
Physics 106a, Caltech 15 November, 2018 Lecture 14: Rotations The motion of solid bodies So far, we have been studying the motion of point particles, which are essentially just translational. Bodies with
More informationIntroduction. How to use this book. Linear algebra. Mathematica. Mathematica cells
Introduction How to use this book This guide is meant as a standard reference to definitions, examples, and Mathematica techniques for linear algebra. Complementary material can be found in the Help sections
More informationMaximal Entropy for Reconstruction of Back Projection Images
Maximal Entropy for Reconstruction of Back Projection Images Tryphon Georgiou Department of Electrical and Computer Engineering University of Minnesota Minneapolis, MN 5545 Peter J Olver Department of
More information