CIVL 8/ D Boundary Value Problems - Rectangular Elements 1/7

Similar documents
INTEGRALS. Chapter 7. d dx. 7.1 Overview Let d dx F (x) = f (x). Then, we write f ( x)

TOPIC 5: INTEGRATION

Instructions for Section 1

Chapter 16. 1) is a particular point on the graph of the function. 1. y, where x y 1

MASTER CLASS PROGRAM UNIT 4 SPECIALIST MATHEMATICS SEMESTER TWO 2014 WEEK 11 WRITTEN EXAMINATION 1 SOLUTIONS

Integration Continued. Integration by Parts Solving Definite Integrals: Area Under a Curve Improper Integrals

( ) Geometric Operations and Morphing. Geometric Transformation. Forward v.s. Inverse Mapping. I (x,y ) Image Processing - Lesson 4 IDC-CG 1

Lecture 11 Waves in Periodic Potentials Today: Questions you should be able to address after today s lecture:

Mathematics. Mathematics 3. hsn.uk.net. Higher HSN23000

Limits Indeterminate Forms and L Hospital s Rule

CONIC SECTIONS. MODULE-IV Co-ordinate Geometry OBJECTIVES. Conic Sections

Miscellaneous open problems in the Regular Boundary Collocation approach

, between the vertical lines x a and x b. Given a demand curve, having price as a function of quantity, p f (x) at height k is the curve f ( x,

COHORT MBA. Exponential function. MATH review (part2) by Lucian Mitroiu. The LOG and EXP functions. Properties: e e. lim.

Section 11.6: Directional Derivatives and the Gradient Vector

SOLVED EXAMPLES. be the foci of an ellipse with eccentricity e. For any point P on the ellipse, prove that. tan

Ch 1.2: Solutions of Some Differential Equations

CONTINUITY AND DIFFERENTIABILITY

Finite element discretization of Laplace and Poisson equations

AS 5850 Finite Element Analysis

CBSE 2015 FOREIGN EXAMINATION

FSA. CmSc 365 Theory of Computation. Finite State Automata and Regular Expressions (Chapter 2, Section 2.3) ALPHABET operations: U, concatenation, *

A general N-dimensional vector consists of N values. They can be arranged as a column or a row and can be real or complex.

This Week. Computer Graphics. Introduction. Introduction. Graphics Maths by Example. Graphics Maths by Example

Chapter 5. Introduction. Introduction. Introduction. Finite Element Modelling. Finite Element Modelling

The Theory of Small Reflections

1 Isoparametric Concept

Continuous Random Variables: Basics

CSE303 - Introduction to the Theory of Computing Sample Solutions for Exercises on Finite Automata

Lecture contents. Bloch theorem k-vector Brillouin zone Almost free-electron model Bands Effective mass Holes. NNSE 508 EM Lecture #9

MASSACHUSETTS INSTITUTE OF TECHNOLOGY HAYSTACK OBSERVATORY WESTFORD, MASSACHUSETTS

Linear Algebra Existence of the determinant. Expansion according to a row.

Section 3: Antiderivatives of Formulas

Self-Adjointness and Its Relationship to Quantum Mechanics. Ronald I. Frank 2016

Introduction to Condensed Matter Physics

Multi-Section Coupled Line Couplers

HIGHER ORDER DIFFERENTIAL EQUATIONS

BSc Engineering Sciences A. Y. 2017/18 Written exam of the course Mathematical Analysis 2 August 30, x n, ) n 2

u 3 = u 3 (x 1, x 2, x 3 )

Finite Element Models for Steady Flows of Viscous Incompressible Fluids

Minimum Spanning Trees

Lie Groups HW7. Wang Shuai. November 2015

(2) If we multiplied a row of B by λ, then the value is also multiplied by λ(here lambda could be 0). namely

Garnir Polynomial and their Properties

Finite Element Method FEM FOR FRAMES

Einstein Equations for Tetrad Fields

VSMN30 FINITA ELEMENTMETODEN - DUGGA

Integration by Parts

ME 522 PRINCIPLES OF ROBOTICS. FIRST MIDTERM EXAMINATION April 19, M. Kemal Özgören

PH427/PH527: Periodic systems Spring Overview of the PH427 website (syllabus, assignments etc.) 2. Coupled oscillations.

Elliptical motion, gravity, etc

However, many atoms can combine to form particular molecules, e.g. Chlorine (Cl) and Sodium (Na) atoms form NaCl molecules.

