Shape functions in 1D

Size: px
Start display at page:

Download "Shape functions in 1D"

Transcription

1 MAE 44 & CIV 44 Introuction to Finite Elements Reaing assignment: ecture notes, ogan.,. Summary: Prof. Suvranu De Shape functions in D inear shape functions in D Quaratic an higher orer shape functions Approimation of strains an stresses in an element Aially loae elastic bar y F A() cross section at b() boy force istribution (force per unit length) E() Young s moulus Potential energy of the aially loae bar corresponing to the eact solution u() Π(u) u EA bu Fu( ) Finite element formulation, takes as its starting point, not the strong formulation, but the Principle of Minimum Potential Energy. Task is to fin the function w that minimizes the potential energy of the system Π(w) EA bw Fw( ) From the Principle of Minimum Potential Energy, that function w is the eact solution.

2 Rayleigh-Ritz Principle Step. Assume a solution w( ) a ϕo ( ) + aϕ ( ) + aϕ ( ) +... Where ϕ o (), ϕ (), are amissible functions an a o, a, etc are constants to be etermine. Step. Plug the approimate solution into the potential energy The approimate solution is u( ) a ϕo ( ) + aϕ ( ) + aϕ ( ) +... Where the coefficients have been obtaine from step Π(w) EA bw Fw( ) Step. Obtain the coefficients a o, a, etc by setting Π(w), a i i,,,... ee to fin a systematic way of choosing the approimation functions. One iea: Choose polynomials! w ( ) a Is this goo? (Is an amissible function?) Finite element iea: Step : Divie the truss into finite elements connecte to each other through special points ( noes ) 4 w ( ) a Is this goo? (Is an amissible function?) El # El # El # Total potential energysum of potential energies of the elements Π(w) EA bw Fw( )

3 Total potential energy Π(w) EA 4 El # El # El # Potential energy of element : Π(w) EA bw Potential energy of element : Π (w) EA bw bw Fw( ) 4 El # El # El # Potential energy of element : 4 4 Π (w) EA bw Fw( ) Total potential energysum of potential energies of the elements Π(w) Π(w) + Π (w) + Π (w) Step : Describe the behavior of each element Recall that in the irect stiffness approach for a bar element, we erive the stiffness matri of each element irectly (See lecture on Trusses) using the following steps: TASK : Approimate the isplacement within each bar as a straight line TASK : Approimate the strains an stresses an realize that a bar (with the approimation state in Task ) is eactly like a spring with kea/ TASK : Use the principle of force equilibrium to generate the stiffness matri ow, we will show you a systematic way of eriving the stiffness matri (sections. an. of ogan). TASK : APPROXIMATE THE DISPACEMET WITHI EACH EEMET TASK : APPROXIMATE THE STRAI an STRESS WITHI EACH EEMET TASK : DERIVE THE STIFFESS MATRIX OF EACH EEMET (net class) USIG THE PRICIPE OF MI. POT EERGY otice that the first two tasks are similar in the two methos. The only ifference is that now we are going to use the principle of minimum potential energy, rather than force equilibrium, to erive the stiffness matri.

4 TASK : APPROXIMATE THE DISPACEMET WITHI EACH EEMET Simplest assumption: isplacement varying linearly insie each bar w() a + a How to obtain a an a? w( ) a + a El # w( ) a + a w( ) a + a w( ) a + a Solve simultaneously a a Hence w() a a + () + () () () Shape functions () an () In matri notation, we write w() () Vector of noal shape functions - - [ () ()] Vector of noal isplacements OTES: PROPERTIES OF THE SHAPE FUCTIOS. Kronecker elta property: The shape function at any noe has a value of at that noe an a value of zero at A other noes. - () - () Check - () - ( ) - an ( ) El # 4

5 . Compatibility: The isplacement approimation is continuous across element bounaries () - - w () + El # El # At () - - w ( ) + () - - w () + () - - w ( ) + Hence the isplacement approimation is continuous across elements. Completeness () + () for all () + () for all - Use the epressions () ; - () An check - - () + () an () + () + Rigi boy moe () + () for all What o we mean by rigi boy moes? Assume that, this means that the element shoul translate in the positive irection by. Hence AY point () on the bar shoul have unit isplacement. et us see whether the isplacement approimation allows this. w() () + () () + () YES! Constant strain states ()+ () at all What o we mean by constant strain states? Assume that an. The strain at AY point () within the bar is ε () et us see whether the isplacement approimation allows this. w() () + () () + () w() Hence, ε () YES! 5

