• Nuclear Engineering and Technology
• /
• v.25 no.1
• /
• pp.51-62
• /
• 1993

A computational program ［TDET］ of the particle transport equation is developed on radiation shielding problem in two-dimensional cartesian geometry based on the discrete element method. Not like the ordinary discrete ordinates method, the quadrature set of angles is not fixed but steered by the spatially dependent angular fluxes. The angular dependence of the scattering source term in the particle transport equation is described by series expansion in spherical harmonics, and the energy dependence of the particles is considered as well. Three different benchmark tests are made for verification of TDET : For the ray effect analysis on a square absorber with a flat isotropic source, the results of TDET calculation are quite well conformed to those of MORSE-CG calculation while TDET ameliorates the ray effect more effectively than S$_{N}$ calculation. In the analysis of the streaming leakage through a narrow vacuum duct in a shield, TDET shows conspicuous and remarkable results of streaming leakage through the duct as well as MORSE-CG does, and quite better than S$_{N}$ calculation. In a realistic reactor shielding situation which treats in two cases of the isotropic scattering and of linearly anisotropic scattering with two groups of energy, TDET calculations show local ray effect between neighboring meshes compared with S$_{N}$ calculations in which the ray effect extends broadly over several meshes.eshes.

• Korea-Australia Rheology Journal
• /
• v.15 no.3
• /
• pp.131-150
• /
• 2003

A new numerical algorithm of finite element methods is presented to solve high Deborah number flow problems with geometric singularities. The steady inertialess planar 4 : 1 contraction flow is chosen for its test. As a viscoelastic constitutive equation, we have applied the globally stable (dissipative and Hadamard stable) Leonov model that can also properly accommodate important nonlinear viscoelastic phenomena. The streamline upwinding method with discrete elastic-viscous stress splitting is incorporated. New interpolation functions classified as rational interpolation, an alternative formalism to enhance numerical convergence at high Deborah number, are implemented not for the whole set of finite elements but for a few elements attached to the entrance comer, where stress singularity seems to exist. The rational interpolation scheme contains one arbitrary parameter b that controls the singular behavior of the rational functions, and its value is specified to yield the best stabilization effect. The new interpolation method raises the limit of Deborah number by 2∼5 times. Therefore on average, we can obtain convergent solution up to the Deborah number of 200 for which the comer vortex size reaches 1.6 times of the half width of the upstream reservoir. Examining spatial violation of the positive definiteness of the elastic strain tensor, we conjecture that the stabilization effect results from the peculiar behavior of rational functions identified as steep gradient on one domain boundary and linear slope on the other. Whereas the rational interpolation of both elastic strain and velocity distorts solutions significantly, it is shown that the variation of solutions incurred by rational interpolation only of the elastic strain is almost negligible. It is also verified that the rational interpolation deteriorates speed of convergence with respect to mesh refinement.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.20 no.3
• /
• pp.291-299
• /
• 2008

The present study is to estimate the effects of diffracted wave fields in a harbor and around harborentrance due to random waves, when a navigation channel is dredged in the vicinity of the a harbor entrance. The cross sections of harbor boundary are considered to be partial or full reflection in this study. The numerical simulation has been performed by the boundary element method, which is to discrete segments of pit- and harbor- boundary with the algorism of auto generated elements. The incident wave conditions are specified using discretized forms of the Mitsuyasu's frequency spectrum and directional function. The results of the present numerical simulation agreed well with those of the published experimental data. It is shown that the ratios of wave height reduction are about 20% for the case of fully reflecting boundary, and 10% for the case of partially reflecting boundary by the effect of placing a pit, respectively.

• Geomechanics and Engineering
• /
• v.12 no.6
• /
• pp.1003-1020
• /
• 2017

