Subdivision Shell Elements with Anisotropic Growth

A thin shell finite element approach based on Loop's subdivision surfaces is proposed, capable of dealing with large deformations and anisotropic growth. To this end, the Kirchhoff-Love theory of thin shells is derived and extended to allow for arbitrary in-plane growth. The simplicity and computational efficiency of the subdivision thin shell elements is outstanding, which is demonstrated on a few standard loading benchmarks. With this powerful tool at hand, we demonstrate the broad range of possible applications by numerical solution of several growth scenarios, ranging from the uniform growth of a sphere, to boundary instabilities induced by large anisotropic growth. Finally, it is shown that the problem of a slowly and uniformly growing sheet confined in a fixed hollow sphere is equivalent to the inverse process where a sheet of fixed size is slowly crumpled in a shrinking hollow sphere in the frictionless, quasi-static, elastic limit.Comment: 20 pages, 12 figures, 1 tabl

Orthotropic rotation-free thin shell elements

A method to simulate orthotropic behaviour in thin shell finite elements is proposed. The approach is based on the transformation of shape function derivatives, resulting in a new orthogonal basis aligned to a specified preferred direction for all elements. This transformation is carried out solely in the undeformed state leaving minimal additional impact on the computational effort expended to simulate orthotropic materials compared to isotropic, resulting in a straightforward and highly efficient implementation. This method is implemented for rotation-free triangular shells using the finite element framework built on the Kirchhoff--Love theory employing subdivision surfaces. The accuracy of this approach is demonstrated using the deformation of a pinched hemispherical shell (with a 18{\deg} hole) standard benchmark. To showcase the efficiency of this implementation, the wrinkling of orthotropic sheets under shear displacement is analyzed. It is found that orthotropic subdivision shells are able to capture the wrinkling behavior of sheets accurately for coarse meshes without the use of an additional wrinkling model.Comment: 10 pages, 8 figure

A variational model of fracture for tearing brittle thin sheets

Tearing of brittle thin elastic sheets, possibly adhered to a substrate, involves a rich interplay between nonlinear elasticity, geometry, adhesion, and fracture mechanics. In addition to its intrinsic and practical interest, tearing of thin sheets has helped elucidate fundamental aspects of fracture mechanics including the mechanism of crack path selection. A wealth of experimental observations in different experimental setups is available, which has been often rationalized with insightful yet simplified theoretical models based on energetic considerations. In contrast, no computational method has addressed tearing in brittle thin elastic sheets. Here, motivated by the variational nature of simplified models that successfully explain crack paths in tearing sheets, we present a variational phase-field model of fracture coupled to a nonlinear Koiter thin shell model including stretching and bending. We show that this general yet straightforward approach is able to reproduce the observed phenomenology, including spiral or power-law crack paths in free standing films, or converging/diverging cracks in thin films adhered to negatively/positively curved surfaces, a scenario not amenable to simple models. Turning to more quantitative experiments on thin sheets adhered to planar surfaces, our simulations allow us to examine the boundaries of existing theories and suggest that homogeneous damage induced by moving folds is responsible for a systematic discrepancy between theory and experiments. Thus, our computational approach to tearing provides a new tool to understand these complex processes involving fracture, geometric nonlinearity and delamination, complementing experiments and simplified theories.Fil: Li, Bin. Universidad Politécnica de Catalunya; España. Sorbonne Université; Francia. Centre National de la Recherche Scientifique; FranciaFil: Millán, Raúl Daniel. Universidad Nacional de Cuyo. Facultad de Ciencias Aplicadas a la Industria; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Mendoza; Argentina. Universidad Politécnica de Catalunya; EspañaFil: Torres Sánchez, Alejandro. Universidad Politécnica de Catalunya; EspañaFil: Roman, Benoît. Centre National de la Recherche Scientifique; Francia. Sorbonne Université; FranciaFil: Arroyo Balaguer, Marino. Universidad Politécnica de Catalunya; Españ

The variational approach to brittle fracture in materials with anisotropic surface energy and in thin sheets

