NEXT-GENERATION SWEEP TOOL: A METHOD FOR GENERATING ALL-HEX MESHES ON TWO-AND-ONE-HALF DIMENSIONAL GEOMTRIES
|
|
- Ralph Davis
- 5 years ago
- Views:
Transcription
1 NEXT-GENERATION SWEEP TOOL: A METHOD FOR GENERATING ALL-HEX MESHES ON TWO-AND-ONE-HALF DIMENSIONAL GEOMTRIES PATRICK M. KNUPP PARALLEL COMPUTING SCIENCES DEPARTMENT SANDIA NATIONAL LABORATORIES M/S 0441, P.O. BOX 5800 ALBUQUERQUE, NM PKNUPP@SANDIA.GOV Abstract. Placement of interior node points is a crucial step in the generation of quality meshes in sweeping algorithms. Two new algorithms were devised for node point placement and implemented in Sweep Tool, the first based on the use of linear transformations between bounding node loops and the second based on smoothing. Examples are given that demonstrate the effectiveness of these algorithms. 1. Introduction It has been demonstrated that geometries that are two-and-one-half dimensional (e.g., generalized cylinders) can be meshed with all-hexahedral finite elements [1], [4]. Because all-hexahedral mesh generation on general three-dimensional geometries remains an elusive goal, algorithms to mesh two-and-one-half dimensional geometries, generally referred to as sweeping or projecting methods, continue to be important. Although mesh sweeping is simple in concept, the first-generation algorithms lacked robustness and often produced poor mesh quality. There are several approaches to mesh generation via sweeping but common to all is the idea of identifying surfaces on a volume to serve as sources or caps and a complemetary set to serve as linking sides. Source surfaces may be arbitrarily meshed with quadrilaterals via paving or multi-block structured meshing and then swept along the linking sides towards a target surface which may or may not be pre-meshed. This is feasible provided the linking side surfaces are meshed with a type of multi-block, structured mesh known as submap [5]. Identification of these source and linking surfaces has been automated in the CUBIT code [6]. Interior mesh connectivity is entirely determined once these surfaces have been identified and meshed. Placement of the interior nodes in space is then a crucial step in generating a quality mesh. For geometries that are sweepable via source mesh translation or rotation, the spatial location of the interior nodes is unambiguous. For more general geometries, however, there is no uniquely correct location for the nodes. The guiding principle in the general case is that nodes are placed so that hexahredral elements of a swept mesh are of good quality. As a minimal requirement, the Jacobian at the eight corners of all hexahedrons should be positive [2]. Another reasonable requirement is that qualities of the source surface mesh such as biasing or the relative areas of quadrilaterals should be transferred, in so far as possible, to the interior layer meshes and to the target. Even with these requirements, flexibility remains in the placement of the interior nodes in sweeping general geometries. This paper reports on a robust method for uniquely specifying the interior node positions in a way that directly addresses the mesh quality issue. SANDIA IS A MULTIPROGRAM LABORATORY OPERATED BY SANDIA CORPORATION, A LOCKHEED MARTIN COMPANY, FOR THE UNITED STATES DEPARTMENT OF ENERGY UNDER CONTRACT DE-ACO4-94AL
2 The outline of this paper is as follows: section 2 describes how point placement was done in the first-generation Project Tool [4] in order to illustrate the challenge posed by general two-and-one-half dimensional geometries. Section 3 treats the point placement algorithm based on linear transformations that is used in Sweep Tool while Section 4 describes results of layer smoothing as a point placement method, both in Sweep and Project Tools. Of course, interior point placement is not the only issue that arises in designing a good sweeping algorithm, so Section 5 briefly covers some additional topics of importance. In section 6 we offer our conclusions. 2. Interior Point Placement in the Project Scheme Interior points are located in the Project scheme by an advancing front mechanism that builds a layer of hexahedrons on the front from the previous layer, begining with the source surface meshes. An individual hexahedron on the front is built using the four nodes x 0, x 1, x 2, and x 3 on the previous layer plus three other nodes x 4, x 5, and x 6 belonging to the hexahedron under construction. The latter are known from the boundary data on the linking surface and from pre-existing nodes on the layer containing the advancing front. In the original algorithm, the final node x 7 of the hexahedron was computed by averaging the three vectors x 4 x 0, x 5 x 1, and x 6 x 2 that approximately point in the direction in which the front advances: (1) x 7 = x {(x 4 x 0 ) + (x 5 x 1 ) + (x 6 x 2 )}. This approach gives the desired nodal position under a translation: if x 4+k = x k +c for some constant vector c R 3 and k = 0,1,2, it immediately follows that x 7 = x 3 + c. However, not all meshes on sweepable geometries can be translated so equation (1) is inadequate for many important cases. Foremost among these are rotatable geometries; rotations are not preserved under (1) and very poor meshes result. An improved planar-intersection algorithm was suggested in [4] in which the point x 7 was determined by the intersection of the three planes containing the points (x 4,x 5,x 6 ), (x 0,x 3,x 4 ), and (x 2,x 3,x 6 ), respectively. While this alternate algorithm improved the quality of the projected meshes, it still did not work effectively on rotatable geometries. Two approaches to solve this difficulty were suggested in [4]. In the first approach, a separate program module (known as scheme Rotate) was developed for geometries whose source meshes can be rotated. The planar-intersection algorithm was discarded in favor of one that works only in the rotational case. The result was a highly efficient module that consistently achieved good quality meshes on a limited set of geometries. A minor disadvantage of adding Rotate to the suite of tools was that there were then two modules to maintain when only one is actually needed. More importantly, however, is that other commonly encountered geometries such as a frustum needed development of additional special modules because the heuristically-based planarintersection algorithm failed to give quality meshes on them as well. To overcome this piecemeal approach, the authors of Project tried mesh smoothing. In some respects this was premature because, as is demonstrated in the next section, the planar-intersection algorithm can be replaced with a unifying technique involving linear transformations that handles all geometries which require translation, rotation, and scaling of layer meshes.
3 3. Interior Point Placement via Linear Transformations Given a layer mesh of points p j R 3 with bounding loop L consisting of the points {x k }, k = 1,2,...,K with K 3, and given a second loop L consisting of K points { x k } in R 3, we seek a vector b R 3 and a 3 3 non-singular linear transformation T between loops such that (2) x k = T x k + b for all k. Then, for any point p j belonging to the layer bounded by loop L, we determine the point p j in the layer bounded by loop L by p j = T p j + b. The vector b can be readily determined from the loop data. Define loop center points c = 1 (3) x k, K (4) c = 1 K k x k. It is not necessary that these center points lie inside the loops. If we sum (2) over all k, we find: (5) thus (2) becomes (6) k b = c T c ; x k c = T (x k c). For convennience, let u k = x k c and ũ k = x k c so that ũ k = T u k. Define K matrices U k = [u k,u k+1,u k+2 ] and Ũk = [ũ k,ũ k+1,ũ k+2 ]. 1 The former matrices will be non-singular provided the three column vectors are linearly independent, i.e., each portion of the loop must be non-planar. If U k is non-singular, then we can uniquely define K linear transformations T k = Ũk U 1 k which satisfy ũ k = T k u k. In general, however, there may not be a single linear transformation T between the loops. If T exists, then T k = T for all k. Since, in general, such a transformation between arbitrary loops does not exist, we perform a least-squares fit to the bounding loop data by minimizing the non-negative function F(T ) = 1 K (7) ũ k T u k 2. 2 k=1 Clearly, if T in (2) exists, then it mimimizes F. To find T, set F/ T ij = 0 to get three uncoupled linear systems that can be written as (8) with (9) M T M = F = k (u k u k ), (10) F = k (u k ũ k ), 1 Since the loops are closed we can achieve periodicity by letting x k+k = x k and x k+k = x k.
4 where is the vector outer product. T is determined provided M is non-singular but, unfortunately, M is singular for the important case of planar loops. Definition: a loop L is planar if there exists n R 3 such that n u k = 0 for all k. Suppose that n exists for loop L. Then Mn = (u k n)u k = 0. Hence the null-space of M is non-trivial and so M is singular. To circumvent the problem of planar loops, redefine the set of vectors u k, ũ k by (11) (12) u k = x k (2c c), ũ k = x k c. It is easy to show that u k = ũ k = K ( c c). Note that u k = (x k c) + ( c c), giving a clear geometric interpretation. Note that for the sweep problem there will not be a vector n that is normal to all the u k unless the geometry is pathologic. We wish to apply the least-squares fit above to find a relation of the form (13) ũ k = T u k. Summing this relation over k we find that T ( c c) = c c. Putting (11-12) into (13), we find that (14) x k c = T (x k c). We have thus arrived at relationship (6), which in turn implies (2), as desired. Note that if T is the identity matrix, then one gets loop translation. Note also that if F is not zero at the minimum, then T does not necessarily send {x k } to { x k }. A new linear transformation is calculated in Sweep Tool for each layer of the advancing front based on consecutive boundary loops from the linking surfaces. 4. Layer Smoothing The algorithm outlined in the previous section was implemented in Sweep Tool and was found effective on a wide variety of problems, including rotatable geometries (see Figure 1). Despite the effectiveness of linear transformations, another algorithm is needed to sweep general two-and-one-half dimensional geometries. To handle general geometries, we resort to layer smoothing techniques in which each layer of quadrilaterals on the advancing front is smoothed independently of connections to nodes on the previous layer. Although fully three-dimensional smoothing is appropriate for the hexahedral meshes generated in sweep, it would be very inefficient to incorporate such a procedure into the advancing front method. An alternative is layer smoothing which provides a relatively efficient, if somewhat less robust, approach. A variety of layer smoothing methods, including Laplacian and isoparametric, were tried in Project with mixed success. Although Project smoothing was generally an improvement over the Project point placement algorithm described in section 2, two difficulties adversely impact mesh quality. The first difficulty is illustrated by the example of sweeping a half-torus with the source surface meshed with a radially-biased circle primitive (Figure 2). Because the Project algorithms of section 2 produce meshes with negative Jacobian elements on rotatable geometries, layer smoothing is necessary. Laplace layer smoothing eliminates the negative Jacobians but the biasing on the source mesh is lost. More importantly, distorted hexahedral elements are created on the first layer of the advancing front
5 Y X Z Fig. 1. Sweep Tool on Rotatable Geometry via Linear Transformation Approach because there is not a smooth transition between the source and target meshes. In this example, the need for the original biasing of the source mesh to be preserved by the layer smoothing is obvious. The second difficulty occurs on volumes having non-convex or non-simply connected source surfaces. Layer smoothing is absolutely necessary on the example in Figure 3 because no linear transformation between the source and target loops exists. Project with Laplacian layer smoothing produces a bad target mesh in this example (note the compression of the quadrilaterals next to the inner boundary). The second difficulty arises in many other important sweepable geometries. Both difficulties are overcome by devising a robust smoother for Sweep Tool that takes mesh biasing into account and performs mesh copying/morphing 2. For example, the mesh in Figure 1 can be obtained with Sweep Tool, with or without layer smoothing. Sweep Tool smoothing also solves the diffulty illustrated in Figure 2. Fig- 2 Space limitations prevent discussion of the smoothing algorithm here. See [3] for a detailed presentation
6 X Z Y Fig. 2. Project Tool: Loss of Mesh Biasing due to Smoothing (Target on Right) ure 4 shows the result of sweeping with Sweep Tool layer smoothing on the example in Figure 3. Experience with a series of realistic problems has shown it is possible to sweep a wide variety of general two-and-one-half dimensional geometries with this smoother. If the curvature of a source surface is significantly different than the target surface, poor mesh quality may result even with layer smoothing. Curvature is somewhat preserved by layer smoothing if the bounding loop reflects the surface curvature but, in general, this is not the case. Research on methods to achieve good mesh quality for curved source and target surfaces is presently being pursued. The approach takes into account both source and target meshes. Sections 3 and 4 give two alternative approaches to placement of the interior nodes, one involving linear transformations and the other layer smoothing. Although linear transformations are less general than layer smoothing, they remain useful because, when they exist, it is computationally much faster to generate them than to perform layer smoothing. Fortunately, (7) gives a direct measure of whether or not such a transformation exists: if F = 0, then the linear transformation exists and smoothing need not be used because good quality is guaranteed, as shown in the following result. Proposition. Assume that a layer mesh has positive Jacobian at all mesh nodes. Let T be a linear transformation with positive determinant. Then T applied to this layer in the manner of (14) will result in another layer mesh which also has positive Jacobians at all of its nodes. Proof. For each quadrilateral of the mesh let U be the matrix [x 1 x 0,x 2 x 0,n] for some unit normal. We have T U = Ũ with the determinant of U positive by assumption. Then the determinant of Ũ is positive because it is the product of the
7 X Z Y Fig. 3. Project Tool: Poor Quality Target Mesh (right) Resulting from Laplacian Layer Smoothing determinant of two matrices which each have positive determinants. There are thus two criteria to determine whether or not smoothing is needed. If (i) F < ǫ where ǫ is some tolerance criterion and (ii) the linear transformation has positive determinant, then smoothing is not needed. Occasionally, layer smoothing will fail to produce adequate mesh quality even with the new smoother because, although the geometry and topology are two-andone-half dimensional, the mesh is fully three-dimensional. When layer smoothing fails, the situation can sometimes be salvaged by invoking the fully three-dimensional smoother as a post-processing procedure. 5. Additional Topics We briefly discuss in this section some Sweep topics not directly related to that of interior node placement. Poor quality meshes on the linking surfaces generally result in poor quality hexadrons generated by sweeping. It is therefore crucial to generate good meshes on the source and linking surfaces. Because surface meshing is outside the control of Sweep Tool, Sweep Tool checks the Jacobians of all the quadrilaterals on the linking surfaces before sweeping. If any Jacobians are negative, then sweep aborts because negative hexahedron Jacobians will surely result. Even when there are no bad quadrilaterals on the linking surfaces, one can still sweep out bad hexahedra due to problems with the linking surface meshes. Highly skewed meshes on the linking surface is a common cause of this problem. Sweep Tool contains a topology-based mesh-matching algorithm because in some
8 X ZY Fig. 4. Sweep Tool: Good Quality Target Mesh (right) Resulting from new Layer Smoothing Algorithm instances it is advantageous to have the target surface meshed prior to sweeping. If the target mesh topology is incompatible with the source mesh topology, the sweep is aborted, otherwise each quadrilateral on the last layer of the advancing front is matched with the appropriate quadrilateral of the target mesh. Multiple source surfaces are allowed in sweep, but multiple targets are not presently allowed. Mesh generation rates of up to 3000 hexahedrons per second have been generated with Sweep Tool on a 300 MHz workstation, compared to 1800 hexahdrons per second for Project Tool. It is possible to achieve this through careful memory-speed tradeoffs to eliminate quadratic-time calculations. 6. Summary and Conclusions. Sweep Tool provides a relatively fast and semi-automatic means of generating volume meshes on general two-and-one-half dimensional geometries. Interior point placement is crucial to the success of the algorithm. Two approaches for the placement of interior node points have been suggested that are robust and highly effective. The linear transformation approach is fast and gives high quality meshes on translatable, rotatable, and scalable geometries. The layer smoothing approach is slower but can give high quallity on more general geometries. Two criteria are given which form the basis for an automatic method for determining if smoothing is needed. Some geometries are border-line sweepable, i.e., even though the geometry is two-and-one-half dimensional, meshes on the linking surfaces conspire to make layer smoothing inadequate. Full 3D smoothing applied as a post-processing step may give
9 the desired mesh quality. Further development of Sweep Tool technology is underway to improve mesh quality in the case of highly curved source and target surfaces. Research on improved algorithms for sweeping to multiple targets is also in progress. Of course, not all geometries are two-and-one-half dimensional so research must continue on meshing the general problem. Acknowledgements The author would like to thank the members of the CUBIT team and the codes users at Sandia National Laboratories for many helpful suggestions. REFERENCES [1] T. Blacker, The Cooper Tool, Proceedings of the 5th International Meshing RoundTable, pp13-29, Pittsburg, Pennsylvania, October 10-11, [2] P. Knupp, On the invertibility of the isoparametric map, Comp. Meth. Appl. Mech. Engr., 78, p , [3] P. Knupp, Applications of Mesh Smoothing: Copy, Morph, and Sweep on Unstructured Quadrilateral Meshes, submitted for publication. [4] L. Mingwu, S. Benzley, G. Sjaardema, T. Tautges, A Multiple Source and Target Sweeping Method for Generating All Hexahedral Finite Element Meshes, Proceedings of the 5th International Meshing RoundTable, p , Pittsburg, Pennsylvania, October 10-11, [5] D. White, L. Mingwu, S. Benzley, G. Sjaardema, Automated Hexahedral Mesh Generation, Proceedings of the 4th International Meshing RoundTable, pp , October 16-17, Albuquerque, New Mexico, [6] David White, private communication.
Hexahedral Mesh Generation for Volumetric Image Data
Hexahedral Mesh Generation for Volumetric Image Data Jason Shepherd University of Utah March 27, 2006 Outline Hexahedral Constraints Topology Boundary Quality Zhang et al. papers Smoothing/Quality Existing
More informationThe Geode Algorithm: Combining Hex/Tet Plastering, Dicing and Transition Elements for Automatic, All-Hex Mesh Generation
The Geode Algorithm: Combining Hex/Tet Plastering, Dicing and Transition Elements for Automatic, All-Hex Mesh Generation Robert W. Leland 1 Darryl J. Melander 1 Ray W. Meyers 1 Scott A. Mitchell 1 Timothy
More informationAUTOMATIC SCHEME SELECTION FOR TOOLKIT HEX MESHING
AUTOMATIC SCHEME SELECTION FOR TOOLKIT HEX MESHING David R. White Sandia National Laboratories, Albuquerque, NM. drwhite@sandia.gov Timothy J. Tautges Sandia National Laboratories, Albuquerque, NM. tjtautg@sandia.gov
More informationA New Smoothing Algorithm for Quadrilateral and Hexahedral Meshes
A New Smoothing Algorithm for Quadrilateral and Hexahedral Meshes Sanjay Kumar Khattri Department of Mathematics, University of Bergen, Norway sanjay@mi.uib.no http://www.mi.uib.no/ sanjay Abstract. Mesh
More informationSURFACE MESH PROJECTION FOR HEXAHEDRAL MESH GENERATION BY SWEEPING
SURFACE MESH PROJECTION FOR HEXAHEDRAL MESH GENERATION BY SWEEPING Xevi Roca 1 Josep Sarrate 1 Antonio Huerta 1 1 Laboratori de Càlcul Numèric (LaCàN), Departament de Matemàtica Aplicada III, Universitat
More informationThe Whisker Weaving Algorithm: A Connectivity-Based Method for Constructing All Hexahedral Finite Element Meshes
The Whisker Weaving Algorithm: A Connectivity-Based Method for Constructing All Hexahedral Finite Element Meshes Timothy J. Tautges Comp Mechanics and Visualization Dept, Sandia National Laboratories,
More informationBasic LOgical Bulk Shapes (BLOBs) for Finite Element Hexahedral Mesh Generation
Basic LOgical Bulk Shapes (BLOBs) for Finite Element Hexahedral Mesh Generation Shang-Sheng Liu and Rajit Gadh Department of Mechanical Engineering University of Wisconsin - Madison Madison, Wisconsin
More information9. F. P. Preparata, M. I. Shamos, Computational Geometry an Introduction, Springer- Verlag, New York, 1985, pp
5. A. P. Gilkey, G. D. Sjaardema, GEN3D: A GENESIS Database 2D to 3D Transformation Program, SAND89-0485, Sandia National Laboratories, Albuquerque, New Mexico, March 1989. 6. M. S. Shephard, M. K. Georges,
More informationMesh Generation. Timothy J. Tautges. Principle Member Technical Staff Sandia National Laboratories
Mesh Generation Timothy J. Tautges Principle Member Technical Staff Sandia National Laboratories Adjunct Professor, Engineering Physics University of Wisconsin-Madison Sandia is a multiprogram laboratory
More informationAn automatic and general least-squares projection procedure for sweep meshing
An automatic and general least-squares projection procedure for sweep meshing Xevi Roca and Josep Sarrate Laboratori de Càlcul Numèric, ETSE de Camins, Canals i Ports de Barcelona, Universitat Politècnica
More informationTopological Issues in Hexahedral Meshing
Topological Issues in Hexahedral Meshing David Eppstein Univ. of California, Irvine Dept. of Information and Computer Science Outline I. What is meshing? Problem statement Types of mesh Quality issues
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 informationUnpublished manuscript.
Unpublished manuscript. Design and Use of Geometry Attributes To Support Mesh Generation and Analysis Applications Timothy J. Tautges Sandia National Laboratories 1, Albuquerque, NM, University of Wisconsin-Madison,
More informationCoupling of Smooth Faceted Surface Evaluations in the SIERRA FEA Code
Coupling of Smooth Faceted Surface Evaluations in the SIERRA FEA Code Timothy J. Tautges Steven J. Owen Sandia National Laboratories University of Wisconsin-Madison Mini-symposium on Computational Geometry
More informationSolids as point set. Solid models. Solid representation schemes (cont d) Solid representation schemes. Solid representation schemes (cont d)
Solid models Solid models developed to address limitations of wireframe modeling. Attempt was to create systems which create only complete representations. Modelers would support direct creation of 3D
More information1. Introduction to Constructive Solid Geometry (CSG)
opyright@010, YZU Optimal Design Laboratory. All rights reserved. Last updated: Yeh-Liang Hsu (010-1-10). Note: This is the course material for ME550 Geometric modeling and computer graphics, Yuan Ze University.
More informationChoosing Corners of Rectangles for Mapped Meshing
Choosing Corners of Rectangles for Mapped Meshing Scott A. Mitchell 1 Abstract 1 Consider mapping a regular i x j quadrilateral mesh of a rectangle onto a surface. The quality of the mapped mesh of the
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 informationWe consider the problem of rening quadrilateral and hexahedral element meshes. For
Rening quadrilateral and hexahedral element meshes R. Schneiders RWTH Aachen Lehrstuhl fur Angewandte Mathematik, insb. Informatik Ahornstr. 55, 5056 Aachen, F.R. Germany (robert@feanor.informatik.rwth-aachen.de)
More informationAutomatic hybrid mesh generation for the boundary face method
Boundary Elements and Other Mesh Reduction Methods XXXVI 139 Automatic hybrid mesh generation for the boundary face method Cheng Huang & Jianming Zhang State Key Laboratory of Advanced Design and Manufacturing
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 informationA Singular Example for the Averaged Mean Curvature Flow
To appear in Experimental Mathematics Preprint Vol. No. () pp. 3 7 February 9, A Singular Example for the Averaged Mean Curvature Flow Uwe F. Mayer Abstract An embedded curve is presented which under numerical
More informationHexahedral Mesh Generation using the Embedded Voronoi Graph
Hexahedral Mesh Generation using the Embedded Voronoi Graph Alla Sheffer, Michal Etzion, Ari Rappoport, Michel Bercovier Institute of Computer Science, The Hebrew University, Jerusalem 91904, Israel. sheffa
More informationHexahedral Meshing of Non-Linear Volumes Using Voronoi Faces and Edges
Hexahedral Meshing of Non-Linear Volumes Using Voronoi Faces and Edges Alla Sheffer and Michel Bercovier Institute of Computer Science, The Hebrew University, Jerusalem 91904, Israel. sheffa berco @cs.huji.ac.il.
More informationAutomatic All-Hex Topology Operations Using Edge Valence Prediction with Application to Localized Coarsening
Brigham Young University BYU ScholarsArchive All Theses and Dissertations 2011-03-17 Automatic All-Hex Topology Operations Using Edge Valence Prediction with Application to Localized Coarsening Timothy
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 informationThe Geometry of Carpentry and Joinery
The Geometry of Carpentry and Joinery Pat Morin and Jason Morrison School of Computer Science, Carleton University, 115 Colonel By Drive Ottawa, Ontario, CANADA K1S 5B6 Abstract In this paper we propose
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 informationProcessing 3D Surface Data
Processing 3D Surface Data Computer Animation and Visualisation Lecture 12 Institute for Perception, Action & Behaviour School of Informatics 3D Surfaces 1 3D surface data... where from? Iso-surfacing
More informationAn Efficient Paving Method of Pure Quad Mesh Generation
2016 International Conference on Intelligent Manufacturing and Materials (ICIMM 2016) ISBN: 978-1-60595-363-2 An Efficient Paving Method of Pure Quad Mesh Generation Yongcai Liu, Wenliang Chen and Yidong
More informationChapter 12 Solid Modeling. Disadvantages of wireframe representations
Chapter 12 Solid Modeling Wireframe, surface, solid modeling Solid modeling gives a complete and unambiguous definition of an object, describing not only the shape of the boundaries but also the object
More information5B.3 An Immersive Topology Environment for Meshing
5B.3 An Immersive Topology Environment for Meshing Steven J. Owen 1, Brett W. Clark 1, Darryl J. Melander 1, Michael Brewer 1, Jason F. Shepherd 1, Karl Merkley 2, Corey Ernst 2, and Randy Morris 2 1 Sandia
More informationShapefile Overlay Using a Doubly-Connected Edge List
Shapefile Overlay Using a Doubly-Connected Edge List Phil Katz and Stephen St.Vincent Swarthmore College 500 College Ave. Swarthmore, PA 19081 [pkatz1, sstvinc2]@swarthmore.edu Abstract We present a method
More informationMETHODS FOR MULTISWEEP AUTOMATION
. { METHODS FOR MULTISWEEP AUTOMATION Jason Shepherd Scott A. Mitchell Patrick Knupp and David White Sandia National Laboratories Albuquerque NM U.S.A. jfshephf2sandia.gov samitch@sandia. gov pknupp @sandia.gov
More information5. Create two planes and split the outer cylinder volume into four quadrants using the menus
Generating a 3D GAMBIT grid for concentric rotating cylinders Prepared by: Jagannath R Nanduri (jaggu.nanduri@mail.wvu.edu) 304-293-3111x2348 1. In order to overcome the problem of generating small size
More informationShape Modeling and Geometry Processing
252-0538-00L, Spring 2018 Shape Modeling and Geometry Processing Discrete Differential Geometry Differential Geometry Motivation Formalize geometric properties of shapes Roi Poranne # 2 Differential Geometry
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 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 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 informationDifferential Geometry: Circle Patterns (Part 1) [Discrete Conformal Mappinngs via Circle Patterns. Kharevych, Springborn and Schröder]
Differential Geometry: Circle Patterns (Part 1) [Discrete Conformal Mappinngs via Circle Patterns. Kharevych, Springborn and Schröder] Preliminaries Recall: Given a smooth function f:r R, the function
More informationTowards a unified measurement of quality for mixed elements
Towards a unified measurement of qualit for mixed elements technical report n 015/01, di utfsm Claudio Lobos Departmento de Informática Universidad Técnica Federico Santa María Santiago, Chile clobos@inf.utfsm.cl
More informationTopic: Orientation, Surfaces, and Euler characteristic
Topic: Orientation, Surfaces, and Euler characteristic The material in these notes is motivated by Chapter 2 of Cromwell. A source I used for smooth manifolds is do Carmo s Riemannian Geometry. Ideas of
More informationLinear Complexity Hexahedral Mesh Generation
Linear Complexity Hexahedral Mesh Generation David Eppstein Department of Information and Computer Science University of California, Irvine, CA 92717 http://www.ics.uci.edu/ eppstein/ Tech. Report 95-51
More informationGeometric Modeling Mortenson Chapter 11. Complex Model Construction
Geometric Modeling 91.580.201 Mortenson Chapter 11 Complex Model Construction Topics Topology of Models Connectivity and other intrinsic properties Graph-Based Models Emphasize topological structure Boolean
More information4 Integer Linear Programming (ILP)
TDA6/DIT37 DISCRETE OPTIMIZATION 17 PERIOD 3 WEEK III 4 Integer Linear Programg (ILP) 14 An integer linear program, ILP for short, has the same form as a linear program (LP). The only difference is that
More informationA NEW TYPE OF SIZE FUNCTION RESPECTING PREMESHED ENTITIES
A NEW TYPE OF SIZE FUNCTION RESPECTING PREMESHED ENTITIES Jin Zhu Fluent, Inc. 1007 Church Street, Evanston, IL, U.S.A. jz@fluent.com ABSTRACT This paper describes the creation of a new type of size function
More informationEngineering designs today are frequently
Basic CAD Engineering designs today are frequently constructed as mathematical solid models instead of solely as 2D drawings. A solid model is one that represents a shape as a 3D object having mass properties.
More informationComputing intersections in a set of line segments: the Bentley-Ottmann algorithm
Computing intersections in a set of line segments: the Bentley-Ottmann algorithm Michiel Smid October 14, 2003 1 Introduction In these notes, we introduce a powerful technique for solving geometric problems.
More informationGuaranteed-Quality All-Quadrilateral Mesh Generation with Feature Preservation
Guaranteed-Quality All-Quadrilateral Mesh Generation with Feature Preservation Xinghua Liang, Mohamed S. Ebeida, and Yongjie Zhang Department of Mechanical Engineering, Carnegie Mellon University, USA
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 informationUntil now we have worked with flat entities such as lines and flat polygons. Fit well with graphics hardware Mathematically simple
Curves and surfaces Escaping Flatland Until now we have worked with flat entities such as lines and flat polygons Fit well with graphics hardware Mathematically simple But the world is not composed of
More information2D rendering takes a photo of the 2D scene with a virtual camera that selects an axis aligned rectangle from the scene. The photograph is placed into
2D rendering takes a photo of the 2D scene with a virtual camera that selects an axis aligned rectangle from the scene. The photograph is placed into the viewport of the current application window. A pixel
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 informationThe goal is the definition of points with numbers and primitives with equations or functions. The definition of points with numbers requires a
The goal is the definition of points with numbers and primitives with equations or functions. The definition of points with numbers requires a coordinate system and then the measuring of the point with
More informationMinkowski Sums. Dinesh Manocha Gokul Varadhan. UNC Chapel Hill. NUS CS 5247 David Hsu
Minkowski Sums Dinesh Manocha Gokul Varadhan UNC Chapel Hill NUS CS 5247 David Hsu Last Lecture Configuration space workspace configuration space 2 Problem Configuration Space of a Translating Robot Input:
More informationSingularity Analysis of an Extensible Kinematic Architecture: Assur Class N, Order N 1
David H. Myszka e-mail: dmyszka@udayton.edu Andrew P. Murray e-mail: murray@notes.udayton.edu University of Dayton, Dayton, OH 45469 James P. Schmiedeler The Ohio State University, Columbus, OH 43210 e-mail:
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 informationCurves and Surfaces. Chapter 7. Curves. ACIS supports these general types of curves:
Chapter 7. Curves and Surfaces This chapter discusses the types of curves and surfaces supported in ACIS and the classes used to implement them. Curves ACIS supports these general types of curves: Analytic
More informationSurface Grid Generation
Applied Mathematical Sciences, Vol. 9, 205, no. 76, 376-3773 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.2988/ams.205.5430 Surface Grid Generation Bashar ogheib and Ali Elsaheli Department of Mathematics
More informationPreferred directions for resolving the non-uniqueness of Delaunay triangulations
Preferred directions for resolving the non-uniqueness of Delaunay triangulations Christopher Dyken and Michael S. Floater Abstract: This note proposes a simple rule to determine a unique triangulation
More informationPlanes Intersecting Cones: Static Hypertext Version
Page 1 of 12 Planes Intersecting Cones: Static Hypertext Version On this page, we develop some of the details of the plane-slicing-cone picture discussed in the introduction. The relationship between the
More informationParallel Computation of Spherical Parameterizations for Mesh Analysis. Th. Athanasiadis and I. Fudos University of Ioannina, Greece
Parallel Computation of Spherical Parameterizations for Mesh Analysis Th. Athanasiadis and I. Fudos, Greece Introduction Mesh parameterization is a powerful geometry processing tool Applications Remeshing
More informationChapter 4 Concepts from Geometry
Chapter 4 Concepts from Geometry An Introduction to Optimization Spring, 2014 Wei-Ta Chu 1 Line Segments The line segment between two points and in R n is the set of points on the straight line joining
More informationBoolean Component. Chapter 1. Boolean Operations
Chapter 1. Boolean Component Component: The Boolean Component (BOOL), in the bool directory, performs Boolean operations on the model topology of bodies, first finding the intersections between bodies,
More informationHexahedral Mesh Refinement Using an Error Sizing Function
Brigham Young University BYU ScholarsArchive All Theses and Dissertations 2011-06-01 Hexahedral Mesh Refinement Using an Error Sizing Function Gaurab Paudel Brigham Young University - Provo Follow this
More informationUsing Isosurface Methods for Visualizing the Envelope of a Swept Trivariate Solid
Using Isosurface Methods for Visualizing the Envelope of a Swept Trivariate Solid Jason Conkey Kenneth I. Joy Center for Image Processing and Integrated Computing Department of Computer Science University
More informationEXTREME POINTS AND AFFINE EQUIVALENCE
EXTREME POINTS AND AFFINE EQUIVALENCE The purpose of this note is to use the notions of extreme points and affine transformations which are studied in the file affine-convex.pdf to prove that certain standard
More informationcoding of various parts showing different features, the possibility of rotation or of hiding covering parts of the object's surface to gain an insight
Three-Dimensional Object Reconstruction from Layered Spatial Data Michael Dangl and Robert Sablatnig Vienna University of Technology, Institute of Computer Aided Automation, Pattern Recognition and Image
More informationAN INTERACTIVE MESH GENERATION ENVIRONMENT FOR GEOMETRY-BASED SIMULATIONS. PROGRESS UPDATE.
6th Workshop on Numerical Methods in Applied Science and Engineering (NMASE 07) Vall de Núria, 9 a 11 de enero de 2007 c LaCàN, www.lacan-upc.es AN INTERACTIVE MESH GENERATION ENVIRONMENT FOR GEOMETRY-BASED
More informationPractical Linear Algebra: A Geometry Toolbox
Practical Linear Algebra: A Geometry Toolbox Third edition Chapter 18: Putting Lines Together: Polylines and Polygons Gerald Farin & Dianne Hansford CRC Press, Taylor & Francis Group, An A K Peters Book
More informationGenerating analysis topology using virtual topology operators
Generating analysis topology using virtual topology operators Tierney, C. M., Sun, L., Robinson, T. T., & Armstrong, C. G. (2015). Generating analysis topology using virtual topology operators. Procedia
More informationOrientation of manifolds - definition*
Bulletin of the Manifold Atlas - definition (2013) Orientation of manifolds - definition* MATTHIAS KRECK 1. Zero dimensional manifolds For zero dimensional manifolds an orientation is a map from the manifold
More informationDigital Geometry Processing Parameterization I
Problem Definition Given a surface (mesh) S in R 3 and a domain find a bective F: S Typical Domains Cutting to a Disk disk = genus zero + boundary sphere = closed genus zero Creates artificial boundary
More information1 Automatic Mesh Generation
1 AUTOMATIC MESH GENERATION 1 1 Automatic Mesh Generation 1.1 Mesh Definition Mesh M is a discrete representation of geometric model in terms of its geometry G, topology T, and associated attributes A.
More informationL1 - Introduction. Contents. Introduction of CAD/CAM system Components of CAD/CAM systems Basic concepts of graphics programming
L1 - Introduction Contents Introduction of CAD/CAM system Components of CAD/CAM systems Basic concepts of graphics programming 1 Definitions Computer-Aided Design (CAD) The technology concerned with the
More informationGeostatistics 2D GMS 7.0 TUTORIALS. 1 Introduction. 1.1 Contents
GMS 7.0 TUTORIALS 1 Introduction Two-dimensional geostatistics (interpolation) can be performed in GMS using the 2D Scatter Point module. The module is used to interpolate from sets of 2D scatter points
More informationUntangling and Smoothing of Quadrilateral and Hexahedral Meshes
Paper 36 Untangling and Smoothing of Quadrilateral and Hexahedral Meshes Civil-Comp Press, 2012 Proceedings of the Eighth International Conference on Engineering Computational Technology, B.H.V. Topping,
More informationCOMPUTING CONSTRAINED DELAUNAY
COMPUTING CONSTRAINED DELAUNAY TRIANGULATIONS IN THE PLANE By Samuel Peterson, University of Minnesota Undergraduate The Goal The Problem The Algorithms The Implementation Applications Acknowledgments
More informationMöbius Transformations in Scientific Computing. David Eppstein
Möbius Transformations in Scientific Computing David Eppstein Univ. of California, Irvine School of Information and Computer Science (including joint work with Marshall Bern from WADS 01 and SODA 03) Outline
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 informationGeometric Templates for Improved Tracking Performance in Monte Carlo Codes
Joint International Conference on Supercomputing in Nuclear Applications and Monte Carlo 2013 (SNA + MC 2013) La Cité des Sciences et de l Industrie, Paris, France, October 27-31, 2013 Geometric Templates
More informationApplying Neural Network Architecture for Inverse Kinematics Problem in Robotics
J. Software Engineering & Applications, 2010, 3: 230-239 doi:10.4236/jsea.2010.33028 Published Online March 2010 (http://www.scirp.org/journal/jsea) Applying Neural Network Architecture for Inverse Kinematics
More informationDesign by Subdivision
Bridges 2010: Mathematics, Music, Art, Architecture, Culture Design by Subdivision Michael Hansmeyer Department for CAAD - Institute for Technology in Architecture Swiss Federal Institute of Technology
More information10.1 Overview. Section 10.1: Overview. Section 10.2: Procedure for Generating Prisms. Section 10.3: Prism Meshing Options
Chapter 10. Generating Prisms This chapter describes the automatic and manual procedure for creating prisms in TGrid. It also discusses the solution to some common problems that you may face while creating
More informationParameterization with Manifolds
Parameterization with Manifolds Manifold What they are Why they re difficult to use When a mesh isn t good enough Problem areas besides surface models A simple manifold Sphere, torus, plane, etc. Using
More informationChapter 18. Geometric Operations
Chapter 18 Geometric Operations To this point, the image processing operations have computed the gray value (digital count) of the output image pixel based on the gray values of one or more input pixels;
More informationUNTANGLING AND OPTIMIZATION OF UNSTRUCTURED HEXAHEDRAL MESHES
UNTANLIN AND OPTIMIZATION OF UNSTRUCTURED HEXAHEDRAL MESHES K. Kovalev (*), M. Delanaye (**), Ch. Hirsch (*) kvk@stro.vub.ac.be (*) Vrije Universiteit Brussel, Pleinlaan,, 1050 Brussels (**) NUMECA Int.,
More informationLesson 2: Wireframe Creation
Lesson 2: Wireframe Creation In this lesson you will learn how to create wireframes. Lesson Contents: Case Study: Wireframe Creation Design Intent Stages in the Process Reference Geometry Creation 3D Curve
More informationLinear Methods for Regression and Shrinkage Methods
Linear Methods for Regression and Shrinkage Methods Reference: The Elements of Statistical Learning, by T. Hastie, R. Tibshirani, J. Friedman, Springer 1 Linear Regression Models Least Squares Input vectors
More informationCOMP 558 lecture 19 Nov. 17, 2010
COMP 558 lecture 9 Nov. 7, 2 Camera calibration To estimate the geometry of 3D scenes, it helps to know the camera parameters, both external and internal. The problem of finding all these parameters is
More information6. Dicretization methods 6.1 The purpose of discretization
6. Dicretization methods 6.1 The purpose of discretization Often data are given in the form of continuous values. If their number is huge, model building for such data can be difficult. Moreover, many
More informationQuadrilateral Meshing by Circle Packing
Quadrilateral Meshing by Circle Packing Marshall Bern 1 David Eppstein 2 Abstract We use circle-packing methods to generate quadrilateral meshes for polygonal domains, with guaranteed bounds both on the
More information3. MODELING A THREE-PIPE INTERSECTION (3-D)
3. MODELING A THREE-PIPE INTERSECTION (3-D) This tutorial employs primitives that is, predefined GAMBIT modeling components and procedures. There are two types of GAMBIT primitives: Geometry Mesh Geometry
More informationA Patch-Based Mesh Optimization Algorithm for Partitioned Meshes
A Patch-Based Mesh Optimization Algorithm for Partitioned Meshes Nicholas Voshell 1, Suzanne Shontz 1, Lori Diachin 2, Patrick Knupp 3, and Todd Munson 4 1 Department of Computer Science and Engineering
More informationLecture 3: Linear Classification
Lecture 3: Linear Classification Roger Grosse 1 Introduction Last week, we saw an example of a learning task called regression. There, the goal was to predict a scalar-valued target from a set of features.
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 informationHexahedral Structured Grid Generation
Hexahedral Structured Grid Generation Javier de Jesús Cortés Aguirre. Pablo Barrera Sánchez. Guilmer González Flores. Universidad Nacional Autónoma de México (UNAM). VII Pan-American Workshop. Applied
More informationThree applications of Euler s formula. Chapter 10
Three applications of Euler s formula Chapter 10 A graph is planar if it can be drawn in the plane R without crossing edges (or, equivalently, on the -dimensional sphere S ). We talk of a plane graph if
More informationPerspective Mappings. Contents
Perspective Mappings David Eberly, Geometric Tools, Redmond WA 98052 https://www.geometrictools.com/ This work is licensed under the Creative Commons Attribution 4.0 International License. To view a copy
More informationInteger Programming Theory
Integer Programming Theory Laura Galli October 24, 2016 In the following we assume all functions are linear, hence we often drop the term linear. In discrete optimization, we seek to find a solution x
More information