• Proceedings of the IEEK Conference
• /
• summer
• /
• pp.263-268
• /
• 2004

The inverse kinematics problem is to find a set of joint variable values that will place the end effector of a robot manipulator into a given pose. Pieper has shown that a sufficient condition for a manipulator to have a closed form solution is that three adjacent joint axes intersects, hence the six axes robot with spherical wrist allows closed form solution. But many industrial robots have a non-spherical wrist to provide a stronger wrist configuration so that they can handle heavy payloads. Also, the use of a non-spherical wrist can result in a cheap and simple wrist arrangement than when all three axes intersect at a common point. In these cases, closed form solutions cannot be found. Therefore numerical technique must be used to solve the inverse kinematics equations. This paper proposes a new algorithm that can be used for finding inverse kinematics solution of the six axes robot with non-spherical wrist. Computer simulations are provided to prove the usefulness of our method.

• Journal of the Korean Institute of Telematics and Electronics S
• /
• v.34S no.3
• /
• pp.24-32
• /
• 1997

This paper is on a discrete controller order reduction with the closed-loop stability and performance guaranteed. to achieve this, after finding the solutionsof two lyapunov inequalities and balancing the full order controller system, we find the reudced order controlers using the balanced truncation (BT) and the balanced singular perturbation approximation (BSPA). When the solutions of the two lyapunov inequalities exist, it is shown that the resulting controllers guarantee the closed-loop stability, and .inf.-norm error bounds are derived for the closed-loop performance region for the BT and in low frequency region for the BSPA. Finally, a numerical example is given to illustrate the validity of the proposed method.

• East Asian mathematical journal
• /
• v.40 no.3
• /
• pp.319-328
• /
• 2024

We establish an existence result for a positive solution to the Schrödinger-type singular semipositone problem: $-{\Delta}u\,=\,V(x)u\,=\,{\lambda}{\frac{f(u)}{u^{\alpha}}}$ in Ω, u = 0 on ∂Ω, where Ω is a bounded domain in ℝ^N , N > 2, λ ∈ ℝ is a positive parameter, V ∈ L^∞(Ω), 0 < α < 1, f ∈ C([0, ∞), ℝ) with f(0) < 0. In particular, when ${\frac{f(s)}{s^{\alpha}}}$ is sublinear at infinity, we establish the existence of a positive solutions for λ ≫ 1. The proofs are mainly based on the sub and supersolution method. Further, we extend our existence result to infinite semipositone problems with mixed boundary conditions.

• 한국해양학회지
• /
• v.30 no.6
• /
• pp.543-549
• /
• 1995

The inverse method has been used to estimate a surface velocity field from sequential AVHRR/SST data. In the model, equation system was composed of heat equation and horizontal divergence minimization and the velocity field contained in the advective term of the heat equation, which was linearized in grid system, was estimated. A constraint was the minimization of horizontal divergence with weighting factor and introduced to compensate the null space(Menke, 1984) of the velocity solutions for the heat equation. The experiments were carried out to set up the range of weighting factor and the matrix equation was solved by SVD(Singular Value Decomposion). In the experiment, the scales of horizontal temperature gradient and divergence of synthetic velocity field were approximated to those of real field. The neglected diffusive effect and the horizontal variation of heat flux in the heat equation were regarded as random temperature errors. According to the result of experiments, the minimum of relative error was more desirable than the minimum of misfit as the criteria of setting up the weighting factor and the error of estimated velocity field became small when the weighting factor was order of $10^{-1}$

• Korea-Australia Rheology Journal
• /
• v.15 no.3
• /
• pp.131-150
• /
• 2003

A new numerical algorithm of finite element methods is presented to solve high Deborah number flow problems with geometric singularities. The steady inertialess planar 4 : 1 contraction flow is chosen for its test. As a viscoelastic constitutive equation, we have applied the globally stable (dissipative and Hadamard stable) Leonov model that can also properly accommodate important nonlinear viscoelastic phenomena. The streamline upwinding method with discrete elastic-viscous stress splitting is incorporated. New interpolation functions classified as rational interpolation, an alternative formalism to enhance numerical convergence at high Deborah number, are implemented not for the whole set of finite elements but for a few elements attached to the entrance comer, where stress singularity seems to exist. The rational interpolation scheme contains one arbitrary parameter b that controls the singular behavior of the rational functions, and its value is specified to yield the best stabilization effect. The new interpolation method raises the limit of Deborah number by 2∼5 times. Therefore on average, we can obtain convergent solution up to the Deborah number of 200 for which the comer vortex size reaches 1.6 times of the half width of the upstream reservoir. Examining spatial violation of the positive definiteness of the elastic strain tensor, we conjecture that the stabilization effect results from the peculiar behavior of rational functions identified as steep gradient on one domain boundary and linear slope on the other. Whereas the rational interpolation of both elastic strain and velocity distorts solutions significantly, it is shown that the variation of solutions incurred by rational interpolation only of the elastic strain is almost negligible. It is also verified that the rational interpolation deteriorates speed of convergence with respect to mesh refinement.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.3 no.3
• /
• pp.91-97
• /
• 1979

