• Proceedings of the Korean Geotechical Society Conference
• /
• 2004.03b
• /
• pp.348-355
• /
• 2004

본 연구에서는 geosynthetics로 보강된 다짐토의 보강 메카니즘을 파악하기 위한 목적으로 실내시험 및 수치계산 수행하였다. 본 연구에서 고려하는 보강 메카니즘은 전단에 의한 다짐토의 체적 팽창(부의 다일렌탄시)을 geosynthetics에 의해 구속 억제하는 과정에서 생성되는 효과로 생각한다. 먼저, 실내실험을 위한 구제직인 방법으로서, geosynthetics의 보강효과를 정량직으로 파악하기 위하여 사질토를 다짐하여 공시체를 만들어 그 주위에 geosynthetics를 설치하여 전체적으로 압축전단 시험을 실시하였다. 또한, 다짐토의 다짐도를 달리 하고 한 가지 종류만의 geosynthetics를 이용하여, 다짐토와 geosynthetics의 상호작용에 따른 압축력 변화, geosynthetics의 인장력 변화 및 공시체의 파괴 진행상황 등을 살펴보았다. 수치계산에서는 다짐토의 다일렌탄시 특성에 대하여 표현 가능한 탄소성 구성모델을 이용하였다. 또한, 탄소성 구성 모델에서의 항복 이전의 탄성영역의 거동을 묘사하기 위하여 Hashiguchi(1989)가 제안한 subloading surface의 개념을 도입하였고, 유한요소(FEM)해석을 통해 얻어진 결과들을 실내시험의 결과와 비교 분석하였으며, 그 결과 양자 양호한 결과를 얻었다.

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

The current design equations for ultimate strength of tubular joints are based on a limited number of experimental results performed on simple joints with simple loading conditions and depend on value of the branch to the chord diameter- ratio $\beta$ too much. Therefore, the purpose of this study is to estimate the ultimate strength of CHS tilbular joints considering the effects of branch inclination angles $\theta$, chord length to diametel ratio $\alpha$ and chord end conditions by finite element analysis. The analyses are performed using finite element software ADINA that is capable of modeling elasto-plastic material behavior as well as geometric nonlinearities. The results show that the current use of sin $\theta$ in normalized design equations for inclined branches is reasonable, but somewhat conservative. When compared with the previous experimental database, the close numerical results are obtained from the parametric studies on the static strength of T-, Y-, DT- and X-joints. Also, a new design equation for ultimate stregth of CHS tubular joints is derived using a modified version of the ring model which can include the effects of $\alpha$ and chord end condtion.

• Journal of the Korean Society for Precision Engineering
• /
• v.15 no.9
• /
• pp.117-126
• /
• 1998

For investigating rubber pad sheet metal forming process, the rubber pad deformation characteristics as well as the contact problem of rubber pad-sheet metal has been analyzed. In this paper, the behavior of the rubber deformation is represented by hyper-elastic constitutive relations based on a generalized Mooney-Rivlin model. Finite element procedures for the two-dimensional responses, employing total Lagrangian formulations are implemented in an implicit form. The volumetric incompressibility condition of the rubber deformation is included in the formulation by using penalty method. The sheet metal is characterized by elasto-plastic material with strain hardening effect and analyzed by a commercial code. The contact procedure and interface program between rubber pad and sheet metal are implemented. Inflation experiment of circular rubber pad identifies the behaviour of the rubber pad deformation during the process. The various form dies and scaled down apparatus of the rubber-pad forming process are fabricated for simulating realistic forming process. The obtaining experimental data and FEM solutions were compared. The numerical solutions illustrate fair agreement with experimental results. The forming pressure distribution according to the dimensions of sheet metal and rubber pads, various rubber models and rubber material are also compared and discussed.

• Journal of the Korean Professional Engineers Association
• /
• v.20 no.3
• /
• pp.5-14
• /
• 1987

