• Coupled systems mechanics
• /
• v.2 no.4
• /
• pp.375-387
• /
• 2013

Soils consist of an assemblage of particles with different sizes and shapes which form a skeleton whose voids are filled with water and air. Hence, soil behaviour must be analyzed by incorporating the effects of the transient flow of the pore-fluid through the voids, and therefore requires a two-phase continuum formulation for saturated porous media. The present paper presents briefly the Biot's basic theory of dynamics of saturated porous media with u-P formulation to determine the responses of pore fluid and soil skeleton during cyclic loading. Kelvin elements are attached to transmitting boundary. The Pastor-Zienkiewicz-Chan model has been used to describe the inelastic behavior of soils under isotropic cyclic loadings. Newmark-Beta method is employed to discretize the time domain. The response of fluid-saturated porous media which are subjected to time dependent loads has been simulated numerically to predict the liquefaction potential of a semi-infinite saturated sandy layer using finite-infinite elements. A settlement of 17.1 cm is observed at top surface. It is also noticed that liquefaction occurs at shallow depth. The mathematical advantage of the coupled finite element analysis is that the excess pore pressure and displacement can be evaluated simultaneously without using any empirical relationship.

• Coupled systems mechanics
• /
• v.3 no.3
• /
• pp.305-318
• /
• 2014

The study deals with the physical modeling of a typical building frame resting on pile foundation and embedded in cohesive soil mass using complete three-dimensional finite element analysis. Two different pile groups comprising four piles ($2{\times}2$) and nine piles ($3{\times}3$) are considered. Further, three different pile diameters along with the various pile spacings are considered. The elements of the superstructure frame and those of the pile foundation are descretized using twenty-node isoparametric continuum elements. The interface between the pile and pile and soil is idealized using sixteen-node isoparametric surface elements. The current study is an improved version of finite element modeling for the soil elements compared to the one reported in the literature (Chore and Ingle 2008). The soil elements are discretized using eight-, nine- and twelve-node continuum elements. Both the elements of superstructure and substructure (i.e., foundation) including soil are assumed to remain in the elastic state at all the time. The interaction analysis is carried out using sub-structure approach in the parametric study. The total stress analysis is carried out considering the immediate behaviour of the soil. The effect of various parameters of the pile foundation such as spacing in a group and number piles in a group, along with pile diameter, is evaluated on the response of superstructure. The response includes the displacement at the top of the frame and bending moment in columns. The soil-structure interaction effect is found to increase displacement in the range of 58 -152% and increase the absolute maximum positive and negative moments in the column in the range of 14-15% and 26-28%, respectively. The effect of the soil- structure interaction is observed to be significant for the configuration of the pile groups and the soil considered in the present study.

• Proceedings of the KSME Conference
• /
• 2000.04a
• /
• pp.351-356
• /
• 2000

Various transition elements are generally used for the effective analysis of a complicated mechanical structure. In this paper, a solid-to-beam transition finite element which connects a continuum element and a $c^1-continuity$ beam element each other is proposed. The shape functions of the transition finite elements, which a 8-noded hexahedral solid element fur 3D analysis and a 4-noded quadrilateral plane element fur 2D analysis are connected to a Euler's beam element, are explicitely formulated. In order to show the effectiveness and convergence characteristics of the proposed transition elements. numerical tests are performed for various examples and their results are compared with those obtained by other methods. As the result of this study. following conclusions are obtained: (1)The proposed transition finite elements show the monotonic convergence characteristics because of having used the compatible displacement folds. (2)As being used the transition element in the finite element analysis, the finite element modelings are more convenient and the analysis results are more accurate because of the formulation characteristies of the Euler's beam element.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.8 no.2
• /
• pp.25-32
• /
• 1988

A formulation of an isoparametric quadrilateral higher-order plate bending finite element is presented. The 8-noded 28-d.o.f. plate element has been degenerated from the three-dimensional continuum by introducing the plate assumptions and considering higher-order in-plane displacement profile. The element characteristics have been derived by the Galerkin's weighted residual method and computed by using the selective reduced integration technique to avoid shear-locking phenomenon. Several numerical examples are given to demonstrate the accuracy and versatility of the proposed quadrilateral higher-order plate bending element over the other existing plate finite elements in both static and dynamic analyses.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.11 no.1
• /
• pp.82-87
• /
• 1987

In this study a method of shape optimization is applied to two dimensional heat transfer system. For this the optimization problem is defined in a functional form including cost, constraints and the system governing equation. Then the material derivative concept in continuum mechanics and the adjoint variable method are employed for the shape design sensitivity analysis. With the sensitivity analysis results, an optimum is sought with the gradient projection optimization algorithm. The two dimensional isoparametric finite elements are used for accurate analysis and sensitivity calculations. The above method is employed to find the boundary shape to achieve a desired temperature distribution along a segment of the boundary subject to the maximum area constraint.