A discontinuous Galerkin formulation of non‐linear Kirchhoff–Love shells

Discontinuous Galerkin (DG) methods provide a means of weakly enforcing the continuity of the unknown‐field derivatives and have particular appeal in problems involving high‐order derivatives. This feature has previously been successfully exploited (Comput. Methods Appl. Mech. Eng. 2008; 197 :2901–2929) to develop a formulation of linear Kirchhoff–Love shells considering only the membrane and bending responses. In this proposed one‐field method—the displacements are the only unknowns, while the displacement field is continuous, the continuity in the displacement derivative between two elements is weakly enforced by recourse to a DG formulation. It is the purpose of the present paper to extend this formulation to finite deformations and non‐linear elastic behaviors. While the initial linear formulation was relying on the direct linear computation of the effective membrane stress and effective bending couple‐stress from the displacement field at the mid‐surface of the shell, the non‐linear formulation considered implies the evaluation of the general stress tensor across the shell thickness, leading to a reformulation of the internal forces of the shell. Nevertheless, since the interface terms resulting from the discontinuous Galerkin method involve only the resultant couple‐stress at the edges of the shells, the extension to non‐linear deformations is straightforward. Copyright © 2008 John Wiley & Sons, Ltd.     

Shear‐flexible subdivision shells

We present a new shell model and an accompanying discretisation scheme that is suitable for thin and thick shells. The deformed configuration of the shell is parameterised using the mid‐surface position vector and an additional shear vector for describing the out‐of‐plane shear deformations. In the limit of vanishing thickness, the shear vector is identically zero and the Kirchhoff–Love model is recovered. Importantly, there are no compatibility constraints to be satisfied by the shape functions used for discretising the mid‐surface and the shear vector. The mid‐surface has to be interpolated with smooth C¹‐continuous shape functions, whereas the shear vector can be interpolated with C⁰‐continuous shape functions. In the present paper, the mid‐surface as well as the shear vector are interpolated with smooth subdivision shape functions. The resulting finite elements are suitable for thin and thick shells and do not exhibit shear locking. The good performance of the proposed formulation is demonstrated with a number of linear and geometrically non‐linear plate and shell examples. Copyright © 2012 John Wiley & Sons, Ltd.     

Elasto‐plastic finite element analysis of shells with damage due to microvoids

This paper presents a non‐linear finite element analysis for the elasto‐plastic behaviour of thick/thin shells and plates with large rotations and damage effects. The refined shell theory given by Voyiadjis and Woelke (Int. J. Solids Struct. 2004; 41 :3747–3769) provides a set of shell constitutive equations. Numerical implementation of the shell theory leading to the development of the C⁰ quadrilateral shell element (Woelke and Voyiadjis, Shell element based on the refined theory for thick spherical shells. 2006, submitted) is used here as an effective tool for a linear elastic analysis of shells. The large rotation elasto‐plastic model for shells presented by Voyiadjis and Woelke (General non‐linear finite element analysis of thick plates and shells. 2006, submitted) is enhanced here to account for the damage effects due to microvoids, formulated within the framework of a micromechanical damage model. The evolution equation of the scalar porosity parameter as given by Duszek‐Perzyna and Perzyna (Material Instabilities: Theory and Applications, ASME Congress, Chicago, AMD‐Vol. 183/MD‐50, 9–11 November 1994; 59–85) is reduced here to describe the most relevant damage effects for isotropic plates and shells, i.e. the growth of voids as a function of the plastic flow. The anisotropic damage effects, the influence of the microcracks and elastic damage are not considered in this paper. The damage modelled through the evolution of porosity is incorporated directly into the yield function, giving a generalized and convenient loading surface expressed in terms of stress resultants and stress couples. A plastic node method (Comput. Methods Appl. Mech. Eng. 1982; 34 :1089–1104) is used to derive the large rotation, elasto‐plastic‐damage tangent stiffness matrix. Some of the important features of this paper are that the elastic stiffness matrix is derived explicitly, with all the integrals calculated analytically (Woelke and Voyiadjis, Shell element based on the refined theory for thick spherical shells. 2006, submitted). In addition, a non‐layered model is adopted in which integration through the thickness is not necessary. Consequently, the elasto‐plastic‐damage stiffness matrix is also given explicitly and numerical integration is not performed. This makes this model consistent mathematically, accurate for a variety of applications and very inexpensive from the point of view of computer power and time. Copyright © 2006 John Wiley & Sons, Ltd.     

