• Structural Engineering and Mechanics
• /
• v.81 no.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.

• Journal of Power System Engineering
• /
• v.15 no.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.

• Journal of Ocean Engineering and Technology
• /
• v.23 no.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.

• Bulletin of the Korean Mathematical Society
• /
• v.51 no.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.

• Computational Structural Engineering
• /
• v.6 no.4
• /
• pp.89-98
• /
• 1993

As the construction of nuclear power plants are increased, nuclear wastes disposal has been faced as a serious problem. If nuclear wastes are to be buried in the underground stratum, thermo-mechanical behavior of stratum must be analyzed, because high temperature distribution has a significant effect on tunnel and surrounding stratum. In this study, in order to analyze the structural behavior of the underground which is subject to concentrated heat sources, a coupling method of nonlinear finite elements and linear boundary elements is proposed. The nonlinear finite elements (NFE) are applied in the vicinity of nuclear depository where thermo-mechanical stress is concentrated. The boundary elements are also used in infinite domain where linear behavior is expected. Using the similar method as for the problem in mechanical field, the coupled nonlinear finite element-boundary element (NFEBE) is developed. It is found that NFEBE method is more efficient than NFE which considers nonlinearity in the whole domain for the nuclear wastes depository that is expected to exhibit local nonlinearity behavior. The effect of coefficients of the rock mass such as Poisson's ratio, elastic modulus, thermal diffusivity and thermal expansion coefficient is investigated through the developed method. As a result, it is revealed that the displacements around tunnel are largely dependent on the thermal expansion coefficients.

• KSTLE International Journal
• /
• v.5 no.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.

• Proceedings of the Korean Society of Machine Tool Engineers Conference
• /
• 1997.04a
• /
• 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.

• Journal of Ocean Engineering and Technology
• /
• v.7 no.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.

• The Journal of the Acoustical Society of Korea
• /
• v.31 no.1
• /
• pp.19-28
• /
• 2012

The time-domain solution from the Kirchhoff integral equation for an exterior problem is not unique at certain eigen-frequencies associated with the fictitious internal modes as happening in frequency-domain analysis. One of the solution methods is the CHIEF (Combined Helmholtz Integral Equation Formulation) approach, which is based on employing additional zero-pressure constraints at some interior points inside the body. Although this method has been widely used in frequency-domain boundary element method due to its simplicity, it was not used in time-domain analysis. In this work, the CHIEF approach is formulated appropriately for time-domain acoustic boundary element method by constraining the unknown surface pressure distribution at the current time, which was obtained by setting the pressure at the interior point to be zero considering the shortest retarded time between boundary nodes and interior point. Sound radiation of a pulsating sphere was used as a test example. By applying the CHIEF method, the low-order fictitious modes could be damped down satisfactorily, thus solving the non-uniqueness problem. However, it was observed that the instability due to high-order fictitious modes, which were beyond the effective frequency, was increased.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2002.04a
• /
• 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.