Search results

The purpose of this paper is to present eight local elasto‐plastic beam element formulations incorporated into the corotational framework for two‐noded three‐dimensional beams. These formulations capture the warping torsional effects of open cross‐sections and are suitable for the analysis of the nonlinear buckling and post‐buckling of thin‐walled frames with generic cross‐sections. The paper highlights the similarities and discrepancies between the different local element formulations. The primary goal of this study is to compare all the local element formulations in terms of accuracy, efficiency and CPU‐running time.

The definition of the corotational framework for a two‐noded three‐dimensional beam element is presented, based upon the works of Battini .The definitions of the local element kinematics and displacements shape functions are developed based on both Timoshenko and Bernoulli assumptions, and considering low‐order as well as higher‐order terms in the second‐order approximation of the Green‐Lagrange strains. Element forces interpolations and generalized stress resultant vectors are then presented for both mixed‐based Timoshenko and Bernoulli formulations. Subsequently, the local internal force vector and tangent stiffness matrix are derived using the principle of virtual work for displacement‐based elements and the two‐field Hellinger‐Reissner assumed stress variational principle for mixed‐based formulations, respectively. A full comparison and assessment of the different local element models are performed by means of several numerical examples.

In this study, it is shown that the higher order elements are more accurate than the low‐order ones, and that the use of the higher order mixed‐based Bernoulli element seems to require the least number of FEs to accurately model the structural behavior, and therefore allows some reduction of the CPU time compared to the other converged solutions; where a larger number of elements are needed to efficiently discretize the structure.

The paper reports computation times for each model in order to assess their relative efficiency. The effect of the numbers of Gauss points along the element length and within the cross‐section are also investigated.

The purpose of this paper is to propose an efficient finite element formulation for nonlinear analysis of clustered tensegrity that consists of classical cables, clustered cables and bars.

The derivation of the finite element formulation is based on the co-rotational approach, which decomposes a geometrically nonlinear deformation into a large rigid body motion and a small-strain deformation. A tangent stiffness matrix of a clustered cable is proposed and the Newton-Raphson scheme is employed to solve the nonlinear equation.

The derived tangent stiffness matrix, including an additional stiffness terms that describes the slide effect of pulleys, can regress to the stiffness matrix of a classical cable, which is convenient for the implementation of finite element procedure. Two typical numerical examples show that the proposed formulation is accurate and requires less iteration than the force density method.

The co-rotational formulation of a clustered cable is originally proposed, although some mature methods, such as the TL, Force Density and Dynamic Relaxation method, have been applied to nonlinear analysis of clustered tensegrity. The proposed co-rotational formulation proved efficient.

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 aim of this paper is to develop a new method for finding multiple bifurcation points in structures.

A brief review of nonlinear analysis is presented. A powerful method (called arc‐length) for tracing nonlinear equilibrium path is described. Techniques for monitoring critical points are discussed to find the rank deficiency of the stiffness matrix. Finally, by using eigenvalue perturbation of tangent stiffness matrix, load parameter associated with multiple bifurcation points is obtained.

Since other methods of finding simple bifurcation points diverge in the neighborhood of critical points, this paper introduces a new method to find multiple bifurcation points. It should be remembered that a simple bifurcation point is a multiple bifurcation point with rank deficiency equal to one. Therefore, the method is applicable to simple critical points as well.

Global buckling of the structures should be considered in design. Many structures (specially symmetric space structures) have multiple bifurcation points, therefore, analyst and designer should be aware of these points and should control them (for example, by changing the geometry or other related factors) for obtaining a safe and optimum design.

In this paper a robust method to find multiple bifurcation points is introduced. By using this method, engineers can be aware of critical load of multiple bifurcation points to control global buckling of related structures.

The general problem of sizing, material and loading parameter sensitivity of non‐linear systems is presented. Both kinematic and path‐dependent material non‐linearities are considered; non‐linear sensitivity path is traced by an incremental solution strategy. The variational approach employed is quite general and can be employed for studying sensitivity of various path‐dependent highly non‐linear phenomena. Both the direct differentiation method (DDM) and adjoint system method (ASM) are discussed in the context of continuum and finite element mechanics. The merits of using the consistent tangent matrix and the necessity of accumulation of design derivatives of stresses and internal parameters are indicated. Aspects of sensitivity problems in metal forming are also discussed. A number of examples illustrate the paper.

