• Computational Structural Engineering
• /
• v.8 no.1
• /
• pp.173-186
• /
• 1995

It is well known that structures constructed by bonding two or more materials and then subjected to temperature change experience thermal stress. This stress results from thermal expansion mismatch of materials. The present paper derives formulas for the stresses in a bimaterial axisymmetric disk which is subjected to a uniform temperature change. First, an approximate solution following strength-of-materials principles is developed. However, the strength-of-materials solution has difficulty in predicting both the peak value of interfacial stresses and its associated distribution. Next, a solution consistent with the theory of elasticity is developed by way of an eigenfunction expansion approach. The eigenfunction analysis is compared with finite element stress analysis results for a specific numerical example. Finite element analysis results show that the interfacial stresses are adequately predicted by eigenfunction solution. Therefore, the method developed in this paper will be useful in determination of the interfacial stress state.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.22 no.12
• /
• pp.1726-1734
• /
• 1998

This paper focuses on the unsteady close-contact melting phenomenon occurring between a phase change material kept at its melting temperature and a flat surface on which constant heat flux is imposed. Based on the same simplifications and framework of analysis as the case of constant surface temperature, an approximate analytical solution which depends only on the liquid-to-solid density ratio is successfully derived. In order to keep consistency with the known solution procedure, both the dimensionless wall heat flux and the Stefan number are properly redefined. The obtained solution proves to agree quite well with the published numerical data and to be capable of resolving the fundamental features of unsteady close-contact melting, especially in the presence of the solid-liquid density difference. The density ratio directly affects the film growth rate and the initial value of solid descending velocity, thereby controlling the duration of unsteady process. The effects of other parameters can be evaluated readily from the steady solution which is implied in the normalized result. Since the dimensionless surface temperature for the present boundary condition increases from zero to unity along the evolution path of the liquid film thickness, the unsteady process lasts longer than that for the case of isothermal heating.

• Proceedings of the Korean Society For Composite Materials Conference
• /
• 2004.04a
• /
• pp.223-227
• /
• 2004

The pultruded structural shapes are usually composed of thin-walled plate elements. Because the composite material has relatively low elastic moduli, the design of pultruded compression members may not be governed by the material strength limit state but by the stability limit state such as the local buckling or the global buckling. Therefore, the stability limit state must be checked to design pultruded columns. In this research, the local buckling analysis of pultruded I-shape column was conducted for various composite materials using the closed-form solution. To establish the design guidelines for the local buckling of pultruded I-shape compression members, the simplified form of equation to find the local buckling coefficient of pultruded I-shape column was proposed as a function of mechanical properties and the width ratio of plate components using the results obtainde by the closed-form solution. In order to verify the validity of proposed solution, the results obtained by the proposed approximate solution were compared with those of the closed-form solution and the experimental results.

• The Journal of the Acoustical Society of Korea
• /
• v.18 no.3
• /
• pp.73-78
• /
• 1999

The paper describes a theoretical study on the flexural vibration of an elastic flat bar with periodically nonuniform material properties. The approximate solution of the natura1 frequency and mode shape has been obtained using the perturbation technique for sinusoidal modulation of the flexural rigidify and mass density. The numerical solution obtained by using the finite element method verifies the trend of the approximate solution. It appears that distributed vibrations exist in the low modes, and this approach can be extended to the vibration analysis of the p1ate in the flat panel speaker.

• Proceedings of the Korean Society of Computational and Applied Mathematics Conference
• /
• 2003.09a
• /
• pp.4.2-4
• /
• 2003

In this paper, we describe how to approximate the solutions of partial differential equations by bicubic spline functions whose interpolation errors at non-grid points are smaller in general than those by linear interpolations of the original solution at grid points. We show that the bicubic spline function can be represented exactly or approximately by a fuzzy system, and that the resulting fuzzy rule table shows the contours of the solution function. Thus, the fuzzy rule table is identified as a digital image and the contours in the rule table provide approximate contours of the solution of partial differential equations. Several illustrative examples are included.

• Bulletin of the Korean Mathematical Society
• /
• v.55 no.3
• /
• pp.881-898
• /
• 2018

This article deals with implementation of a high-order finite difference scheme for numerical solution of the incompressible Navier-Stokes equations on curvilinear grids. The numerical scheme is based on pseudo-compressibility approach. A fifth-order upwind compact scheme is used to approximate the inviscid fluxes while the discretization of metric and viscous terms is accomplished using sixth-order central compact scheme. An implicit Euler method is used for discretization of the pseudo-time derivative to obtain the steady-state solution. The resulting block tridiagonal matrix system is solved by approximate factorization based alternating direction implicit scheme (AF-ADI) which consists of an alternate sweep in each direction for every pseudo-time step. The convergence and efficiency of the method are evaluated by solving some 2D benchmark problems. Finally, computed results are compared with numerical results in the literature and a good agreement is observed.

• Journal of the Korean Society of Manufacturing Process Engineers
• /
• v.19 no.5
• /
• pp.27-37
• /
• 2020

Structural optimization problems with discrete design variables require more function calculations (or finite element analyses) than those in the continuous design space. In this study, a method to find an optimal solution in the discrete design of the truss structure is presented, reducing the number of function calculations. Because a continuous optimal solution is the Karush-Kuhn-Tucker point that satisfies the optimality condition, it is assumed that the discrete optimal solution is around the continuous optimum. Then, response values such as weight, displacement, and stress are predicted using approximate models-referred to as hybrid metamodels-within specified design ranges. The discrete design method using the hybrid metamodels is used as a post-process of the continuous optimization process. Standard truss design problems of 10-bar, 25-bar, 15-bar, and 52-bar are solved to show the usefulness of this method. The results are compared with those of existing methods.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.6 no.4
• /
• pp.453-462
• /
• 1994

An unsteady quasi one-dimensional model of momentum, heat and mass transfer in a falling film of a vertical plate absorber which is cooled by air was developed using the integral method. Energy conservation of the absorber wall is considered in the model. The model can predict absorption rate, film thickness and mean velocity as well as concentration and temperature profiles. Predictions of steady state temperature and concentration profiles for LiBr/water system for constant wall temperature condition are in good agreement with the two-dimensional finite difference method solutions. Effects of operating conditions, such as convective heat transfer coefficient between the cooling air and the absorber wall, cooling air temperature and film thickness at inlet, on absorption rate of water vapor into LiBr/water solution were shown.

• Journal of applied mathematics & informatics
• /
• v.7 no.1
• /
• pp.183-193
• /
• 2000

In this paper we use measure theory to solve a wide range of second-order boundary value ordinary differential equations. First, we transform the problem to a first order system of ordinary differential equations(ODE's)and then define an optimization problem related to it. The new problem in modified into one consisting of the minimization of a linear functional over a set of Radon measures; the optimal measure is then approximated by a finite combination of atomic measures and the problem converted approximatly to a finite-dimensional linear programming problem. The solution to this problem is used to construct the approximate solution of the original problem. Finally we get the error functional E(we define in this paper) for the approximate solution of the ODE's problem.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.22 no.3
• /
• pp.163-177
• /
• 2018

This paper successfully applies the modified Adomian decomposition method to find the approximate solutions of the Caputo fractional integro-differential equations. The reliability of the method and reduction in the size of the computational work give this method a wider applicability. Also, the behavior of the solution can be formally determined by analytical approximation. Moreover, we proved the existence and uniqueness results and convergence of the solution. Finally, an example is included to demonstrate the validity and applicability of the proposed technique.