• Interaction and multiscale mechanics
• /
• 제1권1호
• /
• pp.83-101
• /
• 2008

A derivation is given of two-scale models that are able to describe deformation and flow in a fluid-saturated and progressively fracturing porous medium. From the micromechanics of the flow in the cavity, identities are derived that couple the local momentum and the mass balances to the governing equations for a fluid-saturated porous medium, which are assumed to hold on the macroscopic scale. By exploiting the partition-of-unity property of the finite element shape functions, the position and direction of the fractures are independent from the underlying discretization. The finite element equations are derived for this two-scale approach and integrated over time. The resulting discrete equations are nonlinear due to the cohesive crack model and the nonlinearity of the coupling terms. A consistent linearization is given for use within a Newton-Raphson iterative procedure. Finally, examples are given to show the versatility and the efficiency of the approach.

• 토지주택연구
• /
• 제1권1호
• /
• pp.59-65
• /
• 2010

The paper presents a comparison between a semi-implicit time integration linear finite element implementation and fully-implicit nonlinear Newton-Raphson finite element implementation of a triphasic small strain mixture formulation of an elastic partially saturated porous medium. The pore air phase pressure pa is assumed atmospheric, i.e., $p_a$ = 0, although the formulation and implementation are general to handle increase in pore air pressure as a result of loading, if needed. The solid skeleton phase is assumed linear isotropic elastic and partially saturated 'consolidation' in the presence of surface infiltration and traction is simulated. The verification of the implementation against an analytical solution for partially saturated pore water flow (no deformation) and comparison between the two implementations is presented and the important of the porosity-dependent nature of the partially saturated permeability is assessed on comparison with a commercial code for the partially saturated flow with deformation. As a result, the response of partially saturated permeability subjected to the porosity influences on the saturation of a soil, and the different behaviors of the partially saturated soil between staggered and monolithic coupled programs is worth of attention because the negative pore water pressure in the partially saturated soil depends on the difference.

• Korean Chemical Engineering Research
• /
• 제53권6호
• /
• pp.814-817
• /
• 2015

A new ionic mass transfer correlation is derived for the fluid-saturated, horizontal porous layer. Darcy-Forchheimer model is used to explain characteristics of fluid motion. Based on the microscales of turbulence a backbone mass transfer relation is derived as a function of the Darcy-Rayleigh number, $Ra_D$ and the porous medium Schmidt number, $Sc_p$. For the Darcy's limit of $Sc_p{\gg}Ra_D$, the Sherwood number, Sh is a function of $Ra_D$ only. However, for the region of high $Ra_D$, Sh can be related with $Ra_DSc_p$. Based on the present backbone equation and the electrochemical mass transfer experiments which are electro plating or electroless plating, the new ionic mass transfer correlation is suggested in the porous media.

• 한국지진공학회:학술대회논문집
• /
• 한국지진공학회 2005년도 학술발표회 논문집
• /
• pp.47-51
• /
• 2005

If a structure is founded on the ground saturated with pore water, then the ground should be modeled as a saturated two-phase porous medium for accurate earthquake response analysis. In this study, an axisymmetric transmitting boundary hyperelement is developed for modeling of far field of the ground using u-U formulation for water-saturated transversely isotropic layered stratum. The developed hyperelement is verified by comparing the dynamic stiffness of rigid circular foundation on water-saturated isotropic layered stratum with the case of using equivalent single-phase medium model.

• 한국지진공학회:학술대회논문집
• /
• 한국지진공학회 2006년도 학술발표회 논문집
• /
• pp.345-350
• /
• 2006

If a structure is founded on the ground saturated with pore water, then the ground should be modeled as a saturated two-phase porous medium for accurate earthquake response analysis. In this study, a 3-dimensional transmitting boundary is developed for modeling of far field using u-U formulation for water-saturated transversely isotropic layered stratum. The developed transmitting boundary is verified by comparing the dynamic stiffness of rigid square foundation on water-saturated isotropic layered stratum with the case of using equivalent single-phase medium model.

• Korean Chemical Engineering Research
• /
• 제57권5호
• /
• pp.723-729
• /
• 2019

A theoretical analysis was conducted of convective instability driven by buoyancy forces under transient temperature fields in an annular porous medium bounded by coaxial vertical cylinders. Darcy's law and Boussinesq approximation are used to explain the characteristics of fluid motion and linear stability theory is employed to predict the onset of buoyancy-driven motion. The linear stability equations are derived in a global domain, and then cast into in a self-similar domain. Using a spectral expansion method, the stability equations are reformed as a system of ordinary differential equations and solved analytically and numerically. The critical Darcy-Rayleigh number is founded as a function of the radius ratio. Also, the onset time and corresponding wavelength are obtained for the various cases. The critical time becomes smaller with increasing the Darcy-Rayleigh number and follows the asymptotic relation derived in the infinite horizontal porous layer.

• 한국지반공학회지:지반
• /
• 제13권1호
• /
• pp.5-18
• /
• 1997

