o metod kone n ch prvk lect 01 ppt n.
Download
Skip this Video
Loading SlideShow in 5 Seconds..
O metod ě konečných prvků Lect_01 PowerPoint Presentation
Download Presentation
O metod ě konečných prvků Lect_01

Loading in 2 Seconds...

play fullscreen
1 / 47

O metod ě konečných prvků Lect_01 - PowerPoint PPT Presentation


  • 101 Views
  • Uploaded on

O metod ě konečných prvků Lect_01.ppt. Syllabus and introduction. M. Okrouhlík Ústav termomechaniky, AV ČR , Praha Plze ň , 2010. Syllabus. Deformační varianta MKP Značení Odvození pomocí principu virtuálních prací Diskretizace posuvů, přetvoření a konstitutivních vztahů

loader
I am the owner, or an agent authorized to act on behalf of the owner, of the copyrighted work described.
capcha
Download Presentation

PowerPoint Slideshow about 'O metod ě konečných prvků Lect_01' - julius


An Image/Link below is provided (as is) to download presentation

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.


- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -
Presentation Transcript
o metod kone n ch prvk lect 01 ppt

O metodě konečných prvkůLect_01.ppt

Syllabus and introduction

M. Okrouhlík

Ústav termomechaniky, AV ČR, Praha

Plzeň, 2010

syllabus
Syllabus

Deformační varianta MKP

Značení

Odvození pomocí principu virtuálních prací

Diskretizace posuvů, přetvoření a konstitutivních vztahů

Strukturální prvky – tyč, nosník, membrána, deska, skořepina

Analytický přístup – zobecněné souřadnice

Numerický přístup – isoparametrické prvky

Sestavení matic tuhosti, tlumení a hmotnosti

Předepsání okrajových podmínek

Typy řešených úloh

Řešení statických úloh

Nalezení vlastních čísel a vlastních tvarů kmitu

Řešení nestacionárních úloh – šíření vln

Numerická matematika

Řešení soustav algebraických rovnic

Řešení standardního a zobecněného problému vlastních čísel

Integrace obyčejných diferenciálních rovnic

Metoda konečných prvků pro nelineární úlohy – stručný úvod do problematiky

slide4

http://www.colorado.edu/engineering/CAS/courses.d/IFEM.d/Home.htmlhttp://www.colorado.edu/engineering/CAS/courses.d/IFEM.d/Home.html

jak se nem ps t diserta n pr ce
Jak se nemá psát disertační práce

Při výčtu otců zakladatelů metody konečných prvků autor jedné disertační prácevěnované MKP (práce byla obhajována v roce 2010 a nebyla z Plzně) uvádí, že

Likewise, Argyris a Kesley, publikovali ve roce 1960 …

O muži jménem Likewise jsem měl od počátku velké pochybnosti, přesto jsem šel hledat poučení na internetu. Na adrese

http://books.google.cz/books?id=dQE-aq6JJlQC&pg=PA3&lpg=PA3&dq=Likewise,+Argyris,+Kesley&source=bl&ots=VwjG_lOPZy&sig=Baaka6PDn0hAKBcfcRtzWiJOtVY&hl=cs&ei=_7koS6nGHZPCmgOGteCwDQ&sa=X&oi=book_result&ct=result&resnum=7&ved=0CDYQ6AEwBg#v=onepage&q=&f=false

jsem našel publikaci

Intermediate Finite Element Method: Fluid Flow and Heat Transfer Applications

by Juan C. Heinrich and Darrell W. Proper

jak je dobr um t anglicky
Jak je dobré umět anglicky

kde v oddíle 1.2 na straně 3 a 4 se uvádí

Likewise, Argyris and Kesley published a text in 1960 …

Jenomže, příslovce „likewise“ – tedy „podobně“ – se v angličtině, na rozdíl od češtiny, odděluje čárkou. Je zde však na začátku věty a je tedy s velkým „L“, což by autora – při pečlivém přebírání informací z cizích pramenů – nemělo přimět k víře, že existuje muž jménem Likewise.

Na přiložené reprodukci je inkriminovaná věta podtržena červeně. Je zřejmé, že dává smysl jen s větou předchozí.

Výše zmíněná publikace autorů Juan C. Heinrich a Darrell W. Proper, z níž autor disertační práce doslova převzal a špatně přeložil citovaný text, není uvedena v seznamu literatury.

v tomto kursu mimo jin nab dneme
V tomto kursu,mimo jiné, nabídneme

Porozumění metodě konečných prvků a numerické matematice,

a to prostřednictvím tvorby jednoduchých prográmků na koleně

Například nalezení sil v prutech, reakcí, vlastních frekvencí, vlastních tvarů kmitu, apod

v tomto kursu t nab dneme
V tomto kursu též nabídneme
  • Informaci o práce s daty, jejich zobrazení a seznámíme se se statistickými nástroji pro analýzou výsledků

bio2.ppt

slide13