Estimation of groundwater inflow into underground opening is of critical importance for the design and construction of underground structures. Groundwater inflow into a pilot underground storage facility in China was estimated using analytical equations, numerical modeling and field measurement. The applicability of analytical and numerical methods was examined by comparing the estimated and measured results. Field geological investigation indicated that in local scale the high groundwater inflows are associated with the appearance of open joints, fractured zone or dykes induced by shear and/or tensile tectonic stresses. It was found that 8 groundwater inflow spots with high inflow rates account for about 82% of the total rate for the 9 caverns. On the prediction of the magnitude of groundwater inflow rate, it was found that could both (Finite Element Method) FEM and (Discrete Element Method) DEM perform better than analytical equations, due to the fact that in analytical equations simplified assumptions were adopted. However, on the prediction of the spatial distribution estimation of groundwater inflow, both analytical and numerical methods failed to predict at the present state. Nevertheless, numerical simulations would prevail over analytical methods to predict the distribution if more details in the simulations were taken into consideration.

• Proceedings of the Korean Society for Rock Mechanics Conference
• /
• 2000.09a
• /
• pp.105-115
• /
• 2000

As large underground projects such as radioactive waste disposal, hot water and heat storage, and geothermal energy become influential, the study, which consider all aspects of thermics, hydraulics and mechanics would be needed. Thermo Hydro-Mechanical coupling analysis is one of the most complex numerical technique because it should be implemented with the combined three governing equations to analyze the behavior of rock mass. In this study, finite element code, which is based on Biot's consolidation theory, was developed to analyze the thermo-hydro-mechanical coupling in continuum rock mass. To verify the implemented program, one-dimensional consolidation model under the isothermal and non-isothermal conditions was analyzed and was compared with the analytic solution. The parametric study on two-dimensional consolidation was also performed and the effects of several factors such as poisson's ratio and hydraulic anisotropy on rock mass behavior were investigated. In the future, this program would be revised to be used for analysis of general discontinuous media with incorporating discrete joint model.

• Tunnel and Underground Space
• /
• v.10 no.3
• /
• pp.355-365
• /
• 2000

As large underground projects such as radioactive waste disposal, hot water and heat storage, and geothermal energy become influential, the study, which consider all aspects of thermics, hydraulics and mechanics would be needed. Thermo-Hydro-Mechanical coupling analysis is one of the most complex numerical technique because it should be implemented with the combined three governing equations to analyze the behavior of rock mass. In this study, finite element code, which is based on Biot's consolidation theory, was developed to analyze the thermo-hydro-mechanical coupling in continuum rock mass. To verify the implemented program, one-dimensional consolidation model under the isothermal and non-isothermal conditions was analyzed and was compared with the analytic solution. The parametric study on two-dimensional consolidation was also performed and the effects of several factors such as poisson's ratio and hydraulic anisotropy on rock mass behavior were investigated. In the future, this program would be revised to be used for analysis of general discontinuous media with incorporating discrete joint model.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.38 no.3
• /
• pp.82-91
• /
• 1996

This is a basic study for the static and dynamic analysis on the elasto-plastic and elasto-viscoplastic of an axi-symmetric shell. The objective of this study was to investigate the mechanical characteristics of a nuclear reinforced concrete containment structure, which was selected as a model, by a numerical analysis using a finite element method. The structure was modeled with discrete ring elements of 8-noded isoparametric element rotating against the symmetrical axis, and the interaction between the foundation and the structure was modeled by Winkler's model. Also, the meridional tendon was modeled with 2-node truss elements, and the hoop tendon was done with point elements in two degrees of freedom. The effect of the tendon was considered without the increasement in total degree of freedom as the stiffness matrix of modeled tendon elements was assembled on the stiffness matrix of ring elements linked with the tendon. The results obtained from the analysis of an example were summarized as follows : 1. The stresses in the hoop direction on the interior and exterior surfaces of the structure were shown in changes of similar trend, and high stresses appeared on the structure wall 2. The stresses in the meridional direction on the interior and exterior surfaces were shown in change of different trend. Especially, the stresses at the junctions between the dome and the wall and between the wall and the bottom plate of the structure were very high, compared with those at other parts of the structure. 3. The stress changes in the direction of thickness on the crown of the dome were much linearly distributed. However, as the amount of tendon increased, the stresses in the upper and lower parts of the wall established with the tendon were shown stress concentration. 4. The stress changes in the direction of thickness on the center of the structure wall was linearly distributed in the all cases, and special stress due to the use of the tendon was not shown.

• Steel and Composite Structures
• /
• v.33 no.4
• /
• pp.595-614
• /
• 2019

