• Journal of Korean Society of Steel Construction
• /
• v.11 no.2 s.39
• /
• pp.153-165
• /
• 1999

In order to trace the lateral-torsional post-bucking behaviors of thin-walled space frames with non-symmetric cross sections, a geometrically non-linear finite element formulation is presented by applying incremental equilibrium equations based on the updated Lagrangian formulation and introducing Vlasov's assumption. The improved displacement field for non-symmetric thin-walled cross sections is introduced based on inclusion of second order terms of finite rotations, and the potential energy corresponding to the semitangential rotations and moments is consistently derived. For finite element analysis, tangent stiffness matrices of thin-walled space frame element are derived by using the Hermition polynomials as shape functions. A co-rotational formulation in order to evaluate the unbalanced loads is presented by separating the rigid body rotations and pure deformations from incremental displacements and evaluating the updated direction cosines and incremental member forces.

• Computational Structural Engineering
• /
• v.10 no.1
• /
• pp.201-211
• /
• 1997

A clearly consistent finite element formulation for geometrically non-linear analysis of space frames is presented by applying incremental equilibrium equations based on the updated Lagrangian formulation and introducing Vlasov's assumption. The improved displacement field for symmetric cross sections is introduced based on inclusion of second order terms of finite rotations, and the potential energy corresponding to the semitangential rotations and moments is consistently derived. For finite element analysis, elastic and geometric stiffness matrices of the space frame element are derived by using the Hermitian polynomials as shape functions. A co-rotational formulation in order to evaluate the unbalanced loads is presented by separating the rigid body rotations and pure deformations from incremental displacements and evaluating the updated direction cosines of the frame element due to rigid body rotations and incremental member forces from pure deformaions. Finite element solutions for the spatial buckling and post-buckling analysis of space frames are compared with available solutions and other researcher's results.

• Structural Engineering and Mechanics
• /
• v.67 no.6
• /
• pp.629-641
• /
• 2018

The paper describes the development of a two-dimensional (2D) co-rotational nonlinear beam finite element that includes advanced path-following capabilities for detecting bifurcation instability in elasto-plasticity of steel elements subjected to fire without introducing imperfections. The advantage is twofold: i) no need to assume the magnitude of the imperfections and consequent reduction of the model complexity; ii) the presence of possible critical points is checked at each converged time step based on the actual load and stiffness distribution in the structure that is affected by the temperature field in the elements. In this way, the buckling modes at elevated temperature, that may be different from the ones at ambient temperature, can be properly taken into account. Moreover, an improved displacement predictor for estimating the displacement field allowed significant reduction of the computational cost. A co-rotational framework was exploited for describing the beam kinematic. In order to highlight the potential practical implications of the developed finite element, a parametric analysis was performed to investigate how the beam element compares both with the EN1993-1-2 buckling curve and with experimental tests on axially compressed steel members. Validation against experimental data and numerical outcomes obtained with commercial software is thoroughly described.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.25 no.3
• /
• pp.195-202
• /
• 2012

Two nonlinear frame elements taking into account geometric nonlinearity is presented and compared based on the Lagrangian co-rotational formulation. The first frame element is believed to be geometrically-exact because not only tangent stiffness matrices is exactly evaluated including stiffness matrices due to initial deformation but also total member forces are directly determined from total deformations in the deformed state. Particularly two exact tangent stiffness matrices based on total Lagrangian and updated Lagrangian formulation, respectively, are verified to be identical. In the second frame element, the deformed curved shape is regarded as the polygon and current flexural deformations in iteration process are neglected in evaluating tangent stiffness matrices and total member forces. Two numerical examples are given to demonstrate the accuracy and the good performance of the first frame element compared with the second element. Furthermore it is shown that the first frame element can be used in tracing nonlinear behaviors of cable members.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.48 no.8
• /
• pp.565-571
• /
• 2020

In this paper, precision and efficient nonlinear structural analysis for the aircraft wing-box model is developed. Herein, nonlinear shell element based on the co-rotational (CR) formulation is implemented. Then, parallel computing algorithm, the element-based partitioning technique is developed to accelerate the computational efficiency of the nonlinear structural analysis. Finally, computational performance, i.e., accuracy and efficiency, of the proposed analysis is evaluated by comparing with that of the existing commercial software.

• Computers and Concrete
• /
• v.31 no.3
• /
• pp.223-239
• /
• 2023