연속체의 해석에 있어서, 특별한 경우를 제외하고는, 구조물의 개략적인 거동을 파악해야 될 경우가 종종 있다. 이러한 요구에 부응하기 위해서 강체요소법(Rigid Element Method)이라 불리우는 새로운 해석법이 개발되었다. 강체요소법은 원래 평정연구실에서 벽식프리캐스트 철근콘크리트 구조물의 탄소성해석을 하기 위해서 개발된 해석법에 착안하여, 내수벽과 같은 연속체에 적용함으로서 시작된 수치해석법이다. 그 후 저자들은 도통쉘, 구형쉘 혹은 이들이 조합된 쉘구조물에 적용할 수 있도록 개발 확장하였다. 강체요소법의 기본개념은 연속체의 분해된 각 요소를 강체(rigid body)라고 가정하고, 각 요소들은 요소의 강성으로 치환된 가상스프링으로 서로 연결되어 있다고 가정하여, 이 가상스프링의 거동을 평가함으로서 전체구조물의 거동을 파악하는 해석법이다. 이때 요소의 주변에 취해진 스프링은 해석을 단순화하기 위해서 축력, 면내전단력 및 면외전단력만을 전달한다고 가정하고, 요소의 강체변위(자유도)는 요소내의 임의의 한 점에서 취하며, 이 점에서의 강체변위(rigid displacements)는 요소의 주변에 취해진 스프링을 통하여 다른 요소로 전달된다. 상기와 같은 강체요소법의 개념을 연속체의 탄성 및 탄소성해석에 적용하면, 해석적 개념이 단순할 뿐만 아니라 구조물 전체의 자유도수를 대폭 줄여 컴퓨터 계산시간을 절약할 수 있는 잇점이 있고, 거시적인 모델(macroscopic modeling)과 미시적인 모델 (microscopic modeling)의 중간적인 성격을 가지기 때문에 구조물의 파괴상황에 대해서도 그 개략을 파악할 수 있다. 본 논문에서는 강체요소법을 보다 일반화된 해석법으로 개발, 확장하기 위해서 종전에 단층스프링시스템(single-layer spring system)으로 해석이 어려웠던 문제점들을 보완한 복층프링시스템(double-layer spring system)을 사용함으로서 휨, 비틀림의 효과를 파악할 수 있는 이론적 개념을 적용한 새로운 구요소, 원통요소 및 평면요소를 개발하고, 이러한 강체요소들의 적합매트릭스의 유도 및 해석저긴 방법을 정식화하였다. 또 휨, 비틀림 및 전단력의 효과를 고려한 사각형원통요소 및 능형원 통요소를 이용하여 원통쉘의 탄성 및 탄소성해석할 수 있는 프로그램을 개발하고, 이 프로그램으로 캔틸레버로된 연속형철근콘크리트 원통쉘의 탄성 및 탄소성해석에 적용하여 구조물의 거동에 관한 수치해석의 결과, 즉 내력의 분포, 균열의 진전, 파괴의 상황 및 변형의 상태 등을 파악해 보았다.

• Computational Structural Engineering
• /
• v.10 no.3
• /
• pp.241-250
• /
• 1997

The strain localization of concrete is a phenomenon such that the deformation of concrete is localized in finite region along with softening behavior. The objective of this paper is to develop a plasticity and damage algorithm for the finite element analysis of the strain-localization in concrete. In this paper, concrete member under strain localization is modeled with localized zone and non-localized zone. For modeling of the localized zone in concrete under strain localization, a general Drucker-Prager failure criterion by which the nonlinear strain softening behavior of concrete after peak-stress can be considered is introduced in a thermodynamic formulation of the classical plasticity model. The return-mapping algorithm is used for the integration of the elasto-plastic rate equation and the consistent tangent modulus is also derived. For the modeling of non-localized zone in concrete under strain localization, a consistent nonlinear elastic-damage algorithm is developed by modifying the free energy in thermodynamics. Using finite element program implemented with the developed algorithm, strain localization behaviors for concrete specimens under compression are simulated.

• Journal of the Korean GEO-environmental Society
• /
• v.7 no.6
• /
• pp.133-143
• /
• 2006

A fully coupled effective stress dynamic analysis procedure for modeling seismic liquefaction on slope is presented. An elasto-plastic formulation is used for the constitutive model UBCSAND in which the yield loci are radial lines of constant stress ratio and the flow rule is non-associated. This is incorporated into the 2D version of Fast Lagrangian Analysis of Continua (FLAC) by modifying the existing Mohr-Coulomb model. This numerical procedure is used to simulate centrifuge test data from the Rensselaer Polytechnic Institute (RPI). UBCSAND is first calibrated to cyclic direct simple shear tests performed on Nevada sand. Both pre- and post-liquefaction behaviour is captured. The centrifuge test is then modeled and the predicted accelerations, excess porewater pressures, and displacements are compared with the measurements. The results are shown to be in general agreement. The procedure is currently being used in the design of liquefaction remediation measures for a number of dam, bridge, tunnel, and pipeline projects in Western Canada.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.35 no.9
• /
• pp.1041-1050
• /
• 2011

