• Earthquakes and Structures
• /
• v.17 no.4
• /
• pp.373-385
• /
• 2019

In this article, different frequently adopted modeling aspects of linear and nonlinear dynamic soil-structure interaction (SSI) are studied on a pile-supported integral abutment bridge structure using the open-source platform OpenSees (McKenna et al. 2000, Mazzoni et al. 2007, McKenna and Fenves 2008) for a 2D domain. Analyzed approaches are as follows: (i) free field input at the base of fixed base bridge; (ii) SSI input at the base of fixed base bridge; (iii) SSI model with two dimensional quadrilateral soil elements interacting with bridge and incident input motion propagating upwards at model bottom boundary (with and without considering the effect of abutment backfill response); (iv) simplified SSI model by idealizing the interaction between structural and soil elements through nonlinear springs (with and without considering the effect of abutment backfill response). Salient conclusions of this paper include: (i) free-field motions may differ significantly from those computed at the base of the bridge foundations, thus put a significant bias on the inertial component of SSI; (ii) conventional modeling of SSI through series of soil springs and dashpot system seems to stay on the safer side under dynamic conditions when one considers the seismic actions on the structure by considering a fully coupled SSI model; (iii) consideration of abutment-backfill in the SSI model positively affects the general response of the bridge, as a result of large passive resistance that may develop behind the abutments.

• The Journal of Engineering Geology
• /
• v.25 no.1
• /
• pp.57-66
• /
• 2015

A numerical analysis was performed of the seepage from and stability of a mine waste-dump slope in Imgi, Busan, considering rainfall intensity. The 40-45° slope angle of the waste dump is relatively steep, and the depth of the waste dump down to bedrock is 7-8 m. The groundwater level was 6.6 m below the surface. Various laboratory tests on samples obtained from the waste dump were performed to determine the input data for seepage and stability analyses of the waste-dump slope during rainfall. The results of seepage analysis for various rainfall intensities using the SEEP/W program show that the wetting front moved down with increasing rainfall duration. When the rainfall intensity was > 50 mm/ hour and the duration was > 24 hours, the waste dump became fully saturated because the wetting front reached the groundwater level. The results of slope stability analysis coupled with seepage analysis using the SLOPE/W program show that the safety factor of the slope decreased as the wetting front moved down due to rainfall infiltration. After continuous rainfall for 5-6 hours, the safety factor of the slope suddenly decreased but then recovered and converged. The sudden decrease was induced by an increase in pore-water pressure and a decrease in matric suction down to a certain depth as the wetting front approached the potential sliding surface.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.17 no.1
• /
• pp.137-145
• /
• 2009

The steady-state kinetics of the selective catalytic reduction (SCR) of $NO_X$ with $NH_3$ has been investigated over a commercial ${V_2}{O_5}/TiO_2$ catalyst. In order to account for the influence of transport effects the kinetics are coupled with a fully transient two-phase 1D+1D monolith channel model. The Langmuir-Hinshelwood (L-H) mechanism is adopted to describe the steady-state kinetic behavior of the ${V_2}{O_5}/TiO_2$ catalyst. The reaction rate expressions are based on previously reported papers and are modified to fit the experimental data. The steady-state chemical reaction scheme used in the present mathematical model has been validated extensively with experimental data of selective $NO_X$ reduction efficiency for a wide range of inlet conditions such as space velocity, oxygen concentrations, water concentration, and $NO_2/NO$ ratio. The parametric investigations are performed to examine how the $NH_3$ slip from a SCR $DeNO_X$ catalyst and the conversion of $NO_X$ are affected by the reaction temperature, $NH_3/NO_X$ feed ratio, and space velocity for feed gas compositions with $NO_2/NO_X$ ratios of 0 and 0.5.

• Journal of Ocean Engineering and Technology
• /
• v.19 no.5 s.66
• /
• pp.1-15
• /
• 2005

The motion behaviors including hydrodynamic interaction and mechanical coupling effects on multiple-body floating platforms are simulated by using a time domain hull/mooring/riser coupled dynamics analysis program. The objective of this study is to evaluate off-diagonal hydrodynamic interaction effects and mechanical coupling effects on tandem moored FPSO and shuttle taker motions. In the multiple-body floating platforms interaction, hydrodynamic coupling effects with waves and mechanical coupling effects through the connectors should be considered. Thus, in this study, the multiple-body platform motions are calculated by Combined Matrix Method (CMM) as well as Separated Matrix Method (SMM). The advantage of the combined matrix method is that it can include all the 6Nx6N full hydrodynamic and mechanical interaction effects among N bodies. Whereas, due to the larger matrix size, the calculation time of Combined Matrix Method (CMM) is longer than the Separated Matrix Method (SMM). On the other hand, Separated Matrix Method (SMM) cannot include the off-diagonal 6x6 hydrodynamic interaction coefficients although it can fully include mechanical interactions among N bodies. To evaluate hydrodynamic interaction and mechanical coupling effects, tandem moored FPSO and shuttle tanker is simulated by Combined Matrix Method (CMM) and Separated Matrix Method (SMM). The calculation results give a good agreement between Combined Matrix Method (CMM) and Separated Matrix Method (SMM). The results show that the Separated Matrix Method (SMM) is more efficient for tandem moored FPSO and shuttle tanker. In the numerical calculation, the hydrodynamic coefficients are calculated from a 3D diffraction/radiation panel program WAMIT, and wind and current forces are generated by using the respective coefficients given in the OCIMF data sheet.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.20 no.3
• /
• pp.575-592
• /
• 2018

