• Journal of Mechanical Science and Technology
• /
• v.14 no.6
• /
• pp.666-674
• /
• 2000

Large deflection dynamic responses of laminated composite cylindrical shells under impact are analyzed by the geometrically nonlinear finite element method based on a generalized Sander's shell theory with the first order transverse shear deformation and the von-Karman large deflection assumption. A modified indentation law with inelastic indentation is employed for the contact force. The nonlinear finite element equations of motion of shell and an impactor along with the contact laws are solved numerically using Newmark's time marching integration scheme in conjunction with Akay type successive iteration in each step. The ply failure region of the laminated shell is estimated using the Tsai- Wu quadratic interaction criteria. Numerical results, including the contact force histories, deflections and strains are presented and compared with the ones by linear analysis. The effect of the radius of curvature on the composite shell behaviors is investigated and discussed.

• Structural Engineering and Mechanics
• /
• v.28 no.5
• /
• pp.587-601
• /
• 2008

In this study, a numerical procedure based on the finite element method for materially and geometrically nonlinear analysis of reinforced and prestressed concrete slender columns with arbitrary section subjected to combined biaxial bending and axial load is developed. In order to overcome the low computer efficiency of the conventional section integration method in which the reinforced concrete section is divided into a large number of small areas, an efficient section integration method is used to determine the section tangent stiffness. In this method, the arbitrary shaped cross section is divided into several concrete trapezoids according to boundary vertices, and the contribution of each trapezoid to section stiffness is determined by integrating directly the trapezoid. The space frame flexural theory is utilized to derive the element tangent stiffness matrix. The nonlinear full-range member response is traced by an updated normal plane arc-length solution method. The analytical results agree well with the experimental ones.

• Computational Structural Engineering
• /
• v.11 no.1
• /
• pp.179-190
• /
• 1998

A geometrically nonlinear finite element formulation of spatial cable networks is presented using two cable elements. Firstly, derivation procedures of tangent stiffness and mass matrices for the space truss element and the elastic catenary cable element are summarized. The load incremental method based on Newton-Raphson iteration method and the dynamic relaxation method are presented in order to determine the initial static state of cable nets subjected to self-weights and support motions. Furthermore, static non-linear analysis of cable structures under additional live loads are performed based on the initial configuration. Challenging example problems are presented and discussed in order to demonstrate the feasibility of the present finite element method and investigate static nonlinear behaviors of cable nets.

• Journal of the Society of Naval Architects of Korea
• /
• v.29 no.2
• /
• pp.132-139
• /
• 1992

This paper presents geometrically nonlinear formulation of shell problems using the three-dimensional curved shell element, which includs large displacements and large rotations. Formulations of the geometrically nonlinear problems can be derived in a variety of ways, but most of them have been obtained by assuming that nodal rotations are small. Hence, the tangent stiffness matrix is derived under the assumptions that rotational increments are infinitesimal and the effect of finite rotational increments have to be considered during the equilibrium iterations. To study the large displacement and large rotation problems, the restrictions are removed and the formulations of the curved shell element including the effect of large rotational increments are developed in this paper. The displacement based finite element method using this improved formulation are applied to the analyses of the geometrically nonlinear behaviors of the single and double curved shells, which are compared with the results by others.

• Journal of the Korea institute for structural maintenance and inspection
• /
• v.15 no.2
• /
• pp.125-135
• /
• 2011

The latest study for formulation of finite element method and computation techniques has progressed widely. The classical method in the formulation of frame elements for geometrically nonlinear analysis derives the geometric stiffness directly from the governing differential equation for bending with axial force. From the computational viewpoint of this paper, the most common approach is the finite element method. Commonly, the formulation of frame elements for geometrically nonlinear structures is based on appropriate interpolation functions for the transverse and axial displacements of the member. The formulation of flexibility-based elements, on the other hand, is based on interpolation functions for the internal forces. In this paper, a new method is used to suppose that interpolation functions for the displacements from the curvatures is Lagrangian interpolation. This paper derives flexibility matrix from that displacement functions and is considered the application of it. Using the flexibility matrix, this paper apply the program considered geometrically nonlinear analysis to common problems.

• Structural Engineering and Mechanics
• /
• v.68 no.3
• /
• pp.299-312
• /
• 2018

In this paper, an extended finite element method is proposed to analyze both geometric and material non-linear behavior of general Functionally Graded Material (FGM) plate-shell type structures. A user defined subroutine (UMAT) is developed and implemented in Abaqus/Standard to study the elastoplastic behavior of the ceramic particle-reinforced metal-matrix FGM plates-shells. The standard quadrilateral 4-nodes shell element with three rotational and three translational degrees of freedom per node, S4, is extended in the present study, to deal with elasto-plastic analysis of geometrically non-linear FGM plate-shell structures. The elastoplastic material properties are assumed to vary smoothly through the thickness of the plate-shell type structures. The nonlinear approach is based on Mori-Tanaka model to underline micromechanics and locally determine the effective FGM properties and self-consistent method of Suquet for the homogenization of the stress-field. The elasto-plastic behavior of the ceramic/metal FGM is assumed to follow Ludwik hardening law. An incremental formulation of the elasto-plastic constitutive relation is developed to predict the tangent operator. In order to to highlight the effectiveness and the accuracy of the present finite element procedure, numerical examples of geometrically non-linear elastoplastic functionally graded plates and shells are presented. The effects of the geometrical parameters and the volume fraction index on nonlinear responses are performed.

