TMMI37, vt2, Lecture 8; Introductory 2-dimensional elastostatics; cont.

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

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

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

u x v x dx u x v x v x u x dx d u x v x u x v x dx u x v x dx Integration by Parts Formula

Multiple-Choice Test Introduction to Partial Differential Equations COMPLETE SOLUTION SET

nd the particular orthogonal trajectory from the family of orthogonal trajectories passing through point (0; 1).

Differential Equations

u r du = ur+1 r + 1 du = ln u + C u sin u du = cos u + C cos u du = sin u + C sec u tan u du = sec u + C e u du = e u + C

AS 5850 Finite Element Analysis

a 1and x is any real number.

Higher order derivatives

Integration by Parts

Bifurcation Theory. , a stationary point, depends on the value of α. At certain values

Southern Taiwan University

dx equation it is called a second order differential equation.

(1) Then we could wave our hands over this and it would become:

VSMN30 FINITA ELEMENTMETODEN - DUGGA

Calculus Revision A2 Level

FEM FOR HEAT TRANSFER PROBLEMS دانشگاه صنعتي اصفهان- دانشكده مكانيك

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

INTEGRATION BY PARTS

Image Enhancement in the Spatial Domain

Direct Approach for Discrete Systems One-Dimensional Elements

6.1 Integration by Parts and Present Value. Copyright Cengage Learning. All rights reserved.

Derivation of Eigenvalue Matrix Equations

ELECTRON-MUON SCATTERING

Differential Equations

1 Isoparametric Concept

Finite element discretization of Laplace and Poisson equations


Hydrogen Atom and One Electron Ions

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

Differentiation of Exponential Functions

Dealing with quantitative data and problem solving life is a story problem! Attacking Quantitative Problems

Computer Vision. Fourier Analysis. Computer Science Tripos Part II. Dr Christopher Town. Fourier Analysis. Fourier Analysis

Massachusetts Institute of Technology Department of Mechanical Engineering

PHYSICS 489/1489 LECTURE 7: QUANTUM ELECTRODYNAMICS

Introduction to Condensed Matter Physics

Lecture 37 (Schrödinger Equation) Physics Spring 2018 Douglas Fields

Where k is either given or determined from the data and c is an arbitrary constant.

Calculus II (MAC )

Figure 1: Schematic of a fluid element used for deriving the energy equation.

For more important questions visit :

Addition of angular momentum

Introduction to Arithmetic Geometry Fall 2013 Lecture #20 11/14/2013

Thomas Whitham Sixth Form

Linear-Phase FIR Transfer Functions. Functions. Functions. Functions. Functions. Functions. Let

Numerical methods for PDEs FEM implementation: element stiffness matrix, isoparametric mapping, assembling global stiffness matrix

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

Pipe flow friction, small vs. big pipes

Unit 6: Solving Exponential Equations and More

A Propagating Wave Packet Group Velocity Dispersion

Math 34A. Final Review

First order differential equation Linear equation; Method of integrating factors

ME 321 Kinematics and Dynamics of Machines S. Lambert Winter 2002

SECTION where P (cos θ, sin θ) and Q(cos θ, sin θ) are polynomials in cos θ and sin θ, provided Q is never equal to zero.

NEW APPLICATIONS OF THE ABEL-LIOUVILLE FORMULA

Basic Polyhedral theory

CHAPTER 1. Introductory Concepts Elements of Vector Analysis Newton s Laws Units The basis of Newtonian Mechanics D Alembert s Principle

Partial Derivatives: Suppose that z = f(x, y) is a function of two variables.

Exam 1. It is important that you clearly show your work and mark the final answer clearly, closed book, closed notes, no calculator.

1 General boundary conditions in diffusion

2. Background Material

0WAVE PROPAGATION IN MATERIAL SPACE

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

Phys 344 Lect 17 March 16 th,

ph People Grade Level: basic Duration: minutes Setting: classroom or field site

Review of Exponentials and Logarithms - Classwork

Prelim Examination 2011 / 2012 (Assessing Units 1 & 2) MATHEMATICS. Advanced Higher Grade. Time allowed - 2 hours

Math 120 Answers for Homework 14

First derivative analysis

SAFE HANDS & IIT-ian's PACE EDT-15 (JEE) SOLUTIONS

