• Title/Summary/Keyword: Finite-difference

Search Result 3,280, Processing Time 0.028 seconds

A Numerical Study on the Reduction Effect of Blasting Vibration with Cut Method (심발공법에 따른 발파진동 저감효과에 대한 수치해석적 연구)

  • Son, Ji-Ho;Kim, Byung-Ryeol;Lee, Seung-Joong;Kim, Nam-Soo;Lee, Hyo;Choi, Sung-Oong
    • Explosives and Blasting
    • /
    • v.37 no.1
    • /
    • pp.1-13
    • /
    • 2019
  • The repeated blasting vibration, which is induced commonly in NATM excavation site, can cause a severe damage to the nearby facilities. It is known that the most effective method for reducing blasting vibration includes the use of electronic detonator, deck charge and change of cut method, and so forth. In order to analyze the effect of blasting vibration reduction, in this study, three-dimensional FDM (Finite Difference Method) program FLAC3D has been used for reflecting the blasting hole, delayed time and charging amount. Also the numerical analysis has been performed by applying a dynamic load to each blasting hole. The cut method has been applied with several methods, such as V-cut and Double-drilled parallel cut, which are common in tunnel construction sites. Also, the field test blasting has been carried out in order to compare the measured data with results of numerical analysis. It was shown that the numerical analysis and the field measurement coincide well.

Effect of the lateral earth pressure coefficient on settlements during mechanized tunneling

  • Golpasand, Mohammad-Reza B.;Do, Ngoc Anh;Dias, Daniel;Nikudel, Mohammad-Reza
    • Geomechanics and Engineering
    • /
    • v.16 no.6
    • /
    • pp.643-654
    • /
    • 2018
  • Tunnel excavation leads to a disturbance on the initial stress balance of surrounding soils, which causes convergences around the tunnel and settlements at the ground surface. Considering the effective impact of settlements on the structures at the surface, it is necessary to estimate them, especially in urban areas. In the present study, ground settlements due to the excavation of East-West Line 7 of the Tehran Metro (EWL7) and the Abuzar tunnels are evaluated and the effect of the lateral earth pressure coefficient ($K_0$) on their extension is investigated. The excavation of the tunnels was performed by TBMs (Tunnel Boring Machines). The coefficient of lateral earth pressure ($K_0$) is one of the most important geotechnical parameters for tunnel design and is greatly influenced by the geological characteristics of the surrounding soil mass along the tunnel route. The real (in-situ) settlements of the ground surface were measured experimentally using leveling methods along the studied tunnels and the results were compared with evaluated settlements obtained from both semi-empirical and numerical methods (using the finite difference software FLAC3D). The comparisons permitted to show that the adopted numerical models can effectively be used to predict settlements induced by a tunnel excavation. Then a numerical parametric study was conducted to show the influence of the $K_0$ values on the ground settlements. Numerical investigations also showed that the shapes of settlement trough of the studied tunnels, in a transverse section, are not similar because of their different diameters and depths of the tunnels.

Temperature distribution of ceramic panels of a V94.2 gas turbine combustor under realistic operation conditions

  • Namayandeh, Mohammad Javad;Mohammadimehr, Mehdi;Mehrabi, Mojtaba
    • Advances in materials Research
    • /
    • v.8 no.2
    • /
    • pp.117-135
    • /
    • 2019
  • The lifetime of a gas turbine combustor is typically limited by the durability of its liner, the structure that encloses the high-temperature combustion products. The primary objective of the combustor thermal design process is to ensure that the liner temperatures do not exceed a maximum value set by material limits. Liner temperatures exceeding these limits hasten the onset of cracking which increase the frequency of unscheduled engine removals and cause the maintenance and repair costs of the engine to increase. Hot gas temperature prediction can be considered a preliminary step for combustor liner temperature prediction which can make a suitable view of combustion chamber conditions. In this study, the temperature distribution of ceramic panels for a V94.2 gas turbine combustor subjected to realistic operation conditions is presented using three-dimensional finite difference method. A simplified model of alumina ceramic is used to obtain the temperature distribution. The external thermal loads consist of convection and radiation heat transfers are considered that these loads are applied to flat segmented panel on hot side and forced convection cooling on the other side. First the temperatures of hot and cold sides of ceramic are calculated. Then, the thermal boundary conditions of all other ceramic sides are estimated by the field observations. Finally, the temperature distributions of ceramic panels for a V94.2 gas turbine combustor are computed by MATLAB software. The results show that the gas emissivity for diffusion mode is more than premix therefore the radiation heat flux and temperature will be more. The results of this work are validated by ANSYS and ABAQUS softwares. It is showed that there is a good agreement between all results.

