• Bulletin of the Society of Naval Architects of Korea
• /
• v.24 no.3
• /
• pp.1-8
• /
• 1987

In this paper, the localized finite element method(LFEM) is applied to 3-dimensional ship motion problems in water of infinite depth. The LFEM used here is based on the functional constructed by Bai & Yeung(1974). To test the present numerical scheme, a few vertical axisymmetric bodies are treated by general 3-dimensional formulation. The computed results of hydrodynamic coefficients for a few vertical spheroids and vertical circular cylinders show good agreement with results obtained by others. The advantages of the present numerical method compared with the method of integral equation are as follows; (i) The cumbersome existence of irregular frequencies in the method of conventional integral equation is removed. (ii) The final matrix is banded and symmetric and the computation of the matrix elements is comparatively easier, whereas the size of the matrix in the present scheme is much larger. (iii) In the future research, it is possible to accommodate with the nonlinear exact free surface boundary condition in the localized finite element subdomain, whereas the linear solution is assumed in the truncated(far field) subdomain.

• Structural Engineering and Mechanics
• /
• v.2 no.3
• /
• pp.245-256
• /
• 1994

A new axisymmetric crack model is proposed on the basis of p-version of the finite element method limited to theory of small scale yielding. To this end, axisymmetric stress element is formulated by integrals of Legendre polynomial which has hierarchical nature and orthogonality relationship. The virtual crack extension method has been adopted to calculate the stress intensity factors for 3-D axisymmetric cracked bodies where the potential energy change as a function of position along the crack front is calculated. The sensitivity with respect to the aspect ratio and Poisson locking has been tested to ascertain the robustness of p-version axisymmetric element. Also, the limit value that is an exact solution obtained by FEM when degree of freedom is infinite can be estimated using the extrapolation equation based on error prediction in energy norm. Numerical examples of thick-walled cylinder, axisymmetric crack in a round bar and internal part-thorough cracked pipes are tested with high precision.

• Proceedings of the KSME Conference
• /
• 2000.04a
• /
• pp.176-181
• /
• 2000

Stress intensity factors for a rigid line inclusion tying along a bimaterial interface are calculated by the boundary element method with the multiregion and double-Point techniques. The formula between the stress intensity factors and the inclusion surface stresses are derived. The numerical values of the stress intensity factors for the bimaterial interface rigid line inclusion in the infinite body are proved to be in good agreement within 3% when compared with the previous exact solutions. In the finite bimaterial systems, the stress intensity factors for the center and edge rigid line inclusions at interface are computed with the variation of the rigid line inclusion length and the shear modulus ratio under the biaxial and uniaxial loading conditions.

• Proceedings of the KSME Conference
• /
• 2000.11a
• /
• pp.486-491
• /
• 2000

In this study, finite elements addition and removal method by stress range is applied to optimize shapes in structures, without using classical and numerical optimization methods and search methods. The program based on this algorithm is developed and compared to theoritial results with considerable accuracy. Classical methods need mesh generation for finite element analysis for every iteration, the developed method needs updated mesh data such as coordinates of nodes, elements connectivity, and loads on nodes. And other tools of finite element analysis can be in use as a black box to interface with this program.

• Proceedings of the Korean Society of Machine Tool Engineers Conference
• /
• 2000.10a
• /
• pp.434-439
• /
• 2000

In this paper, in order to study the geometric factor effect of a circular hole near a crack tip in a semi-infinite plate, the Dimensionless Stress Intensity Factor, $F(=\frac K {\sigma {\sqrt{\pi a}}})$ is analyzed at the crack tip using a two Dimensional Boundary Element Method (BEM) program which is known as superior in Fracture Mechanics. Kelvin's solution was used as a fundamental solution in BEM analysis and displacement extrapolation method was used to determine Stress Intensity Factor.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.20 no.4
• /
• pp.309-332
• /
• 2016

Finite element method has been applied to solve the fundamental governing equations of natural convective, electrically conducting, incompressible fluid flow past an infinite vertical plate surrounded by porous medium in presence of thermal radiation, viscous dissipation, Soret and Dufour effects. In this research work, the results of coupled partial differential equations are found numerically by applying finite element technique. The sway of significant parameters such as Soret number, Dufour number, Grashof number for heat and mass transfer, Magnetic field parameter, Thermal radiation parameter, Permeability parameter on velocity, temperature and concentration evaluations in the boundary layer region are examined in detail and the results are shown in graphically. Furthermore, the effect of these parameters on local skin friction coefficient, local Nusselt number and Sherwood numbers is also investigated. A very good agreement is noticed between the present results and previous published works in some limiting cases.

• Structural Engineering and Mechanics
• /
• v.14 no.1
• /
• pp.71-84
• /
• 2002

A finite element formulation for the time-domain analysis of linear transient elastodynamic problems is presented based on the weak form obtained by applying the Galerkin's method to the integro-differential equations which contain the initial conditions implicitly and does not include the inertia terms. The weak form is extended temporally under the assumptions of the constant and linear time variations of field variables, since the time-stepping algorithms such as the Newmark method and the Wilson ${\theta}$-method are not necessary, obtaining two kinds of implicit finite element equations which are tested for numerical accuracy and convergency. Three classical examples having finite and infinite domains are solved and numerical results are compared with the other analytical and numerical solutions to show the versatility and accuracy of the presented formulation.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.17 no.1
• /
• pp.59-65
• /
• 2004

The stress intensity factors have been widely used in numerical studies of crack growth direction. However in many cases, omissive terms of the series expansion are quantitatively significant, so we consider the computation of such terms. For this purpose, we used the finite element method with isoparametric quadratic quarter-point elements. For examples, infinite square plate with a slant crack subjected to a uniaxial load is analyzed. The numerical analysis were performed for the wide range of crack tip element lengths and inclined angles. The numerical results obtained are compared with the theoretical solutions. Also they were accurate and efficient.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.11
• /
• pp.2971-2980
• /
• 1995

A numerical study of fluid flow in driven cavity was carried out using singular finite element method. The driven cavity problem is known to have infinite velocity gradients as well as dual velocity conditions at the singular points. To overcome such difficulties, a finite element method with singular shape functions was used and a special technique was employed to allow multiple values of velocities at the singular points. Application of singular elements in the driven cavity problem has a significant influence on the stability of solution. It was found the singular elements gave a stable solution, especially, for the pressure distribution of the entire flow field by keeping up a large pressure at the singular points. In the existing solutions of driven cavity problem, most efforts were focused on the study of streamlines and vorticities, and pressure were seldom mentioned. In this study, however, more attention was given to the pressure distribution. Computations showed that pressure decreased very rapidly as the distance from the singular point increased. Also, the pressure distribution along the vertical walls showed a smoother transition with singular elements compared to those of conventional method. At the singular point toward the flow direction showed more pressure increase compared with the other side as Reynolds number increased.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2008.04a
• /
• pp.72-77
• /
• 2008

Stability is a very important part which we must consider in structural design. In this paper, we take advantage of finite element method, and study about parametrical instability of star-dome structures, which is subjected to harmonically pulsating load. When calculating stiffness matrix, we consider elastic stiffness and geometrical stiffness simultaneously. In equation of motion, we represent displacements and accelerations by trigonometric series expansions, and then obtain Hill's infinite determinants. After first order approximation, we can get first and second order dynamic instability region finally.