• 한국전산구조공학회논문집
• /
• 제31권2호
• /
• pp.87-95
• /
• 2018

이 논문에서는 스펙트럴 요소법과 외연적 시간적분법을 이용해 SH파의 전파 거동을 계산하는 수치해석 기법을 제시한다. 2차원 영역에서의 탄성파 해석을 위해 해석영역을 유한 영역으로 한정하고 파동이 반사되지 않도록 수치적 파동흡수 경계조건인 perfectly matched layer(PML)를 도입하였다. PML이 포함된 시간영역 파동방정식의 유한요소해법을 위해 스펙트럴 요소법을 적용하였고 Legendre- Gauss-Lobatto 수치적분법을 사용하여 질량행렬을 대각화하였다. 2차 미분방정식 시스템의 파동방정식을 1차 미분방정식 시스템으로 변환하였고 병렬화를 통한 탄성파 해석 성능의 최적화를 위해 외연적 시간적분법인 4차 Runge-Kutta 방법을 이용해 해석영역에서의 변위응답을 계산하였다. 2차원 해석영역에서 SH파의 전파 거동을 계산하는 수치예제를 통해 제시한 외연적 스펙트럴 요소법의 정확성을 검증하였고 PML로 인한 반사파의 감쇠효과를 확인하였다. 외연적 시간적분법을 통한 탄성파 해석 기법은 3차원 영역과 같은 대규모 문제에서의 탄성파 수치해석을 효율적으로 수행할 수 있을 것으로 기대된다.

• KCI Concrete Journal
• /
• 제13권2호
• /
• pp.10-18
• /
• 2001

The spectral analysis of surface waves (SASW) method, which is an in-situ seismic technique, has mainly been developed and used for many years to determine the stiffness profile of layered media (such as asphalt concrete and layered soils) in an infinite half-space. This paper presents a modified experimental technique for nondestructive evaluation of in-place cement mortar compressive strength in single-layer concrete slabs of rather a finite thickness through a correlation to surface wave velocity. This correlation can be used in the quality control of early age cement mortar structures and in evaluating the integrity of structural members where the infinite half space condition is not met. In the proposed SASW field test, the surface of the structural members is subjected to an impact, using a 12 mm steel ball, to generate surface wave energy at various frequencies. Two accelerometer receivers detect the energy transmitted through the medium. By digitizing the analog receiver outputs, and recording the signals for spectral analysis, surface wave velocities can be identified. Modifications to the SASW method includes the reduction of boundary reflections as adopted on the surface waves before the point where the reflected compression waves reach the receivers. In this study, the correlation between the surface wave velocity and the compressive strength of cement mortar is developed using one 36"x36"x4"(91.44$\times$91.44$\times$91.44 cm) cement mortar slab of 2,000 psi (140 kgf/$\textrm{cm}^2$) and two 36"x36"x4"(91.44$\times$91.44$\times$91.44 cm) cement mortar slabs of 3,000 psi (210 kgf/$\textrm{cm}^2$).

• Structural Engineering and Mechanics
• /
• 제14권6호
• /
• pp.625-647
• /
• 2002

Parallel execution of computational mechanics codes requires efficient mesh-partitioning techniques. These mesh-partitioning techniques divide the mesh into specified number of submeshes of approximately the same size and at the same time, minimise the interface nodes of the submeshes. This paper describes a new mesh partitioning technique, employing Genetic Algorithms. The proposed algorithm operates on the deduced graph (dual or nodal graph) of the given finite element mesh rather than directly on the mesh itself. The algorithm works by first constructing a coarse graph approximation using an automatic graph coarsening method. The coarse graph is partitioned and the results are interpolated onto the original graph to initialise an optimisation of the graph partition problem. In practice, hierarchy of (usually more than two) graphs are used to obtain the final graph partition. The proposed partitioning algorithm is applied to graphs derived from unstructured finite element meshes describing practical engineering problems and also several example graphs related to finite element meshes given in the literature. The test results indicate that the proposed GA based graph partitioning algorithm generates high quality partitions and are superior to spectral and multilevel graph partitioning algorithms.

• 전자공학회논문지B
• /
• 제29B권2호
• /
• pp.18-26
• /
• 1992

MR(Model Reduction) is presented in order to estimate the frequency of multiple sinusoids from the finite noisy data with the white or colored noises. MR, using the reduced rank models, is designed, appling the approximation of linear system to LP(Linear Prediction). The MR method is analyzed. Monte-carlo simulations are conducted for MR and Lp. The results are compared with in terms of mean, root-mean square and relative bias. MR eliminates effectevely the extremeous and exceptional poles appearing in LP and improves the accuracy of LP. Especially, MR gives promising results in short noisy measurements, low SNR's and colored noises. Power spectral density and angular frequency position are showed by figures, for examples. Finally, the new method is utilized to the communication and biomedical systems estimating the characteristics of the signal and the system identification modelling the dynamic systems from experimental data.

• 한국지진공학회:학술대회논문집
• /
• 한국지진공학회 2001년도 춘계학술대회 논문집
• /
• pp.108-119
• /
• 2001

