• Structural Engineering and Mechanics
• /
• v.30 no.1
• /
• pp.21-36
• /
• 2008

The static response of a finite beam resting on a tensionless Pasternak foundation and subjected to a concentrated vertical load is assessed in this study. The concentrated vertical load may be applied at the center of the beam, or it may be offset from the center. The tensionless character of the foundation results in the creation of lift-off regions between the beam and the foundation. An analytical/ numerical solution is obtained from the governing equations of the contact and lift-off regions to determine the extent of the contact region. Although there is no nonlinear term in the equations, the problem shows a nonlinear character since the contact region is not known in advance. Due to that nonlinearity, the essentials of the problem (the coordinates of the lift-off points) are calculated numerically using the Newton-Raphson technique. The numerical results are presented in figures to illustrate the behaviours of the free-free and pinned-pinned beams under symmetric or asymmetric loading. The figures illustrate the effects of the shear foundation parameter and the symmetric and asymmetric loading options on the variation of the contact lengths and the displacement of the beam.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.3
• /
• pp.199-209
• /
• 2009

This paper deals with the three-dimensional finite element analysis of failure behavior of UHPFRC I-beam under monotonic load. Different from the constitutive law of normal and high strength concrete, an elastic-plastic fracture model that considers the tensile strain hardening is proposed to describe the material properties of UHPFRC. A multi-directional fixed crack criterion with tensile strain hardening is defined in the tensile region, and Drucker-Prager criterion with an associated flow rule is adopted in the compressive region. The influence of span, prestressing force and section on the behavior of UHPFRC I-beam are investigated. The comparison of the numerical results with the test results indicates a good agreement.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.23 no.7 s.166
• /
• pp.1120-1128
• /
• 1999

An elastic-viscoplastic finite element analysis is performed to investigate detailed growth behavior of creep cracks and the numerical results are compared with experimental results. In Cr-Mo steel stress fields obtained from the crack growth method by mesh translation were compared with both cases that the secondary creep rate is only used as creep material property and the primary creep rate is included. Analytical stress fields, Riedel-Rice(RR) field, Hart-Hui-Riedel(HR) field and Prime(named in here) field, and the results obtained by numerical method were evaluated in details. Time vs. stress at crack tip was showed and crack tip stress fields were plotted. These results were compared with analytical stress fields. There is no difference of stress distribution at remote region between the case of 1st creep rate+2nd creep rate and the case of 2nd creep rate only. In case of slow velocity of crack growth, the effect of 1st creep rate is larger than the one of fast crack growth rate. Stress fields at crack tip region we, in order, Prime field, HR field and RR field from crack tip.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.5 s.236
• /
• pp.568-576
• /
• 2005

A 3-D finite element model combined with a volume tracking method is presented in this work to simulate the mold filling for casting processes. Especially, the analysis involves an adaptive grid method that is created under a criterion of element categorization of filling states and locations in the total region at each time step. By using an adaptive grid wherein the elements, finer than those in internal and external regions, are distributed at the surface region through refinement and coarsening procedures, a more efficient analysis of transient fluid flow with free surface is achieved. Adaptive grid based on VOF method is developed in tetrahedral element system. Through a 3-D analysis of the benchmark test of the casting process, the efficiency of the proposed adaptive grid method is verified. Developed FE code is applied to a typical industrial part of the casting process such as aluminum road wheel.

• Journal of the Korean Society for Precision Engineering
• /
• v.32 no.6
• /
• pp.509-515
• /
• 2015

Recent trends to reduce the size of mobile electronics products have driven miniaturization of various components, including screw parts for assembling components. Considering that the size reduction of screws may degenerate their joining capabilities, the size reduction should not be limited to the thread region but should be extended to its head region. The screw head is usually manufactured by forging in which a profiled punch presses a billet so that plastic deformation occurs to form the desired shape. In this study, finite element (FE) analysis was performed to simulate the forging process of a subminiature screw; a screw head of 1.7 mm diameter is formed out of a 0.82 mm diameter billet. The FE analysis result indicates that this severe forging condition leads to a generation of folding defects. FE analyses were further performed to find appropriate punch design parameters that minimize the amount of folding defects.

• Journal of the Korean Association of Oral and Maxillofacial Surgeons
• /
• v.34 no.2
• /
• pp.151-156
• /
• 2008

