• Structural Engineering and Mechanics
• /
• 제81권2호
• /
• pp.179-194
• /
• 2022

The bending, buckling and free vibration responses of functionally graded material (FGM) beams are investigated semi-analytically by the scaled boundary finite element method (SBFEM) in this paper. In the concepts of the SBFEM, the dimension of computational domain can be reduced by one, therefore only the axial dimension of the beam is discretized using the higher order spectral element, which reduces the amount of calculation and greatly improves the calculation efficiency. The governing equation of FGM beams is derived in detail by the means of the principle of virtual work. Compared with the higher-order beam theory, fewer parameters and simpler control equations are used. And the governing equation is transformed into a first-order ordinary differential equation by introducing intermediate variables. Analytical solutions of the governing equation can be obtained by pade series expansion in the direction of thickness. Numerical example are compared with the numerical solutions provided by the previous researchers to verify the accuracy and applicability of the proposed method. The results show that the proposed formulations can quickly converge to the reference solutions by increasing the order of higher order spectral elements, and high accuracy can be achieved by using a small number of the elements. In addition, the influence of the structural sizes, material properties and boundary conditions on the mechanical behaviors of FG beams subjected to different load types is discussed.

• 동력기계공학회지
• /
• 제15권6호
• /
• pp.27-34
• /
• 2011

The increasing capabilities of the computers enable us to utilize various numerical schemes for the time-domain simulations concerned with 3-dimensional free-surface wave problems. There are still difficulties to solve such kind of problems, however. That's because long time simulations with large computational domain are needed in time-domain analysis. So, we need faster and more efficient numerical schemes to get the solutions practically for these problems. In this paper, a high-order spectral/boundary-element method is used for the numerical investigation of physics involved in wave-body interaction. This method is one of the most efficient numerical methods by which the nonlinear gravity waves can be simulated and hydrodynamic forces also can be calculated in time-domain. To get the robust study in these topics, various numerical tests are performed and compared with others' works.

• 한국해양공학회지
• /
• 제23권6호
• /
• pp.12-16
• /
• 2009

In this study, the wave profiles and kinematics of highly nonlinear waves at various water depths were calculated using a 2D fully nonlinear Numerical Wave Tank (NWT). The NWT was developed based on the Boundary Element Method (BEM) with the potential theory and the mixed Eulerian-Lagrangian (MEL) time marching scheme by 4th-order Runge-Kutta time integration. The spatial variation of intermediate-depth waves along the direction of wave propagation was caused by the unintended generation of 2nd-order free waves, which were originally investigated both theoretically and experimentally by Goda (1998). These free waves were induced by the mismatch between the linear motion of wave maker and nonlinear displacement of water particles adjacent to the maker. When the 2nd-order wave maker motion was applied, the spatial modulation of the waves caused by the free waves was not observed. The respective magnitudes of the nonlinear wave components for various water depths were compared. It was found that the high-order wave components greatly increase as the water depth decreases. The wave kinematics at various locations were calculated and compared with the linear and the Stokes 2nd-order theories.

• 대한수학회보
• /
• 제51권2호
• /
• pp.595-612
• /
• 2014

The spectral collocation method for a second order elliptic boundary value problem on a domain ${\Omega}$ with curved boundaries is studied using the Gordon and Hall transformation which enables us to have a transformed elliptic problem and a square domain S = [0, h] ${\times}$ [0, h], h > 0. The preconditioned system of the spectral collocation approximation based on Legendre-Gauss-Lobatto points by the matrix based on piecewise bilinear finite element discretizations is shown to have the high order accuracy of convergence and the efficiency of the finite element preconditioner.

• 전산구조공학
• /
• 제6권4호
• /
• pp.89-98
• /
• 1993

최근 원자력의 사용이 증가함에 따라 핵폐기물을 효과적으로 처리하는 문제에 관심이 집중되고 있다. 이러한 핵폐기물을 지층내에 저장할 경우 고온의 열에 의해 핵폐기물 구조체에 지대한 영향을 미치므로 지반의 열력학적 거동을 분석할 필요성이 요구된다. 본 연구는 지반내에 처분된 고온의 사용후 핵연료에 의한 열역학적인 응력이 집중되어 비선형 거동이 예상되는 저장구조체 주변에는 비선형 유한요소를 적용하고 선형거동이 예상되는 무한영역에는 선형경계요소를 사용하여, 일반적인 역학적 계와 동일한 방법으로 비선형 유한요소와 경계요소를 조합한 프로그램을 개발하였다. 사용후 핵연료 폐기구조체와 같이 국부적인 비선형거동이 예상되는 구조물에서는 조합방법이 전 영역을 비선형 유한요소로 모형화하여 해석하는 것보다 효율적임을 알 수 있었다. 또한, 지층내 지반에 영향 미치는 주요 지반계수를 변화시킨 경우, 터널경계의 변위에 이러한 계수들이 어떠한 영향을 미치는가를 개발된 방법을 사용하여 검토하였다. 검토결과, 다른 계수들의 변화보다 열팽창계수의 변화가 터널주위의 변위에 상당한 영향을 미침을 알 수 있었다.

