• Journal of the Korean Geotechnical Society
• /
• v.38 no.8
• /
• pp.29-37
• /
• 2022

A three-dimensional analysis is required to interpret the drainage behavior of an improved ground with vertical drains, and the macroelement method enables efficient interpretation considering the three-dimensional drainage effect of vertical drains under two-dimensional plane strain condition. In this study, a novel finite element analysis program was developed by applying the macroelement method to the vacuum consolidation method used in ground improvement practice. The conventional macroelement method was used to calculate the amount of drainage from the vertical drain by setting the excess porewater pressure in the drainage material to zero; however, the program developed in this study was improved to consider negative excess porewater pressure as an actual vacuum consolidation condition. To verify the performance of the program, because of a comparison with the measurement values at the site where the vacuum consolidation method was applied, results predicted by the program and field measurement data showed similar settlement behavior.

• Nuclear Engineering and Technology
• /
• v.41 no.3
• /
• pp.335-346
• /
• 2009

This paper deals with a hydroelastic vibration analysis of two rectangular plates partially coupled with a liquid, which is bounded by two plates and two rigid side walls. The wet displacement of each plate is assumed to be a combination of the modal functions of a dry uniform beam with a clamped boundary condition. As the liquid is assumed to be an ideal liquid, the displacement potential satisfying the Laplace equation is determined so that the liquid boundary conditions can meet the requirements at the rigid surfaces and the free liquid surface. The wet dynamic modal functions of each plate are expanded by using the finite Fourier transform to obtain an appropriate form of the compatibility requirement along the contacting surfaces between the plates and the liquid. The liquid-coupled natural frequencies of the plates are derived by using the Rayleigh-Ritz method. Finite element analyses using commercial software are carried out to verify the proposed theory. It is observed that the theoretical method agrees excellently with the three-dimensional finite element analyses results. The effects of the liquid depth and the liquid thickness on the normalized natural frequencies are investigated to identify the dynamic characteristics of the liquid coupled system.

• Structural Engineering and Mechanics
• /
• v.11 no.6
• /
• pp.651-662
• /
• 2001

In order to investigate the seismic behavior and seismic design methods in the transverse direction of a shield tunnel, a series of model shaking table tests and a two-dimensional finite element dynamic analysis on the tests are carried out. Two kinds of static analytical methods based on ground-tunnel composite finite element model and beam-spring element model are proposed, and the validity of the static analyses is verified by model shaking table tests. The investigation concerns the dynamic response behavior of a tunnel and the ground, the interaction between the tunnel and ground, and an evaluation of different seismic design methods. Results of the investigation indicate that the shield tunnel follows the surrounding ground in displacement and dynamic characteristics in the transverse direction; also, the static analytical methods proposed by the authors can be used directly as the seismic design methods in the transverse direction of a shield tunnel.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.1
• /
• pp.1-10
• /
• 1991

Two-dimensional explicit finite element code was developed. The transient dynamics code can analyse large deformations of non-linear materials subjected to extremely high strain rates. The Lagrangian finite element program uses an explicit time integration operator to integrate the equations of motion, thus the stiffness matrix is not introduced. Cylinder upsetting and ring compression problems are simulated to check the effects of friction and inertia force. It is shown that (1) calculated results agree very well with experimental results, (2) constant shear friction method overestimates the decrease of inner ring radius and then underestimates after on in comparison with the Coulomb friction method, and (3) the effect of the increase in initial strain rate is similar to the effect of higher frictional coefficient.

• Advanced Composite Materials
• /
• v.18 no.1
• /
• pp.43-59
• /
• 2009

The current study assessed the effect of a bolt hole on tensile fatigue properties of CFRP laminates. Two specimens, i.e. $[(0/90)_3]S$, $[(0/45/90/-45)_2]_S$, were analyzed using a finite element method and were experimentally tested for cases, both with and without a hole, whose diameter corresponded to 0.12 times the specimen width. Delamination positions predicted by a 3-dimensional static finite element analysis were matched well to those observed by an ultrasonic imaging system in the middle of fatigue test. A hole whose diameter corresponds to 0.12 times the specimen width caused the fatigue strength to decrease by 9% and 11% under 5 Hz loading frequency, and by 22% and 25% under 10 Hz loading frequency for $[(0/90)_3]_S$ and $[(0/45/90/-45)_2]_S$, respectively. Because the decrease in sectional area due to the hole was normalized in calculation of the tensile strength, a stress concentration around the hole is believed to induce the strength degradation of fatigue specimens. From the finite element analyses, the stress concentration factor around a hole was expected as 8.8 and 9.5 for $[(0/90)_3]_S$ and $[(0/45/90/-45)_2]_S$, respectively.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.209-212
• /
• 2002

