• Structural Engineering and Mechanics
• /
• v.6 no.2
• /
• pp.173-183
• /
• 1998

In the shape design of flexible structures, it is useful to predict the initial shape from the desirable large deformed shapes under some loading conditions. In this paper, we present a numerical procedure of an initial shape determination problem for hyperelastic materials which enables us to calculate an initial shape corresponding to the prescribed deformed shape and boundary condition. The present procedure is based on an Arbitrary Lagrangian-Eulerian (ALE) finite element method for hyperelasticity, in which arbitrary change of shapes in both the initial and deformed states can be treated by considering the variation of geometric mappings in the equilibrium equation. Then the determination problem of the initial shape can be formulated as a nonlinear problem to solve the unknown initial shape for the specified deformed shape that satisfies the equilibrium equation. The present approach can be implemented easily to the finite element method by employing the isoparametric hypothesis. Some basic numerical results are also given to characterize the present procedure.

• Structural Engineering and Mechanics
• /
• v.29 no.2
• /
• pp.203-221
• /
• 2008

In this paper, we present an idea of the geometry-dependent MITC method. The simple concept is exemplified to improve a 2-node iso-beam (isoparametric beam) finite element of varying section. We first study the behavior of a standard 2-node iso-beam finite element of prismatic section, which has been widely used with reduced integration (or the equivalent MITC method) in order to avoid shear locking. Based on analytical studies on cantilever beams of varying section, we propose the axial strain correction (ASC) scheme and the geometry-dependent tying (GDT) scheme for the 2-node iso-beam element. We numerically analyze varying section beam problems and present the improved performance by using both ASC and GDT schemes.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.1
• /
• pp.29-38
• /
• 1994

Linear composite theory as well as a finite element program is developed for axisymmetric elastomeric bearings. This study is limited to axisymmetrically loaded horizontal layered systems with linear, elastic, small' deformation conditions. A multiscale method is used in the development of the composite theory which enables us to model inhomogeneous layered composites as equivalent homogeneous, orthotropic material. Only continuity of the prime variables is required for the finite element analysis, allowing the use of simple $C_o$ elements whereas rather complicated theories presented in the past need more requirements. Four node isoparametric elements are used in the study. The developed theory of this paper is limited to linear conditions, however, the analysis can be extended to nonlinear behavior of flexible material in elastomeric bearing by using multiscale method presented here. Two numerical examples are examined and compared to the results of discrete and previously obtained composite analysis to verify the theory.

• Magazine of the Korea Concrete Institute
• /
• v.8 no.3
• /
• pp.187-195
• /
• 1996

A finite element formulation based on the CFT(Compression Field Theory), considering the effect of compression softening in cracked concrete, and macro-scopic and rotating crack models etc., was presented for the nonlinear behaviour of structural concrete. Considering the computational efficency and the ability of modelling the post-ultimate behaviour as major concerns, the Incremental displacement solution algorithm involving initial material stiffnesses and the relaxation procedure for fast convergence was adopted and formulated in a type of 8-noded quadrilateral isoparametric elements. The analysis program NASCOM(Non1inear Analysis of Structural Concrete by FEM : Monotonic Loading) developed in this way enables the predictions of strength and deformation capacities in a full range, crack patterns and their corresponding widths, and yield extents of reinforcement. As the verification purpose of NASCOM, the predictions were made for Bhide's Panel(PB21) and Leonhardt's deep beam tests. The predicted results shows somewhat stiff behaviour for the panel test, and vice versa for deep beam tests. More refining process would be necessary hereafter in terms of more accurately simulating the effects of tension-stiffening and compression softening in concrete.

• Journal of the Korean Geotechnical Society
• /
• v.36 no.12
• /
• pp.27-34
• /
• 2020