Isogeometric FEM‐BEM coupled structural‐acoustic analysis of shells using subdivision surfaces

We introduce a coupled finite and boundary element formulation for acoustic scattering analysis over thin‐shell structures. A triangular Loop subdivision surface discretisation is used for both geometry and analysis fields. The Kirchhoff‐Love shell equation is discretised with the finite element method and the Helmholtz equation for the acoustic field with the boundary element method. The use of the boundary element formulation allows the elegant handling of infinite domains and precludes the need for volumetric meshing. In the present work, the subdivision control meshes for the shell displacements and the acoustic pressures have the same resolution. The corresponding smooth subdivision basis functions have the C¹ continuity property required for the Kirchhoff‐Love formulation and are highly efficient for the acoustic field computations. We verify the proposed isogeometric formulation through a closed‐form solution of acoustic scattering over a thin‐shell sphere. Furthermore, we demonstrate the ability of the proposed approach to handle complex geometries with arbitrary topology that provides an integrated isogeometric design and analysis workflow for coupled structural‐acoustic analysis of shells.     

Quasi‐static crack propagation in plane and plate structures using set‐valued traction‐separation laws

We introduce a numerical technique to model set‐valued traction‐separation laws in plate bending and also plane crack propagation problems. By using of recent developments in thin (Kirchhoff–Love) shell models and the extended finite element method, a complete and accurate algorithm for the cohesive law is presented and is used to determine the crack path. The cohesive law includes softening and unloading to origin, adhesion and contact. Pure debonding and contact are obtained as particular (degenerate) cases. A smooth root‐finding algorithm (based on the trust‐region method) is adopted. A step‐driven algorithm is described with a smoothed law which can be made arbitrarily close to the exact non‐smooth law. In the examples shown the results were found to be step‐size insensitive and accurate. In addition, the method provides the crack advance law, extracted from the cohesive law and the absence of stress singularity at the tip. Copyright © 2007 John Wiley & Sons, Ltd.     

Non‐linear analysis of shells with arbitrary evolving cracks using XFEM

A new formulation and numerical procedures are developed for the analysis of arbitrary crack propagation in shells using the extended finite element method. The method is valid for completely non‐linear problems. Through‐the‐thickness cracks in sandwich shells are considered. An exact shell kinematics is presented, and a new enrichment of the rotation field is proposed which satisfies the director inextensibility condition. To avoid locking, an enhanced strain formulation is proposed for the 4‐node cracked shell element. A finite strain plane stress constitutive model based on the logarithmic corotational rate is employed. A cohesive zone model is introduced which embodies the special characteristics of the shell kinematics. Stress intensity factors are calculated for selected problems and crack propagation problems are solved. Copyright © 2004 John Wiley & Sons, Ltd.     

A damage to crack transition model accounting for stress triaxiality formulated in a hybrid nonlocal implicit discontinuous Galerkin‐cohesive band model framework