Geometric nonlinear performance simulation and analysis of complex modern buildings and industrial products require high-performance shell elements. Balancing multiple aspects of performance in the one geometric nonlinear analysis element remains challenging. We present a new shell element, flat shell DKMGQ-CR (Co-rotational Discrete Kirchhoff-Mindlin Generalized Conforming Quadrilateral), for linear and geometric nonlinear analysis of both thick and thin shells. The DKMGQ-CR shell element was developed by combining the advantages of high-performance membrane and plate elements in a unified coordinate system and introducing the co-rotational formulation to adapt to large deformation analysis. The effectiveness of linear and geometric nonlinear analysis by DKMGQ-CR is verified through the tests of several classical numerical benchmarks. The computational results show that the proposed new element adapts to mesh distortion and effectively alleviates shear and membrane locking problems in linear and geometric nonlinear analysis. Furthermore, the DKMGQ-CR demonstrates high performance in analyzing thick and thin shells. The proposed element DKMGQ-CR is expected to provide an accurate, efficient, and convenient tool for the geometric nonlinear analysis of shells.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.26 no.3A
• /
• pp.447-455
• /
• 2006

A general purpose of 4-node co-rotational resultant shell element is developed for the solution of nonlinear problems of reinforced concrete, steel and fiber-reinforced composite structures. The formulation of the geometrical stiffness presented here is defined on the mid-surface by using the second order kinematic relations and is efficient for analyzing thick plates and shells by incorporating bending moment and transverse shear resultant forces. The present element is free of shear locking behavior by using the ANS (Assumed Natural Strain) method such that the element performs very well as thin shells. Inelastic behaviour of concrete material is based on the plasticity with strain hardening and elasto-plastic fracture model. The plasticity of steel is based on Von-Mises Yield and Ivanov Yield criteria with strain hardening. The transverse shear stiffness of laminate composite is defined by an equilibrium approach instead of using the shear correction factor. The proposed formulation is computationally efficient and versitile for most civil engineering application and the test results showed good agreement.

• Computational Structural Engineering
• /
• v.10 no.4
• /
• pp.205-216
• /
• 1997

In order to trace the lateral post-buckling behaviors of thin-wafled space frames, a geometrically nonlinear finite element formulation is presented by applying incremental equilibrium equations based on the updated Lagrangian formulation and introducing Vlasov's assumption. The improved displacement field for symmetric thin-walled cross sections is introduced based on inclusion of second order terms of finite rotations, and the potential energy corresponding to the semitangential rotations and moments is consistently derived. For finite element analysis, tangent stiffness matrices of the thinwalled space frame element with 7 degrees of freedom including the restrained warping for each node are derived by using the Hermition polynomials as shape functions. A co-rotational formulation in order to evaluate the unbalanced loads is presented by separating the rigid body rotations and pure deformations from incremental displacements and evaluating the updated direction cosines of the frame element due to rigid body rotations and incremental member forces from pure deformations. Finite element solutions for the spatial buckling and post-buckling analysis of thin-walled space frames are presented and compared with available solutions and other researcher's results.

• Structural Engineering and Mechanics
• /
• v.49 no.6
• /
• pp.727-743
• /
• 2014

Geometrically nonlinear analysis of planar beam and frame structures made of functionally graded material (FGM) by using the finite element method is presented. The material property of the structures is assumed to be graded in the thickness direction by a power law distribution. A nonlinear beam element based on Bernoulli beam theory, taking the shift of the neutral axis position into account, is formulated in the context of the co-rotational formulation. The nonlinear equilibrium equations are solved by using the incremental/iterative procedure in a combination with the arc-length control method. Numerical examples show that the formulated element is capable to give accurate results by using just several elements. The influence of the material inhomogeneity in the geometrically nonlinear behavior of the FGM beam and frame structures is examined and highlighted.

• Journal of Ocean Engineering and Technology
• /
• v.37 no.6
• /
• pp.256-265
• /
• 2023

Many numerical methods have been developed since 1961, but unresolved issues remain. This study developed a numerical method to address these issues and determine the coefficients and properties of rotational waves with a shear current in a finite water depth. The number of unknown constants was reduced significantly by introducing a wavelength-independent coordinate system. The reference depth was calculated independently using the shooting method. Therefore, there was no need for partial derivatives with respect to the wavelength and the reference depth, which simplified the numerical formulation. This method had less than half of the unknown constants of the other method because Newton's method only determines the coefficients. The breaking limit was calculated for verification, and the result agreed with the Miche formula. The water particle velocities were calculated, and the results were consistent with the experimental data. Dispersion relations were calculated, and the results are consistent with other numerical findings. The convergence of this method was examined. Although the required series order was reduced significantly, the total error was smaller, with a faster convergence speed.