Static/Dynamic Deformation with Finite Element Method. Graphics & Media Lab Seoul National University

Similar documents
Fundamentals of Continuum Mechanics. Seoul National University Graphics & Media Lab

From Structural Analysis to FEM. Dhiman Basu

Variational Approach in FEM Part II

8-node quadrilateral element. Numerical integration

BAR & TRUSS FINITE ELEMENT. Direct Stiffness Method

The Hyperelastic material is examined in this section.

Lecture 23 APPLICATIONS OF FINITE ELEMENT METHOD TO SCALAR TRANSPORT PROBLEMS

AS 5850 Finite Element Analysis

Stress-Based Finite Element Methods for Dynamics Analysis of Euler-Bernoulli Beams with Various Boundary Conditions

Math 656 March 10, 2011 Midterm Examination Solutions

Lecture 08 Multiple View Geometry 2. Prof. Dr. Davide Scaramuzza

EE750 Advanced Engineering Electromagnetics Lecture 17

From Structural Analysis to Finite Element Method

PREDICTION OF STRESS CONCENTRATION FACTORS IN UNLAPPED SQUARE HOLLOW "K" JOINTS BY THE FINITE ELEMENT METHOD

Direct Approach for Discrete Systems One-Dimensional Elements

14. MODELING OF THIN-WALLED SHELLS AND PLATES. INTRODUCTION TO THE THEORY OF SHELL FINITE ELEMENT MODELS

CHAPTER 7d. DIFFERENTIATION AND INTEGRATION

Proceedings of the 2nd WSEAS Int. Conference on Applied and Theoretical Mechanics, Venice, Italy, November 20-22,

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

The Finite Element Method. Jerzy Podgórski

Phys 774: Nonlinear Spectroscopy: SHG and Raman Scattering

Grand Canonical Ensemble

Numerical solutions of fuzzy partial differential equations and its applications in computational mechanics. Andrzej Pownuk1

Discrete Shells Simulation

FINITE ELEMENT METHOD II Autumn 2015

