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

Behavior of underground structural systems is usually complicated because of various unknown parameters. In order to construct those structural systems safely and economically, exact identification of the system parameters and accurate analysis of the system behaviors are essentially required. In this study, a forward analysis program, which is able to eliminate numerical errors due to far field boundary effect, is developed by coupling finite and boundary elements. In this coupled analysis, boundary elements are used in the semi-infinite domain where stress variation is small, and finite elements in the stress concentration region where material nonlinearity should be considered. Then, a back analysis program which can identify the system parameters is developed using the direct method to be combined with the forward analysis program. The elastic modulus and initial stress, which are most important in the description of the behavior of underground structures, are taken as the system parameters. A simple example is examined 0 show that the method can be used effectively.

• 대한전기학회논문지:전기기기및에너지변환시스템부문B
• /
• 제55권6호
• /
• pp.306-313
• /
• 2006

A three dimensional adaptive finite element refinement algorithm is developed for non-linear magnetostatic field problems. In the method, the edge elements are used for finite element formulation, and the local error in each element is estimated from the fact that the tangential components of magnetic field intensity and the normal components of magnetic flux density should be continuous at the interface of the two adjacent elements. Based on the estimated error, the elements which have big error are divided into several elements using bisection method. The effectiveness of the developed algorithm is proved through numerical examples.

• Journal of Advanced Marine Engineering and Technology
• /
• 제6권1호
• /
• pp.41-48
• /
• 1982

One of the important subjects in fracture mechanics study is to analyze the stress intensity factor. In this paper, the stress intensity factor in Mode I ($K^{I}$) is determined by J-integral using the finite element method. In this investigation, the values of $K^{I}$ are computed for distorted and undistorted elements of 8-noded isoparametric finite elements. The numerical results obtained are summarized as follows. (1) Through a relatively coarse mesh, the $K^{I}$ values obtained by this method are fairly good accuracy. (2) The $K^{I}$ values for the distorted elements appear to be better than those obtained using the undistorted mesh. (3) Within the limits of these analyses, the solutions obtained through the integral paths in the medium region of elements approach to the analytical solution most closely.

• 한국안전학회지
• /
• 제7권2호
• /
• pp.30-41
• /
• 1992

Analysis of Dynamic Characterisocs of Rectangular Plate by Finite Element Method. Dynamic characteristics of a rectangular plate with opening in it is studied by finite element method. To investigate these characteristics 12 degrees of freedom membrane finite element in used. The rectangular membrane finite elements are defined by specifying geometry, internal displacement functions and strain-displacement relations. Then, the governing equation for the finite element is derived by energy method. To derive the mass matrix and stiffness matrix of the element, expressions for strain and kineic energy in terms of the node displacement are generated. In constructing the overall structure matrix, the matrix of each elements are superposed and partitioned by applying the given boundary condition to obtain a nonslngular matrix. To find the natural freguencies and viration modes, the eigen values and the corresponding eigen vectors are computed by the computer using well known Jacobi power method. In order to verify the capability of the membrane finite element, a flat rectangular plate is analyzed first, and the result is compared with well known analytical results to show the good agreement. A rectangular plate with opening in It is analyzed with the same finite element. The results are presented in this paper. Unfortunately, the literature study could not provide with some results to compare, but the results reveal that the output of this research is phlslcally reasonable. And the results of this research are useful not only in practice but also for the future experimental research in comparison purpose.

• Journal of Mechanical Science and Technology
• /
• 제18권11호
• /
• pp.1989-1995
• /
• 2004

Finite element analysis(FEA) is the most popular numerical method to simulate plasticity-induced fatigue crack closure and can predict fatigue crack closure behavior. Finite element analysis under plane stress state using 4-node isoparametric elements is performed to investigate the detailed closure behavior of fatigue cracks and the numerical results are compared with experimental results. The mesh of constant size elements on the crack surface can not correctly predict the opening level for fatigue crack as shown in the previous works. The crack opening behavior for the size mesh with a linear change shows almost flat stress level after a crack tip has passed by the monotonic plastic zone. The prediction of crack opening level presents a good agreement with published experimental data regardless of stress ratios, which are using the mesh of the elements that are in proportion to the reversed plastic zone size considering the opening stress intensity factors. Numerical interpolation results of finite element analysis can precisely predict the crack opening level. This method shows a good agreement with the experimental data regardless of the stress ratios and kinds of materials.

• Structural Engineering and Mechanics
• /
• 제73권1호
• /
• pp.1-16
• /
• 2020

