Search results

This paper focuses on the development of a new class of eight‐node solid finite elements, suitable for the treatment of volumetric and transverse shear locking problems. Doing so, the proposed elements can be used efficiently for 3D and thin shell applications. The starting point of the work relies on the analysis of the subspace of incompressible deformations associated with the standard (displacement‐based) fully integrated and reduced integrated hexahedral elements. Prediction capabilities for both formulations are defined related to nearly‐incompressible problems and an enhanced strain approach is developed to improve the performance of the earlier formulation in this case. With the insight into volumetric locking gained and benefiting from a recently proposed enhanced transverse shear strain procedure for shell applications, a new element conjugating both the capabilities of efficient solid and shell formulations is obtained. Numerical results attest the robustness and efficiency of the proposed approach, when compared to solid and shell elements well‐established in the literature.

In order to avoid the problem of volumetric locking often encountered when using constant strain tetrahedral finite elements, the purpose of this paper is to present a new composite tetrahedron element which is especially designed for the combined finite-discrete element method (FDEM).

A ten-noded composite tetrahedral (COMPTet) finite element, composed of eight four-noded low order tetrahedrons, has been implemented based on Munjiza’s multiplicative decomposition approach. This approach naturally decomposes deformation into translation, rotation, plastic stretches, elastic stretches, volumetric stretches, shear stretches, etc. The problem of volumetric locking is avoided via a selective integration approach that allows for different constitutive components to be evaluated at different integration points.

A number of validation cases considering different loading and boundary conditions and different materials for the proposed element are presented. A practical application of the use of the COMPTet finite element is presented by quantitative comparison of numerical model results against simple theoretical estimates and results from acrylic fracturing experiments. All of these examples clearly show the capability of the composite element in eliminating volumetric locking.

For this tetrahedral element, the combination of “composite” and “low order sub-element” properties are good choices for FDEM dynamic fracture propagation simulations: in order to eliminate the volumetric locking, only the information from the sub-elements of the composite element are needed which is especially convenient for cases where re-meshing is necessary, and the low order sub-elements will enable robust contact interaction algorithms, which maintains both relatively high computational efficiency and accuracy.

This paper aims to present a 3D static performance analysis model for the gas foil journal bearing to provide better understanding of the gas foil journal bearing and extend the development of the calculation about the static performance.

The foil bearing can be seen as a shell structure, and the mixed interpolation of tensorial components (MITC) element was used to build the shell model. The augmented Lagrange method was used to calculate the contact involving friction between foils and between the foil and the bearing sleeve. A displacement-controlled load scheme was used to calculate the deformation of the foils. A mapping operator was used to map the film pressure from the gas to the surface of the top foil.

This method provides high precision of calculation in the prediction of the static performance. The calculation results were compared with the experimental data, and they show good agreement. Meanwhile, the model can be applied in the prediction of the bearing performance in a broad range of working conditions.

This method extends the calculation of the gas foil journal bearing to a 3D scale and shows good agreement with the experimental data. Meanwhile, the present model has a good adaptability on the revolution speed and can be applied to the predictions in varied working conditions.

As known from nearly incompressible elasticity, selective reduced integration (SRI) is a simple and effective method of overcoming the volumetric locking problem in 2D and 3D solid elements. This method of finite elastoviscoplasticity is discussed as are its well‐known limitations. In this context, an isochoric‐volumetric decoupled material behavior is assumed and thus the additive deviatoric‐volumetric decoupling of the consistent algorithmic moduli tensor is essential. By means of several numerical examples, the performance of elements using selective reduced integration is demonstrated and compared to the performance of other elements such as the enhanced assumed strain elements. It is shown that a minor modification, with little numerical effort, leads to rather robust element behaviour. The application of this process to so‐called solid‐shell elements for thin‐walled structures is also discussed.

The purpose of this paper is to investigate the application of the element‐free Galerkin (EFG) method to the simulation of metal forming processes and to propose a strategy to deal with volumetric locking problem in this context.

The J₂ elastoplastic material model, employed in the work, assumes a multiplicative decomposition of the deformation gradient into an elastic and a plastic part and incorporates a non‐linear isotropic hardening response. The constitutive model is written in terms of the rotated Kirchhoff stress and the logarithmic strain measure. A Total Lagrangian formulation of the problem is considered in order to improve the computational performance of the proposed algorithm. The imposition of the essential boundary conditions and also of the unilateral contact with friction condition are made by the application of the Augmented Lagrangian method. Here, aspects related to the volumetric locking are investigated and an F‐bar approach is applied.

