Scan registration for autonomous mining vehicles using 3D-NDT

Scan registration is an essential subtask when building maps based on range finder data from mobile robots. The problem is to deduce how the robot has moved between consecutive scans, based on the shape of overlapping portions of the scans. This paper presents a new algorithm for registration of 3D data. The algorithm is a generalization and improvement of the normal distributions transform (NDT) for 2D data developed by Biber and Strasser, which allows for accurate registration using a memory-efficient representation of the scan surface. A detailed quantitative and qualitative comparison of the new algorithm with the 3D version of the popular ICP (iterative closest point) algorithm is presented. Results with actual mine data, some of which were collected with a new prototype 3D laser scanner, show that the presented algorithm is faster and slightly more reliable than the standard ICP algorithm for 3D registration, while using a more memory efficient scan surface representation

Learning single-image 3D reconstruction by generative modelling of shape, pose and shading

We present a unified framework tackling two problems: class-specific 3D reconstruction from a single image, and generation of new 3D shape samples. These tasks have received considerable attention recently; however, most existing approaches rely on 3D supervision, annotation of 2D images with keypoints or poses, and/or training with multiple views of each object instance. Our framework is very general: it can be trained in similar settings to existing approaches, while also supporting weaker supervision. Importantly, it can be trained purely from 2D images, without pose annotations, and with only a single view per instance. We employ meshes as an output representation, instead of voxels used in most prior work. This allows us to reason over lighting parameters and exploit shading information during training, which previous 2D-supervised methods cannot. Thus, our method can learn to generate and reconstruct concave object classes. We evaluate our approach in various settings, showing that: (i) it learns to disentangle shape from pose and lighting; (ii) using shading in the loss improves performance compared to just silhouettes; (iii) when using a standard single white light, our model outperforms state-of-the-art 2D-supervised methods, both with and without pose supervision, thanks to exploiting shading cues; (iv) performance improves further when using multiple coloured lights, even approaching that of state-of-the-art 3D-supervised methods; (v) shapes produced by our model capture smooth surfaces and fine details better than voxel-based approaches; and (vi) our approach supports concave classes such as bathtubs and sofas, which methods based on silhouettes cannot learn.Comment: Extension of arXiv:1807.09259, accepted to IJCV. Differentiable renderer available at https://github.com/pmh47/dir

Feature Adaptive Ray Tracing of Subdivision Surfaces

abstract: Subdivision surfaces have gained more and more traction since it became the standard surface representation in the movie industry for many years. And Catmull-Clark subdivision scheme is the most popular one for handling polygonal meshes. After its introduction, Catmull-Clark surfaces have been extended to several eminent ways, including the handling of boundaries, infinitely sharp creases, semi-sharp creases, and hierarchically defined detail. For ray tracing of subdivision surfaces, a common way is to construct spatial bounding volume hierarchies on top of input control mesh. However, a high-level refined subdivision surface not only requires a substantial amount of memory storage, but also causes slow and inefficient ray tracing. In this thesis, it presents a new way to improve the efficiency of ray tracing of subdivision surfaces, while the quality is not as good as general methods.Dissertation/ThesisMasters Thesis Computer Science 201

3-D inelastic analysis methods for hot section components (base program)

A 3-D inelastic analysis methods program consists of a series of computer codes embodying a progression of mathematical models (mechanics of materials, special finite element, boundary element) for streamlined analysis of combustor liners, turbine blades, and turbine vanes. These models address the effects of high temperatures and thermal/mechanical loadings on the local (stress/strain) and global (dynamics, buckling) structural behavior of the three selected components. These models are used to solve 3-D inelastic problems using linear approximations in the sense that stresses/strains and temperatures in generic modeling regions are linear functions of the spatial coordinates, and solution increments for load, temperature and/or time are extrapolated linearly from previous information. Three linear formulation computer codes, referred to as MOMM (Mechanics of Materials Model), MHOST (MARC-Hot Section Technology), and BEST (Boundary Element Stress Technology), were developed and are described

