• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.12 no.4_1
• /
• pp.67-76
• /
• 1992

A hierarchical formulation based on p-version of the finite element method for linear elastic axisymmetric stress analysis is presented. This is accomplished by introducing additional nodal variables in the element displacement approximation on the basis of integrals of Legendre polynomials. Since the displacement approximation is hierarchical, the resulting element stiffness matrix and equivalent nodal load vectors are hierarchical also. The merits of the propoosed element are as follow: i) improved conditioning, ii) ease of joining finite elements of different polynomial order, and iii) utilizing previous solutions and computation when attempting a refinement. Numerical examples are presented to demonstrate the accuracy, efficiency, modeling convenience, robustness and overall superiority of the present formulation. The results obtained from the present formulation are also compared with those available in the literature as well as with the analytical solutions.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.39 no.3
• /
• pp.72-82
• /
• 1997

In all inelastic deformations time rate effects are always present to some degree. Whether or not their exclusion has a significant influence on the prediction of the material behaviour depends upon several factors. In the study of structural components under static loading conditions at normal temperature it is accepted that time rate effects are generally not important. However metals, especially under high temperatures, exhibit simultaneously the phenomena of creep and viscoplasticity. In this study, elastoplastic and elasto-viscoplastic models include nonlinear geometrical effects were developed and several numerical examples are also included to verify the computer programming work developed here in this work. Comparisons of the calculated results, for the elasto-viscoplastic analysis of an internally pressurised thick cylinder under plane strain condition, have shown that the model yields excellent results. The results obtained from the numerical examples for an elasto-viscoplastic analysis of the Nuclear Reinforced Concrete Containment Structure(NRCCS) subjected to an incrementally applied internal pressure were summarized as follows : 1. The steady state hoop stress distribution along the shell layer of dome and dome wall junction part of NRCCS were linearly behave and the stress in interior surfaces was larger than that in exterior. 2.However in the upper part of the wall of NRCCS the steady state hoop stress in creased linearly from its inner to outer surfaces, being the exact reverse to the previous case of dome/dome-wall junction part. 3.At the lower part of wall of NRCCS, the linear change of steady state hoop stress along its wall layer began to disturb above a certain level of load increase.

• The Journal of Korean Academy of Prosthodontics
• /
• v.44 no.1
• /
• pp.10-19
• /
• 2006

Statement of problem : The various kinds of properties of post and core may affect the stress distribution to the root of endodontically treated teeth Purpose: To evaluate the influence of various kinds of properties of post and core to the stress distribution to the root of endodontically treated teeth. Material and methods: Mandibular first premolar, prepared by general shape of post and core with gold crown, was used to two dimensional axisymmetric modeling for finite element analysis. Then property values of 8 different kinds of post and core was substituted for each. Finally, stress distribution shown areas around the root of post and core was analysed after applying 50N of vortical and oblique load. Results: 1. Stress value of oblique load was much higher than the maximum stress value of vertical load. 2. Under oblique load, very concentrated stress was located on post periapical area and variations in stress were very severe. Contrary to this, stress distribution was relatively uniform in vertical load. 3. Post materials with higher elastic modulus showed relatively more apically focused stress, and post materials with lower elastic modulus showed stress focused on cervical area on the axial wall of post. 4. Stress change according to the properties of core was shown only in the cervical area of post and below core as the higher elastic modulus, then increased in stress. 5. Post and core with medium value of elastic modulus showed relatively uniform stress distribution. Conclusions: Post materials with higher elastic modulus showed relatively more apically focused stress, and post materials with lower elastic modulus showed stress focused on cervical area on the axial wall of post. Stress change according to the properties of core was shown only in the cervical area of post and below core.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.30 no.1
• /
• pp.29-38
• /
• 2018

In order to evaluate the performance of the floating wave power, which is an axisymmetric oscillating water column type, linearized free surface boundary condition considering the influence of PTO (power takeoff) was derived and a finite element numerical model was established. Numerical experiments were carried out by varying cylinder length, skirt length, and depth of water, which are design parameters that can change the resonance of water column in cylinder and heave resonance of the float, which is considered to affect the power generation efficiency. Finally, the basic data necessary for the optimum design of the power generation system were obtained. As a result, the efficiency of the power generation system is dominated by the heave motion resonance of the float rather than the water column resonance in the cylinder, and the resonance condition for the heave motion can be changed efficiently by attaching the skirt to the outside of the buoy.

• Structural Engineering and Mechanics
• /
• v.84 no.5
• /
• pp.633-649
• /
• 2022

In this study, the effect of circle tunnel on the force distribution around underground rectangular gallery was investigated using theoretical approach, experimental test and Particle flow code simulation (PFC). Gypsum model with dimension of 1500×1500 mm was built. Tensile strength of material was 1 MPa. Dimension of central gallery was 100 mm×200 mm and diameter of adjacent tunnel in its right side was 20 mm, 40 mm and 60 mm. Horizontal distance between tunnel wall and gallery edge were 25, 50, 75, 100 and 125 mm. using beam theory, the effect of tunnel diameter and distance between tunnel and gallery on the induced force around gallery was analyzed. In the laboratory test, the rate of loading displacement was set to 0.05 millimeter per minute. Also sensitivity analysis has been done. Using PFC2D, interaction between tunnel and gallery was simulated and its results were compared with experimental and theoretical analysis. The results show that the tensile force concentration has maximum value in center of the rectangular space. The tensile force concentration at the right side of the axisymmetric line of gallery has more than its value in the left side of the galleries axisymmetric line. The tensile force concentration was decreased by increasing the distance between tunnel and rectangular space. In whole of the configurations, the angles of micro cracks fluctuated between 75 and 105 degrees, which mean that the variations of tunnel situation have not any influence on the fracture angle.

