• 한국전산유체공학회:학술대회논문집
• /
• 2008.10a
• /
• pp.71-78
• /
• 2008

For analyses of multi-phase flows in a water-cooled nuclear power plant, a three-dimensional SIMPLE-algorithm based hydrodynamic solver CUPID-S has been developed. As governing equations, it adopts a two-fluid three-field model for the two-phase flows. The three fields represent a continuous liquid, a dispersed droplets, and a vapour field. The governing equations are discretized by a finite volume method on an unstructured grid to handle the geometrical complexity of the nuclear reactors. The phasic momentum equations are coupled and solved with a sparse block Gauss-Seidel matrix solver to increase a numerical stability. The pressure correction equation derived by summing the phasic volume fraction equations is applied on the unstructured mesh in the context of a cell-centered co-located scheme. This paper presents the numerical method and the preliminary results of the calculations.

• Structural Engineering and Mechanics
• /
• v.48 no.2
• /
• pp.275-289
• /
• 2013

The cross sections of multi-span beams are sometimes suddenly increased at the interior support of continuous beams to resist high negative moment. An earlier study on elastic lateral torsional buckling of stepped beams was conducted to propose new design equations. This research aims to continue the earlier study by considering the effect of inelastic buckling of stepped beams subjected to pure bending and general loading condition. A three-dimensional finite element-program ABAQUS and a statistical program MINITAB were used in the development of new design equations. The inelastic lateral torsional buckling strengths of 36 and 27 models for singly and doubly stepped beams, respectively, were investigated. The general loading condition consists of 15 loading cases based on the number of inflection point within the unbraced length of the stepped beams. The combined effects of residual stresses and geometrical imperfection were also considered to evaluate the inelastic buckling strengths. The proposed equations in this study will definitely improve current design methods for the inelastic lateral-torsional buckling of stepped beams and will increase efficiency in building and bridge design.

• Journal of the Korean Society of Safety
• /
• v.15 no.3
• /
• pp.14-22
• /
• 2000

A numerical analysis has been performed on the two-dimensional rectangular gallium melting problem using the enthalpy method. The major advantage of this method is that the physical domain is discretized with fixed grids without transforming variables and the interface conditions of phase change are accounted for the definition of suitable source terms in the governing equations. But in the fixed method, there is some ambiguity in defining the porosity constant which has no physical interpretation. If the velocity correction is included in the momentum equation, for the appropriate range of porosity constant, the realistic predictions are obtained. The object of the present work is to predict the phase change within the molten steel with thin riser slab using the modified enthalpy-porosity method. The computational procedures for predicting velocity and temperature are based on the finite volume method and the non-staggered grid system. The influence of natural convection on the melting process is considered. A comparison with the experimental results shows that the modified method is better than the previous one.

• Tunnel and Underground Space
• /
• v.4 no.2
• /
• pp.102-118
• /
• 1994

Thermomechanical analysis is conducted on the radioactive repository in deep rock mass considering the in-situ stress, excavation and thermal loading of a radioactive waste. Thermomechanical properties of a discontinuous rock mass are estimated by a theoretical method so called sequential analysis. Using the estimated properties as input for finite element analysis, the influence on temperature distribution and thermal stress is analyzed within the scope of 2-dimensional steady state and transient heat transfer and coupled thermal elastic plastic behaviour. Granitic rock mass is taken for this analysis. The analysis is done for two different rock mass conditions, i.e. continuous-homogeneous and highly jointed conditions, for the purpose of comparison. In the case of steady state, the extent of disturbed zone around the storage tunnel due to the heat production of the spent-fuel canister varies depending on the thermomechanical properties of the rock mass. In the case of transient analyses, the response of the jointed rock mass to the thermal loading after radioactive waste disposal varies significantly with time, resulting in dramatic changes in the both size and location of disturbed zone.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2008.05a
• /
• pp.231-234
• /
• 2008

This paper is concerned with the analysis of elasto-plastic stress waves by a time discontinuous variational integrator based on Hamiltonian in order to more accurate results in one dimensional dynamic problem. The proposed algorithm adopts both time-discontinuous variational integrator and space-continuous Hamiltonian so as to capture discontinuities of stress waves. This study enables to preserve total mechanical energy such as internal energy, kinetic energy and dissipative energy due to plastic deformation for long integration time. Finite element analysis of elasto-plastic stress waves is carried out in order to demonstrate the accuracy of the proposed algorithm.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.5
• /
• pp.525-533
• /
• 2010

Direct time integration method such as Newmark method is numerically performed under the assumption that continuous load function such as constant amplitude load can be treated as a discontinuous load fuction. It is due that the load can be treated as a constant value at the given time period regardless of variation of load at the time increment interval. It means the numerical results should be accompanied by the error due to approximation of load fuction. In contrast, the load function is calculated by convolution integral for the given time interval at finite element equation based on Gurtin's variation equation. Therefore. precise numerical results can be obtained by Gurtin's method because of convolution integral for the continuous load fuction curve even at the variation of load function in the given time interval. In this study, we prove that Gurtin's method can be more suitable than Newmark method in the problem of constant amplitude loading, using the numerical results for the free end of the one-dimensional rod. This study also shows that Gurtin's method is more effective in constant amplitude loading than in constant loading. The accuracy and the validity are verified by comparison between the results of in-house FORTRAN code and ADINA, a commercial software supporting Newmark method.

