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

Adaptive finite element analysis, in which its solution error meets with the user defined allowable error, is recently used to improve the reliability of finite element analysis results. This adaptive analysis is composed of two procedures; one is the error estimation of an analysis result and the other is the reconstruction of finite elements. In the (p-method, an element size is controlled by relocating of nodal positions (r-method) and the order of an element shape function is determined by the hierarchical polynomial (p-method) corresponding to the clement solution error by the enhanced SPR. In order to show the effectiveness and the accuracy of the suggested rp-method, various numerical examples were analyzed and these analysis results were examined by comparing with those obtained by the existed methods.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1995.10a
• /
• pp.221-226
• /
• 1995

터빈 익렬, 펌프 익차, 원형 냉각탑, 치차 등과 같이 동일한 형상이 원주 방향으로 반복되어 있는 순환 대칭 구조물의 진동특성을 유한 요소법을 사용하여 해석하는 경우에 전체구조를 모델링하는 대신에 구조물을 동일한 형상의 부분구조로 분할하여 부분구조 한개만을 모델링하고 분할된 경계에서 적절한 경계조건을 부과하여 진동해석을 수행함으로서 컴퓨터 기억용량을 절감시키고 계산시간을 단축할 수 있는 방법이 널리 사용되고 있다. Orris and Petyt[1]는 부분구조의 양쪽 분할 경계면, 즉 연결 경계상에 있는 절점변위의 상관관계를 복소파동전파식을 이용해서 구하여 부분구조의 감소된 복소강성행렬 및 질량행렬을 만들고 실수부와 허수부를 분리하여 유한요소해석을 수행하는 방법을 제안하였다. 유한요소 프로그램 ANSYS[2]에서는 이와 같은 방법을 사용하고 있다. Thomas[3]는 순회 정규모드를 이용하였고, 참고문헌[4]에서는 순회행렬을 이용하였다. 또한 유한요소 프로그램 MSC/NASTRAN[5]에서는 푸리에 급수를 이용하고 유한요소 절점의 위치 및 변위를 원통 좌표계를 표현하여 순환대칭구조물의 유한요소해석을 수행할 수 있도록 되어있다. 본 논문에서는 순환 대칭구조물의 형상의 주기성과 순환성을 고려하여 이산퓨리에 변환을 이용함으로써 순환대칭구조물의 유한요소진동해석을 체계적으로 저용량의 컴퓨터에서 신속하고 정확하게 수행할 수 있는 방법을 제안하고자 한다.

• 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.

• Computational Structural Engineering
• /
• v.3 no.2
• /
• pp.77-88
• /
• 1990

Two perforated plates (a square plate and a rectangular plate having an aspect ratio 1.57(L/sub x/=11, L/sub y/=7)) are taken as analysis examples. Each of these plates is given some changes in the boundary conditions. The size of cutouts as well as their locations are also changed in order to examine the variation of two eigenvalues corresponding to the fundamental mode. The relationship between two eigenvalues is established by changing the magnitude of edge thrust.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.21 no.12
• /
• pp.1973-1983
• /
• 1997

Higher order isoparametric elements are usually used in the finite element analysis because they can represent easily the geometric shape of a complex structure ad can improve the solution quality. When these elements are used, the position of internal nodes affects greatly on the solution accuracy. Decreasing of the accuracy related to the position of internal nodes is due to non-conformal mapping often results in an unstable Jacobian value. This paper, in order to remove this difficulty, suggests a modified shape function which can establish conformal mapping between two coordinate systems. Numerical experiments with the proposed shape function show that a stable solution can be obtained without respect to the position of internal nodes, and offer convenience that one can take arbitrarily the position of internal nodes considering only the geometric shape of element boundaries.

• 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.

• Journal of the Society of Naval Architects of Korea
• /
• v.53 no.5
• /
• pp.388-399
• /
• 2016

A motion analysis of an underwater towed cable is a complex task due to its nonlinear nature of the problem. The major source of the nonlinearity of the underwater cable analysis is that the motion of the cable involves large rigid-body motion. This large rigid-body motion makes difficult to use standard displacement-based finite element method. In this paper, the authors apply recently developed nodal position-based finite element method which can deal with the geometric nonlinearity due to the large rigid-body motion. In order to enhance the stability of the large-scale nonlinear cable motion simulation, an efficient time-integration scheme is proposed, namely predictor/multi-corrector Newmark scheme. Three different predictors are introduced, and the best predictor in terms of stability and robustness for impulsive cable motion analysis is proposed. As a result, the nonlinear motion of underwater cable is predicted in a very efficient manner compared to the classical finite element of finite difference methods. The efficacy of the method is demonstrated with several test cases, involving static and dynamic motion of a single cable element, and also under water towed cable composed of multiple cable elements.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.31 no.6
• /
• pp.381-388
• /
• 2018

This paper presents a generalized finite element formulation for plastic hinge modeling based on lumped plasticity in the classical Euler-Bernoulli beam elements. In this approach, the plastic hinges are effectively modeled using proper enrichment functions describing weak discontinuities of the solution. The proposed methodology enables the insertion of plastic hinges at an arbitrary location without modifying the connectivity of elements. The formations of plastic hinges are instead achieved by hierarchically adding degrees of freedom to existing elements. Convergence analyses such as h- and p-extensions are performed to investigate the effectiveness of the proposed method. The analysis results indicate that the proposed generalized finite element method can achieve theoretical convergence rates for both cases where plastic hinges are located at nodes and within an element, thus demonstrating its accuracy.

• 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.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.33 no.4
• /
• pp.310-320
• /
• 2009

In this paper, we estimate area-averaged solution of peening residual stress using a 3-D multi-impact symmetry-cell FE model. The symmetry-cell model includes factors reflecting peening phenomena and plastic shot. Area-averaged solution is much closer to XRD experimental solution than 4-node-averaged solution in plastic shot FE model. We then obtain FE Almen saturation curve corresponding to experimental Almen curve based on area-averaged solution. Using the curve, we obtain FE area-averaged solution in major peening materials, and compare the FE solution with experimental solution. In peening materials, surface, maximum compressive residual stress and deformation depth reach experimental solutions. Thus, FE Almen curve is useful for estimation of residual stress solution and could improve the efficiency of peening process. Consequently, it is confirmed that concept of area-averaged solution is the realistic analytical method for evaluation of peening residual stress.