The purpose of this paper is to describe several novel techniques for implementing a crystal plasticity (CP) material model in a large deformation, implicit finite element framework.

Starting from the key kinematic assumptions of CP, the presentation develops the necessary CP correction terms to several common objective stress rates and the consistent linearization of the stress update algorithm. Connections to models for slip system hardening are isolated from these processes.

A kinematically consistent implementation is found to require a correction to the stress update to include plastic vorticity developed by slip deformation in polycrystals. A simpler, more direct form for the algorithmic tangent is described. Several numerical examples demonstrate the capabilities and computational efficiency of the formulation.

The implementation assumes isotropic slip system hardening. With simple modifications, the described approach extends readily to anisotropic coupled or uncoupled hardening functions.

The modular formulation and implementation support streamlined development of new models for slip system hardening without modifications of the stress update and algorithmic tangent computations. This implementation is available in the open-source code WARP3D.

In the process of developing the CP formulation, this work realized the need for corrections to the Green-Naghdi and Jaumann objective stress rates to account properly for non-zero plastic vorticity. The paper describes fully the consistent linearization of the stress update algorithm and details a new scheme to implement the model with improved efficiency.

The purpose of this paper is to make a theoretical comprehensive efficiency evaluation of a nonlinear analysis method based on the Woodbury formula from the efficiency of the solution of linear equations in each incremental step and the selected iterative algorithms.

First, this study employs the time complexity theory to quantitatively compare the efficiency of the Woodbury formula and the LDLT factorization method which is a commonly used method to solve linear equations. Moreover, the performance of iterative algorithms also significantly effects the efficiency of the analysis. Thus, the three-point method with a convergence order of eight is employed to solve the equilibrium equations of the nonlinear analysis method based on the Woodbury formula, aiming to improve the iterative performance of the Newton–Raphson (N–R) method.

First, the result shows that the asymptotic time complexity of the Woodbury formula is much lower than that of the LDLT factorization method when the number of inelastic degrees of freedom (IDOFs) is much less than that of DOFs, indicating that the Woodbury formula is more efficient for local nonlinear problems. Moreover, the time complexity comparison of the N–R method and the three-point method indicates that the three-point method is more efficient than the N–R method for local nonlinear problems with large-scale structures or a larger ratio of IDOFs number to the DOFs number.

This study theoretically evaluates the efficiency of nonlinear analysis method based on the Woodbury formula, and quantitatively shows the application condition of the comparative methods. The comparison result provides a theoretical basis for the selection of algorithms for different nonlinear problems.

Some frictional contact problems are characterized by significant variations in the location and size of the contact area occurring in the process of deformation. When this feature is combined with strongly non‐linear, path‐dependent material behaviour, difficulties with convergence of the typically used iterative processes can be encountered. Demonstrates this by analysis of press‐fit connection, a typical problem in which both of those characteristics can be present. Offers an explanation as to the possible source of those difficulties. Suggests in support of this explanation, two simple modifications of the usual iterative schemes. In spite of their simplicity, they are found to be more robust than those usual schemes which are normally used in numerical analysis of similar problems.

Recent publications have highlighted the effectiveness of using a consistent tangent modulus when solving elastic‐plastic problems. The formulation of a consistent tangent modulus is closely related to the scheme used to integrate the constitutive equations. Recent work has shown how many of these schemes currently in use can be derived from certain broad classes of algorithms. In this paper these procedures are examined for a number of commonly used yield/failure criteria. For certain cases a remarkably simple formulation results which can lead to considerable savings in computational time.

Presents a robust and unconditionally stable return‐mapping algorithm based on the discrete counterpart of the principle of maximum plastic dissipation. Develops the explicit expression for the consistent elasto‐plastic tangent modulus. All expressions are derived via tensor formulation showing the advantage over the classical matrix notation. The integration algorithm is implemented in the formulation of the four‐node isoparametric assumed‐strain finite‐rotation shell element employing the Mindlin‐Reissner‐type shell model. By applying the layered model, plastic zones can be displayed through the shell thickness. Material non‐linearity described by the von Mises yield criterion and isotropic hardening is combined with a geometrically non‐linear response assuming finite rotations. Numerical examples illustrate the efficiency of the present formulation in conjunction with the standard Newton iteration approach, in which no line search procedures are required. Demonstrates the excellent performance of the algorithm for large time respective load steps.