Computing a Compact Spline Representation of the Medial Axis Transform of a 2D Shape

We present a full pipeline for computing the medial axis transform of an arbitrary 2D shape. The instability of the medial axis transform is overcome by a pruning algorithm guided by a user-defined Hausdorff distance threshold. The stable medial axis transform is then approximated by spline curves in 3D to produce a smooth and compact representation. These spline curves are computed by minimizing the approximation error between the input shape and the shape represented by the medial axis transform. Our results on various 2D shapes suggest that our method is practical and effective, and yields faithful and compact representations of medial axis transforms of 2D shapes.Comment: GMP14 (Geometric Modeling and Processing

A new method for aspherical surface fitting with large-volume datasets

In the framework of form characterization of aspherical surfaces, European National Metrology Institutes (NMIs) have been developing ultra-high precision machines having the ability to measure aspherical lenses with an uncertainty of few tens of nanometers. The fitting of the acquired aspherical datasets onto their corresponding theoretical model should be achieved at the same level of precision. In this article, three fitting algorithms are investigated: the Limited memory-Broyden-Fletcher-Goldfarb-Shanno (L-BFGS), the Levenberg–Marquardt (LM) and one variant of the Iterative Closest Point (ICP). They are assessed based on their capacities to converge relatively fast to achieve a nanometric level of accuracy, to manage a large volume of data and to be robust to the position of the data with respect to the model. Nev-ertheless, the algorithms are first evaluated on simulated datasets and their performances are studied. The comparison of these algorithms is extended on measured datasets of an aspherical lens. The results validate the newly used method for the fitting of aspherical surfaces and reveal that it is well adapted, faster and less complex than the LM or ICP methods.EMR

Automatic fitting of conical envelopes to free-form surfaces for flank CNC machining

We propose a new algorithm to detect patches of free-form surfaces that can be well approximated by envelopes of a rotational cone under a rigid body motion. These conical envelopes are a preferable choice from the manufacturing point of view as they are, by-definition, manufacturable by computer numerically controlled (CNC) machining using the efficient flank (peripheral) method with standard conical tools. Our geometric approach exploits multi-valued vector fields that consist of vectors in which the point-surface distance changes linearly. Integrating such vector fields gives rise to a family of integral curves, and, among them, linear segments that further serve as conical axes are quickly extracted. The lines that additionally admit tangential motion of the associated cone along the reference geometry form a set of candidate lines that are sequentially clustered and ordered to initialize motions of a rigid truncated cone. We validate our method by applying it on synthetic examples with exact envelopes, recovering correctly the exact solutions, and by testing it on several benchmark industrial datasets, detecting manufacturable conical envelope patches within fine tolerances

A sharp interface isogeometric strategy for moving boundary problems

The proposed methodology is first utilized to model stationary and propagating cracks. The crack face is enriched with the Heaviside function which captures the displacement discontinuity. Meanwhile, the crack tips are enriched with asymptotic displacement functions to reproduce the tip singularity. The enriching degrees of freedom associated with the crack tips are chosen as stress intensity factors (SIFs) such that these quantities can be directly extracted from the solution without a-posteriori integral calculation. As a second application, the Stefan problem is modeled with a hybrid function/derivative enriched interface. Since the interface geometry is explicitly defined, normals and curvatures can be analytically obtained at any point on the interface, allowing for complex boundary conditions dependent on curvature or normal to be naturally imposed. Thus, the enriched approximation naturally captures the interfacial discontinuity in temperature gradient and enables the imposition of Gibbs-Thomson condition during solidification simulation. The shape optimization through configuration of finite-sized heterogeneities is lastly studied. The optimization relies on the recently derived configurational derivative that describes the sensitivity of an arbitrary objective with respect to arbitrary design modifications of a heterogeneity inserted into a domain. The THB-splines, which serve as the underlying approximation, produce sufficiently smooth solution near the boundaries of the heterogeneity for accurate calculation of the configurational derivatives. (Abstract shortened by ProQuest.

On initialization of milling paths for 5-axis flank CNC machining of free-form surfaces with general milling tools

We propose a path-planning algorithm for 5-axis flank CNC machining with general tools of varying curvature. Our approach generalizes the initialization strategy introduced for conical tools [Bo et al., 2017] to arbitrary milling tools. Given a free-form (NURBS) surface and a rotational milling tool, we look for its motion in 3D to approximate the input reference surface within a given tolerance. We show that for a general shape of the milling tool, there exist locally and generically four 3D directions in which the point-surface distance follows the shape of the tool up to second order. These directions form a 3D multi-valued vector field and its integration gives rise to a set of integral curves. Among these integral curves, we seek straight line segments that correspond to good initial positions of the axes of the milling tool. We validate our method against synthetic examples with known exact solutions and, on industrial datasets, we detect approximate solutions that meet fine machining tolerances. We also demonstrate applicability of our method for efficient flank milling of convex regions that is not possible using traditional conical tools.RYC-2017-2264

May 12 1997 Cme Event: I. a Simplified Model of the Pre-Eruptive Magnetic Structure

A simple model of the coronal magnetic field prior to the CME eruption on May 12 1997 is developed. First, the magnetic field is constructed by superimposing a large-scale background field and a localized bipolar field to model the active region (AR) in the current-free approximation. Second, this potential configuration is quasi-statically sheared by photospheric vortex motions applied to two flux concentrations of the AR. Third, the resulting force-free field is then evolved by canceling the photospheric magnetic flux with the help of an appropriate tangential electric field applied to the central part of the AR. To understand the structure of the modeled configuration, we use the field line mapping technique by generalizing it to spherical geometry. It is demonstrated that the initial potential configuration contains a hyperbolic flux tube (HFT) which is a union of two intersecting quasi-separatrix layers. This HFT provides a partition of the closed magnetic flux between the AR and the global solar magnetic field. The vortex motions applied to the AR interlock the field lines in the coronal volume to form additionally two new HFTs pinched into thin current layers. Reconnection in these current layers helps to redistribute the magnetic flux and current within the AR in the flux-cancellation phase. In this phase, a magnetic flux rope is formed together with a bald patch separatrix surface wrapping around the rope. Other important implications of the identified structural features of the modeled configuration are also discussed.Comment: 25 pages, 11 figures, to appear in ApJ 200