• Proceedings of the KSME Conference
• /
• 2000.04b
• /
• pp.51-56
• /
• 2000

The Finite Element Solutions Is reported on solid-liquid phase change in porous media with natural convection including freezing. The model is based on volume averaged transport equations, while phase change is assumed to occur over a small temperature range. The FEM (Finite Element Method) algorithm used in this study is 3-step time-splitting method which requires much less execution time and computer storage the velocity-pressure integrated method and the penalty method. And the explicit Lax-Wendroff scheme is applied to nonlinear convective term in the energy equation. For natural convection including melting and solidification the numerical results show reasonable agreement with FDM (Finite Difference Method) results.

• Transactions of Materials Processing
• /
• v.11 no.5
• /
• pp.457-464
• /
• 2002

A finite element inverse analysis is utilized to consider forming effects of an S-rail on the assessment of the crashworthiness with small amount of computation time. A crash analysis can be directly performed after the inverse simulation of a forming process without a smoothing or remeshing scheme. The direct mesh mapping method is used to calculate an initial guess from a sliding constraint surface that is extracted from the die and punch set. Analysis results demonstrate that energy absorption of structures is increased when simulation considers forming effects of thickness variation and work hardening. The finite element inverse analysis is proved to be an effective tool in consideration of forming effects for the crash analysis.

• Journal of Welding and Joining
• /
• v.16 no.4
• /
• pp.109-116
• /
• 1998

The three-dimensional profile of a solder fillet is predicted by minimizing the surface tension and gravity energies of the solder joint using finite element modeling. Geometric complexity stemming from the inclined plane of the gullwing lead is resolved by employing three element types. These element types are used to describe the joint profile formed on the vertical, inclined and interfacial planes. The predicted solder joint profiles show good agreements with the experimental data provided that the solder volume is adjusted considering the wicking effects. Effects of the pad length, inclined lead angle and solder volume on joint profiles are also investigated.

• Nuclear Engineering and Technology
• /
• v.14 no.4
• /
• pp.178-185
• /
• 1982

Albedo-type boundary condition is incorporated into the finite element formulation of the cubic Hermite polynomials for the two-dimensional solution of the two-group diffusion problem. Two modifications are introduced with respect to the conventional expression for the weak form of the group diffusion equation with the zero flux or zero current boundary condition and the cubic element functions over the boundary nodes. The finite element formulations obtained from those modifications are tested with the two-dimensional ZION problem. The numerical effectiveness of the modifications are examined.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.12
• /
• pp.3007-3019
• /
• 1993

For the same configuration, the stiffness of 6-node two dimensional isoparametric is stiffer than that of 8-node two dimensional isoparametric element. This phenomenon may be called 'Relative Stiffness Stiffening Phenomenon.' In this paper, the relative stiffness stiffening phenomenon was studied, and could be corrected by modifying the position of Gauss integral points used in the numerical integration of the stiffness matrix. For the same deformation (bending) energy of 6-node and 8-node two dimensional isoparametric elements, Gauss integral points of 6-node element have to move closer, in comparison with those of 8-node element, in the case of numerical integration along the thickness direction.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2000.10a
• /
• pp.217-221
• /
• 2000

Nonlinear stress analysis of nuclear containment building wall element is carried out using microscopic material model. The present study mainly focuses on the finite element analysis of the nuclear containment building wall element under biaxial tensile stresses and it evaluates the perfomance of adopted microscopic material model in the membrane energy dominant situation. From the numerical analysis, the adopted material model peforms well and has a good agreement with experiment result. Finally, the result of present study can be severed as a benchmark test when concrete material model is in need of evaluation.

• Structural Engineering and Mechanics
• /
• v.5 no.6
• /
• pp.785-801
• /
• 1997

Reinforced Concrete beams with tension and compression softening material constitutive laws are studied. Energy-based and non-local regularisation techniques are presented and applied to a R.C. element. The element characteristics (sectional tangent stiffness matrix, element tangent stiffness matrix restoring forces) are directly derived from their symbolic expressions through numerical integration. In this way the same spatial grid allows us to obtain a non-local strain estimate and also to sample the contributions to the element stiffness matrix. Three examples show the spurious behaviors due to the strain localization and the stabilization effects given by the regularisation techniques, both in the case of tension and compression softening. The possibility to overestimate the ultimate load level when the non-local strain measure is applied to a non softening material is shown.

• KEPCO Journal on Electric Power and Energy
• /
• v.1 no.1
• /
• pp.127-134
• /
• 2015

Recently, the external bonding of carbon fiber reinforced polymer (CFRP) sheets has come to be regarded as a very effective method for strengthening of reinforced concrete structures. The behavior of CFRP-strengthened RC structure is mainly governed by the interfacial behavior, which represents the stress transfer and relative slip between concrete and the CFRP sheet. In this study, the effects of bonded length, width and concrete strength on the interfacial behavior are verified and a bond-slip model is proposed. The proposed bond-slip model has nonlinear ascending regions and exponential descending regions, facilitated by modifying the conventional bilinear bond-slip model. Finite element analysis results of interface element implemented with bond-slip model have shown good agreement with the experimental results performed in this study. It is found that the failure load and strain distribution predicted by finite element analysis with the proposed bond-slip are in good agreement with results of experiments.

• Structural Engineering and Mechanics
• /
• v.87 no.2
• /
• pp.191-200
• /
• 2023

This paper presents a mechanical model of a "tapered composite shaft" rotating at a constant speed around its axis. The spatial equations of motion are solved using the Lagrange technique, and a finite element approach is employed to construct the model. Theoretical analysis is used to compute the kinetic and strain energies. A comparison is made between conventional finite element methods and hierarchical finite element methods, indicating that the former uses fewer elements and provides higher accuracy in determining natural frequencies. Numerical calculations are performed to determine the eigen frequencies and critical speeds of the rotating composite shaft. The critical speeds of composite shaft systems are compared with existing literature to validate the proposed model.

• Journal of Radiation Protection and Research
• /
• v.23 no.3
• /
• pp.123-138
• /
• 1998

Characteristics of element responses of Panasonic UD802 personnel dosimeters in the X, ${\beta}$, ${\gamma}$, ${\gamma}/X$, ${\gamma}/{\beta}$ and ${\gamma}$/neutron mixed fields were assessed. A dose-response algorithm has been developed to decide the high probability of a radiation type and energy by using the distribution in all six ratios of the multi-element TLD. To calculate the 4-element response factors and ratios between the elements of the Panasonic TLDs in the X, $\beta$, and $\gamma$ radiation fields, Panasonic’s UD802 TLDs were irradiated with KINS’s reference irradiation facility. In the photon radiation field, this study confirms that element-3 (E3) and element-4 (E4) of the Panasonic TLDs show energy dependent both in low- and intermediate-energy range, while element-1 (E1) and element-2 (E2) show little energy dependency in the entire whole range. The algorithm, which was developed in this study, was applied to the Panasonic personnel dosimetry system with UD716AGL reader and UD802 TLDs. Performance tests of the algorithm developed was conducted according to the standards and criteria recommended in the ANSI N13.11. The sum of biases and standard deviations was less than 0.232. The values of biases and standard deviations are distributed within a triangle of a lateral value of 0.3 in the ordinate and abscissa, With the above algorithm, Panasonic TLDs satisfactorily perform optimum dose assessment even under an abnormal response of the TLD elements to the energy imparted. This algorithm can be applied to a more rigorous dose assessment by distinguishing an unexpected dose from the planned dose for the most practical purposes, and is useful in conducting an effective personnel dose control program.