Fracture mechanics of brittle materials has focused on bulk materials with isotropic surface energy. In this situation different physical principles for crack path selection are very similar or even equivalent. The situation is radically different when considering crack propagation in brittle materials with anisotropic surface energy. Such materials are important in applications involving single crystals, extruded polymers, or geological and organic materials. When this anisotropy is strong, the phenomenology of crack propagation becomes very rich, with forbidden crack propagation directions or complex sawtooth crack patterns. Thus, this situation interrogates fundamental issues in fracture mechanics, including the principles behind the selection of crack direction. Similarly, tearing of brittle thin elastic sheets, ubiquitous in nature, technology and daily life, challenges our understanding of fracture. Since tearing typically involves large geometric nonlinearity, it is not clear whether the stress intensity factors are meaningful or if and how they determine crack propagation. Geometry, together with the interplay between stretching and bending deformation, leads to complex behaviors, restricting analytical approximate solutions to very simplified settings and specific parameter regimes. In both situations, a rich and nontrivial experimental record has been successfully understood in terms of simple energetic models. However, general modeling approaches to either fracture in the presence of strong surface energy anisotropy or to tearing, capable of exploring new physics, have been lacking. The success of energetic simple models suggests that variational theories of brittle fracture may provide a unifying and general framework capable of dealing with the more general situations considered here. To address fracture in materials with strongly anisotropic surface energy, we propose a variational phase-field model resorting 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. To explore tearing of thin films, we develop a geometrically exact model and a computational framework coupling elasticity (stretching and bending), fracture, and adhesion to a substrate. We numerically implement the model with subdivision surface finite elements. Our simulations qualitatively and quantitatively reproduced the crack patterns observed in tearing experiments. Finally, we examine how shell geometry affects fracture. As suggested by previous results and our own phase-field simulations, shell shape dramatically affects crack evolution and the effective toughness of the shell structure. To gain insight and eventually develop new concepts for optimizing the design of thin shell structures, we derive the configurational force conjugate to crack extension for Koiter's linear thin shell theory. We identify the conservative contribution to this force through an Eshelby tensor, as well as non-conservative contributions arising from curvature.La mécanica de fractura frágil se ha centrado en materiales tridimensionales con una energía de superficie isotrópica. En esta situación, los diferentes principios para la selección del camino de la fisura son muy similares, o incluso equivalentes. La situación es radicalmente opuesta cuando se considera la propagación de fisuras en medios con energía de superficie anisótropa. Estos materiales son importantes en aplicaciones que involucran materiales cristalinos, polímeros extrudidos, o materiales orgánicos y geológicos. Cuando la anisotropía es fuerte, el fenómeno de la propagación de fisuras es muy rico, con direcciones de propagación prohibidas o complejos patrones de ruptura en dientes de sierra. Por tanto, esta situación plantea cuestiones fundamentales en la mecánica de la fractura, incluyendo los principios de selección de la dirección de propagación de la fractura. Igualmente, el proceso de rasgado de láminas delgadas y frágiles, comunes en la naturaleza, la tecnología y la vida diaria, desafía nuestro entendimiento de la fractura. Dado que el rasgado de estas láminas típicamente involucra grandes no linealidades geométricas, no está claro si los factores de intensidad de esfuerzos son válidos o si, y en tal caso cómo determinan la propagación de fisuras. La interacción entre la geometría, las deformaciones y la curvatura da lugar a comportamientos complejos, lo que restringe las soluciones analíticas aproximadas a ejemplos muy simplificados y a regímenes de parámetros limitados. En ambas situaciones, se han podido interpretar experimentos no triviales con modelos energéticos simples. Sin embargo, no se ha profundizado en modelos generales de fractura en presencia de energía de superficie fuertemente anisótropa o en láminas delgadas, ambas interesantes por su capacidad para explorar nueva física. El mencionado éxito de los modelos energéticos simplificados sugiere que las teorías variacionales de fractura en medios frágiles pueden proveer un marco unificador para considerar situaciones más generales, como las que se consideran en este trabajo. Para caracterizar la fractura en materiales con energía de superficie fuertemente anisótropa, proponemos un modelo variacional de campo de fase basado en el modelo extendido de Cahn-Hilliard. Los modelos de campo de fase existentes para la fractura anisótropa fueron formulados en un contexto que sólo admite anisotropía débil. En este trabajo, implementamos numéricamente nuestro modelo de campo de fase de alto orden con aproximantes locales de máxima entropía en un método directo de Garlerkin. Los resultados numéricos muestran todas las características de fractura con anisotropía fuerte, y reproducen llamativamente bien las últimas observaciones experimentales. Para explorar el rasgado de láminas delgadas, desarrollamos un modelo geométricamente exacto y un esquema computacional que acopla elasticidad (estiramiento y flexión), fractura, y la adhesión a un substrato. Implementamos numéricamente el modelo con elementos finitos basados en superficies de subdivisión. Nuestras simulaciones reproducen los patrones de ruptura, tanto cualitativamente como cuantitativamente, observados en los experimentos de rasgado. Finalmente, examinamos cómo la geometría de la lámina afecta la fractura. Como ha sido sugerido en resultados previos y en nuestras propias simulaciones de campo de fase, la forma de la lámina afecta dramáticamente la evolución de fisuras y la resistencia efectiva del material. Para comprender mejor estos fenómenos y con el objetivo de desarrollar nuevos conceptos para la optimización del diseño de estructuras de láminas delgadas, derivamos la fuerza configuracional conjugada a la extensión de la fractura para la teoría lineal de láminas delgadas de Koiter. Identificamos las contribuciones conservativas a esta fuerza a través del tensor de Eshelby, así como las contriuciones no conservativas que aparecen por el efecto de la curvatura

