This work is an overview of available theories and finite elements that have been developed for multilayered, anisotropic, composite plate and shell structures. Although a comprehensive description of several techniques and approaches is given, most of this paper has been devoted to the so called axiomatic theories and related finite element implementations. Most of the theories and finite elements that have been proposed over the last thirty years are in fact based on these types of approaches. The paper has been divided into three parts. Part I, has been devoted to the description of possible approaches to plate and shell structures: 3D approaches, continuum based methods, axiomatic and asymptotic two-dimensional theories, classical and mixed formulations, equivalent single layer and layer wise variable descriptions are considered (the number of the unknown variables is considered to be independent of the number of the constitutive layers in the equivalent single layer case). Complicating effects that have been introduced by anisotropic behavior and layered constructions, such as high transverse deformability, zig-zag effects and interlaminar continuity, have been discussed and summarized by the acronimC
 -Requirements. Two-dimensional theories have been dealt with in Part II. Contributions based on axiomatic, asymtotic and continuum based approaches have been overviewed. Classical theories and their refinements are first considered. Both case of equivalent single-layer and layer-wise variables descriptions are discussed. The so-called zig-zag theories are then discussed. A complete and detailed overview has been conducted for this type of theory which relies on an approach that is entirely originated and devoted to layered constructions. Formulas and contributions related to the three possible zig-zag approaches, i.e. Lekhnitskii-Ren, Ambartsumian-Whitney-Rath-Das, Reissner-Murakami-Carrera ones have been presented and overviewed, taking into account the findings of a recent historical note provided by the author. Finite Element FE implementations are examined in Part III. The possible developments of finite elements for layered plates and shells are first outlined. FEs based on the theories considered in Part II are discussed along with those approaches which consist of a specific application of finite element techniques, such as hybrid methods and so-called global/local techniques. The extension of finite elements that were originally developed for isotropic one layered structures to multilayerd plates and shells are first discussed. Works based on classical and refined theories as well as on equivalent single layer and layer-wise descriptions have been overviewed. Development of available zig-zag finite elements has been considered for the three cases of zig-zag theories. Finite elements based on other approches are also discussed. Among these, FEs based on asymtotic theories, degenerate continuum approaches, stress resultant methods, asymtotic methods, hierarchy-p,_-s global/local techniques as well as mixed and hybrid formulations have been overviewed.

Theories and Finite Elements for Multilayered, Anisotropic, Composite Plates and Shells

This review article concerns a methodology for solving numerically, for engineering purposes, boundary and initial-boundary value problems by a peculiar approach characterized by the following features: the continuous formulation is centered on integral equations based on the combined use of single-layer and double-layer sources, so that the integral operator turns out to be symmetric with respect to a suitable bilinear form. The discretization is performed either on a variational basis or by a Galerkin weighted residual procedure, the interpolation and weight functions being chosen so that the variables in the approximate formulation are generalized variables in Prager’s sense. As main consequences of the above provisions, symmetry is exhibited by matrices with a key role in the algebraized versions; some quadratic forms have a clear energy meaning; variational properties characterize the solutions and other results, invalid in traditional boundary element methods enrich the theory underlying the computational applications. The present survey outlines recent theoretical and computational developments of the title methodology with particular reference to linear elasticity, elastoplasticity, fracture mechanics, time-dependent problems, variational approaches, singular integrals, approximation issues, sensitivity analysis, coupling of boundary and finite elements, and computer implementations. Areas and aspects which at present require further research are identified, and comparative assessments are attempted with respect to traditional boundary integral-elements. This article includes 176 references.

https://hal.archives-ouvertes.fr/hal-00111259/document

Symmetric Galerkin Boundary Element Methods

