• Proceedings of the Korea Society for Simulation Conference
• /
• 2001.10a
• /
• pp.373-379
• /
• 2001

There is a big issue to generate 3D hexahedral finite element (FE) model, since a process to divide the whole domain into several simple-shaped sub-domains is required before generating a continuous mesh with mapped mesh generators. In general, it is nearly impossible to set up proper division numbers interactively to keep mesh connectivity between sub-domains on a complicated arbitrary-shaped domain. If mesh continuity between sub-domains is not required in an analysis, this complicated process can be omitted. Element-free Galerkin method (EFGM) can accept discontinuous meshes, which only requires nodal information. However it is difficult to choose a reasonable influenced domain in moving least squares scheme with non-uniformly distributed nodes in discontinuous FE models. A new FE scheme fur discontinuous mesh is proposed in this paper by applying improved EFGM with some modification to derive FE approximated function in discontinuous parts. Its validity is evaluated on linear elastic problems.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.25 no.11
• /
• pp.1697-1704
• /
• 2001

A finite element model for polarization switching in electro-mechanically coupled materials is proposed and applied to predict the switching behavior of a two-dimensional ferroelectric ceramic. A complicated micro-structure existing in the material is modeled as il continuum body and a simple 3 node triangle finite element with nodal displacement and voltage degrees of freedom is used for a finite element analysis. The elements use nonlinear constitutive equations, switching criterion and kinetic relation, fur representation of material response at strong electric and stress fields. The polarization state of the material is represented by internal variables in each element, which are updated at each simulation step based on the proposed constitutive equations. The model reproduces strain and electric displacement hysteresis loops observed in the material.

• Structural Engineering and Mechanics
• /
• v.63 no.2
• /
• pp.251-257
• /
• 2017

A local damage identification method using measured structural static displacement is proposed in this study. Based on the residual force vector deduced from the static equilibrium equation, residual strain energy (RSE) is introduced, which can localize the damage in the element level. In the case of all the nodal displacements are used, the RSE can localize the true location of damage, while incomplete displacement measurements are used, some suspicious damaged elements can be found. A model updating method based on static displacement response sensitivity analysis is further utilized for accurate identification of damage location and extent. The proposed method is verified by two numerical examples. The results indicate that the proposed method is efficient for damage identification. The advantage of the proposed method is that only limited static displacement measurements are needed in the identification, thus it is easy for engineering application.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.24 no.7 s.178
• /
• pp.1847-1853
• /
• 2000

A computer simulation technique for 2-dimensional laser generated ultrasonic waves was developed for visualization and investigation of ultrasonic propagation in solids. The technique is similar to a finite difference method (FDM) and a mass-particle model method, but uses a new nodal calculation method based on fundamental consideration of an elastic wave equation. By this method, the propagation behavior oflaser generated ultrasonic wave in thermoelastic and ablation mode is visualized and shows good agreement with previous experimental result or the numerical analysis result by Green function.

• Nuclear Engineering and Technology
• /
• v.52 no.11
• /
• pp.2415-2421
• /
• 2020

The axial heterogeneous void distribution in a fuel channel is a relevant and important issue during nuclear reactor analysis for LWR, especially for boiling water channel-type reactors. Variation of the coolant density in fuel channel has an effect on the neutron spectrum that will in turn have an impact on the values of absolute reactivity, the void reactivity coefficient, and the fuel isotopic compositions during irradiation. This effect is referring to as the history effect in light water reactor calculations. As the void reactivity effect is positive in RBMK type reactors, the underestimation of water density heterogeneity in 3D reactor core numerical calculations could cause an uncertainty during assessment of safe operation of nuclear reactor. Thus, this issue is analysed with different cross-section libraries which were generated with WIMS8 code at different reference water densities. The libraries were applied in single fuel model of the nodal code of QUABOX-CUBBOX/HYCA. The thermohydraulic part of HYCA allowed to simulate axial water distribution along fuel assembly model and to estimate water density history effect for RBMK type fuel.

• The Transactions of The Korean Institute of Electrical Engineers
• /
• v.58 no.2
• /
• pp.264-269
• /
• 2009

