• Tribology and Lubricants
• /
• v.38 no.3
• /
• pp.84-92
• /
• 2022

In Part I, developed was a method to obtain the stress field due to an edge dislocation that locates in an elastic half plane beneath the contact edge of an elastically similar square wedge. Essential result was the corrective functions which incorporate a traction free condition of the free surfaces. In the sequel to Part I, features of the corrective functions, F_kij,(k = x, y;i,j = x,y) are investigated in this Part II at first. It is found that F_xxx(ŷ) = F_xyx(ŷ) where ŷ = y/η and η being the location of an edge dislocation on the y axis. When compared with the corrective functions derived for the case of an edge dislocation at x = ξ, analogy is found when the indices of y and x are exchanged with each other as can be readily expected. The corrective functions are curve fitted by using the scatter data generated using a numerical technique. The algebraic form for the curve fitting is designed as F_kij(ŷ) = $\frac{1}{\hat{y}^{1-{\lambda}}I+yp}$$\sum_{q=0}^{m}{\left}$$\left[A_q\left(\frac{\hat{y}}{1+\hat{y}} \right)^q \right]$ where λ_I=0.5445, the eigenvalue of the adhesive complete contact problem introduced in Part I. To investigate the exponent of F_kij, i.e.(1 - λ_I) and p, Log|F_kij|(ŷ)-Log|(ŷ)| is plotted and investigated. All the coefficients and powers in the algebraic form of the corrective functions are obtained using Mathematica. Method of analyzing a surface perpendicular crack emanated from the complete contact edge is explained as an application of the curve-fitted corrective functions.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.4D
• /
• pp.555-561
• /
• 2009

In lieu of section profile data, a fitting of the bored tunnel shape is more significant confirmation for maintenance of a tunnel. Before the permit on the completion of a tunnel, deformation of the completed tunnel with respect to the design model are considered. And deformation can be produced at continuously along the entire of the tunnel section. This study firstly includes an analysis of algebraic approach and test it with an observed field data. And then a number of methods, line search method, genetic algorithm, and pattern search methods, are compared with the 3D tunnel shape fitting. Algebraic methods can solve a simple circular cylinder type as like a railway tunnel. However, a more complex model (compound circular curve and non circular) as like a highway tunnel has to be solved with soft computing tools in the cause of conditional constraints. The genetic algorithm and pattern search methods are computationally more intensive, but they are more flexible at a complex condition. The line search method is fastest, but it needs a narrow bounds of the initial values.

• Journal for History of Mathematics
• /
• v.27 no.4
• /
• pp.271-283
• /
• 2014

This study takes a look at Polya's analysis on Archimedes' "The Method" from a math-historical perspective. We, based on the elaboration of Polya's analysis, investigate the analytic geometric characteristics of Archimedes' "The Method" and discuss the way of using the characteristics in education of school calculus. So this study brings up the educational need of approach of teaching the definite integral by clearly disclosing the transition from length, area, volume etc into the length as an area function under a curve. And this study suggests the approach of teaching both merit and deficiency of the indivisibles method, and the educational necessity of making students realizing that the strength of analytic geometry lies in overcoming deficiency of the indivisibles method by dealing with the relation of variation and rate of change by means of algebraic expression and graph.

• Journal of the Korean Society for Precision Engineering
• /
• v.17 no.12
• /
• pp.121-128
• /
• 2000

The scope of this paper is focused to the systems which have the time period and they should be necessarily studied in the sense of stability and design method of controller to stabilize the orignal unstable systems. In general, the time periodic systems or the systems having same motions during certain time interval are easily found in rotating motion device, i.e., satellite or helicopter and widely used in factory automation systems. The characteristics of the selected dynamic systems are analyzed with the new stability concept and stabilization control method based on Lyapunov direct method. The new method from Lyapunov stability criteria which satisfies the energy convergence is studied with linear algebraic method. And the numerical procedures are developed with computational programming method to apply to the practical linear periodic systems. The results from this paper demonstrate the usefulness in analysis of the asymptotic stability and stabilization of the unstable linear periodic system by using the developed simulation procedures.

