• Title/Summary/Keyword: Finite-difference method

Search Result 2,416, Processing Time 0.037 seconds

Analysis on dynamic numerical model of subsea railway tunnel considering various ground and seismic conditions (다양한 지반 및 지진하중 조건을 고려한 해저철도 터널의 동적 수치모델 분석)

  • Changwon Kwak;Jeongjun Park;Mintaek Yoo
    • Journal of Korean Tunnelling and Underground Space Association
    • /
    • v.25 no.6
    • /
    • pp.583-603
    • /
    • 2023
  • Recently, the advancement of mechanical tunnel boring machine (TBM) technology and the characteristics of subsea railway tunnels subjected to hydrostatic pressure have led to the widespread application of shield TBM methods in the design and construction of subsea railway tunnels. Subsea railway tunnels are exposed in a constant pore water pressure and are influenced by the amplification of seismic waves during earthquake. In particular, seismic loads acting on subsea railway tunnels under various ground conditions such as soft ground, soft soil-rock composite ground, and fractured zones can cause significant changes in tunnel displacement and stress, thereby affecting tunnel safety. Additionally, the dynamic response of the ground and tunnel varies based on seismic load parameters such as frequency characteristics, seismic waveform, and peak acceleration, adding complexity to the behavior of the ground-tunnel structure system. In this study, a finite difference method is employed to model the entire ground-tunnel structure system, considering hydrostatic pressure, for the investigation of dynamic behavior of subsea railway tunnel during earthquake. Since the key factors influencing the dynamic behavior during seismic events are ground conditions and seismic waves, six analysis cases are established based on virtual ground conditions: Case-1 with weathered soil, Case-2 with hard rock, Case-3 with a composite ground of soil and hard rock in the tunnel longitudinal direction, Case-4 with the tunnel passing through a narrow fault zone, Case-5 with a composite ground of soft soil and hard rock in the tunnel longitudinal direction, and Case-6 with the tunnel passing through a wide fractured zone. As a result, horizontal displacements due to earthquakes tend to increase with an increase in ground stiffness, however, the displacements tend to be restrained due to the confining effects of the ground and the rigid shield segments. On the contrary, peak compressive stress of segment significantly increases with weaker ground stiffness and the effects of displacement restrain contribute the increase of peak compressive stress of segment.

Numerical study for Application of H-Pile Connection Plastic Sheet Pile Retaining Wall (HCS) (H-Pile과 Plastic Sheet Pile을 결합한 토류벽체에 대한 수치해석적 연구)

  • Lee, Kyou-Nam;Lim, Hee-Dae
    • The Journal of Engineering Geology
    • /
    • v.27 no.3
    • /
    • pp.331-343
    • /
    • 2017
  • In this study to improve stability, workability and economics of the H-Pile+Earth plate or H-Pile+Earth plate+Cutoff grouting currently in use, we had developed HCS method belonging to the retaining wall which is consisting of a combination H-Pile, Plastic Sheet Pile and Steel Square Pipe for gap maintenance and reinforcement of flexible plastic Sheet Pile, and the behavior of each member composing HCS method is investigated by three-dimensional finite element analysis. To numerically analyze the behavior of the HCS method, we have performed extensive three-dimentional finite element analysis for three kinds of plastic Sheet Pile size, two kinds of H-Pile size and three kinds of H-Pile installation interval, one kinds of Steel Square Pipe and three kinds of Steel Square Pipe installation interval. After analyzing the numerical results, we found that the combinations of $P.S.P-460{\times}131.5{\times}7t$ (PS7) and H-Pile $250{\times}250{\times}9{\times}14$ (H250), $P.S.P473{\times}133.5{\times}9t$ (PS9) and H-Pile $300{\times}200{\times}9{\times}14$ (H300) is the most economical because these combinations are considered to have a stress ratio (=applied stress/allowable stress) close to that as the stiffness of H-Pile, plastic Sheet Pile and Steel Square Pipe composite increased, the horizontal displacement of the retaining wall and the vertical displacement of the upper ground decreased. Especially, due to the arching effects caused by the difference in stiffness between H-Pile and plastic Sheet Pile, a large part of the earth pressure acting on plastic Sheet Pile caused a stress transfer to H-Pile, and the stress and displacement of plastic Sheet Pile were small. Through this study, we can confirm the behavior of each member constituting the HCS method, and based on the confirmed results of this study, it can be used to apply HCS method in reasonable, stable and economical way in the future.