Tunnelling is getting more important solutions for problems induced by the growth population in urban areas. Many studies on tunnelling below existing structure are carried out by many researchers. In general, however, ground water condition is ignored for most of researches using laboratory model test, so far. In case of ground behavior, error can occur if the result of effective stress related to hydraulic condition can't be taken into considerations. In this study, therefore, laboratory model test and the close range photogrammetry were conducted to investigate behaviour of ground and shallow foundation using newly device drainage system which is available to express the ground water condition. Also, numerical analysis was carried out to compare to results from the laboratory model test, and was performed with two methods, one is plastic and the other one is fully coupled analysis. Results from those two methods were compared to that of the laboratory model test.

• Ocean Systems Engineering
• /
• v.5 no.3
• /
• pp.139-160
• /
• 2015

The global performance of the 5MW OC4 semisubmersible floating wind turbine in random waves was numerically simulated by using the turbine-floater-mooring fully coupled and time-domain dynamic analysis program FAST-CHARM3D. There have been many papers regarding floating offshore wind turbines but the effects of second-order wave-body interactions on their global performance have rarely been studied. The second-order wave forces are actually small compared to the first-order wave forces, but its effect cannot be ignored when the natural frequencies of a floating system are outside the wave-frequency range. In the case of semi-submersible platform, second-order difference-frequency wave-diffraction forces and moments become important since surge/sway and pitch/roll natural frequencies are lower than those of typical incident waves. The computational effort related to the full second-order diffraction calculation is typically very heavy, so in many cases, the simplified approach called Newman's approximation or first-order-wave-force-only are used. However, it needs to be justified against more complete solutions with full QTF (quadratic transfer function), which is a main subject of the present study. The numerically simulated results for the 5MW OC4 semisubmersible floating wind turbine by FAST-CHARM3D are also extensively compared with the DeepCWind model test results by Technip/NREL/UMaine. The predicted motions and mooring tensions for two white-noise input-wave spectra agree well against the measure values. In this paper, the numerical static-offset and free-decay tests are also conducted to verify the system stiffness, damping, and natural frequencies against the experimental results. They also agree well to verify that the dynamic system modeling is correct to the details. The performance of the simplified approaches instead of using the full QTF are also tested.

• Tribology and Lubricants
• /
• v.36 no.2
• /
• pp.105-115
• /
• 2020

This paper presents a numerical study on the rotordynamic analysis of a dual-spool turbofan engine in the context of blade defect events. The blades of an axial-type aeroengine are typically well aligned during the compressor and turbine stages. However, they are sometimes exposed to damage, partially or entirely, for several operational reasons, such as cracks due to foreign objects, burns from the combustion gas, and corrosion due to oxygen in the air. Herein, we designed a dual-spool rotor using the commercial 3D modeling software CATIA to simulate blade defects in the turbofan engine. We utilized the rotordynamic parameters to create two finite element Euler-Bernoulli beam models connected by means of an inter-rotor bearing. We then applied the unbalanced forces induced by the mass eccentricities of the blades to the following selected scenarios: 1) fully balanced, 2) crack in the low-pressure compressor (LPC) and high pressure compressor (HPC), 3) burn on the high-pressure turbine (HPT) and low pressure compressor, 4) corrosion of the LPC, and 5) corrosion of the HPC. Additionally, we obtained the transient and steady-state responses of the overall rotor nodes using the Runge-Kutta numerical integration method, and employed model reduction techniques such as component mode synthesis to enhance the computational efficiency of the process. The simulation results indicate that the high-vibration status of the rotor commences beyond 10,000 rpm, which is identified as the first critical speed of the lower speed rotor. Moreover, we monitored the unbalanced stages near the inter-rotor bearing, which prominently influences the overall rotordynamic status, and the corrosion of the HPC to prevent further instability. The high-speed range operation (>13,000 rpm) coupled with HPC/HPT blade defects possibly presents a rotor-case contact problem that can lead to catastrophic failure.

