Abstract

We present a method for the automatic adaption of the support size of meshfree basis functions in the context of the numerical approximation of boundary value problems stemming from a minimum principle. The method is based on a variational approach, and the central idea is that the variational principle selects both the discretized physical fields and the discretization parameters, here those defining the support size of each basis function. We consider local maximum-entropy approximation schemes, which exhibit smooth basis functions with respect to both space and the discretization parameters (the node location and the locality parameters). We illustrate by the Poisson, linear and non-linear elasticity problems the effectivity of the method, which produces very accurate solutions with very coarse discretizations and finds unexpected patterns of the support size of the shape functions.

Abstract

We present a method to process embedded smooth manifolds using sets of points alone. This method avoids any global parameterization and hence is applicable to surfaces of any genus. It combines three ingredients: (1) the automatic detection of the local geometric structure of the manifold by statistical learning methods; (2) the local parameterization of the surface using smooth meshfree (here maximum-entropy) approximants; and (3) patching together the local representations by means of a partition of unity. Mesh-based methods can deal with surfaces of complex topology, since they rely on the element-level parameterizations, but cannot handle high-dimensional manifolds, whereas previous meshfree methods for thin shells consider a global parametric domain, which seriously limits the kinds of surfaces that can be treated. We present the implementation of the method in the context of Kirchhoff–Love shells, but it is applicable to other calculations on manifolds in any dimension. With the smooth approximants, this fourth-order partial differential equation is treated directly. We show the good performance of the method on the basis of the classical obstacle course. Additional calculations exemplify the flexibility of the proposed approach in treating surfaces of complex topology and geometry.

Abstract

Calculations on general point-set surfaces are attractive because of their flexibility and simplicity in the preprocessing but present important challenges. The absence of a mesh makes it nontrivial to decide if two neighboring points in the three-dimensional embedding are nearby or rather far apart on the manifold. Furthermore, the topology of surfaces is generally not that of an open two-dimensional set, ruling out global parametrizations. We propose a general and simple numerical method analogous to the mathematical theory of manifolds, in which the point-set surface is described by a set of overlapping charts forming a complete atlas. We proceed in four steps: (1) partitioning of the node set into subregions of trivial topology; (2) automatic detection of the geometric structure of the surface patches by nonlinear dimensionality reduction methods; (3) parametrization of the surface using smooth meshfree (here maximum-entropy ) approximants; and (4) gluing together the patch representations by means of a partition of unity. Each patch may be viewed as a meshfree macro-element. We exemplify the generality, flexibility, and accuracy of the proposed approach by numerically approximating the geometrically nonlinear Kirchhoff–Love theory of thin-shells. We analyze standard benchmark tests as well as point-set surfaces of complex geometry and topology.

Abstract

We present a phase-field model for fracture in Kirchoff-Love thin shells using the local maximum-entropy (LME) meshfree method. Since the crack is a natural outcome of the analysis it does not require an explicit representation and tracking, which is advantage over techniques as the extended finite element method that requires tracking of the crack paths. The geometric description of the shell is based on statistical learning techniques that allow dealing with general point set surfaces avoiding a global parametrization, which can be applied to tackle surfaces of complex geometry and topology. We show the flexibility and robustness of the present methodology for two examples: plate in tension and a set of open connected pipes.

Abstract

Crack propagation in brittle materials with anisotropic surface energy is important in applications involving single crystals, extruded polymers, or geological and organic materials. Furthermore, when this anisotropy is strong, the phenomenology of crack propagation becomes very rich, with forbidden crack propagation directions or complex sawtooth crack patterns. This problem interrogates fundamental issues in fracture mechanics, including the principles behind the selection of crack direction. Here, we propose a variational phase-field model for strongly anisotropic fracture, which resorts to the extended Cahn-Hilliard framework proposed in the context of crystal growth. Previous phase-field models for anisotropic fracture were formulated in a framework only allowing for weak anisotropy. We implement numerically our higher-order phase-field model with smooth local maximum entropy approximants in a direct Galerkin method. The numerical results exhibit all the features of strongly anisotropic fracture and reproduce strikingly well recent experimental observations. rack propagation in brittle materials with anisotropic surface energy is important in applications involving single crystals, extruded polymers, or geological and organic materials. Furthermore, when this anisotropy is strong, the phenomenology of crack propagation becomes very rich, with forbidden crack propagation directions or complex sawtooth crack patterns. This problem interrogates fundamental issues in fracture mechanics, including the principles behind the selection of crack direction. Here, we propose a variational phase-field model for strongly anisotropic fracture, which resorts to the extended Cahn-Hilliard framework proposed in the context of crystal growth. Previous phase-field models for anisotropic fracture were formulated in a framework only allowing for weak anisotropy. We implement numerically our higher-order phase-field model with smooth local maximum entropy approximants in a direct Galerkin method. The numerical results exhibit all the features of strongly anisotropic fracture and reproduce strikingly well recent experimental observations.