• Computers and Concrete
• /
• 제16권5호
• /
• pp.669-682
• /
• 2015

Degradation of reinforced concrete (RC) structures due to chloride penetration followed by reinforcement corrosion has been a serious problem in civil engineering for many years. The numerical simulation methods at present are mainly finite element method (FEM) and finite difference method (FDM), which are based on mesh. Mesh generation in engineering takes a long time. In the present article, the numerical solution of chloride transport in concrete is analyzed using radial point interpolation method (RPIM) and element-free Galerkin (EFG). They are all meshless methods. RPIM utilizes radial polynomial basis, whereas EFG uses the moving least-square approximation. A Galerkin weak form on global is used to attain the discrete equation, and four different numerical examples are presented. MQ function and appropriate parameters have been proposed in RPIM. Numerical simulation results are compared with those obtained from the finite element method (FEM) and analytical solutions. Two case of chloride transport in full saturated and unsaturated concrete are analyzed to test the practical applicability and performance of the RPIM and EFG. A good agreement is obtained among RPIM, EFG, and the experimental data. It indicates that RPIM and EFG are reliable meshless methods for prediction of chloride concentration in concrete structures.

• 소성∙가공
• /
• 제9권4호
• /
• pp.395-403
• /
• 2000

Formation of Shear Band under the adiabatic condition is widely observed In the engineering materials during rapidly forming process lot a thermally rate-dependent material. The shear band stems from evolution of a narrow region in which an intensive plastic flow occurs. The shear band often plays a role of a precursor of the ductile fracture during a forming process. The objective of this study is to investigate the localization behavior using numerical method. In this work, the implicit finite difference scheme is employed due to the ease of convergence and the numerical stability It is noted that physical and mechanical properties of materials determine how the shear band is formed and then localized. Material properties can be characterized with inertia number dissipation number and diffusion number. It is observed that the dimensionless numbers effect on localization. Using a parametric study, comparison was made between CRS-1018 steel with WHA (tungsten heavy alloy). The deformation behavior of material in this study include an isotropic hardening as well as thermal softening. Moreover, this study suggests that a kinematic hardening constitutive relation be required to predict a more accurate strain level at a shear band.

• 한국지반공학회:학술대회논문집
• /
• 한국지반공학회 2009년도 세계 도시지반공학 심포지엄
• /
• pp.1441-1446
• /
• 2009

Numerical analysis is a powerful method in evaluating the soil-pile-structure interaction under the dynamic loading, and this approach has been applied to the practical area due to the development of computer technology. Finite Difference Method, one of the most popular numerical methods, is sensitive to the shape and the number of mesh. However, the trial and error approach is conducted to obtain the accurate results and the reasonable simulation time because of the lack of researches about mesh size and the number. In this study, FLAC 3D v3.1 program(FDM) is used to simulate the dynamic pile model tests, and the numerical results are compared with the 1G shaking table tests results. With the different size and shape of mesh, the responses of pile behavior and the simulation time are estimated, and the optimum mesh sizes in dynamic analysis of single pile is studied.

• 정보저장시스템학회논문집
• /
• 제5권1호
• /
• pp.19-35
• /
• 2009

The present work is an experimental and analytical study on a flexible disk rotating close to a rigid rotating disk in open air. In the analytical study, the air flow in the gap between the flexible disk and the rigid disk is modeled using Navier-Stokes and continuity equations while the flexible disk is modeled using the linear plate theory. The flow equations are discretized using the cell centered finite volume method (FVM) and solved numerically with semi-implicit pressure-linked equations (SIMPLE algorithm). The spatial terms in the disk equation are discretized using the finite difference method (FDM) and the time integration is performed using fourth-order Runge-Kutta method. An experimental test-rig is designed to investigate the dynamics of the flexible disk when rotating close to a co-rotating, a counter-rotating and a fixed rigid disk, which works as a stabilizer. The effects of rotational speed, initial gap height and inlet-hole radius on the flexible disk displacement and its vibration amplitude are investigated experimentally for the different types of stabilizer. Finally, the analytical and experimental results are compared.

• Ocean Science Journal
• /
• 제40권4호
• /
• pp.201-208
• /
• 2005

An accurate particle tracking method for a finite difference method model is developed using a constant acceleration method. Being assumed constant temporal and spatial gradients, the new method permits temporal-spatial variability of particle velocity. Test results in a solid rotating flow show that the new method has second-order accuracy. The performance of the new method is compared with that of other methods; the first-order Euler forward method, and the second-order Euler predictor-corrector method. The new method is the most efficient method among the three. It is more accurate and efficient than the other two.

• 한국농공학회논문집
• /
• 제57권6호
• /
• pp.1-7
• /
• 2015