Vital signs such as pulse rate and breathing rate are currently measured using contact probes. But, non-contact methods for measuring vital signs are desirable both in hospital settings (e.g. in NICU) and for ubiquitous in-situ health tracking (e.g. on mobile phone and computers with webcams). Recently, camera-based non-contact vital sign monitoring have been shown to be feasible. However, camera-based vital sign monitoring is challenging for people with darker skin tone, under low lighting conditions, and/or during movement of an individual in front of the camera. In this paper, we propose distancePPG, a new camera-based vital sign estimation algorithm which addresses these challenges. DistancePPG proposes a new method of combining skin-color change signals from different tracked regions of the face using a weighted average, where the weights depend on the blood perfusion and incident light intensity in the region, to improve the signal-to-noise ratio (SNR) of camera-based estimate. One of our key contributions is a new automatic method for determining the weights based only on the video recording of the subject. The gains in SNR of camera-based PPG estimated using distancePPG translate into reduction of the error in vital sign estimation, and thus expand the scope of camera-based vital sign monitoring to potentially challenging scenarios. Further, a dataset will be released, comprising of synchronized video recordings of face and pulse oximeter based ground truth recordings from the earlobe for people with different skin tones, under different lighting conditions and for various motion scenarios.

DistancePPG: Robust non-contact vital signs monitoring using a camera

This paper deals with the formulation of finite plate elements for an accurate description of stress and strain fields in multilayered, thick plates subjected to static loadings in the linear, elastic cases. The so-called zig-zag form and interlaminar continuity are addressed in the considered formulations. Two variational statements, the principle of virtual displacements (PVD) and the Reissner mixed variational theorem (RMVT) are employed to derive finite element matrices. Transverse stress assumptions are made in the framework of RMVT and the resulting finite elements describe a priori interlaminar continuous transverse shear and normal stresses. Both modellings which preserve the number of variables independent of the number of layers (equivalent single-layer models, ESLM) and layer-wise models (LWM) in which the same variables are independent in each layer, have been treated. The order N of the expansions assumed for both displacement and transverse stress fields in the plate thickness direction z as well as the number of element nodes Nn have been taken as free parameters of the considered formulations. By varying N, Nn, variable treatment (LW or ESL) as well as variational statements (PVD and RMVT), a large number of newly finite elements have been presented. Finite elements that are based on PVD and RMVT have been called classical and advanced, respectively.

In order to write the matrices related to the considered plate elements in a concise form and to implement them in a computer code (see Part 2), extensive indicial notations have been set out. As a result, all the finite element matrices have been built from only five arrays that were called fundamental nuclei (four are related to RMVT applications and one to PVD cases). These arrays have 3×3 dimensions and are therefore constituted of only nine terms each. The different formulations are then obtained by expanding the indices that were introduced for the N-order expansion, for the number of nodes Nn and for the constitutive layers Nl. Compliances and/or stiffness are accumulated from layer to multilayered level according to the corresponding variable treatment (ESLM or LWM). The numerical evaluations and assessment for the presented plate elements have been provided in the companion paper (Part 2), where it has been concluded that it is convenient to refer to RMVT as a variational tool to formulate multilayered plate elements that are able to give a quasi-three-dimensional description of stress/strain fields in multilayered thick structures. Copyright © 2002 John Wiley & Sons, Ltd.

/pdf/classical-and-advanced-multilayered-plate-elements-based-3ql3nv6vh8.pdf

Classical and advanced multilayered plate elements based upon PVD and RMVT. Part 1: Derivation of finite element matrices

Piezoelectric sensors are increasingly being used in active structural health monitoring, due to their durability, light weight and low power consumption. In the present work damage detection and characterization methodologies based on Lamb waves have been evaluated for aircraft panels. The applicability of various proposed delay-and-sum algorithms on isotropic and composite stiffened panels have been investigated, both numerically and experimentally. A numerical model for ultrasonic wave propagation in composite laminates is proposed and compared to signals recorded from experiments. A modified delay-and-sum algorithm is then proposed for detecting impact damage in composite plates with and without a stiffener which is shown to capture and localize damage with only four transducers.

https://iopscience.iop.org/article/10.1088/0964-1726/23/7/075007/pdf

Assessment of delay-and-sum algorithms for damage detection in aluminium and composite plates