• Steel and Composite Structures
• /
• v.18 no.1
• /
• pp.51-69
• /
• 2015

In recent years, the use of fiber reinforced polymer composites has increased because of their unique features. They have been used widely in the aircraft and space industries, medical and sporting goods and automotive industries. Thanks to their beneficial and various advantages over traditional materials such as high strength, high rigidity, low weight, corrosion resistance, low maintenance cost, aesthetic appearance and easy demountable or moveable construction. In this paper, it is aimed to determine and compare the geometrically nonlinear static and dynamic analysis results of footbridges using steel and glass fiber reinforced polymer composite (GFRP) materials. For this purpose, Halgavor suspension footbridge is selected as numerical examples. The analyses are performed using three identical footbridges, first constructed from steel, second built only with GFRP material and third made of steel- GFRP material, under static and dynamic loadings using finite element method. In the finite element modeling and analyses, SAP2000 program is used. Geometric nonlinearities are taken into consideration in the analysis using P-Delta criterion. The numerical results have indicated that the responses of the three bridges are different and that the response values obtained for the GFRP composite bridge are quite less compared to the steel bridge. It is understood that GFRP material is more useful than the steel for the footbridges.

• Steel and Composite Structures
• /
• v.35 no.1
• /
• pp.77-92
• /
• 2020

The present paper outlined a procedure for geometrically nonlinear dynamic analysis of functionally graded graphene platelets-reinforced (GPLR-FG) nanocomposite cylinder subjected to mechanical shock loading. The governing equation of motion for large deformation problems is derived using meshless local Petrov-Galerkin (MLPG) method based on total lagrangian approach. In the MLPG method, the radial point interpolation technique is employed to construct the shape functions. A micromechanical model based on the Halpin-Tsai model and rule of mixture is used for formulation the nonlinear functionally graded distribution of GPLs in polymer matrix of composites. Energy dissipation in analyses of the structure responding to dynamic loads is considered using the Rayleigh damping. The Newmark-Newton/Raphson method which is an incremental-iterative approach is implemented to solve the nonlinear dynamic equations. The results of the proposed method for homogenous material are compared with the finite element ones. A very good agreement is achieved between the MLPG and FEM with very fine meshing. In addition, the results have demonstrated that the MLPG method is more effective method compared with the FEM for very large deformation problems due to avoiding mesh distortion issues. Finally, the effect of GPLs distribution on strength, stiffness and dynamic characteristics of the cylinder are discussed in details. The obtained results show that the distribution of GPLs changed the mechanical properties, so a classification of different types and volume fraction exponent is established. Indeed by comparing the obtained results, the best compromise of nanocomposite cylinder is determined in terms of mechanical and dynamic properties for different load patterns. All these applications have shown that the present MLPG method is very effective for geometrically nonlinear analyses of GPLR-FG nanocomposite cylinder because of vanishing mesh distortion issue in large deformation problems. In addition, since in proposed method the distributed nodes are used for discretization the problem domain (rather than the meshing), modeling the functionally graded media yields to more accurate results.

• Steel and Composite Structures
• /
• v.15 no.5
• /
• pp.481-505
• /
• 2013

This paper focuses on thermal post-buckling analysis of functionally graded beams with temperature dependent physical properties by using the total Lagrangian Timoshenko beam element approximation. Material properties of the beam change in the thickness direction according to a power-law function. The beam is clamped at both ends. In the case of beams with immovable ends, temperature rise causes compressible forces and therefore buckling and post-buckling phenomena occurs. It is known that post-buckling problems are geometrically nonlinear problems. Also, the material properties (Young's modulus, coefficient of thermal expansion, yield stress) are temperature dependent: That is the coefficients of the governing equations are not constant in this study. This situation suggests the physical nonlinearity of the problem. Hence, the considered problem is both geometrically and physically nonlinear. The considered highly non-linear problem is solved considering full geometric non-linearity by using incremental displacement-based finite element method in conjunction with Newton-Raphson iteration method. In this study, the differences between temperature dependent and independent physical properties are investigated for functionally graded beams in detail in post-buckling case. With the effects of material gradient property and thermal load, the relationships between deflections, critical buckling temperature and maximum stresses of the beams are illustrated in detail in post-buckling case.

• Structural Engineering and Mechanics
• /
• v.44 no.1
• /
• pp.109-125
• /
• 2012

Post-buckling behavior of Timoshenko beams subjected to uniform temperature rising with temperature dependent physical properties are studied in this paper by using the total Lagrangian Timoshenko beam element approximation. The beam is clamped at both ends. In the case of beams with immovable ends, temperature rise causes compressible forces end therefore buckling and post-buckling phenomena occurs. It is known that post-buckling problems are geometrically nonlinear problems. Also, the material properties (Young's modulus, coefficient of thermal expansion, yield stress) are temperature dependent: That is the coefficients of the governing equations are not constant in this study. This situation suggests the physical nonlinearity of the problem. Hence, the considered problem is both geometrically and physically nonlinear. The considered highly non-linear problem is solved considering full geometric non-linearity by using incremental displacement-based finite element method in conjunction with Newton-Raphson iteration method. The beams considered in numerical examples are made of Austenitic Stainless Steel (316). The convergence studies are made. In this study, the difference between temperature dependent and independent physical properties are investigated in detail in post-buckling case. The relationships between deflections, thermal post-buckling configuration, critical buckling temperature, maximum stresses of the beams and temperature rising are illustrated in detail in post-buckling case.