• Transactions of the Korean Society of Mechanical Engineers
• /
• v.14 no.2
• /
• pp.310-323
• /
• 1990

A three-noded, with three degree-of-freedom at each node, in-plane curved beam element is formulated and employed in eigen-analysis of constant curvature beam. The conventional quadratic shape functions used in a three noded C .deg. type curved beam element produce such an undesirable large stiffness that a significant error is introduced in displacements and stresses. These phenomena are called 'Stiffness Locking Phenomena', which result from spurious strain energy due to inappropriate assumptions on independent isoparametric quadratic interpolation functions. Stiffness locking phenomena can be alleviated by using modified interpolation functions which get rid of spurious constraints of conventional interpolation functions. Eigenvalues and their modes as well as displacements and stresses may be locked because they are related to stiffness. Using modified curved beam element in eigenvalue problem of cantilever and arch, the property and performance of modified curved beam element are examined by numerical experimentations. In these eigen-analyses, mass matrices are calculated by using both modified and unmodified curved beam element, are compared with theoretical solutions. These comparisons show that the performance of the modified curved beam element is better than that of the unmodified curved beam element.

• KCI Concrete Journal
• /
• v.13 no.1
• /
• pp.19-25
• /
• 2001

Based on the orthotropic hypoelasticity formulation, a constitutive material model of concrete taking account of triaxial stress state is presented. In this model, the ultimate strength surface of concrete in triaxial stress space is described by the Hsieh's four-parameter surface. On the other hand, the different ultimate strength surface of concrete in strain space is proposed in order to account for increasing ductility in high confinement pressure. Compressive ascending and descending behavior of concrete is considered. Concrete cracking behavior is considered as a smeared crack model, and after cracking, the tensile strain-softening behavior and the shear mechanism of cracked concrete are considered. The proposed constitutive model of concrete is compared with some results obtained from tests under the states of uniaxial, biaxial, and triaxial stresses. In triaxial compressive tests, the peak compressive stress from the predicted results agrees well with the experimental results, and ductility response under high confining pressure matches well the experimental result. The reinforcing bars embedded in concrete are considered as an isoparametric line element which could be easily incorporated into the isoparametric solid element of concrete, and the average stress - average strain relationship of the bar embedded in concrete is considered. From numerical examples for a reinforced concrete simple beam and a structural beam type member, the stress state of concrete in the vicinity of talc critical region is investigated.

• Structural Engineering and Mechanics
• /
• v.86 no.3
• /
• pp.291-309
• /
• 2023

This paper addresses the finite element modeling of functionally graded porous (FGP) beams for free vibration and buckling behaviour cases. The formulated finite element is based on simple and efficient higher order shear deformation theory. The key feature of this formulation is that it deals with Euler-Bernoulli beam theory with only three unknowns without requiring any shear correction factor. In fact, the presented two-noded beam element has three degrees of freedom per node, and the discrete model guarantees the interelement continuity by using both C⁰ and C¹ continuities for the displacement field and its first derivative shape functions, respectively. The weak form of the governing equations is obtained from the Hamilton principle of FGP beams to generate the elementary stiffness, geometric, and mass matrices. By deploying the isoparametric coordinate system, the derived elementary matrices are computed using the Gauss quadrature rule. To overcome the shear-locking phenomenon, the reduced integration technique is used for the shear strain energy. Furthermore, the effect of porosity distribution patterns on the free vibration and buckling behaviours of porous functionally graded beams in various parameters is investigated. The obtained results extend and improve those predicted previously by alternative existing theories, in which significant parameters such as material distribution, geometrical configuration, boundary conditions, and porosity distributions are considered and discussed in detailed numerical comparisons. Determining the impacts of these parameters on natural frequencies and critical buckling loads play an essential role in the manufacturing process of such materials and their related mechanical modeling in aerospace, nuclear, civil, and other structures.

• Proceedings of the Korean Geotechical Society Conference
• /
• 1991.10a
• /
• pp.289-311
• /
• 1991

Analysis of supported excavation system by Elasto-Plastic Isoparametric Finite Element Method and Elasto-Plastic Beam Method have been conducted for the simulation of stepwise underground excavation. Conventional methods, fixed Supported Beam and Spring Supported Beam method, also have been examined and compared with the results of elasto-plastic beam method and field data. Except unavoidable result of upward ground settlement near the top of retaining wall and relatively high bending moment of wall at each excavation level, satisfactory results have been derived using elasto-plastic isopara metric finite element method. The results from elasto-plastic beam analysis program, developed by the author, are proved to be fit field data in acceptable variance as shown in the paper. Displacement and bending moment, of the wall by conventional methods, both fixed supported beam and spring supported beam, are always underestimated than field data, and attention must be given that the diffence increases with deeper excavation depth and lower horizontal subgrade reaction of the ground.