In this contribution we present an application of the lowest order Virtual Element Method (VEM) to the problem of material computational homogenization. Material homogenization allows retrieving material properties through suitable volume averaging procedures, starting from a detailed representation of the micro-constituents of the considered material. The representation of such microstructure constitutes a remarkable effort in terms of data/mesh preparation, especially when there is not evident microstructural regularity. For such a reason, computational micromechanics may represent a challenging benchmark for showing the potential of VEM. In this contribution, polycrystalline materials are considered as an application. The proposed technique constitutes a first step towards modelling of damage processes in micro-structured materials

Virtual Element Method: Micro-Mechanics Applications

In this work a three-dimensional BEM model is used for the analysis of structures with cracks and surface bonded piezoelectric PZT patches used as strain sensors. The cracked structure is modelled by the dual boundary element method, which allows for ac- curate and reliable crack analysis, while the piezoelec- tric patch is analyzed by a finite element state-space ap- proach, that embodies both the full electro-mechanical coupling and the suitable sensor's boundary conditions. The model is used to investigate the strain-transfer mech- anism from an host elastic structure to the piezoelectric layer, taking into account the effect of the adhesive layer, as well as the mechanical interaction between the sen- sor and the measured field. Finally, a fast solver based on the use of the GMRES and the hierarchical matri- ces is employed to speed up the solution of some large scale BEM systems arising from structural health mon- itoring applications. The numerical tool has been used to investigate the sensitivity of sensors arrangements in three-dimensional damaged solids. keyword: Piezoelectric patches, Structural health monitoring, Fast dual boundary element method structural damage by means of the so called damage sig- natures, consisting in adverse changes in some structural features, so to replace scheduled inspection cycles with a continuous monitoring system. In practice, the task of health monitoring is accomplished through a network of sensors and storage/computing units able to track the variables related to the structural response, such as strain, vibration, electrical conductivity and acoustic emission, which are affected by the presence of damage. From this point of view, SHM systems can be collocated in the broad range of the smart structures and their design and implementation rely upon the efficient integration of the sensor network with the structure and on the combination with appropriate electronics, modeling and control algo- rithms. Typical employable sensors for such tasks are strain gauges, accelerometers, fiber optics, piezoelectric films and piezoceramics (Van der Auweraer and Peeters (2003)).

Structures with Surface-Bonded PZT Piezoelectric Patches: a BEM Investigation into the Strain-transfer Mechanism for SHM applications

In this work, the grain-boundary cavitation in polycrystalline aggregates is investigated by means of a grain-scale model. Polycrystalline aggregates are generated using Voronoi tessellations, which have been extensively shown to retain the statistical features of real microstructures. Nucleation, thickening and sliding of cavities at grain boundaries are represented by specific cohesive laws embodying the damage parameters, whose time evolution equations are coupled to the mechanical model. The formulation is presented within the framework of a grain-boundary formulation, which only requires the discretization of the grain surfaces. Some numerical tests are presented to demonstrate the feasibility of the method.

A Micro-Mechanical Model for Grain-Boundary Cavitation in Polycrystalline Materials

Advances in Fracture and Damage Mechanics XII

In this work, a unified scheme for computing the fundamental solutions of a three-dimensional homogeneous elliptic partial differential operator is presented. The scheme is based on the Rayleigh expansion and on the Fourier representation of a homogeneous function. The scheme has the advantage of expressing the fundamental solutions and their derivatives up to the desired order without any term-by-term differentiation. Moreover, the coefficients of the series need to be computed only once, thus making the presented scheme attractive for numerical implementation. The scheme is employed to compute the fundamental solution of isotropic elasticity showing that the spherical harmonics expansions provide the exact expressions. Then, the accuracy of the scheme is assessed by computing the fundamental solutions of a generally anisotropic magneto-electro-elastic material.

Alberto Milazzo

Papers

Virtual Element Method: Micro-Mechanics Applications

Structures with Surface-Bonded PZT Piezoelectric Patches: a BEM Investigation into the Strain-transfer Mechanism for SHM applications

A Micro-Mechanical Model for Grain-Boundary Cavitation in Polycrystalline Materials

Advances in Fracture and Damage Mechanics XII

Spherical Harmonics Expansion of Fundamental Solutions and Their Derivatives for Homogeneous Elliptic Operators