일반적으로 포화지반의 동적거동을 정확히 예편하기 위해서는 다공성(multi-phase) 재료모델과 그 모델을 이용하는 수치해석 프로그램의 개발이 필요하다. 본 논문에서는 높은 동하중을 받는 다공질 재료의 이론적인 거동해석 연구결과와 함께 기존 MPDAP(multi-phase dynamic analysis program)에 JWL(Jones-Withins-Lee) 모델을 삽입시켜 개발한 MPDAP에 대해 다루었다. JWL모델은 기존 모델과는 달리 폭약의 종류 및 특성 등을 고려할 수 있는 모델이다. 또한 본 논문에서는 개발된 프로그램의 적합성을 조사하기 위하여 몇몇 예제에 대한 검증해석을 수행하였다. 검증결과, 단일매체 (single-phase medium)에서의 탄성구형파의 전파특성 해석의 경우 해석결과와 이론해는 거의 일치하는 결과를 나타내었고, 일차원 선형 압밀해석의 경우과잉 간극수압은 Terzaghi의 이론해와 해석된 결과가 비교적 일치하는 경향성을 보여 주었다. 또한 포화지반에서의 평면 압축파 해석의 경우도 해석결과와 이미 검증된 프로그램 또는 완전해의 해석결과는 거의 유사하게 나타나는 것을 알 수 있었다.

• Geomechanics and Engineering
• /
• 제15권2호
• /
• pp.793-803
• /
• 2018

Based on the theory of porous media, an interaction system of a floating pile and a saturated soil in cylindrical coordinates subjected to vertical harmonic load is presented in this paper. The surrounding soil is separated into two distinct layers. The upper soil layer above the level of pile base is described as a saturated viscoelastic medium and the lower soil layer is idealized as equivalent spring-dashpot elements with complex stiffness. Considering the cylindrically symmetry and the pile-soil compatibility condition of the interaction system, a frequency-domain analytical solution for dynamic impedance of the floating pile embedded in saturated viscoelastic soil is also derived, and reduced to verify it with existing solutions. An extensive parametric analysis has been conducted to reveal the effects of the impedance of the lower soil base, the interaction coefficient and the damping coefficient of the saturated viscoelastic soil layer on the vertical vibration of the pile-soil interaction system. It is shown that the vertical dynamic impedance of the floating pile significantly depends on the real stiffness of the impedance of the lower soil base, but is less sensitive to its dynamic damping variation; the behavior of the pile in poro-visco-elastic soils is totally different with that in single-phase elastic soils due to the existence of pore liquid; the effect of the interaction coefficient of solid and liquid on the pile-soil system is limited.

• Coupled systems mechanics
• /
• 제9권2호
• /
• pp.125-145
• /
• 2020

In this paper, we present a 3D thermo-hydro-mechanical coupled discrete beam lattice model of structure built of the nonisothermal saturated poro-plastic medium subjected to mechanical loads and nonstationary heat transfer conditions. The proposed model is based on Voronoi cell representation of the domain with cohesive links represented as inelastic Timoshenko beam finite elements enhanced with additional kinematics in terms of embedded strong discontinuities in axial and both transverse directions. The enhanced Timoshenko beam finite element is capable of modeling crack formation in mode I, mode II and mode III. Mode I relates to crack opening, mode II relates to in-plane crack sliding, and mode III relates to the out-of-plane shear sliding. The pore fluid flow and heat flow in the proposed model are governed by Darcy's law and Fourier's law for heat conduction, respectively. The pore pressure field and temperature field are approximated with linear tetrahedral finite elements. By exploiting nodal point quadrature rule for numerical integration on tetrahedral finite elements and duality property between Voronoi diagram and Delaunay tetrahedralization, the numerical implementation of the coupling results with additional pore pressure and temperature degrees of freedom placed at each node of a Timoshenko beam finite element. The results of several numerical simulations are presented and discussed.

• 한국지진공학회논문집
• /
• 제9권4호
• /
• pp.67-74
• /
• 2005

이 논문에서는 포화된 2상 지반의 동적해석에서 원역을 모형화하기 위한 새로운 무한요소를 제안하였다. 무한요소법은 무한영역 또는 반무한영역을 모형화해야 하는 공학문제에 효과적으로 적용되어 왔다. 그러나 현재까지 개발된 2상지반의 동적해석을 위한 무한요소는 형상함수에 사용될 수 있는 파동성문이 2개(Pl파와 P2파)로 한정되어 있다. 이 논문에서는 이와 같은 제한을 없애고 임의 개수의 파동성분을 고려할 수 있도록 하는 정식화 과정을 제안하였다. 구조물을 포함하는 근역은 유한요소로 나타내며 원역은 평행층상 반무한 지반으로 가정하였다. 제안된 무한요소의 타당성은 1차원 및 2차원 파동전달문제를 해석하고 이를 이론해 및 정밀수지해석 해와 비교하여 검증하였다.