• 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.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2000.10a
• /
• pp.369-376
• /
• 2000

For the general loading condition and boundary condition, it is very difficult to obtain closed-form solutions for buckling loads and natural frequencies of thin-walled structures because its behaviour is very complex due to the coupling effect of bending and torsional behaviour. Consequently most of previous finite element formulations introduced approximate displacement fields using shape functions as Hermitian polynomials, isoparametric interpoation function, and so on. The purpose of this study is to calculate the exact displacement field of a thin-walled straight beam element with the non-symmetric cross section and present a consistent derivation of the exact dynamic stiffness matrix. An exact dynamic element stiffness matrix is established from Vlasov's coupled differential equations for a uniform beam element of non-symmetric thin-walled cross section. This numerical technique is accomplished via a generalized linear eigenvalue problem by introducing 14 displacement parameters and a system of linear algebraic equations with complex matrices. The natural frequencies are evaluated for the non-symmetric thin-walled straight beam structure, and the results are compared with available solutions in order to verify validity and accuracy of the proposed procedures.

• Structural Engineering and Mechanics
• /
• v.22 no.4
• /
• pp.483-510
• /
• 2006

The dynamic instability characteristics of laminated composite stiffened shell panels subjected to in-plane harmonic edge loading are investigated in this paper. The eight-noded isoparametric degenerated shell element and a compatible three-noded curved beam element are used to model the shell panels and the stiffeners respectively. As the usual formulation of degenerated beam element is found to overestimate the torsional rigidity, an attempt has been made to reformulate it in an efficient manner. Moreover the new formulation for the beam element requires five degrees of freedom per node as that of shell element. The method of Hill's infinite determinant is applied to analyze the dynamic instability regions. Numerical results are presented to demonstrate the effects of various parameters like shell geometry, lamination scheme, stiffening scheme, static and dynamic load factors and boundary conditions, on the dynamic instability behaviour of laminated composite stiffened panels subjected to in-plane harmonic loads along the boundaries. The results of free vibration and buckling of the laminated composite stiffened curved panels are also presented.

• Structural Engineering and Mechanics
• /
• v.33 no.4
• /
• pp.447-484
• /
• 2009

The spatially coupled buckling, in-plane, and lateral bucking analyses of thin-walled Timoshenko curved beam with non-symmetric, double-, and mono-symmetric cross-sections resting on elastic foundation are performed based on series solutions. The stiffness matrices are derived rigorously using the homogeneous form of the simultaneous ordinary differential equations. The present beam formulation includes the mechanical characteristics such as the non-symmetric cross-section, the thickness-curvature effect, the shear effects due to bending and restrained warping, the second-order terms of semitangential rotation, the Wagner effect, and the foundation effects. The equilibrium equations and force-deformation relationships are derived from the energy principle and expressions for displacement parameters are derived based on power series expansions of displacement components. Finally the element stiffness matrix is determined using force-deformation relationships. In order to verify the accuracy and validity of this study, the numerical solutions by the proposed method are presented and compared with the finite element solutions using the classical isoparametric curved beam elements and other researchers' analytical solutions.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.7 no.4
• /
• pp.1-13
• /
• 2003

For spatial free vibration of non-symmetric thin-walled curved beams with shear deformation, an improved formulation is proposed in the present study. The elastic strain and the kinetic energies are first derived by considering constant curvature and shear deformation effects due to shear forces and restrained warping torsion. Next equilibrium equations and force-deformation relations are obtained using a stationary condition of total potential energy. And the finite element procedures are developed by using isoparametric curved beam element with arbitray thin-walled sections. Particularly not only shear deformation and thickness-curvature effects on vibration behaviors of curved beams but also mode transition and crossover phenomena with change in curvatures of beams are parametrically investigated. In order to illustrate the accuracy and the reliability of this study, various numerical solutions for spatial free vibration are compared with results by available references and ABAQUS's shell element.

• Computers and Concrete
• /
• v.12 no.4
• /
• pp.443-498
• /
• 2013

A detailed finite element modeling is presented for the simulation of the nonlinear behavior of reinforced concrete structures which manages to predict the nonlinear behavior of four different experimental setups with computational efficiency, robustness and accuracy. The proposed modeling method uses 8-node hexahedral isoparametric elements for the discretization of concrete. Steel rebars may have any orientation inside the solid concrete elements allowing the simulation of longitudinal as well as transverse reinforcement. Concrete cracking is treated with the smeared crack approach, while steel reinforcement is modeled with the natural beam-column flexibility-based element that takes into consideration shear and bending stiffness. The performance of the proposed modeling is demonstrated by comparing the numerical predictions with existing experimental and numerical results in the literature as well as with those of a commercial code. The results show that the proposed refined simulation predicts accurately the nonlinear inelastic behavior of reinforced concrete structures achieving numerical robustness and computational efficiency.