• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.32 no.3
• /
• pp.183-190
• /
• 2019

This paper introduces a near-tip grid refinement and explores its usefulness in the crack analysis by the natural element method(NEM). As a sort of local h-refinement in finite element method(FEM), a NEM grid is locally refined around the crack tip showing high stress singularity. This local grid refinement is completed in two steps in which grid points are added and Delaunay triangles sharing the crack tip node are divided. A plane strain rectangular plate with symmetric edge cracks is simulated to validate the proposed local grid refinement and to examine its usefulness in the crack analysis. The crack analysis is also simulated using a uniform NEM grid for comparison. Unlike the uniform grid, the refined grid provides near-tip stress distributions similar to the analytic solutions and the fine grid. In addition, the refined grid shows higher convergence than the uniform grid, the global relative error to the total number of grid points.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2006.05a
• /
• pp.588-593
• /
• 2006

자연하천은 일반적으로 만곡수로나 사행수로 형태를 보이고 있으며, 직선수로에서와 달리 원심력에 기인한 이차류 영향을 받게 된다. 이차류에 의해서 수면에서는 만곡부 바깥쪽으로, 하상에서는 만곡부 안쪽으로의 흐름특성을 보이게 된다. 만곡부 안쪽으로 가해지는 하상 전단응력에 기인하여 하상에서의 입자가 만곡부 안쪽으로 이송되며, 결과적으로 만곡부 안쪽에는 점사주가, 바깥쪽에는 소(pool)가 생성된다. 또한 지형경사의 생성으로 입자에 가해지는 중력효과도 변화된다. 따라서 이와 같은 자연하천의 흐름과 하상변동을 수치모의 하기 위해서는 만곡부 이차류 특성을 고려한 모형이 필요하다. 본 연구에서는 수심 적분된 흐름방정식과 하상토 보존방정식 (Exner equation)을 이용한 하상변동을 위한 비연계 수치모형을 위해서 하상토 보존방정식의 유한요소 알고리즘을 개발하였다. 하상토 보존방정식은 흐름 특성에 따른 평형 유사량의 공간변화율을 이용하여 일정 기간 동안의 하상 변화량을 계산한다. 이 때 이차류에 의한 하상 전단응력의 편각 및 지형경사 변화에 따른 실제 입자의 이송방향을 보정하여 평형 유사량이 계산된다. 이러한 보정식을 적용시키기 위해서는 유속성분의 공간변화량 및 지형경사의 공간성분이 필요하다. 유한요소법은 연속성 변수를 이산화시켜 근사해를 구하는 수치기법의 일종이기 때문에, 요소망이 불규칙적으로 구성되었을 경우 임의의 절점에서 연속성을 지닌 변수의 공간변화율을 계산하는데 어려움이 있다. 따라서 본 연구에서는 평형 유사량 계산 시에 절점이 아닌 요소 내부에서 평형 유사량을 계산하는, 하상토 보존방정식의 새로운 유한요소 알고리즘을 개발하고, 새로운 알고리즘을 적용시킨 수치모형의 검증을 행하였다. 경계조건 알고리즘의 검증으로 위해서 Soni 등 (1980)이 행한 상류 유입 유사량에 따른 하상변동을 수치 모의하고 실험치와 비교하였으며, Sutmuller와 Glerum (1980)이 수행한 만곡수로에서의 하상변동을 모의하고 실험과 비교하였다. 새로운 알고리즘을 적용시킨 하상토 보존방정식의 유한요소 수치모형의 결과는 매우 안정적이며, 실험과 매우 유사한 결과를 얻을 수 있었다. 본 수치모델은 현재 균일한 입자의 하상토만을 고려하므로, 입자분급이나 하상 장갑화 현상 등은 무시한다.

• Computational Structural Engineering
• /
• v.6 no.3
• /
• pp.125-137
• /
• 1993

Recently, many researches are being dealt with the adaptive method for improving the accuracy of finite element solution. This paper deals with rh-methods that are the combination of r and h-method ; r-method is to relocate the nodes for the grid optimization, h-method is to divide the elements with great error into the equal shape. As a results, rh-method has the same error decrease and convergence as h-method in the same degree of freedom, but it has more exact result of finite element in the state of restraining degrees of freedom than h-method alone.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.28 no.1
• /
• pp.53-61
• /
• 2015