Omezíme se na Newtonovskou fyziku – odhlédneme jak od kvantové fyziky, tak od teorie relativity

asto p id v me omezuj c podm nky
Často přidáváme omezující podmínky

Homogenní kontinuum – identické vlastnosti všech materiálových částic

Isotropní kontinuum – některé vlastnosti jsou nezávislé na směru

slide20

Moudrost vs. Trivialita

Otázka pravdy ve fyzice

what are we contributing to
What are we contributing to?
  • No fundamental laws and principles since Newton’s time
  • Newtonian physics – low velocities
  • Continuum mechanics – no quantum microcosms
  • Rather more sophisticated models, that either work or do not …
  • The question of truth is irrelevant … only the model proved by a proper experiment is acceptable
continuum mechanics 1
Continuum mechanics_1
  • The notion of continuum is one of possible models of matter.
  • The continuity of a structure as we observe is an illusion.
  • In liquids the molecules are loosely bound together by weak electrical forces. The molecules posses a considerable mobility.
  • In gases the intermolecular forces are even weaker
  • In metals there are relatively strong interatomic forces.
continuum mechanics 2
Continuum mechanics_2
  • Continuum mechanics ignores all the five details mentioned above and assumes that the discontinuous structure of real material is considered continuous. So the physical properties of material contained within an infinitesimal element are assumed to be the same as those determined experimentally on samples of finite dimensions.
  • Of course in view of the molecular and atomic structure of the matter the last assumption is false.
  • So the continuum is a model. It could give you the correct results if it is used within the limits of its applicability.
  • Latter on, we will show that FEM is just another model, with its own limits of validity.
  • The question is under what circumstances the continuum model provides a valid description of the flow and deformation of real material.
continuum mechanics 3
Continuum mechanics_3

It is not possible to give a satisfactory mathematical discussion of the validity of the continuum theory.

The ultimate justification of the model is empirical.

So in solid continuum mechanics for metals it is assumed that

if the linear dimension of volume element is greater than

10 000 times the interatomic distance, i.e. 1/1 000 000 [m],

then the continuum theory could still be safety used.

HUNTER, S.C.: Mechanics of Continuous Media,

Ellis Hornwood Ltd, UK, 1983

governing equations of solid continuum mechanics
Governing equations of solid continuum mechanics
  • Cauchy equations of motion
  • Kinematic relations
  • Constitutive relations

3 equations

6 equations

6 equations

pro pro rovnov hu materi lov ho elementu uva ujeme jen 3 rovnice
Proč pro rovnováhu materiálového elementu uvažujeme jen 3 rovnice?

When the stress components – associated with individual cube faces – are being defined and evaluated the material element is considered as a 3D body – the cube.

When equilibrium conditions are considered the material element is considered to be a point.

In developing the partial differential equations of motion, only the equilibrium of the forces was considered. The assumption that the resultant of moments of all forces about the origin must be equal zero can only be used to prove the symmetry of the stress tensor.

and here comes the idea of cosserat brothers
And here comes the idea of Cosserat brothers

Besides the force-stress tensors [Pa]

there are also couple-stress tensors [Pa m]

taken into account when the equilibrium conditions are considered

The Cosserat continuum is usually only effective when there exists a physical motivation for adding couple stresses and microcurvatures as is the case in granular materials.

For a numerical implementation of Cosserat continuum see Sluys and de Borst. See Stein, E., de Borst, R., Hughes, T.J.R.: Encyclopedia of Computational mechanics, Vol. 2, p. 355.

Eugène-Maurice-Pierre Cosserat (4 Mar 1866 – 31 May 1931) was a Frenchmathematician and astronomer. Born in Amiens, he studied at the École Normale Supérieure from 1883 to 1888. He was on Science faculty of Toulouse University from 1889 and director of its observatory from 1908, a position he held for the rest of his life. He was elected to the Académie des Sciences in 1919. His studies included the rings and satellites of Saturn, comets and double stars, but is best remembered for work with his engineer brother François on surface mechanics, particularly problems of elasticity.

co otcov zakladatel
Co otcové zakladatelé?

Teoretické základy mechaniky kontinua jsou známy po více než sto padesát let – Cauchy, Euler, St. Venant, ...

cauchy equations
Cauchy equations

Oeuvres complètes d'Augustin Cauchy. Série 2, tome 8 / publiées...Cauchy, Augustin-Louis (1789-1857), 1882-1974

http://gallica.bnf.fr/ark:/12148/bpt6k90200c.image.f4.langEN

slide31

Strain tensor components

rather body forces

evolution of stress notations
Evolution of stress notations

Todhunter, I. and Pearson, K.:

A History of the Theory of Elasticity, Dover Publications, New York, 1960.

ekvivalence vs rovnov ha
Ekvivalence

Jedna soustava sil = Druhá soustava sil

Rovnováha

Součet sil = 0

Ekvivalence vs. rovnováha

Formálně ano, ale …

Použili jsme d’Alambertova principu,

