Three-dimensional thermal buckling analysis of functionally graded material structures using a modified FSDT-based solid-shell element

On the comparison of two solid-shell formulations based on in-plane reduced and full integration schemes in linear and non-linear applications

On the use of a reduced enhanced solid-shell (RESS) element for sheet forming simulations

The main objective of this paper is to develop a numerical model susceptible to solve the numerical locking problems that may appear when applying the conventional solid and shell finite elements of ABAQUS. This model is based on a hexahedral solid shell element. The formulation of this element relays on the combination of the enhanced assumed strain (EAS) and assumed natural strain (ANS) methods with modified First Shear Deformation Theory (FSDT). The developed element is implemented into the ABAQUS user element (UEL) interface. The performance of this element is demonstrated by different benchmark tests from the literature. Our contribution consists on applying a single solid shell element through the thickness direction to predict the low velocity impact behavior on functionally graded material (FGM) circular plates.

Hygro-thermo-mechanical modeling of thin-walled photovoltaic laminates with polymeric interfaces

This study is focusing on the finite element simulation of low velocity impact applied to glass fiber-reinforced composites. Elasto-visco-plasticity is adopted taking into account strain rate effect, with isotropic hardening. To avoid numerical locking, hexahedral solid shell finite element SS is developed. The formulation is based on the mixed Assumed Natural Strain (ANS) and the Enhanced Assumed Strain (EAS) methods. A geometrical non linear condition is applied. The finite element is implemented in UEL subroutine via the commercial software ABAQUS/Standard. The numerical element proved its efficiency comparing to conventional element and previous experimental results in (Mars et al. 2018). Numerical simulations with solid element are conducted using Abaqus/standard via user material UMAT subroutine. In terms of computation time, the solid shell element is less time consuming when comparing to conventional solid element. This is due to the fact that a single element through the thickness direction is adopted. Element efficiency is assessed for various fibers contents.

In this paper, we present an optimal low-order accurate piezoelectric solid-shell element formulation to model active composite shell structures that can undergo large deformation and large overall motion. This element has only displacement and electric degrees of freedom (dofs), with no rotational dofs, and an optimal number of enhancing assumed strain (EAS) parameters to pass the patch tests (both membrane and out-of-plane bending). The combination of the present optimal piezoelectric solid-shell element and the optimal solid-shell element previously developed allows for efficient and accurate analyses of large deformable composite multilayer shell structures with piezoelectric layers. To make the 3-D analysis of active composite shells containing discrete piezoelectric sensors and actuators even more efficient, the composite solid-shell element is further developed here. Based on the mixed Fraeijs de Veubeke–Hu–Washizu (FHW) variational principle, the in-plane and out-of-plane bending behaviours are improved via a new and efficient enhancement of the strain tensor. Shear-locking and curvature thickness locking are resolved effectively by using the assumed natural strain (ANS) method. We also present an optimal-control design for vibration suppression of a large deformable structure based on the general finite element approach. The linear-quadratic regulator control scheme with output feedback is used as a control law on the basis of the state space model of the system. Numerical examples involving static analyses and dynamic analyses of active shell structures having a large range of element aspect ratios are presented. Active vibration control of a composite multilayer shell with distributed piezoelectric sensors and actuators is performed to test the present element and the control design procedure. Copyright © 2005 John Wiley & Sons, Ltd.

The phase field (PF) approach to fracture has emerged as a promising modeling tool that regularizes the variational fracture theory by Griffith via the introduction of a nonlocal damage-like field variable in the corresponding formulation. In this work, we outline a PF formulation for triggering brittle fracture phenomena in shell structures made of Functionally Graded Materials (FGMs). This model relies on the 6-parameter shell formulation complying with a solid shell kinematic description and incorporates the use of the Enhanced Assumed Strain (EAS) and Assumed Natural Strain (ANS) methods in order to alleviate different locking pathologies. The corresponding multi-field variational formalisms is consistently derived and discretized within the context of the Finite Element Method (FEM). Details regarding the implementation in the general purpose FE packages are outlined. The applicability of this model is demonstrated by means of several numerical applications.

