CHAPTER 4 EVALUATION OF FORCE-CONSTANT MATRIX

Similar documents
PHYS 705: Classical Mechanics. Derivation of Lagrange Equations from D Alembert s Principle

3. A Review of Some Existing AW (BT, CT) Algorithms

Physics 2A Chapter 11 - Universal Gravitation Fall 2017

24-2: Electric Potential Energy. 24-1: What is physics

Integral Vector Operations and Related Theorems Applications in Mechanics and E&M

Scalars and Vectors Scalar

Density Functional Theory I

Set of square-integrable function 2 L : function space F

VEKTORANALYS FLUX INTEGRAL LINE INTEGRAL. and. Kursvecka 2. Kapitel 4 5. Sidor 29 50

Thermodynamics of solids 4. Statistical thermodynamics and the 3 rd law. Kwangheon Park Kyung Hee University Department of Nuclear Engineering

Chapter I Matrices, Vectors, & Vector Calculus 1-1, 1-9, 1-10, 1-11, 1-17, 1-18, 1-25, 1-27, 1-36, 1-37, 1-41.

Review of Vector Algebra and Vector Calculus Operations

Mechanics Physics 151

COMPLEMENTARY ENERGY METHOD FOR CURVED COMPOSITE BEAMS

Machine Learning. Spectral Clustering. Lecture 23, April 14, Reading: Eric Xing 1

Part V: Velocity and Acceleration Analysis of Mechanisms

Contact, information, consultations

19 The Born-Oppenheimer Approximation

