• Structural Engineering and Mechanics
• /
• v.28 no.2
• /
• pp.129-152
• /
• 2008

This paper considers the solution of the stochastic differential equations (SDEs) with random operator and/or random excitation using the spectral SFEM. The random system parameters (involved in the operator) and the random excitations are modeled as second order stochastic processes defined only by their means and covariance functions. All random fields dealt with in this paper are continuous and do not have known explicit forms dependent on the spatial dimension. This fact makes the usage of the finite element (FE) analysis be difficult. Relying on the spectral properties of the covariance function, the Karhunen-Loeve expansion is used to represent these processes to overcome this difficulty. Then, a spectral approximation for the stochastic response (solution) of the SDE is obtained based on the implementation of the concept of generalized inverse defined by the Neumann expansion. This leads to an explicit expression for the solution process as a multivariate polynomial functional of a set of uncorrelated random variables that enables us to compute the statistical moments of the solution vector. To check the validity of this method, two applications are introduced which are, randomly loaded simply supported reinforced concrete beam and reinforced concrete cantilever beam with random bending rigidity. Finally, a more general application, randomly loaded simply supported reinforced concrete beam with random bending rigidity, is presented to illustrate the method.

• Korea-Australia Rheology Journal
• /
• v.16 no.4
• /
• pp.183-191
• /
• 2004

High Deborah or Weissenberg number problems in viscoelastic flow modeling have been known formidably difficult even in the inertialess limit. There exists almost no result that shows satisfactory accuracy and proper mesh convergence at the same time. However recently, quite a breakthrough seems to have been made in this field of computational rheology. So called matrix-logarithm (here we name it tensor-logarithm) formulation of the viscoelastic constitutive equations originally written in terms of the conformation tensor has been suggested by Fattal and Kupferman (2004) and its finite element implementation has been first presented by Hulsen (2004). Both the works have reported almost unbounded convergence limit in solving two benchmark problems. This new formulation incorporates proper polynomial interpolations of the log­arithm for the variables that exhibit steep exponential dependence near stagnation points, and it also strictly preserves the positive definiteness of the conformation tensor. In this study, we present an alternative pro­cedure for deriving the tensor-logarithmic representation of the differential constitutive equations and pro­vide a numerical example with the Leonov model in 4:1 planar contraction flows. Dramatic improvement of the computational algorithm with stable convergence has been demonstrated and it seems that there exists appropriate mesh convergence even though this conclusion requires further study. It is thought that this new formalism will work only for a few differential constitutive equations proven globally stable. Thus the math­ematical stability criteria perhaps play an important role on the choice and development of the suitable con­stitutive equations. In this respect, the Leonov viscoelastic model is quite feasible and becomes more essential since it has been proven globally stable and it offers the simplest form in the tensor-logarithmic formulation.

• Journal of the Korean Institute of Intelligent Systems
• /
• v.24 no.2
• /
• pp.161-165
• /
• 2014

In this paper, we proposed the fuzzy prototype pattern classifier. In the proposed classifier, each prototype is defined to describe the related sub-space and the weight value is assigned to the prototype. The weight value assigned to the prototype leads to the change of the boundary surface. In order to define the prototypes, we use Fuzzy C-Means Clustering which is the one of fuzzy clustering methods. In order to optimize the weight values assigned to the prototypes, we use the Differential Evolutionary Algorithm. We use Linear Discriminant Analysis to estimate the coefficients of the polynomial which is the structure of the consequent part of a fuzzy rule. Finally, in order to evaluate the classification ability of the proposed pattern classifier, the machine learning data sets are used.

• Proceedings of the KSME Conference
• /
• 2007.05b
• /
• pp.2552-2557
• /
• 2007

Shape optimization of an internal cooling passage with staggered dimples on single surface is performed and performances of surrogates are evaluated in this paper. Optimizations are performed so that turbulent heat transfer can be enhanced compromising with pressure loss due to friction. The three-dimensional governing differential equations have been solved to find the overall Nusselt number and friction factor which are related to the objective functions of this problem. Three design variables were selected among the dimensionless geometric variables. Basic surrogate models such as second order polynomial response surface approximation (RSA), Kriging meta-modeling technique, radial basis neural network (RBNN), and derived press based averaged (PBA) surrogate model are constructed. The optimal points are searched from the above constructed surrogates by sequential quadratic programming (SQP). It is shown that use of multiple surrogates can increase the robustness in prediction of better design with minimum computational cost.

• Coupled systems mechanics
• /
• v.1 no.2
• /
• pp.155-163
• /
• 2012

