Search results

The purpose of this paper is to examine quadratic convergence of finite element analysis for hyperelastic material at finite strains via Abaqus‐UMAT as well as classification of the rates of convergence for iterative solutions in regular cases.

Different formulations for stiffness – Hessian form of the free energy functionals – are systematically given for getting the rate‐independent analytical tangent and the numerical tangent as well as rate‐dependent tangents using the objective Jaumann rate of Kirchoff stress tensor as used in Abaqus. The convergence rates for available element types in Abaqus are computed and compared for simple but significant nonlinear elastic problems, such as using the 8‐node linear brick (B‐bar) element – also with hybrid pressure formulation and with incompatible modes – further the 20‐node quadratic brick element with corresponding modifications as well as the 6‐node linear triangular prism element and 4‐node linear tetrahedral element with modifications.

By using the Jaumann rate of Kirchoff stress tensor for both, rate dependent and rate independent problems, quadratic or nearly quadratic convergence is achieved for most of the used elements using Abaqus‐UMAT interface. But in case of using rate independent analytical tangent for rate independent problems, even convergence at all is not assured for all elements and the considered problems.

First time the convergence properties of 3D finite elements available in Abaqus sre systematically treated for elastic material at finite strain via Abaqus‐UMAT.

Several constitutive models are available in the literature to describe the mechanical behaviour of cement stabilized soils. However, difficulties in implementing such models within commercial finite element programs have hindered their application to solve related boundary value problems. Therefore, the aim of this study is to implement a constitutive model, which has the capability to simulate cement stabilized soil behaviour, into the finite element program ABAQUS through the user material subroutine UMAT.

After a detailed review of existing constitutive models for cement stabilized soils, a model based on the elasto‐plastic theory and the extended critical state concept with an associated flow rule is selected for the finite element implementation. A semi‐implicit integration method (cutting plane algorithm) with a continuum elasto‐plastic modulus and path dependent stress prediction strategy has been used in the implementation. The performance of the new finite element formulation of the constitutive model is verified by simulating triaxial test data using the finite element program with the new implementation and predictions from constitutive equations as well as experimental data.

The paper provides the implementation procedure of the constitutive model into ABAQUS but this method is useful for the implementation of any other constitutive model into ABAQUS or any other finite element program. Simulated results for the volumetric deformation of cement stabilized soils show that the cement stabilized soils do not obey the associated flow rule at high confining pressures. The parametric study shows that the influence of cementation increases the brittle nature and the bearing capacity of treated clay. In addition the results show that proposed finite element implementation has the ability to illustrate key features of the cement stabilized clay.

This paper presents an implementation of an elasto‐plastic constitutive model, based on the extended critical state concept, for cement stabilized soils into a finite element programme, which has been identified as an important and challenging topic in computational geomechanics. This implementation is useful in solving boundary value problems in geomechanics involving cement stabilized soils, incorporating key characteristics of these soils.

This paper aims to investigate the sealing characteristic of the rubber X-ring combined seal used in high-pressure hydrogen service and clarify the effect of swelling due to dissolved hydrogen on the sealing behavior.

A finite element analysis method with a user subroutine is proposed to investigate the sealing characteristic of rubber X-ring seals in conjunction with the swelling effect.

The swelling has a noticeable impact on the sealing behavior. The higher peak contact stress suggests that the X-ring seal may be superior to the O-ring seal, while the calculated increase in Mises stress suggests the X-ring exhibits a higher propensity to mechanical damage under low pressure conditions.

This paper describes a method with a user subroutine developed within ABAQUS to simulate the sealing performance coupled hydrogen swelling. The suitability of X-ring seals is evaluated. This work provides insight into the swelling effect on sealing behavior of an X-ring, which is of great value for the design and application of rubber seals used in hydrogen service.

The purpose of this paper is to suggest an implicit integration method for updating the constitutive relationships in the newly proposed anisotropic egg-shaped elastoplastic (AESE) model and to apply it in ABAQUS.

The implicit integration algorithm based on the Newton–Raphson method and the closest point projection scheme containing an elastic predictor and plastic corrector are implemented in the AESE model. Then, the integration code for this model is incorporated into the commercial finite element software ABAQUS through the user material subroutine (UMAT) interface to simulate undrained monotonic triaxial tests for various saturated soft clays under different consolidation conditions.

