In this paper we consider the p-FEM for elliptic boundary value problems on tetrahedral meshes where the entries of the stiffness matrix are evaluated by numerical quadrature. Such a quadrature can be done by mapping the tetrahedron to a hexahedron via the Duffy transformation. We show that for tensor product Gauss-Lobatto-Jacobi quadrature formulas with q+1=p+1 points in each direction and shape functions that are adapted to the quadrature formula, one again has discrete stability for the fully discrete p-FEM. The present error analysis complements the work [Eibner/Melenk 2005] for the p-FEM on triangles/tetrahedra where it is shown that by adapting the shape functions to the quadrature formula, the stiffness matrix can be set up in optimal complexity.

A couple of discretization schemes, based on an FE-like tessellation of the domain and polynomial reproducing, globally smooth shape functions, are considered and numerically explored to a limited extent. The first one among these is an existing scheme, the smooth DMS-FEM, that employs Delaunay triangulation or tetrahedralization (as approximate) towards discretizing the domain geometry employs triangular (tetrahedral) B-splines as kernel functions en route to the construction of polynomial reproducing functional approximations. In order to verify the numerical accuracy of the smooth DMS-FEM vis-à-vis the conventional FEM, a Mindlin-Reissner plate bending problem is numerically solved. Thanks to the higher order continuity in the functional approximant and the consequent removal of the jump terms in the weak form across inter-triangular boundaries, the numerical accuracy via the DMS-FEM approximation is observed to be higher than that corresponding to the conventional FEM. This advantage notwithstanding, evaluations of DMS-FEM based shape functions encounter singularity issues on the triangle vertices as well as over the element edges. This shortcoming is presently overcome through a new proposal that replaces the triangular B-splines by simplex splines, constructed over polygonal domains, as the kernel functions in the polynomial reproduction scheme. Following a detailed presentation of the issues related to its computational implementation, the new method is numerically explored with the results attesting to a higher attainable numerical accuracy in comparison with the DMS-FEM.

Orientador: Marco Lucio Bittencourt / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-08-10T12:57:32Z (GMT). No. of bitstreams: 1 Vazquez_ThaisGodoy_D.pdf: 11719751 bytes, checksum: c6d385d6a6414705c9f468358b8d3bea (MD5) Previous issue date: 2008 / Resumo: Este trabalho tem por objetivo principal o desenvolvimento de funções de interpolaçao e regras de integraçao tensorizaveis para o Metodo dos Elementos Finitos (MEF) de alta ordem hp, considerando os sistemas de referencias locais dos elementos. Para isso, primeiramente, determinam-se ponderaçoes especficas para as bases de funçoes de triangulos e tetraedros, formada pelo produto tensorial de polinomios de Jacobi, de forma a se obter melhor esparsidade e condicionamento das matrizes de massa e rigidez dos elementos. Alem disso, procuram-se novas funçoes de base para tornar as matrizes de massa e rigidez mais esparsas possiveis. Em seguida, escolhe-se os pontos de integraçao que otimizam o custo do calculo dos coeficientes das matrizes de massa e rigidez usando as regras de quadratura de Gauss-Jacobi, Gauss-Radau-Jacobi e Gauss-Lobatto-Jacobi. Por fim, mostra-se a construçao de uma base unidimensional nodal que permite obter uma matriz de rigidez praticamente diagonal para problemas de Poisson unidimensionais. Discute-se ainda extensoes para elementos bi e tridimensionais / Abstract: The main purpose of this work is the development of tensor-based interpolation functions and integration rules for the hp High-order Finite Element Method (FEM), considering the local reference systems of the elements. We first determine specific weights for the shape functions of triangles and tetrahedra, constructed by the tensorial product of Jacobi polynomials, aiming to obtain better sparsity and numerical conditioning for the mass and stiffness matrices of the elements. Moreover, new shape functions are proposed to obtain more sparse mass and stiffness matrices. After that, integration points are chosen that optimize the cost for the calculation of the coefficients of the mass and stiffness matrices using the rules of quadrature of Gauss-Jacobi, Gauss-Radau-Jacobi and Gauss-Lobatto-Jacobi. Finally, we construct an one-dimensional nodal shape function that obtains an almost diagonal stiffness matrix for the 1D Poisson problem. Extensions to two and three-dimensional elements are discussed. / Doutorado / Mecanica dos Sólidos e Projeto Mecanico / Doutor em Engenharia Mecânica

