• Journal of the Korean Geotechnical Society
• /
• v.22 no.9
• /
• pp.37-43
• /
• 2006

Dynamic numerical analysis of geotechnical problems requires a way to simulate the decrease of energy as the domain of interest gets larger. This phenomenon is usually referred to as radiation damping or geometric attenuation and it is distinguished from material damping in which elastic energy is actually dissipated by viscous, hysteretic, or other mechanism. The fact that the domain of analysis in numerical modeling must be chosen, however, causes a need for special attention at the boundary. This observation leads directly to the idea of determining the dynamic response of the interior region from a finite model consisting of the interior region subjected to a boundary condition which ensures that all energy arriving at the boundary is absorbed. This paper presents a simple methodology to simulate transmitting boundaries condition using viscoelastic infinite elements within the recently developed "OpenSees" finite element code. The methodology used here provides that the level of absorption for traveling waves is efficient enough for practical purposes, but unsatisfactory for the case of sharp incident angles. The effectiveness of the infinite elements for the absorption of incident waves at boundaries is evaluated via example analysis.

• Journal of Korea Foundry Society
• /
• v.10 no.4
• /
• pp.299-308
• /
• 1990

Ductile cast iron has a good ductility and ductility and toughness than those of gray cast iron, because the shape of graphite is spheroidal. It has been reported that the strengthening and toughening of the ductile cast iron was resulted from the good modification of various matrix structures obtained by the heat treatment or addition of alloying elements. This study aims to investigate the effect of various special heat treatment(Cyclic Heat Treatment, Intermediate Heat Treatment, Step Quenching), austempering and alloying element(Ni) on the strength and toughness of ductile cast iron. The results obtained from this study are summarized as follows. 1) With addition of Ni, the amount of pearlite or bainite were increased and the morphologies of pearlite or bainite became finer by special heat treatments. 2) As the Ni added and not added ductile cast iron were treated by austenitizing at $900^{\circ}C$ and $840^{\circ}C$, in the latter the austenite was mostly formed in the vicinity of eutectic cell boundary, but in the former on the whole matrix. 3) In cyclic heat treatment, the volume fraction of pearlite was increased and the shape of pearlite was fined with increase of the number of cycle. 4) The shape of pearlite was mostly bar type in the intermediate heat treatment, but spheroidal type in step quenching. 5) The mechanical properties of ductile cast iron containing 1.5%Ni austempered at $400^{\circ}C$ for 25min. after austenitizing at $900^{\circ}C$ for 15min. were a good value of hardness 105(HRB), impact energy 12.5(kg.m), tensile strength 112($kg/mm^2$) and elongation 6.8(%).

• Geomechanics and Engineering
• /
• v.31 no.1
• /
• pp.31-52
• /
• 2022

Investigating and evaluating the long-term creep behavior of historical buildings built on seismic zones is of great importance in terms of transferring these structures to future generations. Furthermore, assessing the earthquake behavior of historical structures such as masonry stone bridges is very important for the future and seismic safety of these structures. For this reason, in this study, earthquake analyses of a masonry stone bridge are carried out considering strong ground motions and various water levels. Tokatli masonry stone arch bridge that was built in the 10th century in Turkey-Karabük is selected for three-dimensional (3D) finite difference analyses and this bridge is modeled using FLAC3D software based on the three-dimensional finite difference method. Firstly, each stone element of the bridge is modeled separately and special stiffness parameters are defined between each stone element. Thanks to these parameters, the interaction conditions between each stone element are provided. Then, the Burger-Creep and Drucker-Prager material models are defined to arch material, rockfill material for evaluating the creep and seismic failure behaviors of the bridge. Besides, the boundaries of the 3D model of the bridge are modeled by considering the free-field and quiet boundary conditions, which were not considered in the past for the seismic behavior of masonry bridges. The bridge is analyzed for 6 different water levels and these water levels are 0 m, 30 m, 60 m, 70 m, 80 m, and 90 m, respectively. A total of 10 different seismic analyzes are performed and according to the seismic analysis results, it is concluded that historical stone bridges exhibit different seismic behaviors under different water levels. Moreover, it is openly seen that the water level is of great importance in terms of earthquake safety of historical stone bridges built in earthquake zones. For this reason, it is strongly recommended to consider the water levels while strengthening and analyzing the historical stone bridges.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.10
• /
• pp.1213-1222
• /
• 1999

A finite element method is developed for calculating the temperature and enthalpy distribution and accordingly the solid, liquid and mushy zone in a three-dimensional body subjected to any heat boundary conditions. The method concurrently consider both temperature and enthalpy for consideration of the latent heat effect, differently from other methods of using a special energy balance equation for solving a mushy zone. The developed brick element has eight nodes with one degree of freedom at each node. The numerical method and procedure are verified using the results of one and two dimensional analytic solutions and by other researchers. It is shown that the present method presents a consistent and stable results in either abrupt or ranged phase change problems. Moreover, the numerical results by the present method are hardly effected by the calculation time steps which otherwise are difficult to determine in most phase change problems. Finally, as a three-dimensional application, a T-shaped body of a phase change is presented and the temperature and enthalpy variation along the time are solved.

• Smart Structures and Systems
• /
• v.6 no.4
• /
• pp.373-388
• /
• 2010