Fundamental Algorithms for System Modeling, Analysis, and Optimization

Functions and Graphs 1. (a) (b) (c) (f) (e) (d) 2. (a) (b) (c) (d)

Background: We have discussed the PIB, HO, and the energy of the RR model. In this chapter, the H-atom, and atomic orbitals.

surface of a dielectric-metal interface. It is commonly used today for discovering the ways in

NONLINEAR ANALYSIS OF PLATE BENDING

Middle East Technical University Department of Mechanical Engineering ME 413 Introduction to Finite Element Analysis

Addition of angular momentum

Middle East Technical University Department of Mechanical Engineering ME 413 Introduction to Finite Element Analysis

ECE COMBINATIONAL BUILDING BLOCKS - INVEST 13 DECODERS AND ENCODERS

6.5 Plate Problems in Rectangular Coordinates

The Matrix Exponential

cycle that does not cross any edges (including its own), then it has at least

3 Finite Element Parametric Geometry

Walk Like a Mathematician Learning Task:

3.4 Repeated Roots; Reduction of Order

The Z transform techniques

Why the Junction Tree Algorithm? The Junction Tree Algorithm. Clique Potential Representation. Overview. Chris Williams 1.

Review of Gaussian Quadrature method

The Matrix Exponential

That is, we start with a general matrix: And end with a simpler matrix:

Lecture 12 Quantum chromodynamics (QCD) WS2010/11: Introduction to Nuclear and Particle Physics

EEO 401 Digital Signal Processing Prof. Mark Fowler

Lecture 4. Conic section

Addition of angular momentum

Fourier Transforms and the Wave Equation. Key Mathematics: More Fourier transform theory, especially as applied to solving the wave equation.


2008 AP Calculus BC Multiple Choice Exam

ME311 Machine Design

IIT JEE MATHS MATRICES AND DETERMINANTS

Notes on Finite Automata Department of Computer Science Professor Goldberg Textbooks: Introduction to the Theory of Computation by Michael Sipser

Basic Polyhedral theory

MATH 319, WEEK 15: The Fundamental Matrix, Non-Homogeneous Systems of Differential Equations

b. How many ternary words of length 23 with eight 0 s, nine 1 s and six 2 s?

Note If the candidate believes that e x = 0 solves to x = 0 or gives an extra solution of x = 0, then withhold the final accuracy mark.

Last time: introduced our first computational model the DFA.

INTRODUCTION TO LINEAR ALGEBRA

10. Limits involving infinity

(a) v 1. v a. v i. v s. (b)

Floating Point Number System -(1.3)

Floating Point Number System -(1.3)

Supplementary Materials

MERIDIAN SURFACE OF WEINGARTEN TYPE IN 4-DIMENSIONAL EUCLIDEAN SPACE E 4 arxiv: v1 [math.dg] 14 May 2013

INTEGRATION BY PARTS

Section 5.1/5.2: Areas and Distances the Definite Integral

1 Introduction to Modulo 7 Arithmetic

Engineering 323 Beautiful HW #13 Page 1 of 6 Brown Problem 5-12

ELASTOPLASTIC ANALYSIS OF PLATE WITH BOUNDARY ELEMENT METHOD

10. The Discrete-Time Fourier Transform (DTFT)

Transcription:

CIVL / -D Boundr Vlu Prolms - Rctngulr Elmnts / RECANGULAR ELEMENS - In som pplictions, it m mor dsirl to us n lmntl rprsnttion of th domin tht hs four sids, ithr rctngulr or qudriltrl in shp. Considr four-nod rctngulr lmnt with prlll sids dfind s: (, ) (, ) (, ) 0 0 RECANGULAR ELEMENS - W will us ilinr distriution of th u ovr th lmnt, thrfor: u (, ) 0 0 RECANGULAR ELEMENS - Whr th function is oth linr in nd nd ilinr in. h four unknowns in th ov prssion m found fitting th curv through th four nodl vlus. u, i,,,or i i i i i i (, ) 0 0 RECANGULAR ELEMENS - Solving th ov four qutions for,,, nd, nd sustituting ths vlus ck into th lmntl distriution, th following st of shp functions m formd: N, N, h pproimtion for u m writtn in mtri form s:,, u un u Nu N i i N N,, RECANGULAR ELEMENS - h lmntl intrpoltions for rctngulr lmnt r: N N RECANGULAR ELEMENS - h shp functions r visull dciving. hr is no curvtur in dirctions prlll to n sid; howvr, thr is twist du to th trm in th lmnt rprsnttion. N N N N N N

