• Structural Engineering and Mechanics
• /
• v.57 no.3
• /
• pp.543-567
• /
• 2016

This study attempts to address the buckling and free vibration characteristics of an isotropic cylindrical panel subjected to non-uniform temperature rise using numerical approach. Finite element analysis has been used in the present study. The approach involves three parts, in the first part non-uniform temperature field is obtained using heat transfer analysis, in the second part, the stress field is computed under the thermal load using static condition and, the last part, the buckling and pre-stressed modal analysis are carried out to compute critical buckling temperature as well as natural frequencies and associated mode shapes. In the present study, the effect of non-uniform temperature field, heat sink temperatures and in-plane boundary constraints are considered. The relation between buckling temperature under uniform and non-uniform temperature fields has been established. Results revealed that decrease (Case (ii)) type temperature variation field influences the fundamental buckling mode shape significantly. Further, it is observed that natural frequencies under free vibration state, decreases as temperature increases. However, the reduction is significantly higher for the lowest natural frequency. It is also found that, with an increase in temperature, nodal and anti-nodal positions of free vibration mode shapes is shifting towards the location where the intensity of the heat source is high and structural stiffness is low.

• Earthquakes and Structures
• /
• v.16 no.2
• /
• pp.221-234
• /
• 2019

Today, many important concrete face rockfill dams (CFRDs) have been built on the world, and some of these important structures are located on the strong seismic regions. In this reason, examination and monitoring of these water construction's seismic behaviour is very important for the safety and future of these dams. In this study, the nonlinear seismic behaviour of Ilısu CFR dam which was built in Turkey in 2017, is investigated for various reservoir water heights taking into account 1995 Kobe near-fault and far-fault ground motions. Three dimensional (3D) finite difference model of the dam is created using the FLAC3D software that is based on the finite difference method. The most suitable mesh range for the 3D model is chosen to achieve the realistic numerical results. Mohr-Coulomb nonlinear material model is used for the rockfill materials and foundation in the seismic analyses. Moreover, Drucker-Prager nonlinear material model is considered for the concrete slab to represent the nonlinearity of the concrete. The dam body, foundation and concrete slab constantly interact during the lifetime of the CFRDs. Therefore, the special interface elements are defined between the dam body-concrete slab and dam body-foundation due to represent the interaction condition in the 3D model. Free field boundary condition that was used rarely for the nonlinear seismic analyses, is considered for the lateral boundaries of the model. In addition, quiet artificial boundary condition that is special boundary condition for the rigid foundation in the earthquake analyses, is used for the bottom of the foundation. The hysteric damping coefficients are separately calculated for all of the materials. These special damping values is defined to the FLAC3D software using the special fish functions to capture the effects of the variation of the modulus and damping ratio with the dynamic shear-strain magnitude. Total 4 different reservoir water heights are taken into account in the seismic analyses. These water heights are empty reservoir, 50 m, 100 m and 130 m (full reservoir), respectively. In the nonlinear seismic analyses, near-fault and far-fault ground motions of 1995 Kobe earthquake are used. According to the numerical analyses, horizontal displacements, vertical displacements and principal stresses for 4 various reservoir water heights are evaluated in detail. Moreover, these results are compared for the near-fault and far-faults earthquakes. The nonlinear seismic analysis results indicate that as the reservoir height increases, the nonlinear seismic behaviour of the dam clearly changes. Each water height has different seismic effects on the earthquake behaviour of Ilısu CFR dam. In addition, it is obviously seen that near-fault earthquakes and far field earthquakes create different nonlinear seismic damages on the nonlinear earthquake behaviour of the dam.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.21 no.11
• /
• pp.2954-2965
• /
• 1996

In this paper, the electromagnetic field characteristics of leaky caxial cable are anlyzed by using the finite difference-time domain(FDTD) technique. Finite difference equations of Maxwell's equations are definedin cylindrical coordinate systems. To simulate the open boundary problem like a free space, the Mur's Absorbing Boundary condition(Mur-ABC) is also used. After modeling the leaky coaxial cable with the three dimensional grid structure, the transient response of the field distribution and the current distribution, the field pattern, the coupling effect are depicted in the time domain.

• 유체기계공업학회:학술대회논문집
• /
• 2005.12a
• /
• pp.670-677
• /
• 2005

In this study, we analyzed the three dimensional unsteady flow field around the motor cooling fan using the unsteady lifting surface theory. We obtained the flow rate for various geometries of fan from the calculated results of velocity field. For the data of design parameter and rotating speed(rpm) of the fan, we can predict the flow rate of the motor cooling fan with thin thickness through numerical analysis without the experimental data of the free stream velocity which is a boundary condition of flow field. the numerical results showed the flow rate within 10% of error in comparison with experimental results. The radial fans, which are often used as internal motor fan were also investigated with the same procedure.

• Earthquakes and Structures
• /
• v.13 no.6
• /
• pp.573-588
• /
• 2017

The earthquake input is required when the soil-structure interaction (SSI) analysis is performed by the direct finite element method. In this paper, the earthquake is considered as the obliquely incident plane body wave arising from the truncated linearly elastic layered half space. An earthquake input method is developed for the time-domain three-dimensional SSI analysis. It consists of a new site response analysis method for free field and the viscous-spring artificial boundary condition for scattered field. The proposed earthquake input method can be implemented in the process of building finite element model of commercial software. It can result in the highly accurate solution by using a relatively small SSI model. The initial condition is considered for the nonlinear SSI analysis. The Daikai subway station is analyzed as an example. The effectiveness of the proposed earthquake input method is verified. The effect of the obliquely incident earthquake is studied.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.33 no.3
• /
• pp.93-100
• /
• 2021

