• Structural Engineering and Mechanics
• /
• v.28 no.2
• /
• pp.221-238
• /
• 2008

Numerical solution to buckling analysis of beams and columns are obtained by the method of differential quadrature (DQ) and harmonic differential quadrature (HDQ) for various support conditions considering the variation of flexural rigidity. The solution technique is applied to find the buckling load of fully or partially embedded columns such as piles. A simple semi- inverse method of DQ or HDQ is proposed for determining the flexural rigidities at various sections of non-prismatic column ( pile) partially and fully embedded given the buckling load, buckled shape and sub-grade reaction of the soil. The obtained results are compared with the existing solutions available from other numerical methods and analytical results. In addition, this paper also uses a recently developed technique, known as the differential transformation (DT) to determine the critical buckling load of fully or partially supported heavy prismatic piles as well as fully supported non-prismatic piles. In solving the problem, governing differential equation is converted to algebraic equations using differential transformation methods (DT) which must be solved together with applied boundary conditions. The symbolic programming package, Mathematica is ideally suitable to solve such recursive equations by considering fairly large number of terms.

• Structural Engineering and Mechanics
• /
• v.7 no.3
• /
• pp.289-302
• /
• 1999

An approximate method for analyzing the bending problems of irregular-shaped plates is proposed. In this paper irregular-shaped plates are such plates as plate with opening, circular plate, semi-circular plate, elliptic plate, triangular plate, skew plate, rhombic plate, trapezoidal plate or the other polygonal plates which are not uniform rectangular plates. It is shown that these irregular-shaped plates can be considered finally as a kind of rectangular plates with non-uniform thickness. An opening in a plate can be considered as an extremely thin part of the plate, and a non-rectangular plate can be translated into a circumscribed rectangular plate whose additional parts are extremely thin or thick according to the boundary conditions of the original plate. Therefore any irregular-shaped plate can be replaced by the equivalent rectangular plate with non-uniform thickness. For various types of irregular-shaped plates the convergency and accuracy of numerical solution by proposed method are investigated.

• Korean Journal of Mathematics
• /
• v.30 no.2
• /
• pp.403-412
• /
• 2022

In this paper, we introduce the following cubic-quartic functional equation of the form $$f(x+4y)+f(x-4y)=16[f(x+y)+f(x-y)]{\pm}30f(-x)+\frac{5}{2}[f(4y)-64f(y)]$$. Further, we investigate the general solution and the Ulam stability for the above functional equation in non-Archimedean spaces by using the direct method.

• 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.

• Structural Engineering and Mechanics
• /
• v.52 no.5
• /
• pp.997-1017
• /
• 2014

According to the structure characteristics of the non-uniform beam bridge, a practical model for calculating the vibration equation of the non-uniform beam bridge is given and the application scope of the model includes not only the beam bridge structure but also the non-uniform beam with added masses and elastic supports. Based on the Bernoulli-Euler beam theory, extending the application of the modal perturbation method and establishment of a semi-analytical method for solving the vibration equation of the non-uniform beam with added masses and elastic supports based is able to be made. In the modal subspace of the uniform beam with the elastic supports, the variable coefficient differential equation that describes the dynamic behavior of the non-uniform beam is converted to nonlinear algebraic equations. Extending the application of the modal perturbation method is suitable for solving the vibration equation of the simply supported and continuous non-uniform beam with its arbitrary added masses and elastic supports. The examples, that are analyzed, demonstrate the high precision and fast convergence speed of the method. Further study of the timesaving method for the dynamic characteristics of symmetrical beam and the symmetry of mode shape should be developed. Eventually, the effects of elastic supports and added masses on dynamic characteristics of the three-span non-uniform beam bridge are reported.

• Journal of Electrical Engineering and Technology
• /
• v.2 no.3
• /
• pp.312-319
• /
• 2007

The solution of the power flow is one of the most important problems in electrical power systems. These traditional methods such as Gauss-Seidel method and Newton-Raphson (NR) method have had drawbacks up to now such as initial values, abnormal operating solutions and divergences in heavy loads. In order to overcome theses problems, the power flow solution incorporating genetic algorithm (GA) is introduced in this paper. General operator of genetic algorithm, arithmetic crossover, and non-uniform mutation operator of GA are suggested to solve the power flow problem. While abnormal solution cannot be obtained by a NR method, multiple power flow solution can be obtained by a GA method. With a heavy load, both normal solution and abnormal solution can be obtained by a proposed method. In this paper, a floating number representation instead of the binary number representation is introduced for accuracy. Simulation results have been compared with traditional methods.

• Journal of applied mathematics & informatics
• /
• v.24 no.1_2
• /
• pp.191-207
• /
• 2007

Based on a mixed Galerkin approximation, we construct the fully discrete approximations of $U_y$ as well as U to a single-phase quasilinear Stefan problem with a forcing term in non-divergence form. We prove the optimal convergence of approximation to the solution {U, S} and the superconvergence of approximation to $U_y$.

• East Asian mathematical journal
• /
• v.33 no.5
• /
• pp.495-509
• /
• 2017

With the accuracy of the nonlinearity guaranteed, plenty of time and large memory space are needed when we solve the finite element numerical solution of nonlinear partial differential equations. In this paper, we use the Group Element Method (GEM) to deal with the non-linearity of the BBM-Burgers Equation with Conservation form and perform a numerical analysis for two particular initial-boundary value (the Dirichlet boundary conditions and Neumann-Dirichlet boundary conditions) problems with the Finite Element Method (FEM). Some numerical experiments are performed to analyze the error between the exact solution and the FEM solution in MATLAB.

• Proceedings of the KSME Conference
• /
• 2000.04b
• /
• pp.144-147
• /
• 2000

A numerical method for the solution of one-dimensional inverse heat conduction problem is established and its performance is demonstrated with computational results. The present work introduces the maximum entropy method in order to build a robust formulation of the inverse problem. The maximum entropy method finds the solution that maximizes the entropy functional under given temperature measurement. The philosophy of the method is to seek the most likely inverse solution. The maximum entropy method converts the inverse problem to a non-linear constrained optimization problem of which constraint is the statistical consistency between the measured temperature and the estimated temperature. The successive quadratic programming facilitates the maximum entropy estimation. The gradient required fur the optimization procedure is provided by solving the adjoint problem. The characteristic feature of the maximum entropy method is discussed with the illustrated results. The presented results show considerable resolution enhancement and bias reduction in comparison with the conventional methods.

• Communications of the Korean Mathematical Society
• /
• v.28 no.1
• /
• pp.191-207
• /
• 2013

In this paper, we introduce the shrinking projection method for hemi-relatively nonexpansive mappings to find a common solution of variational inequality problems and equilibrium problems in uniformly convex and uniformly smooth Banach spaces and prove some strong convergence theorems to the common solution by using the proposed method.