CIVL / -D Boundr Vlu Prolms - Rctngulr Elmnts / RECANGULAR ELEMENS - On vr importnt dvntg of this tp of intrpoltion is tht th drivtivs r no longr constnt ovr n lmnt. h drivtivs in nd r: u u u u u u u u u u u u Rctngulr shp functions lso suffr from intrlmnt discontinuitis, in prticulr, th vlus of norml drivtivs t n dgs. RECANGULAR ELEMENS - h shp function m writtn in trms of th locl coordints nd : N, N, N, N, h shp functions m writtn in trms of nondimnsionl locl coordints s nd t, whr s = /, nd t = /: s t st N s, t N s, t s t st N s, t N s, t RECANGULAR ELEMENS - h dimnsionlss shp functions r dfind ovr unit squr - s, - t. h vrition of th glol coordints nd m prssd in trms of th shp functions s: N N i i i i RECANGULAR ELEMENS - h lmntl mtrics for th Poisson prolm r: k A N N NN da NN ds Lt s first considr how to hndl intgrtion ovr tpicl lmnt r A. For mpl: F, dd A f h A NfdA Nhds F, dd RECANGULAR ELEMENS - his intgrl m trnsformd into th locl coordint spc using = - 0 nd = - 0 : F, dd, F0, 0 dd, A whr (, )/(,) is th Jcoin of th trnsformtion. Using th trnsformtion = s nd = t th intgrl m writtn s: F, dd A, 0, 0 F s t dsdt st, A RECANGULAR ELEMENS - Evlution of k o convrt th intgrl into th dimnsionlss spc (s, t) th following trnsformtion is usd: s 0 Using th chin rul th drivtivs m writtn s: s t s s t t t 0 whr th Jcoin (,)/(s, t) hs vlu of A /.

CIVL / -D Boundr Vlu Prolms - Rctngulr Elmnts / RECANGULAR ELEMENS - Evlution of k In mtri form th ov prssion m writtn s: s s 0 J J s 0 t t h mtri form of th trnsformtion m invrtd. J J s s h mtri form of th trnsformtion m invrtd. J J s s whr J nd J r th first nd scond rows of J - RECANGULAR ELEMENS - Evlution of k hrfor, th prtil drivtivs of th shp functions m writtn s: N N N N J J J J s s Sustituting ll th pics of th trnsformtion in th k trms givs: A k JJ JJ ds dt JJ ds dt whr JJ = (J J + J J )A /. RECANGULAR ELEMENS - Evlution of k h vlus of th mtri nd JJ m computd s: t t t t 0 JJ s s s s 0 h rsulting lmntl stiffnss mtri k is tht is ddd to th glol sstm stiffnss mtri K t th corrsponding nodl loctions. k RECANGULAR ELEMENS - Evlution of k If th lmnt is squr, thn =, thn k coms: k k RECANGULAR ELEMENS - Evlution of f In gnrl, th intgrl is: f For gnrl function f(, ) w will ssum th f vris linrl ovr th lmnt, f(, ) = N f, whr th vctor f contins vlus of th function f t th nod points. With this ssumption th intgrl coms: A NfdA f NN fda da NN f A A RECANGULAR ELEMENS - Evlution of f rnsforming th intgrl into th non-dimnsionl coordints (s, t) ilds: A f NN ds dt f f f f f A f f f f f f f f f f f f f h rsulting lmntl lod vctor contriuts to th glol sstm qutions t thos loctions corrsponding to th four nods dfining th lmnt.

CIVL / -D Boundr Vlu Prolms - Rctngulr Elmnts / RECANGULAR ELEMENS - Evlution of h Considr th intgrl: h Nhds whr th intgrtion is long oundr sgmnt of th lmnt. Sinc, th intgrtion is computd long singl sid of th rctngulr lmnt, th originl shp functions rduc to: j k h k s 0 i l h l l s l h Nk hn k k hn l lds Nl s s Nk Nl l l RECANGULAR ELEMENS - Evlution of h Evluting th intgrls, th h trms rduc to: h l h k hl hk hl h rsulting lmntl lod vctor contriuts to th glol sstm qutions if th lmnt hs sid s prt of th oundr. RECANGULAR ELEMENS - Evlution of h vlution of th is vr similr to h cpt tht thr is n tr N in th intgrnd. h vrition of th function (s) will pproimtd s k N k + l N l. Considr th intgrl : NN ds l k l k l k l k l RECANGULAR ELEMENS Rcll, th glol sstm qutions r composd from th following summtions: K k F f h G G G G G G h rsulting sstm qutions r, in mtri form, givn s: Ku F G G G h rsulting stiffnss mtri contriuts to th glol sstm qutions whn th lmnt hs sid s prt of th oundr. RECANGULAR ELEMENS PROBLEM #9 - Vrif th componnts of th k, f,, nd h mtrics conform to th rsults dvlopd in th clss nots. In othr words, show ll th dtils of th drivtion of th stiffnss mtri nd th loding vctors. RECANGULAR ELEMENS Empl - Considr th non-dimnsionl torsion prolm: with, 0 in 0 on G h strsss nd torqu for th Prndlt strss function r: z G z G G dd

