• 한국압력기기공학회 논문집
• /
• 제13권1호
• /
• pp.84-91
• /
• 2017

This paper is proposed to select the optimal finite element type in finite element analysis. Based on the NUREG reports, static analyses were performed using a commercial analysis program, $ABAQUS^{TM}$. In this study, we used a nonlinear kinematic hardening model proposed by Chaboche. The analysis result of solid elements by inputting the same material constants was different from the results of the NUREG report. This is resulted from the difference between shell element and solid element. Therefore, the material constants that have similar result to the experimental result were determined and compared according to element type. In case of using solid element for efficient finite element analysis, we confirmed that the use of C3D8I element type(incompatible mode 8-node linear brick element) leads the accurate result while reducing the analysis time.

• 한국콘크리트학회:학술대회논문집
• /
• 한국콘크리트학회 1996년도 봄 학술발표회 논문집
• /
• pp.247-252
• /
• 1996

In the concrete structures, cracks occur in various causes and the cracks seriously affect the functions of structures. The analysis techniques of progressive crack in the concrete have been improved with the advance of numerical techniques. The discrete crack model used in finite element program for the analysis of progressive failure is very effective, but it can not be easily implemented into numerical procedures because of difficult handing of nodal points in finite element meshes for crack growth. This paper introduces one of the techniques which skips the difficulty. In this paper, the modeling of progressive failure using finite element formulation is explained for the analysis of concrete fracture. The discontinuous element using the discontinuous shape function and the dual mapping technique in the numerical integration are implemented into finite element code for this purpose. It is shown that developed finite element program can predict the quasi-brittle behavior of concrete including ultimate load. The comparisons of the analysis results with other data are also shown.

• 한국지진공학회논문집
• /
• 제23권1호
• /
• pp.83-88
• /
• 2019

Structural analysis remains as an essential part of any integrated civil engineering system in today's rapidly changing computing environment. Even with enormous advancements in capabilities of computers and mobile tools, enhancing computational efficiency of algorithms is necessary to meet the changing demands for quick real time response systems. The finite element method is still the most widely used method of computational structural analysis; a robust, reliable and automated finite element structural analysis module is essential in a modern integrated structural engineering system. To be a part of an automated finite element structural analysis, an efficient adaptive mesh generation scheme based on R-H refinement for the mesh and error estimates from representative strain values at Gauss points is described. A coefficient that depends on the shape of element is used to correct overly distorted elements. Two simple case studies show the validity and computational efficiency. The scheme is appropriate for nonlinear and dynamic problems in earthquake engineering which generally require a huge number of iterative computations.

• 한국농공학회지
• /
• 제22권2호
• /
• pp.83-100
• /
• 1980

This study was intended to computerize the structural analysis of two-dimensional continuum by finite element method, and to provide a preparatory basis for more sophisticated and more generalized computer programs of this kind. A computer program, applicable to any shape of two-dimensional continuum, was formulated on the basis of 16-degree-of- freedom rectangular element. Various computational aspects pertaining to the implementation of finite element method were reviewed and settled in the course of programming. The validity of the program was checked through several case studies. To assess the accuracy and the convergence characteristics of the method, the results computed by the program were compared with solutions by other methods, namely the analytical Navier's method and the framework method. Through actual programming and analysis of the computed results, the following facts were recognized; 1) The stiffness matrix should necessarily be assembled in a condensed form in order to make it possible to discretize the continuum into practically adequate number of elements without using back-up storage. 2) For minimization of solution time, in-core solution of the equilibrium equation is essential. LDLT decomposition is recommended for stiffness matrices condensed by the compacted column storage scheme. 3) As for rectangular plates, the finite element method shows better performances both in the accuracy and in the rate of convergence than the framework method. As the number of elements increases, the error of the finite element method approaches around 1%. 4) Regardless of the structural shape, there is a uniform tendency in convergence characteristics dependent on the shape of element. Square elements show the best performance. 5) The accuracy of computation is independent of the interpolation function selected.

• Geomechanics and Engineering
• /
• 제12권6호
• /
• pp.1021-1046
• /
• 2017

The assessment of slope stability is an essential task in geotechnical engineering. In this paper, a three-dimensional (3D) finite element analysis (FEA) was employed to investigate the performance of different shear pin arrangements to increase the stability of a soil block resting on an inclined plane with a low-interface friction plane. In the numerical models, the soil block was modeled by volume elements with linear elastic perfectly plastic material in a drained condition, while the shear pins were modeled by volume elements with linear elastic material. Interface elements were used along the bedding plane (bedding interface element) and around the shear pins (shear pin interface element) to simulate the soil-structure interaction. Bedding interface elements were used to capture the shear sliding of the soil on the low-interface friction plane while shear pin interface elements were used to model the shear bonding of the soil around the pins. A failure analysis was performed by means of the gravity loading method. The results of the 3D FEA with the numerical models were compared to those with the physical models for all cases. The effects of the number of shear pins, the shear pin locations, the different shear pin arrangements, the thickness and the width of the soil block and the associated failure mechanisms were discussed.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2003년도 가을 학술발표회 논문집
• /
• pp.3-17
• /
• 2003

