DOCSLIB.ORG

Computational Mechanics Featuring Matlab DRAFT EDITION $Revision: 1.52

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Computational Mechanics Featuring Matlab DRAFT EDITION $Revision: 1.52 Computational Mechanics featuring Matlab DRAFT EDITION $Revision: 1.52 Richard Sonnenfeld November 4, 2013 c 2012 Richard Sonnenfeld Published by: New Mexico Tech Press All rights reserved Contents 0.1 Introduction–forthestudent. i 0.1.1 Howtoreadthisbook..................... i 0.1.2 Extra credit for corrections! . i 0.1.3 Why computer programming in physics . ii 0.1.4 Otherreferences ........................ ii 0.2 Introduction – for the instructor . iii 0.2.1 To the instructor who is new to Matlab or programming . iii 0.2.2 Intended audience . iii 0.2.3 Coursedesign.......................... iv 0.2.3.1 Lecture portion . iv 0.2.3.2 Lab portion . iv 0.2.3.3 Facilities required . v 0.3 Featuresofthisbook.......................... v 0.3.1 Generalfeatures ........................ v 0.3.2 Course timing . vi 0.3.3 Projects............................. vi 0.3.4 Additional topics . vii 0.3.5 Chapterbychaptercontent . vii 1 Hello World! 1 1.1 Stepsto“HelloWorld”......................... 1 1.2 “HelloWorld”.............................. 1 1.3 Gettinghelp............................... 2 1.4 BrowsingHelp ............................. 3 1.5 Matlab interactivemode ........................ 3 1.6 Operatorsandhowtocalculate . 5 1.6.0.1 Exercise 2A . 6 1.6.0.2 Exercise 2B . 6 1.6.0.3 Exercise 2C . 6 1.7 VariablesandMemory ......................... 6 1.7.1 “=”doesnotmean“Equal”. 7 1.7.2 VariableTypes ......................... 8 3 4 CONTENTS 1.7.3 Scientific Notation . 8 1.8 Usingformattedoutput ........................ 8 1.9 Combiningshortscripts . .. .. .. .. .. .. .. .. .. .. 11 1.9.1 script and m-file ........................ 12 1.10 Afirstinteractiveprogram. 12 1.11 ∗Variables,Memory,andbinarycodes . 12 1.12 Review of commands introduced in this chapter . 14 1.13 EndofChapterProblems . 14 2 IWillnotThrowMassiveProjectilesinClass 17 2.1 Getting homework done quickly with loops . 17 2.2 Usingloopstomakeatable . .. .. .. .. .. .. .. .. .. 18 2.3 Projectile motion for heavy projectiles . 19 2.3.1 Using loops to make a table of projectile motion . 20 2.4 Conditionalstatements. 20 2.4.1 Using conditional statements to stop a loop . 20 2.4.2 Using conditional statements to allow choices . 22 2.4.3 Using conditional statements in a while loop ........ 23 2.5 CreatingASCIIdatafiles .. .. .. .. .. .. .. .. .. .. 24 2.6 LogicalorRelationaloperators . 25 2.7 Morebuilt-infunctions. .. .. .. .. .. .. .. .. .. .. 26 2.7.1 Checking the size of an error – Absolute value . 26 2.7.2 Checking for divisibility – fix() and mod() ......... 27 2.7.3 Roundingfunctions.. .. .. .. .. .. .. .. .. .. 27 2.8 What’snext? .............................. 27 2.9 Reviewofcommandsintroducedinthischapter . 28 2.10 EndofChapterProblems . 29 2.11 Codeexamplesforthischapter . 32 3 Arrays, Matrices and Functions 35 3.1 ArraysandMatrices .......................... 35 3.1.1 Arrays.............................. 36 3.1.1.1 Defining and referencing arrays . 36 3.1.1.2 Using the diary() function to log your exploration 37 3.1.2 RowandColumnArrays(1-DMatrices) . 37 3.1.2.1 Managing matrices with the size() function . 38 3.1.3 UsingArraystorepresentvectors. 39 3.1.3.1 Finding the angle between vectors of arbitrary di- mension........................ 40 3.1.4 Using Arrays to represent sets of measurements . 40 3.1.5 Comparing vectorized calculations to element-by-element calculations........................... 42 3.1.5.1 Initialization . 42 3.1.5.2 Assignment...................... 42 3.1.6 More documentation about arrays and matrices . 43 3.2 Userdefinedfunctions . .. .. .. .. .. .. .. .. .. .. 44 3.2.1 HelponUserDefinedfunctions . 45 CONTENTS 5 3.2.2 Differencebetweenascriptandafunction . 45 3.3 Reviewofcommandsintroducedinthischapter . 46 3.4 EndofChapterProblems . .. .. .. .. .. .. .. .. .. 47 3.5 Codeexamplesforthischapter . 50 4 Working with Scientific Data 53 4.1 Workingwithrealdata ........................ 54 4.1.1 Reading and Writing Data files . 54 4.1.1.1 AboutMetadata . .. .. .. .. .. .. .. 56 4.1.2 Working with matrices . 56 4.2 Plotting data with Matlab ....................... 59 4.2.1 Basic plotting . 59 4.2.2 Improvingtheappearanceofbasicplots . 60 4.2.3 Savingfiguresasimages . 61 4.2.4 Advanced plotting: Objects and handle graphics . 62 4.3 NumericalDerivatives . .. .. .. .. .. .. .. .. .. .. 65 4.4 Reviewofcommandsintroducedinthischapter . 66 4.5 EndofChapterProblems . .. .. .. .. .. .. .. .. .. 67 4.6 Codeexamplesforthischapter . 74 5 Projectile Motion with Drag 81 5.1 Trajectoryandrangeofaheavyprojectile . 82 5.2 More Plotting data with Matlab .................... 83 5.2.1 Making publication quality figures . 84 5.2.2 Adding authorship information . 86 5.3 Document, Define, Derive, Display (D4)............... 87 5.4 Simple animation with Matlab ..................... 88 5.5 ViscosityandStoke’sLaw . 90 5.5.1 Defining Viscosity . 90 5.5.2 Examples of Viscosity . 91 5.5.3 Twotypesofviscosity . .. .. .. .. .. .. .. .. 92 5.5.4 WhatcausesViscosity? . 92 5.5.5 Stokes’Law........................... 93 5.6 Linear Drag in One-Dimension, Terminal Velocity . 93 5.7 Linear Drag in Two Dimensions . 96 5.8 Reviewofcommandsintroducedinthischapter . 96 5.9 EndofChapterProblems . .. .. .. .. .. .. .. .. .. 97 5.10 Codeexamplesforthischapter . .100 6 QuadraticDragandtheEulerMethod 103 6.1 InertialDrag ..............................104 6.1.1 Whatcausesinertialdrag? . .104 6.2 ReynoldsNumber............................106 6.2.1 Fluid mechanics and dimensionless parameters . 106 6.2.2 Reynoldsnumberanddragregimes. .107 6.3 One-dimensional (1-D) analytic solution for quadratic drag . 108 6.3.1 Hyperbolic functions . 109 6 CONTENTS 6.4 EulerMethod..............................109 6.5 ApplyingEulerMethodtoQuadraticdrag . .110 6.5.1 Onedimensionalcase . .110 6.5.2 Matlabcodefor1-DEulerMethod . .111 6.5.3 Using sprintf() toannotateplots . .112 6.5.4 Quadratic drag in 2-Dimensions . 112 6.5.5 Getting position from velocity . 113 6.6 EndofChapterProblems . .114 6.7 Codeexamplesforthischapter . .117 7 Newton’s Universal Law of Gravitation 119 7.1 Introduction...............................120 7.1.1 Linear and Angular motion . 120 7.2 Newton’sLawofUniversalGravitation. .120 7.2.1 VectorformofLawofGravitation . .122 7.3 Kepler’sLaws..............................123 7.3.1 Introduction ..........................123 7.3.2 Planets move in ellipses with the sun at one focus . 123 7.3.3 Kepler’s Period Law . 124 7.3.3.1 Fictitious forces – What holds the planets up? . 124 7.3.3.2 Inertialreferenceframes. .126 7.3.3.3 Uniform circular motion . 126 7.3.3.4 Orbits in inertial reference frames . 128 7.3.3.5 Kepler’s period law for circular orbits . 129 7.3.4 Planetssweepout equalareasin equaltimes . 129 7.3.4.1 Rotational form of Newton’s 2nd Law . 130 7.3.4.2 Conservation of Angular Momentum . 130 7.3.4.3 Equivalence of Kepler’s second law and angular momentumconservation. .131 7.4 GravitationalPotential. .132 7.4.1 EscapeVelocity.........................133 7.4.1.1 Example: Calculating escape velocity from Earth 133 7.5 GravitationofExtendedBodies. .134 7.6 CenterofMass .............................137 7.6.1 Example: Calculating center of mass of four point-masses . 137 7.6.2 Correcting Kepler’s period law . 138 7.6.3 Reduced Mass (µ) .......................139 7.6.4 TheBohrAtom.........................140 7.7 Review..................................142 7.8 EndofChapterProblems . .142 8 Runge-KuttaMethodandOrbitalSimulation 149 8.1 Advantagesofhigher-ordermethods . .149 8.2 Firstorder Runge-Kutta(RK1 orEuler method) . 150 8.3 SecondorderRunge-Kutta(RK2) . .151 8.3.1 RK2AlgebraicInterpretation . .152 8.3.2 RK2GraphicalInterpretation. .153 CONTENTS 7 8.4 FourthorderRunge-Kutta(RK4). .153 8.5 Applying Runge-Kutta methods to quadratic drag . 154 8.6 Applying Runge-Kutta methods to gravitation . 157 8.6.1 Sample listing for Runge-Kutta2 applied to gravitation . 157 8.7 EndofChapterProblems . .161 8.8 PossibleProjectsforPhysics241 . .168 8.8.1 Preproposal...........................168 8.8.2 Proposal.............................168 8.8.3 Difficulty and Grading . 168 8.8.4 TheProjects ..........................168 8.8.4.1 Modeling . 168 8.8.4.2 Data Analysis . 171 8 CONTENTS 0.1. INTRODUCTION – FOR THE STUDENT i 0.1 Introduction – for the student 0.1.1 How to read this book In this age of hyperlinks and non-linear thinking, I believe learning can still be best served by an engaging story that starts at the beginning and progresses logically to the end. Hyperlinks are good when you are already an expert in an area and you know just what you are looking for. The story is best when exploring a world that is new to you; you do not yet know what you are looking for. While some of the students taking this course may have computer programming experience, none is expected in advance. The best way to read a physics (or math, or engineering text) is with pencil and paper, so that you repeat (or fill in) all the details of derivations and examples that are given in the book. Psychologists have long known that “active learning” is the only way to keep your mind engaged in a new and difficult subject. Just scanning the words on the page is fine for a novel, but not for learning a technical subject. When the subject is explicitly computer programming, one should have, in lieu of pencil and paper, an open terminal, or in this case, and open Matlab command line. This book in no way tries to be exhaustive in its coverage of the Matlab language itself. In fact Matlab is huge, and it grows every year. Since it is a commercial program, it is in the interests
Load more

