16.810 (16.682) Engineering Design and Rapid Prototyping Instructor(s) Finite Element Method January 12, 2004 Prof. Olivier de Weck Dr. Il Yong Kim 16.810 (16.682) 2 Plan for Today ? FEM Lecture (ca. 50 min) ? FEM fundamental concepts, analysis procedure ? Errors, Mistakes, and Accuracy ? Cosmos Introduction (ca. 30 min) ? Follow along step-by-step ? Conduct FEA of your part (ca. 90 min) ? Work in teams of two ? First conduct an analysis of your CAD design ? You are free to make modifications to your original model 16.810 (16.682) 3 Course Concept today 16.810 (16.682) 4 Course Flow Diagram CAD/CAM/CAE Intro FEM/Solid Mechanics Overview Manufacturing Training Structural Test “Training” Design Optimization Hand sketching CAD design FEM analysis Produce Part 1 Test Produce Part 2 Optimization Problem statement Final Review Test Learning/Review Deliverables Design Sketch v1 Analysis output v1 Part v1 Experiment data v1 Design/Analysis output v2 Part v2 Experiment data v2 Drawing v1 Design Intro due today today Wednesday 16.810 (16.682) 5 Numerical Method Finite Element Method Boundary Element Method Finite Difference Method Finite Volume Method Meshless Method 16.810 (16.682) 6 What is the FEM? Description - FEM cuts a structure into several elements (pieces of the structure). - Then reconnects elements at “nodes” as if nodes were pins or drops of glue that hold elements together. - This process results in a set of simultaneous algebraic equations. FEM: Method for numerical solution of field problems. Number of degrees-of-freedom (DOF) Continuum: Infinite FEM: Finite (This is the origin of the name, Finite Element Method) 16.810 (16.682) 7 Fundamental Concepts (1) Elastic problems Thermal problems Fluid flow Electrostatics etc. Many engineering phenomena can be expressed by “governing equations” and “boundary conditions” Governing Equation (Differential equation) () 0LfI  Boundary Conditions () 0BgI  16.810 (16.682) 8 Elastic deformation Thermal behavior etc. Governing Equation: () 0LfI  Boundary Conditions: () 0BgI  []{}{} Ku F A set of simultaneous algebraic equations FEM Approximate! Fundamental Concepts (2) Example: Vertical machining center Geometry is very complex! You know all the equations, but you cannot solve it by hand 16.810 (16.682) 9 []{}{} Ku F 1 {} [ ]{}  uKF Property Behavior Action Elastic Thermal Fluid Electrostatic Behavior {}uProperty []K Action {}F stiffness displacement force conductivity temperature heat source viscosity velocity body force dialectri permittivity electric potential charge Unknown Fundamental Concepts (3) 16.810 (16.682) 10 It is very difficult to make the algebraic equations for the entire domain Divide the domain into a number of small, simple elements Adjacent elements share the DOF at connecting nodes Fundamental Concepts (4) Finite element: Small piece of structure A field quantity is interpolated by a polynomial over an element 16.810 (16.682) 11 Obtain the algebraic equations for each element (this is easy!) Put all the element equations together []{}{} Ku F []{}{} E EE Ku F []{}{} E EE Ku F []{}{} E EE Ku F []{}{} E EE Ku F []{}{} E EE Ku F []{}{} E EE Ku F []{}{} E EE Ku F []{}{} E EE Ku F []{}{} E EE Ku F Fundamental Concepts (5) 16.810 (16.682) 12 []{}{} Ku F 1 {} [ ]{}  uKF Solve the equations, obtaining unknown variables at nodes. Fundamental Concepts (6) 16.810 (16.682) 13 Concepts - Summary - FEM uses the concept of piecewise polynomial interpolation. - By connecting elements together, the field quantity becomes interpolated over the entire structure in piecewise fashion. - A set of simultaneous algebraic equations at nodes. []{}{} Ku F Property Behavior Action K: Stiffness matrix x: Displacement F: Load KxF= K F x 16.810 (16.682) 14 Brief History - The term finite element was first coined by clough in 1960. In the early 1960s, engineers used the method for approximate solutions of problems in stress analysis, fluid flow, heat transfer, and other areas. - The first book on the FEM by Zienkiewicz and Chung was published in 1967. - In the late 1960s and early 1970s, the FEM was applied to a wide variety of engineering problems. - Most commercial FEM software packages originated in the 1970s. (Abaqus, Adina, Ansys, etc.) - Klaus-Jurgen Bathe in ME at MIT Reference [2] 16.810 (16.682) 15 Can readily handle very complex geometry: - The heart and power of the FEM Can handle a wide variety of engineering problems - Solid mechanics - Dynamics - Heat problems - Fluids - Electrostatic problems Can handle complex restraints - Indeterminate structures can be solved. Can handle complex loading - Nodal load (point loads) - Element load (pressure, thermal, inertial forces) - Time or frequency dependent loading Advantages of the FEM 16.810 (16.682) 16 Disadvantages of the FEM A general closed-form solution, which would permit one to examine system response to changes in various parameters, is not produced. The FEM obtains only "approximate" solutions. The FEM has "inherent" errors. Mistakes by users can be fatal. 16.810 (16.682) 17 Typical FEA Procedure by Commercial Software Preprocess Process Postprocess User User Computer Build a FE model Conduct numerical analysis See results 16.810 (16.682) 18 [1] Select analysis type - Structural Static Analysis - Modal Analysis - Transient Dynamic Analysis -Buckling Analysis - Contact - Steady-state Thermal Analysis - Transient Thermal Analysis [2] Select element type 2-D 3-D Linear Quadratic Beam Truss Shell Solid Plate [3] Material properties ,,,,E QUD" Preprocess (1) 16.810 (16.682) 19 Preprocess (2) [4] Make nodes [5] Build elements by assigning connectivity [6] Apply boundary conditions and loads 16.810 (16.682) 20 Process and Postprocess - Solve the boundary value problem [7] Process - See the results [8] Postprocess Displacement Stress Strain Natural frequency Temperature Time history 16.810 (16.682) 21 Responsibility of the user Results obtained from ten reputable FEM codes and by users regarded as expert can differ greatly.* BC: Hinged supports Load: Pressure pulse * R. D. Cook, Finite Element Modeling for Stress Analysis, John Wiley & Sons, 1995 Fancy, colorful contours can be produced by any model, good or bad!! 1 ms pressure pulse 200 mm Unknown: Lateral mid point displacement in the time domain 16.810 (16.682) 22 Errors Inherent in FEM Formulation Quadratic element Cubic element - Field quantity is assumed to be a polynomial over an element. (which is not true) True deformation - Geometry is simplified. Domain Approximated domain FEM Linear element FEM - Use very simple integration techniques (Gauss Quadrature) x f(x) 1 -1 1 1 11 Area: ( ) 33 fxdx f f  §·§ · | ¨?¨ ? ?1? 1 3 16.810 (16.682) 23 - The computer carries only a finite number of digits. - Numerical Difficulties e.g.) Very large stiffness difference e.g.) 2 1.41421356, 3.14159265S 12 2 ,0kk k| 12 22 2 2 [( ) ] 0 PP kk ku P u k  ? | Errors Inherent in Computing u 2 u 1 k 1 k 2 P 12 16.810 (16.682) 24 Mistakes by Users - Elements are of the wrong type e.g) Shell elements are used where solid elements are needed - Distorted elements - Supports are insufficient to prevent all rigid-body motions - Inconsistent units (e.g. E=200 GPa, Force = 100 lbs) - Too large stiffness differences ? Numerical difficulties 16.810 (16.682) 25 Plan for Today ? FEM Lecture (ca. 50 min) ? FEM fundamental concepts, analysis procedure ? Errors, Mistakes, and Accuracy ? Cosmos Introduction (ca. 30 min) ? Follow along step-by-step ? Conduct FEA of your part (ca. 90 min) ? Work in teams of two ? First conduct an analysis of your CAD design ? You are free to make modifications to your original model 16.810 (16.682) 26 References Glaucio H. Paulino, Introduction to FEM (History, Advantages and Disadvantages), http://cee.ce.uiuc.edu/paulino Robert Cook et al., Concepts and Applications of Finite Element Analysis, John Wiley & Sons, 1989 Robert Cook, Finite Element Modeling For Stress Analysis, John Wiley & Sons, 1995 Introduction to Finite Element Method, http://210.17.155.47 (in Korean) J. Tinsley Oden et al., Finite Elements – An Introduction, Prentice Hall, 1981