Dg And Hdg Methods For Curved Structures

We introduce and analyze discontinuous Galerkin methods for a Naghdi type arch model. We prove that, when the numerical traces are properly chosen, the methods display optimal convergence uniformly with respect to the thickness of the arch. These methods are thus free from membrane and shear locking. We also prove that, when polynomials of degree $k$ are used, {\em all} the numerical traces superconverge with a rate of order h 2k+1. Based on the superconvergent phenomenon and we show how to post-process them in an element-by-element fashion to obtain a far better approximation. Indeed, we prove that, if polynomials of degree k are used, the post-processed approximation converges with order 2k+1 in the L2-norm throughout the domain. This has to be contrasted with the fact that before post-processing, the approximation converges with order k+1 only. Moreover, we show that this superconvergence property does not deteriorate as the thickness of the arch becomes extremely small. Since the DG methods suffer from too many degree of freedoms we introduce and analyze a class of hybridizable discontinuous Galerkin (HDG) methods for Naghdi arches. The main feature of these methods is that they can be implemented in an efficient way through a hybridization procedure which reduces the globally coupled unknowns to approximations to the transverse and tangential displacement and bending moment at the element boundaries. The error analysis of the methods is based on the use of a projection especially designed to fit the structure of the numerical traces of the method. This property allows to prove in a very concise manner that the projection of the errors is bounded in terms of the distance between the exact solution and its projection. The study of the influence of the stabilization function on the approximation is then reduced to the study of how they affect the approximation properties of the projection in a single element. Consequently, we prove that HDG methods have the same result as DG methods. At the end of the thesis, we talk a little bit of shell problems

Exploiting Superconvergence Through Smoothness-Increasing Accuracy-Conserving (SIAC) Filtering

There has been much work in the area of superconvergent error analysis for finite element and discontinuous Galerkin (DG) methods. The property of superconvergence leads to the question of how to exploit this information in a useful manner, mainly through superconvergence extraction. There are many methods used for superconvergence extraction such as projection, interpolation, patch recovery and B-spline convolution filters. This last method falls under the class of Smoothness-Increasing Accuracy-Conserving (SIAC) filters. It has the advantage of improving both smoothness and accuracy of the approximation. Specifically, for linear hyperbolic equations it can improve the order of accuracy of a DG approximation from k + 1 to 2k + 1, where k is the highest degree polynomial used in the approximation, and can increase the smoothness to k − 1. In this article, we discuss the importance of overcoming the mathematical barriers in making superconvergence extraction techniques useful for applications, specifically focusing on SIAC filtering

Interlaminar Stresses by Refined Beam Theories and the Sinc Method Based on Interpolation of Highest Derivative

Computation of interlaminar stresses from the higher-order shear and normal deformable beam theory and the refined zigzag theory was performed using the Sinc method based on Interpolation of Highest Derivative. The Sinc method based on Interpolation of Highest Derivative was proposed as an efficient method for determining through-the-thickness variations of interlaminar stresses from one- and two-dimensional analysis by integration of the equilibrium equations of three-dimensional elasticity. However, the use of traditional equivalent single layer theories often results in inaccuracies near the boundaries and when the lamina have extremely large differences in material properties. Interlaminar stresses in symmetric cross-ply laminated beams were obtained by solving the higher-order shear and normal deformable beam theory and the refined zigzag theory with the Sinc method based on Interpolation of Highest Derivative. Interlaminar stresses and bending stresses from the present approach were compared with a detailed finite element solution obtained by ABAQUS/Standard. The results illustrate the ease with which the Sinc method based on Interpolation of Highest Derivative can be used to obtain the through-the-thickness distributions of interlaminar stresses from the beam theories. Moreover, the results indicate that the refined zigzag theory is a substantial improvement over the Timoshenko beam theory due to the piecewise continuous displacement field which more accurately represents interlaminar discontinuities in the strain field. The higher-order shear and normal deformable beam theory more accurately captures the interlaminar stresses at the ends of the beam because it allows transverse normal strain. However, the continuous nature of the displacement field requires a large number of monomial terms before the interlaminar stresses are computed as accurately as the refined zigzag theory