Mini-Workshop: Analytical and Numerical Methods in Image and Surface Processing

The workshop successfully brought together researchers from mathematical analysis, numerical mathematics, computer graphics and image processing. The focus was on variational methods in image and surface processing such as active contour models, Mumford-Shah type functionals, image and surface denoising based on geometric evolution problems in image and surface fairing, physical modeling of surfaces, the restoration of images and surfaces using higher order variational formulations

Fracture of solar-grade anisotropic polycrystalline Silicon: A combined phase field–cohesive zone model approach

Artículo Open Access en el sitio web del editor. Pago por publicar en abierto.This work presents a novel computational framework to simulate fracture events in brittle anisotropic polycrystalline materials at the microscopical level, with application to solar-grade polycrystalline Silicon. Quasi-static failure is modeled by combining the phase field approach of brittle fracture (for transgranular fracture) with the cohesive zone model for the grain boundaries (for intergranular fracture) through the generalization of the recent FE-based technique published in [M. Paggi, J. Reinoso, Comput. Methods Appl. Mech. Engrg., 31 (2017) 145–172] to deal with anisotropic polycrystalline microstructures. The proposed model, which accounts for any anisotropic constitutive tensor for the grains depending on their preferential orientation, as well as an orientation-dependent fracture toughness, allows to simulate intergranular and transgranular crack growths in an efficient manner, with or without initial defects. One of the advantages of the current variational method is the fact that complex crack patterns in such materials are triggered without any user-intervention, being possible to account for the competition between both dissipative phenomena. In addition, further aspects with regard to the model parameters identification are discussed in reference to solar cells images obtained from transmitted light source. A series of representative numerical simulations is carried out to highlight the interplay between the different types of fracture occurring in solar-grade polycrystalline Silicon, and to assess the role of anisotropy on the crack path and on the apparent tensile strength of the material.Unión Europea FP/2007–2013/ERC 306622Ministerio de Economía y Competitividad MAT2015–71036-P y MAT2015–71309-PJunta de Andalucía P11-TEP-7093 y P12-TEP- 105

Approximation of tensor fields on surfaces of arbitrary topology based on local Monge parametrizations

We introduce a new method, the Local Monge Parametrizations (LMP) method, to approximate tensor fields on general surfaces given by a collection of local parametrizations, e.g.~as in finite element or NURBS surface representations. Our goal is to use this method to solve numerically tensor-valued partial differential equations (PDE) on surfaces. Previous methods use scalar potentials to numerically describe vector fields on surfaces, at the expense of requiring higher-order derivatives of the approximated fields and limited to simply connected surfaces, or represent tangential tensor fields as tensor fields in 3D subjected to constraints, thus increasing the essential number of degrees of freedom. In contrast, the LMP method uses an optimal number of degrees of freedom to represent a tensor, is general with regards to the topology of the surface, and does not increase the order of the PDEs governing the tensor fields. The main idea is to construct maps between the element parametrizations and a local Monge parametrization around each node. We test the LMP method by approximating in a least-squares sense different vector and tensor fields on simply connected and genus-1 surfaces. Furthermore, we apply the LMP method to two physical models on surfaces, involving a tension-driven flow (vector-valued PDE) and nematic ordering (tensor-valued PDE). The LMP method thus solves the long-standing problem of the interpolation of tensors on general surfaces with an optimal number of degrees of freedom.Comment: 16 pages, 6 figure

Improved variational description of the Wick-Cutkosky model with the most general quadratic trial action

We generalize the worldline variational approach to field theory by introducing a trial action which allows for anisotropic terms to be induced by external 4-momenta of Green's functions. By solving the ensuing variational equations numerically we demonstrate that within the (quenched) scalar Wick-Cutkosky model considerable improvement can be achieved over results obtained previously with isotropic actions. In particular, the critical coupling associated with the instability of the model is lowered, in accordance with expectations from Baym's proof of the instability in the unquenched theory. The physical picture associated with a different quantum mechanical motion of the dressed particle along and perpendicular to its classical momentum is discussed. Indeed, we find that for large couplings the dressed particle is strongly distorted in the direction of its four-momentum. In addition, we obtain an exact relation between the renormalized coupling of the theory and the propagator. Along the way we introduce new and efficient methods to evaluate the averages needed in the variational approach and apply them to the calculation of the 2-point function.Comment: 32 pages, 4 figures, Latex. Some typos corrected and expanded discussion of the instability of the model provided. Accepted in Eur. Phys. J.