Search results

The purpose of this paper is to evaluate the capabilities of the generalized finite element method (GFEM) under the context of the geometrically nonlinear analysis. The effect of large displacements and deformations, typical of such analysis, induces a significant distortion of the element mesh, penalizing the quality of the standard finite element method approximation. The main concern here is to identify how the enrichment strategy from GFEM, that usually makes this method less susceptible to the mesh distortion, may be used under the total and updated Lagrangian formulations.

An existing computational environment that allows linear and nonlinear analysis, has been used to implement the analysis with geometric nonlinearity by GFEM, using different polynomial enrichments.

The geometrically nonlinear analysis using total and updated Lagrangian formulations are considered in GFEM. Classical problems are numerically simulated and the accuracy and robustness of the GFEM are highlighted.

This study shows a novel study about GFEM analysis using a complete polynomial space to enrich the approximation of the geometrically nonlinear analysis adopting the total and updated Lagrangian formulations. This strategy guarantees the good precision of the analysis for higher level of mesh distortion in the case of the total Lagrangian formulation. On the other hand, in the updated Lagrangian approach, the need of updating the degrees of freedom during the incremental and iterative solution are for the first time identified and discussed here.

A simple beam element is developed for the solution of large deflection problems. The total Lagrangian formulation is based on the kinematic relations proposed by Reissner for finite rotations and stretching as well as shearing of plane beams. The motion is discretized by linear expansions of the global displacement components and the cross‐sectional rotation in two‐dimensional Euclidean space yielding a simple beam element with three degrees of freedom at the two nodes. The shear locking is reduced by selective integration in order to eliminate the spurious shear constraint similar to interdependent variable interpolation. The large rotation formulation is compared with two forms of moderate rotation theories which have been used in the past to develop the geometric stiffness properties for linear stability analysis of the so‐called Mindlin plate elements. The predictive value of different geometric stiffness approximations is assessed with several examples which range from the static and kinetic stability analysis of the classical Euler‐column to the large deflection problem of a clamped beam.

A generalized displacement method has been previously presented for the analysis of thin plate‐shell structures with the use of bilinear 4‐node isoparametric shell elements. Following this approach, a procedure for the geometrically non‐linear analysis of thin plates and shells based on both updated and total Lagrangian formulations is developed. The results of some numerical examples are presented to show the versatility and effectiveness of the method.

Large strain model can be formulated in terms of the Lagrangian or the Eulerian frame. In this paper, the Eulerian type large strain models are studied. Numerical examples on the Lagrangian and Eulerian types large strain models are investigated and compared. It is found that the differences in the choice of large strain model under large strain and rotation problems are noticeable but not significant if small load step is used for analysis. Furthermore, we have also found that unsymmetrical formulation instead of symmetrical formulation should be adopted for Eulerian type large strain models.

The purpose of this paper is to present an efficient smoothed particle hydrodynamics (SPH) method particularly adapted for the geometrically nonlinear analysis of structures.

In order to resolve the inconsistency phenomenon which systematically occurs in the standard SPH method at the domain’s boundaries of the studied structure, the classical kernel function and its spatial derivatives were modified by the use of Taylor series expansion. The well-known tensile instabilities inherent to the Eulerian SPH formulation were attenuated by the use of the Total Lagrangian Formulation (TLF).

In order to demonstrate the effectiveness of the present improved SPH method, several numerical applications involving geometrically nonlinear behaviors were carried out using the explicit dynamics scheme for the time integration of the PDEs. Comparisons of the obtained results using the present SPH model with analytical reference solutions and with those obtained using ABAQUS finite element (FE) commercial software, show its good accuracy and robustness.

An additional application including a multilayered composite structure and involving buckling and delamination was investigated using the present improved SPH model and the results are compared to the FE results, they confirmed both the efficiency and the accuracy of the proposed method.