The results show that the proposed approach presents no volumetric locking phenomenon when using the mean dilation approach. Moreover, differently from finite element approximations, no hour‐glass instabilities in the deformation pattern are observed, avoiding in this way the need to devise additional stabilization procedures in the proposed procedure.

This paper demonstrates the implementation and validation of the mean dilation approach, in the scope of the EFG, which was successful in coping with the volumetric locking phenomena and presented no hour‐glass instabilities in the problem cases considered in this work.

An important characteristic of many soil models is a volume change during plastic flow. In computations, this plastic volume change is expressed via a kinematic constraint on the possible deformations. Due to this constraint the plane‐strain three‐noded triangular element exhibits locking when plastic deformations occur, under dilatant, contractant and isochoric conditions. It is demonstrated that using the method of enhanced assumed strains by Simol this locking cannot be remedied. For six‐noded wedges and four‐noded and five‐noded pyramids the same conclusion is obtained.

Considers the problem of stability of the enhanced strain elements in the presence of large deformations. The standard orthogonality condition between the enhanced strains and constant stresses ensures satisfaction of the patch test and convergence of the method in case of linear elasticity. However, this does not hold in the case of large deformations. By analytic derivation of the element eigenvalues in large strain states additional orthogonality conditions can be derived, leading to a stable formulation, regardless of the magnitude of deformations. Proposes a new element based on a consistent formulation of the enhanced gradient with respect to new orthogonality conditions which it retains with four enhanced modes volumetric and shear locking free behaviour of the original formulation and does not exhibit hour‐glassing for large deformations.

The purpose of this paper is to present a new stabilised low-order finite element methodology for large strain fast dynamics.

The numerical technique describing the motion is formulated upon the mixed set of first-order hyperbolic conservation laws already presented by Lee et al. (2013) where the main variables are the linear momentum, the deformation gradient tensor and the total energy. The mixed formulation is discretised using the standard explicit two-step Taylor-Galerkin (2TG) approach, which has been successfully employed in computational fluid dynamics (CFD). Unfortunately, the results display non-physical spurious (or hourglassing) modes, leading to the breakdown of the numerical scheme. For this reason, the 2TG methodology is further improved by means of two ingredients, namely a curl-free projection of the deformation gradient tensor and the inclusion of an additional stiffness stabilisation term.

A series of numerical examples are carried out drawing key comparisons between the proposed formulation and some other recently published numerical techniques.

Both velocities (or displacements) and stresses display the same rate of convergence, which proves ideal in the case of industrial applications where low-order discretisations tend to be preferred. The enhancements introduced in this paper enable the use of linear triangular (or bilinear quadrilateral) elements in two dimensional nearly incompressible dynamics applications without locking difficulties. In addition, an artificial viscosity term has been added into the formulation to eliminate the appearance of spurious oscillations in the vicinity of sharp spatial gradients induced by shocks.

Currently, there are some finite element head models developed by research groups all around the world. Nevertheless, the majority are not geometrically accurate. One of the problems is the brain geometry, which usually resembles a sphere. This may raise problems when reconstructing any event that involves brain kinematics, such as accidents, affecting the correct evaluation of resulting injuries. Thus, the purpose of this study is to develop a new finite element head model more accurate than the existing ones.

In this work, a new and geometrically detailed finite element brain model is proposed. Special attention was given to sulci and gyri modelling, making this model more geometrically accurate than currently available ones. In addition, these brain features are important to predict specific injuries such as brain contusions, which usually involve the crowns of gyri.

The model was validated against experimental data from impact tests on cadavers, comparing the intracranial pressure at frontal, parietal, occipital and posterior fossa regions.

As this model is validated, it can be now used in accident reconstruction and injury evaluation and even as a design tool for protective head gear.

An 8‐node solid element applicable for thin structures is presented. The element employs eighteen assumed stress modes and the conventional displacement interpolation. The formulation starts with the hybrid stress element proposed by Pian and Tong. The higher order stress modes are first decomposed into the ones which do and do not lead to thin‐element locking. The recently established methodology of admissible matrix formulation allows the decoupling of the above two categories of stress modes in the flexibility matrix without triggering element instability or failure of the patch test. The element stiffness can thus be decomposed into a series of matrices. Locking can be eliminated by adjusting the magnitude of the pertinent matrices. Accuracy and convergence rate of the present element are found to be competent to many of the existing plate and shell models.