La découpe virtuelle d'objets déformables est au coeur de nombreuses applications pour la simulation interactive. Nous présentons un nouvel algorithme de remaillage permettant la simulation de découpes avec la méthode des éléments finis. Nous avons combiné notre algorithme avec la méthode du «snapping» déplaçant les noeuds la surface de coupe, pour des tétradres linéaires. Notre approche permet de maîtriser le nombre de noeuds et de la qualité numérique du maillage durant les coupes. Elle donne des résultats similaires pour les fonctions de forme quadratiques. Dans ce cadre, nous avons évalué le «snapping» pour simuler la fracture de surfaces triangulées. Nous avons appliqué nos résultats en 3D pour l'assistance aux gestes chirurgicaux, en étant les premiers présenter des résultats sur la détection de déchirures dans un flux vidéo monoculaire. La robustesse de notre algorithme et l'augmentation des structures internes des organes. / Virtual cutting of deformable objects is at the core of many applications in interactive simulation for computational medicine. We present a new remeshing algorithm simulating cuts based on the finite element method. For tetrahedral elements with linear shape tunctions, we combined our algorithm with the movement of the nodes on the cutting surface, called snapping. Our approach shows benefits when evaluating the impact of cuts on the number of nodes and the numerical quality of the mesh. The remeshing algorithm yields similar results for quadratic elements. However, the snapping of nodes entails higher challenges and has been evaluated on triangular elements, simulating fractures. For augmented reality applications, we are the first to present results on the detection of fractures, tearing and cutting from a monocular video stream. Examples in different contexts show the robustness of our algorithm and the augmentation of internai organ structures highlights the clinical interest.

In this paper we consider the p-FEM for elliptic boundary value problems on tetrahedral meshes where the entries of the stiffness matrix are evaluated by numerical quadrature. Such a quadrature can be done by mapping the tetrahedron to a hexahedron via the Duffy transformation. We show that for tensor product Gauss-Lobatto-Jacobi quadrature formulas with q+1=p+1 points in each direction and shape functions that are adapted to the quadrature formula, one again has discrete stability for the fully discrete p-FEM. The present error analysis complements the work [Eibner/Melenk 2005] for the p-FEM on triangles/tetrahedra where it is shown that by adapting the shape functions to the quadrature formula, the stiffness matrix can be set up in optimal complexity.