Z hlediska inerciálního systému je ta síla je fiktivní …

Rovnice platí jen právě teď

Takže bych tomu neříkal podmínky rovnováhy, ale pohybové rovnice

back to fem t oday approximate methods of solution s prevail
Back to FEMToday, approximate methods of solutions prevail

They are based on discretization in spaceand time and have numerous variants

  • Finite difference method
  • Transfer matrix method
  • Matrix methods
  • Finite element method
    • Displacement formulation
    • Force formulation
    • Hybrid formulation
  • Boundary element method
  • Meshless element method
  • From CAD and FEA to Isogeometric Analysis. By Cottrell, J.A., Hughes, TJ.R., Bazilevs, Z.: Isogeometric Analysis. Towards Integration of CAD and FEA, Wiley, Chichester, 2009.
numerical methods in finite element analysis
Numerical methods in Finite Element Analysis

Equilibrium problems

Space discretization only – Solution of systems of algebraic equations

Steady state vibration problems

Generalized eigenvalue problem

Propagation problems

Space and time discretization – step by step integration in time, for example CD (central differences) or NM (Newmark).

Timestep of integration corresponds to sampling increment in experiment.

Nyquist frequency plays the same role both in experiment and in computation.

Robust procedures and their efficient implementation are crucial for solving ‘large sized’ tasks typical for transient FE analysis

computability limits of are based on
Computability limits of are based on
  • Limits of physics
    • Limits of technology
    • Instrumental limits of
      • computers
      • experiments
  • But first of all on
    • Validity limits of employed models
      • Continuum mechanics is a model
      • Computational treatment is another model
      • Experiment is a tool for observing the nature

– but not the nature itself

experiments and axioms
Experiments and axioms

The body of theory furnishes the concept and formulae by means of which the experiment can be conceived and interpreted.

From experiment we may find agreement which develops confidence in the theory – but establish a theory by experiment we never can.

Experiment is a necessary adjunct to a physical theory – but it is an adjunct, not the master.

No experiment can be interpreted without recourse to ideas that are a part of the theory under examination

Adjunct – asistent, výpomocná síla

přídavek, doplněk

Fluegge, S.: (Editor) Encyclopedia of Physics, Vol. III, Principles of Classical Mechanics and Field Theory, Springer, Berlin, 1960 – Truesdell, C. and Toupin, R.: The Classical Field Theories, p. 228

slide41

Roger Bacon

On Experimental Science (1268)

Experimental science does not receive truth from superior science. She is the mistress and the other sciences are her servants.

… experimental science is a study entirely unknown by the common people …

…no science can be known without mathematics …

Roger Bacon, (c. 1214–1294), also known as Doctor Mirabilis (wonderful teacher), was an Englishphilosopher and Franciscan friar who placed considerable emphasis on empirical methods. He is sometimes credited as one of the earliest European advocates of the modern scientific method inspired by the works of Plato and Aristotle via early Islamic scientists such as Avicenna and Averroes.

From Wikipedia

slide43

Not a vector in a physical or mathematical sense, ie.

the quantity defined by the magnitude and the direction

details about tensor and matrix notation rules and terminology
Details about tensor and matrix notation, rules and terminology
  • Tensor rank
  • Kronecker delta
  • Summation rule
  • Orthogonal transformation
  • Addition, subtraction
  • Contraction
  • Outer and inner products
  • Scalar and dyadic products

See also

cm_part_1.ppt

Okrouhlík, M., Pták, S.: Počítačová mechanika kontinua I,

Základy nelineární mechaniky kontinua, Česká technika, Nakladatelství ČVUT, 2006

basic principles of solid continuum mechanics
Basic principles of solid continuum mechanics
  • Preliminaries
  • Gradient
  • Gauss divergence theorem
  • The generalization of ‘per partes’ integration (integration by parts)
  • Kinetic and strain energies
  • Material derivative
  • Conservation laws
  • Equilibrium
  • Cauchy equation of motion

See also

all_together_01_05_c4.doc

Okrouhlík, M., Pták, S.: Počítačová mechanika kontinua I,

Základy nelineární mechaniky kontinua, Česká technika, Nakladatelství ČVUT, 2006

doporu en literatura
Doporučená literatura
  • Bathe, K.-J.: Finite Element Procedures,

Prentice-Hall, Inc., Englewood Cliffs, 1996

  • Belytschko T., Liu, W.K., Moran, K.:

Nonlinear Finite Elements for Continua and Structures,

John Wiley, Chichester, 2000

  • Fung, Z.C., Tong, P.:

Classical and computational solid mechanics,

World Scientific, Singapore, 2001

  • Okrouhlík, M., Pták, S.:

Počítačová mechanika I, Základy nelineární mechaniky kontinua,

Nakladatelství ČVUT, 2006

  • Stejskal, V., Okrouhlík, M.:

Kmitání s Matlabem, (Vibration with Matlab),

Vydavatelství ČVUT, Praha 2002, ISBN 80-01-02435-0