Modelling the entire ductile fracture process remains a challenge. On the one hand, continuous damage models succeed in capturing the initial diffuse damage stage but are not able to represent discontinuities or cracks. On the other hand, discontinuous methods, as the cohesive zones, which model the crack propagation behaviour, are suited to represent the localised damaging process. However, they are unable to represent diffuse damage. Moreover, most of the cohesive models do not capture triaxiality effect. In this paper, the advantages of the two approaches are combined in a single damage to crack transition framework. In a small deformation setting, a nonlocal elastic damage model is associated with a cohesive model in a discontinuous Galerkin finite element framework. A cohesive band model is used to naturally introduce a triaxiality‐dependent behaviour inside the cohesive law. Practically, a numerical thickness is introduced to recover a 3D state, mandatory to incorporate the in‐plane stretch effects. This thickness is evaluated to ensure the energy consistency of the method and is not a new numerical parameter. The traction‐separation law is then built from the underlying damage model. The method is numerically shown to capture the stress triaxiality effect on the crack initiation and propagation.     

New development in extended finite element modeling of large elasto‐plastic deformations

This paper presents new achievements in the extended finite element modeling of large elasto‐plastic deformation in solid problems. The computational technique is presented based on the extended finite element method (X‐FEM) coupled with the Lagrangian formulation in order to model arbitrary interfaces in large deformations. In X‐FEM, the material interfaces are represented independently of element boundaries, and the process is accomplished by partitioning the domain with some triangular sub‐elements whose Gauss points are used for integration of the domain of elements. The large elasto‐plastic deformation formulation is employed within the X‐FEM framework to simulate the non‐linear behavior of materials. The interface between two bodies is modeled by using the X‐FEM technique and applying the Heaviside‐ and level‐set‐based enrichment functions. Finally, several numerical examples are analyzed, including arbitrary material interfaces, to demonstrate the efficiency of the X‐FEM technique in large plasticity deformations. Copyright © 2008 John Wiley & Sons, Ltd.     

Three-dimensional numerical simulations of dynamic fracture in silicon carbide reinforced aluminum