• KSTLE International Journal
• /
• 제5권2호
• /
• pp.52-56
• /
• 2004

The hard disk drive performance depends strongly on air bearing characterisitcs of the head slider. The objective of the slider design is to provide accurate positioning of the magnetic read/write element at the very small height above the disk. Application of the numerical methods is required due to complexity of the air bearing surface of the slider. The Boundary-Fitted Coordinate System Divergence Formulation method can be used for calculation of pressure distribution in the case of steep film thickness gradients. In the present work, the interpolating functions used in the expression for the Couette flow are modified in order to improve the solution characteristics in the extremely high compressibility number region. The advantages of the modified method are demonstrated on example of the flat skewed slider. Finally, the modi.ed method is applied to analysis of the static characteristics of the femto-slider. The analysis results indicate the effect of the silder's air bearing surface crown on the flying height and the pitching angle in steady state position.

• 한국공작기계학회:학술대회논문집
• /
• 한국공작기계학회 1997년도 춘계학술대회 논문집
• /
• pp.156-163
• /
• 1997

Recently, attempts have been made to be join ceramics to metals in order to make up for the brittleness of ceramics. The difference in the coefficients of linear expansion of the two materials joined at high temperature will cause residual stress, which has a strong influence on the strength of the bonded joints. In this paper, the residual stress distribution and stress intensity factors of the ceramic/metal bonded joints were analyzed by 2-dimensional elastic boundary element method. Fracture toughness tests of ceramic/metal bonded joints with an interface crack were carried out. So the advanced method of quantitative strength evaluation for ceramic/metal bonded joints is to be suggested. Fracture surface and crack propagation path were observed using scanning electron microscope.

• 한국해양공학회지
• /
• 제7권2호
• /
• pp.9-18
• /
• 1993

Numerical wave tank is a robust tool by which the nonlinear interactions between the body and the free-surface can be treated in time-domain. In this paper, a two-dimensional numerical wave tank based on the Spectral/Boundary-Element Method is developed, and applied successfully to the study of nonlinear wave diffraction around a submerged circular cylinder. Particularly, it is shown that the high-order wave components of significant wave height are developed in the lee-side of the cylinder and that these waves result in a negative drift force on the circular cylider.

• 한국음향학회지
• /
• 제31권1호
• /
• pp.19-28
• /
• 2012

Kirchhoff 적분식을 이용하여 외부 음향 문제의 시간 영역 응답을 계산하는 경우, 주파수영역 해석과 마찬가지로 가상적인 내부 음향 모드에 기인한 비유일성 문제가 발생한다. 이를 해결하는 방법들 중의 하나로서 CHIEF(Combined Helmholtz Integral Equation Formulation) 방법이 쓰이는데, 이는 몇몇 내부 수음점의 응답을 0으로 추가하여 구속하는 조건을 부가하는 기법이다. 이 기법은 주파수 영역 경계요소법에서는 간편한 수식 때문에 많이 사용되고 있지만, 시간 영역에서는 사용된 예가 없다. 본 연구에서는 대상체 내부의 가상 수음점과 경계 표면의 절점들간의 최소 거리에 대한 지연시간을 고려하여, 계산하고자 하는 미지수인 현재 시간의 경계 표면 음장을 구속함으로써, 시간 영역 해석에 적합하도록 CHIEF 방법을 수식화하였다. 예제로서, 반지름 방향으로 진동하는 구의 음향 방사 문제를 다루었다. CHIEF 방법을 적용함에 따라 저차의 내부 음향 모드에 기인한 비유일성 문제를 해결할 수 있었고, 비요동 모드에 의한 수치적 불안정성을 피할 수 있었다. 그러나, 유효주파수 밖에 남은 내부 음향의 고차모드들에 의한 수치적 불안정성은 증가하였다.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2002년도 봄 학술발표회 논문집
• /
• pp.12-18
• /
• 2002

For the analysis of Kirchhoff plate bending problems, a new meshless method is implemented. For the satisfaction of the C¹ continuity condition in which the first derivative is treated as another primary variable, Hermite interpolation is enforced on standard reproducing kernel particle method. In order to impose essential boundary conditions on solving C¹ continuity problems, shape function modifications are adopted. Through numerical tests, the characteristics and accuracy of the HRKPM are investigated and compared with the finite element analysis. By this implementation, it is shown that high accuracy is achieved by using HRKPM fur solving Kirchhoff plate bending problems.