The parametrization of shell structures using the so-called solid shell concept has been widely exploited in the last decades. This trend is mainly attributed to the relatively simple kinematic treatment of solid shells in the corresponding finite element formulation in conjunction with the use of unmodified three-dimensional material laws, among other aspects. In the present investigation, we provide a comprehensive finite element implementation of solid shells incorporating: (i) the use of Enhanced Assumed Strain (EAS) and the Assumed Natural Strain (ANS) methods to prevent locking issues, (ii) the phase-field approach for triggering fracture events, and (iii) some representative inelastic material models. The current modular implementation has been integrated into the FE package ABAQUS via the user-defined routine UEL. Several representative examples demonstrate the applicability of the present formulation.

The principle of robot based incremental sheet metal forming is based on flexible shaping by means of a freely programmable path-synchronous movement of two tools, which are operated by two industrial robots. The final shape is produced by the incremental infeed of the forming tool in depth direction and its movement along the geometry’s contour in lateral direction. The main problem during the forming process is the influence on the dimensional accuracy resulting from the compliance of the involved machine structures and the springback effects of the workpiece. The project aims to predict these deviations caused by resiliences and to carry out a compensative path planning based on this prediction. Therefore a planning tool is implemented which compensates the robot’s compliance and the springback effects of the sheet metal. Finite element analysis using a material model developed at the Institute of Applied Mechanics (IFAM) [1] has been used for the simulation of the forming process. The finite strain constitutive model combines nonlinear kinematic and isotropic hardening and is derived in a thermodynamical setting. It is based on the multiplicative split of the deformation gradient in the context of hyperelasticity. The kinematic hardening component represents a continuum extension of the classical rheological model of Armstrong–Frederick kinematic hardening which is widely adopted as capable of representing the above metal hardening effects. The major problem of low-order finite elements used to simulate thin sheet structures, such as used for the experiments, is locking, a non-physical stiffening effect. Recent research focuses on the large deformation version of a new eight-node solid-shell finite element based on reduced integration with hourglass stabilization. In the solid-shell formulation developed at IFAM ([2], [3]) the enhanced assumed strain (EAS) concept as well as the assumed natural strain (ANS) concept are implemented to circumvent locking. These tools are very important to obtain a good correlation between experiment and simulation.

ABSTRACTThermal buckling and vibration of functionally graded (FG) sinusoidal microbeams with temperature-dependent properties and three kinds of temperature distributions are investigated in this article. As one material length scale is introduced, the modified couple stress theory is capable of predicting the small-scale effects. Material properties of FG microbeams are calculated using the Mori–Tanaka method. Furthermore, temperature-dependent properties are taken into account to investigate the mechanical characteristics of FG microbeams in high–thermal-gradient environment. Motion equations and the associated boundary conditions are obtained simultaneously through variational principle. Then Navier procedure and the differential quadrature method incorporating an iterative procedure are used to solve the governing differential equations with temperature-dependent properties and general boundary conditions. Numerical examples are performed for demonstrating the influences of temperature distribution, beam thickness, material length scale, slenderness ratio, shear deformation, functionally graded index, boundary conditions, and temperature-dependent/independent properties on thermal buckling and free vibration behaviors of FG microbeams.

Author index of volume 70

Optimum distribution of materials for functionally graded rectangular plates considering thermal buckling

Nonlinear numerical dissipative elastodynamics of an optimal solid shell element

The paper deals with the validation of a recently proposed hexahedral solid-shell finite element in the buckling analysis of a laminated composite plate with delaminations. The object is to study the buckling behavior of structures with delaminations using the enhanced assumed strain (EAS) solid shell element with 5, 7 and 9 parameters. The EAS three-dimensional finite element formulation presented in this paper is free from shear locking and leads to accurate results for distorted element shapes. The developed FE model is used to study the effects of some parameters in the buckling load, such as the stacking sequences, delamination size, aspect ratio, width-to-thickness ratio. The feasibility of the proposed method is confirmed by numerical examples. Results show that using hexahedral solid-shell finite element in the buckling analysis is more efficient than using the enhanced solid finite element.

On the finite element analysis of functionally graded sandwich curved beams via a new refined higher order shear deformation theory