ON THE IMPLEMENTATION OF THE P-ADAPTIVE FINITE ELEMENT METHOD USING THE OBJECT
|
|
- Shawn Evans
- 6 years ago
- Views:
Transcription
1 ON THE IMPLEMENTATION OF THE P-ADAPTIVE FINITE ELEMENT METHOD USING THE OBJECT ORIENTED PROGRAMMING PHILOSOFY Philippe R.B. Devloo Faculdade de Engenharia Civil - UNICAMP Carlos Armando Magalhaes Duarte & Andre Tamagnini Noel Departamento de Engenharia Mec^anica - UFSC phil@fec.unicamp.br SUMMARY In this paper, an implementation of the p-adaptive nite element method is described using the object oriented programming philosofy. It is shown that, using the object oriented programming language C++, the p-adaptive method can be implemented with relative ease. The structure of the resulting program permits the same nite element class to be applied to a variety of computational mechanics problems. The blockoriented banded matrix implementation allows to increase and/or decrease the polynomial order of a given element with high eciency. The concepts oered by the object oriented programming philosofy can be used to extend the nite element program to include new elds of computational mechanics and implement dierent algorithms such as h-p adaptivity, multigrid, etc. The authors hope that their eort will form the basis of future scientic cooperation projects about the development of adaptive nite element software. 1 INTRODUCTION In recent years, the object oriented programming philosofy has established itself as the programming methodology of choice for such areas as operating systems, graphical user interfaces, graphics programming and databases. ([1]). In this work, an object oriented implementation of a nite element program is presented, bringing the benets of object oriented programming to scientic computing. Rather than proposing an object oriented implementation of the traditional nite element data structures, a new structure is developed which is particularly suited for implementing adaptive nite element algorithms. In a parallel eort, an object oriented system-independent graphics environment is being developed which facilitates the development of interactive programs. In the near future, both eorts will be combined to form a p-adaptive nite element program with a graphical interface. The Program Development Cycle During the development of an object oriented program, it is common to work with both a bottom up and top down design. With top down design, we denote the description of the global organization of the program. At this level, the interaction of the dierent classes is described without considering the details of their implementation. In the bottom up design, the ner granularity of the
2 program is discussed. In a typical nite element program, this concerns the implementation of an element, a node, an integration rule, etc. These individual modules then interact to form a complete program. Both aproaches are usually combined to form a functional program. In the rst section, a global description of the nite element program is given. In the second and further sections, individual building blocks are described, motivating the choices which were made when dening their behaviour. 2 GENERAL ORGANIZATION The main purpose of a nite element program is to perform a nite element analysis. Therefore, the top level class in the program is the analysis class. The analysis class combines a datale, a geometric grid, a computable grid, a matrix and a solution vector in order to perform the global steps of a nite element analysis : data initialization, matrix assembly and decomposition. Using these four objects, the analysis class reads a datale and builds the geometric and computable grid. The matrix object stored in the analysis class allows for a choice of various equation solvers. The geometric grid contains the data items needed to describe the geometry of the domain of the problem. As such the geometric grid functions as a database of geometric nodes, geometric elements and geometric descriptions of boundary conditions (elements and sides, nodes, etc). Putting the information in a database like structure facilitates the use of out-of-core memory schemes and save/restore operations. The computable grid contains the data items necessary to dene the dierential equation, its boundary conditions and the nite element interpolation. These items are dened through lists of materials, computable elements, computable nodes and boundary conditions. Like the geometric grid, the computable grid functions as a database of these items rather than as an object which does the actual computations. The nite element operations such a computation of the stiness matrix, assembly and decomposition of the global matrix are done by the analysis class. The matrix class consists of a family of classes which represent the global stiness matrix of the nite element approximation. Depending on the class being used, the matrix object uses a banded, block banded or element by element storage scheme. The implementation of the block banded storage scheme will be discussed in a later section of this paper. Finally, the datale consists of a family of classes, each being specialized in reading a particular datale format. In this sense, classes derived from the base class datale can be developed to read NASTRAN input les, ANSYS les or any format corresponding to a traditional nite element program. It should be noted that, at this level, only the global functionality of the program is described, without mentioning the data structures which will sustain its implementation. This is typical for object oriented software, which is based on the interaction of objects rather than on the judicious choice of data structures.
3 3 ELEMENT AND NODE IMPLEMENTATION As mentioned before, in the proposed program structure, a separatation is made between the geometric description of the domain and the description of the differential equation, its boundary conditions and interpolation scheme. At the element and node level, this separation is evidenced by the existence of geometric and computable elements and geometric and computable nodes. The geometric element and node are responsible for the geometric description of a nite element. Its only responsibility towards the nite element computations is the generation of the jacobian matrix at a given integration point. This approach allows for the easy implementation of blending functions within the given program structure and for the dierential use of interpolation functions for geometric modeling and interpolation. Blending functions are particularly suitable for mapping analytically described curved boundary segments onto a square master region, a feature which is essential for the correct implementation of the p-adaptive method. The dierential use of families of shapefunctions for geometric modeling and interpolation is important because the geometric mapping best uses Lagrangian type shapefunctions whereas the better families for interpolation functions are hierarchical, such as Chebyshev or Legendre polynomials. ([6], [5]). The geometric node has an id, a coordinate system and the coordinates of the point associated with it. Its main function is to be used as a reference by the geometric elements. The computable element denes the element interpolation and implements the computation of the stiness matrix and boundary condition contribution. In its current implementation, the quadrilateral element implements arbitrary interpolation orders in the and directions of the master element. Corresponding to the interpolation orders the computable element contains two pointers to objects containing the Gauss integration points and weights. Each quadrilateral element has nine computable nodes, some of which may have zero degrees of freedom (e.g. in the case of linear interpolation). When an element is enriched in either of both directions, the number of degrees of freedom of the corresponding computable nodes is increased. The computable nodes have a unique id, an equation number and a vector of degrees of freedom. In the case of linear analysis, the values of the degrees of freedom are only used during post-processing. 4 ELEMENT INTEGRATION In almost all nite element programs, the element stiness is computed within an appropriate stiness routine. Each element implemented in the program has a separate stiness routine. Therefore, the completeness of a nite element program is often judged by the number of elements which are available. In the program presented, a dierent approach is taken, using the resources oered by the object oriented programming philosofy : a single stiness routine is written which is aplicable to any linear problem using a single integration rule. Such approach has the advantage that, once the p-adaptive algorithm has been implemented for one particular problem, it is available for all similar problems,
4 irrespective of the number of degrees of freedom per node. In a traditional nite element program, the material properties are passed on to the stiness routine as dierential equation coecients. These coecients are then used to integrate the stiness matrix and right hand side. In the approach presented, the stiness routine only performs the integration of a function returning matrix values. The function in this case is a method of the material class which, as a function of the shape function values and their derivatives, returns the corresponding matrix value. As the matrix class is an abstract concept independent of its size, the stiness routine presented is independent of the number of degrees of freedom per node and is therefore aplicable to any linear dierential equation. It is now irrelevant to ask how many elements are implemented. One should rather ask which material classes are available. Another advantage of this approach is that, whenever a new material is implemented, its approximations are automatically available in the p-adaptive context and, in the near future, the h-p-adaptive context as well. 5 BOUNDARY CONDITION IMPLEMENTA- TION In a traditional nite element program, the boundary conditions are invariably implemented on the global stiness matrix. Such approach is valid because all degrees of freedom have a physical meaning (e.g. displacement, temperature, pressure, etc.) and because the program performs only one type of simulation (heat transfer, elasticity, etc.). In p-adaptive methods, the shape functions are generally hierarchical and therefore the degrees of freedom lose their physical context. It is also impossible to apply Neumann or mixed boundary conditions at the level of the global stiness matrix. All boundary conditions (Dirichlet, Neumann and Mixed) are therefore implemented at the element level : the local stiness matrix and/or the right hand side are modied so that the global solution satises the given boundary conditions. The stiness matrix and right hand side are modied by integrating the proper coecients over the element sides. As in the stiness integration, the contributions to the stiness matrix are computed by the object of the material class, which, given the values of the shape functions and their derivatives on the boundary, computes the contribution to the element stiness matrix. This approach allows for the specication of boundary condition for elements with arbitrary number of shape functions and arbitrary sets of shape functions (e.g. Chebyshev, Legendre, etc.). It is also possible to apply this methodology for most problems in computational mechanics, linear or non-linear.
5 6 NODE NUMBERING, THE IMPLEMEN- TATION OF BLOCKED MATRICES As mentioned before, geometric nodes and computable nodes are implemented as separate classes. The computable nodes have a unique id, an equation number and a vector of degrees of freedom. Each computable node can have an arbitrary number of degrees of freedom, determined by the size of the vector of degrees of freedom. The notion of equation and equation number has been slightly extended in that an equation now refers to a group of algebraic equations corresponding to the degrees of freedom of a computable node. The concept of matrix was also extended to include blocked matrices. A blocked matrix can be viewed as rectangular matrix of submatrices. Each submatrix can have an arbitrary number of rows and columns, although it is assumed that the number of rows of the submatrices of a given row of the blocked matrix are identical and that the number of columns of the submatrices of a given column of a blocked matrix are also identical. The number of rows and/or columns of a submatrix can be zero, as is the case when the corresponding node has zero degrees of freedom. The data of the blocked matrix can be accessed in the traditional fashion as (i; j) denoting the element of the ith row and jth column, and also as (i; j; k; l) denoting the (k; l) element of submatrix (i; j). The blocked matrix includes a method to increase/decrease the number rows of the submatrices of a given row and/or increase/decrease the number columns of the submatrices of a given column. The underlying data structure of the blocked matrix is hidden from the user. In fact, two implementations were tested : a rst implementation in which the blocked matrix actually stores a matrix of objects of the class matrix and a second implementation in which the behaviour of the blocked matrix was implemented on top of the traditional storage scheme. The rst implementation is ecient for submatrix operations and is slower for indexing operations, whereas the second implementation is ecient for indexing and slower for submatrix operations. In the nite element program presented, the indexing operations clearly dominated and the second implementation was preferred. 7 THE IMPLEMENTATION OF THE P-RE- FINEMENT METHOD Given the structure of the program described in the previous sections, the p- renement algorithm is almost trivially implemented. As mentioned before, the quadrilateral computable elements have nine nodes, 4 nodes corresponding to the corners, 4 midside nodes and one node corresponding to the internal degrees of freedom. Each node can have an arbitrary number of degrees of freedom. The equations corresponding to a node are stored on a per node basis in a blocked matrix. When increasing/decreasing the polynomial order of either direction of the quadrilateral element, it is only needed to change the number of degrees of freedom of the corresponding nodes, without modifying the structure of the
6 blocked stiness matrix. Therefore, the equation (block) number of the computable nodes remains unchanged during a p-adaptive cycle, resulting in a very ecient adaptive code. This methodology is also trivially extended to triangles. The p-adaptive technique has been applied to several test cases described in [7] and yielded identical results as the examples mentioned in that work. No numerical results are presented here because the scope of this work is to describe the implementation of the p-adaptive element rather than to show numerical results. The real value of adaptive algorithms will only become aparent in the context of an interactive graphics environment. In a parallel eort, the rst author of this work is developing a graphics environment based on the object oriented programming philosofy with the perspective of incorporating the p-adaptive algorithm presented here. 8 CONCLUSIONS The object oriented programming philosofy enables the development of an innovative approach to nite element programming. The new structure presented is ideally suited for implementing the p-adaptive nite element method. It also allows for a single nite element program to be applied to virtually any linear simulation in the area of computational mechanics, irrespective of the number of degrees of freedom per node. It is believed that the data structure presented above takes advantage of the benets oered by the object oriented programming phylosofy and will form the basis of a powerful nite element program to be developed in the near future. References [1] Encarnac~ao J.L. and Lockemann, P.C. - Engeneering Databases. Springer Verlag, Berlin, [2] Forde, B. W. R., Foschi, R. O. and Stiemer, S. F. - Object Oriented Finite Element Analysis, Computers and Structures, 34, No 1, , (1990). [3] Alves Filho, J.S.R., Devloo, P.R.B. - Object Oriented Programming in Scientic Computing, The Beginning of a New Era, Engineering Computations, 8, (1991). [4] Devloo, P.R.B. and Alves Filho, J.S.R. - An Object Oriented Approach to Finite Element Programming (Phase 1) : A System Independent Windowing Environment for Developing Scientic Programs, Advances in Engineering Software and Workstations, to appear. [5] Devloo, P.R.B., Oden, J.T., Pattani, P. - An h-p Adaptive Finite Element Method for the Numerical Simulation of Compressible Flow, em Computer Methods in Applied Mechanics and Engineering, 70, , (1988). [6] Zienkiewicz, O.C., Craig, A. - Adaptive Renement, Error Estimates, Multigrid Solution, and Hierarchical Finite Element Method Concept, in Accuracy
7 Estimates and Adaptive Renements in Finite Element Computations, eds. Babuska, I., Zienkiewicz, O.C., Gago, J., Oliveira, E.R.de A., John Wiley & Sons, Chisester, [7] Duarte, C.A.M. Estudo da Vers~ao P do Metodo de Elementos Finitos para Problemas da Elasticidade e de Potencial, Master's Thesis, UFSC, Florianopolis, Brasil, November [8] Thomas, D. - What's in an object Byte, March, (1989). [9] Wiener, R.S., Pinson, L.J. An Introduction to Object Oriented Programming and C++. Addison Wesley, [10] Stroustrup, B. The C++ Programming Language. Addison Wesley, [11] Meyer, B. Object Oriented Software Construction. Prentice Hall, [12] CNS, Inc. C++/Views User Guide, software manual, 1991.
TICAM - Texas Institute for Computational and Applied Mathematics. The University of Texas at Austin. Taylor Hall Abstract
A New Cloud-Based hp Finite Element Method J. T. Oden, C. A. M. Duarte y and O. C. Zienkiewicz z TICAM - Texas Institute for Computational and Applied Mathematics The University of Texas at Austin Taylor
More informationJ. TINSLEY ODEN, and C. ARMANDO DUARTE y. Texas Inst. for Comput. and Appl. Math. The University of Texas at Austin. Taylor Hall 2.
CLOUDS, CRACKS AND FEM'S J. TINSLEY ODEN, and C. ARMANDO DUARTE y Texas Inst. for Comput. and Appl. Math. The University of Texas at Austin Taylor Hall 2.400 Austin, Texas, 78712, U.S.A. Dedicated to Professor
More informationDEVELOPMENT OF A PARAMETRIC PROGRAM FOR SIMULATION OF THE STAMPING PROCESS USING FINITE ELEMENTS METHOD
DEVELOPMENT OF A PARAMETRIC PROGRAM FOR SIMULATION OF THE STAMPING PROCESS USING FINITE ELEMENTS METHOD Erick Siqueira Guidi, eguidi@uol.com.br Fernando de Azevedo Silva, fazevedo@feg.unesp.br UNESP São
More informationContents. I The Basic Framework for Stationary Problems 1
page v Preface xiii I The Basic Framework for Stationary Problems 1 1 Some model PDEs 3 1.1 Laplace s equation; elliptic BVPs... 3 1.1.1 Physical experiments modeled by Laplace s equation... 5 1.2 Other
More informationJournal of Engineering Research and Studies E-ISSN
Journal of Engineering Research and Studies E-ISS 0976-79 Research Article SPECTRAL SOLUTIO OF STEADY STATE CODUCTIO I ARBITRARY QUADRILATERAL DOMAIS Alavani Chitra R 1*, Joshi Pallavi A 1, S Pavitran
More informationcourse outline basic principles of numerical analysis, intro FEM
idealization, equilibrium, solutions, interpretation of results types of numerical engineering problems continuous vs discrete systems direct stiffness approach differential & variational formulation introduction
More informationreasonable to store in a software implementation, it is likely to be a signicant burden in a low-cost hardware implementation. We describe in this pap
Storage-Ecient Finite Field Basis Conversion Burton S. Kaliski Jr. 1 and Yiqun Lisa Yin 2 RSA Laboratories 1 20 Crosby Drive, Bedford, MA 01730. burt@rsa.com 2 2955 Campus Drive, San Mateo, CA 94402. yiqun@rsa.com
More informationf ( 1,η ) 1,2 2,2 2,1
Chapter 5 Mesh Generation and Assembly 5.1 Introduction There are several reasons for the popularity of nite element methods. Large code segments can be implemented for a wide class of problems. The software
More informationSpot Weld. Sheet1 Sheet2. diameter
High Performance, Process Oriented, Weld Spot Approach Dipl.-Phys. Daniel Heiserer 1, Mladen Chargin 2, Dipl.-Ing. Juergen Siela 3 1 Abstract Simulation of automotive vehicles has become a more important
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 informationH (div) approximations based on hp-adaptive curved meshes using quarter point elements
Trabalho apresentado no CNMAC, Gramado - RS, 2016. Proceeding Series of the Brazilian Society of Computational and Applied Mathematics H (div) approximations based on hp-adaptive curved meshes using quarter
More information(x,y ) 3 (x,y ) (x,y ) φ 1
Chapter 4 Finite Element Approximation 4. Introduction Our goal in this chapter is the development of piecewise-polynomial approximations U of a two- or three-dimensional function u. For this purpose,
More informationROTATIONAL DEPENDENCE OF THE SUPERCONVERGENT PATCH RECOVERY AND ITS REMEDY FOR 4-NODE ISOPARAMETRIC QUADRILATERAL ELEMENTS
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng, 15, 493±499 (1999) ROTATIONAL DEPENDENCE OF THE SUPERCONVERGENT PATCH RECOVERY AND ITS REMEDY FOR 4-NODE ISOPARAMETRIC QUADRILATERAL
More informationAPPLICATION OF ALGORITHMS FOR AUTOMATIC GENERATION OF HEXAHEDRAL FINITE ELEMENT MESHES
MESTRADO EM ENGENHARIA MECÂNICA November 2014 APPLICATION OF ALGORITHMS FOR AUTOMATIC GENERATION OF HEXAHEDRAL FINITE ELEMENT MESHES Luís Miguel Rodrigues Reis Abstract. The accuracy of a finite element
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Lecture - 24
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Lecture - 24 So in today s class, we will look at quadrilateral elements; and we will
More informationA class of parallel multiple-front algorithms on subdomains
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2003; 56:1569 1592 (DOI: 10.1002/nme.627) A class of parallel multiple-front algorithms on subdomains A. Bose 1, G.
More informationThe driving motivation behind the design of the Janus framework is to provide application-oriented, easy-to-use and ecient abstractions for the above
Janus a C++ Template Library for Parallel Dynamic Mesh Applications Jens Gerlach, Mitsuhisa Sato, and Yutaka Ishikawa fjens,msato,ishikawag@trc.rwcp.or.jp Tsukuba Research Center of the Real World Computing
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 informationObject-oriented programming in boundary element methods using C11
Advances in Engineering Software 30 (1999) 127±132 Object-oriented programming in boundary element methods using C11 Wenqing Wang*, Xing Ji, Yuangong Wang Department of Engineering Mechanics and Technology,
More informationW 1. x x L p p
Introduction to Finite Element Methods Paul S. Heckbert Carnegie Mellon University Global Illumination Course SIGGRAPH '93 Most simulation problems cannot be solved symbolically, and must instead be solved
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 informationAMS527: Numerical Analysis II
AMS527: Numerical Analysis II A Brief Overview of Finite Element Methods Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao SUNY Stony Brook AMS527: Numerical Analysis II 1 / 25 Overview Basic concepts Mathematical
More information1.2 Numerical Solutions of Flow Problems
1.2 Numerical Solutions of Flow Problems DIFFERENTIAL EQUATIONS OF MOTION FOR A SIMPLIFIED FLOW PROBLEM Continuity equation for incompressible flow: 0 Momentum (Navier-Stokes) equations for a Newtonian
More informationCONSTRUCTIONS OF QUADRILATERAL MESHES: A COMPARATIVE STUDY
South Bohemia Mathematical Letters Volume 24, (2016), No. 1, 43-48. CONSTRUCTIONS OF QUADRILATERAL MESHES: A COMPARATIVE STUDY PETRA SURYNKOVÁ abstrakt. Polygonal meshes represent important geometric structures
More informationHowever, m pq is just an approximation of M pq. As it was pointed out by Lin [2], more precise approximation can be obtained by exact integration of t
FAST CALCULATION OF GEOMETRIC MOMENTS OF BINARY IMAGES Jan Flusser Institute of Information Theory and Automation Academy of Sciences of the Czech Republic Pod vodarenskou vez 4, 82 08 Prague 8, Czech
More informationAdaptive Surface Modeling Using a Quadtree of Quadratic Finite Elements
Adaptive Surface Modeling Using a Quadtree of Quadratic Finite Elements G. P. Nikishkov University of Aizu, Aizu-Wakamatsu 965-8580, Japan niki@u-aizu.ac.jp http://www.u-aizu.ac.jp/ niki Abstract. This
More informationCURVILINEAR MESH GENERATION IN 3D
CURVILINEAR MESH GENERATION IN 3D Saikat Dey, Robert M. O'Bara 2 and Mark S. Shephard 2 SFA Inc. / Naval Research Laboratory, Largo, MD., U.S.A., dey@cosmic.nrl.navy.mil 2 Scientific Computation Research
More informationAS 5850 Finite Element Analysis
AS 5850 Finite Element Analysis INTRODUCTION Prof. IIT Madras Goal of engineering computations Perform analysis and design of physical systems and processes subjected to imposed conditions (or loads) and
More informationSecond-order shape optimization of a steel bridge
Computer Aided Optimum Design of Structures 67 Second-order shape optimization of a steel bridge A.F.M. Azevedo, A. Adao da Fonseca Faculty of Engineering, University of Porto, Porto, Portugal Email: alvaro@fe.up.pt,
More informationLecture 3.2 Methods for Structured Mesh Generation
Lecture 3.2 Methods for Structured Mesh Generation 1 There are several methods to develop the structured meshes: Algebraic methods, Interpolation methods, and methods based on solving partial differential
More informationElement Quality Metrics for Higher-Order Bernstein Bézier Elements
Element Quality Metrics for Higher-Order Bernstein Bézier Elements Luke Engvall and John A. Evans Abstract In this note, we review the interpolation theory for curvilinear finite elements originally derived
More informationlinearize discretize Galerkin optimize sample
Fairing by Finite Dierence Methods Leif Kobbelt Abstract. We propose an ecient and exible scheme to fairly interpolate or approximate the vertices of a given triangular mesh. Instead of generating a piecewise
More informationParameterization of triangular meshes
Parameterization of triangular meshes Michael S. Floater November 10, 2009 Triangular meshes are often used to represent surfaces, at least initially, one reason being that meshes are relatively easy to
More informationDevelopment of a user-friendly interface for the creation of user elements
Development of a user-friendly interface for the creation of user elements Kinshuk German National Research Center for Information Technology GMD-FIT, Schloss Birlinghoven D-53754 Sankt Augustin, Germany
More informationAn Object Oriented Finite Element Library
An Object Oriented Finite Element Library Release 3.1.0 Rachid Touzani Laboratoire de Mathématiques Blaise Pascal Université Clermont Auvergne 63177 Aubière, France e-mail: Rachid.Touzani@univ-bpclermont.fr
More informationJDEP 384H: Numerical Methods in Business
Instructor: Thomas Shores Department of Mathematics Lecture 1, January 9, 2007 110 Kaufmann Center Outline 1 2 3 Solving Systems Matrix and Vector Algebra Welcome to Matlab! Rational File Management File
More information13.472J/1.128J/2.158J/16.940J COMPUTATIONAL GEOMETRY
13.472J/1.128J/2.158J/16.940J COMPUTATIONAL GEOMETRY Lecture 23 Dr. W. Cho Prof. N. M. Patrikalakis Copyright c 2003 Massachusetts Institute of Technology Contents 23 F.E. and B.E. Meshing Algorithms 2
More informationA complete methodology for the implementation of XFEM models. Team: - Carlos Hernán Villanueva - Kai Yu
A complete methodology for the implementation of XFEM models Team: - Carlos Hernán Villanueva - Kai Yu Problem Statement QUAD4 DISCONTINUITY Problem Statement Traditional FEM Mesh generates discrete representation
More informationA numerical grid and grid less (Mesh less) techniques for the solution of 2D Laplace equation
Available online at www.pelagiaresearchlibrary.com Advances in Applied Science Research, 2014, 5(1):150-155 ISSN: 0976-8610 CODEN (USA): AASRFC A numerical grid and grid less (Mesh less) techniques for
More informationA ow-condition-based interpolation nite element procedure for triangular grids
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2006; 51:673 699 Published online in Wiley InterScience (www.interscience.wiley.com).1246 A ow-condition-based interpolation
More informationΩ2
CACHE BASED MULTIGRID ON UNSTRUCTURED TWO DIMENSIONAL GRIDS CRAIG C. DOUGLAS, JONATHAN HU y, ULRICH R UDE z, AND MARCO BITTENCOURT x. Abstract. High speed cache memory is commonly used to address the disparity
More informationu 0+u 2 new boundary vertex
Combined Subdivision Schemes for the design of surfaces satisfying boundary conditions Adi Levin School of Mathematical Sciences, Tel-Aviv University, Tel-Aviv 69978, Israel. Email:fadilev@math.tau.ac.ilg
More information2) For any triangle edge not on the boundary, there is exactly one neighboring
Triangulating Trimmed NURBS Surfaces Chang Shu and Pierre Boulanger Abstract. This paper describes techniques for the piecewise linear approximation of trimmed NURBS surfaces. The problem, called surface
More informationA Boolean Expression. Reachability Analysis or Bisimulation. Equation Solver. Boolean. equations.
A Framework for Embedded Real-time System Design? Jin-Young Choi 1, Hee-Hwan Kwak 2, and Insup Lee 2 1 Department of Computer Science and Engineering, Korea Univerity choi@formal.korea.ac.kr 2 Department
More informationMeshing of flow and heat transfer problems
Meshing of flow and heat transfer problems Luyao Zou a, Zhe Li b, Qiqi Fu c and Lujie Sun d School of, Shandong University of science and technology, Shandong 266590, China. a zouluyaoxf@163.com, b 1214164853@qq.com,
More informationSeed Point. Agglomerated Points
AN ASSESSMENT OF LINEAR VERSUS NON-LINEAR MULTIGRID METHODS FOR UNSTRUCTURED MESH SOLVERS DIMITRI J. MAVRIPLIS Abstract. The relative performance of a non-linear FAS multigrid algorithm and an equivalent
More informationAalborg Universitet. Object-Oriented Structuring of Finite Elements Hededal, O. Publication date: 1994
Aalborg Universitet Object-Oriented Structuring of Finite Elements Hededal, O. Publication date: 1994 Document Version Publisher's PDF, also known as Version of record Link to publication from Aalborg
More informationImprovement of semi-analytical design sensitivities of non-linear structures using equilibrium relations
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2001; 50:2127 2142 Improvement of semi-analytical design sensitivities of non-linear structures using equilibrium relations
More informationStatic and Dynamic Analysis Of Reed Valves Using a Minicomputer Based Finite Element Systems
Purdue University Purdue e-pubs International Compressor Engineering Conference School of Mechanical Engineering 1980 Static and Dynamic Analysis Of Reed Valves Using a Minicomputer Based Finite Element
More informationFinite element algorithm with adaptive quadtree-octree mesh refinement
ANZIAM J. 46 (E) ppc15 C28, 2005 C15 Finite element algorithm with adaptive quadtree-octree mesh refinement G. P. Nikishkov (Received 18 October 2004; revised 24 January 2005) Abstract Certain difficulties
More informationInvestigation of a Robust Method for Connecting Dissimilar 3D Finite Element Models. David M. Trujillo 1. December 2005
Investigation of a Robust Method for Connecting Dissimilar 3D Finite Element Models by David M. Trujillo 1 December 2005 1 Consultant, TRUCOMP, Fountain Valley, California trucomp@earthlink.net Abstract
More informationCHAPTER 5 USE OF STL FILE FOR FINITE ELEMENT ANALYSIS
CHAPTER 5 USE OF STL FILE FOR FINITE ELEMENT ANALYSIS 5.1 Introduction: Most CAD software in the market can generate STL files, and these are generally used for prototyping and rendering purposes. These
More informationParallel hp-finite Element Simulations of 3D Resistivity Logging Instruments
Parallel hp-finite Element Simulations of 3D Resistivity Logging Instruments M. Paszyński 1,3, D. Pardo 1,2, L. Demkowicz 1, C. Torres-Verdin 2 1 Institute for Computational Engineering and Sciences 2
More informationJassim Katwan, Osama (2010) Discrete modelling of heat transfer. MSc(R) thesis, University of Glasgow.
Jassim Katwan, Osama (2010) Discrete modelling of heat transfer. MSc(R) thesis, University of Glasgow. http://theses.gla.ac.uk/1602/ Copyright and moral rights for this thesis are retained by the author
More informationOperator-Based Backward Motion Estimation. Aria Nosratinia and Michael T. Orchard. Beckman Institute for Advanced Science and Technology.
Operator-Based Backward Motion Estimation Aria Nosratinia and Michael T. Orchard Beckman Institute for Advanced Science and Technology University of Illinois at Urbana-Champaign 405 N. Mathews Ave., Urbana,
More informationAN OBJECT ORIENTED CLASS ORGANIZATION FOR DYNAMIC GEOMETRICALLY NON-LINEAR FEM ANALYSIS
CMNE/CILAMCE 2007 Porto, 13 a 15 de Junho, 2007 c APMTAC, Portugal 2007 AN OBJECT ORIENTED CLASS ORGANIZATION FOR DYNAMIC GEOMETRICALLY NON-LINEAR FEM ANALYSIS Flavio T. Fonseca, Roque L. S. Pitangueira
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 informationObject Oriented Finite Element Modeling
Object Oriented Finite Element Modeling Bořek Patzák Czech Technical University Faculty of Civil Engineering Department of Structural Mechanics Thákurova 7, 166 29 Prague, Czech Republic January 2, 2018
More informationReduced Order Modeling Techniques for Mesh Movement Strategies as Applied to Fluid Structure Interactions. Alfred E.J. Bogaers
Reduced Order Modeling Techniques for Mesh Movement Strategies as Applied to Fluid Structure Interactions Alfred E.J. Bogaers University of Pretoria Reduced order modeling techniques for mesh movement
More informationParameterization. Michael S. Floater. November 10, 2011
Parameterization Michael S. Floater November 10, 2011 Triangular meshes are often used to represent surfaces, at least initially, one reason being that meshes are relatively easy to generate from point
More informationFully-automated hex-dominant mesh generation with directionality control via packing rectangular solid cells
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2003; 57:2099 2129 (DOI: 10.1002/nme.754) Fully-automated hex-dominant mesh generation with directionality control
More informationMid-Year Report. Discontinuous Galerkin Euler Equation Solver. Friday, December 14, Andrey Andreyev. Advisor: Dr.
Mid-Year Report Discontinuous Galerkin Euler Equation Solver Friday, December 14, 2012 Andrey Andreyev Advisor: Dr. James Baeder Abstract: The focus of this effort is to produce a two dimensional inviscid,
More information(Refer Slide Time: 00:02:24 min)
CAD / CAM Prof. Dr. P. V. Madhusudhan Rao Department of Mechanical Engineering Indian Institute of Technology, Delhi Lecture No. # 9 Parametric Surfaces II So these days, we are discussing the subject
More informationTilings. Mark McClure. July 16, Self-similarity is a concept often associated with fractal geometry. There
Digraph Self-Similar Sets and Aperiodic Tilings Mark McClure July 16, 2001 1 Introduction Self-similarity is a concept often associated with fractal geometry. There are many interesting self-similar sets
More informationUNSTRUCTURED GRIDS ON NURBS SURFACES. The surface grid can be generated either in a parameter. surfaces. Generating grids in a parameter space is
UNSTRUCTURED GRIDS ON NURBS SURFACES Jamshid Samareh-Abolhassani 1 Abstract A simple and ecient computational method is presented for unstructured surface grid generation. This method is built upon an
More informationLab 2: Support Vector Machines
Articial neural networks, advanced course, 2D1433 Lab 2: Support Vector Machines March 13, 2007 1 Background Support vector machines, when used for classication, nd a hyperplane w, x + b = 0 that separates
More informationA Finite Element Method for Deformable Models
A Finite Element Method for Deformable Models Persephoni Karaolani, G.D. Sullivan, K.D. Baker & M.J. Baines Intelligent Systems Group, Department of Computer Science University of Reading, RG6 2AX, UK,
More informationS 1,3 S 3,3 S 3,2 S 1,2 S 1,1 S 3,1
New Approaches for Particle Tracing on Sparse Grids Christian Teitzel and Thomas Ertl Computer Graphics Group, University of Erlangen Am Weichselgarten 9, 91058 Erlangen, Germany fteitzel,ertlg@informatik.uni-erlangen.de
More informationRELIABILITY OF THE FEM CALCULATIONS OF THE FRACTURE MECHANICS PARAMETERS
International Conference on Economic Engineering and Manufacturing Systems Braşov, 26 27 November 2009 RELIABILITY OF THE FEM CALCULATIONS OF THE FRACTURE MECHANICS PARAMETERS Galina TODOROVA, Valentin
More informationStatic analysis of eolicblade through finite element method and OOP C++
International Conference on Control, Engineering & Information echnology (CEI 4) Proceedings Copyright IPCO-4 ISSN 56-568 Static analysis of eolicblade through finite element method and OOP C++ MateusDantas,
More informationOn Partitioning Dynamic Adaptive Grid Hierarchies. Manish Parashar and James C. Browne. University of Texas at Austin
On Partitioning Dynamic Adaptive Grid Hierarchies Manish Parashar and James C. Browne Department of Computer Sciences University of Texas at Austin fparashar, browneg@cs.utexas.edu (To be presented at
More informationFinite element solution of multi-scale transport problems using the least squares based bubble function enrichment
Finite element solution of multi-scale transport problems using the least squares based bubble function enrichment A. Yazdani a, V. Nassehi b1 a Cranfield University, School of Applied Sciences, Cranfield,
More informationFOR ALL GRID SIZES. Thor Gjesdal. Christian Michelsen Research A/S. N-5036 Fantoft, Norway SUMMARY
A CELL-CENTERED MULTIGRID ALGORITHM FOR ALL GRID SIZES Thor Gjesdal Christian Michelsen Research A/S N-5036 Fantoft, Norway SUMMARY Multigrid methods are optimal; that is, their rate of convergence is
More informationTAU mesh deformation. Thomas Gerhold
TAU mesh deformation Thomas Gerhold The parallel mesh deformation of the DLR TAU-Code Introduction Mesh deformation method & Parallelization Results & Applications Conclusion & Outlook Introduction CFD
More informationRefine boundary at resolution r. r+1 r. Update context information CI(r) based on CI(r-1) Classify at resolution r, based on CI(r), update CI(r)
Context Based Multiscale Classication of Images Jia Li Robert M. Gray EE Department EE Department Stanford Univ., CA 94305 Stanford Univ., CA 94305 jiali@isl.stanford.edu rmgray@stanford.edu Abstract This
More informationTopology Optimization and JuMP
Immense Potential and Challenges School of Engineering and Information Technology UNSW Canberra June 28, 2018 Introduction About Me First year PhD student at UNSW Canberra Multidisciplinary design optimization
More informationPartial Differential Equations
Simulation in Computer Graphics Partial Differential Equations Matthias Teschner Computer Science Department University of Freiburg Motivation various dynamic effects and physical processes are described
More informationExtra-High Speed Matrix Multiplication on the Cray-2. David H. Bailey. September 2, 1987
Extra-High Speed Matrix Multiplication on the Cray-2 David H. Bailey September 2, 1987 Ref: SIAM J. on Scientic and Statistical Computing, vol. 9, no. 3, (May 1988), pg. 603{607 Abstract The Cray-2 is
More informationConditional Branching is not Necessary for Universal Computation in von Neumann Computers Raul Rojas (University of Halle Department of Mathematics an
Conditional Branching is not Necessary for Universal Computation in von Neumann Computers Raul Rojas (University of Halle Department of Mathematics and Computer Science rojas@informatik.uni-halle.de) Abstract:
More informationA. Atamturk. G.L. Nemhauser. M.W.P. Savelsbergh. Georgia Institute of Technology. School of Industrial and Systems Engineering.
A Combined Lagrangian, Linear Programming and Implication Heuristic for Large-Scale Set Partitioning Problems 1 A. Atamturk G.L. Nemhauser M.W.P. Savelsbergh Georgia Institute of Technology School of Industrial
More informationFlow simulation. Frank Lohmeyer, Oliver Vornberger. University of Osnabruck, D Osnabruck.
To be published in: Notes on Numerical Fluid Mechanics, Vieweg 1994 Flow simulation with FEM on massively parallel systems Frank Lohmeyer, Oliver Vornberger Department of Mathematics and Computer Science
More informationCHAPTER 5 FINITE ELEMENT METHOD
CHAPTER 5 FINITE ELEMENT METHOD 5.1 Introduction to Finite Element Method Finite element analysis is a computer based numerical method to deduce engineering structures strength and behaviour. Its use can
More informationLS-DYNA 980 : Recent Developments, Application Areas and Validation Process of the Incompressible fluid solver (ICFD) in LS-DYNA.
12 th International LS-DYNA Users Conference FSI/ALE(1) LS-DYNA 980 : Recent Developments, Application Areas and Validation Process of the Incompressible fluid solver (ICFD) in LS-DYNA Part 1 Facundo Del
More informationAvailable from Deakin Research Online:
Deakin Research Online Deakin University s institutional research repository DDeakin Research Online Research Online This is the published version (version of record) of: Xiang, Yang and Zhou, Wanlei 2002,
More informationFEM Convergence Requirements
19 FEM Convergence Requirements IFEM Ch 19 Slide 1 Convergence Requirements for Finite Element Discretization Convergence: discrete (FEM) solution approaches the analytical (math model) solution in some
More informationParallel Clustering on a Unidirectional Ring. Gunter Rudolph 1. University of Dortmund, Department of Computer Science, LS XI, D{44221 Dortmund
Parallel Clustering on a Unidirectional Ring Gunter Rudolph 1 University of Dortmund, Department of Computer Science, LS XI, D{44221 Dortmund 1. Introduction Abstract. In this paper a parallel version
More informationIntroduction to Finite Element Method
Guest Lecture in Prodi Teknik Sipil Introduction to Finite Element Method Wong Foek Tjong, Ph.D. Petra Christian University Surabaya Lecture Outline 1. Overview of the FEM 2. Computational steps of the
More information10/26/ Solving Systems of Linear Equations Using Matrices. Objectives. Matrices
6.1 Solving Systems of Linear Equations Using Matrices Objectives Write the augmented matrix for a linear system. Perform matrix row operations. Use matrices and Gaussian elimination to solve systems.
More informationWei Shu and Min-You Wu. Abstract. partitioning patterns, and communication optimization to achieve a speedup.
Sparse Implementation of Revised Simplex Algorithms on Parallel Computers Wei Shu and Min-You Wu Abstract Parallelizing sparse simplex algorithms is one of the most challenging problems. Because of very
More informationCHAPTER 6 EXPERIMENTAL AND FINITE ELEMENT SIMULATION STUDIES OF SUPERPLASTIC BOX FORMING
113 CHAPTER 6 EXPERIMENTAL AND FINITE ELEMENT SIMULATION STUDIES OF SUPERPLASTIC BOX FORMING 6.1 INTRODUCTION Superplastic properties are exhibited only under a narrow range of strain rates. Hence, it
More informationAdaptive numerical methods
METRO MEtallurgical TRaining On-line Adaptive numerical methods Arkadiusz Nagórka CzUT Education and Culture Introduction Common steps of finite element computations consists of preprocessing - definition
More informationA new cloud-based hp finite element method
A new cloud-based hp finite element method J.T. Oden, C.A.M. Duarte, O.C. Zienkiewicz To cite this version: J.T. Oden, C.A.M. Duarte, O.C. Zienkiewicz. A new cloud-based hp finite element method. Computer
More informationA spectral boundary element method
Boundary Elements XXVII 165 A spectral boundary element method A. Calaon, R. Adey & J. Baynham Wessex Institute of Technology, Southampton, UK Abstract The Boundary Element Method (BEM) is not local and
More informationEFFICIENT SOLVER FOR LINEAR ALGEBRAIC EQUATIONS ON PARALLEL ARCHITECTURE USING MPI
EFFICIENT SOLVER FOR LINEAR ALGEBRAIC EQUATIONS ON PARALLEL ARCHITECTURE USING MPI 1 Akshay N. Panajwar, 2 Prof.M.A.Shah Department of Computer Science and Engineering, Walchand College of Engineering,
More informationN. Hitschfeld. Blanco Encalada 2120, Santiago, CHILE.
Generalization of modied octrees for geometric modeling N. Hitschfeld Dpto. Ciencias de la Computacion, Univ. de Chile Blanco Encalada 2120, Santiago, CHILE E-mail: nancy@dcc.uchile.cl Abstract. This paper
More informationFinite Element Implementation
Chapter 8 Finite Element Implementation 8.1 Elements Elements andconditions are the main extension points of Kratos. New formulations can be introduced into Kratos by implementing a new Element and its
More informationWhat is Multigrid? They have been extended to solve a wide variety of other problems, linear and nonlinear.
AMSC 600/CMSC 760 Fall 2007 Solution of Sparse Linear Systems Multigrid, Part 1 Dianne P. O Leary c 2006, 2007 What is Multigrid? Originally, multigrid algorithms were proposed as an iterative method to
More informationFinite Element Model of Fracture Formation on Growing Surfaces
Finite Element Model of Fracture Formation on Growing Surfaces Pavol Federl and Przemyslaw Prusinkiewicz Department of Computer Science, University of Calgary Calgary, Alberta, Canada T2N 1N4 e mail: federl
More informationCHAPTER 1. Introduction
ME 475: Computer-Aided Design of Structures 1-1 CHAPTER 1 Introduction 1.1 Analysis versus Design 1.2 Basic Steps in Analysis 1.3 What is the Finite Element Method? 1.4 Geometrical Representation, Discretization
More informationADVANCES IN ADAPTIVE METHODS IN COMPUTATIONAL FLUID MECHANICS. J. Tinsley Oden. Abstract
ADVANCES N ADAPTVE METHODS N COMPUTATONAL FLUD MECHANCS J. Tinsley Oden Texas nstitute for Computational Mechanics The University of Texas at Austin Austin, Texas 78712 Abstract Recent developments in
More information