• Steel and Composite Structures
• /
• v.19 no.3
• /
• pp.679-695
• /
• 2015

A numerical solution using finite difference method to evaluate the thermal buckling of simply supported FGM plate with variable thickness is presented in this research. First, the governing differential equation of thermal stability under uniform temperature through the plate thickness is derived. Then, the governing equation has been solved using finite difference method. After validating the presented numerical method with the analytical solution, the finite difference formulation has been extended in order to include variable thickness. The accuracy of the finite difference method for variable thickness plate has been also compared with the literature where a good agreement has been found. Furthermore, a parametric study has been conducted to analyze the effect of material and geometric parameters on the thermal buckling resistance of the FGM plates. It was found that the thickness variation affects isotropic plates a bit more than FGM plates.

• 제어로봇시스템학회:학술대회논문집
• /
• 1994.10a
• /
• pp.156-161
• /
• 1994

In this paper, we consider an optimal control problem of a nonlinear stochastic system. Dynamic programming approach is employed for the formulation of a stochastic optimal control problem. As an optimality condition, dynamic programming equation so called the Bellman equation is obtained, which seldom yields an analytical solution, even very difficult to solve numerically. We obtain the numerical solution of the Bellman equation using an algorithm based on the finite difference approximation and the contraction mapping method. Optimal controls are constructed through the solution process of the Bellman equation. We also construct a test case in order to investigate the actual performance of the algorithm.

• Journal of the Society of Naval Architects of Korea
• /
• v.44 no.2 s.152
• /
• pp.130-138
• /
• 2007

A coupled variational equation for fluid-structure interaction (FSI) problems is derived from a steady state Navier-Stokes equation for incompressible Newtonian fluid and an equilibrium equation for geometrically nonlinear structures. For a fully coupled FSI formulation, between fluid and structures, a traction continuity condition is considered at interfaces where a no-slip condition is imposed. Under total Lagrange formulation in the structural domain, finite rotations are well described by using the second Piola-Kirchhoff stress and Green-Lagrange strain tensors. An adjoint shape design sensitivity analysis (DSA) method based on material derivative approach is applied to the FSI problem to develop a shape design optimization method. Demonstrating some numerical examples, the accuracy and efficiency of the developed DSA method is verified in comparison with finite difference sensitivity. Also, for the FSI problems, a shape design optimization is performed to obtain a maximal stiffness structure satisfying an allowable volume constraint.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.13 no.5
• /
• pp.2403-2407
• /
• 2012

Finite number of pile elements are considered in load transfer method. And section force and movement of each pile element are computed by considering compatibilities between pile displacement and the load transfer along a pile and between displacement and resistance at the tip of the pile. For the conventional load transfer method, large amount of computations due to iterations are needed. Formulation of finite difference equation from the differential equation which depicts pile behavior under axial loading was accomplished in order to simplify the computation for obtaining pile section forces and displacements. By comparing the results between the simplified computation method and the reported data, there was no difference between the two results.

• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 2000.10a
• /
• pp.81-88
• /
• 2000

The elastic wave equation is solved using the finite-difference method in 3D space to simulate the seismic wave propagation. It is based on the velocity-stress formulation of the equation of motion on a staggered grid. The nonreflecting boundary conditions are used to attenuate the wave field close to the numerical boundary. To satisfy the stress-free conditions at the free-surface boundary, a new formulation combining the zero-stress formalism with the vacuum one is applied. The effective media parameters are employed to satisfy the traction continuity condition across the media interface. With use of the moment-tensor components, the wide range of source mechanism parameters can be specified. The numerical experiments are carried out in order to test the applicability and accuracy of this scheme and to understand the fundamental features of the wave propagation under the generalized elastic media structure. Computational results show that the scheme is sufficiently accurate for modeling wave propagation in 3D elastic media and generates all the possible phases appropriately in under the given heterogeneous velocity structure. Also the characteristics of the ground motion in an sedimentary basin such as the amplification, trapping, and focusing of the elastic wave energy are well represented. These results demonstrate the use of this simulation method will be helpful for modeling the ground motion of seismological and engineering purpose like earthquake hazard assessment, seismic design, city planning, and etc..

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.24 no.6 s.177
• /
• pp.1520-1528
• /
• 2000

