Isogeometric Analysis Application to Car Crash Simulation
|
|
- Maurice Carson
- 5 years ago
- Views:
Transcription
1 Isogeometric Analysis Application to Car Crash Simulation S. Bouabdallah 2, C. Adam 2,3, M. Zarroug 2, H. Maitournam 3 De Vinci Engineering Lab, École Supérieure d Ingénieurs Léonard de Vinci 2 Direction Scientifique et Technologies Futures, PSA Peugeot Citroën 3 Laboratoire de Mécanique des Solides, École Polytechnique Forum Teratec 24 - July 2, 24
2 Outline Introduction Industrial context Numerical locking 2 Timoshenko beams Pure bending problem Membrane and transverse shear locking D B-Spline reduced quadrature rules 3 Reissner-Mindlin plates Pure bending problem 2D B-Spline reduced quadrature rules 4 Reissner-Mindlin shells Shell model Shell obstacle course 5 Non-linear solid shells Benchmark problems 6 Conclusions S. Bouabdallah, C. Adam, M. Zarroug, H. Maitournam Teratec 24 - IGA Application to Car Crash Simulation / 26
3 Industrial context Good predictive capability expected. Complete mesh size: 3 million elements. = Explicit simulations performed on the whole model. Stable time step: µs. hours to simulate 3ms on 96 CPUs. = Excessive computational cost. 9% of the geometry meshed with linear Reissner-Mindlin shell elements. Figure: Model size evolution S. Bouabdallah, C. Adam, M. Zarroug, H. Maitournam Teratec 24 - IGA Application to Car Crash Simulation 2 / 26
4 Industrial context Geometry has proved to be very significant. Need to ine some components. = Stable time step greatly penalized. Accurate and smooth geometry perred for contact problems. Large time used for preparing and meshing parts Figure: Estimation of relative time costs for model generation and analysis process S. Bouabdallah, C. Adam, M. Zarroug, H. Maitournam Teratec 24 - IGA Application to Car Crash Simulation 3 / 26
5 Industrial context = IsoGeometric Analysis for automobile crashworthiness. Parametrized CAD geometry for optimisation process. Analysis and Design share the same model. Higher regularity between elements within a patch. Figure: Numerical simulation process in isogeometric analysis. S. Bouabdallah, C. Adam, M. Zarroug, H. Maitournam Teratec 24 - IGA Application to Car Crash Simulation 4 / 26
6 Numerical locking Fully integrated low order elements, based on Timoshenko/Mindlin hypothesis, present poor performances with thin beams, plates and shells. Transverse shear and membrane locking appear. Classical Lagrange-based reduced and selective integration rules do not remove these pathologies. = Need to define B-Spline/NURBS-based integration schemes. Framework: Multi-patch problems with uniform regularity within each patch. Quadratic and cubic polynomials, regularity r from to p. C r C S. Bouabdallah, C. Adam, M. Zarroug, H. Maitournam Teratec 24 - IGA Application to Car Crash Simulation 5 / 26
7 Introduction Industrial context Numerical locking 2 Timoshenko beams Pure bending problem Membrane and transverse shear locking D B-Spline reduced quadrature rules 3 Reissner-Mindlin plates Pure bending problem 2D B-Spline reduced quadrature rules 4 Reissner-Mindlin shells Shell model Shell obstacle course 5 Non-linear solid shells Benchmark problems 6 Conclusions S. Bouabdallah, C. Adam, M. Zarroug, H. Maitournam Teratec 24 - IGA Application to Car Crash Simulation 6 / 26
8 Pure bending problem Thick curved beam of length L and radius R clamped at one end. Moment of inertia for a rectangular section: I = bh3 2. Distributed loading moment proportional to the bending stiffness: m(s) = EI ( π 2L) 2 sin ( π 2L s ). = Exact solution known, depends only of L and R: u n (s) = 2L ( ( π s ) ( π )) ( 2L cos cos πr )2 R 2L s. Exact membrane and shear strains are zero. Severe test for membrane and shear locking. S. Bouabdallah, C. Adam, M. Zarroug, H. Maitournam Teratec 24 - IGA Application to Car Crash Simulation 7 / 26
9 Membrane and transverse shear locking One element with B-Splines of degree p from to 5. Shear locking for slenderness ratio L/h and membrane locking for L/h R/L. Both shear and membrane locking appear for all degrees p. Rotation relative error with R/L = 3 Rotation relative error with R/L = 3 2 θ b ex θb h L 2 / θ b ex L SHEAR LOCKING MEMBRANE LOCKING Linear Quad. Cubic Quartic Quintic 2 R/L L/h θ b ex θb h L 2 / θ b ex L 2 SHEAR LOCKING MEMBRANE LOCKING Linear 2 Quad. Cubic 3 Quartic 4 Quintic 2 R/L L/h (a) Linear scale (b) Logarithmic scale Figure: One element curved thick beam: rotation relative error in L 2 norm. S. Bouabdallah, C. Adam, M. Zarroug, H. Maitournam Teratec 24 - IGA Application to Car Crash Simulation 8 / 26
10 S. Bouabdallah, C. Adam, M. Zarroug, H. Maitournam 2 Teratec 24 - IGA Application to Car Crash Simulation 9 / 26 Transverse shear locking Different effects on convergence according to the couple (p,r). Degree : no pre-asymptotic convergence. Degree 2: slow pre-asymptotic convergence when C and no pre-asymptotic convergence with few elements when C. Spurious constraints strengthened by the high regularity. = Excessive bending stiffeness of the structure. θ b ex θb h L 2 / θ b ex L Rotation relative error with L/h = 3 Linear Quad. C Quad C Cubic C Cubic C Number of elements / l (a) Pre-asymptotic convergence with L/h = θ b ex θb h L 2 / θ b ex L Rotation relative error with L/h = Linear Quad. C Quad C Cubic C Cubic C Number of elements / l (b) Pre-asymptotic and asymptotic convergence with L/h = 2 4 3
11 Mathematical framework Based on the work in FEA of [Prathap,93]. Field-consistency paradigm in thick structural elements. Transverse shear strain in straight beam: γ s (s) = w,s (s) θ(s). In the limit of an infinitely thin beam, the shear strain energy must vanish. = Different physical and spurious constraints are obtained according to the couple (p,r). Mathematical induction is performed adding elements one-by-one. New unknowns and conditions are numbered. An over-constrained linear system leads to numerical locking whereas an under-constrained system produces zero energy modes. S. Bouabdallah, C. Adam, M. Zarroug, H. Maitournam Teratec 24 - IGA Application to Car Crash Simulation / 26
12 D B-Spline reduced quadrature rules Optimal selective reduced integration schemes adapted to high regularity basis functions. The higher the regularity is, the fewer Gauss points are needed. Shear and membrane locking are completely removed. Improved accuracy of curved beam B-Spline elements. Order C Order 2 C Order 2 C Order 3 C Order 3 C Order 3 C2 Reduced integration for transverse shear: Gauss points ξ = s / L θ ex θ h L 2 / θ ex L C full C selec. C full C selec. Quadratic B spline 3 2 Number of elements 2 4 θ ex θ h L 2 / θ ex L Cubic B spline C full C selec. C full C selec. C 2 full C 2 selec. 2 Number of elements Figure: Rotation relative error in L 2 norm with full and selective reduced integrations S. Bouabdallah, C. Adam, M. Zarroug, H. Maitournam Teratec 24 - IGA Application to Car Crash Simulation / 26
13 Introduction Industrial context Numerical locking 2 Timoshenko beams Pure bending problem Membrane and transverse shear locking D B-Spline reduced quadrature rules 3 Reissner-Mindlin plates Pure bending problem 2D B-Spline reduced quadrature rules 4 Reissner-Mindlin shells Shell model Shell obstacle course 5 Non-linear solid shells Benchmark problems 6 Conclusions S. Bouabdallah, C. Adam, M. Zarroug, H. Maitournam Teratec 24 - IGA Application to Car Crash Simulation 2 / 26
14 Pure bending problem Same behaviours as with beams. High regularity wanted for more accurate results with fewer control points. Need to unlock B-Splines especially when C p regularity. Possibility to extend integration rules from beams to plates and shells. Vertical displacement of point C: concentrated load 8 w c / w c Linear C L/h = 3 Linear C L/h = 5 Quad. C L/h = 3 Quad. C L/h = 5 Quad. C L/h = 3 Quad. C L/h = Number of elements in each direction Figure: Central deflection of a simply supported plate (L/h = 3, 5 ). S. Bouabdallah, C. Adam, M. Zarroug, H. Maitournam Teratec 24 - IGA Application to Car Crash Simulation 3 / 26
15 Pure bending problem Same behaviours as with beams. High regularity wanted for more accurate results with fewer control points. Need to unlock B-Splines especially when C p regularity. Possibility to extend integration rules from beams to plates and shells. Vertical displacement of point C: uniform pressure w c / w c Linear C L/h = 3 Linear C L/h = 5 Quad. C L/h = 3 Quad. C L/h = 5 Quad. C L/h = 3 Quad. C L/h = Number of elements in each direction Figure: Central deflection of a simply supported plate (L/h = 3, 5 ). S. Bouabdallah, C. Adam, M. Zarroug, H. Maitournam Teratec 24 - IGA Application to Car Crash Simulation 3 / 26
16 2D B-Spline reduced quadrature rules Uni-dimensional B-Spline-based reduced quadrature rules extended to bi-dimensional rules by tensor product. The high regularity lowers the required number of Gauss points. Hourglass modes are suppressed using additional quadrature points in boundary elements, where the control points accumulate. B spline selective integration: Gauss points B spline selective integration: Gauss points B spline selective integration: Gauss points η.5 η.5 η ξ (a) Quadratic C ξ (b) Cubic C ξ (c) Figure: B-Spline selective integration quadrature rules and central deflection. S. Bouabdallah, C. Adam, M. Zarroug, H. Maitournam Teratec 24 - IGA Application to Car Crash Simulation 4 / 26
17 2D B-Spline reduced quadrature rules Resulting under-integrated elements are free from transverse shear locking. Numerical solutions independent of the thickness. Improved accuracy of plate isogeometric elements. No spurious zero energy modes. Vertical displacement of point C: concentrated load Vertical displacement of point C: concentrated load 8 8 w c / w c Linear C L/h = 3 Linear C L/h = 5 Quad. C L/h = 3 Quad. C L/h = 5 Quad. C L/h = 3 Quad. C L/h = Number of elements in each direction w c / w c Linear C L/h = 3 Linear C L/h = 5 Quad. C L/h = 3 Quad. C L/h = 5 Quad. C L/h = 3 Quad. C L/h = Number of elements in each direction (a) Full integration (b) Reduced integration Figure: Central deflection of a simply supported plate (L/h = 3, 5 ). S. Bouabdallah, C. Adam, M. Zarroug, H. Maitournam Teratec 24 - IGA Application to Car Crash Simulation 5 / 26
18 2D B-Spline reduced quadrature rules Resulting under-integrated elements are free from transverse shear locking. Numerical solutions independent of the thickness. Improved accuracy of plate isogeometric elements. No spurious zero energy modes. Vertical displacement of point C: uniform pressure Vertical displacement of point C: uniform pressure w c / w c Linear C L/h = 3 Linear C L/h = 5 Quad. C L/h = 3 Quad. C L/h = 5 Quad. C L/h = 3 Quad. C L/h = Number of elements in each direction w c / w c Linear C L/h = 3 Linear C L/h = 5 Quad. C L/h = 3 Quad. C L/h = 5 Quad. C L/h = 3 Quad. C L/h = Number of elements in each direction (c) Full integration (d) Reduced integration Figure: Central deflection of a simply supported plate (L/h = 3, 5 ). S. Bouabdallah, C. Adam, M. Zarroug, H. Maitournam Teratec 24 - IGA Application to Car Crash Simulation 5 / 26
19 Introduction Industrial context Numerical locking 2 Timoshenko beams Pure bending problem Membrane and transverse shear locking D B-Spline reduced quadrature rules 3 Reissner-Mindlin plates Pure bending problem 2D B-Spline reduced quadrature rules 4 Reissner-Mindlin shells Shell model Shell obstacle course 5 Non-linear solid shells Benchmark problems 6 Conclusions S. Bouabdallah, C. Adam, M. Zarroug, H. Maitournam Teratec 24 - IGA Application to Car Crash Simulation 6 / 26
20 Shell model Reissner-Mindlin shell model obtained from a degenerated three-dimensional model. First-order kinematic description through the thickness with transverse shear. Exact geometry of the shell defined by x(ξ) = nm A= R A (ξ, η)x A + h ζn(ξ, η). 2 Interpolated displacements described by u(ξ) = nm A= R A (U A + h ) 2 ζθ A n. S. Bouabdallah, C. Adam, M. Zarroug, H. Maitournam Teratec 24 - IGA Application to Car Crash Simulation 7 / 26
21 Characterization of the normal Several possibilities: Exact normal: n = x,ξ x,η x,ξ x,η 2. 2 Collocated normal: n = n A A [, nm ]. a. n A defined by the orthogonal projection of the associated control point X A onto the mid-surface of the shell, b. n A defined by a uniform distribution of the shell normals in the parametric space, c. n A defined at the Greville abscissae in the parametric space. Exact normal formulation is incompatible with reduced integration. Collocation at Greville abscissae results in good compromise between accuracy and computation time. S. Bouabdallah, C. Adam, M. Zarroug, H. Maitournam Teratec 24 - IGA Application to Car Crash Simulation 8 / 26
22 Shell obstacle course problems Evaluate time efficiency and accuracy of the proposed under-integrated elements. Linear elasticity benchmark problems. Bending dominated problems: transverse shear and membrane locking. (a) Scordelis-Lo roof (b) Pinched cylinder (c) Pinched hemisphere Figure: The shell obstacle course: problem descriptions and data. S. Bouabdallah, C. Adam, M. Zarroug, H. Maitournam Teratec 24 - IGA Application to Car Crash Simulation 9 / 26
23 Accuracy Vertical displacement of point B : full integration Vertical displacement of point C: full integration Horizontal displacement of point A: full integration w B / w B Linear C Quad. C Quad. C Cubic C Cubic C w C / w C Linear C Quad. C Quad. C Cubic C Cubic C u A / u A Linear C Quad. C Quad. C Cubic C Cubic C Number of elements per side Number of elements per side Number of elements per side (a) Roof: full (b) Cylinder: full (c) Hemisphere: full Vertical displacement of point B : reduced integration Vertical displacement of point C: reduced integration Horizontal displacement of point A: reduced integration w B / w B Linear C Quad. C Quad. C Cubic C Cubic C w C / w C Linear C Quad. C Quad. C Cubic C Cubic C u A / u A Linear C Quad. C Quad. C Cubic C Cubic C Number of elements per side Number of elements per side Number of elements per side (d) Roof: reduced (e) Cylinder: reduced (f) Hemisphere: reduced Figure: Shell obstacle course: displacement convergence S. Bouabdallah, C. Adam, M. Zarroug, H. Maitournam Teratec 24 - IGA Application to Car Crash Simulation 2 / 26
24 Time efficiency Vertical displacement of point B : full integration Vertical displacement of point C: full integratione Horizontal displacement of point A: full integration w B / w B Linear C Quad. C Quad. C Cubic C Cubic C CubicC 2 w C / w C Linear C Quad. C Quad. C Cubic C Cubic C u A / u A Linear C Quad. C Quad. C Cubic C Cubic C 2 Computational time (s) 2 Computational time (s) 2 Computational time (s) (a) Roof: full (b) Cylinder: full (c) Hemisphere: full Vertical displacement of point B : reduced integration Vertical displacement of point C: reduced integration Horizontal displacement of point A: reduced integration w B / w B 4 2 Linear C Quad. C Quad. C Cubic C Cubic C w C / w C 4 2 Linear C Quad. C Quad. C Cubic C Cubic C u A / u A 4 2 Linear C Quad. C Quad. C Cubic C Cubic C 2 Computational time (s) 2 Computational time (s) 2 Computational time (s) (d) Roof: reduced (e) Cylinder: reduced (f) Hemisphere: reduced Figure: Shell obstacle course: displacement convergence S. Bouabdallah, C. Adam, M. Zarroug, H. Maitournam Teratec 24 - IGA Application to Car Crash Simulation 2 / 26
25 Total number of Gauss points Lagrange-based quadrature: n bending GP = ( e ( )) 2p+ 2 2 B-Spline-based quadrature: and n shear GP = ( e ( )) 2p n bending GP = ngp shear = ( e ( ) ) 2p+ 2 2 r + 2r. Computational effort depends on the difference p r when performing B-Spline-based reduced integration. Speed up factor greater than 5 (resp. 6) is reached with the under-integrated quadratic (resp. cubic) C p elements. Computation time: full integration Computation time: selective integration Computation time (s) Linear C Quad. C Quad. C Cubic C Cubic C Computation time (s) Linear C Quad. C Quad. C Cubic C Cubic C Number of elements per side (a) Full integration Number of elements per side (b) Selective integration Figure: Computation time for the pinched cylinder problem. S. Bouabdallah, C. Adam, M. Zarroug, H. Maitournam Teratec 24 - IGA Application to Car Crash Simulation 22 / 26
26 Introduction Industrial context Numerical locking 2 Timoshenko beams Pure bending problem Membrane and transverse shear locking D B-Spline reduced quadrature rules 3 Reissner-Mindlin plates Pure bending problem 2D B-Spline reduced quadrature rules 4 Reissner-Mindlin shells Shell model Shell obstacle course 5 Non-linear solid shells Benchmark problems 6 Conclusions S. Bouabdallah, C. Adam, M. Zarroug, H. Maitournam Teratec 24 - IGA Application to Car Crash Simulation 23 / 26
27 Benchmark problems Solid-shell elements with one quadratic element in the thickness. Geometric non-linear elasticity benchmark problems. Bi-variate B-Spline-based quadrature rules applied. (a) Pinched hemisphere (b) Stretched cylinder Figure: Geometric non-linear analysis of shells: problem descriptions and data. S. Bouabdallah, C. Adam, M. Zarroug, H. Maitournam Teratec 24 - IGA Application to Car Crash Simulation 24 / 26
28 Accuracy and time efficiency Bi-quadratic C shape functions in the in-plane directions. Improved accuracy and time efficiency. Speed up factor of approximately 3. Load P in % of P max (N) Load displacement curve 2 2x2 full 2x2 reduced Shell erence S4R Displacement (m) (a) Pinched hemisphere Load P in % of P max (N) Load displacement curve 2x2 full 2x2 reduced Shell erence S4R Displacement (m) (b) Stretched cylinder Figure: Displacement convergence with Total Lagrangian Formulation S. Bouabdallah, C. Adam, M. Zarroug, H. Maitournam Teratec 24 - IGA Application to Car Crash Simulation 25 / 26
29 Conclusions Severe numerical locking appear in quadratic and cubic C p elements. = Reduced integration to improve the performances of IG elements. Classical reduced/selective quadrature rules remove numerical locking in C elements only. Mathematical framework to define B-Spline-based reduced D quadrature rules. Extension to 2D quadrature rules by tensor product. = Efficient thick shell elements obtained as in FEA. Improved accuracy and time efficiency. Computational cost depends on the difference p r. S. Bouabdallah, C. Adam, M. Zarroug, H. Maitournam Teratec 24 - IGA Application to Car Crash Simulation 26 / 26
Isogeometric analysis of thin Reissner-Mindlin plates and shells: Locking phenomena and generalized local B method
Isogeometric analysis of thin Reissner-Mindlin plates and shells: Locking phenomena and generalized local B method Qingyuan Hu a,b, Yang Xia c, Sundararajan Natarajan d, Andreas Zilian b, Ping Hu c, Stéphane
More informationCurrent Status of Isogeometric Analysis in LS-DYNA
Current Status of Isogeometric Analysis in LS-DYNA Stefan Hartmann Developer Forum, September 24 th, 2013, Filderstadt, Germany Development in cooperation with: D.J. Benson: Professor of Applied Mechanics,
More informationAN ISOGEOMETRIC REISSNER-MINDLIN SHELL WITH LAGRANGE BASIS
11th World Congress on Computational Mechanics (WCCM XI) 5th European Conference on Computational Mechanics (ECCM V) 6th European Conference on Computational Fluid Dynamics (ECFD VI) E. Oñate, J. Oliver
More informationANSYS Element. elearning. Peter Barrett October CAE Associates Inc. and ANSYS Inc. All rights reserved.
ANSYS Element Selection elearning Peter Barrett October 2012 2012 CAE Associates Inc. and ANSYS Inc. All rights reserved. ANSYS Element Selection What is the best element type(s) for my analysis? Best
More informationLOCAL STRESS ANALYSIS OF STIFFENED SHELLS USING MSC/NASTRAN S SHELL AND BEAM p-elements
LOCAL STRESS ANALYSIS OF STIFFENED SHELLS USING MSC/NASTRAN S SHELL AND BEAM p-elements Sanjay Patel, Claus Hoff, Mark Gwillim The MacNeal-Schwendler Corporation Abstract In large finite element models
More informationSolid and shell elements
Solid and shell elements Theodore Sussman, Ph.D. ADINA R&D, Inc, 2016 1 Overview 2D and 3D solid elements Types of elements Effects of element distortions Incompatible modes elements u/p elements for incompressible
More informationTesting of Shell Elements using Challenging Benchmark Problems
Testing of Shell Elements using Challenging Benchmark Problems F.T. Wong Assistant Professor, Department of Civil Engineering, Petra Christian University, Surabaya ABSTRACT This paper presents the author
More informationElement Order: Element order refers to the interpolation of an element s nodal results to the interior of the element. This determines how results can
TIPS www.ansys.belcan.com 鲁班人 (http://www.lubanren.com/weblog/) Picking an Element Type For Structural Analysis: by Paul Dufour Picking an element type from the large library of elements in ANSYS can be
More informationRecent Advances on Higher Order 27-node Hexahedral Element in LS-DYNA
14 th International LS-DYNA Users Conference Session: Simulation Recent Advances on Higher Order 27-node Hexahedral Element in LS-DYNA Hailong Teng Livermore Software Technology Corp. Abstract This paper
More informationPATCH TEST OF HEXAHEDRAL ELEMENT
Annual Report of ADVENTURE Project ADV-99- (999) PATCH TEST OF HEXAHEDRAL ELEMENT Yoshikazu ISHIHARA * and Hirohisa NOGUCHI * * Mitsubishi Research Institute, Inc. e-mail: y-ishi@mri.co.jp * Department
More informationModeling of Shells with Three-dimensional Finite Elements
Universität Stuttgart Fakultät für Bau und- Umweltingenieurwissenschaften Baustatik und Baudynamik Modeling of Shells with Three-dimensional Finite Elements Manfred Bischoff Institute of Structural Mechanics
More informationFigure 30. Degrees of freedom of flat shell elements
Shell finite elements There are three types of shell finite element; 1) flat elements, 2) elements based on the Sanders-Koiter equations and 3) elements based on reduction of a solid element. Flat elements
More informationGuidelines for proper use of Plate elements
Guidelines for proper use of Plate elements In structural analysis using finite element method, the analysis model is created by dividing the entire structure into finite elements. This procedure is known
More informationIntroduction to 2 nd -order Lagrangian Element in LS-DYNA
Introduction to 2 nd -order Lagrangian Element in LS-DYNA Hailong Teng Livermore Software Technology Corporation Nov, 2017 Motivation Users are requesting higher order elements for implicit. Replace shells.
More informationFinite Element Method. Chapter 7. Practical considerations in FEM modeling
Finite Element Method Chapter 7 Practical considerations in FEM modeling Finite Element Modeling General Consideration The following are some of the difficult tasks (or decisions) that face the engineer
More informationSet No. 1 IV B.Tech. I Semester Regular Examinations, November 2010 FINITE ELEMENT METHODS (Mechanical Engineering) Time: 3 Hours Max Marks: 80 Answer any FIVE Questions All Questions carry equal marks
More informationRevised Sheet Metal Simulation, J.E. Akin, Rice University
Revised Sheet Metal Simulation, J.E. Akin, Rice University A SolidWorks simulation tutorial is just intended to illustrate where to find various icons that you would need in a real engineering analysis.
More informationRevision of the SolidWorks Variable Pressure Simulation Tutorial J.E. Akin, Rice University, Mechanical Engineering. Introduction
Revision of the SolidWorks Variable Pressure Simulation Tutorial J.E. Akin, Rice University, Mechanical Engineering Introduction A SolidWorks simulation tutorial is just intended to illustrate where to
More informationMODELLING STIFFENED LIGHTWEIGHT STRUCTURES WITH ISOGEOMETRIC ANALYSIS VIA MORTAR METHODS
ECCOMAS Congress 2016 VII European Congress on Computational Methods in Applied Sciences and Engineering M. Papadrakakis, V. Papadopoulos, G. Stefanou, V. Plevris (eds.) Crete Island, Greece, 5 10 June
More informationIntroduction to Finite Element Analysis using ANSYS
Introduction to Finite Element Analysis using ANSYS Sasi Kumar Tippabhotla PhD Candidate Xtreme Photovoltaics (XPV) Lab EPD, SUTD Disclaimer: The material and simulations (using Ansys student version)
More informationCOMPUTER AIDED ENGINEERING DESIGN (BFF2612)
COMPUTER AIDED ENGINEERING DESIGN (BFF2612) BASIC MATHEMATICAL CONCEPTS IN CAED by Dr. Mohd Nizar Mhd Razali Faculty of Manufacturing Engineering mnizar@ump.edu.my COORDINATE SYSTEM y+ y+ z+ z+ x+ RIGHT
More informationChapter 3 Analysis of Original Steel Post
Chapter 3. Analysis of original steel post 35 Chapter 3 Analysis of Original Steel Post This type of post is a real functioning structure. It is in service throughout the rail network of Spain as part
More informationCOMPUTER AIDED ENGINEERING. Part-1
COMPUTER AIDED ENGINEERING Course no. 7962 Finite Element Modelling and Simulation Finite Element Modelling and Simulation Part-1 Modeling & Simulation System A system exists and operates in time and space.
More informationFlexible multibody systems - Relative coordinates approach
Computer-aided analysis of multibody dynamics (part 2) Flexible multibody systems - Relative coordinates approach Paul Fisette (paul.fisette@uclouvain.be) Introduction In terms of modeling, multibody scientists
More informationNon-Linear Finite Element Methods in Solid Mechanics Attilio Frangi, Politecnico di Milano, February 3, 2017, Lesson 1
Non-Linear Finite Element Methods in Solid Mechanics Attilio Frangi, attilio.frangi@polimi.it Politecnico di Milano, February 3, 2017, Lesson 1 1 Politecnico di Milano, February 3, 2017, Lesson 1 2 Outline
More informationStatic analysis of an isotropic rectangular plate using finite element analysis (FEA)
Journal of Mechanical Engineering Research Vol. 4(4), pp. 148-162, April 2012 Available online at http://www.academicjournals.org/jmer DOI: 10.5897/JMER11.088 ISSN 2141-2383 2012 Academic Journals Full
More informationMeshless Modeling, Animating, and Simulating Point-Based Geometry
Meshless Modeling, Animating, and Simulating Point-Based Geometry Xiaohu Guo SUNY @ Stony Brook Email: xguo@cs.sunysb.edu http://www.cs.sunysb.edu/~xguo Graphics Primitives - Points The emergence of points
More informationMotivation Patch tests Numerical examples Conclusions
Weakening the tight coupling between geometry and simulation in isogeometric analysis: from sub- and super- geometric analysis to Geometry Independent Field approximation (GIFT) Elena Atroshchenko, Gang
More informationSmooth finite elements
Smooth finite elements seamless handling of incompressibility, distorted and polygonal meshes; links with equilibrium methods Stéphane Bordas * Nguyen-Xuan Hung ** Nguyen-Dang Hung *** * University of
More informationCITY AND GUILDS 9210 UNIT 135 MECHANICS OF SOLIDS Level 6 TUTORIAL 15 - FINITE ELEMENT ANALYSIS - PART 1
Outcome 1 The learner can: CITY AND GUILDS 9210 UNIT 135 MECHANICS OF SOLIDS Level 6 TUTORIAL 15 - FINITE ELEMENT ANALYSIS - PART 1 Calculate stresses, strain and deflections in a range of components under
More informationGenerative Part Structural Analysis Fundamentals
CATIA V5 Training Foils Generative Part Structural Analysis Fundamentals Version 5 Release 19 September 2008 EDU_CAT_EN_GPF_FI_V5R19 About this course Objectives of the course Upon completion of this course
More informationModelling Flat Spring Performance Using FEA
Modelling Flat Spring Performance Using FEA Blessing O Fatola, Patrick Keogh and Ben Hicks Department of Mechanical Engineering, University of Corresponding author bf223@bath.ac.uk Abstract. This paper
More informationThe Application of. Interface Elements to Dissimilar Meshes in. Global/Local Analysis
The Application of Interface Elements to Dissimilar Meshes in Global/Local Analysis John E Schiermeier Senior Development Engineer The MacNeal Schwendler Corporation Los Angeles, California Jerrold M Housner
More informationElfini Solver Verification
Page 1 Elfini Solver Verification Preface Using this Guide Where to Find More Information Conventions What's new User Tasks Static Analysis Cylindrical Roof Under its Own Weight Morley's Problem Twisted
More informationFINITE ELEMENT ANALYSIS OF A COMPOSITE CATAMARAN
NAFEMS WORLD CONGRESS 2013, SALZBURG, AUSTRIA FINITE ELEMENT ANALYSIS OF A COMPOSITE CATAMARAN Dr. C. Lequesne, Dr. M. Bruyneel (LMS Samtech, Belgium); Ir. R. Van Vlodorp (Aerofleet, Belgium). Dr. C. Lequesne,
More informationApplication of Shell elements to buckling-analyses of thin-walled composite laminates
Application of Shell elements to buckling-analyses of thin-walled composite laminates B.A. Gӧttgens MT 12.02 Internship report Coach: Dr. R. E. Erkmen University of Technology Sydney Department of Civil
More informationA Multiple Constraint Approach for Finite Element Analysis of Moment Frames with Radius-cut RBS Connections
A Multiple Constraint Approach for Finite Element Analysis of Moment Frames with Radius-cut RBS Connections Dawit Hailu +, Adil Zekaria ++, Samuel Kinde +++ ABSTRACT After the 1994 Northridge earthquake
More informationIsogeometric Analysis of Solids & Structures Small and Finite Deformation Applications
Trends & Challenges in Computational Mechanics Padua, Italy - September 12-14, 2011 Isogeometric Analysis of Solids & Structures Small and Finite Deformation Applications Robert L. Taylor Department of
More informationReduction of Finite Element Models for Explicit Car Crash Simulations
Reduction of Finite Element Models for Explicit Car Crash Simulations K. Flídrová a,b), D. Lenoir a), N. Vasseur b), L. Jézéquel a) a) Laboratory of Tribology and System Dynamics UMR-CNRS 5513, Centrale
More informationElastic Analysis of a Bending Plate
analys: linear static. constr: suppor. elemen: plate q12pl. load: elemen face force. materi: elasti isotro. option: direct units. post: binary ndiana. pre: dianai. result: cauchy displa extern force green
More informationModule 1: Introduction to Finite Element Analysis. Lecture 4: Steps in Finite Element Analysis
25 Module 1: Introduction to Finite Element Analysis Lecture 4: Steps in Finite Element Analysis 1.4.1 Loading Conditions There are multiple loading conditions which may be applied to a system. The load
More informationAnalysis and Design of Cantilever Springs
Analysis and Design of Cantilever Springs Hemendra Singh Shekhawat, Hong Zhou Department of Mechanical Engineering Texas A&M University-Kingsville Kingsville, Texas, USA Abstract Cantilever springs are
More informationEXACT BUCKLING SOLUTION OF COMPOSITE WEB/FLANGE ASSEMBLY
EXACT BUCKLING SOLUTION OF COMPOSITE WEB/FLANGE ASSEMBLY J. Sauvé 1*, M. Dubé 1, F. Dervault 2, G. Corriveau 2 1 Ecole de technologie superieure, Montreal, Canada 2 Airframe stress, Advanced Structures,
More informationCHAPTER 2 LITERATURE REVIEW. The finite element method today, has become a powerful tool in engineering analysis
CHAPTER 2 LITERATURE REVIEW 2.1 General The finite element method today, has become a powerful tool in engineering analysis and design. The ease of application, reliability of solutions and ability to
More informationINSTRUCTIONAL PLAN L( 3 ) T ( ) P ( ) Instruction Plan Details: DELHI COLLEGE OF TECHNOLOGY & MANAGEMENT(DCTM), PALWAL
DELHI COLLEGE OF TECHNOLOGY & MANAGEMENT(DCTM), PALWAL INSTRUCTIONAL PLAN RECORD NO.: QF/ACD/009 Revision No.: 00 Name of Faculty: Course Title: Theory of elasticity L( 3 ) T ( ) P ( ) Department: Mechanical
More informationFB-MULTIPIER vs ADINA VALIDATION MODELING
FB-MULTIPIER vs ADINA VALIDATION MODELING 1. INTRODUCTION 1.1 Purpose of FB-MultiPier Validation testing Performing validation of structural analysis software delineates the capabilities and limitations
More informationMAE 323: Lab 7. Instructions. Pressure Vessel Alex Grishin MAE 323 Lab Instructions 1
Instructions MAE 323 Lab Instructions 1 Problem Definition Determine how different element types perform for modeling a cylindrical pressure vessel over a wide range of r/t ratios, and how the hoop stress
More informationModule: 2 Finite Element Formulation Techniques Lecture 3: Finite Element Method: Displacement Approach
11 Module: 2 Finite Element Formulation Techniques Lecture 3: Finite Element Method: Displacement Approach 2.3.1 Choice of Displacement Function Displacement function is the beginning point for the structural
More informationModeling lattice structured materials with micropolar elasticity. Accuracy of the micropolar model. Marcus Yoder. CEDAR presentation Spring 2017
Modeling lattice structured materials with micropolar elasticity Accuracy of the micropolar model CEDAR presentation Spring 2017 Advisors: Lonny Thompson and Joshua D. Summers Outline 2/25 1. Background
More informationModeling Skills Stress Analysis J.E. Akin, Rice University, Mech 417
Introduction Modeling Skills Stress Analysis J.E. Akin, Rice University, Mech 417 Most finite element analysis tasks involve utilizing commercial software, for which you do not have the source code. Thus,
More informationThe effect of parameterization on isogeometric analysis of free-form curved beams
The effect of parameterization on isogeometric analysis of free-form curved beams Hosseini, SF, Moetakef-Imani, B, Hadidi-Moud, S & Hassani, B Author post-print (accepted) deposited by Coventry University
More informationA Numerical Form Finding Method for Minimal Surface of Membrane Structure
10 th World Congress on Structural and Multidisciplinary Optimization May 19-24, 2013, Orlando, Florida, US Numerical Form Finding Method for Minimal Surface of Membrane Structure Koichi Yamane 1, Masatoshi
More informationTABLE OF CONTENTS SECTION 2 BACKGROUND AND LITERATURE REVIEW... 3 SECTION 3 WAVE REFLECTION AND TRANSMISSION IN RODS Introduction...
TABLE OF CONTENTS SECTION 1 INTRODUCTION... 1 1.1 Introduction... 1 1.2 Objectives... 1 1.3 Report organization... 2 SECTION 2 BACKGROUND AND LITERATURE REVIEW... 3 2.1 Introduction... 3 2.2 Wave propagation
More informationA Comprehensive Introduction to SolidWorks 2011
A Comprehensive Introduction to SolidWorks 2011 Godfrey Onwubolu, Ph.D. SDC PUBLICATIONS www.sdcpublications.com Schroff Development Corporation Chapter 2 Geometric Construction Tools Objectives: When
More informationScientific Manual FEM-Design 17.0
Scientific Manual FEM-Design 17. 1.4.6 Calculations considering diaphragms All of the available calculation in FEM-Design can be performed with diaphragms or without diaphragms if the diaphragms were defined
More informationIntroduction. Section 3: Structural Analysis Concepts - Review
Introduction In this class we will focus on the structural analysis of framed structures. Framed structures consist of components with lengths that are significantly larger than crosssectional areas. We
More informationWP1 NUMERICAL BENCHMARK INVESTIGATION
WP1 NUMERICAL BENCHMARK INVESTIGATION 1 Table of contents 1 Introduction... 3 2 1 st example: beam under pure bending... 3 2.1 Definition of load application and boundary conditions... 4 2.2 Definition
More informationBuckling Analysis of a Thin Plate
Buckling Analysis of a Thin Plate Outline 1 Description 2 Modeling approach 3 Finite Element Model 3.1 Units 3.2 Geometry definition 3.3 Properties 3.4 Boundary conditions 3.5 Loads 3.6 Meshing 4 Structural
More informationLS-DYNA s Linear Solver Development Phase 1: Element Validation
LS-DYNA s Linear Solver Development Phase 1: Element Validation Allen T. Li 1, Zhe Cui 2, Yun Huang 2 1 Ford Motor Company 2 Livermore Software Technology Corporation Abstract LS-DYNA is a well-known multi-purpose
More informationNumerical simulation of CAD thin structures using the extended Finite Element Method and Level Sets
Numerical simulation of CAD thin structures using the extended Finite Element Method and Level Sets G. Legrain a,, C. Geuzaine b, J.F. Remacle c, N. Moës a, P. Cresta d, J. Gaudin e a LUNAM Université,
More informationME 475 FEA of a Composite Panel
ME 475 FEA of a Composite Panel Objectives: To determine the deflection and stress state of a composite panel subjected to asymmetric loading. Introduction: Composite laminates are composed of thin layers
More informationMANUAL COMPUTERS & STRUCTURES INC.
SAP2000 Integrated Finite Element Analysis and Design of Structures VERIFICATION MANUAL COMPUTERS & STRUCTURES INC. Computers and Structures, Inc. Berkeley, California, USA Version 6.1 Revised July 1997
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Lecture - 36
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Lecture - 36 In last class, we have derived element equations for two d elasticity problems
More informationIsogeometric fluid structure interaction analysis with emphasis on non-matching discretizations, and with application to wind turbines
Isogeometric fluid structure interaction analysis with emphasis on non-matching discretizations, and with application to wind turbines Yuri Bazilevs a,, Ming-Chen Hsu a, and Michael A. Scott b a Department
More informationExercise 1. 3-Point Bending Using the GUI and the Bottom-up-Method
Exercise 1 3-Point Bending Using the GUI and the Bottom-up-Method Contents Learn how to... 1 Given... 2 Questions... 2 Taking advantage of symmetries... 2 A. Preprocessor (Setting up the Model)... 3 A.1
More informationMetafor FE Software. 2. Operator split. 4. Rezoning methods 5. Contact with friction
ALE simulations ua sus using Metafor eao 1. Introduction 2. Operator split 3. Convection schemes 4. Rezoning methods 5. Contact with friction 1 Introduction EULERIAN FORMALISM Undistorted mesh Ideal for
More informationCHAPTER 4. Numerical Models. descriptions of the boundary conditions, element types, validation, and the force
CHAPTER 4 Numerical Models This chapter presents the development of numerical models for sandwich beams/plates subjected to four-point bending and the hydromat test system. Detailed descriptions of the
More informationPro MECHANICA STRUCTURE WILDFIRE 4. ELEMENTS AND APPLICATIONS Part I. Yves Gagnon, M.A.Sc. Finite Element Analyst & Structural Consultant SDC
Pro MECHANICA STRUCTURE WILDFIRE 4 ELEMENTS AND APPLICATIONS Part I Yves Gagnon, M.A.Sc. Finite Element Analyst & Structural Consultant SDC PUBLICATIONS Schroff Development Corporation www.schroff.com
More informationLarge Rotation Analysis of thin-walled shells
Large Rotation Analysis of thin-walled shells Y. Başar, O. Kintzel Ruhr-University Bochum Bochum, Germany Summary The objective of this contribution is to develop an isoparametric -node shell element by
More informationMAE Advanced Computer Aided Design. 01. Introduction Doc 02. Introduction to the FINITE ELEMENT METHOD
MAE 656 - Advanced Computer Aided Design 01. Introduction Doc 02 Introduction to the FINITE ELEMENT METHOD The FEM is A TOOL A simulation tool The FEM is A TOOL NOT ONLY STRUCTURAL! Narrowing the problem
More informationEngineering Effects of Boundary Conditions (Fixtures and Temperatures) J.E. Akin, Rice University, Mechanical Engineering
Engineering Effects of Boundary Conditions (Fixtures and Temperatures) J.E. Akin, Rice University, Mechanical Engineering Here SolidWorks stress simulation tutorials will be re-visited to show how they
More informationUtilizing ABAQUS 10-Node Modified Tet for Analyzing Impact Problems Involving Thin-Walled Structures.
Utilizing ABAQUS 10-Node Modified Tet for Analyzing Impact Problems Involving Thin-Walled Structures. Abstract Ted Diehl* and Doug Carroll** *Mechanical Technology Center, PCS **Advanced Technology, Smart
More informationMultilevel optimization of. of Composite panels under complex load and boundary conditions.
Loughborough University Institutional Repository Multilevel optimization of composite panels under complex load and boundary conditions This item was submitted to Loughborough University's Institutional
More informationTHREE DIMENSIONAL ACES MODELS FOR BRIDGES
THREE DIMENSIONAL ACES MODELS FOR BRIDGES Noel Wenham, Design Engineer, Wyche Consulting Joe Wyche, Director, Wyche Consulting SYNOPSIS Plane grillage models are widely used for the design of bridges,
More informationME Optimization of a Frame
ME 475 - Optimization of a Frame Analysis Problem Statement: The following problem will be analyzed using Abaqus. 4 7 7 5,000 N 5,000 N 0,000 N 6 6 4 3 5 5 4 4 3 3 Figure. Full frame geometry and loading
More informationNEW WAVE OF CAD SYSTEMS AND ITS APPLICATION IN DESIGN
Vol 4 No 3 NEW WAVE OF CAD SYSTEMS AND ITS APPLICATION IN DESIGN Ass Lecturer Mahmoud A Hassan Al-Qadisiyah University College of Engineering hasaaneng@yahoocom ABSTRACT This paper provides some lighting
More informationIn this course we will need a set of techniques to represent curves and surfaces in 2-d and 3-d. Some reasons for this include
Parametric Curves and Surfaces In this course we will need a set of techniques to represent curves and surfaces in 2-d and 3-d. Some reasons for this include Describing curves in space that objects move
More informationShell-to-Solid Element Connector(RSSCON)
WORKSHOP 11 Shell-to-Solid Element Connector(RSSCON) Solid Shell MSC.Nastran 105 Exercise Workbook 11-1 11-2 MSC.Nastran 105 Exercise Workbook WORKSHOP 11 Shell-to-Solid Element Connector The introduction
More informationMultiframe May 2010 Release Note
Multiframe 12.02 18 May 2010 Release Note This release note describes the version 12.02 release of Multiframe, Steel Designer and Section Maker. This release will run on Windows XP/2003/Vista/7. Contents
More informationixcube 4-10 Brief introduction for membrane and cable systems.
ixcube 4-10 Brief introduction for membrane and cable systems. ixcube is the evolution of 20 years of R&D in the field of membrane structures so it takes a while to understand the basic features. You must
More informationChapter 7 Practical Considerations in Modeling. Chapter 7 Practical Considerations in Modeling
CIVL 7/8117 1/43 Chapter 7 Learning Objectives To present concepts that should be considered when modeling for a situation by the finite element method, such as aspect ratio, symmetry, natural subdivisions,
More informationCOSMOSWorks Verification
2006 COSMOSWorks Verification COSMOSWorks Verification Problems This document contains verification problems to demonstrate the accuracy of the COMSOSWorks software in comparison to analytical results.
More informationA review of the state-of-the-art in vehicle modeling for crashworthiness analysis using LSDYNA
A review of the state-of-the-art in vehicle modeling for crashworthiness analysis using LSDYNA We have about 15 years of experience and every crashworthiness simulation is a compromise between quality
More informationIntroduction to FEM Modeling
Total Analysis Solution for Multi-disciplinary Optimum Design Apoorv Sharma midas NFX CAE Consultant 1 1. Introduction 2. Element Types 3. Sample Exercise: 1D Modeling 4. Meshing Tools 5. Loads and Boundary
More informationVIRTUAL TESTING OF AIRCRAFT FUSELAGE STIFFENED PANELS
4 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES VIRTUAL TESTING OF AIRCRAFT FUSELAGE STIFFENED PANELS Peter Linde*, Jürgen Pleitner*, Wilhelm Rust** *Airbus, Hamburg, Germany, **CAD-FE GmbH, Burgdorf
More informationIndustrial Applications of Computational Mechanics Plates and Shells Mesh generation static SSI. Prof. Dr.-Ing. Casimir Katz SOFiSTiK AG
Industrial Applications of Plates and Shells Mesh generation static SSI Prof. Dr.-Ing. Casimir Katz SOFiSTiK AG FEM - Reminder A mathematical method The real (continuous) world is mapped on to a discrete
More information1. INTRODUCTION. A shell may be defined as the solid material enclosed. between two closely spaced curved surfaces (1), the distance
1 1. INTRODUCTION 1.1. GENERAL A shell may be defined as the solid material enclosed between two closely spaced curved surfaces (1), the distance between these two surfaces being the thickness of the shell.
More informationKU Leuven vibro-acoustics activities in an Industry 4.0 context
KU Leuven vibro-acoustics activities in an Industry 4.0 context Wim Desmet KU Leuven Department of Mechanical Engineering Flanders Make - Virtual Department Mechatronics & Design overview KU Leuven team
More informationSETTLEMENT OF A CIRCULAR FOOTING ON SAND
1 SETTLEMENT OF A CIRCULAR FOOTING ON SAND In this chapter a first application is considered, namely the settlement of a circular foundation footing on sand. This is the first step in becoming familiar
More informationDevelopment of Eight Node Isoparametric Quadrilateral (QUAD8) Elements in Java
709 Development of Eight Node Isoparametric Quadrilateral (QUAD8) Elements in Java I.E Umeonyiagu 1, N.P Ogbonna 2 1 Department of Civil Engineering, Chukwuemeka Odumegwu Ojukwu University, Uli. Nigeria
More informationExercise 2: Mesh Resolution, Element Shapes, Basis Functions & Convergence Analyses
Exercise 2: Mesh Resolution, Element Shapes, Basis Functions & Convergence Analyses Goals In this exercise, we will explore the strengths and weaknesses of different element types (tetrahedrons vs. hexahedrons,
More informationCourse in. FEM ANSYS Classic
Course in Geometric modeling Modeling Programme for Lesson: Modeling considerations Element Type Real Constants Material Properties Sections Geometry/Modeling WorkPlane & Coordinate systems Keypoints Lines
More informationENGINEERING TRIPOS PART IIA FINITE ELEMENT METHOD
ENGINEERING TRIPOS PART IIA LOCATION: DPO EXPERIMENT 3D7 FINITE ELEMENT METHOD Those who have performed the 3C7 experiment should bring the write-up along to this laboratory Objectives Show that the accuracy
More informationFinite Element Course ANSYS Mechanical Tutorial Tutorial 3 Cantilever Beam
Problem Specification Finite Element Course ANSYS Mechanical Tutorial Tutorial 3 Cantilever Beam Consider the beam in the figure below. It is clamped on the left side and has a point force of 8kN acting
More informationOptimization in the Abaqus Environment Using TOSCA
Optimization in the Abaqus Environment Using TOSCA Luca Furbatto, Giovanni Di Lorenzo and Claus B.W. Pedersen Luca Furbatto and Giovanni Di Lorenzo (McLaren Racing), Claus B.W. Pedersen (FE-Design GmbH)
More informationRefinable C 1 spline elements for irregular quad layout
Refinable C 1 spline elements for irregular quad layout Thien Nguyen Jörg Peters University of Florida NSF CCF-0728797, NIH R01-LM011300 T. Nguyen, J. Peters (UF) GMP 2016 1 / 20 Outline 1 Refinable, smooth,
More informationFinite Element Modeling Techniques (2) دانشگاه صنعتي اصفهان- دانشكده مكانيك
Finite Element Modeling Techniques (2) 1 Where Finer Meshes Should be Used GEOMETRY MODELLING 2 GEOMETRY MODELLING Reduction of a complex geometry to a manageable one. 3D? 2D? 1D? Combination? Bulky solids
More informationHourglass (HG) Modes. Hourglass (HG) Modes
Hourglass (HG) Modes Hourglass modes are nonphysical modes of deformation that occur in underintegrated elements and produce no stress. By underintegrated, we mean Solid elements with a single integration
More informationGeneral modeling guidelines
General modeling guidelines Some quotes from industry FEA experts: Finite element analysis is a very powerful tool with which to design products of superior quality. Like all tools, it can be used properly,
More informationComputations of stresses with volume-elements in rectangular and HE sections
CT3000: Bachelor Thesis Report, Izik Shalom (4048180) Computations of stresses with volume-elements in rectangular and HE sections Supervisors: dr. ir. P.C.J. Hoogenboom en Ir. R. Abspoel June 2013 Preface
More information