Different kinds of forces can be applied to the biological tissue. The analysis of the applied force is highly important to explain the mechanism of cellular response. In this study, the applied force to the collagen gel was analyzed by the finite elements analysis. The model received two different kinds of static force (compression and tension). The force range was 50g to 400g. In results, von Mises stress was concentrated in the peripheral region in the compression model. It was concentrated in the central area in the tension model. However, the compressive force was high in the peripheral area of the compression model and the tensional force was also high in the same area of the tension model. In conclusion, the applied force could be different to the region and it should be considered in the experiment to analyze the effects of the mechanical force on the cells.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2004.10a
• /
• pp.75-78
• /
• 2004

A mesh generation algorithm adapted to the mesh density map using the Delaunay mesh generation technique is developed. In the finite element analyses of the forging processes, the numerical error increases as the process goes on because of discrete property of the finite elements or severe distortion of elements. Especially, in the region where stresses and strains are concentrated, the numerical discretization error will be highly increased. However, it is too time consuming to use a uniformly fine mesh in the whole domain to reduce the expected numerical error. Therefore, it is necessary to construct locally refined mesh at the region where the error is concentrated such as at the die corner. In this study, the point insertion algorithm is used and the mesh size is controlled by moving nodes to optimized positions according to a mesh density map constructed with a posteriori error estimation. An optimization technique is adopted to obtain a good position of nodes. And optimized smoothing techniques are also adopted to have smooth distribution of the mesh and improve the mesh element quality.

• Ocean Systems Engineering
• /
• v.2 no.3
• /
• pp.229-238
• /
• 2012

We discuss the bi - stability that is possibly exhibited by a liquid free surface in a parametrically - driven two-dimensional (2D) rectangular tank with finite liquid depth. Following the method of adaptive mode ordering, assuming two dominant modes and retaining polynomial nonlinearities up to third-order, a nonlinear finite-dimensional nonlinear modal system approximation is obtained. A "continuation method" of nonlinear dynamics is then used in order to elicit efficiently the instability boundary in parameters' space and to predict how steady surface elevation changes as the frequency and/or the amplitude of excitation are varied. Results are compared against those of the linear version of the system (that is a Mathieu-type model) and furthermore, against an intermediate model also derived with formal mode ordering, that is based on a second - order ordinary differential equation having nonlinearities due to products of elevation with elevation velocity or acceleration. The investigation verifies that, in parameters space, there must be a region, inside the quiescent region, where liquid surface instability is exhibited. There, behaviour depends on initial conditions and a wave form would be realised only if the free surface was substantially disturbed initially.

• Journal of Mechanical Science and Technology
• /
• v.20 no.10
• /
• pp.1741-1752
• /
• 2006

A combined Streamline Upwind Petrov-Galerkin method (SUPG) and segregated finite element algorithm for solving conjugate heat transfer problems where heat conduction in a solid is coupled with heat convection in viscous fluid flow is presented. The Streamline Upwind Petrov-Galerkin method is used for the analysis of viscous thermal flow in the fluid region, while the analysis of heat conduction in solid region is performed by the Galerkin method. The method uses the three-node triangular element with equal-order interpolation functions for all the variables of the velocity components, the pressure and the temperature. The main advantage of the presented method is to consistently couple heat transfer along the fluid-solid interface. Four test cases, which are the conjugate Couette flow problem in parallel plate channel, the counter-flow in heat exchanger, the conjugate natural convection in a square cavity with a conducting wall, and the conjugate natural convection and conduction from heated cylinder in square cavity, are selected to evaluate efficiency of the presented method.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.26 no.2
• /
• pp.108-120
• /
• 2022

We considered qualitative behaviour of the generalization of Einstein's model of Brownian motion when the key parameter of the time interval of free jump degenerates. Fluids will be characterised by number of particles per unit volume (density of fluid) at point of observation. Degeneration of the phenomenon manifests in two scenarios: a) flow of the fluid, which is highly dispersing like a non-dense gas and b) flow of fluid far away from the source of flow, when the velocity of the flow is incomparably smaller than the gradient of the density. First, we will show that both types of flows can be modeled using the Einstein paradigm. We will investigate the question: What features will particle flow exhibit if the time interval of the free jump is inverse proportional to the density and its gradient ? We will show that in this scenario, the flow exhibits localization property, namely: if at some moment of time t₀ in the region, the gradient of the density or density itself is equal to zero, then for some T during time interval [t₀, t₀ + T] there is no flow in the region. This directly links to Barenblatt's finite speed of propagation property for the degenerate equation. The method of the proof is very different from Barenblatt's method and based on the application of Ladyzhenskaya - De Giorgi iterative scheme and Vespri - Tedeev technique. From PDE point of view it assumed that solution exists in appropriate Sobolev type of space.