• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1998.10a
• /
• pp.463-470
• /
• 1998

The boundary/finite element analysis for the seismic wave amplifications due to nonhomogeneous alluvial deposits was performed in this study. For numerical analysis, the homogeneous linear elastic soil half-space was modeled by using the 3-node isoparametric boundary elements and the inhomogeneous alluvial deposit was modeled by using the 8-node isoparametric finite elements. The two elements at interface were coupled together by the equilibrium condition of the tractions and the compatibility condition of the displacements. As a prarmetric variable, the incident angle and the dimensionsless frequency of the SH, P and SV-waves and the shear wave velocity ratio and the mass density ratio between the half-space and the alluvial deposit were selected.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.25 no.5
• /
• pp.413-420
• /
• 2012

In this study, the Kernel integration scheme for 2D linear elastic direct boundary element method has been discussed on the basis of subparametric element. Usually, the isoparametric based boundary element uses same polynomial order in the both basis function and mapping function. On the other hand, the order of mapping function is lower than the order of basis function to define displacement field when the subparametric concept is used. While the logarithmic numerical integration is generally used to calculate Kernel integration as well as Cauchy principal value approach, new formulation has been derived to improve the accuracy of numerical solution by algebraic modification. The subparametric based direct boundary element has been applied to 2D elliptical partial differential equation, especially for plane stress/strain problems, to demonstrate whether the proposed algebraic expression for integration of singular Kernel function is robust and accurate. The problems including cantilever beam and square plate with a cutout have been tested since those are typical examples of simple connected and multi connected region cases. It is noted that the number of DOFs has been drastically reduced to keep same degree of accuracy in comparison with the conventional isoparametric based BEM. It is expected that the subparametric based BEM associated with singular Kernel function integration scheme may be extended to not only subparametric high order boundary element but also subparametric high order dual boundary 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.73 no.6
• /
• pp.715-724
• /
• 2020

The influence on flexural strength of Glass/Epoxy laminated composite curved panels of different geometries (cylindrical, spherical, elliptical, hyperboloid and flat) due to inclusion of nano cenosphere filler examined in this research article. The deflection responses of the hybrid structure are evaluated numerically using the isoparametric finite element technique and modelled mathematically via higher-order displacement structural kinematics. To predict the deflection values, a customised in-house computer code in MATLAB environment is prepared using the higher-order isoparametric formulation. Subsequently, the numerical model validity has been established by comparing with those of available benchmark solution including the convergence characteristics of the finite element solution. Further, a few cenosphere filled hybrid composite are prepared for different volume fractions for the experimental purpose, to review the propose model accuracy. The experimental deflection values are compared with the finite element solutions, where the experimental elastic properties are adopted for the computation. Finally, the effect of different variable design dependent parameter and the percentages of nano cenosphere including the geometrical shapes obtained via a set of numerical experimentation.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1993.10a
• /
• pp.1-7
• /
• 1993

Most of wall or plate structures subjected to nonuniform earth or water pressure has variable thickness. These problems were generally solved by models with uniform thickness and pressure. To obtain more accurate and economic solution for this type of problem. efficient isoparametric plate element considering variable thickness and nonuniform pressure were developed. Some example problems demonstrated efficiency of the proposed element.

• Journal of Mechanical Science and Technology
• /
• v.15 no.11
• /
• pp.1499-1506
• /
• 2001

Various transition elements are used in general for the effective finite element analysis of complicated mechanical structures. In this paper, a solid-to-beam transition finite element, which can b e used for connecting a C1-continuity beam element to a continuum solid element, is proposed. The shape functions of the transition finite element are derived to meet the compatibility condition, and a transition element equation is formulated by the conventional finite element procedure. In order to show the effectiveness and convergence characteristics of the proposed transition element, numerical tests are performed for various examples. As a result of this study, following conclusions are obtained. (1) The proposed transition element, which meets the compatibility of the primary variables, exhibits excellent accuracy. (2) In case of using the proposed transition element, the number of nodes in the finite element model may be considerably reduced and the model construction becomes more convenient. (3) This formulation method can be applied to the usage of higher order elements.

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

• Transactions of Materials Processing
• /
• v.4 no.1
• /
• pp.69-78
• /
• 1995

The arbitrary Lagrangian-Eulerian(ALE) finite element analysis is applied to the axisymmetric hot extrusion through continuous dies. In order to simulate hot forming problems, an ALE scheme for temperature analysis is proposed. The computed results are compared with experimental results as with those by pure Lagrangian method. In the present study mesh control is accomplished by the use of isoparametric mapping of quadrilaterals.

• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 1999.04a
• /
• pp.292-299
• /
• 1999

In order to investigate dynamic behaviors of marine cables under wave and earthquake forces a geometric nonlinear. F, E formulation of marine cables is presented and tangent stiffness and mass matrices for the isoparametric cable element are derived, The initial equilibrium state of cables subjected to self -weights and current forces is determined and free vibration and dynamic nonlinear analysis of cable structures under additional environmental loads are performed based on the initial configuration Challenging examples are presented and discussed in order to demonstrate the feasibility of the present finite element method and investigate dynamic nonlinear behaviors of marine cables.

• Journal of the Korean Mathematical Society
• /
• v.53 no.6
• /
• pp.1411-1429
• /
• 2016

We consider the nite volume element method for elliptic problems using isoparametric bilinear elements on quadrilateral meshes. A gradient recovery method is presented by using the patch interpolation technique. Based on some superclose estimates, we prove that the recovered gradient $R({\nabla}u_h)$ possesses the superconvergence: ${\parallel}{\nabla}u-R({\nabla}u_h){\parallel}=O(h^2){\parallel}u{\parallel}_3$. Finally, some numerical examples are provided to illustrate our theoretical analysis.