Search results

The purpose of this study is to demonstrate that using appropriate values for fictitious parameters is very important in dynamic relation methods. It will be shown that a better scheme can be made by modifying these terms.

Former research studies have proposed diverse values for fictitious parameters. These factors are very essential and highly affect structural analyses’ abilities. In this paper, the fictitious masses in ten previous well-known schemes are replaced with each other. These formulations lead to the extra 41 different new procedures.

To compare the skills of the created processes with those of the ten previous ones, 14 benchmark problems with geometrical nonlinear behaviour are analysed. The performances’ evaluations are based on the number of iterations and analysis time. Considering these two criteria, the score of each technique is found for the ranking assessments.

To solve a static problem by using a dynamic relaxation (DR) scheme, it should be first converted to a dynamic space. Using the appropriate values for fictitious terms is very important in this approach. The fictitious mass matrix and damping factor play the most effective role in the process stability. Besides, the fictitious time step is necessary for improving the method convergence rate.

Different famous DR procedures were compared with each other previously. These solvers used their original assumptions for the imaginary mass and damping. So far, no attempt has been made to change the fictitious parameters of the well-known DR methods. As these fictitious factors highly affect structural analyses’ efficiencies, these solvers are formulated again by using new parameters. In this study, the fictitious masses of ten previous famous methods are replaced with each other. These substitutions give 51 different procedures.

It is concluded that the present formulations lead to more effective and favourable methods than the solvers with previous assumptions.

This paper aims to propose a new robust membrane finite element for the analysis of plane problems. The suggested element has triangular geometry. Four nodes and 11 degrees of freedom (DOF) are considered for the element. Each of the three vertex nodes has three DOF, two displacements and one drilling. The fourth node that is located inside the element has only two translational DOF.

The suggested formulation is based on the assumed strain method and satisfies both compatibility and equilibrium conditions within each element. This establishment results in higher insensitivity to the mesh distortion. Enforcement of the equilibrium condition to the assumed strain field leads to considerably high accuracy of the developed formulation.

To show the merits of the suggested plane element, its different properties, including insensitivity to mesh distortion, particularly under transverse shear forces, immunities to the various locking phenomena and convergence of the element are studied. The obtained results demonstrate the superiority of the suggested element compared with many of the available robust membrane elements.

According to the attained results, the proposed element performs better than the well-known displacement-based elements such as linear strain triangular element, Q4 and Q8 and even is comparable with robust modified membrane elements.

In this paper, quasi‐Tresca yield surfaces are reviewed. In order to do elasto‐plastic analysis, a new yield criterion is presented. The proposed yield surface can be used in nonlinear three‐dimensional analysis of structures. Function of the yield surface is presented in principal stress space and also Cartesian one. A computer program has been developed for nonlinear analysis in C++. Numerical examples have been solved by the proposed yield surface and good results have been obtained.

Numerous finite elements are proposed based on analytical solutions. However, it is difficult to find the solutions for complicated governing equations. This paper aims to present a novel formulation in the framework of assumed stress quasi-conforming method for the static and free vibration analysis of anisotropic and symmetric laminated plates.

Firstly, an initial stress approximation ruled by 17 parameters, which satisfies the equilibrium equations is derived to improve the performance of the constructed element. Then the stress matrix is treated as the weighted function to weaken the strain-displacement equations. Finally, the Timoshenko’s laminated composite beam functions are adopted as boundary string-net functions for strain integration.

Several numerical examples are presented to show the performance of the new element, and the results obtained are compared with other available ones. Numerical results have proved that the new element is free from shear locking and possesses high accuracy for the analysis of anisotropic and symmetric laminated plates.

This paper proposes a new QC element for the static and free vibration analysis of anisotropic and symmetric laminated plates. In contrast with the complicated analytical solutions of the equilibrium equations, an initial stress approximation ruled by 17 parameters is adopted here. The Timoshenkos laminated composite beam functions are introduced as boundary string-net functions for strain integration. Numerical results demonstrate the new element is free from shear locking and possesses high accuracy for the analysis of anisotropic and symmetric laminated plates.

Numerical experiences reveal that the performances of the dynamic relaxation (DR) method are related to the structural types. This paper is devoted to compare the DR schemes for geometric nonlinear analysis of shells. To achieve this task, 12 famous approaches are briefly introduced. The differences among these schemes are between the estimation of the time step, the mass and the damping matrices. In this study, several benchmark structures are analyzed by using these 12 techniques. Based on the number of iterations and the analysis duration, their performances are graded. Numerical findings reveal the high efficiency of the kinetic DR (kdDR) approach and Underwood’s strategy.

Up to now, the performances of various DR algorithms for geometric nonlinear analysis of thin shells have not been investigated. In this paper, 12 famous DR methods have been used for solving these structures. It should be noted that the difference between these approaches is in the estimation of the fictitious parameters. The aforementioned techniques are used to solve several numerical samples. Then, the performances of all schemes are graded based on the number of iterations and the analysis duration.