Putting holes in holey geometry: Topology change for arbitrary surfaces

This paper presents a method for computing topology changes for triangle meshes in an interactive geometric modeling environment. Most triangle meshes in practice do not exhibit desirable geometric properties, so we develop a solution that is independent of standard assumptions and robust to geometric errors. Specifically, we provide the first method for topology change applicable to arbitrary non-solid, non-manifold, non-closed, self-intersecting surfaces. We prove that this new method for topology change produces the expected conventional results when applied to solid (closed, manifold, non-self-intersecting) surfaces---that is, we prove a backwards-compatibility property relative to prior work. Beyond solid surfaces, we present empirical evidence that our method remains tolerant to a variety of surface aberrations through the incorporation of a novel error correction scheme. Finally, we demonstrate how topology change applied to non-solid objects enables wholly new and useful behaviors

Computational Gradient Elasticity and Gradient Plasticity with Adaptive Splines

Classical continuum mechanics theories are largely insufficient in capturing size effects observed in many engineering materials: metals, composites, rocks etc. This is attributed to the absence of a length scale that accounts for microstructural effects inherent in these materials. Enriching the classical theories with an internal length scale solves this problem. One way of doing this, in a theoretically sound manner, is introducing higher order gradient terms in the constitutive relations. In elasticity, introducing a length scale removes the singularity observed at crack tips using the classical theory. In plasticity, it eliminates the spurious mesh sensitivity observed in softening and localisation problems by defining the width of the localisation zone thereby maintaining a well-posed boundary value problem. However, this comes at the cost of more demanding solution techniques. Higher-order continuity is usually required for solving gradient-enhanced continuum theories, a requirement difficult to meet using traditional finite elements. Hermitian finite elements, mixed methods and meshless methods have been developed to meet this requirement, however these methods have their drawbacks in terms of efficiency, robustness or implementational convenience. Isogeometric analysis, which exploits spline-based shape functions, naturally incorporates higher-order continuity, in addition to capturing the exact geometry and expediting the design-through-analysis process. Despite its potentials, it is yet to be fully explored for gradient-enhanced continua. Hence, this thesis develops an isogeometric analysis framework for gradient elasticity and gradient plasticity. The linearity of the gradient elasticity formulation has enabled an operator-split approach so that instead of solving the fourth-order partial differential equation monolithically, a set of two partial differential equations is solved in a staggered manner. A detailed convergence analysis is carried out for the original system and the split set using NURBS and T-splines. Suboptimal convergence rates in the monolithic approach and the limitations of the staggered approach are substantiated. Another advantage of the spline-based approach adopted in this work is the ease with which different orders of interpolation can be achieved. This is useful for consistency, and relevant in gradient plasticity where the local (explicit formulation) or nonlocal (implicit formulation) effective plastic strain needs to be discretised in addition to the displacements. Using different orders of interpolation, both formulations are explored in the second-order and a fourth-order implicit gradient formulation is proposed. Results, corroborated by dispersion analysis, show that all considered models give good regularisation with mesh-independent results. Compared with finite element approaches that use Hermitian shape functions for the plastic multiplier or mixed finite element approaches, isogeometric analysis has the distinct advantage that no interpolation of derivatives is required. In localisation problems, numerical accuracy requires the finite element size employed in simulations to be smaller than the internal length scale. Fine meshes are also needed close to regions of geometrical singularities or high gradients. Maintaining a fine mesh globally can incur high computational cost especially when considering large structures. To achieve this efficiently, selective refinement of the mesh is therefore required. In this context, splines need to be adapted to make them analysis-suitable. Thus, an adaptive isogeometric analysis framework is also developed for gradient elasticity and gradient plasticity. The proposed scheme does not require the mesh size to be smaller than the length scale, even during analysis, until a localisation band develops upon which adaptive refinement is performed. Refinement is based on a multi-level mesh with truncated hierarchical basis functions interacting through an inter-level subdivision operator. Through Bezier extraction, truncation of the bases is simplified by way of matrix multiplication, and an element-wise standard finite element data structure is maintained. In sum, a robust computational framework for engineering analysis is established, combining the flexibility, exact geometry representation, and expedited design-through analysis of isogeometric analysis, size-effect capabilities and mesh-objective results of gradient-enhanced continua, the standard convenient data structure of finite element analysis and the improved efficiency of adaptive hierarchical refinement