Stress analysis of Multiloop Edgewise Arch Wire with various degree of tip back bend : a study using the finite element method (Multiloop Edgewise Arch Wire의 tip back 정도에 따른 응력 분포에 관한 유한요소법적 연구)

  • Lee, Young-Il;Cha, Kyung-Suk;Ju, Jin-Won;Lee, Jin-Woo
    • The korean journal of orthodontics
    • /
    • v.30 no.2 s.79
    • /
    • pp.127-142
    • /
    • 2000
  • This study have been carried out to find out the mechnical effect of Multiloop Edgewise Arch Wire(MEAW) making use of the finite element method. The tip back bend of MEAW taken in this analysis is $5^{\circ},\;10{\circ}\;and\;15{\circ}$. In addition, Class II or up & down elastic is applied to find out stress distribution and their values in PDL. A adult male of normal occlusion was selected to create the models of teeth and PDL. And the model of MEAW was also created using commercial finite element code (ANSYS version 5.2). The MEAW is forcibly engaged with a class II or up & down elastic, to determine the initial stress generated in PDL. Comparing the compressive and tensile stress at each reference-planes, following results are obtained. 1. When a MEAW of $5^{\circ},\;10{\circ}\;15{\circ}$ tip back bend was engaged with Class II or up & down elastic, the distribution of compressive, tensile stress in entire PDL is similar in each case. 2. The values of compressive and tensile stress in PDL is higher in $15{\circ}$ tip back bend case than in $10{\circ}\;or\;15{\circ}$ tip back bend case. 3. In the distal PDL of 1st and 2nd molar, compressive stress appears. The compressive area is more wide and its values is higher in PDL of 2nd molar than those in 1st molar. The compressive area and its values become more wide and higher according to the increase of the tip back bend. 4. The values of compressive stress are comparatively smaIIer in PDL of molars than those in premolars. 5. Comparing class II and up & down elastic case, tensile stress values in anterior teeth PDL are smaller md their distribution is more wide in up & down elastic case than class If elastic case. On another hand, there is no difference in distribution and stress values in PDL of posterior teeth between two cases. 6. Comparing the tensile area in PDL of anterior teeth, tensile stress values are maximum in PDL of canine.

  • PDF

A numerical simulation of propagating turbidity currents using the ULTIMATE scheme (ULTIMATE 기법을 이용한 부유사 밀도류 전파 수치모의)

  • Choi, Seongwook;Choi, Sung-Uk
    • Journal of Korea Water Resources Association
    • /
    • v.50 no.1
    • /
    • pp.55-64
    • /
    • 2017
  • This study presents a numerical model for simulating turbidity currents using the ULTIMATE scheme. For this, the layer-averaged model is used. The model is applied to laboratory experiments, where the flume is composed of sloping and flat parts, and the characteristics of propagating turbidity currents are investigated. Due to the universal limiter of the ULTIMATE scheme, the frontal part of the turbidity currents at a sharp gradient without numerical oscillations is computed. Simulated turbidity currents propagate super-critically to the end of the flume, and internal hydraulic jumps occur at the break-in-slope after being affected by the downstream boundary. It is found that the hydraulic jumps are computed without numerical oscillations if Courant number is less than 1. In addition, factors that affect propagation velocity of turbidity currents is studied. The particle size less than $9{\mu}m$ does not affect propagation velocity but the buoyancy flux affects clearly. Finally, it is found that the numerical model computes the bed elevation change due to turbidity currents properly. Specifically, a discontinuity in the bed elevation, arisen from the hydraulic jumps and resulting difference in sediment entrainment, is observed.

A Study on Heat Transfer Analysis around the Square Heat Source of Interior Solid by Using Finite Element Method (유한요소법에 의한 고체내부의 사각열원 주위 열전도 특성연구)

  • Jang, Jae-Eun;Hong, Bong-Gi
    • Journal of the Korean Society of Fisheries and Ocean Technology
    • /
    • v.18 no.2
    • /
    • pp.101-108
    • /
    • 1982
  • In this paper the result of heat transfer analysis around the square heat source of interior solid by using the F. E. M is reported. Calculation for temperature distribution and each element was used by F. E. M. the solid is sub-divided into system of equal size triangular shape. These values of temperature distribution will valuable for design of jet engine and steam generator and the results gained are as follow; 1. Calculation by F. E. M is identified with the experiment. 2. Temperature distribution on the horizontal surface is $\theta$=0.698 in model 4 and the other hand $\theta$=0.401 in model 6 for X=16cm. In intermediates surface between heat source and bottom surface, the influence of L is more greater than that of height in the temperature difference. 3. Temperature distribution on the vertical surface for model 2 is resulted strong influence by K. In the case of Y=4cm is identified with $\theta$=0.0790 for K=7 and also $\theta$=0.0036 for K=0.3. In the difference of temperature distribution.

  • PDF

