• Journal of the Korean Geotechnical Society
• /
• v.39 no.5
• /
• pp.51-63
• /
• 2023

Numerical analysis was conducted to determine the effect of soil behavior by thawing and freezing of seasonal frozen soil on pile foundations. The analysis was performed using the finite element method (FEM) to simulate soil-pile interaction based on the atmosphere temperature change. Thermomechanical coupled modeling using FEM was applied with the temperature-dependent nonlinear properties of the frozen soil. The analysis model cases were applied to the MCR and HDP models to simulate the elastoplastic behavior of soil. The numerical analysis results were analyzed and compared with various conditions having different length and width sizes of the pile. The results of the numerical analysis showed t hat t he HDP model was relat ively passive, and t he aspect and magnit ude of t he bearing capacit y and displacement of the pile head were similar depending on the length and width of the pile conditions. The vertical displacement of the pile head by thawing and freezing of the ground showed a large variation in displacement for shorter length conditions. In the MCR model, the vertical displacement appeared in the maximum thaw settlement and frost heaving of 0.0387 and 0.0277 m, respectively. In the HDP model, the vertical displacement appeared in the maximum thaw settlement and frost heaving of 0.0367 and 0.0264 m, respectively. The results of the pile bearing capacity for the two elastoplastic models showed a larger difference in the width condition than the length condition of the pile, with a maximum of about 14.7% for the width L condition, a maximum of about 5.4% for M condition, and a maximum of about 5.3% for S condition. The significance of the effect on the displacement of the pile head and the bearing capacity depended on the pile-soil contact area, and the difference depended on the presence or absence of an active layer in the soil and its thickness.

• Journal of the Korean Society of Safety
• /
• v.11 no.4
• /
• pp.143-150
• /
• 1996

This is an explicit-Implicit, finite element analysis for linear as well as nonlinear hygrothermal stress problems. Additional features, such as moisture diffusion equation, crack element and virtual crack extension(VCE ) method for evaluating J-integral are implemented in this program. The Linear Elastic Fracture Mechanics(LEFM) Theory is employed to estimate the crack driving force under the transient condition for and existing crack. Pores in materials are assumed to be saturated with moisture in the liquid form at the room temperature, which may vaporize as the temperature increases. The vaporization effects on the crack driving force are also studied. The Ideal gas equation is employed to estimate the thermodynamic pressure due to vaporization at each time step after solving basic nodal values. A set of field equations governing the time dependent response of porous media are derived from balance laws based on the mixture theory Darcy's law Is assumed for the fluid flow through the porous media. Perzyna's viscoplastic model incorporating the Von-Mises yield criterion are implemented. The Green-Naghdi stress rate is used for the invariant of stress tensor under superposed rigid body motion. Isotropic elements are used for the spatial discretization and an iterative scheme based on the full newton-Raphson method is used for solving the nonlinear governing equations.

• Geomechanics and Engineering
• /
• v.23 no.2
• /
• pp.151-163
• /
• 2020

Grouting method is an effective way of reinforcing cracked rock masses and plugging water gushing. Current grouting diffusion models are generally developed for horizontal cracks, which is contradictory to the fact that the crack generally occurs in rock masses with irregular spatial distribution characteristics in real underground environments. To solve this problem, this study selected a cement-sodium silicate slurry (C-S slurry) generally used in engineering as a fast-curing grouting material and regarded the C-S slurry as a Bingham fluid with time-varying viscosity for analysis. Based on the theory of fluid mechanics, and by simultaneously considering the deadweight of slurry and characteristics of non-uniform spatial distribution of viscosity of fast-curing grouts, a theoretical model of slurry diffusion in an oblique crack in rock masses at constant grouting rate was established. Moreover, the viscosity and pressure distribution equations in the slurry diffusion zone were deduced, thus quantifying the relationship between grouting pressure, grouting time, and slurry diffusion distance. On this basis, by using a 3-d finite element program in multi-field coupled software Comsol, the numerical simulation results were compared with theoretical calculation values, further verifying the effectiveness of the theoretical model. In addition, through the analysis of two engineering case studies, the theoretical calculations and measured slurry diffusion radius were compared, to evaluate the application effects of the model in engineering practice. Finally, by using the established theoretical model, the influence of cracking in rock masses on the diffusion characteristics of slurry was analysed. The results demonstrate that the inclination angle of the crack in rock masses and azimuth angle of slurry diffusion affect slurry diffusion characteristics. More attention should be paid to the actual grouting process. The results can provide references for determining grouting parameters of fast-curing grouts in engineering practice.

