• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.4 no.2
• /
• pp.13-24
• /
• 1984

The main part of this thesis deals with a three-dimensional transient seepage flow, discussing its fundamental conception into Which the mathematical analysis was introduced based upon the finite element method, propriety of geometrical analysis and flow movement of seepage water. As is commented in the body of the report, the saturated-unsaturated method with respect to the free water surface method is revealed to have the relative merits and demerits. Generally, the unknown sections in the form of free water surface matters in case that tunnels or pores appear. On the other hand, the applicable range is quite large and variable to the methods accordingly. Also, an example as to the seepage analysis in the rock bed section is described. An unfavorable time-consuming problem is inevitale because the complicated and widely multifarious joints and faults are responsible for the intensity of seepage flow.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.18 no.2
• /
• pp.101-108
• /
• 1982

In this paper the result of heat transfer analysis around the square heat source of interior solid by using the F. E. M is reported. Calculation for temperature distribution and each element was used by F. E. M. the solid is sub-divided into system of equal size triangular shape. These values of temperature distribution will valuable for design of jet engine and steam generator and the results gained are as follow; 1. Calculation by F. E. M is identified with the experiment. 2. Temperature distribution on the horizontal surface is $\theta$=0.698 in model 4 and the other hand $\theta$=0.401 in model 6 for X=16cm. In intermediates surface between heat source and bottom surface, the influence of L is more greater than that of height in the temperature difference. 3. Temperature distribution on the vertical surface for model 2 is resulted strong influence by K. In the case of Y=4cm is identified with $\theta$=0.0790 for K=7 and also $\theta$=0.0036 for K=0.3. In the difference of temperature distribution.

• Journal of the Korea Concrete Institute
• /
• v.18 no.1 s.91
• /
• pp.109-116
• /
• 2006

An incremental approach to predict the time dependent flexural behavior of composite girder is presented in the framework of incremental finite element method. Age dependent nature of creep, shrinkage, and maturing of elastic modulus of concrete is prescribed in the incremental tangent description of constitutive relation derived based on the first order Taylor series expansion applying to the total from of stress-strain relation. The loop phenomenon in which age dependent nature of concrete causes stress redistribution and it causes creep in turn is taken into account in the formulation through the incremental representation of constitutive relation. The developed algorithm predicts the time dependent deflections of 4.8m long two span double composite box girder subjected to shrinkage, maturing of elastic modulus, and creep initially induced by self weight. Comparison shows a good agreement between the predicted and measured results.

• Journal of Biomedical Engineering Research
• /
• v.21 no.4
• /
• pp.373-384
• /
• 2000

The goal of this paper is to develop a computational method to predict cancellous bone density distributions based upon continuum levels of volumetric strain. Volumetric strain is defined as the summation of normal strains, excluding shear strains, within an elastic range of loadings. Volumetric strain at a particular location in a cancellous structure changes with changes of the boundary conditions (prescribed displacements, tractions, and pressure). This change in the volumetric strain is postulated to predict the adaptive change in the bone apparent density. This bone remodeling theory based on volumetric strain is then used with the finite element method to compute the apparent density distribution for cancellous bone in both lumbar spine and proximal femur using an iterative algorithm, considering the dead zone of strain stimuli. The apparent density distribution of cancellous bone predicted by this method has the same pattern as experimental data reported in the literature (Wolff 1892, Keller et al. 1989, Cody et al. 1992). The resulting bone apparent density distributions predict Young's modulus and strength distributions throughout cancellous bone in agreement with the literature (Keller et al. 1989, Carter and Hayes 1977). The method was convergent and sensitive to changes in boundary conditions. Therefore, the computational algorithm of the present study appears to be a useful approach to predict the apparent density distribution of cancellous bone (i.e. a numerical approximation for Wolff's Law)

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.11 no.9
• /
• pp.3176-3183
• /
• 2010

A multi-scale fatigue life prediction methodology of composite pressure vessels subjected to multi-axial loading has been proposed in this paper. The multi-scale approach starts from the constituents, fiber, matrix and interface, leading to predict behavior of ply, laminates and eventually the composite structures. The multi-scale fatigue life prediction methodology is composed of two steps: macro stress analysis and micro mechanics of failure based on fatigue analysis. In the macro stress analysis, multi-axial fatigue loading acting at laminate is determined from finite element analysis of composite pressure vessel, and ply stresses are computed using a classical laminate theory. The micro stresses are calculated in each constituent from ply stresses using a micromechanical model. Three methods are employed in predicting fatigue life of each constituent, i.e. a maximum stress method for fiber, an equivalent stress method for multi-axially loaded matrix, and a critical plane method for the interface. A modified Goodman diagram is used to take into account the generic mean stresses. Damages from each loading cycle are accumulated using Miner's rule. Monte Carlo simulation has been performed to predict the overall fatigue life of a composite pressure vessel considering statistical distribution of material properties of each constituent, fiber volume fraction and manufacturing winding angle.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.16 no.12
• /
• pp.8654-8664
• /
• 2015