Analysis of Subwavelength Metal Hole Array Structure for the Enhancement of Quantum Dot Infrared Photodetectors

  • Ha, Jae-Du;Hwang, Jeong-U;Gang, Sang-U;No, Sam-Gyu;Lee, Sang-Jun;Kim, Jong-Su;Krishna, Sanjay;Urbas, Augustine;Ku, Zahyun
    • Proceedings of the Korean Vacuum Society Conference
    • /
    • 2013.02a
    • /
    • pp.334-334
    • /
    • 2013
  • In the past decade, the infrared detectors based on intersubband transition in quantum dots (QDs) have attracted much attention due to lower dark currents and increased lifetimes, which are in turn due a three-dimensional confinement and a reduction of scattering, respectively. In parallel, focal plane array development for infrared imaging has proceeded from the first to third generations (linear arrays, 2D arrays for staring systems, and large format with enhanced capabilities, respectively). For a step further towards the next generation of FPAs, it is envisioned that a two-dimensional metal hole array (2D-MHA) structures will improve the FPA structure by enhancing the coupling to photodetectors via local field engineering, and will enable wavelength filtering. In regard to the improved performance at certain wavelengths, it is worth pointing out the structural difference between previous 2D-MHA integrated front-illuminated single pixel devices and back-illuminated devices. Apart from the pixel linear dimension, it is a distinct difference that there is a metal cladding (composed of a number of metals for ohmic contact and the read-out integrated circuit hybridization) in the FPA between the heavily doped gallium arsenide used as the contact layer and the ROIC; on the contrary, the front-illuminated single pixel device consists of two heavily doped contact layers separated by the QD-absorber on a semi-infinite GaAs substrate. This paper is focused on analyzing the impact of a two dimensional metal hole array structure integrated to the back-illuminated quantum dots-in-a-well (DWELL) infrared photodetectors. The metal hole array consisting of subwavelength-circular holes penetrating gold layer (2DAu-CHA) provides the enhanced responsivity of DWELL infrared photodetector at certain wavelengths. The performance of 2D-Au-CHA is investigated by calculating the absorption of active layer in the DWELL structure using a finite integration technique. Simulation results show the enhanced electric fields (thereby increasing the absorption in the active layer) resulting from a surface plasmon, a guided mode, and Fabry-Perot resonances. Simulation method accomplished in this paper provides a generalized approach to optimize the design of any type of couplers integrated to infrared photodetectors.

  • PDF

Assessment of Composite Material Flaws on the Type III Cylinders for Compressed Natural Gas Vehicles (압축천연가스자동차용 Type III 용기의 복합재 결함 평가)

  • Kim, Young-Seob;Kim, Lae-Hyun;Yang, Dong-Ju
    • Journal of Energy Engineering
    • /
    • v.20 no.2
    • /
    • pp.90-95
    • /
    • 2011
  • This study was conducted to judge requalification of cylinders by assessing composite flaws such as scratches, cuts, and gouges damaging on the composite of Type III cylinders for compressed natural gas vehicles. As a result of the flaw tolerance test, all specimens have satisfied with minimum requirement cycles according to damage levels based on ISO 19078 and cyclic performance for pressure showed beyond twenty thousands in damage level 1 and 2, and did eighteen thousands to twenty-one thousands in damage level 3. Eight of twelve specimens failed the test due to composite flaws and the rest of the cylinders failed regardless of flaws. The results of Finite Element Method followed by the computer simulation indicated that the stress of 79.5 MPa calculated on the flaw model of $1.25\;mm{\times}200\;mm$ and the stress of 66.6 MPa on the non-flaw model when the service pressure applied to inside of cylinder. The difference between the models is about 19.37%. We concluded that this difference influences fatigue life and this flaw model is a critical value affecting cyclic performance of cylinders.