Lateral torsional buckling of steel I-beams: Effect of initial geometric imperfection

  • Bas, Selcuk
    • Steel and Composite Structures
    • /
    • v.30 no.5
    • /
    • pp.483-492
    • /
    • 2019
  • In the current study, the influence of the initial lateral (sweep) shape and the cross-sectional twist imperfection on the lateral torsional buckling (LTB) response of doubly-symmetric steel I-beams was investigated. The material imperfection (residual stress) was not considered. For this objective, standard European IPN 300 beam with different unbraced span was numerically analyzed for three imperfection cases: (i) no sweep and no twist (perfect); (ii) three different shapes of global sweep (half-sine, full-sine and full-parabola between the end supports); and (iii) the combination of three different sweeps with initial sinusoidal twist along the beam. The first comparison was done between the results of numerical analyses (FEM) and both a theoretical solution and the code lateral torsional buckling formulations (EC3 and AISC-LRFD). These results with no imperfection effects were then separately compared with three different shapes of global sweep and the presence of initial twist in these sweep shapes. Besides, the effects of the shapes of initial global sweep and the inclusion of sinusoidal twist on the critical buckling load of the beams were investigated to unveil which parameter was considerably effective on LTB response. The most compatible outcomes for the perfect beams was obtained from the AISC-LRFD formulation; however, the EC-3 formulation estimated the $P_{cr}$ load conservatively. The high difference from the EC-3 formulation was predicted to directly originate from the initial imperfection reduction factor and high safety factor in its formulation. Due to no consideration of geometric imperfection in the AISC-LFRD code solution and the theoretical formulation, the need to develop a practical imperfection reduction factor for AISC-LRFD and theoretical formulation was underlined. Initial imperfections were obtained to be more influential on the buckling load, as the unbraced length of a beam approached to the elastic limit unbraced length ($L_r$). Mode-compatible initial imperfection shapes should be taken into account in the design and analysis stages of the I-beam to properly estimate the geometric imperfection influence on the $P_{cr}$ load. Sweep and sweep-twist imperfections led to 10% and 15% decrease in the $P_{cr}$ load, respectively, thus; well-estimated sweep and twist imperfections should considered in the LTB of doubly-symmetric steel I-beams.

Physical and numerical modelling of the inherent variability of shear strength in soil mechanics

  • Chenari, Reza Jamshidi;Fatahi, Behzad;Ghoreishi, Malahat;Taleb, Ali
    • Geomechanics and Engineering
    • /
    • v.17 no.1
    • /
    • pp.31-45
    • /
    • 2019
  • In this study the spatial variability of soils is substantiated physically and numerically by using random field theory. Heterogeneous samples are fabricated by combining nine homogeneous soil clusters that are assumed to be elements of an adopted random field. Homogeneous soils are prepared by mixing different percentages of kaolin and bentonite at water contents equivalent to their respective liquid limits. Comprehensive characteristic laboratory tests were carried out before embarking on direct shear experiments to deduce the basic correlations and properties of nine homogeneous soil clusters that serve to reconstitute the heterogeneous samples. The tests consist of Atterberg limits, and Oedometric and unconfined compression tests. The undrained shear strength of nine soil clusters were measured by the unconfined compression test data, and then correlations were made between the water content and the strength and stiffness of soil samples with different consistency limits. The direct shear strength of heterogeneous samples of different stochastic properties was then evaluated by physical and numerical modelling using FISH code programming in finite difference software of $FLAC^{3D}$. The results of the experimental and stochastic numerical analyses were then compared. The deviation of numerical simulations from direct shear load-displacement profiles taken from different sources were discussed, potential sources of error was introduced and elaborated. This study was primarily to explain the mathematical and physical procedures of sample preparation in stochastic soil mechanics. It can be extended to different problems and applications in geotechnical engineering discipline to take in to account the variability of strength and deformation parameters.