The three-point bending test by Kolsky-bar apparatus is a convenient technique to test the dynamic fracture properties of materials. This paper presents detailed three-dimensional finite element simulations of a silicon particle reinforced aluminum (SiC_p/Al) experiment (Li et al., [Proceedings of the US Army Symposium on Solid Mechanics]. In the simulations, the interaction between the input bar and the specimen is modeled by coupled boundary conditions. The material model includes large plastic deformations, strain-hardening and strain-rate hardening mechanisms. Furthermore, crack initiation and propagation processes are simulated by a cohesive element model. The simulation results quantitatively agree with the experimental measurements on three fronts: (1) the structural response of the specimen, (2) the time of unstable crack propagation, and (3) the local deformations at the crack-tip zone. The simulations reveal crack propagation characteristics, including crack-tip plastic deformation, crack front curving, and crack velocity profile. The effectiveness of Kolsky-bar type fracture tests is verified. It is shown that a rate-independent cohesive model can describe the complicated dynamic elastic–plastic fracture process in the SiC_p/Al material.     

Lebensdauerberechnung gekerbter Bauteile auf Basis der elastisch‐plastischen Schwingbruchmechanik

Fatigue life calculation of notched components based on the elastic‐plastic fatigue fracture mechanics The life of notched components is subdivided into the pre‐crack, or crack‐initiation, and crack propagation phases within and outside notch area. It is known that a major factor governing the service life of notched components under cyclic loading is fatigue crack growth in notches. Therefore a uniform elastic‐plastic crack growth model, based on the J‐Integral, was developed which especially considers the crack opening and closure behaviour and the effect of residual stresses for the determination of crack initiation and propagation lives for cracks in notches under constant and variable‐amplitude loading. The crack growth model will be introduced and verified by experiments.     

An explicit discontinuous Galerkin method for non‐linear solid dynamics: Formulation,parallel implementation and scalability properties

An explicit‐dynamics spatially discontinuous Galerkin (DG) formulation for non‐linear solid dynamics is proposed and implemented for parallel computation. DG methods have particular appeal in problems involving complex material response, e.g. non‐local behavior and failure, as, even in the presence of discontinuities, they provide a rigorous means of ensuring both consistency and stability. In the proposed method, these are guaranteed: the former by the use of average numerical fluxes and the latter by the introduction of appropriate quadratic terms in the weak formulation. The semi‐discrete system of ordinary differential equations is integrated in time using a conventional second‐order central‐difference explicit scheme. A stability criterion for the time integration algorithm, accounting for the influence of the DG discretization stability, is derived for the equivalent linearized system. This approach naturally lends itself to efficient parallel implementation. The resulting DG computational framework is implemented in three dimensions via specialized interface elements. The versatility, robustness and scalability of the overall computational approach are all demonstrated in problems involving stress‐wave propagation and large plastic deformations. Copyright © 2007 John Wiley & Sons, Ltd.     

A rational elasto‐plastic spatially curved thin‐walled beam element

Torsion is one of the primary actions in members curved in space, and so an accurate spatially curved‐beam element needs to be able to predict the elasto‐plastic torsional behaviour of such members correctly. However, there are two major difficulties in most existing finite thin‐walled beam elements, such as in ABAQUS and ANSYS, which may lead to incorrect predictions of the elasto‐plastic behaviour of members curved in space. Firstly, the integration sample point scheme cannot capture the shear strain and stress information resulting from uniform torsion. Secondly, the higher‐order twists are ignored which leads to loss of the significant effects of Wagner moments on the large twist torsional behaviour. In addition, the initial geometric imperfections and residual stresses are significant for the elasto‐plastic behaviour of members curved in space. Many existing finite thin‐walled beam element models do not provide facilities to deal with initial geometric imperfections. Although ABAQUS and ANSYS have facilities for the input of residual stresses as initial stresses, they cannot describe the complicated distribution patterns of residual stresses in thin‐walled members. Furthermore, external loads and elastic restraints may be applied remote from shear centres or centroids. The effects of the load (and restraint) positions are important, but are not considered in many beam elements. This paper presents an elasto‐plastic spatially curved element with arbitrary thin‐walled cross‐sections that can correctly capture the uniform shear strain and stress information for integration, and includes initial geometric imperfections, residual stresses and the effects of the load and restraint positions. The element also includes elastic restraints and supports, which have to be modelled separately as spring elements in some other finite thin‐walled beam elements. Comparisons with existing experimental and analytical results show that the elasto‐plastic spatially curved‐beam element is accurate and efficient. Copyright © 2006 John Wiley & Sons, Ltd.     

Ductile damage analysis of elasto‐plastic shells at large inelastic strains

The objective of this contribution is to model ductile damage phenomena under consideration of large inelastic strains, to couple the corresponding constitutive law with a multi‐layer shell kinematics and to give finally an adequate finite element formulation. An elastic–plastic constitutive law is formulated by using a spatial hyperelasto‐plastic formulation based on the multiplicative decomposition of the deformation gradient. To include isotropic ductile damage the continuum damage model of Rousselier is modified so as to consider large strains and additionally extended by various void nucleation and macro‐crack criteria. In order to achieve numerical efficiency, elastic strains are supposed to be sufficiently small providing a numerical effective integration based on the backward Euler rule. Finite element formulation is enriched by means of the enhanced strain concept. Thus the well‐known deficiencies due to incompressible deformations and the inclusion of transverse strains are avoided. Several examples are given to demonstrate the performance of the algorithms developed concerning large inelastic strains of shells and ductile damage phenomena. Copyright © 2000 John Wiley & Sons, Ltd.     

Nonlinear rotation‐free three‐node shell finite element formulation

A new triangular thin‐shell finite element formulation is presented, which employs only translational degrees of freedom. The formulation allows for large deformations, and it is based on the nonlinear Kirchhoff thin‐shell theory. A number of static and dynamic test problems are considered for which analytical or benchmark solutions exist. Comparisons between the predictions of the new model and these solutions show that the new model accurately reproduces complex nonlinear analytical solutions as well as solutions obtained using existing, more complex finite element formulations. Copyright © 2013 John Wiley & Sons, Ltd.     

Mesh‐independent matrix cracking and delamination modeling in laminated composites

The initiation and evolution of transverse matrix cracks and delaminations are predicted within a mesh‐independent cracking (MIC) framework. MIC is a regularized extended finite element method (x‐FEM) that allows the insertion of cracks in directions that are independent of the mesh orientation. The Heaviside step function that is typically used to introduce a displacement discontinuity across a crack surface is replaced by a continuous function approximated by using the original displacement shape functions. Such regularization allows the preservation of the Gaussian integration schema regardless of the enrichment required to model cracking in an arbitrary direction. The interaction between plies is anchored on the integration point distribution, which remains constant through the entire simulation. Initiation and propagation of delaminations between plies as well as intra‐ply MIC opening is implemented by using a mixed‐mode cohesive formulation in a fully three‐dimensional model that includes residual thermal stresses. The validity of the proposed methodology was tested against a variety of problems ranging from simple evolution of delamination from existing transverse cracks to strength predictions of complex laminates withouttextita priori knowledge of damage location or initiation. Good agreement with conventional numerical solutions and/or experimental data was observed in all the problems considered. Published 2011. This article is a US Government work and is in the public domain in the USA.     

Experimentelle Charakterisierung und zweidimensionale Simulation der mikrostrukturbestimmten Ermüdungsrissausbreitung in einer β‐Titanlegierung

In this project the initiation and propagation of short fatigue cracks in the metastable β‐titanium alloy TIMETAL®LCB is investigated. By means of an interferometric strain/displacement gauge system (ISDG) to measure the crack opening displacement (COD) and the electron back scattered diffraction technique (EBSD) to determine the orientation of individual grains the microstructural influence on short crack initiation and growth can be characterized. Finite element calculations show a high influence of the elastic anisotropy on the initiation sites of cracks. Crack propagation takes place transgranulary along slip planes as well as intergranulary along grain boundaries. The crack growth rate depends strongly on the active mechanism at the crack tip which in turn is influenced by crack length, the applied stress and the orientation of the grains involved. The value of the steady state crack closure stress changes from a positive value at low applied stresses (roughness induced) to a negative one at higher applied stresses (due to plastic deformations at the crack tip). The crack growth simulation is realised by a two‐dimensional boundary element technique, which contains the ideas of Navarro und de los Rios. The model includes the sequence of the applied stress amplitude as well as the experimental measured roughness induced crack closure.     

Weight‐adjusted discontinuous Galerkin methods: Matrix‐valued weights and elastic wave propagation in heterogeneous media

Weight‐adjusted inner products are easily invertible approximations to weighted L² inner products. These approximations can be paired with a discontinuous Galerkin (DG) discretization to produce a time‐domain method for wave propagation which is low storage, energy stable, and high‐order accurate for arbitrary heterogeneous media and curvilinear meshes. In this work, we extend weight‐adjusted DG methods to the case of matrix‐valued weights, with the linear elastic wave equation as an application. We present a DG formulation of the symmetric form of the linear elastic wave equation, with upwind‐like dissipation incorporated through simple penalty fluxes. A semidiscrete convergence analysis is given, and numerical results confirm the stability and high‐order accuracy of weight‐adjusted DG for several problems in elastic wave propagation.     

A statically determinate truss model for thin shells: Two-surface analysis (bending theory)

This paper presents a simple scheme for analysing shells subject to the Kirchhoff–Love hypothesis. The method combines a previously developed, statically determinate truss model with the static-geometric analogy of elastic thin shell theory. In this way, both in-plane and flexural mid-surface actions are obtained by simply repeating the elementary computations performed on an essentially ‘membrane’ truss.     

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, quasistatic, elastic limit. Copyright © 2013 John Wiley & Sons, Ltd.