CSU ATS601 Fall Other reading: Vallis 2.1, 2.2; Marshall and Plumb Ch. 6; Holton Ch. 2; Schubert Ch r or v i = v r + r (3.

Chapter 23: Electric Potential

CSJM University Class: B.Sc.-II Sub:Physics Paper-II Title: Electromagnetics Unit-1: Electrostatics Lecture: 1 to 4

Complex atoms and the Periodic System of the elements

Rigid Bodies: Equivalent Systems of Forces

Physics 11b Lecture #2. Electric Field Electric Flux Gauss s Law

Stellar Astrophysics. dt dr. GM r. The current model for treating convection in stellar interiors is called mixing length theory:

Energy in Closed Systems

Hybridization Based Reachability of Uncertain Planar Affine Systems

Multistage Median Ranked Set Sampling for Estimating the Population Median

Instantaneous velocity field of a round jet

Inseting this into the left hand side of the equation of motion above gives the most commonly used algoithm in classical molecula dynamics simulations

4.4 Continuum Thermomechanics

anubhavclasses.wordpress.com CBSE Solved Test Papers PHYSICS Class XII Chapter : Electrostatics

3.1 Electrostatic Potential Energy and Potential Difference

Physical & Interfacial Electrochemistry 2013

Remember: When an object falls due to gravity its potential energy decreases.

Analysis of Truss Structures with Uncertainties: From Experimental Data to Analytical Responses

Multipole Radiation. March 17, 2014

Electrostatics. 1. Show does the force between two point charges change if the dielectric constant of the medium in which they are kept increase?

gravity r2,1 r2 r1 by m 2,1

Engineering Mechanics. Force resultants, Torques, Scalar Products, Equivalent Force systems

UNIT10 PLANE OF REGRESSION

Molecular Dynamic Simulations of Nickel Nanowires at Various Temperatures

Using DP for hierarchical discretization of continuous attributes. Amit Goyal (31 st March 2008)

Correspondence Analysis & Related Methods

Physics 202, Lecture 2. Announcements

Chapter 8. Linear Momentum, Impulse, and Collisions

Chapter 12: Kinematics of a Particle 12.8 CURVILINEAR MOTION: CYLINDRICAL COMPONENTS. u of the polar coordinate system are also shown in

Test 1 phy What mass of a material with density ρ is required to make a hollow spherical shell having inner radius r i and outer radius r o?

V. Principles of Irreversible Thermodynamics. s = S - S 0 (7.3) s = = - g i, k. "Flux": = da i. "Force": = -Â g a ik k = X i. Â J i X i (7.

As is natural, our Aerospace Structures will be described in a Euclidean three-dimensional space R 3.

Amplifier Constant Gain and Noise

MULTIPOLE FIELDS. Multipoles, 2 l poles. Monopoles, dipoles, quadrupoles, octupoles... Electric Dipole R 1 R 2. P(r,θ,φ) e r

P 365. r r r )...(1 365

Tian Zheng Department of Statistics Columbia University

A. Thicknesses and Densities

Unit_III Complex Numbers: Some Basic Results: 1. If z = x +iy is a complex number, then the complex number z = x iy is

Chapter Fifiteen. Surfaces Revisited

A New Approach for Deriving the Instability Potential for Plates Based on Rigid Body and Force Equilibrium Considerations

The Forming Theory and the NC Machining for The Rotary Burs with the Spectral Edge Distribution

Photodisintegration of light nuclei

COORDINATE SYSTEMS, COORDINATE TRANSFORMS, AND APPLICATIONS

Dynamics of Rigid Bodies

Summer Workshop on the Reaction Theory Exercise sheet 8. Classwork

ISOPARAMETRIC ELEMENTS

8 Baire Category Theorem and Uniform Boundedness

Qualifying Examination Electricity and Magnetism Solutions January 12, 2006

Physics Exam II Chapters 25-29

A New Fuzzy Control Model for Kidney Patients

Capítulo. Three Dimensions

ˆ x ESTIMATOR. state vector estimate

Interatomic Forces. Overview

PHYS Week 5. Reading Journals today from tables. WebAssign due Wed nite

Solving the Dirac Equation: Using Fourier Transform

Asymptotic Solutions of the Kinetic Boltzmann Equation and Multicomponent Non-Equilibrium Gas Dynamics

DYNAMICS VECTOR MECHANICS FOR ENGINEERS: Kinematics of Rigid Bodies in Three Dimensions. Seventh Edition CHAPTER

Solutions to selected problems from homework 1.

On the Essence of Electric Charge Part 1

Supersymmetry in Disorder and Chaos (Random matrices, physics of compound nuclei, mathematics of random processes)

Physics 235 Chapter 5. Chapter 5 Gravitation

Review: Electrostatics and Magnetostatics

Physica A 392 (2013) Contents lists available at SciVerse ScienceDirect. Physica A. journal homepage:

J. N. R E DDY ENERGY PRINCIPLES AND VARIATIONAL METHODS APPLIED MECHANICS

Some Approximate Analytical Steady-State Solutions for Cylindrical Fin

SOME NEW SELF-DUAL [96, 48, 16] CODES WITH AN AUTOMORPHISM OF ORDER 15. KEYWORDS: automorphisms, construction, self-dual codes

Slide 1. Quantum Mechanics: the Practice

Exam 1. Sept. 22, 8:00-9:30 PM EE 129. Material: Chapters 1-8 Labs 1-3

To Feel a Force Chapter 7 Static equilibrium - torque and friction

ABSTRACT. JIN, YUN Atomistic Simulations of Fracture of 2D Graphene Systems and the Elastic

Rotating Disk Electrode -a hydrodynamic method

arxiv: v1 [physics.pop-ph] 3 Jun 2013

Fundamental principles

Chapter 5 Linear Equations: Basic Theory and Practice

The geometric construction of Ewald sphere and Bragg condition:

Minimization of the Free Energy under a Given Pressure by Natural Iteration Method

Chapter IV Vector and Tensor Analysis IV.2 Vector and Tensor Analysis September 29,

2 Governing Equations

The Peak Shape of the Pair Distribution Function

An Approach to Inverse Fuzzy Arithmetic

Chapter 2. A Brief Review of Electron Diffraction Theory

Chapter 13 - Universal Gravitation

Transcription:

CHAPTER 4 EVALUATION OF FORCE-CONSTANT MATRIX 4.- AIM OF THE WORK As antcpated n the ntodcton the am of the pesent ok s to obtan the nmecal vale of the foce-constant matx fo tantalm. In a fst step expesson of the total enegy of the elemental system s taken nto accont. It conssts of pa nteacton potentals embedded potentals and angla dependent potentals. etals abot the expesson fo the potentals sch as vales of paametes and lattce popetes ae pesented n the paagaph 4.. Statng fom the expesson of the total enegy of the system t s possble to evalate the fst and second devatve of sch enegy. Examples and explcatons abot ho to evalate these devatves ae pesented n the paagaph 4.3. Once evalated the tems of the devatves t s possble to eoganze them to evalate the foce-constant matx both n the eal doman that n the complex doman va Foe tansfom. These evalatons ae obtaned thanks to the se of the comptatonal softae pogam "Wolfam Mathematca 9.0". Moe detals abot ho the foce-constant matx as evalated and ts nmecal vale ae pesented n paagaph 4.4. 4.- TOTAL POTENTIAL ENERGY OF TANTALUM Tantalm a chemcal element hch symbol s T a and atomc nmbe 73 s a body-contended cbc BBC tanston metal.

Fge 4. - Tantalm The mechancal and physcal popetes depend to a geat extent on the metal pty. Pty s manly elated to the pesence of metallods sch as cabon oxygen ntogen and hydogen of hch even lo concentaton case a consdeable changes n the popetes. One of the man popetes of tantalm s ts hgh esstance to cooson. The alloys podced th tantalm have hgh meltng ponts ae stong and have good dctlty. It s sed ethe n a pe fom o as a component n alloys n a nmbe of technologcal applcatons ncldng hgh-tempeate stctal mateals and dffson baes n ntegated ccts. In aeospace applcatons t s manly sed n the podcton of spealloys fo et engne components. In ecent yeas sgnfcant effects have been made to mpove the basc ndestandng of the mechancal behavo of tantalm by means of compte modellng.

3 In the pesent ok to evalate the foce-constant matx the statng pont s the expesson of the total potental enegy fo tantalm pesented n scvee femento and hee shong: E tot φ F µ λ 6 ν 4. hee ndces and enmeate atoms and the spescpts 3 efe to a Catesan components of vectos and tensos. As befly ntodced n chapte 3.4. bt th dffeent notaton the fst tem n eqaton 4. epesents pa nteactons beng the pa-nteacton potental beteen an atom located at postons and an atom at poston. The fncton F s the enegy of embeddng atom n the host electon densty ndced at ste by all othe atoms of the system. The host electon densty s gven by : 4. hee s the electon densty fncton assgned to an atom. The fst to tems n 4. epesent the egla EAM pesented n chapte 3.4. and have a cental foce chaacte. They depend only on nteatomc dstances bt not on bond angles. The non-cental component of bondng s ntodced thogh the dpole vectos µ 4.3 and qadpole tensos

4 λ 4.4 hee ν s the tace of λ ν λ 4.5 The angla tems n 4.3-4.5 penalze the total enegy fo devatons of local atomc envonments fom cbc symmety. The expesson of the fnctons shoed above ee gven n a nmecal fom. A table th to colmns as sed fo each fncton n hch n one colmn thee ae the vales of the ndependent vaable and n the second thee ae the coespondng vales of the fncton. These tables can be fond n the C-ROM of the ok as a "Mcosoft Excel" fle named "modelo Tanalo Mshn.xlsx" Mshn Y. Elasse C. Gmbsch P. Lozovo Y. Hashbon A. 008. The man popetes of the BCC Tantalm ae shon n the follong fge :

5 Fge 4. - Popetes of BCC Tantalm calclated th the angla dependent potental n compason th expemental fst-pncple data It s mpotant to emak that the angla dependent potental s an empcal potental sed to evalate the total enegy of the system. It ncldes also the EAM potentals ntodced n chapte 3.4. and t s one of the many potentals avalable to model the total potental enegy. Each model can be bette than the othes dependng by the specfc applcatons. The potental sed fo ths ok s the one pesented n Mshn Y. Lozovo A. Y. 006. 4.3- FIRST AN SECON ERIVATIVE OF THE TOTAL POTENTIAL ENERGY

6 As pevosly ntodced n chapte 4. the foce constant matx can be constcted statng fom the second devatve of the total potental enegy expessed by 4.. In fact developng n sees of Taylo ntl the second ode the expesson of the potental of the cystal t s possble to obtan an expesson n hch appea the elements of the foce constant matx elated to the second devatve of the total potental enegy thogh the elaton : Φ k l m E l ToT k m 0 4.6 l Whee s the dsplacement of the atom m decton k s the dsplacement of the atom decton and E ToT s the total potental enegy The elaton 4. s fo a geneal case. l n the th coodnates m n the k th coodnate Applyng the dscete lattce model pesented n chapte 3 to the total potental enegy expessed by 4. the nteatomc dstance ae consdeed -cell elements an elementay segment hle the lattce ste poston ae consdeed 0-cell elements o nteconnected vetces.

7 Once made these assmptons t s necessay to choose the nmbe of neghbos consdeed n the expesson of the potental. In the pesent model ae taken nto accont the fst and second neghbos hch mply 4 atoms of efeence as schematcally shon n the fge.9. It s mpotant to note that n the expesson of the total potental enegy thee ae tems hch depend explctly fom the nteatomc dstance o elementay segment e sch as the pa nteacton potentals; and tems that depends explctly by the ste poston bt mplctly by the nteatomc dstance. Consdeng that to evalate the foce-constant matx t s necessay to deve the total potental enegy th espect to the -cell elementay segment hen devng the tems hch depends fom the 0-cell lattce poston sch as the angla-dependent nteatomc potentals o the EAM potental the chan le has to be appled. In the pesent ok to bette manage all the tems geneated by the second devatve of the total potental enegy sch devatves ae evalated ndvdal by fo each of the fve tems hch consttte the total potental enegy. Ths fomal expesson takng nto accont the chan le ae shon belo: E φ φ 4.7

8 F E k k k k F F F ] [ 4.8 E 3 k k k k k k k k k k k k k k

9 [ ] µ µ 4.9 4 E σ σ σ σ σ σ k k k k k k k k k

0 λ λ σ σ σ σ σ σ } σ σ 4.0 E 6 5 k 6 ν ν 6

4. The notaton sed s the one ntodced n chapte 4. epesented the nteatomc dstance o -cell elementay segment denotes the Konecke delta and the symbol denotes the dyadc podct. Wth the symbol ae denoted the atoms that ae the second exteme of a -cell elementay segment hch has as fst exteme the atom. Some smplfcatons ee sed to smplfy these expessons sch as : 4. λ λ 4.3 Whee λ s a vale obtaned hen the 3x3 matx λ s constcted µ 0 4.4 4.5 4.6 ν 4.7

The fst and second devatve of the fnctons φ F and hch appea nto the expessons 4.7-4. ae evalated nmecally as the fomat as nmecal. The fomlas sed to evalate these nmecal devatves ae: f t 6 f t 3 f t f t f ' t 6h 4 5 4 f '' t [ f t.. f t f t f t f t h 3 3 The nmecal vales can be seen n C-ROM nde the name "Mshn Tantalo.xlsx". 4.4- EVALUATION OF THE FORCE CONSTANT MATRIX Once evalated the second devatve of the fve tems of the total potental enegy chapte 4.3 t s possble to constct the foce-constant matx B. Fo each of the fve tems deved t s possble to constct a foceconstant matx egopng coectly the tems of each devatve. The fomal expesson of the foce-constant matx B of each tem s: B e e B P e e 0 Q 4. 4.0

3 B e e P Q 4. B e e R k k 4.3 B 3 e e P 3 Q 3 4.4 B 3 k k k e R3 T3 M 3 3 k k e N 4.5 B 4 σ σ e e P4 σ Q4 S4 4 E4 σ σ F4 G4λ 4.6 B 4 σ σ σ k k k e R4 T4 M 4 4 k k e N σ L 4 σ k k 4.7 B 5 e e P 5 Q 5 λ 4.8

4 B 5 e e R5 k k 4.9 Whee the geneal 3x3 elements B e e epesent the elements of the B e e man dagonal and the 3x3 elements ae the othes elements. The tems P Q R P... ae nmbes that mltples the dffeent 3x3 elements of B and ae obtaned egopng all the tems that fo each second devatve mltply the 3x3 element of B consdeed. The fnal foce-constant matx s obtaned by addng the fve matces: B e e B e e B e e B3 e e B4 e e B5 e e 4.30 B 3 4 B 5 e e e B e e B e e B e e e 4.3 As an example the elements B and B ae shon : 0.649868 B 0.397763 0.397763 0.397763 0.649868 0.397763 0.397763 0.397763 0.649868 0.009089 B 0.0487 0.06 0.0487.009089 0.06 0.06 0.06 0.057455 Whee epesents the -cell elements and shon n fge 3.. The fll expesson of the matx B can be fond n the C-ROM nde the name "vefca del fle nmec secondo vecno.nb"

5 Usng the Foe tansfom t s possble to evalate the matx Ψˆ θ l Ψ l m e θl 4.3 At last sng the elaton 3. t s possble to evalate the matx Φˆ hose nmecal expesson s: 0.07907 ˆΦ 0.08899 0.00363806 0.08899 0.33309 0.0056 0.00363806 0.0056 0.0665869 hee the expesson of Φˆ s evalated fst applyng the le : θ4 θ θ θ3 ; θ5 θ θ3 ; θ 6 θ θ3 ; θ 7 θ θ and then gvng the vales : θ 0. ; θ 0. ; θ 3 0. othese the matx Φˆ depends by θ and has a vey lage expesson. All these evalatons ee made sng the softae Wolfam Mathematca 9.0 a comptatonal softae pogam.