An efficient 2D-continuum SPH model for the geometrically nonlinear analysis of thin and thick structures is proposed. Contrarily to the classical SPH approaches, here the constitutive material relations are used to link naturally the stresses and strains. The Total Lagrangian approach is investigated to alleviate the tensile instabilities problem, allowing at the same time to avoid the updating procedure of the neighboring particles search and therefore reducing CPU usage. The proposed approach is valid for isotropic and multilayered composites structures undergoing large transformations. CPU time savings and better results with the new 2D-continuum SPH formulation compared to the classical continuum SPH. The explicit dynamic scheme was used for time integration allowing a fast resolution algorithm even for highly nonlinear problems.

A previously presented finite element shell formulation is extended to the application of large strains. The finite elements are those based on the concept of ‘the degenerated solids’, which are widely used in non‐linear finite element programs. The constitutive equation of hyperelastic incompressible material is adopted and specialized to the Mooney‐Rivlin law. The additional state variable, the hydrostatic pressure, which occurs for incompressible materials, is eliminated on element level using the plane stress condition. Attention is drawn to the efficient calculation of the element matrices by applying a layer concept. The effectiveness of the proposed total Lagrangian formulation is demonstrated in a number of example problems.

This article presents a geometrically non-linear finite element formulation for the analysis of planar two-layer beam-columns taking into account the inter-layer slip and uplift.

The co-rotational method is adopted, in which the motion of the element is decomposed into a rigid body motion and a small deformational one. The geometrically linear formulation can be used in the local frame and automatically be transformed into a geometrically nonlinear one. In co-rotational frame, both layers are assumed to be discretely connected at the element ends. Slips and uplifts are assumed to be small. Consequently, the condition of non interpenetration between the layers can be treated using a node-to-node contact algorithm. The resolution methods such as penalty (PM) and augmented Lagrangian method (ALM) with Uzawa updating scheme can be used.

The non-penetration condition between the layers of composite beams can be formulated by using contact law. It is found that despite a low convergence rate of augmented Lagrangian method compared to penalty method, the former prevents the unrealistic penetration. Besides, it is shown that the buckling load of the composite beam-column is largely affected by the uplift stiffness of the connectors.

The proposed finite element model is capable of simulating accurately the geometrically non-linear behavior of planar two-layer beam-columns taking into account the inter-layer slip and uplift. Regarding uplift, the non-penetration condition is strictly enforced by considering rigorous contact conditions at the interface. The constraint problem is solved using the penalty method or the augmented Lagrangian method with the Uzawa updating scheme.

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.

Gives a bibliographical review of the finite element methods (FEMs) applied for the linear and nonlinear, static and dynamic analyses of basic structural elements from the theoretical as well as practical points of view. The range of applications of FEMs in this area is wide and cannot be presented in a single paper; therefore aims to give the reader an encyclopaedic view on the subject. The bibliography at the end of the paper contains 2,025 references to papers, conference proceedings and theses/dissertations dealing with the analysis of beams, columns, rods, bars, cables, discs, blades, shafts, membranes, plates and shells that were published in 1992‐1995.

nonlinear dynamic analysis of triangular and quadrilateral membrane elements with in-plane drilling rotational degree of freedom.

The nonlinear analysis is carried out using the updated co-rotational Lagrangian description. In this purpose, in-plane co-rotational formulation that considers the in-plane drilling rotation is developed and presented for triangular and quadrilateral elements, and a tangent stiffness matrix is derived. Furthermore, a simple and effective in-plane mass matrix that takes into account the in-plane rotational inertia, which permit true representation of in-plane vibrational modes is adopted for dynamic analysis, which is carried out using the Newmark direct time integration method.

The proposed numerical tests show that the presented elements exhibit very good performances and could return true in-plane rotational vibrational modes. Also, when using a well-chosen co-rotational formulation these elements shows good results for nonlinear static and dynamic analysis.

Publications that describe geometrical nonlinearity of the in-plane behaviour of membrane element with rotational d.o.f are few, and often they are based on the total Lagrangian formulation or on the rate form. Also these elements, at the author knowledge, have not been extended to the nonlinear dynamic analysis. Thus, an appropriate extension of triangular and quadrilateral membrane elements with drilling rotation to nonlinear dynamic analysis is required.