An elasto-plastic finite element method using the theory of strain gradient plasticity is proposed to evaluate the size dependency of structural plasticity that occurs when the configuration size decreases to micron scale. For this method, we suggest a low-order plane and three-dimensional displacement-based elements, eliminating the need for a high order, many degrees of freedom, a mixed element, or super elements, which have been considered necessary in previous researches. The proposed method can be performed in the framework of nonlinear incremental analysis in which plastic strains are calculated and averaged at nodes. These strains are then interpolated and differentiated for gradient calculation. We adopted a strain-gradient-hardening constitutive equation from the Taylor dislocation model, which requires the plastic strain gradient. The developed finite elements are tested numerically on the basis of typical size-effect problems such as micro-bending, micro-torsion, and micro-voids. With respect to the strain gradient plasticity, i.e., the size effects, the results obtained by using the proposed method, which are simple in their calculation, are in good agreement with the experimental results cited in previously published papers.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.10 no.4
• /
• pp.349-360
• /
• 2008

In this paper, the modeling procedure of random field with an elasto-plastic finite element algorithm and probability of failure on closely-spaced tunnel were investigated. Local average subdivision (LAS) method which can generate discrete random variables fast and accurately as well as change the resolution in certain region was used. And correlated value allocating and weighted average method were suggested to implement geometrical characteristics of tunnel. After the probability of failure on the test problem was thoroughly investigated using random finite element method, the results were compared with the deterministic strength reduction factor method and single random variable method. Of particular importance in this work, is the conclusion that the probability of failure determined by simplified probabilistic analysis, in which spatial variability is ignored by assuming perfect correlation, can be estimated from the safety factor determined by strength reduction factor method. Also, single random variable method can lead to unconservative estimates of the probability of failure.

• Geomechanics and Engineering
• /
• v.27 no.5
• /
• pp.465-479
• /
• 2021

One of the important causes of building and infrastructure failure, such as bridges on pile foundations, is the placement of the piles in liquefiable soil that can become unstable under seismic loads. Therefore, the overarching aim of this study is to investigate the seismic behavior of a soil-pile system in liquefiable soil using three-dimensional numerical FEM analysis, including soil-pile interaction. Effective parameters on concrete pile response, involving the pile diameter, pile length, soil type, and base acceleration, were considered in the framework of finite element non-linear dynamic analysis. The constitutive model of soil was considered as elasto-plastic kinematic-isotropic hardening. First, the finite element model was verified by comparing the variations on the pile response with the measured data from the centrifuge tests, and there was a strong agreement between the numerical and experimental results. Totally 64 non-linear time-history analyses were conducted, and the responses were investigated in terms of the lateral displacement of the pile, the effect of the base acceleration in the pile behavior, the bending moment distribution in the pile body, and the pore pressure. The numerical analysis results demonstrated that the relationship between the pile lateral displacement and the maximum base acceleration is non-linear. Furthermore, increasing the pile diameter results in an increase in the passive pressure of the soil. Also, piles with small and big diameters are subjected to yielding under bending and shear states, respectively. It is concluded that an effective stress-based ground response analysis should be conducted when there is a liquefaction condition in order to determine the maximum bending moment and shear force generated within the pile.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.10 no.3 s.49
• /
• pp.101-111
• /
• 2006

This paper proposes a method for modeling new fault ground motion due to moderate size earthquakes in Stable Continental Regions (SCRs) for the first time. The near fault ground motion is characterized by a single long period velocity pulse of large amplitude. In order to model the velocity pulse, its period and peak amplitude need be determined in terms of earthquake magnitude and distance from the causative fault. Because there have been observed very few new fault ground motions, it is difficult to derive the model directly from the recorded data in SCRs. Instead an indirect approach is adopted in this work. The two parameters, the period and peak amplitude of the velocity pulse, are known to be functions of the rise time and the slip velocity. For Western United States (WUS) that belongs active tectonic regions, there art empirical formulas for these functions. The relations of rise time and slip velocity on the magnitude in SCRs are derived by comparing related data between Western United States and Central-Eastern United States that belongs to SCRs. From these relations, the functions of these pulse parameters for NFGM in SCRs can be expressed in terms of earthquake magnitude and distance. A time history of near fault ground motion of moderate magnitude earthquake in stable continental regions is synthesized by superposing the velocity pulse on the for field ground motion that is generated by stochastic method. As an demonstrative application, the response of a single degree of freedom elasto-plastic system is studied.