CIVL 8/ D Boundary Value Problems - Quadrilateral Elements (Q4) 1/8

Similar documents
Elementary Differential Equations and Boundary Value Problems

DEPARTMENT OF ELECTRICAL &ELECTRONICS ENGINEERING SIGNALS AND SYSTEMS. Assoc. Prof. Dr. Burak Kelleci. Spring 2018

2. The Laplace Transform

CIVL 8/ D Boundary Value Problems - Triangular Elements (T6) 1/8

Boyce/DiPrima 9 th ed, Ch 2.1: Linear Equations; Method of Integrating Factors

Transfer function and the Laplace transformation

Chapter 12 Introduction To The Laplace Transform

Midterm exam 2, April 7, 2009 (solutions)

Boyce/DiPrima 9 th ed, Ch 7.8: Repeated Eigenvalues

Wave Equation (2 Week)

Chap.3 Laplace Transform

S.Y. B.Sc. (IT) : Sem. III. Applied Mathematics. Q.1 Attempt the following (any THREE) [15]

UNIT #5 EXPONENTIAL AND LOGARITHMIC FUNCTIONS

Lecture 26: Leapers and Creepers

LaPlace Transform in Circuit Analysis

EXERCISE - 01 CHECK YOUR GRASP

Chapter 5 The Laplace Transform. x(t) input y(t) output Dynamic System

Instructors Solution for Assignment 3 Chapter 3: Time Domain Analysis of LTIC Systems

4.i. 4 Numerical Model Development 4.1

Lecture 1: Numerical Integration The Trapezoidal and Simpson s Rule

REPETITION before the exam PART 2, Transform Methods. Laplace transforms: τ dτ. L1. Derive the formulas : L2. Find the Laplace transform F(s) if.

CHAPTER. Linear Systems of Differential Equations. 6.1 Theory of Linear DE Systems. ! Nullcline Sketching. Equilibrium (unstable) at (0, 0)

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

Copyright 2012 Pearson Education, Inc. Publishing as Prentice Hall.

Logistic equation of Human population growth (generalization to the case of reactive environment).

Final Exam : Solutions

PWM-Scheme and Current ripple of Switching Power Amplifiers

Spring 2006 Process Dynamics, Operations, and Control Lesson 2: Mathematics Review

1 Finite Automata and Regular Expressions

Let s look again at the first order linear differential equation we are attempting to solve, in its standard form:

AR(1) Process. The first-order autoregressive process, AR(1) is. where e t is WN(0, σ 2 )

Transmission Line Theory

where: u: input y: output x: state vector A, B, C, D are const matrices

The transition:transversion rate ratio vs. the T-ratio.

14.02 Principles of Macroeconomics Fall 2005 Quiz 3 Solutions

SOLUTIONS. 1. Consider two continuous random variables X and Y with joint p.d.f. f ( x, y ) = = = 15. Stepanov Dalpiaz

fiziks Institute for NET/JRF, GATE, IIT JAM, JEST, TIFR and GRE in PHYSICAL SCIENCES MATEMATICAL PHYSICS SOLUTIONS are

CSE 245: Computer Aided Circuit Simulation and Verification

DE Dr. M. Sakalli

On the Speed of Heat Wave. Mihály Makai

Applied Statistics and Probability for Engineers, 6 th edition October 17, 2016

Lecture 2: Current in RC circuit D.K.Pandey

Lecture 4: Laplace Transforms

Ma/CS 6a Class 15: Flows and Bipartite Graphs

15th European Signal Processing Conference (EUSIPCO 2007), Poznan, Poland, September 3-7, 2007, copyright by EURASIP

NARAYANA I I T / P M T A C A D E M Y. C o m m o n P r a c t i c e T e s t 1 6 XII STD BATCHES [CF] Date: PHYSIS HEMISTRY MTHEMTIS

Math 266, Practice Midterm Exam 2

On the Derivatives of Bessel and Modified Bessel Functions with Respect to the Order and the Argument