:2;$-$(01*%<*=,-./-*=0;"%/;"-*

Linear Algebra Provides a Basis for Elasticity without Stress or Strain

Phy213: General Physics III 4/10/2008 Chapter 22 Worksheet 1. d = 0.1 m

Hydrogen Atom and One Electron Ions

Finite element discretization of Laplace and Poisson equations

3.4 Properties of the Stress Tensor

Finite Element Based Implementation of Fiala s Thermal Manikin in THESEUS-FE

External Equivalent. EE 521 Analysis of Power Systems. Chen-Ching Liu, Boeing Distinguished Professor Washington State University

Polytropic Process. A polytropic process is a quasiequilibrium process described by

Journal of Chemical and Pharmaceutical Research, 2014, 6(5): Research Article

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

Introduction to elastic wave equation. Salam Alnabulsi University of Calgary Department of Mathematics and Statistics October 15,2012

Work is the change in energy of a system (neglecting heat transfer). To examine what could

ANALYSIS: The mass rate balance for the one-inlet, one-exit control volume at steady state is

SERBIATRIB th International Conference on Tribology. Belgrade, Serbia, May 2015

A Note on Estimability in Linear Models

HORIZONTAL IMPEDANCE FUNCTION OF SINGLE PILE IN SOIL LAYER WITH VARIABLE PROPERTIES

IV. First Law of Thermodynamics. Cooler. IV. First Law of Thermodynamics

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

PHYS 1443 Section 004 Lecture #12 Thursday, Oct. 2, 2014

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

Three-Node Euler-Bernoulli Beam Element Based on Positional FEM

Finite Element Models for Steady Flows of Viscous Incompressible Fluids

Jones vector & matrices

MECH321 Dynamics of Engineering System Week 4 (Chapter 6)

First looking at the scalar potential term, suppose that the displacement is given by u = φ. If one can find a scalar φ such that u = φ. u x.

FEFF and Related Codes

Modeling of the Through-the-Thickness Electric Potentials of a Piezoelectric Bimorph Using the Spectral Element Method

September 27, Introduction to Ordinary Differential Equations. ME 501A Seminar in Engineering Analysis Page 1. Outline

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

Potential energy of a structure. Vforce. joints j

PHYS 1441 Section 002 Lecture #15

Division of Mechanics Lund University MULTIBODY DYNAMICS. Examination Name (write in block letters):.

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

PHYS 1443 Section 002 Lecture #20

The gravitational field energy density for symmetrical and asymmetrical systems

Electrochemical Equilibrium Electromotive Force. Relation between chemical and electric driving forces

Nonlinear Behavior of Pile-Soil Subjected to Torsion due to Environmental Loads on Jacket Type Platforms

STABILIZED FINITE ELEMENTS IN GEOMECHANICAL APPLICATIONS

horizontal force output data block Hankel matrix transfer function complex frequency response function impedance matrix

The van der Waals interaction 1 D. E. Soper 2 University of Oregon 20 April 2012

Physics 181. Particle Systems

Finite Element Method for Turbomachinery Flows

1 Introduction. M.P. Santisi d Avila 1 and F. Lopez Caballero 2. Abstract

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

CHAPTER 8 Potential Energy and Conservation of Energy

Section 11.6: Directional Derivatives and the Gradient Vector

Heat/Di usion Equation. 2 = 0 k constant w(x; 0) = '(x) initial condition. ( w2 2 ) t (kww x ) x + k(w x ) 2 dx. (w x ) 2 dx 0.

Einstein Equations for Tetrad Fields

Humanistic, and Particularly Classical, Studies as a Preparation for the Law

TRANSIENT RESPONSE OF A THREE LAYER SANDWICH PLATE WITH ELECTRORHEOLOGICAL CORE AND ORTHOTROPIC FACES

(A) the function is an eigenfunction with eigenvalue Physical Chemistry (I) First Quiz

ACOUSTIC WAVE EQUATION. Contents INTRODUCTION BULK MODULUS AND LAMÉ S PARAMETERS

System variables. Basic Modeling Concepts. Basic elements single and. Power = effort x flow. Power = F x v. Power = V x i. Power = T x w.

Distortional Analysis of Thin-Walled Box Girder Structure: A Comparative Study

COMPUTATIONAL MECHANICS (2º cycle, Master of Science Degree)

Physics 256: Lecture 2. Physics

Conservation of Energy

Problem 22: Journey to the Center of the Earth

C PLANE ELASTICITY PROBLEM FORMULATIONS

LARGE DISPLACEMENT ANALYSIS OF SLENDER ARCHES

Dynamic Modeling and Vibration Control for Spacecraft s Solar Array Jian-Ping JIANG 1,a, *, Rui XU 2,b

1. Basic Concepts. 1.1 What is Kinematics 1.2 Solution Methods 1.3 Choice of Coordinates 1.4 2D versus 3D. Professor Wan-Suk Yoo

Review - Probabilistic Classification

,. *â â > V>V. â ND * 828.

A Probabilistic Characterization of Simulation Model Uncertainties

Constitutive behaviour of masonry

Stretching and bending deformations due to normal and shear tractions of doubly curved shells using third-order shear and normal deformable theory

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

Ερωτήσεις και ασκησεις Κεφ. 10 (για μόρια) ΠΑΡΑΔΟΣΗ 29/11/2016. (d)

(most) due to long range e m forces i.e. via atomic collisions or due to short range nuclear collisions or through decay ( = weak interactions)

r/lt.i Ml s." ifcr ' W ATI II. The fnncrnl.icniccs of Mr*. John We mil uppn our tcpiiblicnn rcprc Died.

H A M M IG K S L IM IT E D, ' i. - I f

.-I;-;;, '.-irc'afr?*. P ublic Notices. TiffiATRE, H. aiety

W i n t e r r e m e m b e r t h e W O O L L E N S. W rite to the M anageress RIDGE LAUNDRY, ST. H E LE N S. A uction Sale.

Transcription:

Statc/Dynamc Dormaton wth Fnt Elmnt Mthod Graphcs & Mda Lab Sol Natonal Unvrsty

Statc/Dynamc Dormaton Statc dormaton Dynamc dormaton ndormd shap ntrnal + = nrta = trnal dormd shap statc qlbrm dynamc qlbrm

How to Smlat Statc Dormaton?

Govrnng Eqaton W s th dsplacmnt ld Contnos D vctor ncton a : R R a At vry matral pont, nd sch that lastc potntal nrgy U t trnal orc Lagrang Eqaton

Ovrvw o Statc Dormaton U t. Drv th lastc potntal nrgy n trms o stran and strss. Solv th partal drntal qaton wth nt lmnt dscrtzaton

For Now, Us Cachy s n. Stran nsor Cachy Lt s dn a colmn vctor... ε Cachy Cachy ε

Hook s Law n a Matr Form z yz y zz yy z yz y zz yy z z yz yz y y G G G yy zz zz zz yy yy zz yy E E E

Hook s Law n a Matr Form ε E σ z yz y zz yy z yz y zz yy G G G E E E E E E E E E Not that th stran nrgy dnsty U can now b dnd as ε σ U yy zz zz zz yy yy zz yy E E E z z yz yz y y G G G

Now, th Elastc Potntal Enrgy U V U dv dv ε σ ε E ε V V dv yy zz y yz z yy zz y yz z Now, th govrnng qaton: U t ε V Eε dv t

Ovrvw o Statc Dormaton U t. Drv th lastc potntal nrgy n trms o stran and strss. Solv th partal drntal qatons wth nt lmnt dscrtzaton

Ovrvw o Statc Dormaton U t. Drv th lastc potntal nrgy n trms o stran and strss. Solv th partal drntal qatons wth nt lmnt dscrtzaton

Fnt Elmnt Dscrtzaton h abov govrnng qaton s or th contnm. t holds or ach matral pont. For smlaton, w nd to prorm nt lmnt dscrtzaton. W s lnar tntrahdral lmnts ndormd dormd t V Eε dv ε n n n n

Eprssng n trms o Nodal Ponts Eprss an arbtrary pont n trms o,,, = b + b + b + b, wth b + b + b + b = hy ar calld th barycntrc coordnats b b b b b X sx X sx A HW: Drv ths X s ndormd dormd

Eprssng n trms o b b b b vctor o nodal dsplacmnts b b b b H matr H th lnar shap ncton H b b b b b b b b b b b b ndormd dormd

HW: Calclat th matr B dpnds on th orgnal ttrahdron only! Now, th Stran nsor 6 B R H B H ε z yz y zz yy

Strss and Potntal Enrgy Strss as a ncton o th dsplacmnts σ Eε EB ε B Enrgy as a ncton o th dsplacmnts U V ε σ dv V B EB dv V B EB V B EB dv

Stnss Matr and Elastc Forc Potntal nrgy U V B EB K stnss matr Forc d to lastc potntal nrgy U K K HW: Drv ths Wth ttrahdral lmnts and Cachy stran mtrc h stnss matr dpnds only on th rst stat msh h lastc orc s lnar wth th nodal dsplacmnt Do th abov hold whn Grn s sd or th stran mtrc? HW: stdy ths and rport n th nt class

Statc Eqlbrm ntrnal orc rsltng rom lastc potntal nrgy o th lmnt shold balanc wth smmaton o all th trnal orcs actng on that lmnt K trnal orcs Gravtatonal orc actng on an lmnt Dormd or ndormd, an lmnt s ndr gravty. So RHS o th qlbrm qaton shold ncld th trms contrbtd rom gravty. How mch dos a dormd lmnts gt th gravtatonal potntal?

Gravtatonal Forc Frst, th gravtatonal potntal nrgy o a dormd lmnt Now, th gravtatonal orc dv G V g g H dv G V g H 4 V V dv HW: Drv ths HW: Drv ths g H g dv dv V V

Wth Cachy Stran nsor For an lmnt U t U K

Assmblng Elmnt Stnss Matrcs For an lmnt U t U K ndormd dormd

Assmblng Elmnt Stnss Matrcs For a sngl lmnt... For th ntr body... K K n n K 9 7 4 9 7 4 K n n K j K j

Assmblng Elmnt Stnss Matrcs For a sngl lmnt... For th ntr body... K K n n K 9 7 4 9 7 4 K n n K j K j Chck Ponts What dos K j mans? How many non-zro blocks ar thr?

How to Smlat Dynamc Dormaton?

Dynamc Dormaton Govrnng qatons: Elr-Lagrang qatons kntc nrgy lastc potntal nrgy d dt U t trnal orc, t, t ndormd dormd

Kntc Enrgy h Kntc Enrgy V dv t dt d t dt d,,, t, t ndormd dormd M V dv H H th mass matr HW: Drv ths

Forc d to Kntc Enrgy h kntc nrgy Forc d to kntc nrgy dt d dt d M M V dv t dt d t dt d,, M V dv H H th mass matr HW: Drv ths M V HW: Drv ths

What s th nttv Manng? Govrnng qatons: Elr-Lagrang qatons t U dt d K M

Dampng Matr Govrnng qatons wth dampng trms d dt U C t M C K dampng matr Adopt proportonal Raylgh dampng: Or, ncorporat vsco-lastc ormlaton C M K

Assmblng Eqs or All th Elmnts Sngl lmnt Entr body K C M K C M n n C K M,, n n K C M,, 9 7 4 9 7 4 K C M,,

Solvng th Lnar Systm M, C, K matrcs ar tm ndpndnt! Convrt t to a lnar systm Eplct ntgraton mplct ntgraton K C M h h K Cv M v v v K C M v v h dt d h h h h h Kv K Cv K C M v v v HW: Drv ths

Statc/Dynamc Dormaton wth Fnt Elmnt Mthod HE END