• 대한기계학회논문집B
• /
• 제21권9호
• /
• pp.1230-1243
• /
• 1997

A numerical technique for simulating incompressible viscous flow with free surface is presented. The flow field is obtained by penalty finite element formulation. In this work, a modified volume of fluid (VOF) method which is compatible with 4-node element is proposed to track the moving free surface. This scheme can be applied to irregular mesh system, and can be easily extended to three dimensional geometries. Numerical analyses were done for two benchmark examples, namely the broken dam problem and the solitary wave propagation problem. The numerical results were in close agreement with the existing data. Illustrative examples were studied to show the effectiveness of the proposed numerical scheme.

• 비파괴검사학회지
• /
• 제34권1호
• /
• pp.1-9
• /
• 2014

A two-dimensional numerical model based on the finite element method was built to simulate the wave propagation phenomena that occur during the ultrasonic time of flight diffraction (TOFD) process. First, longitudinal-wave TOFD was simulated, and the numerical results agreed well with the theoretical results. Shear-wave TOFD was also investigated because shear waves have higher intensity and resolution. The shear wave propagation was studied using three models with different boundary conditions, and the tip-diffracted shear-to-longitudinal wave was extracted from the A-scan signal difference between the cracked and non-cracked specimens. This signal showed very good agreement between the geometrical and numerical arrival times. The results of this study not only provide better understanding of the diffraction phenomena in TOFD, but also prove the potential of shear-wave TOFD for practical application.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2007년도 춘계학술대회B
• /
• pp.2947-2950
• /
• 2007

A 3D human ventricular model is proposed to simulate an integrative analysis of heart physiology and blood hemodynamics. This consists of the models of electrophysiology of human cells, electric wave propagation of tissue, heart solid mechanics, and 3D blood hemodynamics. The 3D geometry of human heart is discretized to a finite element mesh for the simulation of electric wave propagation and mechanics of heart. In cellular level, excitations by action potential are simulated using the existing human model. Then the contraction mechanics of a whole cell is incorporated to the excitation model. The excitation propagation to ventricular cells are transiently computed in the 3D cardiac tissue using a mono-domain method of electric wave propagation in cardiac tissue. Blood hemodynamics in heart is also considered and incorporated with muscle contraction. We use a PISO type finite element method to simulate the blood hemodynmaics in the human ventricular model.

• 동력기계공학회지
• /
• 제14권5호
• /
• pp.48-55
• /
• 2010

When a large ship is advancing in waves, ship undergoes the hydroelastic response, which has influences on structural stability and the fatigue destruction etc. of the ship. Therefore, to predict accurate hydroelastic response, it is necessary to analyze hydroelastic response including fluid-structure interaction. In this research, a ship is divided into many hull elements to calculate the fluid forces and wave exciting forces on each elements using three-dimensional source distribution method. The calculated fluid forces and wave exciting forces are assigned to nodes of hull elements. The neighbor nodes are connected with elastic beam elements. We analyzed hydroelastic responses, and those are formulated by using finite element method. Particularly, to estimate the influence of forward speed on the hydroelastic responses, we use two different methods : Full Hull Rotation Method(FHRM) and Sectional Hull Rotation Method(SHRM).

• 한국해안해양공학회지
• /
• 제10권4호
• /
• pp.174-186
• /
• 1998