• Journal of Electrical Engineering and Technology
• /
• v.13 no.6
• /
• pp.2402-2411
• /
• 2018

Fluoronitriles-$CO_2$ gas mixtures are promising alternatives to $SF_6$ in environmentally-friendly gas-insulated transmission lines (GILs). Insulating gas heat transfer characteristics are of major significance for the current-carrying capacity design and operational state monitoring of GILs. In this paper, a three-dimensional calculation model was established for a GIL using the thermal-fluid coupled finite element method. The calculated results showed close agreement with experimentally measured data. The temperature distribution of a GIL filled with the Fluoronitriles-$CO_2$ mixture was obtained and compared with those of GILs filled with $CO_2$ and $SF_6$. Furthermore, the effects of the mixture ratio of the component gases and the gas pressure on the temperature rise and current-carrying capacity of the GIL were analyzed. Results indicated that the heat transfer performance of the Fluoronitriles-$CO_2$ gas mixture was better than that of $CO_2$ but worse than that of $SF_6$. When compared with $SF_6$, use of the Fluoronitriles-$CO_2$ gas mixture caused a reduction in the GIL's current-carrying capacity. In addition, increasing the Fluoronitriles gas component ratio or increasing the pressure of the insulating gas mixture could improve the heat dissipation and current-carrying capacity of the GIL. These research results can be used to design environmentally-friendly GILs containing Fluoronitriles-$CO_2$ gas mixtures.

• Journal of the Korean Society of Groundwater Environment
• /
• v.7 no.2
• /
• pp.89-96
• /
• 2000

Slug tests are the most widely used field method for quantification of hydraulic conductivity of porous media. Well recovery is affected by well casing, borehole radii, screened length, hydraulic conductivity, and specific storage of porous media. In this study, a new slug tests model was developed through finite element approximation and the validity and usefulness of the model were tested in various ways. Water level fluctuation in a well under slug test and cons-equent groundwater flow in the surrounding porous medium were appropriately coupled through estimation of well-flux using an iteration technique. Numerical accuracy of the model was verified using the Cooper et al. (1967) solution. The model has advantages in simulations for monitored slug tests, partial penetration, and inclusion of storage factor. Volume coverage of slug tests is significantly affected by storage factor. Magnitude and speed of propagation of head changes from a well increases as storage factor becomes low. It will be beneficial to use type curves of monitored head transients in the surrounding porous formation for estimation of specific storage. As the vertical component of groundwater flow is enhanced, the influence of storage factor on well recovery decreases. For a radial-vertical flow around a partially penetrated well, deviations between hydraulic estimates by various methods and data selection of recovery curve are negligible on practical purposes, whereas the deviations are somewhat significant for a radial flow.

• 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 the Microelectronics and Packaging Society
• /
• v.16 no.2
• /
• pp.43-50
• /
• 2009

Although thin PCBs(Printed Circuit Boards) have recently been required for high density interconnection, high electrical performance, and low manufacturing cost, the utilization of thin PCBs is severely limited by warpage and reliability issues. Warpage of the thin PCB leads to failure in solder-joints and chip. The $B^2it$(Buried Bump Interconnection Technology) for PCB has been developed to achieve a competitive manufacturing price. In this study, chip temperature, package warpage, chip stress and solder-joints stress characteristics of the PCB prepared with $B^2it$ process have been calculated using thermo-mechanical coupled analysis by the FEM(Finite Element Method). FEM computation was carried out with the variations in bump shapes and kinds of materials under 1.5W power of chip and constant convection heat transfer. The results show that chip temperature distribution reached more quickly steady-state status with PCB prepared with $B^2it$ process than PCB prepared with conventional via interconnection structure. Although $B^2it$ structures are effective on low package warpage and chip stress, with high strength bump materials arc disadvantage for low stress of solder-joints. Therefore, it is recommended that optimized bump shapes and materials in PCB design should be considered in terms of reliability characteristics in the packaging level.