The comparison between the simulated results from ABAQUS and the experimental results demonstrates the satisfactory performance of this implicit integration algorithm in terms of effectiveness and robustness and the ability of the proposed model to predict the characteristics of soft clay.

The rotational hardening rule in the AESE model together with the implicit integration algorithm cannot be considered.

The singularity problem existing in most elastoplastic models is eliminated by the closed, smooth and flexible anisotropic egg-shaped yield surface form in the AESE model. In addition, this notion leads to an efficient implicit integration algorithm for updating the highly nonlinear constitutive equations for unsaturated soft clay.

This study aims to explain the characterization of cyclic behavior of a tube made of functionally graded material (FGM) under different combinations of internal pressure and cyclic through-thickness temperature gradients.

The normality rule, nonlinear kinematic hardening Chaboche model and Von Mises yield criterion were used to model the constitutive behavior of an FG tube in the incremental form. The material properties and hardening parameters of the Chaboche model vary according to the power-law function in the radial direction. The backward Euler integration scheme combined with return mapping algorithm which relies on the solution of a nonlinear equation performs the numerical procedure. The algorithm is implemented within the user subroutine UMAT in ABAQUS/standard.

The published works on FG components considering only the mechanical and physical properties as a function of spatial coordinate and nonlinear kinematic hardening parameters have not been considered to be changed continuously from one surface to another. Motivated by this, the present paper has deliberately been targeted to tackle this kind of problem to simulate the cyclic behavior of an FG tube as accurately as possible. In addition, to classify various behaviors the FG tube under cyclic thermomechanical loadings, Bree’s interaction diagram as an essential tool in designing of the FG pressure vessels in many engineering sectors is presented.

Provides a detailed description of the FG parameters of Chaboche kinematic hardening parameters in the adopted constitutive equations. In this paper, the significant effects of internal pressure values, kinematic hardening models and also FG inhomogeneity index related to the hardening rule parameters on plastic deformation of the FG tube are illustrated. Finally, the various cyclic behaviors of the FG tube under different combinations of thermomechanical loading are fully explored.

The purpose of this paper is to develop a progressive damage framework to predict the fatigue life of cord-reinforced rubber composite under cyclic loadings. Special attention has been paid to failure mechanisms, like cord–rubber interfacial debonding, and rubber matrix damage.

The constitutive modeling is based on the continuum damage mechanics (CDMs) and the thermodynamics of irreversible process. The damage in rubber is described by an istropic law, whereas elasto-plastic continuum model has been proposed for cord–rubber interphase layer. The numerical framework is implemented into commercial finite element code Abaqus/Standard via user subroutine (UMAT).

One of the most important findings obtained from reviewing various techniques is that meso-level fatigue damage modeling based on developed framework can simulate competitive damage scenarios, e.g. debonding, delamination or matrix failure.

A systematic framework for predicting failure in cord-reinforced rubber composite is formulated within the context of CDMs that can also be applied for industrial components, such as tires and airsprings.

This study aims to develop an innovative actuator for improving the performance of future aircraft, by adapting the airfoil shape according to the flight conditions. The flap’s camber of a civil regional transportation aircraft’s trailing edge actuated and morphed with the use of shape memory alloys (SMA) actuator technology, instead of the conventional split flap mechanism is studied.

For the flap’s members sizing an efficient methodology is utilised based on finite element (FE) stress analysis combined to analytically formulated design criteria. A mechanical simulation within an FE approach simulated the performance of the moving rib, integrating both aerodynamic loads and SMA phenomenology, implementing Lagouda’s constitutive model. Aim of this numerical simulation is to provide guidelines for further development of the flap. A three-dimensional assembly of the flap is constructed to produce manufacturing drawing and to ensure that during its morphing no interference between the members occurrs. Eventually, the manufactured flap is integrated on a test rig and the experimental characterisations under no and static loads, and dynamic excitation are performed.

Experimental results showed that the rib’s SMA mechanism can adequate function under load providing satisfactory morphing capabilities.

The investigated approach is an internal into the flap mechanism based on the shape memory effect of thin wires. In the developed mechanism, SMA wires are attached to the wing structure, where they function as actuating elements.

