50 research outputs found

    A Hermite interpolatory subdivision scheme for C2C^2-quintics on the Powell-Sabin 12-split

    Get PDF
    In order to construct a C1C^1-quadratic spline over an arbitrary triangulation, one can split each triangle into 12 subtriangles, resulting in a finer triangulation known as the Powell-Sabin 12-split. It has been shown previously that the corresponding spline surface can be plotted quickly by means of a Hermite subdivision scheme. In this paper we introduce a nodal macro-element on the 12-split for the space of quintic splines that are locally C3C^3 and globally C2C^2. For quickly evaluating any such spline, a Hermite subdivision scheme is derived, implemented, and tested in the computer algebra system Sage. Using the available first derivatives for Phong shading, visually appealing plots can be generated after just a couple of refinements.Comment: 17 pages, 7 figure

    SUBDIVIDE AND CONQUER RESOLUTION

    Get PDF
    This contribution will be freewheeling in the domain of signal, image and surface processing and touch briefly upon some topics that have been close to the heart of people in our research group. A lot of the research of the last 20 years in this domain that has been carried out world wide is dealing with multiresolution. Multiresolution allows to represent a function (in the broadest sense) at different levels of detail. This was not only applied in signals and images but also when solving all kinds of complex numerical problems. Since wavelets came into play in the 1980's, this idea was applied and generalized by many researchers. Therefore we use this as the central idea throughout this text. Wavelets, subdivision and hierarchical bases are the appropriate tools to obtain these multiresolution effects. We shall introduce some of the concepts in a rather informal way and show that the same concepts will work in one, two and three dimensions. The applications in the three cases are however quite different, and thus one wants to achieve very different goals when dealing with signals, images or surfaces. Because completeness in our treatment is impossible, we have chosen to describe two case studies after introducing some concepts in signal processing. These case studies are still the subject of current research. The first one attempts to solve a problem in image processing: how to approximate an edge in an image efficiently by subdivision. The method is based on normal offsets. The second case is the use of Powell-Sabin splines to give a smooth multiresolution representation of a surface. In this context we also illustrate the general method of construction of a spline wavelet basis using a lifting scheme

    Non-Uniform Rational B-Splines and Rational Bezier Triangles for Isogeometric Analysis of Structural Applications

    Full text link
    Isogeometric Analysis (IGA) is a major advancement in computational analysis that bridges the gap between a computer-aided design (CAD) model, which is typically constructed using Non-Uniform Rational B-splines (NURBS), and a computational model that traditionally uses Lagrange polynomials to represent the geometry and solution variables. In IGA, the same shape functions that are used in CAD are employed for analysis. The direct manipulation of CAD data eliminates approximation errors that emanate from the process of converting the geometry from CAD to Finite Element Analysis (FEA). As a result, IGA allows the exact geometry to be represented at the coarsest level and maintained throughout the analysis process. While IGA was initially introduced to streamline the design and analysis process, this dissertation shows that IGA can also provide improved computational results for complex and highly nonlinear problems in structural mechanics. This dissertation addresses various problems in structural mechanics in the context of IGA, with the use of NURBS and rational BĂ©zier triangles for the description of the parametric and physical spaces. The approaches considered here show that a number of important properties (e.g., high-order smoothness, geometric exactness, reduced number of degrees of freedom, and increased flexibility in discretization) can be achieved, leading to improved numerical solutions. Specifically, using B-splines and a layer-based discretization, a distributed plasticity isogeometric frame model is formulated to capture the spread of plasticity in large-deformation frames. The modeling approach includes an adaptive analysis where the structure of interest is initially modeled with coarse mesh and knots are inserted based on the yielding information at the quadrature points. It is demonstrated that improvement on efficiency and convergence rates is attained. With NURBS, an isogeometric rotation-free multi-layered plate formulation is developed based on a layerwise deformation theory. The derivation assumes a separate displacement field expansion within each layer, and considers transverse displacement component as C0-continuous at dissimilar material interfaces, which is enforced via knot repetition. The separate integration of the in-plane and through-thickness directions allows to capture the complete 3D stresses in a 2D setting. The proposed method is used to predict the behavior of advanced materials such as laminated composites, and the results show advantages in efficiency and accuracy. To increase the flexibility in discretizing complex geometries, rational BĂ©zier triangles for domain triangulation is studied. They are further coupled with a Delaunay-based feature-preserving discretization algorithm for static bending and free vibration analysis of Kirchhoff plates. Lagrange multipliers are employed to explicitly impose high-order continuity constraints and the augmented system is solved iteratively without increasing the matrix size. The resulting discretization is geometrically exact, admits small geometric features, and constitutes C1-continuity. The feature-preserving rational BĂ©zier triangles are further applied to smeared damage modeling of quasi-brittle materials. Due to the ability of Lagrange multipliers to raise global continuity to any desired order, the implicit fourth- and sixth-order gradient damage models are analyzed. The inclusion of higher-order terms in the nonlocal Taylor expansion improves solution accuracy. A local refinement algorithm that resolves marked regions with high resolution while keeping the resulting mesh conforming and well-conditioned is also utilized to improve efficiency. The outcome is a unified modeling framework where the feature-preserving discretization is able to capture the damage initiation and early-stage propagation, and the local refinement technique can then be applied to adaptively refine the mesh in the direction of damage propagation.PHDCivil EngineeringUniversity of Michigan, Horace H. Rackham School of Graduate Studieshttps://deepblue.lib.umich.edu/bitstream/2027.42/147668/1/ningliu_1.pd

    Splines for damage and fracture in solids

    Get PDF
    This thesis addresses different aspects of numerical fracture mechanics and spline technology for analysis. An energy-based arc-length control for physically non-linear problems is proposed. It switches between an internal energy-based and a dissipation-based arc-length method. The arc-length control allows to trace an equilibrium path with multiple snap-through and/or snap-back phenomena and only requires two parameters. Phase field models for brittle and cohesive fracture are numerically assessed. The impact of different parameters and boundary conditions on the phase field model for brittle fracture is investigated. It is demonstrated that Γ-convergence is not attained numerically for the phase field model for brittle fracture and that the phase field model for cohesive fracture does not pass a two-dimensional patch test when using an unstructured mesh. The properties of the Bézier extraction operator for T-splines are exploited for the determination of linear dependencies, partition of unity properties, nesting behaviour and local refinement. Unstructured T-spline meshes with extraordinary points are modified such that the blending functions fulfil the partition of unity property and possess a higher continuity. Bézier extraction for Powell-Sabin B-splines is introduced. Different spline technologies are compared when solving Kirchhoff-Love plate theory on a disc with simply supported and clamped boundary conditions. Powell-Sabin B-splines are utilised for smeared and discrete approaches to fracture. Due to the higher continuity of Powell-Sabin B-splines, the implicit fourth order gradient damage model for quasi-brittle materials can be solved and stresses can be computed directly at the crack tip when considering the cohesive zone method

    Computation of the crack opening displacement in the phase-field model

    Get PDF
    The phase-field model treats discrete cracks in a smeared sense by a regularisation technique. It holds attractive properties: there is no need to consider cracks as geometric discontinuities, and it avoids remeshing around crack tips. The method has been employed in the analysis of brittle and cohesive fracture problems. In the brittle fracture model, a Griffith-like energy functional is used in the simulation, while in the cohesive fracture model, the fracturing problem exploits a displacement jump governed energy functional. Obviously, the displacement jump is crucial in the cohesive fracture model and in certain other applications, e.g. for hydraulic fracturing. In the current study, the approximated form of the crack opening displacement is derived for the brittle and cohesive fracture models. In both models, the crack opening displacement is associated with a line integral normal to the crack, but different factors in front of the integral apply. The derived integrals are verified analytically in a one-dimensional setting and numerically in multi-dimensional set-ups, featuring straight and curved cracks

    New Techniques for the Modeling, Processing and Visualization of Surfaces and Volumes

    Get PDF
    With the advent of powerful 3D acquisition technology, there is a growing demand for the modeling, processing, and visualization of surfaces and volumes. The proposed methods must be efficient and robust, and they must be able to extract the essential structure of the data and to easily and quickly convey the most significant information to a human observer. Independent of the specific nature of the data, the following fundamental problems can be identified: shape reconstruction from discrete samples, data analysis, and data compression. This thesis presents several novel solutions to these problems for surfaces (Part I) and volumes (Part II). For surfaces, we adopt the well-known triangle mesh representation and develop new algorithms for discrete curvature estimation,detection of feature lines, and line-art rendering (Chapter 3), for connectivity encoding (Chapter 4), and for topology preserving compression of 2D vector fields (Chapter 5). For volumes, that are often given as discrete samples, we base our approach for reconstruction and visualization on the use of new trivariate spline spaces on a certain tetrahedral partition. We study the properties of the new spline spaces (Chapter 7) and present efficient algorithms for reconstruction and visualization by iso-surface rendering for both, regularly (Chapter 8) and irregularly (Chapter 9) distributed data samples
    corecore