Purpose of the Course
Lecturer: Lars Andersen, MSc, PhD, Associate Professor Department of Civil Engineering, Division of Structural Mechanics A lb Aalborg U University, i it S Sohngaardsholmsvej h d h l j 57 57, DK DK-9000 9000 A Aalborg lb Phone: 9940 8455 | E-mail:
[email protected] Homepage of the course: www.wind.civil.aau.dk → Teaching Activities → Finite Element Design After the course, the student must be able to: Understand the basic concepts in Finite Element Analysis (FEA) Use a FEA program (STAAD.Pro 2007) Know K the th basic b i tterms and d algorithms l ith b behind hi d th the analysis l i Be able to analyse large complex structures Come up with realistic dimensions of structural elements.
1
Contents of the lecture
Introduction to the Finite-Element Method (FEM) ﻩWhat is Finite-Element Analysis (FEA)? ﻩHistoric overview ﻩWhy use FEA? p and output p from an FEA ﻩInput ﻩHow does FEM work? ﻩExample in STAAD.Pro ﻩExercise: Learn to use STAAD.Pro (plane-frame problem)
2
Short Historic Overview of FEM
FEM is tied with the development of computer technology Approximately 40 years old NASA developed NASTRAN in the 1960s First College Course in FEM was offered in 1970 In the 1970s, FEM was limited to large corporations with expensive mainframe computers In the 1980s, “powerful” desktop computers made FEM an indispensable engineering tool In the 1990s, more complex elements are introduced, optimization capabilities are integrated integrated, and CAD programs are used for modelling complex structures g y the method was developed p for the analysis y of Originally, stresses in structures – but today FEM is used to analyse heat transfer, fluid flow, electric and magnetic fields etc.
3
Why use FEA? Simple Statically Determinate System
Freeway-crossing north of Aalborg ...
... approximated i t db by statically t ti ll d determinate t i t system t
5 reaction forces determined by 5 equilibrium equations
4
Why use FEA? Simple Statically Indeterminate System
Freeway-crossing north of Aalborg ...
... approximated i t db by statically t ti ll iindeterminate d t i t system t
6 reaction forces but only 3 equilibrium equations
5
Why use FEA? Complex Structure
6
David Fay Custom Chair
Golden Gate Bridge, San Francisco
Antwerp Railway Station, Belgium
Kandahar Airport, Afghanistan
Why use FEA? Load combinations
7
Discrete Versus Continuous System
8
Prototype = reality
Continuous system
Discrete system = FE model
Input and Output from an FEA
Input: ﻩDimensions Di i off the h structure ﻩCross-section types (circular, rectangular, I-profile, ...) ﻩMaterial properties (wood, steel, concrete, glass, ...) ﻩSupports (fixed, free, moving, ...) ﻩLoads (concentrated, line, surface, combinations) Output: ﻩDeformation components (translation, rotation) ﻩSection force curves, reactions (shear, normal, moment) ﻩStrains and stresses (shear, normal) ﻩFulfilment of design criteria (Eurocode, ...) ﻩEigenmodes (dynamic resonance risk, ...)
9
Eigenmodes/Eigenfrequencies
Tacoma Narrows Bridge, 1940
10
Eigenmodes/Eigenfrequencies
Millenium Bridge London, 2000/2002
11
Eigenmodes/Eigenfrequencies
12
FEA Programs
Commercial Finite-Element Programs: ﻩABAQUS (www.simulia.com) ﻩCOSMOSWorks (SolidWorks) (www.cosmosm.com) ﻩFEMLAB (www.comsol.dk) ﻩSTAAD.Pro (www.bentley.com) ﻩANSYS Structural (www.ansys.com) ﻩetc. etc
13
How Does FEA work? Dividing the structure into elements
14
ﻩFrom the user input, a given structure is divided into small elements (finite elements) (done partly by user and partly by program) ﻩEach element is assigned material properties (done by user) ﻩEach element’s element s mechanical behaviour is defined by a set of differential equations from the choice of element type and material properties (done by program)
How Does FEA work? Matrix Equations for the Elements are Found
ﻩThe differential equations for each element are solved and arranged d iinto t a matrix t i fformulation l ti suitable it bl ffor computert aided solutions (done by program)
15
How Does FEA work? Matrix Equation for the Global System is Assembled
ﻩThe element matrices are combined into a global system of equations ti d defined fi d ffrom th the placement l t off th the respective ti elements (done by program) ﻩFrom this the global structural equation is obtained (done by program)
16
How Does FEA work? Load and Boundary Conditions are Applied
ﻩThe boundary conditions (loads and supports) are specified (d (done b by user)) ﻩThe boundary conditions are incorporated into the system of differential equations (done by program)
17
How Does FEA work? The Structural Matrix Equation is Solved
ﻩThe displacement (and rotations) of all nodes are found from solving th system the t off equations ti (done (d b by program)) ﻩDisplacements at intermediate points are found from interpolation of nodal values (done by program)
18
How Does FEA work? Stresses and Strains are found
ﻩThe strains are found from the displacements (done by program) ﻩStresses are found from a constitutive relation (done by program)
19
How Does FEA work? The design criteria are checked
ﻩThe design criteria are checked (done by user/program) ﻩThe structure is modified to fulfil criteria and a new analysis is made (done by user)
20
Example – Analytical Solution
Structural system:
Analytical solution:
21
STAAD.Pro – Overview
22
Editor
Add beam Mark beam
Hold down ctrl to move the starting point of a beam
Menu
S Snap node d
STAAD.Pro – Overview
General, General property, support, l d load
23
D t b Database
STAAD.Pro – Overview
Analyse/print All
24
STAAD.Pro – Script file
25
Node coordinates Member definition M t i ld Material definition fi iti Section assignment Material assignment Support assignment Load assignment Analysis definition
STAAD.Pro – Analysis
Analyze
Mode
26
STAAD.Pro – Results
Double click gives Double-click section displacement
Node displacement Reactions Section forces
27
Example – Numerical Result from STAAD.Pro
Analytical solution:
28
Today's Problem
29
Plane bridge
Determine ete e tthe ep profile o e types from o tthe e de deformation o at o ccriteria te a ((1/200 / 00 o of spa span). ) Change the supports and determine the profile types in the same manner. Get familiar with the p program. g ((Use the menu Geometry/Split y p Beam to divide the vertical beam into 3 for easy applying the load)
Today's Problem
30
Support types
Profile types
Maximum deformation
(F)ixed, (P)ined, (F)ixed (B)ut (direction)
(Global deformation), local (x,y)-coordinates (x,y) coordinates
+ FB (Fx,Mz)
P
F
FB (Fx,Mz)
P
F
F
F
F
P
P
P
50 mm
IPE160
HE200B
y-dir di 139.2 mm x=4.167
y-dir di 165.6 mm x=5.833
x-dir di 33.8 mm y=2.85+1.45
50 mm 50 mm
Maximum allowed deformation 1/200 of span = 50 mm