• Structural Engineering and Mechanics
• /
• v.63 no.4
• /
• pp.551-565
• /
• 2017

In this paper, an analytical approach is proposed for determining vibration characteristics of cracked non-uniform continuous Timoshenko beam carrying an arbitrary number of spring-mass systems. This method is based on the Timoshenko beam theory, transfer matrix method and numerical assembly method to obtain natural frequencies and mode shapes. Firstly, the beam is considered to be divided into several segments by spring-mass systems and support points, and four undetermined coefficients of vibration modal function are contained in each sub-segment. The undetermined coefficient matrices at spring-mass systems and pinned supports are obtained by using equilibrium and continuity conditions. Then, the overall matrix of undetermined coefficients for the whole vibration system is obtained by the numerical assembly technique. The natural frequencies and mode shapes of a cracked non-uniform continuous Timoshenko beam carrying an arbitrary number of spring-mass systems are obtained from the overall matrix combined with half-interval method and Runge-Kutta method. Finally, two numerical examples are used to verify the validity and reliability of this method, and the effects of cracks on the transverse vibration mode shapes and the rotational mode shapes are compared. The influences of the crack location, depth, position of spring-mass system and other parameters on natural frequencies of non-uniform continuous Timoshenko beam are discussed.

• International Journal of Aeronautical and Space Sciences
• /
• v.3 no.1
• /
• pp.74-85
• /
• 2002

The numerical experiment has been conducted to investigate the unsteady shock wave reflecting phenomena. The cell-vertex finite-volume, Roe's upwind flux difference splitting method with unstructured grid is implemented to solve unsteady Euler equations. The $4^{th}$-order Runge-Kutta method is applied for time integration. A linear reconstruction of the flux vector using the least-square method is applied to obtain the $2^{nd}$-order accuracy for the spatial derivatives. For a better resolution of the shock wave and slipline, the dynamic grid adaptation technique is adopted. The new concept of grid adaptation technique, which is much simpler than that of conventional techniques, is introduced for the current study. Three error indicators (divergence and curl of velocity, and gradient of density) are used for the grid adaptation procedure. Considering the quality of the solution and the numerical efficiency, the grid adaptation procedure was updated up to $2^{nd}$ level at every 20 time steps. For the convenience of comparison with other experimental and analytical results, the case of interaction between the straight incoming shock wave and a sharp wedge is simulated for various flow conditions. The numerical results show good agreement with other experimental and analytical results, in the shock wave reflecting structure, slipline, and the trajectory of the triple points. Some critical cases show disagreement with the analytical results, but these cases also have been proven to show hysteresis phenomena.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.24 no.1
• /
• pp.45-53
• /
• 2014

For the time integration solution of the impulsive dynamic contact problem of high-speed rotating disks formulated by the finite element technique, the velocity and acceleration contact constraints as well as the displacement contact constraint are imposed for the numerical stability without spurious oscillations. The solution of the present technique is checked by the numerical simulation using the concentric high-speed rotating disks with the clearance and impulsive loading. It is shown that the almost steady state solution agrees with the corresponding analytical solution of the elasticity and that the differentiated constraints are crucial for the numerical stability of such high-speed contact problems of the disks under impulsive loading.

• Journal of Navigation and Port Research
• /
• v.43 no.5
• /
• pp.320-327
• /
• 2019

This study examined the port expansion plan for the fishery port at the east coast of Korea, in accordance with permission conditions for coastal ports such as a limit on the cargo volume and passenger demand for the coastal tourism belt formation. The site was chosen as a municipal coastal port attracting the new ocean industry and building waterfront zone as a hub of new marine tourism. Two different numerical models (Swan and Bouss 2D) were used. Before applying to the target sea area, some numerical tests were conducted for the variation according to Bouss-2D's strong/weak and nonlinear technique compared to the irregular diffraction of semi-infinite breakwater with a theoretical solution. As a result, there was a difference in strong nonlinearity with breaking waves and it was necessary to experiment with a strong nonlinear analysis technique for the actual site. Two numerical models were applied to the fishery port site and the tranquility of some alternatives were analyzed. The numerical results show the most suitable plan was ALT-1, with satisfied harbor tranquility and reasonable economic sense. The extension of the east breakwater and enlarged turning basin of the F-Land plan have brought generally more stable harbor tranquility than the ALT-1. The result can be used as a reference for the port expansion plan in the future.

• Geomechanics and Engineering
• /
• v.16 no.5
• /
• pp.539-545
• /
• 2018

