• Steel and Composite Structures
• /
• 제19권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년도 Proceedings of the Korea Automatic Control Conference, 9th (KACC) ; Taejeon, Korea; 17-20 Oct. 1994
• /
• 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.

• 대한조선학회논문집
• /
• 제44권2호
• /
• 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.

• 한국산학기술학회논문지
• /
• 제13권5호
• /
• pp.2403-2407
• /
• 2012

압축하중을 받는 말뚝의 하중-침하관계를 구하기 위한 하중전이법에서는 말뚝을 유한개의 요소로 분할하고 각각의 요소들에 대하여 말뚝주면변위와 주면 전이하중 그리고 말뚝선단변위와 선단저항력 사이의 적합조건으로부터 말뚝 요소의 단면력과 변위를 계산하게 된다. 통상적인 하중전이법에 의할 경우 모든 말뚝요소들에 대하여 하중-변위 사이의 적합조건을 만족하도록 반복계산을 수행해야 하므로 계산량이 과다한 단점이 있다. 본 연구에서는 압축하중을 받는 말뚝에 대한 미분방정식을 유한차분식으로 정식화하여 말뚝 바닥요소를 제외한 요소에 대해서는 반복계산 없이 직접 말뚝 단면력과 변위량을 계산할 수 있는 간편 계산법을 도입하였다. 간편 계산법에 의한 계산결과를 기존의 해석자료와 비교하여 보았는데 하중-침하관계와 말뚝 단면력 분포가 일치함을 알 수 있었다.

• 한국지진공학회:학술대회논문집
• /
• 한국지진공학회 2000년도 추계 학술발표회 논문집 Proceedings of EESK Conference-Fall 2000
• /
• 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..

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

• 소음진동
• /
• 제6권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.

• 한국산학기술학회논문지
• /
• 제13권4호
• /
• pp.1895-1899
• /
• 2012

유한차분법을 적용한 압밀해석을 수행하였는데 순간하중이 재하되는 경우 유한차분법에 의해 예측되는 시간별 침하량과 Terzaghi 방법에 의한 침하량 사이의 차이는 시간격자간격을 충분히 작게 하여 해결할 수 있음을 알 수 있었다. 점증하중에 대한 압밀해석을 위한 유한차분식을 유도하였는데 해석결과에 따른 과잉간극수압의 분포가 Olson의 이론해와 일치하였다. 점증하중이 작용하는 경우에 대해 예측한 시간-침하거동에 있어 유도된 유한차분식에 의한 결과와 Terzaghi 및 Olson 에 의한 결과 또한 거의 일치하였다. 다단 점증하중에 대한 해석결과 또한 신뢰성이 높은 것으로 보인다.

• Tribology and Lubricants
• /
• 제18권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.

• 소음진동
• /
• 제10권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.