The authors examined the numerical natural frequency analysis of a 2D functionally graded (FG) truncated thick hollow cone using 3D elasticity theory.

The material properties of the 2D-FGM (two dimensional-functionally graded materials) cone are graded along the radial and axial axes of the cone using a power–law distribution. The eigenvalue problem was solved using finite element analysis (FEA) employing graded hexahedral elements, and the verification of the finite element approach was assessed by comparing the current solution to earlier experimental studies.

The effects of semivertex angle, material distribution and the cone configuration on the natural frequencies have been analyzed. For various semivertex angles, thickness, length and power law exponents, many results in the form of natural frequencies and mode shapes are presented for the 2D-FGM cone. As a result, the effects of the given parameters were addressed, and the results were compared, demonstrating the direct efficiency of raising the power–law exponents and cone thickness on the rise of natural frequencies.

For the first time, the numerical natural frequency analysis of a 2D-FG truncated thick hollow truncated cone based on 3D equations of elasticity has been investigated. The material properties of the truncated cone have been distributed along two directions, which has not been considered before in any research for the truncated thick cone. The reason for using these innovative volume fraction functions is the lack of accurate coverage by functions that are available in the literature (Asemi et al., 2011; Babaei et al. 2021).

The purpose of this study, which is designed for the implementation of models in the implicit finite element framework, is to propose a robust, stable and efficient explicit integration algorithm for rate-independent elasto-plastic constitutive models.

The proposed automatic substepping algorithm is founded on an explicit integration scheme. The estimation of the maximal subincrement size is based on the stability analysis.

In contrast to other explicit substepping schemes, the algorithm is self-correcting by definition and generates no cumulative drift. Although the integration proceeds with maximal possible subincrements, high level of accuracy is attained. Algorithmic tangent stiffness is calculated in explicit form and optionally no analytical second-order derivatives are needed.

The algorithm is convenient for elasto-plastic constitutive models, described with an algebraic constraint and a set of differential equations. This covers a large family of materials in the field of metal plasticity, damage mechanics, etc. However, it cannot be directly used for a general material model, because the presented algorithm is convenient for solving a set of equations of a particular type.

The estimation of the maximal stable subincrement size is computationally cheap. All expressions in the algorithm are in explicit form, thus the implementation is simple and straightforward. The overall performance of the approach (i.e. accuracy, time consumption) is fully comparable with a default (built-in) ABAQUS/Standard algorithm.

The estimated maximal subincrement size enables the algorithm to be stable by definition. Subincrements are much larger than those in conventional substepping algorithms. No error control, error correction or local iterations are required even in the case of large increments.

This paper aims to first apply more advanced anisotropic yield criterions as Yld91 and Yld2004 to spherical indentation simulations, and investigate plastic anisotropy identified from indentation simulations following different yield criterions (Hill48, Yld91, Yld2004) to discover laws. It also aims to compare the difference in plastic anisotropy identified from indentation on three yield criterions and evaluate the applicability of plastic anisotropy.

This paper uses indentation simulations on different yield criterions to identify plastic anisotropy. First, the trust-region techniques based on the nonlinear least-squares method are used to determine anisotropy coefficients of Yld91 and Yld2004. Then, Yld91 and Yld2004 are implemented into ABAQUS software using user-defined material (UMAT) subroutines with the proposed universal structure. Finally, through considering comprehensively the key factors, the locations of the optimal data acquisition points in indentation simulations on different yield criterions are determined. And, the identified stress–strain curves are compared with experimental data.

This paper discovers that indentation on Yld2004 is able to fully identify difference in equivalent plastic strain between 0° and 90° directions when indentation depth h_t is relatively smaller. And, this research demonstrates conclusively that plastic anisotropy identified from indentation on Yld2004 and Yld91 is more applicable at larger strains than that on Hill48, and that on Yld2004 is more applicable than that on Yld91, overall. In addition, the method on the determination of the locations of the optimal data acquisition points is demonstrated to be also valid for anisotropic material.

This paper first investigates plastic anisotropic properties and laws identified from indentation simulations following more advanced anisotropic yield criterions and provides reference for later research.