Since Berkhoff proposed the mild-slope equation in 1972, it has widely been used for calculation of shallow water wave transformation. Recently, it was extended to give an extended mild-slope equation, which includes the bottom slope squared term and bottom curvature term so as to be capable of modeling wave transformation on rapidly varying topography. These equations were derived by integrating the Laplace equation vertically. In the present study, we develop a finite element model to solve the Laplace equation directly while keeping the same computational efficiency as the mild-slope equation. This model assumes the vertical variation of wave potential as a cosine hyperbolic function as done in the derivation of the mild-slope equation, and the Galerkin method is used to discretize . The computational domain was discretized with proper finite elements, while the radiation condition at infinity was treated by introducing the concept of an infinite element. The upper boundary condition can be either free surface or a solid structure. The applicability of the developed model was verified through example analyses of two-dimensional wave reflection and transmission. .

• Geomechanics and Engineering
• /
• v.34 no.4
• /
• pp.341-380
• /
• 2023

This study evaluates the seismic fragilities in fill slopes in South Korea through parametric finite element analyses that have been barely investigated thus far. We consider three slope geometries for a slope of height 10 m and three slope angles, and two soil types, namely frictional and frictionless, associated with two soil states, loose and dense for frictional soils and soft and stiff for frictionless soils. The input ground motions accounting for four site conditions in South Korea are obtained from one-dimensional site response analyses. By comparing the numerical modeling of slopes using PLAXIS^2D against the previous studies, we compiled suites of the maximum permanent slope displacement (D_max) against two ground motion parameters, namely, peak ground acceleration (PGA) and Arias Intensity (I_A). A probabilistic seismic demand model is adopted to compute the probabilities of exceeding three limit states (minor, moderate, and extensive). We propose multiple seismic fragility curves as functions of a single ground motion parameter and numerous seismic fragility surfaces as functions of two ground motion parameters. The results show that soil type, slope angle, and input ground motion influence these probabilities, and are expected to help regional authorities and engineers assess the seismic fragility of fill slopes in the road systems in South Korea.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.14 no.5
• /
• pp.1-12
• /
• 2010

This study presents a finite element analysis method that can accurately evaluate the nonlinear behaviour of structures affected by shallow soft subsoils and the soil-structure interaction. A two-dimensional finite element model that consists of a structure and shallow soft subsoil was used. The finite element model was used for a nonlinear time domain analysis of the OpenSees program. A parametric study was performed to investigate the effects of soil shear velocities, earthquake input motions, soft soil depth, and soil-structure interaction. The result of the proposed nonlinear finite element analysis method was compared with the result of an existing frequency domain analysis method, which is frequently used for addressing nonlinear soil behavior. The result showed that the frequency domain analysis, which uses equivalent secant soil stiffness and does not address the soil-structure interaction, significantly overestimated the response of the structures with short dynamic periods. The effect of the soil-structure interaction on the response spectrum did not significantly vary with the foundation dimensions and structure mass.

• Journal of Ocean Engineering and Technology
• /
• v.12 no.2 s.28
• /
• pp.57-64
• /
• 1998

This study is concerned with criterion for membrane to shell conversion in two-dimensional elastic-plastic finite element analysis using membrane/shell mixed element. It is well known that in the sheet metal forming some parts of the sheet deform under almost pure stretching (membrane) conditions, whereas other parts in contact with sharp tooling surfaces can develop significant bending strains. The membrane analysis has a short computational time however, in the membrane analysis the bending effects can not be condidered at all. On the other hand, the shell analysis allows the consideration of bending effects, but involves too much computational time. So Onatel),2), Yang et al3),4) developed the membrane/shell mixed element. Onate introduced the energy ratio parameter and Yang et al introduced the ratio of thickness to radius of curvature as the criterion. In the present study we propose a new criterion by using the angle between both side elements in the nodal point.

• International Journal of Aeronautical and Space Sciences
• /
• v.16 no.2
• /
• pp.206-222
• /
• 2015

This work uses different finite element approaches to the free vibration analysis of reinforced shell structures, and a simplified model of a typical launcher with two boosters is used as an example. The results obtained using a refined one-dimensional (1D) beam model are compared to those obtained with commercial finite element software. The 1D models that are used in the present work are based on the Carrera Unified Formulation (CUF), which assumes a variable kinematic displacement field over the cross-sections of the beam. Two different sets of polynomials that correspond to Taylor (TE) or Lagrange (LE) expansions were used. The analyses focused on three reinforced structures: a stiffened panel, a reinforced cylinder and the complete structure of the launcher. The frequencies and natural modes obtained using one-dimensional models are compared to those obtained from classical finite element analysis. The classical FE models were built using a beam-shell or solid elements, and the results indicate that the refined beam models can in fact be used to investigate the behavior of very complex reinforced structures. These models can predict the shell-like modes that are typical of thin-walled structures that cannot be detected using classical beam models. The refined 1D models used in the present work provide results that are as accurate as those from solid FE models, but the 1D models have a much lower computational cost.