• Structural Engineering and Mechanics
• /
• v.44 no.1
• /
• pp.109-125
• /
• 2012

Post-buckling behavior of Timoshenko beams subjected to uniform temperature rising with temperature dependent physical properties are studied in this paper by using the total Lagrangian Timoshenko beam element approximation. The beam is clamped at both ends. In the case of beams with immovable ends, temperature rise causes compressible forces end therefore buckling and post-buckling phenomena occurs. It is known that post-buckling problems are geometrically nonlinear problems. Also, the material properties (Young's modulus, coefficient of thermal expansion, yield stress) are temperature dependent: That is the coefficients of the governing equations are not constant in this study. This situation suggests the physical nonlinearity of the problem. Hence, the considered problem is both geometrically and physically nonlinear. The considered highly non-linear problem is solved considering full geometric non-linearity by using incremental displacement-based finite element method in conjunction with Newton-Raphson iteration method. The beams considered in numerical examples are made of Austenitic Stainless Steel (316). The convergence studies are made. In this study, the difference between temperature dependent and independent physical properties are investigated in detail in post-buckling case. The relationships between deflections, thermal post-buckling configuration, critical buckling temperature, maximum stresses of the beams and temperature rising are illustrated in detail in post-buckling case.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03b
• /
• pp.264-267
• /
• 2008

Dual-cooled fuel with inner and outer flow channel was proposed for high burup, next generation nuclear fuel design. The annular cylinder of dual cooled fuel has higher structural strength compared to the conventional one, but also have concerns about flow induced vibration due to an additional flow of inner channel and the difference of flow velocity in between inner and outer channel. In this study, the dynamic stability of flexible, annular cylinder was evaluated according to the flow variation and compared to the that of the conventional PWR fuel rod. Centrifugal and Coriolis force by the additional flow in the inner channel were added in the dynamic equation of flexible beam in uniform, external, and axial flow. Complex eigenfrequency was calculated by the finite element method. Stability margin of annular cylinder compared to the solid cylinder and change of the dynamic characteristic are presented and discussed as a analysis results.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.145-155
• /
• 1998

A computer code has been developed for the computation of the viscous flow around a ship model with the free surface. In this code, the incompressible Reynolds-averaged Navier-Stokes equations are solved numerically by a finite difference method which employes second-order finite differences for the spatial discretization and a four-stage Runge-Kutta scheme for the temporal integration of the governing equations. For the turbulence closure, a modified version of the Baldwin-Lomax model is exploited. The location of the free surface is determined by solving the equation of the kinematic free-surface condition using the Lax-Wendroff scheme and the boundary-fitted grid is generated at each time step so that one of the grid surfaces always coincides with the free surface. An inviscid approximation of the dynamic free-surface boundary condition is applied as the boundary conditions for the velocity and pressure on the free surface. To validate the computational method and the computer code developed in the present study, the numerical computations are carried out for both Wigley parabolic hull and Series 60 $C_B=0.6$ ship model and the computational results are compared with the experimental data.

• Structural Engineering and Mechanics
• /
• v.54 no.3
• /
• pp.433-451
• /
• 2015

This paper focuses on large deflection static behavior of edge cracked simple supported beams subjected to a non-follower transversal point load at the midpoint of the beam by using the total Lagrangian Timoshenko beam element approximation. The cross section of the beam is circular. The cracked beam is modeled as an assembly of two sub-beams connected through a massless elastic rotational spring. It is known that large deflection problems are geometrically nonlinear problems. The considered highly nonlinear problem is solved considering full geometric non-linearity by using incremental displacement-based finite element method in conjunction with Newton-Raphson iteration method. There is no restriction on the magnitudes of deflections and rotations in contradistinction to von-Karman strain displacement relations of the beam. The beams considered in numerical examples are made of Aluminum. In the study, the effects of the location of crack and the depth of the crack on the non-linear static response of the beam are investigated in detail. The relationships between deflections, end rotational angles, end constraint forces, deflection configuration, Cauchy stresses of the edge-cracked beams and load rising are illustrated in detail in nonlinear case. Also, the difference between the geometrically linear and nonlinear analysis of edge-cracked beam is investigated in detail.

• Journal of Korean Port Research
• /
• v.14 no.4
• /
• pp.451-461
• /
• 2000

The design and maintenance of navigation channel and water facilities of an harbor which is located at the mouth of river or at the estuary area are difficult due to the complexity of estuarial water and sediment circulation. Effects of deepening navigable waterways, of changing coastline configurations, or of discharging dredged material to the open sea are necessary to be investigated and predicted in terms of water quality and possible physical changes to the coastal environment. A borad analysis of the transport mechanism in the estuary area was made in terms of sediment property, falling velocity, concentration and flow characteristics. In order to simulate the transport processes, a two-dimensional finite element model is developed, which includes erosion, transport and deposition mechanism of suspended sediments. Galerkin’s weighted residual method is used to solve the transient convection-diffusion equation. The fluid domain is subdivided into a series of triangular elements in which a quadratic approximation is made for suspended sediment concentration. Model could deal with a continuous aggregation by stipulating the settling velocity of the flocs in each element. The model provides suspended sediment concentration, bed shear stress, erosion versus deposition rate and bed profile at the given time step.

• Proceedings of the Korean Society of Tribologists and Lubrication Engineers Conference
• /
• 2002.10b
• /
• pp.103-104
• /
• 2002

The influence of the slip flow on the MEMS/MST based gas-lubricated proceeding bearing is investigated. Based on the modified Reynolds equation, the numerical analysis of the finite difference method was developed by applying the first order slip flow approximation. The numerical prediction of bearing performance provides the significant results concerning the slip flow effect in micro scale gas-lubricated proceeding bearing. The result indicates that the load-carrying capacity as well as the rotordynamic coefficients were significantly reduced due to the slip flow. Through this work, it is concluded that the slip flow effect could not be ignored in the micro gas-lubricated proceeding bearing.

• Journal of Institute of Control, Robotics and Systems
• /
• v.5 no.4
• /
• pp.481-488
• /
• 1999

A thermal model based on the chamber geometry of the industry-standard AST SHS200MA rapid thermal processing system has been developed for the study of thermal uniformity and process repeatability thermal model combines radiation energy transfer directly from the tungsten-halogen lamps and the steady-state thermal conducting equations. Because of the difficulties of solving partial differential equation, calculation of wafer temperature was performed by using finite-difference approximation. The proposed thermal model was verified via titanium silicidation experiments. As a result, we can conclude that the thermal model show good estimation of wafer surface temperature distribution.

• 제어로봇시스템학회:학술대회논문집
• /
• 2004.08a
• /
• pp.990-994
• /
• 2004

This paper presents a method to predict the transient characteristic of the air lubricated slider head in a hard disk drive by using optimization technique. The time dependent modified Reynolds equation based on the molecular slip flow approximation equations was used to describe the fluid flow within the air bearing and the implicit finite difference scheme is applied to calculate the pressure distribution under the slider head. The exhaustive search combined with the Broyden-Fletcher-Goldfarb-Shanno method were employed to obtain optimum design variables which are taper angle, rail width and taper length in order to keep the forces and moments acting on the slider head in dynamic equilibrium. The results show that the optimal head slider of the magnetic head has good stability characteristic that can reach the steady state within 0.5 microsecond.

• Journal of Industrial Technology
• /
• v.17
• /
• pp.199-204
• /
• 1997

In this paper, the result of application of vibration method to the orthotropic plates with free edges supported on elastic foundation and with a pair of opposite edges under axial forces is presented. Such plates represent the concrete highway slab and hybrid composite pavement of bridges. The reinforced concrete slab can be assumed as a special orthotropic plate, as a close approximation. The highway slab is supported on elastic foundation, with free boundaries. Sometimes, the pair of edges perpendicular to the traffic direction may be subject to the axial forces. The plate is subject to the concentrated load/loads, in the form of traffic loads, or the test equipments. Finite difference method is used to obtain the deflection influence surfaces needed for vibration analysis. The influence of the modulus of the foundation, the aspect ratio of the plate, and the magnitudes of the axial forces and the concentrated attached mass on the plate, under the natural frequency is thoroughly studied.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.3
• /
• pp.342-351
• /
• 1999

A jet injected normally into a cross flow has been found to have the cross section of a horseshoe shape. It occurs by a twin vortex motion in the region downstream of the jet injection. Such a flow is inherently and highly three-dimensional and numerical calculations should play an important role. The three-dimensional momentum equations with buoyancy effect and energy equation are solved to obtain the velocity distributions, center-line trajectories, cross sectional shape and entrainment. The density difference is sufficiently small, so that the Boussinesq approximation is considered to be valid. The SIMPLE algorithm is applied in a staggered grid system of a calculational domain for the numerical method.