• The Journal of Korean Academy of Prosthodontics
• /
• v.38 no.3
• /
• pp.348-357
• /
• 2000

The purpose of this study is to examine the effect of partial osseointegration situation on bone loading patterns around two different free-standing screw shaped implants (Nobel Biocare, Gothenburg, Sweden and Degussa-Huls, Hanau, German). Two dimensional axisymmetric Finite element models of two implants(10mm length and 4mm diameter) were created according to different bone quantity, quality and osseointegration ratio in maxilla and mandible bone. At the same time uni-cortical and hi-cortical fixation were analyzed. Generally, full bond case showed less stress than partial bond case in overall area and mandibular model showed less amount of stress than that of maxilla model. Maximum stress of the Branemark implant is higher than that of ANKYLOS regardless of bonding ratio at crestal and apex region. However, more stress concentration was noted in ANKYLOS implant at screw body area especially in mandible. The effect of bicortical fixation on crestal bone stress reduction is dramatical in mandible however, there was no significant effect in maxillary case. The effect of partial bond on stress distribution was more significant at screw body and apex region than in crestal region. Partial bond cases demonstrated greater stress accumulation in trabecular bone than cortical bone. It is concluded that the more accurate model of implant and bone which affects stress and strain distribution is needed to mimic in vivo behavior of implants.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2006.10a
• /
• pp.291-311
• /
• 2006

The acceleration of consolidation by stone columns was mostly analysed within the framework of a basic unit cell model (i.e. a cylindrical soil body around a column). A method of converting the axisymmetric unit cell into the equivalent plane-strain model would be required for two-dimensional numerical modelling of multi-column field applications. This paper proposes two practical simplified conversion methods to obtain the equivalent plane-strain model of the unit cell, and investigates their applicability to multi-column reinforced ground. In the first conversion method, the soil permeability is matched according to an analytical equation, whereas in the second method, the column width is matched based on the equivalence of column area. The validity of these methods is tested by comparison with the numerical results of unit-cell simulations and with the field data from an embankment case history. The results show that for the case of linear-elastic material modelling, both methods produce reasonably accurate long-term consolidation settlements, whereas for the case of elasto-plastic material modelling, the second method is preferable as the first one gives erroneously lower long-term settlements, where plastic yielding of stone column are ignored.

• Proceedings of the KSME Conference
• /
• 2007.05a
• /
• pp.683-688
• /
• 2007

Surface heat treatment of SM45C round bar by diode laser was simulated to find it's condition by using commercial finite element code MARC. Due to axisymmetric geometry, a quarter of model for SM45C round bar was considered and user subroutines were applied to boundary condition for the heat transfer. Material properties such as conductivity, specific heat and mass density were given as a function of temperature. Rotation speed of round bar and feed rate of beam were considered to design heat source model. Shape parameter values of heat source were determined by beam profile. As results, Three dimensional heat source model for diode laser beam conditions of surface hardening has been designed by the comparison between the finite element analysis results and experimental data on SM45C round bar. Diode laser surface hardening for SM45C round bar was successfully simulated and it should be useful to determine optimal heat treatment condition.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2011.11a
• /
• pp.321-324
• /
• 2011

A computational analysis of nozzle flow characteristics and plume structure using Reynolds-averaged Navier-Stokes equations with $k-{\omega}$ SST turbulence model was conducted to examine performance of the supersonic nozzle employed in a small liquid-rocket engine for ground firing test. Computed results and experimental outcome of 2-D converging-diverging nozzle flow were compared for verifying the computational capability as well as the turbulence model validity. Numerical computations of 2-D axisymmetric nozzle flow was carried out with the selected model. As a result, flow separation with backflow appeared around the nozzle exit. This investigation was reported as a background data for the optimal nozzle design of small liquid-propellant rocket engine for ground test.

• Journal of computational fluids engineering
• /
• v.11 no.2 s.33
• /
• pp.8-18
• /
• 2006

A CFD benchmark calculation for a steam blowdown test was performed for 30 seconds to develop the methodology of numerical analysis for the thermal mixing between steam and subcooled water. In the CFD analysis, the grid model simulating the sparger and the IRWST pool were developed by the axisymmetric condition and then the steam condensation phenomena by a direct contact was modelled by the so-called condensation region model. Thermal mixing phenomenon in the subcooled water tank was treated as an incompressible flow, a free surface flow between the air and the water, a turbulent flow, and a buoyancy flow. The comparison of the CFD results with the test data showed a good agreement as a whole, but a small temperature difference was locally found at some locations. The commercial CFD code of CFX4.4 together with the condensation region model can simulate the thermal mixing behavior reasonably well when a sufficient number of mesh distribution and a proper numerical method are adopted.