CIVL / -D Boundr Vlu Prolms - Rctngulr Elmnts / RECANGULAR ELEMENS Empl Discrtiztion - h simplst modl for torsion of squr r, utilizing smmtr is singl rctngulr lmnt. h gnrl prolm domin nd th FEM msh r shown low. Lins of Smmtr 0 (,) (0,) 0 0 n (0,0) (,0) cross-sctionl r 0 n RECANGULAR ELEMENS Empl - Elmntl Formultion - Using linr rctngulr lmnt th lmntl stiffnss mtri componnts r: k A loding function of f = givs lod vctor of: f f f f A f f f f f f f f f f f f f RECANGULAR ELEMENS Empl - Assml - Sinc thr is onl on lmnt in th modl th ssml is simpl: Constrints - For this modl, = 0 on th oundr, thrfor,, nd = 0. 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 RECANGULAR ELEMENS Empl - Solution - In this cs, th solution is quit simpl: G Computtion of Drivd Vrils - h prtil drivtivs with rspct to nd tht dfin th strss componnts r: G A G dd G.G G ct.g RECANGULAR ELEMENS Empl - Solution - In this cs, th solution is quit simpl: G (0,0) (0,) A (,0) (,) 0 (,) (0,) 0 0 n (0,0) (,0) 0 n RECANGULAR ELEMENS Empl Rpt th prvious non-dimnsionl torsion prolm using four rctngulr lmnts:, 0 in 0 on with G h strsss nd torqu for th Prndlt strss function r: z G z G G dd

CIVL / -D Boundr Vlu Prolms - Rctngulr Elmnts / RECANGULAR ELEMENS Empl Discrtiztion - Discrtiztion - h simplst modl for torsion of squr r, utilizing smmtr nd using four rctngulr lmnts is shown low. Lins of Smmtr 0 9 (0,) (,) 0 0 n (0,0) (,0) cross-sctionl r 0 n RECANGULAR ELEMENS Empl - Elmntl Formultion - Using linr rctngulr lmnt th lmntl stiffnss mtri componnts r: k For lmnt & : k k RECANGULAR ELEMENS Empl - Elmntl Formultion - Using linr rctngulr lmnt th lmntl stiffnss mtri componnts r: k For lmnt & : k 9 k 9 RECANGULAR ELEMENS Empl - Elmntl Formultion - A loding function of f = givs lod vctor of: f f f f A f f f f f f f f f f f f f f f f f 9 RECANGULAR ELEMENS Empl - Assml - Compiling ch lmntl mtri into th sstm qutions ilds: 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 9 RECANGULAR ELEMENS Empl - Constrints - For this modl, = 0 on th oundr, thrfor,,,,, nd 9 = 0. 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 9

CIVL / -D Boundr Vlu Prolms - Rctngulr Elmnts / RECANGULAR ELEMENS Empl - Solution - In this cs, th solution is quit simpl: 0 0 0 0 0 G G 0G 0G G 0 0 0 0 RECANGULAR ELEMENS Empl - Computtion of Drivd Vrils - h prtil drivtivs with rspct to nd tht dfin th strss componnts r: G dd ct,9g,0 G A.0G.G A RECANGULAR ELEMENS Empl - Computtion of Drivd Vrils - h prtil drivtivs with rspct to nd tht dfin th strss componnts r: (0,) 9 (,) (0,0) (,0) 0 n (0,) 0 9 (0,0) (,0) 0 n (,) 0 RECANGULAR ELEMENS PROBLEM #0 - Considr th non-dimnsionl torsion prolm:, 0 in 0 on G Lins of Smmtr cross-sctionl r 0 n (0,) 0 9 (,) (0,0) (,0) 0 n 0 End of Rctngulr Elmnts