Mandatory programming project

Steak, shake or break –
and other applications of
FEM
TMA4220
Outline





Present the governing equations (3 of them)
Building your FEM code
Mesh generation
Experimentation
Integration with other courses
Time dependent Poisson


“Steak”: Heat transfer
The way heat propagates through objects
u
 u 0
t
Linear elasticity


“Break”: Stress analysis
How things bend and deform
 ( u )  f
Vibration analysis


“Shake”: Eigenmode analysis
Natural frequency of which things oscillate
 u   (u ) 0
Building your FEM code
Series of problems 1a) – 3b)
 Gauss quadrature





FEM theory
Sample 2d problem
Sample 3d problem
Your choice and experimentation (Problem 4)
How to work

On problem 1-3: answer short and answer right

On problem 4: Elaborate as much as possible
Mesh generation


Lots of issues (probably the most understated
part of FEM training)
3 options
Pre-made by myself
 GMSH (free, easy-to-use mesh tool, readable output)


I’ll get back to #3 later
Experimentation

Find your own problem
Choose an equation
 Choose a geometry
 Choose an objective
 Analyze the results


Experiment with geometry, material properties
and/or boundary conditions
Example 1: beef

How to cook the optimal beef




Choose an equation
Choose a geometry
Choose an objective
 Optimality criterium?
Experimentation



Too short cooking – raw and unedible
Too long cooking – dry and tasteless
Too varying temperature – all of the above
uu()225
 u 0
t
Example 2: Structural analysis

Linear elasticity



Choose an equation
Choose a geometry
 http://www.youtube.com/watch?v=v5fS4RF3ht0
Choose
an objective


Max stress lower than material critical point
Will it withstand gravity?
Will it withstand external forces? (creepers?)
 ( u )  f

Example 3: free vibration

What are the natural frequencies of your geometry (and
material)?



Choose an equation
Choose an objective
 Find the natural frequencies
Choose
a geometry



Resonance with enviromental frequencies (wind, noise etc)
Sound/music (harmonic overtones)
Purely visual
 u   (u ) 0
Integration with other courses



Numerical linear algebra (TMA4205)
Mathematical Modeling (TMA4195)
Experts in Team (TMA4850)
Practicalities






2 persons to each group
Deadline: 18th Nov.
Written report and presentation
No min/max pages
Access to computer lab
I’m away next week, Anne away the next 2
weeks
Summary



Choose your own problem
Be creative
Play on your strengths
Geometry modeling
 Visualization
 Physical Experiments
 CPU speed