The final ranking of each strategy will be obtained after studying all numerical examples. It is worth emphasizing that the number of iterations and that of convergence points of the arc length algorithms are dependent on the value of the initial arc length. In other words, a slight change in the magnitude of the arc length may lead to the wrong responses. Contrary to this behavior, the analyzer’s role in the dynamic relaxation techniques is considerably less than the arc length method. In the DR strategies when the answer approaches the limit points, the iteration number increases automatically. As a result, this algorithm can be used to analyze the structures with complex equilibrium paths.

Numerical experiences reveal that the DR method performances are related to the structural types. This paper is devoted to compare the DR schemes for geometric nonlinear analysis of shells.

Geometric nonlinear analysis of shells is a sophisticated procedure. Consequently, extensive research studies have been conducted to analyze the shells efficiently. The most important characteristic of these structures is their high resistance against pressure. This study demonstrates the performances of various DR methods in solving shell structures.

Up to now, the performances of various DR algorithms for geometric nonlinear analysis of thin shells are not investigated.

The purpose of this study is dedicated to use an efficient mixed strain finite element approach to develop a three-node triangular shell element. Moreover, large deformation analysis of the functionally graded material shells is the main contribution of this research. These target structures include thin or moderately thick panels.

Due to reach these goals, Green–Lagrange strain formulation with respect to small strains and large deformations with finite rotations is used. First, an efficient three-node triangular degenerated shell element is formulated using tensorial components of two-dimensional shell theory. Then, the variation of Young’s modulus through the thickness of shell is formulated by using power function. Note that the change of Poisson’s ratio is ignored. Finally, the governing linearized incremental relation was iteratively solved using a high potential nonlinear solution method entitled generalized displacement control.

Some well-known problems are solved to validate the proposed formulations. The suggested triangular shell element can obtain the exact responses of functionally graded (FG) shell structures, without any shear locking, instabilities and ill-conditioning, even by using fewer numbers of the elements. The obtained outcomes are compared with the other reference solutions. All findings demonstrate the accuracy and capability of authors’ element for analyzing FG shell structures.

A mixed strain finite element approach is used for nonlinear analysis of FG shells. These structures are curved thin and moderately thick shells. Small strains and large deformations with finite rotations are assumed.

FG shells are mostly made curved thin or moderately thick, and these structures have a lot of applications in the civil and mechanical engineering.

The social implication of this study is concerned with how technology impacts the world. In short, the presented scheme can improve structural analysis ways.

Developing an efficient three-node triangular element, for geometrically nonlinear analysis of FG doubly-curved thin and moderately thick shells, is the main contribution of the current research. Finite rotations are considered by using the Taylor’s expansion of the rotation matrix. Mixed interpolation of strain fields is used to alleviate the locking phenomena. Using fewer numbers of shell elements with fewer numbers of degrees of freedom can reduce the computational costs and errors significantly.

The purpose of this paper is to suggest a robust hybrid method for updating the stress and plastic internal variables in plasticity considering damage mechanics.

By benefiting the properties of the well-known explicit and implicit integrations, a new mixed method is derived. In fact, the advantages of the mentioned techniques are used to achieve an efficient integration.

The numerical studies demonstrate the high precision and robustness of the suggested algorithm.

The perfect von-Mises plasticity together with Lemaitre damage model is considered within the realm of small deformations.

Updating stress and plastic internal variables are of utmost importance in elastoplastic analyses of structures. The accuracy and efficiency of stress-updating methods significantly affect the final outcomes of nonlinear analyses.

The idea which is used to derive the hybrid method leads to an efficient integration method for updating the constitutive equations of the damage mechanics.

The purpose of this paper is to present a new effective integration method for cyclic plasticity models.

By defining an integrating factor and an augmented stress vector, the system of differential equations of the constitutive model is converted into a nonlinear dynamical system, which could be solved by an exponential map algorithm.

The numerical tests show the robustness and high efficiency of the proposed integration scheme.

The von‐Mises yield criterion in the regime of small deformation is assumed. In addition, the model obeys a general nonlinear kinematic hardening and an exponential isotropic hardening.

Integrating the constitutive equations in order to update the material state is one of the most important steps in a nonlinear finite element analysis. The accuracy of the integration method could directly influence the result of the elastoplastic analyses.

The paper deals with integrating the constitutive equations in a nonlinear finite element analysis. This subject could be interesting for the academy as well as industry. The proposed exponential‐based integration method is more efficient than the classical strategies.

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.

This paper aims to provide a simple and accurate higher order predictor‐corrector integration which can be used in dynamic analysis and to compare it with the previous works.

The predictor‐corrector integration is defined by combining the higher order explicit and implicit integrations in which displacement and velocity are assumed to be functions of accelerations of several previous time steps. By studying the accuracy and stability conditions, the weighted factors and acceptable time step are determined.

Simplicity and vector operations plus accuracy and stability are the main specifications of the new predictor‐corrector method. This procedure can be used in linear and nonlinear dynamic analysis.

In the proposed integration, time step is assumed to be constant.

The numerical integration is the heart of a dynamic analysis. The result's accuracy is strongly influenced by the accuracy and stability of the numerical integration.

This paper presents simple and accurate predictor‐corrector integration based on accelerations of several previous time steps. This may be used as a routine in any dynamic analysis software to enhance accuracy and reduce computational time.