Since the IB (Immersed Boundary) method, which can perform coupling analysis with objects and fluids having an impermeable boundary of arbitrary shape on a fixed grid system, has been developed, the IB method in various CFD models is increasing. The representative IB methods are the directing-forcing method and the ghost cell method. The directing-forcing type method numerically satisfies the boundary condition from the fluid force calculated at the boundary surface of the structure, and the ghost-cell type method is a computational method that satisfies the boundary condition through interpolation by placing a virtual cell inside the obstacle. These IB methods have a disadvantage in that the computational algorithm is complex. In this study, the simplified immersed boundary (SIB) method enables the analysis of temporary structures on a fixed grid system and is easy to expand to three proposed dimensions. The SIB method proposed in this study is based on a one-field model for immiscible two-phase fluid that assumes that the density function of each phase moves with the center of local mass. In addition, the volume-weighted average method using the density function of the solid was applied to handle moving solid structures, and the CIP method was applied to the advection calculation to prevent numerical diffusion. To examine the analysis performance of the proposed SIB method, a numerical simulation was performed on an object falling to the free water surface. The numerical analysis result reproduced the object falling to the free water surface well.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.04a
• /
• pp.116-119
• /
• 2007

In this study the unsteady aerodynamic analysis of a hovering helicopter rotor is performed. For the accurate flow field analysis Euler equations and the free wake method are coupled. The Euler equations are solved to find the pressure distribution around the rotor, and free wake method is used to give the boundary condition for the solution of Euler equations. Also, vortex strength and wake motion after the rotor are simulated by the free wake method. The accuracy of the present method is compared with the source sink model. The present method is applied to the hovering Caradonna-Tung rotor and compared with experimental results.

• The Magazine of the Society of Air-Conditioning and Refrigerating Engineers of Korea
• /
• v.14 no.2
• /
• pp.108-117
• /
• 1985

In this study, a numerical analysis was performed for the natural convection heat transfer in a vertical channel which was consisted of two finite-thickness vertical walls with heat source. The ratio of the thermal conductivity of wall to air played an important role in the analysis. The case for which one side wall has protrusion resistances was also examined. The governing equations for the system was discretized by control volume formulation and solved by SIMPLE method. As the result of this study, it was found that the uniform heat flux boundary condition could be applied when the conductivity ratio was below approximately 50 and the uniform temperature boundary condition could be used when the conductivity rat io was over approximately 15,000. However, when the conductivity ratio was between 50 and 15,000, the thermal conductivity ratio value should be considered for the analysis. It was also found that the existence of protrusion resistance influenced the thermal field up to the distance of 3-4 times of the protrusion length.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.9 no.2 s.42
• /
• pp.71-80
• /
• 2005

The solution to a liquid sloshing problem is challenge to the field of engineering. This is not only because the dynamic boundary condition at the free surface is nonlinear, but also because the position of the free surface varies with time in a manner not known a priori. Therefore, this nonlinear phenomenon, which is characterized by the oscillation of the unrestrained free surface of the fluid, is a difficult mathematical problem to solve numerically and analytically. In this study, three-dimensional boundary element method(BEM), which is based on the so-called an arbitrary Lagrangian-Eulerian(ALE) approach for the fluid flow problems with a free surface, was formulated to solve the behavior of the nonlinear free surface motion. An ALE-BEM has the advantage to track the free surface along any prescribed paths by using only one displacement variable, even for a three-dimensional problem. Also, some numerical examples were presented to demonstrate the validity and the applicability of the developed procedure.

• Journal of Ocean Engineering and Technology
• /
• v.34 no.6
• /
• pp.406-418
• /
• 2020

In this study, a nonlinear wave simulation code is developed using a higher-order spectral (HOS) method. The HOS method is very efficient because it can determine the solution of the boundary value problem using fast Fourier transform (FFT) without matrix operation. Based on the HOS order, the vertical velocity of the free surface boundary was estimated and applied to the nonlinear free surface boundary condition. Time integration was carried out using the fourth order Runge-Kutta method, which is known to be stable for nonlinear free-surface problems. Numerical stability against the aliasing effect was guaranteed by using the zero-padding method. In addition to simulating the initial wave field distribution, a nonlinear adjusted region for wave generation and a damping region for wave absorption were introduced for wave generation simulation. To validate the developed simulation code, the adjusted simulation was carried out and its results were compared to the eighth order Stokes theory. Long-time simulations were carried out on the irregular wave field distribution, and nonlinear wave propagation characteristics were observed from the results of the simulations. Nonlinear adjusted and damping regions were introduced to implement a numerical wave tank that successfully generated nonlinear regular waves. According to the variation in the mean wave steepness, irregular wave simulations were carried out in the numerical wave tank. The simulation results indicated an increase in the nonlinear interaction between the wave components, which was numerically verified as the mean wave steepness. The results of this study demonstrate that the HOS method is an accurate and efficient method for predicting the nonlinear interaction between waves, which increases with wave steepness.