This thesis deals with smooth discretization schemes and inverse problems, the former used in efficient yet accurate numerical solutions to forward models required in turn to solve inverse problems. The aims of the thesis include, (i) development of a stabilization techniques for a class of forward problems plagued by unphysical oscillations in the response due to the presence of jumps/shocks/high gradients, (ii) development of a smooth hybrid discretization scheme that combines certain useful features of Finite Element (FE) and Mesh-Free (MF) methods and alleviates certain destabilizing factors encountered in the construction of shape functions using the polynomial reproduction method and, (iii) a first of its kind attempt at the joint inversion of both static and dynamic source parameters of the 2004 Sumatra-Andaman earthquake using tsunami sea level anomaly data. Following the introduction in Chapter 1 that motivates and puts in perspective the work done in later chapters, the main body of the thesis may be viewed as having two parts, viz., the first part constituting the development and use of smooth discretization schemes in the possible presence of destabilizing factors (Chapters 2 and 3) and the second part involving solution to the inverse problem of tsunami source recovery (Chapter 4). In the context of stability requirements in numerical solutions of practical forward problems, Chapter 2 develops a new stabilization scheme. It is based on a stochastic representation of the discretized field variables, with a view to reduce or even eliminate unphysical oscillations in the MF numerical simulations of systems developing shocks or exhibiting localized bands of extreme plastic deformation in the response. The origin of the stabilization scheme may be traced to nonlinear stochastic filtering and, consistent with a class of such filters, gain-based additive correction terms are applied to the simulated solution of the system, herein achieved through the Element-Free Galerkin (EFG) method, in order to impose a set of constraints that help arresting the spurious oscillations. The method is numerically illustrated through its application to a gradient plasticity model whose response is often characterized by a developing shear band as the external load is gradually increased. The potential of the method in stabilized yet accurate numerical simulations of such systems involving extreme gradient variations in the response is thus brought forth. Chapter 3 develops the MF-based discretization motif by balancing this with the widespread adoption of the FE method. Thus it concentrates on developing a 'hybrid' scheme that aims at the amelioration of certain destabilizing algorithmic issues arising from the necessary condition of moment matrix invertibility en route to the generation of smooth shape functions. It sets forth the hybrid discretization scheme utilizing bivariate simplex splines as kernels in a polynomial reproducing approach adopted over a conventional FE-like domain discretization based on Delaunay triangulation. Careful construction of the simplex spline knotset ensures the success of the polynomial reproduction procedure at all points in the domain of interest, a significant advancement over its precursor, the DMS-FEM. The shape functions in the proposed method inherit the global continuity ( C p 1 ) and local supports of the simplex splines of degree p . In the proposed scheme, the triangles comprising the domain discretization also serve as background cells for numerical integration which here are near-aligned to the supports of the shape functions (and their intersections), thus considerably ameliorating an oft-cited source of inaccuracy in the numerical integration of MF-based weak forms. Numerical experiments establish that the proposed method can work with lower order quadrature rules for accurate evaluation of integrals in the Galerkin weak form, a feature desiderated in solving nonlinear inverse problems that demand cost-effective solvers for the forward models. Numerical demonstrations of optimal convergence rates for a few test cases are given and the hybrid method is also implemented to compute crack-tip fields in a gradient-enhanced elasticity model. Chapter 4 attempts at the joint inversion of earthquake source parameters for the 2004 Sumatra-Andaman event from the tsunami sea level anomaly signals available from satellite altimetry. Usual inversion for earthquake source parameters incorporates subjective elements, e.g. a priori constraints, posing and parameterization, trial-and-error waveform fitting etc. Noisy and possibly insufficient data leads to stability and non-uniqueness issues in common deterministic inversions. A rational accounting of both issues favours a stochastic framework which is employed here, leading naturally to a quantification of the commonly overlooked aspects of uncertainty in the solution. Confluence of some features endows the satellite altimetry for the 2004 Sumatra-Andaman tsunami event with unprecedented value for the inversion of source parameters for the entire rupture duration. A nonlinear joint inversion of the slips, rupture velocities and rise times with minimal a priori constraints is undertaken. Large and hitherto unreported variances in the parameters despite a persistently good waveform fit suggest large propagation of uncertainties and hence the pressing need for better physical models to account for the defect dynamics and massive sediment piles. Chapter 5 concludes the work with pertinent comments on the results obtained and suggestions for future exploration of some of the schemes developed here.

The presented dissertation thesis is focused, as the title suggests, on the analysis of stress state aspects of quasi-brittle fracture. That means the fracture of composite materials with cement matrix (such as concrete, mortar, plaster, etc.), ceramics and other composites. Used methods are based on the theory of multi-parameter linear elastic fracture mechanics, which highlights the importance of considering of several initial terms of Williams power series, approximating the stress and displacement fields in a cracked body, within conducted fracture analyses. Determination of values of coefficients of terms of this series, recalculated into the shape functions serving in most of the conducted stress state analyses, is performed via the so called over-deterministic method. Another tool for the problem solving is nonlinear fracture mechanics, represented primarily by the cohesive crack model, namely the crack band model implemented in the used ATENA software. For the backward reconstruction of stress field in the cracked bodies the application ReFraPro is used. The analytical part deals with various aspects of wedge-splitting test – from the boundary conditions, though various possibilities of nodal selection (required as input variables for the over-deterministic method) up to the advanced (automated) analysis of numerical model. Special chapter includes atypical test specimens designed for adjusting of various levels of constraint of stress and deformation at the propagating crack tip. The study of this geometry and also the subsequent detail analysis reveals important information for real experiments. Backward reconstruction of stress field presents analysis on suitable possibilities of nodal selections as inputs into the procedure of approximation of the crack tip fields and answers the question of the necessity of application of the multi-parameter linear elastic fracture mechanics for certain fracture analyses of specimens from quasi-brittle materials. The th