• Structural Engineering and Mechanics
• /
• v.29 no.3
• /
• pp.289-300
• /
• 2008

The use of metric trial functions to represent the real stress field in what is called the unsymmetric finite element formulation is an effective way to improve predictions from distorted finite elements. This approach works surprisingly well because the use of parametric functions for the test functions satisfies the continuity conditions while the use of metric (Cartesian) shape functions for the trial functions attempts to ensure that the stress representation during finite element computation can retrieve in a best-fit manner, the actual variation of stress in the metric space. However, the issue of how to handle situations where there is locking along with mesh distortion has never been addressed. In this paper, we show that the use of a consistent definition of the constrained strain field in the metric space can ensure a lock-free solution even when there is mesh distortion. The three-noded Timoshenko beam element is used to illustrate the principles. Some significant conclusions are drawn regarding the optimal strategy for finite element modelling where distortion effects and field-consistency requirements have to be reconciled simultaneously.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.18 no.4 s.70
• /
• pp.385-393
• /
• 2005

This paper established the dynamic model of a flexible Timoshenko beam capable of geometrical nonlinearities subject to large overall motions by using the finite element method. Equations of motion are derived by using Hamilton principle and are formulated in terms of finite elements in which the nonlinear constraint equations are adjoined to the system using Lagrange multipliers. The Newmark direct integration method and the Newton-Raphson iteration are employed here for the numerical study which is to demonstrate the efficiency of the proposed formulation.

• Journal of Korean Society of Steel Construction
• /
• v.19 no.6
• /
• pp.599-610
• /
• 2007

In this paper, we described a co-rotational formulation for the geometrical nonlinear analysis of three-dimensional frames. We suggested a new concept called the Zero-Twist-Section Condition (ZTSC) to decide the element coordinate system consistently. According to the ZTSC procedure, it is possible to obtain an element coordinate system and natural deformations consistently when finite displacements and rotations are induced in an element. Based on the developed procedure, numerical examples are investigated to calculate natural rotations while finite displacements are imposed on an element. Also, the developed co-rotational procedure gives accurate results in the analysis of post-buckling problems with finite rotations.

• International Journal of Aerospace System Engineering
• /
• v.3 no.1
• /
• pp.10-16
• /
• 2016

A trigonometric Layerwise higher order shear deformation theory (TLHSDT) is developed and implemented for free vibration and buckling analysis of laminated composite and sandwich plates by analytical and finite element formulation. The present model assumes parabolic variation of out-plane stresses through the depth of the plate and also accomplish the zero transverse shear stresses over the surface of the plate. Thus a need of shear correction factor is obviated. The present zigzag model able to meet the transverse shear stress continuity and zigzag form of in-plane displacement continuity at the plate interfaces. Hence, botheration of shear correction coefficient is neglected. In the case of analytical method, the governing differential equation and boundary conditions are obtained from the principle of virtual work. For the finite element formulation, an efficient eight noded $C^0$ continuous isoparametric serendipity element is established and employed to examine the dynamic analysis. Like FSDT, the considered mathematical model possesses similar number of variables and which decides the present models computationally more effective. Several numerical predictions are carried out and results are compared with those of other existing numerical approaches.

• Nuclear Engineering and Technology
• /
• v.50 no.6
• /
• pp.971-980
• /
• 2018

Using the concept of quasi-static decomposition and using three-noded isoparametric locking-free element, this article presents a formulation of the finite element method for Timoshenko beam subjected to spatially different time-dependent motions at supports. To verify the validity of the formulation, three fixed-hinged beams excited by the real seismic motions are examined; one is a slender beam, another is a stocky one, and the other is an intermediate one. The numerical results of time histories of motions of the three beams are compared with corresponding analytical solutions. The internal loads such as bending moment and shearing force at a specific time are also compared with analytic solutions. These comparisons show good agreements. The comparisons between static components of the internal loads and the corresponding total internal loads show that the static components predominate in the stocky beam, whereas the dynamic components predominate in the slender one. Thus, the total internal loads of the stocky beam, which is governed by static components, can be predicted simply by static analysis. Careful numerical experiments indicate that the fundamental frequency of a beam can be used as a parameter identifying such a stocky beam.