• Journal of Korean Society of Steel Construction
• /
• v.26 no.4
• /
• pp.323-334
• /
• 2014

In this paper, experimental program and associated numerical study were carried out to evaluate the fire resistance of unprotected concrete-filled rectangular steel tubular (CFT) columns subjected to the standard fire. The key testing parameters included the length effect, the load ratio, and the sectional dimensions of the CFT columns. Temperature distribution and axial deformation of the CFT column specimens were measured and analyzed. Rather early local buckling of steel tubes was observed in all the specimens. This caused subsequent load transfer from steel tube to concrete, and eventually triggered concrete crushing, or complete loss of the load bearing capacity of the column. This implies that the limit state of local buckling as well as overall flexural buckling should be incorporated in fire design procedure. As expected, the fire resistance time of specimen with higher load ratio consistently lessened. The prediction of fire resistance time of unprotected CFT columns based on the limiting steel temperature in current design codes or the formula proposed by previous studies is slightly conservative compared to the fire test results available. To establish the finite element analysis model that can be used to predict the thermal and structural behaviour of unprotected CFT columns in fire, the fully coupled thermal-stress analysis was also tried by using the commercial code ABAQUS. The numerical results showed a reasonable global correlation with the experimental results.

• Journal of Nuclear Fuel Cycle and Waste Technology(JNFCWT)
• /
• v.17 no.2
• /
• pp.183-202
• /
• 2019

For design and performance assessment of a high-level radioactive waste (HLW) disposal system, it is necessary to understand the characteristics of coupled thermo-hydro-mechanical (THM) behavior. However, in previous studies for the Korean Reference HLW Disposal System (KRS), thermal analysis was performed to determine the spacing of disposal tunnels and interval of disposition holes without consideration of the coupled THM behavior. Therefore, in this study, TOUGH2-MP/FLAC3D is used to conduct THM modeling for performance assessment of the Korean Reference HLW Disposal System (KRS). The peak temperature remains below the temperature limit of $100^{\circ}C$ for the whole period. A rapid rise of temperature caused by decay heat occurs in the early years, and then temperature begins to decrease as decay heat from the waste decreases. The peak temperature at the bentonite buffer is around $96.2^{\circ}C$ after about 3 years, and peak temperature at the rockmass is $68.2^{\circ}C$ after about 17 years. Saturation of the bentonite block near the canister decreases in the early stage, because water evaporation occurs owing to temperature increase. Then, saturation of the bentonite buffer and backfill increases because of water intake from the rockmass, and bentonite buffer and backfill are fully saturated after about 266 years. The stress is calculated to investigate the effect of thermal stress and swelling pressure on the mechanical behavior of the rockmass. The calculated stress is compared to a spalling criterion and the Mohr-Coulumb criterion for investigation of potential failure. The stress at the rockmass remains below the spalling strength and Mohr-Coulumb criterion for the whole period. The methodology of using the TOUGH2-MP/FLAC3D simulator can be applied to predict the long-term behavior of the KRS under various conditions; these methods will be useful for the design and performance assessment of alternative concepts such as multi-layer and multi-canister concepts for geological spent fuel repositories.

• Proceedings of the Korean Geotechical Society Conference
• /
• 1998.05a
• /
• pp.35-81
• /
• 1998

Distinct Element Method(DEM) has a great advantage to model the discontinuous behaviour of jointed rock masses such as rotation, sliding, and separation of rock blocks. Geometrical data of joints by a field monitoring is not enough to model the jointed rock mass though the results of DE analysis for the jointed rock mass is most sensitive to the distributional properties of joints. Also, it is important to use a properly joint law in evaluating the stability of a jointed rock mass because the joint is considered as the contact between blocks in DEM. In this study, a stochastic modelling technique is developed and the dilatant rock joint is numerically modelled in order to consider th geometrical and mechanical properties of joints in DE analysis. The stochastic modelling technique provides a assemblage of rock blocks by reproducing the joint distribution from insufficient joint data. Numerical Modelling of joint dilatancy in a edge-edge contact of DEM enable to consider not only mechanical properties but also various boundary conditions of joint. Preprocess Procedure for a stochastic DE model is composed of a statistical process of raw data of joints, a joint generation, and a block boundary generation. This stochastic DE model is used to analyze the effect of deviations of geometrical joint parameters on .the behaviour of jointed rock masses. This modelling method may be one tool for the consistency of DE analysis because it keeps the objectivity of the numerical model. In the joint constitutive law with a dilatancy, the normal and shear behaviour of a joint are fully coupled due to dilatation. It is easy to quantify the input Parameters used in the joint law from laboratory tests. The boundary effect on the behaviour of a joint is verified from shear tests under CNL and CNS using the numerical model of a single joint. The numerical model developed is applied to jointed rock masses to evaluate the effect of joint dilation on tunnel stability.