This study presents a method for determining bearing stiffness and damping coefficients of air-lubricated slider bearing, and shows influences of air-bearing surface geometry(recess depth, crown an d pivot location) on flying attitude and dynamic characteristics. To derive the dynamic lubrication equation, the perturbation method is applied to the generalized lubrication equation which based on linearized Boltzmann equation. The generalized lubrication equation and the dynamic lubrication equation are converted to a control volume formulation, and then, the static and dynamic pressure distributions are calculated by finite difference method. The recess depth and crown of the slider show significantly influence on flying attitude and dynamic characteristics comparing with those of pivot location.

• Journal of KSNVE
• /
• v.6 no.5
• /
• pp.607-616
• /
• 1996

A combined computational fluid dynamics(CFD)-Kirchhoff method is presented for predicting high-speed impulsive noise generated by a hovering blade. Two types of Kirchhoff integral formula are used; one for the classical linear Kirchhoff formulation and the other for the nonlinear Kirchhoff formulation. An Euler finite difference solver is solved first to obtain the flow field close to the blade, and then this flow field is used as an input to a Kirchhoff formulation to predict the acoustic far-field. These formulas are used at Mach numbers of 0.90 and 0.95 to investigate the effectiveness of the linear and nonlinear Kirchhoff formulas for delocalized flow. During these calculiations, the retarded time equation is also carefully examined, in particular, for the cases of the control surface located outside of the sonic cylinder, where multiple roots are obtained. Predicted results of acoustic far-field pressure with the linear Kirchhoff formulation agree well with experimental data when the control surface is at the certain location(R=1.46), but the correlation is getting worse before or after this specific location of the control surface due to the delocalized nonlinear aerodynamic flow field. Calculations based on the nonlinear Kirchhoff equation using a linear sonic cylinder as a control surface show a reasonable agreement with experimental data in negative amplitudes for both tip Mach numbers of 0.90 and 0.95, except some computational integration problems over a shock. This concliudes that a nonlinear formulation is necessary if the control surface is close to the blade and the flow is delocalized.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.13 no.4
• /
• pp.1895-1899
• /
• 2012

Formulation of finite difference method for analyzing consolidation were carried out. It can be seen that the differences in settlement with time obtained by FDM and Terzaghi method are diminished by fine discretization of time increment. Excess pore pressures predicted by the derived finite difference equation were same as those calculated by Olson's method. Predicted time-settlement behavior from the derived finite difference method were almost same as those calculated by Terzaghi's method and Olson's method. Analysis results obtained from the assumed multi-step time dependent loading are thought to be reasonable.

• Tribology and Lubricants
• /
• v.18 no.6
• /
• pp.402-410
• /
• 2002

In this study, a hydrodynamic analysis of the reciprocating compressor crankshaft considering a finite bearing modelling of the journal bearings used in small refrigeration compressors is performed. In the problem formulation of the compression mechanism dynamics, all corresponding hydrodynamic forces and moments are considered using the finite bearing analysis in order to determine the crankshaft trajectory at each step. The solution of the Reynolds' equation is determined numerically using a finite difference method and a Newton-Raphson procedure was employed in solving the dynamic equations of the crankshaft. The crankshaft orbits fur the finite bearing model and short bearing theory were used to compare the effect of the hydrodynamic farces of the journal bearings on the dynamic and lubrication characteristics of the crankshaft-journal bearing system. Results show that the finite bearing model for the journal bearings must be considered in calculating for the accurate dynamic characteristics of the reciprocating compressor crankshaft.

• Journal of KSNVE
• /
• v.10 no.1
• /
• pp.97-106
• /
• 2000

The dynamic behavior of crankshaft-bearing system in scroll compressor has been investigated using the combined methodologies of finite elements and transfer matrices. The finite element formulation is proposed including the field element for a shaft section and the point element at balancer weight locations, bearing locations, etc., whereas the conventional method is used with the elements. The Houbolt method is used to consider the time march for the integration of the system equations. The linear stiffness and damping coefficients are calculated for a finite cylindrical fluid-film bearing by solving the Reynolds equation, using finite difference method. The orbital response of crankshaft supported on the linear bearing model is obtained, considering balancer weights of motor rotor. And, the steady state displacement of crankshaft are compared with a variation in balancer weight. The loci of crankshaft at bearing locations are composed of the synchronous whirl component and the non-synchronous whirl component.