Longitudinal Displacement Profile (LDP) is an appropriate tool for determination of the displacement magnitude of the tunnel walls as a function of the distance to the tunnel face. Some useful formulations for calculation of LDP have been developed based on the monitoring data on site or by 3D numerical simulations. However, the presented equations are only based on the tunnel dimensions and for different quality of rock masses proposed a unique LDP. In the present study, it is tried to present a new formulation, for calculation of LDP, on the basis of Rock mass quality. For this purpose, a comprehensive numerical simulation program was developed to investigate the effect of rock mass quality on the LDP. Results of the numerical modelling were analyzed and the least square technique was used for fitting an appropriate curve on the derived data from the numerical simulations. The proposed formulation in the present study, is a logistic function and the constants of the logistic function were predicted by rock mass quality index (GSI). Results of this study revealed that, the LDP curves of the tunnel surrounded by rock masses with high quality (GSI>60) match together; because the rock mass deformation varies over an elastic range.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2006.11a
• /
• pp.373-376
• /
• 2006

In general, a method for generating irregular wave by combination of component waves obtained from linear wave theory is widely used. In these method, however, mean water surface elevation is rising from time to time because of nonlinear effect of wave. In this study, for the rising problem of mean water surface elevation and stabilization of calculation from time to time, mass transport velocity for horizontal velocity at wave source position is considered. The rising problem of mean water surface elevation is checked by comparing calculated wave profile from numerical technique proposed in this study with target wave profile at wave source position in numerical wave tank by using CADMAS-SURF code. And, by generating irregular wave, the validity of wave overtopping rate estimated from this numerical analysis is discussed by comparing computed results with measured results in hydraulic model experiments for vertical seawall located on a sloping sea bottom. As a results, the computations are validated against the previously experimental results by hydraulic model test and numerical results of this study and a good agreement is observed. Therefore, numerical technique of this study is a powerful tool for estimating wave overtopping rate over the crest of coastal structure.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2010.09b
• /
• pp.171-180
• /
• 2010

In this study, applicability of software used to analyze the lateral deformations and its shape of braced cuts, executed to excavate the ground for constructing the underground structures, was assessed by performing field measurements and numerical analyses and their behaviors were also compared with results of previous studies. Three typical sections, located at the construction site where the subway was on the process of construction at Suwon city in Gyeonggido, were chosen and the data of field measurement at those sections were collected. Numerical analyses with FEM technique of using the software PLAXIS and elasto-plastic approach of using the software MIDAS were performed. In general, the deformed shapes of braced cut, obtained from numerical approaches, were in relatively good agreements with results from field measurements. For sections of A-A and B-B, measured values were greater than analyzed ones whereas they were in relatively good agreements in the section C-C. As results of comparing the values from the measurements and the estimations, they were found to be close to each other so that numerical approaches were assessed to be appropriate to estimate the lateral deformation. The numerical technique with FEM was preferred to use because it estimated closer to the measurements than the elasto-plastic approach.

• Journal of computational fluids engineering
• /
• v.13 no.4
• /
• pp.96-106
• /
• 2008

We review a series of studies on turbulent skin friction drag reduction in wall-turbulence recently conducted in Japan. First, an identity equation relating the skin friction drag and the Reynolds shearstress (the FIK identity) is introduced. Based on the implication of the FIK identity, a new analytical suboptimal feedback control law requiring the streamwise wall-shear stress only is introduced and direct numerical simulation (DNS) results of turbulent pipe flow with that control is reported. We also introduce DNS of an anisotropic compliant surface and parameter optimization using an evolutionary optimization technique.

• Structural Engineering and Mechanics
• /
• v.4 no.2
• /
• pp.177-190
• /
• 1996

A method to numerically evaluate the behaviour of single span beams, prestressed with external tendons and symmetrically loaded is presented. This algorithm, based on the Finite Difference Method, includes second order effects and large displacements in an attempt to more fully understand the behaviour of the beam up to collapse. The numerical technique discussed is particularly appropriate for the analysis of R.C. and P.C. beams rehabilitated or strengthened by means of external prestressing but it is reliable for the analysis of new beams as well.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03a
• /
• pp.347-354
• /
• 2008

We introduce a series of studies on turbulent skin friction drag reduction in wall-turbulence. First, an identity equation relating the skin friction drag and the Reynolds shear stress (the FIK identity) is introduced. Based on the implication of the FIK identity, a new analytical suboptimal feedback control law requiring the streamwise wall-shear stress only is introduced and direct numerical simulation (DNS) results of turbulent pipe flow with that control is reported. We also introduce DNS of an anisotropic compliant surface and parameter optimization using an evolutionary optimization technique.