whereby we can express the phase by any one of the formulas cos ( 3 whereby we can express the phase by any one of the formulas

3(8 ) (8 x x ) 3x x (8 )

2.1. Differential Equations and Solutions #3, 4, 17, 20, 24, 35

A THREE COMPARTMENT MATHEMATICAL MODEL OF LIVER

Institute of Actuaries of India

XV Exponential and Logarithmic Functions

H is equal to the surface current J S

Chapter 3: Fourier Representation of Signals and LTI Systems. Chih-Wei Liu

Poisson process Markov process

Midterm Examination (100 pts)

[ ] 1+ lim G( s) 1+ s + s G s s G s Kacc SYSTEM PERFORMANCE. Since. Lecture 10: Steady-state Errors. Steady-state Errors. Then

4.1 The Uniform Distribution Def n: A c.r.v. X has a continuous uniform distribution on [a, b] when its pdf is = 1 a x b

SME 3033 FINITE ELEMENT METHOD. Bending of Prismatic Beams (Initial notes designed by Dr. Nazri Kamsah)

Decline Curves. Exponential decline (constant fractional decline) Harmonic decline, and Hyperbolic decline.

3.4 Repeated Roots; Reduction of Order

Heat flow in composite rods an old problem reconsidered

( ) ( ) + = ( ) + ( )

Module 1-2: LTI Systems. Prof. Ali M. Niknejad

Review Lecture 5. The source-free R-C/R-L circuit Step response of an RC/RL circuit. The time constant = RC The final capacitor voltage v( )

Control System Engineering (EE301T) Assignment: 2

CPSC 211 Data Structures & Implementations (c) Texas A&M University [ 259] B-Trees

1. Inverse Matrix 4[(3 7) (02)] 1[(0 7) (3 2)] Recall that the inverse of A is equal to:

Lecture 1: Growth and decay of current in RL circuit. Growth of current in LR Circuit. D.K.Pandey

THE INVOLUTE-EVOLUTE OFFSETS OF RULED SURFACES *

Fourier. Continuous time. Review. with period T, x t. Inverse Fourier F Transform. x t. Transform. j t

FIRST-ORDER SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS I: Introduction and Linear Systems

Laplace Transforms recap for ccts

Phys463.nb Conductivity. Another equivalent definition of the Fermi velocity is

Why Laplace transforms?

Math 3301 Homework Set 6 Solutions 10 Points. = +. The guess for the particular P ( ) ( ) ( ) ( ) ( ) ( ) ( ) cos 2 t : 4D= 2

Partial Fraction Expansion

Part I: Short Answer [50 points] For each of the following, give a short answer (2-3 sentences, or a formula). [5 points each]

Part 3 System Identification

Chapter 10. The singular integral Introducing S(n) and J(n)

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

Discussion 06 Solutions

Numerical Simulation for the 2-D Heat Equation with Derivative Boundary Conditions

y z P 3 P T P1 P 2. Werner Purgathofer. b a

Charging of capacitor through inductor and resistor

Lecture 4: Parsing. Administrivia

16.512, Rocket Propulsion Prof. Manuel Martinez-Sanchez Lecture 3: Ideal Nozzle Fluid Mechanics

Source code. where each α ij is a terminal or nonterminal symbol. We say that. α 1 α m 1 Bα m+1 α n α 1 α m 1 β 1 β p α m+1 α n

Revisiting what you have learned in Advanced Mathematical Analysis

EE 434 Lecture 22. Bipolar Device Models

with Dirichlet boundary conditions on the rectangle Ω = [0, 1] [0, 2]. Here,

ANALYSIS OF LAMINATED CONICAL SHELL STRUCTURES USING HIGHER ORDER MODELS. and J. Herskovits c

AN INTRODUCTION TO FOURIER ANALYSIS PROF. VEDAT TAVSANOĞLU

Fourier Series and Parseval s Relation Çağatay Candan Dec. 22, 2013

Voltage v(z) ~ E(z)D. We can actually get to this wave behavior by using circuit theory, w/o going into details of the EM fields!

5. An object moving along an x-coordinate axis with its scale measured in meters has a velocity of 6t

ERROR ANALYSIS A.J. Pintar and D. Caspary Department of Chemical Engineering Michigan Technological University Houghton, MI September, 2012

EAcos θ, where θ is the angle between the electric field and

Transcription:

CIVL 8/7111 -D Boundar Vau Prom - Quadriara Emn (Q) 1/8 ISOPARAMERIC ELEMENS h inar rianguar mn and h iinar rcanguar mn hav vra imporan diadvanag. 1. Boh mn ar una o accura rprn curvd oundari, and. h provid poor approimaion for an drivd varia. B uiizing highr-ordr mn ioparamric mn h dficinci of h inar mn ar ovrcom. ISOPARAMERIC ELEMENS h rm ioparamric man ha an qua numr of paramr ar ud o rprn h gomr and h dpndn varia. In ohr word, h am hap funcion ud o inrpoa h dpndn varia ar ud o rprn h gomr. A uparamric mn i on whr h gomr i dfind fwr paramr han ud o inrpoa h ouion. A uprparamric mn i on whr h gomr i dfind mor paramr han ud o inrpoa h ouion. ISOPARAMERIC ELEMENS h rm ioparamric man ha an qua numr of paramr ar ud o rprn h gomr and h dpndn varia. u un i i 6-nod rianguar mn (6) N N i i i i 8-nod quadriara mn (Q8) h four-nod quadriara mn i vr imiar o h four-nod rcanguar mn. Boh provid iinar inrpoaion ovr h mn for h gnra ouion and a inar approimaion of h fir drivaiv. Howvr, h quadriara mn i no rricd o rcanguar hap; in fac i ma rprn an numr of rciinar hap. h gnra approimaion ovr a quadriara mn i aumd o : u, un, un, un, un, u NuN,,,, h mna coordina ar and. h inrpoaion funcion N i (, ) ar idnica o h rcanguar mn inrpoaion funcion. N, N, 1 1 1 1 1 N, N, 1 1 1 1

CIVL 8/7111 -D Boundar Vau Prom - Quadriara Emn (Q) /8 h ranformaion from goa coordina (, ) o mna or parna coordina (, ) i accompihd uing h foowing coordina mapping:,,,,, N N N N,,,,, N N N N N N N N h inrpoaion funcion for a quadriara mn ar rfrrd a hap funcion inc h prion dfin h hap of h mn. L ook a how h in and poin in h parn mn (, ) ar mappd or ranformd o h goa m (, ).,, * (,) 1 1,, * (,) 1 1 1, 1, (,) 1 a 1 1 1, 1, (,) 1 a 1 1 In gnra raigh in in h - coordina m ar mappd ino raigh in in h quadriara mn. For amp, poin a and ar mappd o poin a* and * in h quadriara. h raigh in = 1 in h parn mn, conncing nod and, ma mappd ino h quadriara a:,, * (,) 1 1,, * (,) 1 1 1, 1, (,) 1 a 1 1 1, 1, (,) 1 a 1 1 In ordr o undrand h naur of h ranformaion, h chain ru i ud o form h diffrnia raionhip: In mari noaion h drivaiv ma wrin a: J whr J J h drminan of h Jacoian mari, J, i a of h invriii of h ranformaion = (, ) and = (, ). Whn J i poiiv vrwhr in a rgion, h ranformaion ma invrd o drmin = (,) and = (, ). hi man ha for a givn poin (, ) in h quadriara hr i a uniqu corrponding poin (, ) in h parn mn.

CIVL 8/7111 -D Boundar Vau Prom - Quadriara Emn (Q) /8 h mn of h Jacoian mari ma compud a: 11 1 11 1 1 11 1 11 hrfor, J ha h form: J j1 j j j whr h conan j i ar funcion of h coordina of h nod (funcion of h gomr of h mn). h conan j i acua zro o ha J i a inar funcion of and. h drminan J gomrica rprn h raionhip wn h diffrnia ara da(, ) = d d and h diffrnia ara da(, ). J da, da, J d d PROBLEM #1 - Show ha h conan j dfining h drminan of h Jacoian J i zro. Eamp - Conidr h quadriara mn givn ow. 1, 0, 10,10 0, 11 h mn of h Jacoian mari ma cacuad a: 1 1 1 1 18 Eamp - Conidr h quadriara mn givn ow. 1, 0, 10,10 0, 11 h mn of h Jacoian mari ma cacuad a: 1 1 1 19 1 1 1 1 Eamp - h drminan i: 10 8 10 9 178 188 170 J 8 Conidr h ocaion of nod hr a = 19 and =, hn h ranformaion i givn : 67 6 1 98 111 66 79 1 177 J 8

CIVL 8/7111 -D Boundar Vau Prom - Quadriara Emn (Q) /8 Eamp - Ovr h inrva -1 1 and -1 1 h drminan of h Jacoian, J, i poiiv, indicaing ha hr ar no prom wih h mapping from - o -. h ara of h quadriara ma compud a: 98 A A d d J d d 1 177 d d 8 66 79 177 Eamp - h parn mn in h - m and h hap of h quadriara mn ar hown in h figur ow: 1, 0 10,10 * 19, (,) (,) 0, 11 1 a 1 1 1 1 111 Eamp - Conidr h ocaion of nod hr a = 16 and = 1, hn h ranformaion i givn : 9 9 7 5 7 11 9 7 J 0 1 5 J Eamp - Ovr h inrva -1 1 and -1 1 h drminan of h Jacoian, J, i no awa poiiv. In fac, J i ngaiv for ( + ) 5/. For poin (, ) who vau aif ( + ) 5/, h mapping oca hm ouid of h quadriara. h prom i aociad wih h fac h rgion i conv aong h id dfind nod,, and. 9 7 J 0 8 1 8 1 Eamp - h parn mn in h - m and h hap of h quadriara mn ar hown in h figur ow: 1, 0 10,10 (,) 16, 1 (,) 0, 11 1 a 1 1 1 1 h mna maric for h Poion prom ar: NN NN k da f f da N A A a NN d h Nh d

CIVL 8/7111 -D Boundar Vau Prom - Quadriara Emn (Q) 5/8 h dvopmn of oh h k and f rm i idnica o ha prnd for rcanguar mn: k JJ d d f Nf J d d whr JJ = (J 1 J 1 + J J ) J. h vau of h mari ma compud a: Rca, J 1 and J ar h fir and cond row of h invr of h Jacoian mari. 1 1 J J1 J J J J h mn of h invr of h Jacoian ar: 11 1 11 1 1 11 1 11 J For a gnra quadriara mn, h mari JJ i a funcion of and wih h Jacoian J in h dnominaor. h ruing prion for JJ i vr difficu o vaua ac and h ingraion ar uua don numrica. hrfor, h rm k and f ma ca in h foowing form: G, dd whr G(,) i a compicad funcion of h varia and. In princip, i ma poi o vaua h f rm, howvr, numrica ingraion i pica mor pracica. For h k rm, h apparanc of h Jacoian J in h ingrand gnra indica h u of numrica quadraur. I i foruna ha h imi of h ingra conform o h forma ud in Gau-Lgndr quadraur. Sinc h ingraion ar in oh h and -dircion, on approach i o app Gauian quadraur para in ach dircion. hrfor, h gnra prion for k and f an N-rm Gauian quadraur i: k JJ d d ww JJ f Nf Jdd i n i n 1 j1 1 j1 whr h ucrip ij rfr o h vauaion of h quani a i and j. n n i j ij ij ij,,,, ww N f J i j i j i j i j ij Conidr h ingra a and h : a NN d h Nh d whr h ingraion i aong a oundar gmn of h mn. Sinc, h ingraion i compud aong a ing id of h quadriara mn, h quadriara hap funcion rduc o inar funcion.

CIVL 8/7111 -D Boundar Vau Prom - Quadriara Emn (Q) 6/8 h ru for hi ingraion ar idnica o prion oaind for inar rianguar and rcanguar mn. j a h k h k 0 i Nk N k k k Nk hn k k hn d N h N N N N d Nk 1 N Evauaing h ingra, h a and h rm rduc o: j a k h k 0 i k k 1 k k h Nk 1 N h h k h 6 hk h h ruing a mna iffn mari and h 1 h mna oad vcor conriu o h goa m quaion if h mn ha a id a par of h oundar. h goa m quaion ar compod from h foowing ummaion: K k a F f h G G G G G G h ruing m quaion ar, in mari form, givn a: Ku F G G G Eamp - Conidr h quadriara mn givn ow. 1, 0, 10,10 0, 11 Cacua h vau of h iffn mari k for h aov mn whn = 19 and = 1. k JJ d d f Nf J d d 1 1 Eamp - Rca JJ = (J 1 J 1 + J J ) J. h vau of h mari ma compud a: whr J 1 and J ar h fir and cond row of h invr of h Jacoian mari. 1 1 1 J J J 1 J J1 J Eamp - h mn of h invr of h Jacoian ar: 1 1 1 1 1 1 1 1 1 1 1 1

CIVL 8/7111 -D Boundar Vau Prom - Quadriara Emn (Q) 7/8 Eamp - h cacuaion ar prformd compuing h vau of ach componn of,, J 1, J, and J maric a ach Gau poin, prforming h mari muipicaion, and fina umming h ru for a h Gau poin. h u of numrica quadraur a hi poin in h formuaion ma gnra om rror. Howvr, if nough Gau poin ar ud h numrica ru ar qui accura whn compard o h ac vau for h ingraion. Eamp - A wo-dimniona quadraur uing N = i acua Gau poin and uing N = i 16 quadraur poin. 1 N = 1 Gau poin 1 N = 1 N = Eamp - hi fac i vidn in h foowing cacuaion uing Gauian quadraur wih N = 1,,, and. N 1 0.09605 0.101695 0.09605 0.101695 0.6559 0.101695 0.6559 k 0.09605 0.101695 mmric 0.6559 Eac k 0.5519 0.0106 0.781 0.0975 0.7519 0.05195 0.6917 0.6771 0.878 mmric 0.91507 Eamp - hi fac i vidn in h foowing cacuaion uing Gauian quadraur wih N = 1,,, and. N Eac k k 0.55081 0.1989 0.86 0.0966 0.790 0.006 0.705 0.669 0.71 mmric 0.9171 0.5519 0.0106 0.781 0.0975 0.7519 0.05195 0.6917 0.6771 0.878 mmric 0.91507 Eamp - hi fac i vidn in h foowing cacuaion uing Gauian quadraur wih N = 1,,, and. N Eac k k 0.551 0.0101 0.790 0.097 0.7515 0.05188 0.695 0.6770 0.876 mmric 0.9159 0.5519 0.0106 0.781 0.0975 0.7519 0.05195 0.6917 0.6771 0.878 mmric 0.91507 Eamp - hi fac i vidn in h foowing cacuaion uing Gauian quadraur wih N = 1,,, and. N Eac k k 0.5519 0.0106 0.781 0.0975 0.7519 0.05195 0.6917 0.6771 0.878 mmric 0.91507 0.5519 0.0106 0.781 0.0975 0.7519 0.05195 0.6917 0.6771 0.878 mmric 0.91507

CIVL 8/7111 -D Boundar Vau Prom - Quadriara Emn (Q) 8/8 Eamp - h ru for N = 1 ar no vr accura. Sinc on on Gau poin i ud in h quadraur, hi ru i no unpcd. h vau for k uing Gauian quadraur wih N = and indica ha h accurac of h vauaion incra a h numr of Gau poin incra. Eamp - Drivd Varia - h drivd varia for h Poion or Lapac quaion ar h paria drivaiv u/ and u/. h rm ar compud a funcion of poiion in an mn uing h hap funcion and h noda vau. h ru uing N = ar ac o fiv or i dcima pac. Eamp - Drivd Varia - h paria drivaiv ma comind o giv h norma or dirciona drivaiv: nin jn u u n un u u n n u u J u u J1 u 1 u J u u J Eamp - Drivd Varia - Rca J 1 and J ar h fir and cond row of J -1, rpciv. Rca ha u / i: u [ Nu] u u u [ Nu] hrfor: u u J1 u J u h dirciona or norma drivaiv i: u nj n J u n 1 End of Quadriara mn (Q)