Recommended publications
  • Verification and Validation in Computational Fluid Dynamics1
    SAND2002 - 0529 Unlimited Release Printed March 2002 Verification and Validation in Computational Fluid Dynamics1 William L. Oberkampf Validation and Uncertainty Estimation Department Timothy G. Trucano Optimization and Uncertainty Estimation Department Sandia National Laboratories P. O. Box 5800 Albuquerque, New Mexico 87185 Abstract Verification and validation (V&V) are the primary means to assess accuracy and reliability in computational simulations. This paper presents an extensive review of the literature in V&V in computational fluid dynamics (CFD), discusses methods and procedures for assessing V&V, and develops a number of extensions to existing ideas. The review of the development of V&V terminology and methodology points out the contributions from members of the operations research, statistics, and CFD communities. Fundamental issues in V&V are addressed, such as code verification versus solution verification, model validation versus solution validation, the distinction between error and uncertainty, conceptual sources of error and uncertainty, and the relationship between validation and prediction. The fundamental strategy of verification is the identification and quantification of errors in the computational model and its solution. In verification activities, the accuracy of a computational solution is primarily measured relative to two types of highly accurate solutions: analytical solutions and highly accurate numerical solutions. Methods for determining the accuracy of numerical solutions are presented and the importance of software testing during verification activities is emphasized. The fundamental strategy of 1Accepted for publication in the review journal Progress in Aerospace Sciences. 3 validation is to assess how accurately the computational results compare with the experimental data, with quantified error and uncertainty estimates for both.
    [Show full text]
  • Computational Physics * Cruz a Editorial Office, Pure and Applied Physics, India
    Research & Reviews: Journal of Pure and Applied Physics e-ISSN:2320-2459 p-ISSN:2347-2316 Computational Physics * Cruz A Editorial Office, Pure and Applied Physics, India Editorial Note Received date: 21/02/2021 *For Correspondence Accepted date: 23/02/2021 Published date: 28/02/2021 Andrew Cruz, Computational Physics, Editorial Office, Pure and Applied Physics, India. E-mail: [email protected] EDITORIAL NOTE Computational physics is that the study and implementation of numerical analysis to unravel problems in physics that a quantitative theory already exists. Historically, computational physics was the primary application of recent computers in science, and is now a subset of computational science. It is sometimes considered a sub discipline (or offshoot) of theoretical physics, but others consider it an intermediate branch between theoretical and experimental physics-a neighbourhood of study which supplements both theory and experiment. Computational physics problems are generally very difficult to unravel exactly. This is often thanks to several (mathematical) reasons: lack of algebraic and/or analytic solvability, complexity, and chaos. For example, even apparently simple problems, like calculating the wave function of an electron orbiting an atom during a strong field (Stark effect), may require great effort to formulate a practical algorithm (if one are often found); other cruder or brute-force techniques, like graphical methods or root finding, could also be required. On the more advanced side, mathematical perturbation theory is additionally sometimes used (a working is shown for this particular example here). In addition, the computational cost and computational complexity for many-body problems (and their classical counterparts) tend to grow quickly.
    [Show full text]
  • Fast and Efficient MATLAB-Based MPM Solver
    Fast and efficient MATLAB-based MPM solver (fMPMM-solver v1.1) Emmanuel Wyser1, Yury Alkhimenkov1,2,3, Michel Jaboyedoff1,2, and Yury Y. Podladchikov1,2,3 1Institute of Earth Sciences, University of Lausanne, 1015 Lausanne, Switzerland 2Swiss Geocomputing Centre, University of Lausanne, 1015 Lausanne, Switzerland 3Department of Mechanics and Mathematics, Moscow State University, Moscow, Russia Correspondence: Emmanuel Wyser ([email protected]) Abstract. We present an efficient MATLAB-based implementation of the material point method (MPM) and its most recent variants. MPM has gained popularity over the last decade, especially for problems in solid mechanics in which large defor- mations are involved, i.e., cantilever beam problems, granular collapses, and even large-scale snow avalanches. Although its numerical accuracy is lower than that of the widely accepted finite element method (FEM), MPM has been proven useful 5 in overcoming some of the limitations of FEM, such as excessive mesh distortions. We demonstrate that MATLAB is an ef- ficient high-level language for MPM implementations that solve elasto-dynamic and elasto-plastic problems. We accelerate the MATLAB-based implementation of MPM method by using the numerical techniques recently developed for FEM opti- mization in MATLAB. These techniques include vectorisation, the usage of native MATLAB functions, the maintenance of optimal RAM-to-cache communication, and others. We validate our in-house code with classical MPM benchmarks including 10 i) the elastic collapse under self-weight of a column, ii) the elastic cantilever beam problem, and iii) existing experimental and numerical results, i.e., granular collapses and slumping mechanics respectively. We report a performance gain by a factor of 28 for a vectorised code compared to a classical iterative version.
    [Show full text]
  • Curriculum Vitae Boyce E
    Curriculum Vitae Boyce E. Griffith Assistant Professor of Mathematics Adjunct Assistant Professor of Biomedical Engineering Member, McAllister Heart Institute The University of North Carolina at Chapel Hill Contact Information Department of Mathematics Phillips Hall, Campus Box 3250 University of North Carolina Chapel Hill, NC 27599-3250 email: [email protected] web: http://griffith.web.unc.edu phone: 919.962.7110 fax: 919.962.2568 Research Interests Mathematical and computational methods in medicine and biology; computer simulation in physiology, especially cardiovascular mechanics and fluid-structure interaction, cardiac electrophysiology, and cardiac electro-mechanical coupling; multiscale physiological models; medical devices; adaptive numerical methods; and scientific computing Education Courant Institute of Mathematical Sciences, New York University, New York, New York Ph.D. in Mathematics, September 2005 Rice University, Houston, Texas B.A. in Computational and Applied Mathematics and in Mathematics, May 2000 B.S. in Computer Science, May 2000 Professional Experience Employment 2014–present Assistant Professor of Mathematics (tenure track), University of North Carolina at Chapel Hill 2012–2014 Assistant Professor of Medicine and Mathematics (tenure track), New York University 2008–2012 Assistant Professor of Medicine (tenure track), New York University 2006–2008 American Heart Association Postdoctoral Research Fellow, Courant Institute of Mathematical Sciences, New York University 2005–2006 Courant Instructor, Department of Mathematics,
    [Show full text]
  • Application of Monte Carlo Method in Grid Computing and Allied Fields
    International Journal of Advanced Research in Computer Engineering & Technology (IJARCET) Volume 3, Issue 2, February 2014 Application of Monte Carlo method in Grid Computing and Allied Fields Manas kumar Yogi, Vijaya kranthi chinthala Abstract— This paper presents the strategy involved in Monte Carlo Method and how it impacts the functionality of many dependable application domains like Stock Markets, (1) Nuclear reactors, Gaming Industry, Biological fields. The huge benefits outweigh the small number of challenges posed while This approximation tends towards the exact value of the area simulation with this method. Monte Carlo methods are of the circle as the number of points N tends to infinity, i.e., generally useful for tools to model empirical studies where complex input nature demands precision handling. Index Terms—Monte Carlo, Random ,Generator, (2) 2 simulation, Grid. The value of π is then obtained with the relation π = Acircle/r . I. INTRODUCTION The Monte Carlo method is intensively used in various II. MECHANISM FOR APPLYING MONTE CARLO METHOD IN areas of scientific research. From computational physics to GRID COMPUTING fluid dynamics, it has seen exponential growth with the A. Random Number Generator introduction of computationally powerful computers. Truly, the Monte Carlo method is of great interest for solving The Monte Carlo method uses the stochastic landing of systems with unknown analytical solutions. In the real world, points in the square to find the area of the circle. This area more often than not, straight forward analytical solutions are tends towards the exact value of the area of the circle as the number of points tends towards infinity.
    [Show full text]
  • Stochastic Computational Mechanics
    AMAS COURSE ON RELIABILITY-BASED DESIGN AND OPTIMISATION RB0'03- (PP.99-156) -WARSAW, SEPTEMBER 15-18, 2003 Stochastic computational mechanics C. BUCHER and M. MACKE Bauhaus-University Weimar Institute of Structural Mechanics Marienstrasse. 15, D-99421 Weimar, Germany christian.bucher~bauing.uni-weimar.de Recent methods in computational structural mechanics put more emphasis on the consideration of uncertainties of loads and structural properties within the analysis. In most cases, these methods rely heavily on probability theory and statistics. Computational procedures for an efficient handling of random phe­ nomena in structural mechanics have been developed over the last few decades. This includes specific concepts for stochastic finite elements which are based on the possibly very close relation between structural and statistical/stochastic data. The paper starts with a description of essential elements from probability the­ ory. This includes random variables, random vectors, and random fields. Special emphasis is then put on outlining methods to determine structural failure prob­ abilities under static and dynamic loading conditions. These methods cover ac­ curate methods such as Monte Carlo simulation as well as approximate methods such as the first order reliability method or the response surface method. Several numerical examples from different areas of structural mechanics demonstrate the application of these methods. Key words: Reliability analysis, nonlinear structu.rol analysis, stochastic mechan­ ics, rondom fields, computational mechanics, stochastic finite elements. 1. Introduction Numerical methods for structural analysis have been developed quite sub­ stantially over the last decades. In particular, finite element methods and closely related approximations have become state of the art. The modeling capabilities and the solution possibilities lead, on the one hand, to an in­ creasing refinement allowing for more and more details to be captured in the analysis.
    [Show full text]
  • Quantum Mechanics Computational Physics
    Quantum Mechanics_Computational physics Computational physics is the study and implementation of Numerical analysis to solve problems in physics for which a quantitative theory already exists.[1] Historically, computational physics was the first application of modern computers in science, and is now a subset of computational science.[2] It is sometimes regarded as a subdiscipline (or offshoot) of theoretical physics, but others consider it an intermediate branch between theoretical and experimental physics, a third way that supplements theory and experiment. Overview In physics, different theories based on mathematical models provide very precise predictions on how systems behave. Unfortunately, it is often the case that solving the mathematical model for a particular system in order to produce a useful prediction is not feasible. This can occur, for instance, when the solution does not have a closed- form expression, or is too complicated. In such cases, numerical approximations are required. Computational physics is the subject that deals with these numerical approximations: the approximation of the solution is written as a finite (and typically large) number of simple mathematical operations (algorithm), and a computer is used to perform these operations and compute an approximated solution and respective error.[1] Challenges in computational physics Physics problems are in general very difficult to solve exactly. This is due to several (mathematical) reasons: lack of algebraic and/or analytic solubility, complexity and chaos. For example - even apparently simple problems, such as calculating the wavefunction of an electron orbiting an atom in a strong electric field (Stark effect), may require great effort to formulate a practical algorithm (if one can be found); other cruder or brute-force techniques, such as graphical methods or root finding, may be required.
    [Show full text]
  • The Multiscale Damage Mechanics in Objected-Oriented Fortran Framework
    The Multiscale Damage Mechanics in Objected-oriented Fortran Framework Zifeng Yuan Submitted in Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy in the Graduate School of Arts and Sciences Columbia University 2016 ©2016 Zifeng Yuan All Rights Reserved ABSTRACT The Multiscale Damage Mechanics in Objected-oriented Fortran Framework Zifeng Yuan We develop a dual-purpose damage model (DPDM) that can simultaneously model intralayer damage (ply failure) and interlayer damage (delamination) as an alternative to conventional practices that models ply failure by continuum damage mechanics (CDM) and delamination by cohesive elements. From purely computational point of view, if successful, the proposed approach will significantly reduce computational cost by eliminating the need for having double nodes at ply interfaces. At the core, DPDM is based on the regularized continuum damage mechanics approach with vectorial representation of damage and ellipsoidal damage surface. Shear correction factors are introduced to match the mixed mode fracture toughness of an analytical cohesive zone model. A predictor-corrector local-nonlocal regularization scheme, which treats intralayer portion of damage as nonlocal and interlayer damage as local, is developed and verified. Two variants of the DPDM are studied: a single- and two- scale DPDM. For the two- scale DPDM, reduced-order-homogenization (ROH) framework is employed with matrix phase modeled by the DPDM while the inclusion phase modeled by the CDM. The proposed DPDM is verified on several multi-layer laminates with various ply orientations including double-cantilever beam (DCB), end-notch-flexure (ENF), mixed-mode-bending (MMB), and three-point-bending (TPB). The simulation is executed in the platform of FOOF (Finite element solver based on Object-Oriented Fortran).
    [Show full text]
  • Discontinuous Galerkin Methods for Computational Fluid Dynamics
    Discontinuous Galerkin Methods for Computational Fluid Dynamics B. Cockburn School of Mathematics, University of Minnesota, Minneapolis, Minnesota 55455, U.S.A. ABSTRACT The discontinuous Galerkin methods are locally conservative, high-order accurate, and robust methods that can easily handle elements of arbitrary shapes, irregular triangulations with hanging nodes, and polynomial approximations of different degrees in different elements. These properties, which render them ideal for hp-adaptivity in domains of complex geometry, have brought them to the main stream of computational fluid dynamics. In this paper, we study the properties of the DG methods as applied to awidevarietyofproblemsincludinglinear,symmetrichyperbolicsystems,theEulerequationsofgas dynamics, purely elliptic problems and the incompressible and compressible Navier-Stokes equations. In each instance, we discuss the main properties of the methods, display the mechanisms that make them work so well, and present numerical experiments showing their performance. key words: computational fluid dynamics, discontinuous finite elements 1. Introduction This paper is a short essay on discontinuous Galerkin (DG) methods for fluid dynamics. The DG methods provide discontinuous approximations defined by using a Galerkin method element by element, the connection between the values of the approximation in different elements being established by the so-called numerical traces.Sincethemethodsuse discontinuous approximations, they can easily handle elements of arbitrary shapes, irregular triangulations with hanging nodes, and polynomial approximations of different degrees in different elements. The methods are thus ideally suited for hp-adaptivity in domains of complex geometry. Moreover, since they use a Galerkin method on each element, they can easily achieve high-order accuracy when the exact solution is smooth and high-resolution when it is not.
    [Show full text]
  • D7.5 Application Performance on Accelerators
    E-Infrastructures H2020-EINFRA-2014-2015 EINFRA-4-2014: Pan-European High Performance Computing Infrastructure and Services PRACE-4IP PRACE Fourth Implementation Phase Project Grant Agreement Number: EINFRA-653838 D7.5 Application performance on accelerators Final Version: 1.0 Author(s): Victor Cameo Ponz, CINES Date: 24.03.2017 D7.5 Application performance on accelerators Project and Deliverable Information Sheet PRACE Project Project Ref. №: EINFRA-653838 Project Title: PRACE Fourth Implementation Phase Project Project Web Site: http://www.prace-project.eu Deliverable ID: < D7.5> Deliverable Nature: <DOC_TYPE: Report / Other> Dissemination Level: Contractual Date of Delivery: PU 31 / 03 / 2017 Actual Date of Delivery: 31 / 03 / 2017 EC Project Officer: Leonardo Flores Añover * - The dissemination level are indicated as follows: PU – Public, CO – Confidential, only for members of the consortium (including the Commission Services) CL – Classified, as referred to in Commission Decision 2991/844/EC. Document Control Sheet Title: Application performance on accelerators Document ID: D7.5 Version: <1.0> Status: Final Available at: http://www.prace-project.eu Software Tool: Microsoft Word 2010 File(s): D7.5 Written by: Victor Cameo Ponz, CINES Authorship Contributors: Adem Tekin, ITU Alan Grey, EPCC Andrew Emerson, CINECA Andrew Sunderland, STFC Arno Proeme, EPCC Charles Moulinec, STFC Dimitris Dellis, GRNET Fiona Reid, EPCC Gabriel Hautreux, INRIA Jacob Finkenrath, CyI James Clark, STFC Janko Strassburg, BSC Jorge Rodriguez, BSC Martti Louhivuori,
    [Show full text]
  • COMPUTATIONAL MECHANICS (2º Cyy,Cle, Master of Science Degree)
    LECTURE (CU) COMPUTATIONAL MECHANICS (2º cyy,cle, Master of Science Degree) Paulo Piloto Applied Mechanics Department School of Technology and Management Polytechnic Institute of Braganza CURRICULAR UNIT ‐ SYLLABUS • Subject: Computational Mechanics. • Course Speciality: Construction and Industrial Engineering. • Main Scientific area: Structures and Solid Mechanics (MSE). • Classes: 60 h/Semester: T (Theoretical) PL (Practice and Laboratory). • Cycle: 2º (Master degree of Science). • Year / Semester: 1º year/ 2nd Semester. • Learning outcomes and competences: – Understand and apply Finite Element Method formulation. – To be aware of beam and bar finite element formulation. – Understand and appl y two and three dimensional elasticity formulation. – Understand and apply plate and shell finite element formulation. – Understand and apply the finite element method and the numerical solutions. – Learn to use commercial finite element software. CURRICULAR UNIT ‐ SYLLABUS • Cour se conte nts ((ededExtended veesorsion) : – Chapter 1 ‐ Stages of the FEM. Bar finite element: • Introduction, advantages and applications of the finite element method (FEM). Basic concepts in matrix analysis of structures. Types of analysis. Fundamental steps in the FEM. Phases of the method. Mathematical model formulation. Discrete mathematical models. Static and dynamic formulations. Stiffness matrix and element assembly. Continuous mathematical models. Variational formulation. Bar element formulation. Matrix formulation of the element equations. Isoparametric formulation and numerical integration. – Chapter 2 ‐ Finite element formulation: • Standard flowchart of a finite element code. FEM general methodology. Shape functions. Interpolation of displacements. Displacement and strain fields. Stress field. Constitutive models. Solution of the FEM equations. FEM convergence requirements and error types. OiOptima lpoints for stress callilculations. – Chapter 3 – Beam finite elements: • Euler‐Bernoulli beam finite element.
    [Show full text]
  • Eighty Years of the Finite Element Method: Birth, Evolution, and Future
    Eighty Years of the Finite Element Method: Birth, Evolution, and Future Wing Kam Liu, Northwestern University Shaofan Li, UC Berkeley Harold S. Park, Boston University This year marks the eightieth anniversary of the invention of the finite element method (FEM). FEM has become the computational workhorse for engineering design analysis and scientific modeling of a wide range of physical processes, including material and structural mechanics, fluid flow and heat conduction, various biological processes for medical diagnosis and surgery planning, electromagnetics and semi-conductor circuit and chip design and analysis, additive manufacturing, i.e. virtually every conceivable problem that can be described by partial differential equations (PDEs). FEM has fundamentally revolutionized the way we do scientific modeling and engineering design, ranging from automobiles, aircraft, marine structures, bridges, highways, and high-rise buildings. Associated with the development of finite element methods has been the concurrent development of an engineering science discipline called computational mechanics, or computational science and engineering. In this paper, we present a historical perspective on the developments of finite element methods mainly focusing on its applications and related developments in solid and structural mechanics, with limited discussions to other fields in which it has made significant impact, such as fluid mechanics, heat transfer, and fluid-structure interaction. To have a complete storyline, we divide the development of the finite element method into four time periods: I. (1941-1965) Early years of FEM; II. (1966-1991) Golden age of FEM; III. (1992-2017) Large scale, industrial applications of FEM and development of material modeling, and IV (2018-) the state-of-the-art FEM technology for the current and future eras of FEM research.
    [Show full text]