Most of chloride diffusion models based on finite difference method (FDM) could not express the diffusion in horizontal direction at each elevation. To overcome these weakness, two dimensional chloride ion penetration model based on finite element method (FEM) to be able to combine various multi-physics simultaneously was suggested by introducing moving mesh technique. To avoid the generation of mesh being able to be distorted depending on the relative movement of water level to static concrete, a rectangular type of mesh was intentionally adopted and the total number of meshes was empirically selected. The simulated results showed that the contents of surface chloride decreased following to the increase of elevation in the top part of low sea level, whereas there were no changes in the bottom part of low level. In the DuraCrete model, the diffusion coefficient of splashed zone is generally smaller than submerged zone, whereas the trend of Life365 model is reverse. Therefore, it could be understood that the developed model using moving mesh technique effectively reflects $DuraCrete^{TM}$ model rather than $Life365^{TM}$ model. In the future, the model will be easily expanded to be combined with various multi-physics models considering water evaporation, heat of hydration, irradiation effect of sun and so on because it is based on FEM.

• 한국소음진동공학회논문집
• /
• 제15권2호
• /
• pp.206-212
• /
• 2005

In many machines handling lightweight and flexible media, such as magnetic tape drives, xerographic copiers and sewing machines, the media must transit an open space. It is important to predict the static and dynamic behavior of the sheets with a high degree of reliability. The nonlinear theory of the dynamic elastica has often been used to a nonlinear dynamic deflection model. In this paper, the governing equation is derived and simulated by the finite difference method. The parametric cubic curve is applied for defining the guide shape. The dynamic contact conditions suggested by Klarbring is used to predict the direction of the flexible media according to the initial velocity and the friction coefficient. The analysis is also compared to the conventional model, showing that after contacting a $45^{\circ}$ wall, the directions of flexible media of two models are different.

• Geomechanics and Engineering
• /
• 제5권1호
• /
• pp.37-55
• /
• 2013

An important issue in the design of soil-nailing systems, as long-term retaining walls, is to assess their stability during seismic events. As such, this study is aimed at simulating the dynamic behavior and failure pattern of nailed structures using two series of numerical analyses, namely dynamic time history and pseudo-static. These numerical simulations are performed using the Finite Difference Method (FDM). In order to consider the actual response of a soil-nailed structure, nonlinear soil behaviour, soil-structure interaction effects, bending resistance of structural elements and construction sequences have been considered in the analyses. The obtained results revealed the efficiency of both analysis methods in simulating the seismic failure mechanism. The predicted failure pattern consists of two sliding blocks enclosed by three slip surfaces, whereby the bottom nails act as anchors and the other nails hold a semi-rigid soil mass. Moreover, it was realized that an increase in the length of the lowest nails is the most effective method to improve seismic stability of soil-nailed structures. Therefore, it is recommended to first estimate the nails pattern for static condition with the minimum required static safety factor. Then, the required seismic stability can be obtained through an increase in the length of the lowest nails. Moreover, placement of additional long nails among lowest nails in existing nailed structures can be considered as a simple retrofitting technique in seismic prone areas.

• Nuclear Engineering and Technology
• /
• 제53권5호
• /
• pp.1403-1415
• /
• 2021

The transient capability of a SP₃ based pin-wise core analysis code SPHINCS is developed and verified through the analyses of the C5G7-TD benchmark. Spatial discretization is done by the fine mesh finite difference method (FDM) within the framework of the coarse mesh finite difference (CMFD) formulation. Pin size fine meshes are used in the radial fine mesh kernels. The time derivatives of the odd moments in the time-dependent SP₃ equations are neglected. The pin homogenized group constants and Super Homogenization (SPH) factors generated from the 2D single assembly calculations at the unrodded and rodded conditions are used in the transient calculations via proper interpolation involving the approximate flux weighting method for the cases that involve control rod movement. The simplifications and approximations introduced in SPHINCS are assessed and verified by solving all the problems of C5G7-TD and then by comparing with the results of the direct whole core calculation code nTRACER. It is demonstrated that SPHINCS yields accurate solutions in the transient behaviors of core power and reactivity.

• Corrosion Science and Technology
• /
• 제18권3호
• /
• pp.102-109
• /
• 2019

The hydrogen diffusion and trapping model with a numerical finite difference method (FDM) was modified and extended to accommodate $H_2S$ corrosion and scale forming processes of high-strength steel under tensile stress condition. The newly proposed diffusion model makes it possible to clearly understand combined effect of tensile stress and $H_2S$ corrosion process on hydrogen diffusion behaviors. The core concept of this theoretical approach is that overall diffusion behavior is separated into diffusion process through two respective layers: an outer sulfide scale and an inner steel matrix. Diffusion coefficient values determined by curve-fitting permeation data reported previously with the newly proposed diffusion model indicate that the application of tensile stress can contribute to continual increase in the diffusivity in the sulfide scale with a high density of defect. This suggests that the scale with a lower stability under the stress condition can be a key parameter to enhance hydrogen influx in the steel matrix. Consequently, resistance to hydrogen assisted cracking of the steel under tensile stress can be decreased significantly.