Use of axisymmetric shell element for the structure increases the efficiency and accuracy in finite element analysis of the interaction between the ground and the structure. This paper derived the force balance equation and the moment balance equation for an axisymmetric shell element based on Kirchhoff's theory. The governing equation for the axial deformation used the isoparametric shape function in the Galerkin formulation, and the governing equation for the shell bending used the higher-order shape function. The developed axisymmetric shell element was combined with Geo-COUS, a geotechnical finite element program for the coupled analysis with the ground. The accuracy of the developed element was confirmed through the example analyses of the circular plate and the liquid storage tank. And the energy balance equation for the axisymmetric shell element is presented.

• 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.38 no.3
• /
• pp.261-281
• /
• 2011

A new method of formulation of a class of elements that are immune to mesh distortion effects is proposed here. The simple three-noded bar element with an offset of the internal node from the element center is employed here to demonstrate the method and the principles on which it is founded upon. Using the function space approach, the modified formulation is shown here to be superior to the conventional isoparametric version of the element since it satisfies the completeness requirement as the metric formulation, and yet it is in agreement with the best-fit paradigm in both the metric and the parametric domains. Furthermore, the element error is limited to only those that are permissible by the classical projection theorem of strains and stresses. Unlike its conventional counterpart, the modified element is thus not prone to any errors from mesh distortion. The element formulation is symmetric and thus satisfies the requirement of the conservative nature of problems associated with all self-adjoint differential operators. The present paper indicates that a proper mapping set for distortion immune elements constitutes geometric and displacement interpolations through parametric and metric shape functions respectively, with the metric components in the displacement/strain replaced by the equivalent geometric interpolation in parametric co-ordinates.

• Structural Engineering and Mechanics
• /
• v.19 no.5
• /
• pp.477-501
• /
• 2005

In this study, a new smart beam finite element is proposed for the finite element modeling of beam-type smart structures that are equipped with bonded plate-type piezoelectric sensors and actuators. Constitutive equations for the direct piezoelectric effect and converse piezoelectric effect of piezoelectric materials are considered in the formulation. By using a variational principle, the equations of motion for the smart beam finite element are derived. The proposed 2-node beam finite element is an isoparametric element based on Timoshenko beam theory. The proposed smart beam finite element is applied to the free vibration control adopting a constant gain feedback scheme. The electrical force vector, which is obtained in deriving an equation of motion, is the control force equivalent to that in existing literature. Validity of the proposed element is shown through comparing the analytical results of the verification examples with those of other previous researchers. With the use of smart beam finite elements, simulation of free vibration control is demonstrated by sensing the voltage of the piezoelectric sensors and by applying the voltages to the piezoelectric actuators.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.26 no.5
• /
• pp.904-913
• /
• 2002

In this paper, the dynamic stress intensity factors of fracture mechanics are numerically computed in time domain using the FEM. For which the finite element formulations are derived applying the Galerkin method to the equations of motion in convolution integral as has been presented in the previous paper. To assure the strain fields of r$^{-1}$ 2/ singularity near the crack tip, the triangular quarter-point singular elements are imbedded in the finite element mesh discretized by the isoparametric quadratic quadrilateral elements. Two-dimensional problems of the elastodynamic fracture mechanics under the impact load are solved and compared with the existing numerical and analytical solutions, being shown that numerical results of good accuracy are obtained by the presented method.

• Structural Engineering and Mechanics
• /
• v.6 no.7
• /
• pp.817-826
• /
• 1998

A finite element model based on the incremental endochronic theory for flexible pavement materials was developed in this study. Three grid systems with eight-node cubic isoparametric elements, and different loading steps were used to perform the calculations for a specimen of circular cylinder. The uniaxial stress experimental results on an asphalt mixture at $60^{\circ}C$ in SHRP conducted by University of California at Berkeley were used to check the ability of the derived numerical model. Then, the numerical results showed isotropic response and deviatoric response on the specimen in a three dimensional manner, which provided a better understanding for a deformed flexible material under the specified loading conditions.