• Proceedings of the KSR Conference
• /
• 2000.05a
• /
• pp.471-478
• /
• 2000

For many decades, the railway was constructed using tracks with jointed rails of relatively short lengths. The joints cause many drawbacks in the track and lead to signeficant maintenance cost. so, railroad engineers became interested in eliminating joints. Continuous welded rail(CWR) track has many advantages over the conventional jointed-rail track. but, in the case of the elimination of rail joints, it may cause the track to be suddenly buckled laterally by thermal loads. In this paper, firstly, 3-D CWR track model and CWRB program exactly considering the influence of tie are developed far linear static and buckling analysis using finite element method. Characteristics of CWR track model are using 7-dofs beam element as rail, Offset technic exactly considering centroid axies difference of track components(rail, rail-pad-fastener, tie), and Thermal gradient considering thermal difference of top flange and bottom flange in rail section.. second,, Through the static and linear buckling analysis by CWRB, Influences of various track components (rail, ballast, fastener, tie and so on..) on CWR track behavior and stability was characterized.

• Journal of the Korean Geotechnical Society
• /
• v.16 no.1
• /
• pp.51-64
• /
• 2000

This paper presents a procedure of calculating a safety factor of the unsaturated slope suffering from the rainfall infiltration. The process of infiltration into a slope due to rainfall and its effect on the behavior of the soil slope are examined by using a two dimensional finite element flow-deformation coupled analysis. A factor of safety is calculated at various elapsed times after the commencement of rainfall as in the following procedure. First, stresses are estimated at each Gaussian point from the coupled finite element analysis. Then, the global stress smoothing method is applied to get a continuous stress field. Based on this stress field, a factor of safety is calculated for a specified slip surface by a stress integration scheme. Then, a search strategy is used to find out a critical slip surface which is associated with the minimum factor of safety. Some numerical examples are analyzed in order to study the effect of hydraulic conductivity on the slope stability during rain-induced infiltration. According to the results, local failure zone can be formed near the slope surface due to inhomogeneous distribution of hydraulic conductivity If the failure zone is once formed, then the region extends until a large amount of slide activates. Therefore the local failure can be neglected no longer in the stability analysis.

• Steel and Composite Structures
• /
• v.37 no.6
• /
• pp.663-678
• /
• 2020

This paper proposes a novel topology optimization method generating multiple materials for external linear plane crack structures based on the combination of IsoGeometric Analysis (IGA) and eXtended Finite Element Method (X-FEM). A so-called eXtended IsoGeometric Analysis (X-IGA) is derived for a mechanical description of a strong discontinuity state's continuous boundaries through the inherited special properties of X-FEM. In X-IGA, control points and patches play the same role with nodes and sub-domains in the finite element method. While being similar to X-FEM, enrichment functions are added to finite element approximation without any mesh generation. The geometry of structures based on basic functions of Non-Uniform Rational B-Splines (NURBS) provides accurate and reliable results. Moreover, the basis function to define the geometry becomes a systematic p-refinement to control the field approximation order without altering the geometry or its parameterization. The accuracy of analytical solutions of X-IGA for the crack problem, which is superior to a conventional X-FEM, guarantees the reliability of the optimal multi-material retrofitting against external cracks through using topology optimization. Topology optimization is applied to the minimal compliance design of two-dimensional plane linear cracked structures retrofitted by multiple distinct materials to prevent the propagation of the present crack pattern. The alternating active-phase algorithm with optimality criteria-based algorithms is employed to update design variables of element densities. Numerical results under different lengths, positions, and angles of given cracks verify the proposed method's efficiency and feasibility in using X-IGA compared to a conventional X-FEM.

• Steel and Composite Structures
• /
• v.6 no.6
• /
• pp.459-477
• /
• 2006

A theoretical research project is undertaken to develop integrated analysis and design tools for long span composite beams in modern high-rise buildings, and it aims to develop non-linear finite element models for practical design of composite beams. As the first paper in the series, this paper presents the development study as well as the calibration exercise of the proposed finite element models for simply supported composite beams. Other practical issues such as continuous composite beams, the provision of web openings for passage of building services, the partial continuity offered by the connections to columns as well as the behaviour of both unprotected and protected composite beams under fires will be reported separately. In this paper, details of the finite elements and the material models for both steel and reinforced concrete are first described, and finite element studies of composite beams with full details of test data are then presented. It should be noted that in the proposed finite element models, both steel beams and concrete slabs are modelled with two dimensional plane stress elements whose widths are assigned to be equal to the widths of concrete flanges, and the flange widths and the web thicknesses of steel beams as appropriate. Moreover, each shear connector is modelled with one horizontal spring and one vertical spring to simulate its longitudinal shear and pull-out actions based on measured load-slippage curves of push-out tests of shear connectors. The numerical results are then carefully analyzed and compared with the corresponding test results in terms of load mid-span deflection curves as well as load end-slippage curves. Other deformation characteristics of the composite beams such as stress and strain distributions across the composite cross-sections as well as distributions of shear forces and slippages in shear connectors along the beam spans are also examined in details. It is shown that the numerical results of the composite beams compare well with the test data in terms of various load-deformation characteristics along the entire deformation ranges. Hence, the proposed analysis and design tools are considered to be simple and yet effective for composite beams with practical geometrical dimensions and arrangements. Structural engineers are strongly encouraged to employ the models in their practical work to exploit the full advantages offered by composite construction.