Numerical modelling of bottom-hole rock in underbalanced drilling using thermo-poroelastoplasticity model

  • Liu, Weiji;Zhou, Yunlai;Zhu, Xiaohua;Meng, Xiannan;Liu, Mei;Wahab, Magd Abdel
    • Structural Engineering and Mechanics
    • /
    • v.69 no.5
    • /
    • pp.537-545
    • /
    • 2019
  • Stress analysis of bottom-hole rock has to be considered with much care to further understand rock fragmentation mechanism and high penetration rate. This original study establishes a fully coupled simulation model and explores the effects of overburden pressure, horizontal in-situ stresses, drilling mud pressure, pore pressure and temperature on the stress distribution in bottom-hole rock. The research finds that in air drilling, as the well depth increases, the more easily the bottom-hole rock is to be broken. Moreover, the mud pressure has a great effect on the bottom-hole rock. The bigger the mud pressure is, the more difficult to break the bottom-hole rock is. Furthermore, the maximum principal stress of the bottom-hole increases as the mud pressure, well depth and temperature difference increase. The bottom-hole rock can be divided into three main regions according to the stress state, namely a) three directions tensile area, b) two directions compression areas and c) three directions compression area, which are classified as a) easy, b) normal and c) hard, respectively, for the corresponding fragmentation degree of difficulty. The main contribution of this paper is that it presents for the first time a thorough study of the effect of related factors, including stress distribution and temperature, on the bottom-hole rock fracture rather than the well wall, using a thermo-poroelastoplasticity model.

Numerical Study on Package Warpage as Structure Modeling Method of Materials for a PCB of Semiconductor Package (반도체 패키지용 PCB의 구조 모델링 방법에 따른 패키지의 warpage 수치적 연구)

  • Cho, Seunghyun;Ceon, Hyunchan
    • Journal of the Microelectronics and Packaging Society
    • /
    • v.25 no.4
    • /
    • pp.59-66
    • /
    • 2018
  • In this paper, we analyzed the usefulness of single-structured printed circuit board (PCB) modeling by using numerical analysis to model the PCB structure applied to a package for semiconductor purposes and applying modeling assuming a single structure. PCBs with circuit layer of 3rd and 4th were used for analysis. In addition, measurements were made on actual products to obtain material characteristics of a single structure PCB. The analysis results showed that if the PCB was modeled in a single structure compared to a multi-layered structure, the warpage analysis results resulting from modeling the PCB structure would increase and there would be a significant difference. In addition, as the circuit layer of the PCB increased, the mechanical properties of the PCB, the elastic coefficient and inertia moment of the PCB increased, decreasing the package's warpage.

Analysis of the Deformed Unit Cell by Clamping Force Through the FEM and CFD Interaction (FEM과 CFD 연동을 통한 스택 체결 시 압력에 의해 변형된 단위 전지 해석)

  • YOO, BIN;LIM, KISUNG;JU, HYUNCHUL
    • Journal of Hydrogen and New Energy
    • /
    • v.32 no.4
    • /
    • pp.228-235
    • /
    • 2021
  • Polymer electrolyte membrane fuel cells (PEMFC) are currently being used in various transport applications such as drones, unmanned aerial vehicles, and automobiles. The power required is different according to the type of use, purpose, and the conditions adjusted using a cell stack. The fuel cell stack is compressed to reduce the size and prevent fuel leakage. The unit cells that make up the cell stack are subjected to compression by clamping force, which makes geometrical changes in the porous media and it impacts on cell performance. In this study, finite elements method (FEM) and computational fluid dynamics (CFD) analysis for the deformed unit cell considering the effects of clamping force is performed. First, structural analysis using the FEM technique over the deformed gas diffusion layer (GDL) considering compression is carried out, and the resulting porosity changed in the GDL is calculated. The PEMFC model is then verified by a three-dimensional, two-phase fuel cell simulation applying the physical properties and geometry obtained before and after compression. The detailed simulation results showed different concentration distributions of fuel between the original and deformed geometry, resulting in the difference in the distribution of current density is represented at compressed GDL region with low oxygen concentration.