This paper deals with the problems and their possible solutions in the development of finite element for analysis of shell. Based on these solution schemes, a series of flat shell elements are established which show no signs of membrane locking and other defects even though the coarse meshes are used. In the element formulation, non-conforming displacement modes are extensively used for improvement of element behaviors. A number of numerical tests are performed to prove the validity of the solutions to the problems involved in establishing a series of high performance flat shell elements. The test results reveal among others that the high accuracy and fast convergence characteristics of the elements are obtainable by the use of various non-conforming modes and that the ‘Direct Modification Method’ is a very useful tool for non-conforming elements to pass the patch tests. Furthermore, hierarchical and higher order non-conforming modes are proved to be very efficient not only to make an element insensitive to the mesh distortion but also to remove the membrane locking. Some numerical examples are solved to demonstrate the validity and applicability of the presented elements to practical engineering shell problems.

• Coupled systems mechanics
• /
• 제3권3호
• /
• pp.305-318
• /
• 2014

The study deals with the physical modeling of a typical building frame resting on pile foundation and embedded in cohesive soil mass using complete three-dimensional finite element analysis. Two different pile groups comprising four piles ($2{\times}2$) and nine piles ($3{\times}3$) are considered. Further, three different pile diameters along with the various pile spacings are considered. The elements of the superstructure frame and those of the pile foundation are descretized using twenty-node isoparametric continuum elements. The interface between the pile and pile and soil is idealized using sixteen-node isoparametric surface elements. The current study is an improved version of finite element modeling for the soil elements compared to the one reported in the literature (Chore and Ingle 2008). The soil elements are discretized using eight-, nine- and twelve-node continuum elements. Both the elements of superstructure and substructure (i.e., foundation) including soil are assumed to remain in the elastic state at all the time. The interaction analysis is carried out using sub-structure approach in the parametric study. The total stress analysis is carried out considering the immediate behaviour of the soil. The effect of various parameters of the pile foundation such as spacing in a group and number piles in a group, along with pile diameter, is evaluated on the response of superstructure. The response includes the displacement at the top of the frame and bending moment in columns. The soil-structure interaction effect is found to increase displacement in the range of 58 -152% and increase the absolute maximum positive and negative moments in the column in the range of 14-15% and 26-28%, respectively. The effect of the soil- structure interaction is observed to be significant for the configuration of the pile groups and the soil considered in the present study.

• Structural Engineering and Mechanics
• /
• 제28권5호
• /
• pp.603-633
• /
• 2008

A four-node plate finite element for the analysis of laminated composites which is developed using strain gradient notation is presented. The element is based on a first-order shear deformation theory and on the equivalent lamina assumption. Strains and stresses can be calculated at different points through the thickness of the plate. They are averaged values due to the equivalent lamina assumption. A shear correction factor is used as the transverse shear strain is taken to be constant over the plate thickness while its actual variation is parabolic. Strain gradient notation, which is physically interpretable, allows for the detailed a-priori analysis of the finite element model. The polynomial expansions are inspected and spurious terms responsible for modeling errors are identified in the shear strains polynomial expansions. The element is corrected by simply removing the spurious terms from the shear strains expansions. The element is implemented into a FORTRAN finite element code in two versions; namely, with and without spurious terms. Results are compared to show the effects of the spurious terms on the solutions. It is also shown that a refined mesh composed of corrected elements provides solutions which approximate very well the analytical solutions, validating the procedure.

• 한국해안해양공학회지
• /
• 제14권1호
• /
• pp.86-93
• /
• 2002

항만부진등 해석에 사용되는 유한요소의 적정 크기를 검토하기 위해 완전개방 직사각형 항만에 대해서 다양한 크기의 요소를 사용한 수치실험을 실시하였다. 수치실험 결과, 공진주기와 증폭비의 허용 오차율을 모두 2%로 하는 경우 파장당 9개, 1%로 하는 경우 12개의 요소가 필요한 것으로 나타났다. 또한, 요소 수의 증가에 따라 공진주기의 오차율은 선형적으로, 증폭비의 오차는 진동하면서 감소하였으며, 요소 수를 9개보다 적게 하는 경우 증폭비의 오차가 상대적으로 크게 증가하였다

• 대한전기학회:학술대회논문집
• /
• 대한전기학회 1989년도 추계학술대회 논문집 학회본부
• /
• pp.41-44
• /
• 1989

This paper describes an analysis of the three-dimensional eddy current problems by the finite element method using magnetic vector potential and electric scalar potential. The finite element formulation uses first-order shape functions and tetrahedral elements. The validity of this formalation is ensured as the result of the sphere conductor model problem in a sinusoidal magnetic field.