• Journal of the Korean Mathematical Society
• /
• v.52 no.5
• /
• pp.1069-1096
• /
• 2015

In this article we introduce a numerical method, named Gegenbauer wavelets method, which is derived from conventional Gegenbauer polynomials, for solving fractional initial and boundary value problems. The operational matrices are derived and utilized to reduce the linear fractional differential equation to a system of algebraic equations. We perform the convergence analysis for the Gegenbauer wavelets method. We also combine Gegenbauer wavelets operational matrix method with quasilinearization technique for solving fractional nonlinear differential equation. Quasilinearization technique is used to discretize the nonlinear fractional ordinary differential equation and then the Gegenbauer wavelet method is applied to discretized fractional ordinary differential equations. In each iteration of quasilinearization technique, solution is updated by the Gegenbauer wavelet method. Numerical examples are provided to illustrate the efficiency and accuracy of the methods.

• East Asian mathematical journal
• /
• v.33 no.1
• /
• pp.53-65
• /
• 2017

Numerical solutions of the Rosenau-Burgers equation based on the cubic B-spline finite element method are introduced. The backward Euler method is used for discretization in time, and the obtained nonlinear algebraic system is changed to a linear system by the Newton's method. We show that those methods are unconditionally stable. Two test problems are studied to demonstrate the accuracy of the proposed method. The computational results indicate that numerical solutions are in good agreement with exact solutions.

• The Journal of Natural Sciences
• /
• v.8 no.2
• /
• pp.9-11
• /
• 1996

We apply a full mutigrid scheme to computing eigenvalues of the Laplace eigenvalue problem with Dirichlet boundary condition. We use finite difference method to get an algebraic equation and apply inverse power method to estimating the smallest eigenvalue. Our result shows that combined method of inverse power method and full multigrid scheme is very effective in calculating eigenvalue of the eigenvalue problem.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2007.04a
• /
• pp.495-500
• /
• 2007

A Meshfree is a method used to establish algebraic equations of system for the whole problem domain without the use of a predefined mesh for the domain discretization. A point interpolation method is based on combining radial and polynomial basis functions. Involvement of radial basis functions overcomes possible singularity. Furthermore, the interpolation function passes through all scattered points in an influence domain and thus shape functions are of delta function property. This makes the implementation of essential boundary conditions much easier than the meshfree methods based on the moving least-squares approximation. This study aims to investigate a stress analysis of structural element between a meshfree method and the finite element method. Examples on cantilever type plate and stress concentration problems show that the accuracy and convergence rate of the meshfree methods are high.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.21 no.8
• /
• pp.1311-1321
• /
• 1997

In multibody dynamics, differential and algebraic equations which can satisfy both equation of motion and kinematic constraint equation should be solved. To solve these equations, coordinate partitioning method and constraint stabilization method are commonly used. In the coordinate partitioning method, the coordinates are divided into independent and dependent and coordinates. The most typical coordinate partitioning method are LU decomposition, QR decomposition, and SVD (singular value decomposition). The objective of this research is to find an efficient coordinate partitioning method in the dynamic analysis of flexible multibody systems. Comparing two coordinate partitioning methods, i.e. LU and QR decomposition in the flexible multibody systems, a new hybrid coordinate partitioning method is suggested for the flexible multibody analysis.

• International Journal of Automotive Technology
• /
• v.2 no.3
• /
• pp.117-122
• /
• 2001

For a dynamic analysis of a constrained multibody system, it is necessary to have a routine for satisfying kinematic constraints. LU decomposition scheme, which is used to divide coordinates into dependent and independent coordinates, is efficient but has great difficulty near the singular configuration. Other method such as the projection matrix, which is more stable near a singular configuration, takes longer simulation time due to the large amount of calculation for decomposition. In this paper, the row space and the null space of the Jacobian matrix are proposed by using the pseudo-inverse method and the projection matrix. The equations of the motion of a system are replaced with independent acceleration components using the null space of the Jacobian matrix. Also a new hybrid method is proposed, combining the LU decomposition and the projection matrix. The proposed hybrid method has following advantages. (1) The simulation efficiency is preserved by the LU method during the simulation. (2) The accuracy of the solution is also achieved by the projection method near the singular configuration.