The updated finite element model of K$\ddot{o}$m$\ddot{u}$rhan Highway Bridge on the Firat River located on the $51^{st}$ km of Elazi$\breve{g}$-Malatya highway is obtained by using analytical and experimental results. The 2D and 3D finite element model of the bridge is created by using SAP2000 structural analyses software, and the dynamic characteristics of the bridge are determined analytically. The experimental measurements are carried out by Operational Modal Analysis Method under traffic induced vibrations and the dynamic characteristics are obtained experimentally. The vibration data are gathered from the both box girder and the deck of the bridge, separately. Due to the expansion joint in the middle of the bridge, special measurement points are selected when experimental test setups constitute. Measurement duration, frequency span and effective mode number are determined by considering similar studies in literature. The Peak Picking method in the frequency domain is used in the modal identification. At the end of the study, analytical and experimental dynamic characteristic are compared with each other and the finite element model of the bridge is updated by changing some uncertain parameters such as material properties and boundary conditions. Maximum differences between the natural frequencies are reduced from 10% to 2%, and a good agreement is found between natural frequencies and mode shapes after model updating.

• Journal of the Society of Naval Architects of Korea
• /
• v.37 no.1
• /
• pp.67-81
• /
• 2000

An accurate and efficient numerical method for two-dimensional nonlinear radiation problem has been developed. The wave motion due to a moving body is described by the assumption of ideal fluid flow, and the governing Laplace equation can be effectively solved by the higher-order boundary element method with the help of the GMRES (Generalized Minimal RESidual) algorithm. The intersection or corner problem is resolved by utilizing the so-called discontinuous elements. The implicit trapezoidal rule is used in updating solutions at new time steps by considering stability and accuracy. Traveling waves caused by the oscillating body are absorbed downstream by the damping zone technique. It is demonstrated that the present method for time marching and radiation condition works efficiently for nonlinear radiation problem. To avoid the numerical instability enhanced by the local gathering of grid points, the regriding technique is employed so that all the grids on the free surface may be distributed with an equal distance. This makes it possible to reduce time interval and improve numerical stability. Special attention is paid to the local flow around the body during time integration. The nonlinear radiation force is calculated by the "acceleration potential technique". Present results show good agreement with other numerical computations and experiments.

• Journal of the Microelectronics and Packaging Society
• /
• v.18 no.2
• /
• pp.11-15
• /
• 2011

Heat transfer equation is first reviewed and then governing equation of moisture diffusion. Analogy scheme is applied to analysis the moisture absorption problem of polymers. It make possible to numerically analyze the diffusion problem for single medium by using commercial finite element code if it is under the isothermal loading condition. It is extended to special multimaterial system by introducing pressure ratio function, whose moisture characteristics of materials are proportional to temperature only. The weight changes of silicon-nonconductive-polymer joint model due to moisture absorption is measured and been very close to the numerical results as for single media with boundary condition with zero concentration, but yields numerical errors as for multisystem media.

• Journal of Ocean Engineering and Technology
• /
• v.23 no.5
• /
• pp.39-44
• /
• 2009

The hydrostatic slipper bearing is generally used in high pressure axial piston pumps to support the load generated from two surfaces which are sliding relatively at low speed. The object of the bearing is to remove the possibility of direct contact by maintenance of an adequate oil film thickness between two metal surfaces. Because the bearing performance is influenced by the bearing deformation, it is highly dependent on the injection pressure, the bearing surface profile and so on. In this study, the deformation characteristics of a hydrostatic slipper bearing is investigated according to the injection pressure by the finite element analysis. In the analysis, the special boundary condition to take the fluid-structure interaction (FSI) into account is used on the interactive surface. The results, such as bearing deformation, stress and lifting force, obtained from the fully coupled analysis are compared with those from the single step sequential method.

• Special Issue of the Society of Naval Architects of Korea
• /
• 2015.09a
• /
• pp.56-62
• /
• 2015

In ship structure, many parts are in contact with inner or outer fluid as stern, ballast and oil tanks. Fatigue damages are sometimes observed in these tanks which seem to be caused by resonance with exciting force of engine and propeller. Vibration characteristics of these tanks in contact with fluid are significantly affected by fluid coupling effect. Therefore it is important to exactly predict vibration characteristics of tank structure. In order to estimate the vibration characteristics, the fluid-structure interaction(FSI) problem should be solved precisely. But it is difficult to estimate exactly the magnitude of the fluid coupling effect because it has some problems such as a fluid-structure interaction, influence by the free surface, vibration modes of structural panels and depth of water. In this paper, with fluid coupling effect, the effect of structural constraint between panels on the vibration characteristics are investigated numerically and discussed.

• Structural Engineering and Mechanics
• /
• v.6 no.6
• /
• pp.583-612
• /
• 1998

A unified third-order laminate plate theory that contains classical, first-order and third-order theories as special cases is presented. Analytical solutions using the Navier and L$\acute{e}$vy solution procedures are presented. The Navier solutions are limited to simply supported rectangular plates while the L$\acute{e}$vy solutions are restricted to rectangular plates with two parallel edges simply supported and other two edges having arbitrary combination of simply supported, clamped, and free boundary conditions. Numerical results of bending and vibration for a number of problems are discussed in the second part of the paper.