Beam-type structures with non-uniform cross sections are widely used in mechanical, architectural, and civil engineering fields. This paper deals with dynamic characteristics and vibration problems. Governing equations are first derived by using local coordinates. Their solutions are then assumed by using Galerkin's mode summation method. Bisection method is also applied in solving the determinant of the matrix which can provide natural frequencies. Whereas finite element methods adopt admissible functions satisfying only geometric boundary condition, in this study we apply Galerkin's mode summation method which uses eigen-functions satisfying both governing equations and boundary conditions. Modal analysis and experimental tests are finally performed using simply-supported beams with four different non-uniform cross-sections. Our analytical results then show good agreement with experimental ones.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.12 no.8
• /
• pp.3374-3383
• /
• 2011

In this study, numerical analyses for slamming impact phenomena have been carried out using a 2-dimensional wedge shaped structure having finite deadrise angles. Fluid is assumed incompressible and entry speed of the structure is kept constant. Geo-reconstruct(or PLIC-VOF) scheme is used for the tracking of the deforming free surface. Numerical analyses are carried out for the deadrise angles of $10^{\circ}$, $20^{\circ}$ and $30^{\circ}$. For each deadrise angle, variations are made for the grid size on the wedge bottom and for the entry speed. The magnitude and the location of impact pressure and the total drag force, which is the summation of pressure distributed at the bottom of the structure, are analyzed. Results of the analyses are compared with the results of the Dobrovol'skaya similarity solutions, the asymptotic solution based on the Wagner method and the solution of Boundary Element Method(BEM).

• Journal of Energy Engineering
• /
• v.7 no.1
• /
• pp.113-121
• /
• 1998

In natural convection flows, the fluid velocities are highly dependent on the thermal field and property variations can have a strong effect on both the flow and thermal fields. To examine the effect of property variations, at first, numerical analyses covering wide range of the Prandtl number under the same Rayleigh numbers have been carried out. Next, we have modeled the viscosity and thermal conductivity as parabolic functions of temperature and a comprehensive set of numerical solutions have been obtained to understand the effect. The Prandtl number dependence of Nusselt number is fairly strong even though the effect is still weak compared to the Rayleigh number dependence. When thermophysical properties are dependent on temperature, the flow field showed a fairly weak variation except near boundaries, whereas the temperature field is strongly affected, especially by the temperature dependent thermal conductivity.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.3A
• /
• pp.367-373
• /
• 2008

In the design of horizontally curved plate girder web panels, it is required to evaluate accurately the elastic buckling strength under pure shear. Currently, elastic shear buckling coefficients of curved web panels stiffened by transverse intermediate stiffeners are determined by assuming conservatively that straight web panels without curvature are simply supported at the juncture between the flange and web. However, depending upon the geometry and the properties of the curved plate girder, the elastically restrained support may behave rather closer to a fixed support. The buckling strength of curved girder web is much greater (maximum 38%) than that of a straight girder calculated under the assumption that all four edges are simply supported in Lee and Yoo (1999). In the present study, a series of numerical analyses based on a 3D finite element modeling is carried out to investigate the effects of geometric parameters on both the boundary condition at the juncture and the horizontal curvature of web panel, and the resulting data are quantified in a simple design equation.

• The Journal of the Acoustical Society of Korea
• /
• v.41 no.4
• /
• pp.375-388
• /
• 2022

This paper deals with the process of obtaining sound transmission loss of a cylindrical silencer lined with multi-layered poroelastic sound absorbing materials. The Biot model and the Johnson-Champoux-Allard-Lafarge (JCAL) model were used to deal with waves propagating in multi-layered poroelastic materials. The boundary conditions required for analysis of the silencer were obtained and the numerical process of finding modes was explained. A numerical experiment was conducted on the 2-layered silencer using the modes and the transmission loss converged with the first 12 modes. Finally, the mode matching method proposed in this research was validated by being compared with the results calculated from Finite Element Method (FEM) about different kinds of sound absorbing materials.