Heisenberg Model. Sayed Mohammad Mahdi Sadrnezhaad. Supervisor: Prof. Abdollah Langari

MA 262, Spring 2018, Final exam Version 01 (Green)

Atomic Physics. Final Mon. May 12, 12:25-2:25, Ingraham B10 Get prepared for the Final!

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

8-node quadrilateral element. Numerical integration

Section 11.6: Directional Derivatives and the Gradient Vector

ME311 Machine Design

Addition of angular momentum

MCE503: Modeling and Simulation of Mechatronic Systems Discussion on Bond Graph Sign Conventions for Electrical Systems

Linear Non-Gaussian Structural Equation Models

Thomas Whitham Sixth Form

Probability and Stochastic Processes: A Friendly Introduction for Electrical and Computer Engineers Roy D. Yates and David J.

Chapter 8: Electron Configurations and Periodicity

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

This test is for two independent Populations. The test is sometimes called the Mann-Whitney U test or the Rank Sum Wilcoxon. They are equivalent.

Lecture 4: Parsing. Administrivia

y cos x = cos xdx = sin x + c y = tan x + c sec x But, y = 1 when x = 0 giving c = 1. y = tan x + sec x (A1) (C4) OR y cos x = sin x + 1 [8]

Introduction to the Fourier transform. Computer Vision & Digital Image Processing. The Fourier transform (continued) The Fourier transform (continued)

( ) Differential Equations. Unit-7. Exact Differential Equations: M d x + N d y = 0. Verify the condition

Derivation of Electron-Electron Interaction Terms in the Multi-Electron Hamiltonian

Potential energy of a structure. Vforce. joints j

5.80 Small-Molecule Spectroscopy and Dynamics

Slide 1. Slide 2. Slide 3 DIGITAL SIGNAL PROCESSING CLASSIFICATION OF SIGNALS

UNTYPED LAMBDA CALCULUS (II)

High Energy Physics. Lecture 5 The Passage of Particles through Matter

Einstein Equations for Tetrad Fields

Sara Godoy del Olmo Calculation of contaminated soil volumes : Geostatistics applied to a hydrocarbons spill Lac Megantic Case

y = 2xe x + x 2 e x at (0, 3). solution: Since y is implicitly related to x we have to use implicit differentiation: 3 6y = 0 y = 1 2 x ln(b) ln(b)

Transcription:

Lctr 8; ntrodctor 2-dimnsional lastostatics; cont. (modifid 23--3)

ntrodctor 2-dimnsional lastostatics; cont. W will now contin or std of 2-dim. lastostatics, and focs on a somwhat mor adancd lmnt thn th CST-lmnt, naml th -nodd bilinar (and isoparamtric) lmnt. (, ) (, 3 3 ) n (, ) n 2 (, 2 2 ) 2 Howr, lt s first start with a qick riw of th CST-lmnt! 2 (modifid 23--3)

Riw- th CST-lmnt (, 3 3 ) W do not s an lmntwis/local coordinat sstm; onl th global on! n (, ) n 2 ( 2, 2 ) 2 Contrclockwis nmbring of th nods 3 (modifid 23--3)

Th CST-lmnt; cont. On th prios lctr w fond th following transformation diagram T V f C T F E d k K C D n n 2 2 (modifid 23--3)

Th CST-lmnt; cont. E 2 ( ) / 2 T V f C T F E d k K C D n n 2 2 5 (modifid 23--3)

Th CST-lmnt; cont. T V f C T F E d k K C D n n 2 2 6 (modifid 23--3)

Th CST-lmnt; cont. 2 T V 2 f 3 C 3 T 2 2 3 3 F Th rst of th diagram as sal! ow, what do w gt for or - nodd lmnt? E d k K C D n n 2 2 7 (modifid 23--3)

Th bilinar -nodd lmnt; cont. As for CST-lmnt! E T _ tda f C d T k K C F D n n 2 2 8 (modifid 23--3)

Th bilinar -nodd lmnt; cont. Sinc th strss/strain filds will not b constant in a -nodd lmnt, th strain nrg mst b writtn as an intgral. From this it follows (b MPE) that w gt an intgral in th transformation diagram! E T _ tda f d C T k K C F D n n 2 2 9 (modifid 23--3)