This paper presented the nonlinear behaviors of the single-layered lattice dome, which is widely used for the long-span structure system. The behaviors were analysed through the classical shell buckling theory as the single-layered lattice dome behaves like continum thin shell due to its geometric characteristics, and finite element analysis method using the software program Nastran. Shell buckling theory provides two types of buckling loads, the global- and member buckling, and finite element analysis provides the ultimate load of geometric nonlinear analysis as well as the buckling load of Eigen value solution. Two types of models for the lattice dome were analysed, that is rigid- and pin-jointed structure. Buckling load using the shell buckling theory for each type of lattice dome, governed by the minimum value of global buckling or member buckling load, resulted better estimation than the buckling load with Eigen value analysis. And it is useful to predict the buckling pattern, that is global buckling or member buckling.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.16 no.2
• /
• pp.207-216
• /
• 2003

In this paper, the finite element method for the elastodynamic axisymmetric fracture analysis is presented in matrix form through the application of the Galerkin method to the time integral equations of motion with no inertia forces. Isoparametric quadratic quadrilateral element and triangular crack tip singular elements with one-quarter node are used in the mesh division of the finite element model. To show the validity and accuracy of the proposed method, the infinite elastic medium with the penny shaped crack is solved first and compared with the analytical solution and the numerical results by the finite difference method and the boundary element method existing in the published literatures, and then the dynamic stress intensity factors of solid and hollow cylinders of finite dimensions haying penny-shaped cracks and internal and external circumferential tracks are computed in detail.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.12 no.1
• /
• pp.1-8
• /
• 1992

A method analizing contact pressure between plate and elastic half space is presented by using F.E.M. With the method, the pressure intensities at surface nodes of half space cae be directly calculated by using flexibility matrix of half space. The method is originally presented by Y.K. Cheung et al.(3) Insted of Y.K. Cheung's method, which use a conception of equi-contact pressure area around each surface nodes of half space in the noded rectanqular element area. We use the equi-contact pressure area around the Gaussian integration points of half space surface in the noded isoparametric element area. Numarical examples are presented and compared with other's studies.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.17 no.3
• /
• pp.279-293
• /
• 2004

The analysis of linear static and free vibration problems of isotropic and laminated composite plates and shells is performed by the improved 9-node shell element with the new strain displacement relationship. In that relationship, the effect of new additional terms between the bending strain and displacement has been investigated in the warping problem. Natural co ordinate based strains, stresses and constitutive equations are used. The assumed natural strain method is used to alleviate both membrane and shear locking behavior from the element. The Lanczos method is employed in the calculation of the eigenvalues of laminated composite structures and the Gauss integration rule is adopted to evaluate the mass matrix. The numerical examples are compared with the analytical solutions to validate the current formulation and the results presented could be useful for the understanding of the behaviour of laminates under free vibration conditions.

• Journal of the Korean Society of Manufacturing Technology Engineers
• /
• v.7 no.6
• /
• pp.110-116
• /
• 1998

For the vibration analysis of 3-dimensional piping system containing fluid flow, a transfer matrix method is presented. The fluid velocity and pressure were considered, that coupled to longitudinal and flexural vibrations. Transfer matrices and point matrices were derived from direct solutions of the differential equations of motion of pipe conveying fluids, and the variations of natural frequency with flow velocity for 3-dimensional piping system were investigated.

• Proceedings of the IEEK Conference
• /
• 1998.06a
• /
• pp.359-362
• /
• 1998

2차원 혼합 모드 소자-회로 시뮬레이터를 이용한 과도상태 해석의 알고리즘을 제시한다. 1변수 muller 및 regular falsi법을 회로의 절점 전압과 분기(branch) 전류를 계산하는데 적용하였다. 제안된 알고리즘의 정확도와 유호성을 검증하기 위해 PN 다이오드의 양극(anode)에 저항이 직렬로 연결된 회로의 모의실험을 수행한 결과, MEDICI의 모의실험 결과에 비해 과도상태에서 전류 및 전압 특성은 각각 0.06%, 0.2% 오차 범위 한도 내에서 일치함을 보였다.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 1997.10a
• /
• pp.659-663
• /
• 1997

In this paper, a simplified vibration analysis model for bellows was presented to avoid excessive effort required for shell model. To reduce degree of freedom, bellows was modelled using one dimensional beam element. The equivalent mass and stiffness matrices were obtained based upon Guyan reduction process. The results were compared with the confirmed results, which were in good agreement.