The Spectral Analysis of surface Waves(SASW) Method has a great has a great potential for rapid determination of shear wave velocity profile of ground. However, it has an inherent limitation in the interpretation of test results due to the assumption that the ground is layered horizontally. The reason of the assumption is that difficulties exist in obtaining analytical solutions of wave equation when a soil system is composed of inclined soil layer. In this study, a finite-element method has been employed to assess the effects of dip angle and stiffness contrast of inclined soil layers and the testing direction on the dispersion curve. The propagation of wave front in the inclined soil layer was also investigated. The results indicated that the influence of dip angle on the dispersion curve is getting obvious as the dip angle increases and the propagation of wave front in the inclined layer also entirely different compared with the case of the horizontal layer.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2004년도 춘계학술대회논문집
• /
• pp.681-686
• /
• 2004

Identification of asymmetry and anisotropy of rotor system is important for diagnosis of rotating machinery. Directional frequency response functions (dFRFs) are known to be powerful tool in effectively detecting the presence of asymmetry or anisotropy. In this paper, an estimation method of dFRFs for rotors is newly developed, when both asymmetry and anisotropy are present. The method transforms the finite degrees-of-freedom time-varying linear differential equation of motion to an infinite degree-of-freedom time-invariant linear one, employing the modulated coordinates. The validity of the method is demonstrated by numerical simulation with a simple rotor model.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2001년도 추계 학술대회논문집
• /
• pp.84-90
• /
• 2001

Secondary instability of flow past a circular cylinder is examined using direct numerical simulation at Reynolds number 220 and 250. The higher-order finite difference scheme is employed for the spatial distributions along with the second order Adams-Bashforth and the first order backward-Euler time integration. In x-y plane, the convection term is applied by the 5th order upwind scheme, and the pressure and viscosity terms are applied by the 4th order central difference. In spanwise, Navier-Stokes equation is distributed using Spectral Method. The critical Reynolds number for this instability is found to be about Re=190. The secondary instability leads re three-dimensionality with a spanwise wavelength about 4 cylinder diameters at onset (A-mode). Results of three-dimensional effect in wake of a circular cylinder are represented with spanwise and streamwise vorticity contours as Reynolds numbers.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2001년도 추계학술대회논문집B
• /
• pp.570-577
• /
• 2001

Three-dimensional time-dependent flow past a circular cylinder is numerically investigated using direct numerical simulation for Reynolds number 280 and 300. The higher-order finite difference scheme is employed for the spatial distributions along with the second order Adams-Bashforth and the first order backward-Euler time integration. In x-y plane, the convection term is applied by the 5th order upwind scheme and the pressure and viscosity terms are applied by the 4th order central difference. And in spanwise, Navier-Stokes equation is distributed using of Spectral Method. At Reynolds number 259 the two-dimensional wake becomes linearly unstable to a second branch of modes with wavelength about 1.0 diameters at onset (B-mode). Present results of three-dimensional effects of in wake of a circular cylinder is represented with spanwise and streamwise vorticity contours as Reynolds numbers.

• 한국전자파학회논문지
• /
• 제21권7호
• /
• pp.805-814
• /
• 2010

기존 Crank-Nicolson FDTD 기법(CN FDTD 기법)의 비등방성 분산 특성을 개선하기 위한 CN ID-FDTD 기법을 제안하였다. 제안한 CN ID-FDTD 기법은 공간 미분 연산을 위해 기존 CN FDTD 기법의 centered 유한 차분식 (Finite Difference equation: FD 연산식)이 아닌 isotropic-dispersion 유한 차분식(ID-FD 연산식)$^{[1],[2]}$을 이용한다. 본 논문에서는 손실 매질에 대한 CN ID-FDTD 기법의 분산 관계식을 유도하였고, 이 분산 관계식을 이용해 ID-FD 연산식에서 분산 오차(dispersion error)를 줄이는 가중치(weighting factor)와 보정값(scaling factor)을 제시하였다. 그리고 해석 결과의 정확성 비교를 통해 CN ID-FDTD 기법에서는 기존 CN FDTD 기법의 단점이었던 비등방성 분산 오차가 확연하게 감소하는 것을 확인하였다.

• Smart Structures and Systems
• /
• 제21권5호
• /
• pp.591-600
• /
• 2018

The probabilistic characterization of wind field characteristics is a significant task for fatigue reliability assessment of long-span railway bridges in wind-prone regions. In consideration of the effect of wind direction, the stochastic properties of wind field should be represented by a bivariate statistical model of wind speed and direction. This paper presents the construction of the bivariate model of wind speed and direction at the site of a railway arch bridge by use of the long-term structural health monitoring (SHM) data. The wind characteristics are derived by analyzing the real-time wind monitoring data, such as the mean wind speed and direction, turbulence intensity, turbulence integral scale, and power spectral density. A sequential quadratic programming (SQP) algorithm-based finite mixture modeling method is proposed to formulate the joint distribution model of wind speed and direction. For the probability density function (PDF) of wind speed, a double-parameter Weibull distribution function is utilized, and a von Mises distribution function is applied to represent the PDF of wind direction. The SQP algorithm with multi-start points is used to estimate the parameters in the bivariate model, namely Weibull-von Mises mixture model. One-year wind monitoring data are selected to validate the effectiveness of the proposed modeling method. The optimal model is jointly evaluated by the Bayesian information criterion (BIC) and coefficient of determination, $R^2$. The obtained results indicate that the proposed SQP algorithm-based finite mixture modeling method can effectively establish the bivariate model of wind speed and direction. The established bivariate model of wind speed and direction will facilitate the wind-induced fatigue reliability assessment of long-span bridges.