Finite elements based on the partition of unity (PU) approximation have powerful capabilities for p-adaptivity and solutions with high smoothness without remeshing of the domain. Recently, the PU approximation was successfully applied to the three-node shell finite element, properly eliminating transverse shear locking and showing excellent convergence properties and solution accuracy. However, the enrichment with the PU approximation results in a significant increase in the number of degrees of freedom; therefore, it requires greater computational cost, thus making it less suitable for practical engineering. To circumvent this disadvantage, we propose a new strategy to decrease the total number of degrees of freedom in the existing PU-based shell element, without loss of optimal convergence and accuracy. To alleviate the locking phenomenon, we use the method of mixed interpolation of tensorial components and perform convergence studies to show the accuracy and capability of the proposed shell element. The excellent performances of the new shell elements are illustrated in three benchmark problems.

• 한국전산구조공학회논문집
• /
• 제19권1호
• /
• pp.39-46
• /
• 2006

본 연구는 파 해석에 있어서 공간-시간 분할 개념을 도입하여 켈러킨 방법으로 해석하였다. 공간-시간 유한요소법은 오직 공간에 대해서만 분할하는 일반적인 유한요소법보다 간편하다. 비교적 큰 시간간격에 대해서 공간과 시간을 동시에 분할하는 방법을 제시하며 가중잔차법이 공간-시간 영역에서 유한요소 정식화에 이용되었다. 큰 시간 간격으로 인하여 문제의 해가 발산하는 경우가 동적인 문제에서 흔히 발생한다. 이러한 결점을 보완한 사각형 공간-시간 요소를 취하여 문제를 해석하고 해의 안정에 대해 기술하였다. 다수의 수치해석을 통하여 이 방법이 효과적 임을 알 수 있었다.

• Structural Engineering and Mechanics
• /
• 제53권2호
• /
• pp.205-226
• /
• 2015

Finite element method (FEM) is an effective quantitative method to solve complex engineering problems. The basic idea of FEM for a complex problem is to be able to find a solution by reducing the problem made simple. If mathematical tools are inadequate to obtain precise result, even approximate result, FEM is the only method that can be used for structural analyses. In FEM, the domain is divided into a large number of simple, small and interconnected sub-regions called finite elements. FEM has been used commonly for linear and nonlinear analyses of different types of structures to give us accurate results of plane stress and plane strain problems in civil engineering area. In this paper, FEM is used to investigate stress analysis of a shear wall which is subjected to concentrated loads and fundamental principles of stress analysis of the shear wall are presented by using matrix displacement method in this paper. This study is consisting of two parts. In the first part, the shear wall is discretized with constant strain triangular finite elements and stiffness matrix and load vector which is attained from external effects are calculated for each of finite elements using matrix displacement method. As to second part of the study, finite element analysis of the shear wall is made by ANSYS software program. Results obtained in the second part are presented with tables and graphics, also results of each part is compared with each other, so the performance of the matrix displacement method is demonstrated. The solutions obtained by using the proposed method show excellent agreements with the results of ANSYS. The results show that this method is effective and preferable for the stress analysis of shell structures. Further studies should be carried out to be able to prove the efficiency of the matrix displacement method on the solution of plane stress problems using different types of structures.

• 동력기계공학회지
• /
• 제11권1호
• /
• pp.92-97
• /
• 2007

It is important to compute the structural analysis of plate structures in structural design. In this paper, the author uses the finite element-transfer stiffness coefficient method (FE-TSCM) for the structural analysis of plate structures. The FE-TSCM is based on the concept of the successive transmission of the transfer stiffness coefficient method and the modeling technique of the finite element method (FEM). The algorithm for in-plane structural analysis of a rectangular plate structure is formulated by using the FE-TSCM. In order to confirm the validity of the FE-TSCM for structural analysis of plate structures, two numerical examples for the in-plane structural analysis of a plate with triangular elements and the bending structural analysis of a plate with rectangular elements are computed. The results of the FE-TSCM are compared with those of the FEM on a personal computer.

• Advances in Computational Design
• /
• 제5권3호
• /
• pp.277-289
• /
• 2020

This article describes the technology for constructing of a multiply-connected three-dimensional area's finite element representation. Representation of finite-element configuration of an area is described by a discrete set that consist of the number of nodes and elements of the finite-element grid, that are orderly set of nodes' coordinates and numbers of finite elements. Corresponding theorems are given, to prove the correctness of the solution method. The adequacy of multiply-connected area topology's finite element model is shown. The merging of subareas is based on the criterion of boundary nodes' coincidence by establishing a simple hierarchy of volumes, surfaces, lines and points. Renumbering nodes is carried out by the frontal method, where nodes located on the outer edges of the structure are used as the initial front.