• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.5 no.3
• /
• pp.39-47
• /
• 1985

In this paper, a FEM analysis model was developed to solve the consolidation phenomena of embankment on the soft foundation. The developed FEM model was based on the Biot's consolidation equation which was coupled with one of three stress-strain constitutive relationships. In order to check the validity of the newly developed FEM model, the program input data were used by a test embankment which had been already constructed at Cubzac-les-ponts in France by Magnan et al. The FEM results compared to the experimental and analytical results which were obtained by the Magnan's group at Cubzac-les-ponts. The results compared showed that the consolidation phenomena were well explained by the author's FEM model which results were more accurate than the others. As for the pore water pressure, Christian-Boehmer's method used in this paper was considered preferable to Sandhu-Wilson's used by Magnan.

• Nuclear Engineering and Technology
• /
• v.20 no.2
• /
• pp.88-104
• /
• 1988

The FUel Rod Analysis(FURA) code is developed using two-dimensional finite element methods for axisymmetric and plane stress analysis of fuel rod. It predicts the thermal and mechanical behavior of fuel rod during normal and load follow operations. To evaluate the exact temperature distribution and the inner gas pressure, the radial deformation of pellet and clad, the fission gas release are considered over the full-length of fuel rod. The thermal element equation is derived using Galerkin's techniques. The displacement element equation is derived using the principle of virtual works. The mechanical analysis can accommodate various components of strain: elastic, plastic, creep and thermal strain as well as strain due to swelling, relocation and densification. The 4-node quadratic isoparametric elements are adopted, and the geometric model is confined to a half-pellet-height region with the assumption that pellet-pellet interaction is symmetrical. The pellet cracking and crack healing, pellet-cladding interaction are modelled. The Newton-Raphson iteration with an implicit algorithm is applied to perform the analysis of non-linear material behavior accurately and stably. The pellet and cladding model has been compared with both analytical solutions and experimental results. The observed and predicted results are in good agreement. The general behavior of fuel rod is calculated by axisymmetric system and the cladding behavior against radial crack is used by plane stress system. The sensitivity of strain aging of PWR fuel cladding tube due to load following is evaluated in terms of linear power, load cycle frequency and amplitude.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.39 no.8
• /
• pp.807-812
• /
• 2015

Li-ion batteries are energy sources that are widely used in applications such as notebooks, cellular phones, power tools, and vehicles. They are devices in which stored chemical energy is changed to electrical energy by electrochemical reactions. They have a high energy density, small size, and are lightweight. In particular, power tools and vehicles require high charge/discharge rates. Therefore, in this paper, we perform electrochemical simulations using a commercial finite-element analysis program to determine the high discharge-rate characteristics of Li-ion cells. In addition, by performing high discharge-rate simulations, we found that the limited discharge current was 63 A. Based on the results obtained, we investigate the behavior of Li-ion cells with a high rate of discharge.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.22 no.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.

• Composites Research
• /
• v.15 no.4
• /
• pp.23-31
• /
• 2002

The two dimensional size effect of specimen gauge section ($length{\;}{\times}{\;}width$) was investigated on the compressive behavior of a T300/924 $\textrm{[}45/-45/0/90\textrm{]}_{3s}$, carbon fiber-epoxy laminate. A modified ICSTM compression test fixture was used together with an anti-buckling device to test 3mm thick specimens with a $30mm{\;}{\times}{\;}30mm,{\;}50mm{\;}{\times}{\;}50mm,{\;}70mm{\;}{\times}{\;}70mm{\;}and{\;}90mm{\;}{\times}{\;}90mm$ gauge length by width section. In all cases failure was sudden and occurred mainly within the gauge length. Post failure examination suggests that $0^{\circ}$ fiber microbuckling is the critical damage mechanism that causes final failure. This is the matrix dominated failure mode and its triggering depends very much on initial fiber waviness. It is suggested that manufacturing process and quality may play a significant role in determining the compressive strength. When the anti-buckling device was used on specimens, it was showed that the compressive strength with the device was slightly greater than that without the device due to surface friction between the specimen and the device by pretoque in bolts of the device. In the analysis result on influence of the anti-buckling device using the finite element method, it was found that the compressive strength with the anti-buckling device by loaded bolts was about 7% higher than actual compressive strength. Additionally, compressive tests on specimen with an open hole were performed. The local stress concentration arising from the hole dominates the strength of the laminate rather than the stresses in the bulk of the material. It is observed that the remote failure stress decreases with increasing hole size and specimen width but is generally well above the value one might predict from the elastic stress concentration factor. This suggests that the material is not ideally brittle and some stress relief occurs around the hole. X-ray radiography reveals that damage in the form of fiber microbuckling and delamination initiates at the edge of the hole at approximately 80% of the failure load and extends stably under increasing load before becoming unstable at a critical length of 2-3mm (depends on specimen geometry). This damage growth and failure are analysed by a linear cohesive zone model. Using the independently measured laminate parameters of unnotched compressive strength and in-plane fracture toughness the model predicts successfully the notched strength as a function of hole size and width.