Relative Geologic Time By Dynamic Time Warping

This thesis considers an approach to tackle a core problem within seismic interpretation, which is bringing an autonomously generated interpretation of the seismic data, which is now known as a Relative Geologic Time. The proposed method readily utilizes the method of Dynamic Time Warping, which is an established method within signal processing. Using Dynamic Time Warping is thought to replicate similar interpretations an interpreter would conduct when fulfilling an interpretation of the subsurface. Utilizing Dynamic Time Warping to seismic data results in a fully autonomous interpretation of the subsurface, conducted in minutes and seconds. The method is simple and extendable, which can easily be further expanded. The workflow established during the thesis work results in a method that successfully produces an RGT volume. However, problems related to the method must be improved to enhance the outcome further and diminish errors present in the result. Furthermore, even with problems associated with the method, potential solutions are described in detail in the discussion and appendix. Discussion affiliated with previous attempts in solving Relative Geologic Time volumes is emphasized. The research conducted in Dynamic Time Warping is promising and emits potential for further research. LaTeX setup by Gunn and Patel (2017)

Lagrangian-on-Lagrangian Garment Design

Since the discovery of elastomeric materials, such as spandex or lycra, skintight clothing has revolutionized many different areas of the clothing industry, such as body-shaping clothing, athletic wear, and medical garments, among others. Often, this kind of clothing is designed to fulfill a given purpose, such as providing comfort, mobility, or improving recovery in the case of an athlete, provide support or exert some desired pressure in the case of medical garments, or actively deform the body to acquire some desired shape. Additionally, some designs aim to improve the life of the garment by, for example, minimizing tractions across the seams. While many tight-skin garments are sold in the market for generic body shapes, many of the purposes here mentioned are only achievable through a personalized fitting. To this end, we introduce a novel model, where the cloth is modeled as a membrane, parameterized as a function of the body. The cloth, is then able to slide on the body and deform it while staying always in contact. We call this model Lagrangian-on-Lagrangian. Based on this model, we develop an optimization framework, based on sensitivity analysis, capable of developing sewable patterns such that, when worn by a person, satisfy a given design target. With the framework, we include several design targets such as, body shape, stretch, pressure, sliding under motion, and seam traction. We evaluate our method on a variety of applications, as well as body shapes

Development of an integrated BEM approach for hot fluid structure interaction: BEST-FSI: Boundary Element Solution Technique for Fluid Structure Interaction

As part of the continuing effort at NASA LeRC to improve both the durability and reliability of hot section Earth-to-orbit engine components, significant enhancements must be made in existing finite element and finite difference methods, and advanced techniques, such as the boundary element method (BEM), must be explored. The BEM was chosen as the basic analysis tool because the critical variables (temperature, flux, displacement, and traction) can be very precisely determined with a boundary-based discretization scheme. Additionally, model preparation is considerably simplified compared to the more familiar domain-based methods. Furthermore, the hyperbolic character of high speed flow is captured through the use of an analytical fundamental solution, eliminating the dependence of the solution on the discretization pattern. The price that must be paid in order to realize these advantages is that any BEM formulation requires a considerable amount of analytical work, which is typically absent in the other numerical methods. All of the research accomplishments of a multi-year program aimed toward the development of a boundary element formulation for the study of hot fluid-structure interaction in Earth-to-orbit engine hot section components are detailed. Most of the effort was directed toward the examination of fluid flow, since BEM's for fluids are at a much less developed state. However, significant strides were made, not only in the analysis of thermoviscous fluids, but also in the solution of the fluid-structure interaction problem