본 연구에서는 지배방정식으로 확장형 완경사방정식을 사용하고 무한요소를 이용하여 방사조건을 처리하는 Galerkin 유한요소 모형을 수립하였다. 수립된 모형의 타당성과 적용성을 입증하기 위하여 Ippen and Goda((1963)의 완전개방 직사각형 모형항만에서의 항만 공진과 Sharp(1968) 및 Chandrasekera and Cheung(1997)의 원형 천뢰상을 전파하는 파랑 변형에 대한 수치해석을 실시하였다. 수리모형실험 및 복합요소 모형에 의한 결과와의 비교를 통하여 본 모형이 급경사 지형에도 매우 양호한 결과를 제시함을 확인하였다. 마지막으로 방파제의 대안으로 고려될 수 있는 원형 해저 우물을 설정하고 이를 지나는 파의 변형 특성을 검토하였다.

• 비파괴검사학회지
• /
• 제28권2호
• /
• pp.199-204
• /
• 2008

The generation of axisymmetric Lamb waves and interaction with wall thinning (corrosion) defects in hollow cylinders are simulated using the finite element method. Guided wave interaction with defects in cylinders is challenged by the multi-mode dispersion and the mode conversion. In this paper, two longitudinal, axisymmetric modes are generated using the concept of a time-delay periodic ring arrays (TDPRA), which makes use of the constructive/destructive interference concept to achieve the unidirectional emission and reception of guided waves. The axisymmetric scattering by the wall thinning extending in full circumference of a cylinder is studied with a two-dimensional FE simulation. The effect of wall thinning depth, axial extension, and the edge shape on the reflections of guided waves is discussed.

• 대한기계학회논문집A
• /
• 제35권8호
• /
• pp.957-964
• /
• 2011

본 연구에서는 압전 밴드 갭 구조물(포논 결정) 에 대한 체적 탄성파의 전파 특성을 주파수 및 모드 별로 파악하기 위한 유한 요소법의 적용 방안을 제안하였다. 이를 위해 체적 탄성 진행파의 면내 모드 뿐만 아니라 면외 모드를 포함하도록 3 차원 주기 경계 조건을 고려하였다. 특히, 체적 탄성파 모드 간의 비연성 특성을 전기 분극 방향에 따라 유도한 다음, 그 결과를 유한 요소 모델링에 반영하였다. 제안된 방법은 실제 시뮬레이션을 통해 다양한 형태의 압전 밴드 갭 구조물의 파동 특성 분석에 적용될 수 있는 일반적이고 효율적인 방법임을 확인하였다.

• 대한기계학회논문집A
• /
• 제25권4호
• /
• pp.574-581
• /
• 2001

Transient linear elastodynamic problems are numerically analyzed in a time-domain by the Finite Element Method, for which the variational formulation based upon the equations of motion in convolution integral is newly derived. This formulation is implicit and does not include the time derivative terms so that the computation procedure is simple and less assumptions are required comparing to the conventional time-domain dynamic numerical algorithms, being able to get the improved numerical accuracy and stability. That formulation is expanded using the semi-discrete approximation to obtain the finite element equations. In the temporal approximation, the time axis is divided equally and constant and linear time variations are assumed in those intervals. It is found that unconditionally stable numerical results are obtained in case of the constant time variation. Some numerical examples are given to show the versatility of the presented formulation.

• 한국압력기기공학회 논문집
• /
• 제14권1호
• /
• pp.38-47
• /
• 2018

For safe operation of nuclear power plants, a loose-part monitoring system (LPMS) is used to detect and locate loose-parts within the reactor coolant system, and to estimate their mass and damage potential. There are several methods to estimate mass, such as the center frequency method based on the Hertz's impact theory, a frequency ratio method and so on, but it is known that these methods cannot provide accurate information on impact response for identifying the impact source. Thanks to increasing computing power, finite element analysis (FEA) method recently become an available option to calculate reliably impact response behavior. In this paper, a finite element analysis model to simulate the propagation behavior of the bending wave, generated by a metal ball impact, is validated by performing a series of impact tests and the corresponding finite element analyses for flat plate and shell structures. Also, a FEA-based metal sphere signal map is developed, and then blind tests are performed to verify the map. This study provides an accurate simulation method for predicting the metal impact behavior and for building a metal sphere signal map, which can be used to estimate the mass of loose-parts on site in nuclear power plants.

• 한국방재학회:학술대회논문집
• /
• 한국방재학회 2007년도 정기총회 및 학술발표대회
• /
• pp.137-139
• /
• 2007

한반도의 지각에 대한 지진파의 파형들을 피킹하기 위하여 한반도의 남동부지역에서 시행된 폭파로부터 각 관측소에서 관측된 지진자료를 2 차원 수치 순산모델링으로 계산된 이론적인 신호와 비교하였다. 실제 지진자료와 같은 완전한 파형을 얻기 위해서는 실체파뿐만 아니라 표면파를 포함한 모든 파형을 이론적으로 생성하여야 한다. 또한, 각 파형의 Q 값을 고려하여 그 감쇄 정도를 파악하여야 실제 파형과 똑 같은 이론적인 파형을 만들 수 있다. 본 연구에서는 FEM(Finite Element Method)과 FDM(Finite Difference Method) 같은 전형적인 수치 모델링 기법 대신에 수도스펙트럼기법(pseudo-spectral method)을 사용하여 이론적인 파형을 계산하였다. Q 값을 고려하지 않았기 때문에 완전한 파형을 얻기에는 어려움이 있었다. 그러나 각 관측소에 도달하는 실체파의 초기 파형들의 주시를 실제 파형과 맞추었다는 데에 그 성과가 있다고 할 수 있다. 본 연구의 궁극적인 목표는 수치 모델링을 통하여 지진원의 위치와 종류를 밝혀 내는데 그 정확성을 높이는 것이다.