Numerical Study on Vertical Stress Estimation for Panel Pillars at Room and Pillar Mines (주방식 광산의 패널 광주 수직응력 추정을 위한 수치해석 연구)

  • Yoon, Dong-Ho;Song, Jae-Joon
    • Tunnel and Underground Space
    • /
    • v.30 no.5
    • /
    • pp.473-483
    • /
    • 2020
  • This paper examines the vertical stress change concentrated on mine pillar which occurs due to the stress disturbance from opening excavation at room and pillar mine by FLAC3D, a finite difference method (FDM) software. The mesh size combination is decided with a careful consideration of relative error and run-time, then its performance is verified. A series of numerical analyses is conducted and the vertical stress at central pillar was observed for the test cases of 1×1 to 11×11 mine pillars, 40 m to 320 m depth with 40 m difference. The results show that the vertical stress of pillar approaches to the similar value with the value estimated by tributary area theory(TAT) when the development area (NP) is increased or the height of overburden (HOB) is decreased, while it is overestimated in the opposite case. Furthermore, it also represents that the vertical stress factor (VSF) converges to a specific value when the depth is increased whille keeping the development area identical.

Moisture Content Change of Korean Red Pine Logs During Air Drying: II. Prediction of Moisture Content Change of Korean Red Pine Logs under Different Air Drying Conditions (소나무 원목의 천연건조 중 함수율 변화: II. 소나무 원목의 천연건조 중 함수율 변화 예측)

  • HAN, Yeonjung;CHANG, Yoon-Seong;EOM, Chang-Deuk;LEE, Sang-Min
    • Journal of the Korean Wood Science and Technology
    • /
    • v.47 no.6
    • /
    • pp.732-750
    • /
    • 2019
  • Air drying was carried out on 15 Korean red pine logs to provide a prediction model of the moisture content (MC) change in the wood during drying. The final MC was 17.4% after 880 days since the beginning of air drying in the summer for 6 Korean red pine logs with 68.7% initial MC. The final MC was 16.0% after 760 days since the beginning of air drying in the winter for 9 Korean red pine logs with 35.8% initial MC. A regression model with R-squared of 0.925 was obtained as a result of multiple regression analyses with initial MC, top diameter, temperature, relative humidity, and wind speed as independent variable and and MC change during air drying as dependent variable. The initial MC and top diameter, which is the characteristic of Korean red pine, have greater effect on the MC decrease during air drying compared to meteorological factors such as the temperature, relative humidity, and wind speed. Two-dimensional mass transfer analysis was performed to predict the MC distribution of Korean red pine logs during air drying. Two prediction models with different air drying days and different meteorological factors for the determination of the diffusion coefficient and surface emission coefficient were presented. The error between the different two methods ranged from 0.1 to 0.8% and the difference from the measured value ranged from 2.2 to 3.6%. By measuring the internal MC during air drying of Korean pine logs with various initial MC and diameter, and calculating the moisture transfer coefficient in wood for each meteorological condition, the error of the prediction model can be reduced.

Seismic Performance Evaluation of the Underground Utility Tunnel by Response Displacement Method and Response History Analysis (응답변위법과 응답이력해석법을 이용한 지중 공동구의 내진성능 평가)

  • Kwon, Ki-Yong;Lee, Jin-Sun;Kim, Yong-Kyu;Youn, Jun-Ung;Jeong, Soon-Yong
    • Journal of the Korean Geotechnical Society
    • /
    • v.36 no.11
    • /
    • pp.119-133
    • /
    • 2020
  • Underground utility tunnel, the most representative cut and cover structure, is subjected to seismic force by displacement of the surrounding soil. In 2020, Korea Infrastructure Safety Corporation has published "Seismic Performance Evaluation Guideline for Existing Utility Tunnel." This paper introduces two seismic evaluation methods, RDM (Response Displacement Method) and RHA (Response History Analysis) adopted in the guide and compares the methods for an example of an existing utility tunnel. The test tunnel had been constructed in 1988 and seismic design was not considered. RDM is performed by single and double cosine methods based on the velocity response spectrum at the base rock. RHA is performed by finite difference analysis that is able to consider nonlinear behavior of soil and structure together in two-dimensional plane strain condition. The utility tunnel shows elastic behavior for RDM, but shows plastic hinge for RHA under the collapse prevention level earthquake.