6 Completeness Rigi boy moes + Constant Strain states Compatibility + Completeness Convergence Ensure that the solution gets better as more elements are introuce an, in the limit, approaches the eact answer. 4. How to write the epressions for the shape functions easily (without having to erive them each time): Start with the Kronecker elta property (the shape function at any noe has value of at that noe an a value of zero at all other noes) - () - () () () ( - ) ( - ) ( - ) ( - ) El # oe at which is otice that the length of the element - ( - ) The enominator is ( - ) the numerator evaluate at the noe itself A slightly fancier assumption: isplacement varying quaratically insie each bar () () () () () () El # ( - )( - ) ( - )( - ) w() () + () + () ( - )( - ) ( - )( - ) This is a quaratic finite element in ( - )( - ) D an it has three noes an three ( - )( - ) associate shape functions per element. TASK : APPROXIMATE THE STRAI an STRESS WITHI EACH EEMET From equation (), the isplacement within each element w() Recall that the strain in the bar Hence w ε ε B () The matri B is known as the strain-isplacement matri B 6

7 For a linear finite element Hence - - [ () ()] - B - ε B - [ ] Hence, strain is a constant within each element (only for a linear element)! Displacement is linear w() a + a El # Strain is constant - ε El # Recall that the stress in the bar u σ Eε E Hence, insie the element, the approimate stress is σ EB () For a linear element the stress is also constant insie each element. This has the implication that the stress (an strain) is iscontinuous across element bounaries in general. Summary Insie an element, the three most important approimations in terms of the noal isplacements () are: Displacement approimation in terms of shape functions u() Strain approimation in terms of strain-isplacement matri ε() B () Stress approimation in terms of strain-isplacement matri an Young s moulus σ EB () () 7

8 Summary For a linear element Displacement approimation in terms of shape functions - - u() Strain approimation Stress approimation ε [ ] E σ [ ] 8

dw 2 3(w) x 4 x 4 =L x 3 x 1 =0 x 2 El #1 El #2 El #3 Potential energy of element 3: Total potential energy Potential energy of element 1:

dw 2 3(w) x 4 x 4 =L x 3 x 1 =0 x 2 El #1 El #2 El #3 Potential energy of element 3: Total potential energy Potential energy of element 1: MAE 44 & CIV 44 Introduction to Finite Elements Reading assignment: ecture notes, ogan. Summary: Pro. Suvranu De Finite element ormulation or D elasticity using the Rayleigh-Ritz Principle Stiness matri

More information

Chapter 9 Method of Weighted Residuals

Chapter 9 Method of Weighted Residuals Chapter 9 Metho of Weighte Resiuals 9- Introuction Metho of Weighte Resiuals (MWR) is an approimate technique for solving bounary value problems. It utilizes a trial functions satisfying the prescribe

More information

Lecture 3: Development of the Truss Equations.

Lecture 3: Development of the Truss Equations. 3.1 Derivation of the Stiffness Matrix for a Bar in Local Coorinates. In 3.1 we will perform Steps 1-4 of Logan s FEM. Derive the truss element equations. 1. Set the element type. 2. Select a isplacement

More information

THE ACCURATE ELEMENT METHOD: A NEW PARADIGM FOR NUMERICAL SOLUTION OF ORDINARY DIFFERENTIAL EQUATIONS

THE ACCURATE ELEMENT METHOD: A NEW PARADIGM FOR NUMERICAL SOLUTION OF ORDINARY DIFFERENTIAL EQUATIONS THE PUBISHING HOUSE PROCEEDINGS O THE ROMANIAN ACADEMY, Series A, O THE ROMANIAN ACADEMY Volume, Number /, pp. 6 THE ACCURATE EEMENT METHOD: A NEW PARADIGM OR NUMERICA SOUTION O ORDINARY DIERENTIA EQUATIONS

More information

MECh300H Introduction to Finite Element Methods. Finite Element Analysis (F.E.A.) of 1-D Problems

MECh300H Introduction to Finite Element Methods. Finite Element Analysis (F.E.A.) of 1-D Problems MECh300H Introduction to Finite Element Methods Finite Element Analysis (F.E.A.) of -D Problems Historical Background Hrenikoff, 94 frame work method Courant, 943 piecewise polynomial interpolation Turner,

More information

Section 7.1: Integration by Parts

Section 7.1: Integration by Parts Section 7.1: Integration by Parts 1. Introuction to Integration Techniques Unlike ifferentiation where there are a large number of rules which allow you (in principle) to ifferentiate any function, the

More information

Calculus in the AP Physics C Course The Derivative

Calculus in the AP Physics C Course The Derivative Limits an Derivatives Calculus in the AP Physics C Course The Derivative In physics, the ieas of the rate change of a quantity (along with the slope of a tangent line) an the area uner a curve are essential.

More information

FE Modeling and Analysis of Isotropic and Orthotropic Beams Using First Order Shear Deformation Theory

FE Modeling and Analysis of Isotropic and Orthotropic Beams Using First Order Shear Deformation Theory Materials Sciences an Applications,, 4, 77- http://.oi.org/.46/msa..4 Publishe Online January (http://www.scirp.org/journal/msa) 77 FE Moeling an Analysis of Isotropic an Orthotropic Beams Using First

More information

The thin plate theory assumes the following deformational kinematics:

The thin plate theory assumes the following deformational kinematics: MEG6007 (Fall, 2017) - Solutions of Homework # 8 (ue on Tuesay, 29 November 2017) 8.0 Variational Derivation of Thin Plate Theory The thin plate theory assumes the following eformational kinematics: u

More information

Work and Kinetic Energy

Work and Kinetic Energy Work Work an Kinetic Energy Work (W) the prouct of the force eerte on an object an the istance the object moes in the irection of the force (constant force only). W = " = cos" (N " m = J)! is the angle

More information

Using Quasi-Newton Methods to Find Optimal Solutions to Problematic Kriging Systems

Using Quasi-Newton Methods to Find Optimal Solutions to Problematic Kriging Systems Usin Quasi-Newton Methos to Fin Optimal Solutions to Problematic Kriin Systems Steven Lyster Centre for Computational Geostatistics Department of Civil & Environmental Enineerin University of Alberta Solvin

More information

EXPONENTIAL FOURIER INTEGRAL TRANSFORM METHOD FOR STRESS ANALYSIS OF BOUNDARY LOAD ON SOIL

EXPONENTIAL FOURIER INTEGRAL TRANSFORM METHOD FOR STRESS ANALYSIS OF BOUNDARY LOAD ON SOIL Tome XVI [18] Fascicule 3 [August] 1. Charles Chinwuba IKE EXPONENTIAL FOURIER INTEGRAL TRANSFORM METHOD FOR STRESS ANALYSIS OF BOUNDARY LOAD ON SOIL 1. Department of Civil Engineering, Enugu State University

More information

WJEC Core 2 Integration. Section 1: Introduction to integration

WJEC Core 2 Integration. Section 1: Introduction to integration WJEC Core Integration Section : Introuction to integration Notes an Eamples These notes contain subsections on: Reversing ifferentiation The rule for integrating n Fining the arbitrary constant Reversing

More information

Strength Analysis of CFRP Composite Material Considering Multiple Fracture Modes

Strength Analysis of CFRP Composite Material Considering Multiple Fracture Modes 5--XXXX Strength Analysis of CFRP Composite Material Consiering Multiple Fracture Moes Author, co-author (Do NOT enter this information. It will be pulle from participant tab in MyTechZone) Affiliation

More information

Formal Methods for Deriving Element Equations

Formal Methods for Deriving Element Equations Formal Methods for Deriving Element Eqations And the importance of Shape Fnctions Formal Methods In previos lectres we obtained a bar element s stiffness eqations sing the Direct Method to obtain eact

More information

SYNCHRONOUS SEQUENTIAL CIRCUITS

SYNCHRONOUS SEQUENTIAL CIRCUITS CHAPTER SYNCHRONOUS SEUENTIAL CIRCUITS Registers an counters, two very common synchronous sequential circuits, are introuce in this chapter. Register is a igital circuit for storing information. Contents

More information

Summary: Differentiation

Summary: Differentiation Techniques of Differentiation. Inverse Trigonometric functions The basic formulas (available in MF5 are: Summary: Differentiation ( sin ( cos The basic formula can be generalize as follows: Note: ( sin

More information

Table of Common Derivatives By David Abraham

Table of Common Derivatives By David Abraham Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec

More information

05 The Continuum Limit and the Wave Equation

05 The Continuum Limit and the Wave Equation Utah State University DigitalCommons@USU Founations of Wave Phenomena Physics, Department of 1-1-2004 05 The Continuum Limit an the Wave Equation Charles G. Torre Department of Physics, Utah State University,

More information

Unit #6 - Families of Functions, Taylor Polynomials, l Hopital s Rule

Unit #6 - Families of Functions, Taylor Polynomials, l Hopital s Rule Unit # - Families of Functions, Taylor Polynomials, l Hopital s Rule Some problems an solutions selecte or aapte from Hughes-Hallett Calculus. Critical Points. Consier the function f) = 54 +. b) a) Fin

More information

Lectures HD#14 Dynamic response of continuous systems. Date: April

Lectures HD#14 Dynamic response of continuous systems. Date: April ectures -3 Date: April 4 7 Toay: Vibrations of continuous systems HD#4 Dynamic response of continuous systems Free vibrations of elastic bars an beams. Properties of normal moe functions. Force response

More information

Structural Analysis of Truss Structures using Stiffness Matrix. Dr. Nasrellah Hassan Ahmed

Structural Analysis of Truss Structures using Stiffness Matrix. Dr. Nasrellah Hassan Ahmed Structural Analysis of Truss Structures using Stiffness Matrix Dr. Nasrellah Hassan Ahmed FUNDAMENTAL RELATIONSHIPS FOR STRUCTURAL ANALYSIS In general, there are three types of relationships: Equilibrium

More information

Lecture 20: Isoparametric Formulations.

Lecture 20: Isoparametric Formulations. Chapter #0 Isoparametric Formulation. Isoparametric formulations help us solve two problems. Help simplify the definition of the approimate displacement field for more comple planar elements (4-sided elements,

More information

Calculus of Variations

Calculus of Variations Calculus of Variations Lagrangian formalism is the main tool of theoretical classical mechanics. Calculus of Variations is a part of Mathematics which Lagrangian formalism is base on. In this section,

More information

Conservation laws a simple application to the telegraph equation

Conservation laws a simple application to the telegraph equation J Comput Electron 2008 7: 47 51 DOI 10.1007/s10825-008-0250-2 Conservation laws a simple application to the telegraph equation Uwe Norbrock Reinhol Kienzler Publishe online: 1 May 2008 Springer Scienceusiness

More information

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent

More information

Final Exam Study Guide and Practice Problems Solutions

Final Exam Study Guide and Practice Problems Solutions Final Exam Stuy Guie an Practice Problems Solutions Note: These problems are just some of the types of problems that might appear on the exam. However, to fully prepare for the exam, in aition to making

More information

Development of Truss Equations

Development of Truss Equations MANE & CIV Introuction to Finite Eeents Prof. Suvranu De Deveopent of Truss Equations Reaing assignent: Chapter : Sections.-.9 + ecture notes Suar: Stiffness atri of a bar/truss eeent Coorinate transforation

More information

Summary of the Class before Exam1

Summary of the Class before Exam1 uar o the lass beore Ea Builing a FEA Moel Ingreients o a FEA sotware pacage teps in builing a FEA oel Moeling consierations D pring/truss Eleents ingle D spring/truss eleent Global stiness atri; properties

More information

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University

More information

five moments Moments Moments Moments 8 a force acting at a different point causes a different moment:

five moments Moments Moments Moments 8 a force acting at a different point causes a different moment: ELEENTS O RCHITECTURL STRUCTURES: OR, EHVIOR, ND DESIGN DR. NNE NICHOLS SPRING 2014 forces have the tenency to make a boy rotate about an ais lecture five moments http://www.physics.um.eu same translation

More information

Nonlinear Dielectric Response of Periodic Composite Materials

Nonlinear Dielectric Response of Periodic Composite Materials onlinear Dielectric Response of Perioic Composite aterials A.G. KOLPAKOV 3, Bl.95, 9 th ovember str., ovosibirsk, 639 Russia the corresponing author e-mail: agk@neic.nsk.su, algk@ngs.ru A. K.TAGATSEV Ceramics

More information

Planar sheath and presheath

Planar sheath and presheath 5/11/1 Flui-Poisson System Planar sheath an presheath 1 Planar sheath an presheath A plasma between plane parallel walls evelops a positive potential which equalizes the rate of loss of electrons an ions.

More information

160 CHAPTER 2 Axially Loaded Numbers. Stress Concentrations

160 CHAPTER 2 Axially Loaded Numbers. Stress Concentrations 160 HTR xially oae Numers Stress oncentrations The prolems for Section.10 are to e solve y consiering the stress-concentration factors an assuming linearly elastic ehavior. rolem.10-1 The flat ars shown

More information

Integration Review. May 11, 2013

Integration Review. May 11, 2013 Integration Review May 11, 2013 Goals: Review the funamental theorem of calculus. Review u-substitution. Review integration by parts. Do lots of integration eamples. 1 Funamental Theorem of Calculus In

More information

Chapter 2 Lagrangian Modeling

Chapter 2 Lagrangian Modeling Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie

More information

Chapter 1 Overview: Review of Derivatives

Chapter 1 Overview: Review of Derivatives Chapter Overview: Review of Derivatives The purpose of this chapter is to review the how of ifferentiation. We will review all the erivative rules learne last year in PreCalculus. In the net several chapters,

More information

3.7 Implicit Differentiation -- A Brief Introduction -- Student Notes

3.7 Implicit Differentiation -- A Brief Introduction -- Student Notes Fin these erivatives of these functions: y.7 Implicit Differentiation -- A Brief Introuction -- Stuent Notes tan y sin tan = sin y e = e = Write the inverses of these functions: y tan y sin How woul we

More information

Implicit Differentiation

Implicit Differentiation Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,

More information

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations Lecture XII Abstract We introuce the Laplace equation in spherical coorinates an apply the metho of separation of variables to solve it. This will generate three linear orinary secon orer ifferential equations:

More information

2.13 Variation and Linearisation of Kinematic Tensors

2.13 Variation and Linearisation of Kinematic Tensors Section.3.3 Variation an Linearisation of Kinematic ensors.3. he Variation of Kinematic ensors he Variation In this section is reviewe the concept of the variation, introce in Part I, 8.5. he variation

More information

Math 3A Midterm 1 Solutions

Math 3A Midterm 1 Solutions Math 3A Miterm Solutions Rea all of the following information before starting the exam: 0/0/00 Check your exam to make sure all pages are present. When you use a major theorem (like the inermeiate value

More information

23 Implicit differentiation

23 Implicit differentiation 23 Implicit ifferentiation 23.1 Statement The equation y = x 2 + 3x + 1 expresses a relationship between the quantities x an y. If a value of x is given, then a corresponing value of y is etermine. For

More information

Mechanics Physics 151

Mechanics Physics 151 Mechanics Physics 151 Lecture 3 Continuous Systems an Fiels (Chapter 13) Where Are We Now? We ve finishe all the essentials Final will cover Lectures 1 through Last two lectures: Classical Fiel Theory

More information

Chapter 4. Electrostatics of Macroscopic Media

Chapter 4. Electrostatics of Macroscopic Media Chapter 4. Electrostatics of Macroscopic Meia 4.1 Multipole Expansion Approximate potentials at large istances 3 x' x' (x') x x' x x Fig 4.1 We consier the potential in the far-fiel region (see Fig. 4.1

More information

Construction of the Electronic Radial Wave Functions and Probability Distributions of Hydrogen-like Systems

Construction of the Electronic Radial Wave Functions and Probability Distributions of Hydrogen-like Systems Construction of the Electronic Raial Wave Functions an Probability Distributions of Hyrogen-like Systems Thomas S. Kuntzleman, Department of Chemistry Spring Arbor University, Spring Arbor MI 498 tkuntzle@arbor.eu

More information

Vibration Analysis of Railway Tracks Forced by Distributed Moving Loads

Vibration Analysis of Railway Tracks Forced by Distributed Moving Loads IJR International Journal of Railway Vol. 6, No. 4 / December 13, pp. 155-159 The Korean Society for Railway Vibration Analysis of Railway Tracks Force by Distribute Moving Loas Sinyeob Lee*, Dongkyu Kim*,

More information

Code_Aster. Detection of the singularities and computation of a card of size of elements

Code_Aster. Detection of the singularities and computation of a card of size of elements Titre : Détection es singularités et calcul une carte [...] Date : 0/0/0 Page : /6 Responsable : Josselin DLMAS Clé : R4.0.04 Révision : 9755 Detection of the singularities an computation of a car of size

More information

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21 Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting

More information

CHAPTER 1 INTRODUCTION TO THE FINITE ELEMENT METHOD AND BAR ELEMENTS...

CHAPTER 1 INTRODUCTION TO THE FINITE ELEMENT METHOD AND BAR ELEMENTS... able of Contents CHPER INRODUCION O HE FINIE EEMEN MEHOD ND BR EEMENS.... Simple springs.... Bar elements....3 System analysis... 4.3. Eample ially loae bars with varying cross-sections... 5.4 Properties

More information

Calculus Class Notes for the Combined Calculus and Physics Course Semester I

Calculus Class Notes for the Combined Calculus and Physics Course Semester I Calculus Class Notes for the Combine Calculus an Physics Course Semester I Kelly Black December 14, 2001 Support provie by the National Science Founation - NSF-DUE-9752485 1 Section 0 2 Contents 1 Average

More information

Math 1271 Solutions for Fall 2005 Final Exam

Math 1271 Solutions for Fall 2005 Final Exam Math 7 Solutions for Fall 5 Final Eam ) Since the equation + y = e y cannot be rearrange algebraically in orer to write y as an eplicit function of, we must instea ifferentiate this relation implicitly

More information

Homework 7 Due 18 November at 6:00 pm

Homework 7 Due 18 November at 6:00 pm Homework 7 Due 18 November at 6:00 pm 1. Maxwell s Equations Quasi-statics o a An air core, N turn, cylinrical solenoi of length an raius a, carries a current I Io cos t. a. Using Ampere s Law, etermine

More information

Integration: Using the chain rule in reverse

Integration: Using the chain rule in reverse Mathematics Learning Centre Integration: Using the chain rule in reverse Mary Barnes c 999 University of Syney Mathematics Learning Centre, University of Syney Using the Chain Rule in Reverse Recall that

More information

The Sokhotski-Plemelj Formula

The Sokhotski-Plemelj Formula hysics 25 Winter 208 The Sokhotski-lemelj Formula. The Sokhotski-lemelj formula The Sokhotski-lemelj formula is a relation between the following generalize functions (also calle istributions), ±iǫ = iπ(),

More information

Experimental Robustness Study of a Second-Order Sliding Mode Controller

Experimental Robustness Study of a Second-Order Sliding Mode Controller Experimental Robustness Stuy of a Secon-Orer Sliing Moe Controller Anré Blom, Bram e Jager Einhoven University of Technology Department of Mechanical Engineering P.O. Box 513, 5600 MB Einhoven, The Netherlans

More information

Optimal LQR Control of Structures using Linear Modal Model

Optimal LQR Control of Structures using Linear Modal Model Optimal LQR Control of Structures using Linear Moal Moel I. Halperin,2, G. Agranovich an Y. Ribakov 2 Department of Electrical an Electronics Engineering 2 Department of Civil Engineering Faculty of Engineering,

More information

Cable holds system BUT at t=0 it breaks!! θ=20. Copyright Luis San Andrés (2010) 1

Cable holds system BUT at t=0 it breaks!! θ=20. Copyright Luis San Andrés (2010) 1 EAMPLE # for MEEN 363 SPRING 6 Objectives: a) To erive EOMS of a DOF system b) To unerstan concept of static equilibrium c) To learn the correct usage of physical units (US system) ) To calculate natural

More information

A simple model for the small-strain behaviour of soils

A simple model for the small-strain behaviour of soils A simple moel for the small-strain behaviour of soils José Jorge Naer Department of Structural an Geotechnical ngineering, Polytechnic School, University of São Paulo 05508-900, São Paulo, Brazil, e-mail:

More information

First Order Linear Differential Equations

First Order Linear Differential Equations LECTURE 8 First Orer Linear Differential Equations We now turn our attention to the problem of constructing analytic solutions of ifferential equations; that is to say,solutions that can be epresse in

More information

Mathcad Lecture #5 In-class Worksheet Plotting and Calculus

Mathcad Lecture #5 In-class Worksheet Plotting and Calculus Mathca Lecture #5 In-class Worksheet Plotting an Calculus At the en of this lecture, you shoul be able to: graph expressions, functions, an matrices of ata format graphs with titles, legens, log scales,

More information

5.4 Fundamental Theorem of Calculus Calculus. Do you remember the Fundamental Theorem of Algebra? Just thought I'd ask

5.4 Fundamental Theorem of Calculus Calculus. Do you remember the Fundamental Theorem of Algebra? Just thought I'd ask 5.4 FUNDAMENTAL THEOREM OF CALCULUS Do you remember the Funamental Theorem of Algebra? Just thought I' ask The Funamental Theorem of Calculus has two parts. These two parts tie together the concept of

More information

18 EVEN MORE CALCULUS

18 EVEN MORE CALCULUS 8 EVEN MORE CALCULUS Chapter 8 Even More Calculus Objectives After stuing this chapter you shoul be able to ifferentiate an integrate basic trigonometric functions; unerstan how to calculate rates of change;

More information

Lecture notes Models of Mechanics

Lecture notes Models of Mechanics Lecture notes Models of Mechanics Anders Klarbring Division of Mechanics, Linköping University, Sweden Lecture 7: Small deformation theories Klarbring (Mechanics, LiU) Lecture notes Linköping 2012 1 /

More information

Damage identification based on incomplete modal data and constrained nonlinear multivariable function

Damage identification based on incomplete modal data and constrained nonlinear multivariable function Journal of Physics: Conference Series PAPER OPEN ACCESS Damage ientification base on incomplete moal ata an constraine nonlinear multivariable function To cite this article: S S Kourehli 215 J. Phys.:

More information

3.2 Differentiability

3.2 Differentiability Section 3 Differentiability 09 3 Differentiability What you will learn about How f (a) Might Fail to Eist Differentiability Implies Local Linearity Numerical Derivatives on a Calculator Differentiability

More information

LINEAR DIFFERENTIAL EQUATIONS OF ORDER 1. where a(x) and b(x) are functions. Observe that this class of equations includes equations of the form

LINEAR DIFFERENTIAL EQUATIONS OF ORDER 1. where a(x) and b(x) are functions. Observe that this class of equations includes equations of the form LINEAR DIFFERENTIAL EQUATIONS OF ORDER 1 We consier ifferential equations of the form y + a()y = b(), (1) y( 0 ) = y 0, where a() an b() are functions. Observe that this class of equations inclues equations

More information

3.6. Let s write out the sample space for this random experiment:

3.6. Let s write out the sample space for this random experiment: STAT 5 3 Let s write out the sample space for this ranom experiment: S = {(, 2), (, 3), (, 4), (, 5), (2, 3), (2, 4), (2, 5), (3, 4), (3, 5), (4, 5)} This sample space assumes the orering of the balls

More information

Crack onset assessment near the sharp material inclusion tip by means of modified maximum tangential stress criterion

Crack onset assessment near the sharp material inclusion tip by means of modified maximum tangential stress criterion Focuse on Mechanical Fatigue of Metals Crack onset assessment near the sharp material inclusion tip by means of moifie maximum tangential stress criterion Onřej Krepl, Jan Klusák CEITEC IPM, Institute

More information

Mathematics 116 HWK 25a Solutions 8.6 p610

Mathematics 116 HWK 25a Solutions 8.6 p610 Mathematics 6 HWK 5a Solutions 8.6 p6 Problem, 8.6, p6 Fin a power series representation for the function f() = etermine the interval of convergence. an Solution. Begin with the geometric series = + +

More information

Quantum Mechanics in Three Dimensions

Quantum Mechanics in Three Dimensions Physics 342 Lecture 20 Quantum Mechanics in Three Dimensions Lecture 20 Physics 342 Quantum Mechanics I Monay, March 24th, 2008 We begin our spherical solutions with the simplest possible case zero potential.

More information

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013 Survey Sampling Kosuke Imai Department of Politics, Princeton University February 19, 2013 Survey sampling is one of the most commonly use ata collection methos for social scientists. We begin by escribing

More information

Average value of position for the anharmonic oscillator: Classical versus quantum results

Average value of position for the anharmonic oscillator: Classical versus quantum results verage value of position for the anharmonic oscillator: Classical versus quantum results R. W. Robinett Department of Physics, The Pennsylvania State University, University Park, Pennsylvania 682 Receive

More information

Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum

Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum October 6, 4 ARDB Note Analytic Scaling Formulas for Crosse Laser Acceleration in Vacuum Robert J. Noble Stanfor Linear Accelerator Center, Stanfor University 575 San Hill Roa, Menlo Park, California 945

More information

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5

More information

Momentum and Energy. Chapter Conservation Principles

Momentum and Energy. Chapter Conservation Principles Chapter 2 Momentum an Energy In this chapter we present some funamental results of continuum mechanics. The formulation is base on the principles of conservation of mass, momentum, angular momentum, an

More information

Calculus of Variations

Calculus of Variations 16.323 Lecture 5 Calculus of Variations Calculus of Variations Most books cover this material well, but Kirk Chapter 4 oes a particularly nice job. x(t) x* x*+ αδx (1) x*- αδx (1) αδx (1) αδx (1) t f t

More information

8.022 (E&M) Lecture 19

8.022 (E&M) Lecture 19 8. (E&M) Lecture 19 Topics: The missing term in Maxwell s equation Displacement current: what it is, why it s useful The complete Maxwell s equations An their solution in vacuum: EM waves Maxwell s equations

More information

Static analysis of beams on elastic foundation by the method of discrete singular convolution

Static analysis of beams on elastic foundation by the method of discrete singular convolution International Journal of Engineering & Applie Sciences (IJEAS) Beir Agöz, Kair Mercan, Çiğem Demir, Ömer Civale Vol.8, Issue 3(016) 67-73 Static analysis of beams on elastic founation by the metho of iscrete

More information

Chapter 2 Governing Equations

Chapter 2 Governing Equations Chapter 2 Governing Equations In the present an the subsequent chapters, we shall, either irectly or inirectly, be concerne with the bounary-layer flow of an incompressible viscous flui without any involvement

More information

Thermal Modulation of Rayleigh-Benard Convection

Thermal Modulation of Rayleigh-Benard Convection Thermal Moulation of Rayleigh-Benar Convection B. S. Bhaauria Department of Mathematics an Statistics, Jai Narain Vyas University, Johpur, Inia-3400 Reprint requests to Dr. B. S.; E-mail: bsbhaauria@reiffmail.com

More information

Polynomial Inclusion Functions

Polynomial Inclusion Functions Polynomial Inclusion Functions E. e Weert, E. van Kampen, Q. P. Chu, an J. A. Muler Delft University of Technology, Faculty of Aerospace Engineering, Control an Simulation Division E.eWeert@TUDelft.nl

More information

ACCURATE BENDING STRENGTH ANALYSIS OF ASYMMETRIC GEARS USING THE NOVEL ES-PIM WITH TRIANGULAR MESH

ACCURATE BENDING STRENGTH ANALYSIS OF ASYMMETRIC GEARS USING THE NOVEL ES-PIM WITH TRIANGULAR MESH International Journal of Automotive an Mechanical Engineering (IJAME) ISSN: 2229-8649 (Print); ISSN: 2180-1606 (Online); Volume 4, pp. 373-396, July-December 2011 Universiti Malaysia Pahang DOI: http://x.oi.org/10.15282/ijame.4.2011.1.0031

More information

PART I. Basic Concepts

PART I. Basic Concepts PART I. Basic Concepts. Introduction. Basic Terminology of Structural Vibration.. Common Vibration Sources.. Forms of Vibration.3 Structural Vibration . Basic Terminology of Structural Vibration The term

More information

TAYLOR S POLYNOMIAL APPROXIMATION FOR FUNCTIONS

TAYLOR S POLYNOMIAL APPROXIMATION FOR FUNCTIONS MISN-0-4 TAYLOR S POLYNOMIAL APPROXIMATION FOR FUNCTIONS f(x ± ) = f(x) ± f ' (x) + f '' (x) 2 ±... 1! 2! = 1.000 ± 0.100 + 0.005 ±... TAYLOR S POLYNOMIAL APPROXIMATION FOR FUNCTIONS by Peter Signell 1.

More information

Chapter 14 Truss Analysis Using the Stiffness Method

Chapter 14 Truss Analysis Using the Stiffness Method Chapter 14 Truss Analsis Using the Stiffness Method Structural Mechanics 2 ept of Arch Eng, Ajou Univ Outline undamentals of the stiffness method Member stiffness matri isplacement and force transformation

More information

CS9840 Learning and Computer Vision Prof. Olga Veksler. Lecture 2. Some Concepts from Computer Vision Curse of Dimensionality PCA

CS9840 Learning and Computer Vision Prof. Olga Veksler. Lecture 2. Some Concepts from Computer Vision Curse of Dimensionality PCA CS9840 Learning an Computer Vision Prof. Olga Veksler Lecture Some Concepts from Computer Vision Curse of Dimensionality PCA Some Slies are from Cornelia, Fermüller, Mubarak Shah, Gary Braski, Sebastian

More information

Chapter 3 Variational Formulation & the Galerkin Method

Chapter 3 Variational Formulation & the Galerkin Method Institute of Structural Engineering Page 1 Chapter 3 Variational Formulation & the Galerkin Method Institute of Structural Engineering Page 2 Today s Lecture Contents: Introduction Differential formulation

More information

Crack-tip stress evaluation of multi-scale Griffith crack subjected to

Crack-tip stress evaluation of multi-scale Griffith crack subjected to Crack-tip stress evaluation of multi-scale Griffith crack subjecte to tensile loaing by using periynamics Xiao-Wei Jiang, Hai Wang* School of Aeronautics an Astronautics, Shanghai Jiao Tong University,

More information

Finite 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 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 information

3-D FEM Modeling of fiber/matrix interface debonding in UD composites including surface effects

3-D FEM Modeling of fiber/matrix interface debonding in UD composites including surface effects IOP Conference Series: Materials Science an Engineering 3-D FEM Moeling of fiber/matrix interface eboning in UD composites incluing surface effects To cite this article: A Pupurs an J Varna 2012 IOP Conf.

More information

Linear First-Order Equations

Linear First-Order Equations 5 Linear First-Orer Equations Linear first-orer ifferential equations make up another important class of ifferential equations that commonly arise in applications an are relatively easy to solve (in theory)

More information

Physics 505 Electricity and Magnetism Fall 2003 Prof. G. Raithel. Problem Set 3. 2 (x x ) 2 + (y y ) 2 + (z + z ) 2

Physics 505 Electricity and Magnetism Fall 2003 Prof. G. Raithel. Problem Set 3. 2 (x x ) 2 + (y y ) 2 + (z + z ) 2 Physics 505 Electricity an Magnetism Fall 003 Prof. G. Raithel Problem Set 3 Problem.7 5 Points a): Green s function: Using cartesian coorinates x = (x, y, z), it is G(x, x ) = 1 (x x ) + (y y ) + (z z

More information

A new identification method of the supply hole discharge coefficient of gas bearings

A new identification method of the supply hole discharge coefficient of gas bearings Tribology an Design 95 A new ientification metho of the supply hole ischarge coefficient of gas bearings G. Belforte, F. Colombo, T. Raparelli, A. Trivella & V. Viktorov Department of Mechanics, Politecnico

More information

Designing Information Devices and Systems II Fall 2017 Note Theorem: Existence and Uniqueness of Solutions to Differential Equations

Designing Information Devices and Systems II Fall 2017 Note Theorem: Existence and Uniqueness of Solutions to Differential Equations EECS 6B Designing Information Devices an Systems II Fall 07 Note 3 Secon Orer Differential Equations Secon orer ifferential equations appear everywhere in the real worl. In this note, we will walk through

More information

Calculus and optimization

Calculus and optimization Calculus an optimization These notes essentially correspon to mathematical appenix 2 in the text. 1 Functions of a single variable Now that we have e ne functions we turn our attention to calculus. A function

More information

Development of Truss Equations

Development of Truss Equations CIVL 7/87 Chapter 3 - Truss Equations - Part /53 Chapter 3a Development of Truss Equations Learning Objectives To derive the stiffness matri for a bar element. To illustrate how to solve a bar assemblage

More information

Static Equilibrium. Theory: The conditions for the mechanical equilibrium of a rigid body are (a) (b)

Static Equilibrium. Theory: The conditions for the mechanical equilibrium of a rigid body are (a) (b) LPC Physics A 00 Las Positas College, Physics Department Staff Purpose: To etermine that, for a boy in equilibrium, the following are true: The sum of the torques about any point is zero The sum of forces

More information

Chapter Primer on Differentiation

Chapter Primer on Differentiation Capter 0.01 Primer on Differentiation After reaing tis capter, you soul be able to: 1. unerstan te basics of ifferentiation,. relate te slopes of te secant line an tangent line to te erivative of a function,.

More information

ARCH 614 Note Set 5 S2012abn. Moments & Supports

ARCH 614 Note Set 5 S2012abn. Moments & Supports RCH 614 Note Set 5 S2012abn Moments & Supports Notation: = perpenicular istance to a force from a point = name for force vectors or magnitue of a force, as is P, Q, R x = force component in the x irection

More information