Search results

The present paper is directed towards elasto‐plastic large deformation analysis of thin shells based on the concept of degenerated solids. The main aspect of the paper is the derivation of an efficient computational strategy placing emphasis on consistent elasto‐plastic tangent moduli and stress integration with the radial return method under the restriction of ‘zero normal stress condition’ in thickness direction. The advantageous performance of the standard Newton iteration using a consistent tangent stiffness matrix is compared to the classical scheme with an iteration matrix based on the infinitesimal elasto‐plastic constitutive tensor. Several numerical examples also demonstrate the effectiveness of the standard Newton iteration with respect to modified and quasi‐Newton methods like BFGS and others.

The present study focusses on algorithmic aspects related to deformation dependent loads in non‐linear static finite element analysis. If the deformation dependency is considered only on the right hand side, a considerable increase in the number of iterations follows. It may also cause failure of convergence in the proximity of critical points. If in turn the deformation dependent loading is included within the consistent linearization, an additional left hand side term emerges, the so‐called load stiffness matrix. In this paper several numerical test cases are used to show and quantify the influence of the two different approaches on the iteration process. Consideration of the complete load stiffness matrix may result in a cumbersome coding effort, different for each load case, and in certain cases its derivation is even not practicable at all. Therefore also several formulations for approximated load stiffness matrices are presented. It is shown that these simplifications not only reduce the additional effort for linearization and implementation, but also keep the iterative costs relatively small and still allow the calculation of the entire equilibrium path.

In this work we present interrelations between different finite rotation parametrizations for geometrically exact classical shell models (i.e. models without drilling rotation). In these kind of models the finite rotations are unrestricted in size but constrained in the 3‐d space. In the finite element approximation we use interpolation that restricts the treatment of rotations to the finite element nodes. Mutual relationships between different parametrizations are very clearly established and presented by informative commutative diagrams. The pluses and minuses of different parametrizations are discussed and the finite rotation terms arising in the linearization are given in their explicit forms.

In this paper, a new homogenization technique for the determination of dynamic and kinematic quantities of representative elementary volumes (REVs) in granular assemblies is presented. Based on the definition of volume averages, expressions for macroscopic stress, couple stress, strain and curvature tensors are derived for an arbitrary REV. Discrete element model simulations of two different test set‐ups including cohesionless and cohesive granular assemblies are used as a validation of the proposed homogenization technique. A non‐symmetric macroscopic stress tensor, as well as couple stresses are obtained following the proposed procedure, even if a single particle is described as a standard continuum on the microscopic scale.

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.

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.

The present contribution deals with the sensitivity analysis and optimization of structures for path‐dependent structural response. Geometrically as well as materially non‐linear behavior with hardening and softening is taken into account. Prandtl‐Reuss‐plasticity is adopted so that not only the state variables but also their sensitivities are path‐dependent. Because of this the variational direct approach is preferred for the sensitivity analysis. For accuracy reasons the sensitivity analysis has to be consistent with the analysis method evaluating the structural response. The proposed sensitivity analysis as well as its application in structural optimization is demonstrated by several examples.

The use of enhanced strains leads to an improved performance of low order finite elements. A modified Hu‐Washizu variational formulation with orthogonal stress and strain functions is considered. The use of orthogonal functions leads to a formulation with B (overline) ‐strain matrices which avoids numerical inversion of matrices. Depending on the choice of the stress and strain functions in Cartesian or natural element coordinates one can recover, for example, the hybrid stress element P‐S of Pian‐Sumihara or the Trefftz‐type element QE2 of Piltner and Taylor. With the mixed formulation discussed in this paper a simple extension of the high precision elements P‐S and QE2 to general non‐linear problems is possible, since the final computer implementation of the mixed element is very similar to the implementation of a displacement element. Instead of sparse B‐matrices, sparse B (overline) ‐matrices are used and the typical matrix inversions of hybrid and mixed methods can be avoided. The two most efficient four‐node B (overline) ‐elements for plane strain and plane stress in this study are denoted B (overline)(x, y)‐QE4 and B (overline)(ξ, η)‐QE4.

This study aims to develop a new analysis approach devised by incorporating a gradient-enhanced microplane damage model (GeMpDM) into isogeometric analysis (IGA), which shows computational stability and capability in accurately predicting crack propagations in structures with complex geometries.

For the non-local microplane damage modeling, the maximum modified von-Mises equivalent strain among all microplanes is regularized as a representative quantity. This characterization implies that only one additional governing equation is considered, which improves computational efficiency dramatically. By combined use of GeMpDM and IGA, quasi-static and dynamic numerical analyses are conducted to demonstrate the capability in predicting crack paths of complex geometries in comparison to FEM and experimental results.

The implicit scheme with the adopted damage model shows favorable numerical stability and the numerical results exhibit appropriate convergence characteristics concerning the mesh size. The damage evolution is successfully controlled by a tension-compression damage factor. Thanks to the advanced geometric design capability of IGA, the details of crack patterns can be predicted reliably, which are somewhat difficult to be acquired by FEM. Additionally, the damage distribution obtained in the dynamic analysis is in close agreement with experimental results.

The paper originally incorporates GeMpDM into IGA. Especially, only one non-local variable is considered besides the displacement field, which improves the computational efficiency and favorable convergence characteristics within the IGA framework. Also, enjoying the geometric design ability of IGA, the proposed analysis method is capable of accurately predicting crack paths reflecting the complex geometries of target structures.

In this paper we derive a simple finite element formulation for geometrical nonlinear shell structures. The formulation bases on a direct introduction of the isoparametric finite element formulation into the shell equations. The element allows the occurrence of finite rotations which are described by two independent angles. A layerwise linear elastic material model for composites has been chosen. A consistent linearization of all equations has been derived for the application of a pure Newton method in the nonlinear solution process. Thus a quadratic convergence behaviour can be achieved in the vicinity of the solution point. Examples show the applicability and effectivity of the developed element.