Marginal Loss Factor(MLF) is represented as the sensitivity of transmission loss, which is computed from the change of the generation at slack bus by the change of the load at the arbitrary bus. The MLF dependent on the selection of slack bus is one of the key factors affecting nodal pricing, Genco's profits, social welfare(SW) and Nash Equilibrium in a competitive electricity market. This paper addresses the methodology of slack bus selection by using Cournot model of Cost Based Pool market. Numerical results from sample cases show that the slack bus of MLF of the highest average is beneficial from the view points of SW.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2002.04a
• /
• pp.3-11
• /
• 2002

A moment-curvature relationship to simulate the behavior of reinforced concrete (RC) columns under cyclic loading is introduced. Unlike previous moment-curvature models and the layered section approach, the unposed model takes into account the bond-slip effect by using a monotonic moment-curvature relationship constructed on the basis of the bond-slip relation and corresponding equilibrium equation at each nodal point. In addition, the use of curved unloading and reloading branches inferred from the stress-strain relation of steel gives more exact numerical result. The pinching enact caused by axial force is considered with an assumption that the absorbing energy corresponding to any deformation level maintains constant regardless of the magnitude of applied axial force. The advantages of the proposed model, comparing tn layered section approach, may be on the reduction in calculation time and memory space in case of its application to large structures.. Finally, correlation studies between analytical results and experimental studies are conducted to establish the validity of the proposed mood.

• Proceedings of the KSR Conference
• /
• 2008.11b
• /
• pp.1379-1388
• /
• 2008

In this paper, we present the equations of motion by which the natural vibration of a rotating annular disk can be accurately analyzed. These equations are derived from the theory of finite deformation and the principle of virtual work. The radial displacements of annular disk which is rotating at constant angular velocity are determined by non-linear equations formulated using 1-dimensional finite elements in radial direction. The equations of the in-plane vibrations at disturbed state are also formulated using 1-dimensional finite elements in radial direction along the number of nodal diameters. They are expressed as in functions of the radial displacements at the steady state and the disturbed displacements about the steady state. In-plane static deformation modes of the annular disk are used as the interpolation functions of 1-dimensional finite elements in radial direction. The natural vibrations of an annular disk with different boundary conditions are analyzed by using the presented model and the 3-dimensional finite element model to verify accuracy of the presented equations of motion. Its results are compared and discussed.

• Proceedings of the Korean Nuclear Society Conference
• /
• 1996.05a
• /
• pp.40-45
• /
• 1996

A new concept of equivalent homogenization is proposed. The concept employs new set of homogenized parameters: homogenized cross sections (XS) and interface matrix (IM), which relates partial currents at the celt interfaces. The idea of interface matrix generalizes the idea of discontinuity factors (DFs), proposed and developed by K.Koebke and K.Smith. The method of K.Smith can be simulated within framework of new method, while the new method approximates heterogeneous cell better in case of the steep flux gradients at the cell interfaces. The attractive shapes of new concept are: improved accuracy, simplicity of incorporation in the existing codes, equal numerical expenses in comparison to the K.Smith's approach. The new concept is useful for: (a) explicit reflector/baffle simulation; (b) control blades simulation; (c) mixed UO2/MOX core simulation, The offered model has been incorporated in the finite difference code and in the nodal code PANBOX. The numerical results show good accuracy of core calculations and insensitivity of homogenized parameters with respect to in- core conditions.

• Computers and Concrete
• /
• v.28 no.6
• /
• pp.621-634
• /
• 2021

For the reliable strut-and-tie model (STM) design of disturbed regions of concrete members, structural designers must accurately determine the strength of concrete struts to check the strength conditions of a selected STM el and the anchorage of reinforcing bars in nodal zones. In this study, the author proposed a consistent numerical method for strut strength, applicable to all two-dimensional STMs. The proposed method includes the effects of a biaxial stress state associated with tensile strains in reinforcing bars crossing a strut, deviation angle between strut orientation and compressive principal stress flow, and degree of confinement provided by reinforcement. The author examined the method's validity through the STM prediction of the ultimate strengths of 517 reinforced concrete (RC) deep beams, 24 RC panels, and 258 RC corbels, all tested to failure.