In cold-formed steel (CFS) structures, such as trusses, transmission towers and portal frames, the use of back-to-back built-up CFS unequal angle sections are becoming increasingly popular. In such an arrangement, intermediate welds or screw fasteners are required at discrete points along the length, preventing the angle sections from buckling independently. Limited research is available in the literature on axial strength of back-to-back built-up CFS unequal angle sections. The issue is addressed herein. This paper presents an experimental investigation on both the welded and screw fastened back-to-back built-up CFS unequal angle sections under axial compression. The load-axial shortening and the load verses lateral displacement behaviour along with the deformed shapes at failure are reported. A nonlinear finite element (FE) model was then developed, which includes material non-linearity, geometric imperfections and modelling of intermediate fasteners. The FE model was validated against the experimental test results, which showed good agreement, both in terms of failure loads and deformed shapes at failure. The validated FE model was then used for the purpose of a parametric study to investigate the effect of different thicknesses, lengths and, yield stresses of steel on axial strength of back-to-back built-up CFS unequal angle sections. Five different thicknesses and seven different lengths (stub to slender columns) with two different yield stresses were investigated in the parametric study. Axial strengths obtained from the experimental tests and FE analyses were used to assess the performance of the current design guidelines as per the Direct Strength Method (DSM); obtained comparisons show that the current DSM is conservative by only 7% on average, while predicting the axial strengths of back-to-back built-up CFS unequal angle sections.

• Geomechanics and Engineering
• /
• v.11 no.1
• /
• pp.1-39
• /
• 2016

The finite element method (FEM), discrete element method (DEM), and Discontinuous deformation analysis (DDA) are among the standard numerical techniques applied in computational geo-mechanics. However, in some cases there no possibility for modelling by traditional finite analytical techniques or other mesh-based techniques. The solution presented in the current study as a completely Lagrangian and mesh-free technique is smoothed particle hydrodynamics (SPH). This method was basically applied for simulation of fluid flow by dividing the fluid into several particles. However, several researchers attempted to simulate soil-water interaction, landslides, and failure of soil by SPH method. In fact, this method is able to deal with behavior and interaction of different states of materials (liquid and solid) and multiphase soil models and their large deformations. Soil indicates different behaviors when interacting with water, structure, instrumentations, or different layers. Thus, study into these interactions using the mesh based grids has been facilitated by mesh-less SPH technique in this work. It has been revealed that the fast development, computational sophistication, and emerge of mesh-less particle modeling techniques offer solutions for problems which are not modeled by the traditional mesh-based techniques. Also it has been found that the smoothed particle hydrodynamic provides advanced techniques for simulation of soil materials as compared to the current traditional numerical methods. Besides, findings indicate that the advantages of applying this method are its high power, simplicity of concept, relative simplicity in combination of modern physics, and particularly its potential in study of large deformations and failures.

• Structural Engineering and Mechanics
• /
• v.89 no.2
• /
• pp.181-197
• /
• 2024

A numerical method is presented in this paper, for buckling analysis of thin arbitrary stiffened composite cylindrical shells under axial compression. The stiffeners can be placed inside and outside of the shell. The shell and stiffeners are operated as discrete elements, and their interactions are taking place through the compatibility conditions along their intersecting lines. The governing equations of motion are obtained based on Koiter's theory and solved by utilizing the principle of the minimum potential energy. Then, the buckling load coefficient and the critical buckling load are computed by solving characteristic equations. In this formulation, the elastic and geometric stiffness matrices of a single curved strip of the shell and stiffeners can be located anywhere within the shell element and in any direction are provided. Moreover, five stiffened composite shell specimens are made and tested under axial compression loading. The reliability of the presented method is validated by comparing its numerical results with those of commercial software, experiments, and other published numerical results. In addition, by using the ANSYS code, a 3-D finite element model that takes the exact geometric arrangement and the properties of the stiffeners and the shell into consideration is built. Finally, the effects of Poisson's ratio, shell length-to-radius ratio, shell thickness, cross-sectional area, angle, eccentricity, torsional stiffness, numbers and geometric configuration of stiffeners on the buckling of stiffened composite shells with various end conditions are computed. The results gained can be used as a meaningful benchmark for researchers to validate their analytical and numerical methods.