54 research outputs found
The BEM with graded meshes for the electric field integral equation on polyhedral surfaces
We consider the variational formulation of the electric field integral
equation on a Lipschitz polyhedral surface . We study the Galerkin
boundary element discretisations based on the lowest-order Raviart-Thomas
surface elements on a sequence of anisotropic meshes algebraically graded
towards the edges of . We establish quasi-optimal convergence of
Galerkin solutions under a mild restriction on the strength of grading. The key
ingredient of our convergence analysis are new componentwise stability
properties of the Raviart-Thomas interpolant on anisotropic elements
The hp-BEM with quasi-uniform meshes for the electric field integral equation on polyhedral surfaces: a priori error analysis
This paper presents an a priori error analysis of the hp-version of the
boundary element method for the electric field integral equation on a piecewise
plane (open or closed) Lipschitz surface. We use H(div)-conforming
discretisations with Raviart-Thomas elements on a sequence of quasi-uniform
meshes of triangles and/or parallelograms. Assuming the regularity of the
solution to the electric field integral equation in terms of Sobolev spaces of
tangential vector fields, we prove an a priori error estimate of the method in
the energy norm. This estimate proves the expected rate of convergence with
respect to the mesh parameter h and the polynomial degree p
A note on the polynomial approximation of vertex singularities in the boundary element method in three dimensions
We study polynomial approximations of vertex singularities of the type on three-dimensional surfaces. The analysis focuses on the case when . This assumption is a minimum requirement to guarantee that the above singular function is in the energy space for boundary integral equations with hypersingular operators. Thus, the approximation results for such singularities are needed for the error analysis of boundary element methods on piecewise smooth surfaces. Moreover, to our knowledge, the approximation of strong singularities () by high-order polynomials is missing in the existing literature. In this note we prove an estimate for the error of polynomial approximation of the above vertex singularities on quasi-uniform meshes discretising a polyhedral surface. The estimate gives an upper bound for the error in terms of the mesh size and the polynomial degree
Optimal error estimation for H(curl)-conforming p-interpolation in two dimensions
In this paper we prove an optimal error estimate for the H(curl)-conforming
projection based p-interpolation operator introduced in [L. Demkowicz and I.
Babuska, p interpolation error estimates for edge finite elements of variable
order in two dimensions, SIAM J. Numer. Anal., 41 (2003), pp. 1195-1208]. This
result is proved on the reference element (either triangle or square) K for
regular vector fields in H^r(curl,K) with arbitrary r>0. The formulation of the
result in the H(div)-conforming setting, which is relevant for the analysis of
high-order boundary element approximations for Maxwell's equations, is provided
as well
SOLID-SHELL FINITE ELEMENT MODELS FOR EXPLICIT SIMULATIONS OF CRACK PROPAGATION IN THIN STRUCTURES
Crack propagation in thin shell structures due to cutting is conveniently simulated
using explicit finite element approaches, in view of the high nonlinearity of the problem. Solidshell
elements are usually preferred for the discretization in the presence of complex material
behavior and degradation phenomena such as delamination, since they allow for a correct
representation of the thickness geometry. However, in solid-shell elements the small thickness
leads to a very high maximum eigenfrequency, which imply very small stable time-steps. A new
selective mass scaling technique is proposed to increase the time-step size without affecting
accuracy. New âdirectionalâ cohesive interface elements are used in conjunction with selective
mass scaling to account for the interaction with a sharp blade in cutting processes of thin ductile
shells
Schnelle Löser fĂŒr Partielle Differentialgleichungen
The workshop Schnelle Löser fĂŒr partielle Differentialgleichungen, organised by Randolph E. Bank (La Jolla), Wolfgang Hackbusch (Leipzig), and Gabriel Wittum (Frankfurt am Main), was held May 22ndâMay 28th, 2011. This meeting was well attended by 54 participants with broad geographic representation from 7 countries and 3 continents. This workshop was a nice blend of researchers with various backgrounds
Failure processes in soft and quasi-brittle materials with nonhomogeneous microstructures
Material failure pervades the fields of materials science and engineering; it occurs at various scales and in various contexts. Understanding the mechanisms by which a material fails can lead to advancements in the way we design and build the world around us. For example, in structural engineering, understanding the fracture of concrete and steel can lead to improved structural systems and safer designs; in geological engineering, understanding the fracture of rock can lead to increased efficiency in oil and gas extraction; and in biological engineering, understanding the fracture of bone can lead to improvements in the design of bio-composites and medical implants. In this thesis, we numerically investigate a wide spectrum of failure behavior; in soft and quasi-brittle materials with nonhomogeneous microstructures considering a statistical distribution of material properties.
The first topic we investigate considers the influence of interfacial interactions on the macroscopic constitutive response of particle reinforced elastomers. When a particle is embedded into an elastomer, the polymer chains in the elastomer tend to adsorb (or anchor) onto the surface of the particle; creating a region in the vicinity of each particle (often referred to as an interphase) with distinct properties from those in the bulk elastomer. This interphasial region has been known to exist for many decades, but is primarily omitted in computational investigations of such composites. In this thesis, we present an investigation into the influence of interphases on the macroscopic constitutive response of particle filled elastomers undergoing large deformations. In addition, at large deformations, a localized region of failure tends to accumulate around inclusions. To capture this localized region of failure (often referred to as interfacial debonding), we use cohesive zone elements which follow the Park-Paulino-Roesler traction-separation relation. To account for friction, we present a new, coupled cohesive-friction relation and detail its formulation and implementation. In the process of this investigation, we developed a small library of cohesive elements for use with a commercially available finite element analysis software package.
Additionally, in this thesis, we present a series of methods for reducing mesh dependency in two-dimensional dynamic cohesive fracture simulations of quasi-brittle materials. In this setting, cracks are only permitted to propagate along element facets, thus a poorly designed discretization of the problem domain can introduce artifacts into the fracture behavior. To reduce mesh induced artifacts, we consider unstructured polygonal finite elements. A randomly-seeded polygonal mesh leads to an isotropic discretization of the problem domain, which does not bias the direction of crack propagation. However, polygonal meshes tend to limit the possible directions a crack may travel at each node, making this discretization a poor candidate for dynamic cohesive fracture simulations. To alleviate this problem, we propose two new topological operators. The first operator we propose is adaptive element-splitting, and the second is adaptive mesh refinement. Both operators are designed to improve the ability of unstructured polygonal meshes to capture crack patterns in dynamic cohesive fracture simulations. However, we demonstrate that element-splitting is more suited to pervasive fracture problems, whereas, adaptive refinement is more suited to problems exhibiting a dominant crack.
Finally, we investigate the use of geometric and constitutive design features to regularize pervasive fragmentation behavior in three-dimensions. Throughout pervasive fracture simulations, many cracks initiate, propagate, branch and coalesce simultaneously. Because of the cohesive element method's unique framework, this behavior can be captured in a regularized manner. In this investigation, unstructuring techniques are used to introduce randomness into a numerical model. The behavior of quasi-brittle materials undergoing pervasive fracture and fragmentation is then examined using three examples. The examples are selected to investigate some of the significant factors influencing pervasive fracture and fragmentation behavior; including, geometric features, loading conditions, and material gradation
Mathematical foundations of adaptive isogeometric analysis
This paper reviews the state of the art and discusses recent developments in
the field of adaptive isogeometric analysis, with special focus on the
mathematical theory. This includes an overview of available spline technologies
for the local resolution of possible singularities as well as the
state-of-the-art formulation of convergence and quasi-optimality of adaptive
algorithms for both the finite element method (FEM) and the boundary element
method (BEM) in the frame of isogeometric analysis (IGA)
- âŠ