This paper presents analytical approximate solutions for the initial post-buckling deformation of the pipes in oil and gas wells. The governing differential equation with sinusoidal nonlinearity can be reduced to form a third-order-polynomial nonlinear equation, by coupling of the well-known Maclaurin series expansion and orthogonal Chebyshev polynomials. Analytical approximations to the resulting boundary condition problem are established by combining the Newton's method with the method of harmonic balance. The linearization is performed prior to proceeding with harmonic balancing thus resulting in a set of linear algebraic equations instead of one of non-linear algebraic equations, unlike the classical method of harmonic balance. We are hence able to establish analytical approximate solutions. The approximate formulae for load along axis, and periodic solution are established for derivative of the helix angle at the end of the pipe. Illustrative examples are selected and compared to "reference" solution obtained by the shooting method to substantiate the accuracy and correctness of the approximate analytical approach.

• Structural Engineering and Mechanics
• /
• v.60 no.3
• /
• pp.455-470
• /
• 2016

In this paper, we investigate the free vibration of axially loaded non-uniform Rayleigh cantilever beams. The Rayleigh beams account for the rotary inertia effect which is ignored in Euler-Bernoulli beam theory. Using an inverse problem approach we show, that for certain polynomial variations of the mass per unit length and the flexural stiffness, there exists a fundamental closed form solution to the fourth order governing differential equation for Rayleigh beams. The derived property variation can serve as test functions for numerical methods. For the rotating beam case, the results have been compared with those derived using the Euler-Bernoulli beam theory.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.6 no.4
• /
• pp.390-396
• /
• 1982

The onset of longitudinal vortices in horizontal Blasius flow isothermally heated from below is studied analytically. The assumption that at the onset of thermal instability the thermal disturbances are confined within the thermal boundary layer is employed for the limiting case of large Prandtl number. Polynomial representations for the basic quantities obtained by the integral method of the boundary layer analysis have been used. Then the system of differential equations and boundary conditions for disturbance quantities is reformulated in a convenient form so that the solutions may be constructed as rapidly convergent power series. The critical buoyancy parameter G $r_{x}$ $^{*}$ /R $e^{*1.5}$ falls between 2 and 6, which is about one order of magnitude lower than the existing experimental values. It is also shown that the positions of the onset of instability can be closely predicted by the present theory.y.y.

• Communications of the Korean Mathematical Society
• /
• v.18 no.3
• /
• pp.569-580
• /
• 2003

The Lotka-McKendrick equation which describes the evolution of a single population under the phenomenological conditions is developed from the well-known Malthus’model. In this paper, we introduce the Lotka-McKendrick equation for the description of the dynamics of a population. We apply a discontinuous Galerkin finite element method in age-time domain to approximate the solution of the system. We provide some numerical results. It is experimentally shown that, when the mortality function is bounded, the scheme converges at the rate of $h^2$ in the case of piecewise linear polynomial space. It is also shown that the scheme converges at the rate of $h^{3/2}$ when the mortality function is unbounded.

• Steel and Composite Structures
• /
• v.9 no.5
• /
• pp.473-497
• /
• 2009

For the spatially coupled free vibration analysis of composite box beams resting on elastic foundation under the axial force, the exact solutions are presented by using the power series method based on the homogeneous form of simultaneous ordinary differential equations. The general vibrational theory for the composite box beam with arbitrary lamination is developed by introducing Vlasov°Øs assumption. Next, the equations of motion and force-displacement relationships are derived from the energy principle and explicit expressions for displacement parameters are presented based on power series expansions of displacement components. Finally, the dynamic stiffness matrix is calculated using force-displacement relationships. In addition, the finite element model based on the classical Hermitian interpolation polynomial is presented. To show the performances of the proposed dynamic stiffness matrix of composite box beam, the numerical solutions are presented and compared with the finite element solutions using the Hermitian beam elements and the results from other researchers. Particularly, the effects of the fiber orientation, the axial force, the elastic foundation, and the boundary condition on the vibrational behavior of composite box beam are investigated parametrically. Also the emphasis is given in showing the phenomenon of vibration mode change.

• Journal of Institute of Control, Robotics and Systems
• /
• v.5 no.8
• /
• pp.969-976
• /
• 1999

This paper is concerned with the generation of a balancing trajectory for improving the walking performance. The balancing motion has been determined by solving a second -order differential equation. However, this method caused some difficulties in linearizing and approximating the equation and had restrictions on using various balancing trajectories. The proposed difficulties in linearizing and approximating the equation and had restrictions on using various balancing trajectories. The proposed method i this paper is based on the genetic algorithm for minimizing the motins of balancing joints, whose trajectories are generated by the fifth-order polynomial interpolation after planning leg trajectories. The real walking experiments are made on the biped robot IWR-III, developed by our Automatic Control Laboratory. The system has 8 degrees of freedom and the structure of three pitches in each leg, and one roll and one prismatic joint in the balancing joints. The experimental result shows the validity and applicability of the new proposed algorithm.