• Structural Engineering and Mechanics
• /
• 제30권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.

• 한국전산구조공학회논문집
• /
• 제22권3호
• /
• pp.199-209
• /
• 2009

본 논문은 단조하중을 받는 초고강도 섬유보강 콘크리트 I형보의 파괴거동에 대한 3차원 유한요소해석을 수행하였다. 보통 또는 고강도 콘크리트의 구성방정식과 달리 초고강도 섬유보강 콘크리트의 재료적 특성을 나타내기 위해 인장변형률 경화관계를 고려한 탄소성 파괴역학 모델을 제안하였다. 인장영역에서는 인장경화 변형률을 고려한 다차원적 균열기준을 정의하였고, 압축영역에서는 associated flow rule을 고려한 Drucker-Prager기준을 채택하여 해석을 수행하였다. UHPFRCI형보의 지간, 프리스트레스 하중 및 단면의 영향에 관한 수치해석 결과를 실험 거동와 비교한 결과 정확한 해석 결과를 보여주었다.

• 대한기계학회논문집A
• /
• 제23권7호
• /
• 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.

• 대한기계학회논문집B
• /
• 제29권5호
• /
• 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.

• 한국정밀공학회지
• /
• 제32권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
• /
• 제34권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.

• 한국소성가공학회:학술대회논문집
• /
• 한국소성가공학회 2004년도 추계학술대회논문집
• /
• 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
• /
• 제2권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
• /
• 제20권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
• /
• 제26권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.