• Smart Structures and Systems
• /
• v.13 no.4
• /
• pp.517-546
• /
• 2014

This paper proposes a refined electro-mechanical beam formulation. Lagrange-type polynomials are used to interpolate the unknowns over the beam cross section. Three- (L3), four- (L4), and nine-point(L9) polynomials are considered which lead to linear, bi-linear, and quadratic displacement field approximations over the beam cross-section. Finite elements are obtained by employing the principle of virtual displacements in conjunction with the Carrera Unified Formulation (CUF). The finite element matrices and vectors are expressed in terms of fundamental nuclei whose forms do not depend on the assumptions made. Additional refined beam models are implemented by introducing further discretizations, over the beam cross-section. Some assessments from bibliography have been solved in order to validate the electro-mechanical formulation. The investigations conducted show that the present formulation is able to detect the electro-mechanical interaction.

• Structural Engineering and Mechanics
• /
• v.62 no.1
• /
• pp.33-42
• /
• 2017

We present in this paper a finite element formulation for nonlinear torsional analysis of 3D beams with arbitrary composite cross-sections. Since the proposed formulation employs a continuum mechanics based beam element with kinematics enriched by the extended St. Venant solutions, it can precisely account higher order warping effect and its 3D couplings. We propose a numerical procedure to calculate the extended St. Venant equation and the twisting center of an arbitrary composite cross-section simultaneously. The accuracy and efficiency of the proposed formulation are thoroughly investigated through representative numerical examples.

• Journal of KSNVE
• /
• v.6 no.3
• /
• pp.287-296
• /
• 1996

A virtual work form of flexible multibody dynamic formulation with rotary inertia has been derived. For the analysis of large flexible multibody systems, deformation modal coordinates have been employed to represent coupled motion between gross and vibrational motion. For the efficient evaluation of the entries in the mass matrix, a flexible body has been treated as a collection of mass points. The rotary inertia was generated from the consistent mass matrix in a finite element model. Deformation mode shapes were obtained from finite element analysis. Bending and twisting vibration analyses of a cantilever have been carried out to see rotary inertia effects. A space flexible robot simulation has been also carried out to show effectiveness of the proposed formulation. This formulation is effective to the model that consists of beam, plate, or shell element that contains rotational degree of freedom at the nodal point. It is also effective to the flexible body model to which a large lumped rotary inertia is attached.

• Coupled systems mechanics
• /
• v.6 no.3
• /
• pp.351-367
• /
• 2017

In this article, the multiphysics response of magneto-electro-elastic (MEE) cantilever beam subjected to thermo-mechanical loading is analysed. The equilibrium equations of the system are obtained with the aid of the principle of total potential energy. The constitutive equations of a MEE material accounting the thermal fields are used for analysis. The corresponding finite element (FE) formulation is derived and model of the beam is generated using an eight noded 3D brick element. The 3D FE formulation developed enables the representation of governing equations in all three axes, achieving accurate results. Also, geometric, constitutive and loading assumptions required to dimensionality reduction can be avoided. Numerical evaluation is performed on the basis of the derived formulation and the influence of various mechanical loading profiles and volume fractions on the direct quantities and stresses is evaluated. In addition, an attempt has been made to compare the individual effect of thermal and mechanical loading with the combined effect. It is believed that the numerical results obtained helps in accurate design and development of sensors and actuators.

• Journal of The Korean Society of Agricultural Engineers
• /
• v.47 no.4
• /
• pp.15-24
• /
• 2005

This study is focused on modeling to predict the behavior of two-way RC slabs. A new finite element model will be presented to analyze the nonlinear behavior of RC slabs. The numerical approach is based on the p-version degenerate shell element including theory of anisotropic laminated composites, theory of materially and geometrically nonlinear plates. In the nonlinear formulation of this model, the total Lagrangian formulation is adopted with large deflections and moderate rotations being accounted for in the sense of von Karman hypothesis. The material model is based on the Kuper's yield criterion, hardening rule, and crushing condition. The validity of the proposed p-version nonlinear RC finite element model is demonstrated through the load-deflection curves and the ultimate loads. It is shown that the proposed model is able to adequately predict the deflection and ultimate load of two-way slabs with respect to steel arrangements and steel ratios.