The effect of curvature at the bottom of a soft ground tunnel by numerical analysis (수치해석에 의한 연약지반 터널의 바닥부 곡률의 영향 분석)

  • You, Kwangho;Kim, Kangsan
    • Journal of Korean Tunnelling and Underground Space Association
    • /
    • v.23 no.2
    • /
    • pp.107-118
    • /
    • 2021
  • Due to the acceleration of road construction, the number and extension of tunnels are increasing every year. A lot of research has been done on the collapse of tunnels, but research on the invert heaving is insufficient. Therefore, in this study, a sensitivity analysis was performed using a geotechnical general-purpose program to analyze the effect of the invert curvature of a tunnel excavated on the soft ground. As a result, it was quantitatively confirmed that the stability of a tunnel was increased as the curvature of the tunnel invert was increased so that the safety factor was calculated to be large regardless of the ground conditions and the thickness of the support. In addition, it was confirmed that the stability of the tunnel was increased by reducing the convergence of the tunnel and the maximum bending stress supported by shotcrete. Therefore, when a tunnel is excavated on soft ground, it is believed that applying a curvature to the invert will increase the stability of the tunnel.

Development of 2D Urban Inundation Analysis Model using Adaptive Mesh Refinement Method (메쉬 세분화 기법을 이용한 2차원 침수해석 모형의 개발)

  • Lee, Seung-soo
    • Proceedings of the Korea Water Resources Association Conference
    • /
    • 2016.05a
    • /
    • pp.93-93
    • /
    • 2016
  • 최근 증가하고 있는 기후변화에 의해 설계빈도를 상회하는 강우의 발생빈도가 증가하고 있으며, 이로 인한 도시유역의 내수범람 피해가 증가하고 있다. 도시유역에서 발생하는 침수 피해의 경우 인적 물적 자원이 집중되어 있는 도시의 특성으로 인해 침수로 인한 직접적 피해 규모가 상당할 뿐만 아니라 침수 발생 후 세균 및 박테리아에 의해 발생하는 수인성 전염병의 유행 등과 같은 2차적 피해 또한 심각한 사회적 비용을 초래할 수 있어 도시유역의 침수 피해를 저감시키기 위한 대책이 절실히 요구되어지고 있다. 도시유역의 침수를 예방하기 위한 대책은 구조적 비구조적 대책으로 구분되어 질 수 있으며 구조적 대책의 경우 침수 피해 예방에 직접적인 효과를 낼 수 있다는 장점이 있으나 대규모 사업예산 및 사업 기간으로 인해 직접적 효과를 보기까지 상대적으로 긴 시간이 필요할 뿐만 아니라 사업 진행 중 대상지역 거주민들의 민원으로 인한 갈등 조정 등으로 인해 사업실행에 어려움을 겪고 있다. 이러한 측면에서 비구조적 대책의 일환인 수치해석을 통한 침수피해 재현 및 침수원인 파악을 통한 구조개선 제안은 구조적 대책의 단점을 보완할 수있는 좋은 대안이 될 수 있다. 도시유역의 경우 비도시유역과 대조적인 차이점으로는 높은 비율의 불투수층, 복잡한 지형, 다수의 인공 구조물 및 배수관망 시스템 등을 들 수 있으며, 침수해석 모형의 정확도를 높이기 위해서는 복잡한 지형의 효율적인 처리가 무엇보다 중요하다. 일반적으로 이용되는 2차원 침수해석 모형들은 직교구조 격자 또는 비구조 격자를 이용하여 지형을 묘사하고 있으며 DEM 자료를 직접 사용하는 직교구조 격자의 경우 지형 데이터 생성이 상대적으로 쉽다는 장점이 있으나 복잡한 지형을 표현하기 위해서는 불필요한 지역까지 높은 해상도를 이용해야 하며 이로 인하여 모의시간이 지나치게 길어지는 문제점이 발생한다. 비구조 격자의 경우 상대적으로 복잡한 도시 유역을 잘 묘사할 수 있다는 장점이 있으나 격자망 생성에 필요한 데이터가 많고 격자망 생성에 지나치게 많은 시간과 노력이 소요된다는 단점이 있다. 따라서 본 연구에서는 위에서 언급한 두 가지 방법의 장점만을 취할 수 있도록 메쉬 세분화 기법을 이용한 2차원 침수해석 모형을 개발 하여 복잡한 지형은 고해상도 메쉬를 이용하여 보다 자세히 묘사하고 상대적으로 복잡하지 않은 지형은 저해상도 메쉬를 이용하여 계산시간을 단축시킬 수 있도록 하였다. 수치해석 기법으로는 엇갈림 격자를 이용하는 Leap-Frog 기법과 유한차분 (Finite difference Method)기법을 이용하였다.

  • PDF