Boundary layer equations (partial differential equations) can be transformed to ordinary diffential equations with constant coeffieients in terms of similarity transformed to ordinary differential equations with constant coeffieients in terms of similarity transformations in the heat tranfer analysis on the surface of any axiaymmetric boiles. The azimuth functions can not be uniquely determined because of the singular behavior at the stagnation point(X=0.deg.).In spite of the azimuth functions behaving singularly, many of researchers have analyzed the heat transfer problem on a horizontal chlinder or a sphere, supposing the set of solutions( $H_{1}$ & G$_{1}$) of being yieled from the simple differential equation to be unique solution of therazimuth functions. In order to ascertain whether mathematical incompatibility as mentioned above can be admitted in the viewpoint of enginerring or not, condensation heat transfer coefficients on a sphere are computed for all azimuth functions( $H_{1}$ G$_{1}$ & $H_{2}$ G$_{2}$) and comparisons with the experimental result are discussed.

• Journal of Institute of Control, Robotics and Systems
• /
• v.5 no.3
• /
• pp.324-337
• /
• 1999

In this work, a new 3-PSP type spatial 3-degree-of-freedom parallel mechanism is proposed. And a 6 DOF hybrid manipulator which consists of a 3-PPR type planar 3 DOF parallel mechanism and a new 3-PSP type spatial 3-degree-of-freedom parallel mechanism is proposed. Both 3 DOF mechanism modules have closed-form forward position solutions and particularly, 3-PSP spatial module has unique forward position solution. Firstly, the closed-form position analysis and first-order kinematic analysis for the proposed 3-PSP type module are carried out, and the first-order kinematic characteristics are examined via maximum singular value and the isotropic index of the mechanism. It is shown through these analyses that the mechanism has excellent isotrpic property throughout the workspace. Secondly, position and kinematic analysis of the 3-PPR planar module are briefly described. Thirdly, the forward position analysis for the 3-PPR 3-PSP type 6 degree-of-freedom hybrid mechanism consisting of a 3-PPR planar module and a 3-PSP spatial module is performed along with the analysis of the workspace size and first-order kinematic characteristics. The kinematic characteristics of the proposed hybrid manipulator are compared to those of geometrically similar Stewart manipulator.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.17 no.1
• /
• pp.59-65
• /
• 2004

The stress intensity factors have been widely used in numerical studies of crack growth direction. However in many cases, omissive terms of the series expansion are quantitatively significant, so we consider the computation of such terms. For this purpose, we used the finite element method with isoparametric quadratic quarter-point elements. For examples, infinite square plate with a slant crack subjected to a uniaxial load is analyzed. The numerical analysis were performed for the wide range of crack tip element lengths and inclined angles. The numerical results obtained are compared with the theoretical solutions. Also they were accurate and efficient.

• Structural Engineering and Mechanics
• /
• v.54 no.1
• /
• pp.69-85
• /
• 2015

In this paper, a receding contact problem for an elastic layer resting on two quarter planes is considered. The layer is pressed by a stamp and distributed loads. It is assumed that the contact surfaces are frictionless and only compressive traction can be transmitted through the contact surfaces. In addition the effect of body forces are neglected. Firstly, the problem is solved analytically based on theory of elasticity. In this solution, the problem is reduced into a system of singular integral equations in which contact areas and contact stresses are unknowns using boundary conditions and integral transform techniques. This system is solved numerically using Gauss-Jacobi integral formulation. Secondly, two dimensional finite element analysis of the problem is carried out using ANSYS. The dimensionless quantities for the contact areas and the contact pressures are calculated under various distributed load conditions using both solutions. It is concluded that the position and the magnitude of the distributed load have an important role on the contact area and contact pressure distribution between layer and quarter plane contact surface. The analytic results are verified by comparison with finite element results.

• Advances in Computational Design
• /
• v.2 no.4
• /
• pp.313-331
• /
• 2017

In this study, teaching-learning based optimization (TLBO) is improved by incorporating model of multiple teachers, adaptive teaching factor, self-motivated learning, and learning through tutorial. Modified TLBO (MTLBO) is applied for simultaneous topology, shape, and size optimization of space and planar trusses to study its effectiveness. All the benchmark problems are subjected to stress, displacement, and kinematic stability constraints while design variables are discrete and continuous. Analyses of unacceptable and singular topologies are prohibited by seeing element connectivity through Grubler's criterion and the positive definiteness. Performance of MTLBO is compared to TLBO and state-of-the-art algorithms available in literature, such as a genetic algorithm (GA), improved GA, force method and GA, ant colony optimization, adaptive multi-population differential evolution, a firefly algorithm, group search optimization (GSO), improved GSO, and intelligent garbage can decision-making model evolution algorithm. It is observed that MTLBO has performed better or found nearly the same optimum solutions.