Th bilinar -nodd lmnt; cont. E 2 T _ tda 2 3 f d C 3 T k K C F D 2 2 3 3 n ot as, as in th CSTcas, to find th shap fnctions to! n 2 2 (modifid 23--3)

Th bilinar -nodd lmnt; cont. Th basic mthod to st p th shap fnctions for a -nodd lmnt is to introdc an additional lmntwis/local (so calld natral) coordinat sstm. Howr, bfor that, w will look a littl bit at th lowr lft part of th transformation diagram! E T _ tda f d C T k K C F D n n 2 2 (modifid 23--3)

Th bilinar -nodd lmnt; cont. Mor spcificall, lt s split th matri [ ε] in two parts (which is not against th law :) d n n 2 2 2 (modifid 23--3)

(modifid 23--3) 3 d W ha i.. Th bilinar -nodd lmnt; cont.

Th bilinar -nodd lmnt; cont. Th lowr lft part of th transformation diagram ma ths b rwrittn as shown blow. And wh is that fin, on ma ask? Wll, as w soon will s, it opns th wa to mak an adantag of an additionall introdcd natral coordinat sstm! Priosl w had [ ε] hr! d n n 2 2 (modifid 23--3)

Th bilinar -nodd lmnt; cont. Lt s now introdc th natral coordinat sstm according to th illstration, whr th smbol n in th transformation diagram indicats that an ntit/qantit dpnds on th natral coordinats. T n n n n n d n n 2 2 n n 2 2 5 (modifid 23--3)

Th bilinar -nodd lmnt; cont. ilinar shap fnctions ar as to st p in th natral coordinats! n n 2 T n n n n n d n n 2 2 6 (modifid 23--3)

Th bilinar -nodd lmnt; cont. 2 2 n n 2 T n n n n n d n n 2 2 7 (modifid 23--3)

Th bilinar -nodd lmnt; cont. 3 3 n n 2 T n n n n n d n n 2 2 8 (modifid 23--3)

Th bilinar -nodd lmnt; cont. n n 2 T n n n n n d n n 2 2 9 (modifid 23--3)

Th bilinar -nodd lmnt; cont. n n 2 As can b sn, w will for simplicit skip th ind n whn writing th natral shap fnctions! Howr, thr is no risk of confsing things, sinc w onl will st p th shap fnctions in th natral coordinats. n n 2 2 2 (modifid 23--3)

(modifid 23--3) 2 2 3 (,) (,) 2 (,) 3 (,) 3 3 2 2 3 2 3 2 3 2 3 2 (,) (,) W not Th bilinar -nodd lmnt; cont.

Th bilinar -nodd lmnt; cont. T n n n n n Plas obsr! Ths displacmnt componnts dpnd on th natral coordinats! d n n 2 2 22 (modifid 23--3)

(modifid 23--3) 23 d n n n n T n Th sb-matri fond in th transformation matri [Tn] is calld th Jacoban matri of th coordinat chang, and will b dnotd [J]! Th chain rl implis [Tn] Th bilinar -nodd lmnt; cont.

Th bilinar -nodd lmnt; cont. Ths, th strain-displacmnt matri [] is gin b th following prssion T n n n Sinc w alrad ar in fll control of [], [ n] and [n], it jst rmains to find th inrs of th transformation matri [Tn]! W ha T n n n n n d J Tn Tn J J 22 J J J 2 J 2 J J J J 22 2 J J 2 2 (modifid 23--3)

Th bilinar -nodd lmnt; cont. Ths, w mst fi [J]! T Tn T n n J n J n n n n n f w choos to intrpolat th coordinats in th sam wa as th displacmnts (so calld isoparamtr), w ha d n : nod and can find [J] b diffrntiating th prssion for {} in an appropriat wa- ALL DOE! 25 (modifid 23--3)

Th bilinar -nodd lmnt; cont. n ordr to calclat th lmnt stiffnss, w mst tak car of th intgral k T E tda This intgration (in th natral coordinats) is don nmricall b so calld Gass-qadratr, which basicall mans that w alat th intgrand at a nmbr of points, and mltipl with an associatd ara. Dpnding on lmnt tp, a crtain nmbr of alation points ar ndd, which for or -nodd lmnt i. t mst finall b notd that also th abo lmnt somtims bhas poorl, and that a highr ordr lmnt thn mst b sd. 26 (modifid 23--3)