Meshless methods for shear-deformable beams and plates based on mixed weak forms

Thin structural theories such as the shear-deformable Timoshenko beam and Reissner-Mindlin plate theories have seen wide use throughout engineering practice to simulate the response of structures with planar dimensions far larger than their thickness dimension. Meshless methods have been applied to construct numerical methods to solve the shear deformable theories. Similarly to the finite element method, meshless methods must be carefully designed to overcome the well-known shear-locking problem. Many successful treatments of shear-locking in the finite element literature are constructed through the application of a mixed weak form. In the mixed weak form the shear stresses are treated as an independent variational quantity in addition to the usual displacement variables. We introduce a novel hybrid meshless-finite element formulation for the Timoshenko beam problem that converges to the stable first-order/zero-order finite element method in the local limit when using maximum entropy meshless basis functions. The resulting formulation is free from the effects shear-locking. We then consider the Reissner-Mindlin plate problem. The shear stresses can be identified as a vector field belonging to the Sobelov space with square integrable rotation, suggesting the use of rotated Raviart-Thomas-Nedelec elements of lowest-order for discretising the shear stress field. This novel formulation is again free from the effects of shear-locking. Finally we consider the construction of a generalised displacement method where the shear stresses are eliminated prior to the solution of the final linear system of equations. We implement an existing technique in the literature for the Stokes problem called the nodal volume averaging technique. To ensure stability we split the shear energy between a part calculated using the displacement variables and the mixed variables resulting in a stabilised weak form. The method then satisfies the stability conditions resulting in a formulation that is free from the effects of shear-locking.Open Acces

HDGlab: An Open-Source Implementation of the Hybridisable Discontinuous Galerkin Method in MATLAB

This paper presents HDGlab, an open source MATLAB implementation of the hybridisable discontinuous Galerkin (HDG) method. The main goal is to provide a detailed description of both the HDG method for elliptic problems and its implementation available in HDGlab. Ultimately, this is expected to make this relatively new advanced discretisation method more accessible to the computational engineering community. HDGlab presents some features not available in other implementations of the HDG method that can be found in the free domain. First, it implements high-order polynomial shape functions up to degree nine, with both equally-spaced and Fekete nodal distributions. Second, it supports curved isoparametric simplicial elements in two and three dimensions. Third, it supports non-uniform degree polynomial approximations and it provides a flexible structure to devise degree adaptivity strategies. Finally, an interface with the open-source high-order mesh generator Gmsh is provided to facilitate its application to practical engineering problems

MATHICSE Technical Report : Fast and accurate elastic analysis of laminated composite plates via isogeometric collocation and an equilibrium-based stress recovery approach

A novel approach which combines isogeometric collocation and an equilibriumbased stress recovery technique is applied to analyze laminated composite plates. Isogeometric collocation is an appealing strong form alternative to standard Galerkin approaches, able to achieve high order convergence rates coupled with a significantly reduced computational cost. Laminated composite plates are herein conveniently modeled considering only one element through the thickness with homogenized material properties. This guarantees accurate results in terms of displacements and in-plane stress components. To recover an accurate out-of-plane stress state, equilibrium is imposed in strong form as a post-processing correction step, which requires the shape functions to be highly continuous. This continuity demand is fully granted A novel approach which combines isogeometric collocation and an equilibriumbased stress recovery technique is applied to analyze laminated composite plates. Isogeometric collocation is an appealing strong form alternative to standard Galerkin approaches, able to achieve high order convergence rates coupled with a significantly reduced computational cost. Laminated composite plates are herein conveniently modeled considering only one element through the thickness with homogenized material properties. This guarantees accurate results in terms of displacements and in-plane stress components. To recover an accurate out-of-plane stress state, equilibrium is imposed in strong form as a post-processing correction step, which requires the shape functions to be highly continuous. This continuity demand is fully granted