• Journal of applied mathematics & informatics
• /
• v.18 no.1_2
• /
• pp.171-182
• /
• 2005

An analysis of the non-homogeneous term involved in the free surface condition for second order wave diffraction on a pair of cylinders is presented. In the computations of the nonlinear loads on offshore structures, the most challenging task is the computation of the free surface integral. The main contribution to this integrand is due to the non-homogeneous term present in the free surface condition for second order scattered potential. In this paper, the free surface condition for the second order scattered potential is derived. Under the assumption of large spacing between the two cylinders, waves scattered by one cylinder may be replaced in the vicinity of the other cylinder by equivalent plane waves together with non-planner correction terms. Then solving a complex matrix equation, the first order scattered potential is derived and since the free surface term for second order scattered potential can be expressed in terms of the first order potentials, the free surface term can be obtained using the knowledge of first order potentials only.

• Structural Engineering and Mechanics
• /
• v.43 no.2
• /
• pp.163-177
• /
• 2012

This paper presents a LMI(linear matrix inequality)-based fuzzy approach of modeling and active vibration control of geometrically nonlinear flexible plates with piezoelectric materials as actuators and sensors. The large-amplitude vibration characteristics and dynamic partial differential equation of a piezoelectric flexible rectangular thin plate structure are obtained by using generalized Fourier series and numerical integral. Takagi-Sugeno (T-S) fuzzy model is employed to approximate the nonlinear structural system, which combines the fuzzy inference rule with the local linear state space model. A robust fuzzy dynamic output feedback control law based on the T-S fuzzy model is designed by the parallel distributed compensation (PDC) technique, and stability analysis and disturbance rejection problems are guaranteed by LMI method. The simulation result shows that the fuzzy dynamic output feedback controller based on a two-rule T-S fuzzy model performs well, and the vibration of plate structure with geometrical nonlinearity is suppressed, which is less complex in computation and can be practically implemented.

• The Transactions of the Korean Institute of Electrical Engineers
• /
• v.39 no.1
• /
• pp.22-28
• /
• 1990

This paper is a study on the reduction of computation time in case of nonlinear magnetic field analysis by finite element method and Newton-Raphson method. For the purpose, the nonlinear convergence equation is computed by the conjugate gradient method which is known to be applicable to symmetric positive definite matrix equations only. As the results, we can not prove mathematically that the system Jacobian is positive definite, but when we applied this method, the diverging case did not occur. And the computation time is reduced by 25-55% and 15-45% in comparison with the case of direct and successive over-relaxation method, respectively. Therefore, we proved the utility of conjugate gradient method.

• Journal of the Korean Society for Precision Engineering
• /
• v.12 no.6
• /
• pp.145-151
• /
• 1995

The first-order Taylor series expansion can be evaluated analytically from the formulated symbolic nonlinear dynamic equations. A closed-form linear dynamic euation is derived about a nominal trajectory. The state space representation of the linearized dynamics can be derived easily from the closed-form linear dynamic equations. But manual symbolic expansion of dynamic equations and linearization is tedious, time-consuming and error-prone. So it is desirable to manipulate the procedures using a computer. In this paper, the analytic linearization is performed using the symbolic language MATHEMATICA. Two examples are given to illustrate the approach anbd to compare nonlinear model with linear model.

• Journal of Institute of Control, Robotics and Systems
• /
• v.11 no.11
• /
• pp.895-900
• /
• 2005

We derive a linear matrix inequality(LMI) condition guaranteeing that any invariant zeros of a triple (A, B, C) lie in the open left half plane of the complex plane, i.e. $C(sI-A)^{-1}B$ is minimum phase. The LMI condition is equivalent to a certain constrained Lyapunov matrix equation which can be found in many results relating to stability analysis or control design. We show that the LMI condition can be used to simplify various control engineering problems such as a dynamic output feedback control problem, a variable structure static output feedback control problem, and a nonlinear system observer design problem. Finally, we give some numerical examples.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.19 no.2
• /
• pp.608-616
• /
• 2018

If a general nonlinear system is linearized by the successive multiplication of the 1st and 2nd order systems, then there are four types of poles in this linearized system: the pole of the 1st order system and the equal poles, two distinct real poles, and complex conjugate pair of poles of the 2nd order system. Linear Quadratic (LQ) control is a method of designing a control law that minimizes the quadratic performance index. It has the advantage of ensuring the stability of the system and the pole placement of the root of the system by weighted matrix adjustment. LQ control by the weighted matrix can move the position of the pole of the system arbitrarily, but it is difficult to set the weighting matrix by the trial and error method. This problem can be solved using the characteristic equations of the Hamiltonian system, and if the control weighting matrix is a symmetric matrix of constants, it is possible to move several poles of the system to the desired closed loop poles by applying the control law repeatedly. The paper presents a method of calculating the state weighting matrix and the control law for moving the equal poles with Jordan blocks to two real poles using the characteristic equation of the Hamiltonian system. We express this characteristic equation with a state weighting matrix by means of a trigonometric function, and we derive the relation function (${\rho},\;{\theta}$) between the equal poles and the state weighting matrix under the condition that the two real poles are the roots of the characteristic equation. Then, we obtain the moving-range of the two real poles under the condition that the state weighting matrix becomes a positive semi-finite matrix. We calculate the state weighting matrix and the control law by substituting the two real roots selected in the moving-range into the relational function. As an example, we apply the proposed method to a simple example 3rd order system.

• Journal of the Korean Society for Precision Engineering
• /
• v.29 no.9
• /
• pp.986-994
• /
• 2012

This paper is concerned with the application of the vision control algorithm with weighting matrix in robot point placement task. The proposed vision control algorithm involves four models, which are the robot kinematic model, vision system model, the parameter estimation scheme and robot joint angle estimation scheme. This proposed algorithm is to make the robot move actively, even if relative position between camera and robot, and camera's focal length are unknown. The parameter estimation scheme and joint angle estimation scheme in this proposed algorithm have form of nonlinear equation. In particular, the joint angle estimation model includes several restrictive conditions. For this study, the weighting matrix which gave various weighting near the target was applied to the parameter estimation scheme. Then, this study is to investigate how this change of the weighting matrix will affect the presented vision control algorithm. Finally, the effect of the weighting matrix of robot vision control algorithm is demonstrated experimentally by performing the robot point placement.

• Structural Engineering and Mechanics
• /
• v.54 no.4
• /
• pp.809-827
• /
• 2015

In this paper, seismic energy response of inelastic steel structures under earthquake excitations is investigated. For this purpose, a numerical procedure based on nonlinear dynamic analysis is developed by considering material, geometric and connection nonlinearities. Material nonlinearity is modeled by the inversion of Ramberg-Osgood equation. Nonlinearity caused by the interaction between the axial force and bending moment is also defined considering stability functions, while the geometric nonlinearity caused by axial forces is described using geometric stiffness matrix. Cyclic behaviour of steel connections is taken into account by employing independent hardening model. Dynamic equation of motion is solved by Newmark's constant acceleration method in the time history domain. Energy response analysis of space frames is performed by using this proposed numerical method. Finally, for the first time, the distribution of the different energy types versus time at the duration of the earthquake ground motion is obtained where in addition error analysis for the numerical solutions is carried out and plotted depending on the relative error calculated as a function of energy balance versus time.

• Honam Mathematical Journal
• /
• v.43 no.3
• /
• pp.373-383
• /
• 2021

In this study, a numerical method deals with the Chebyshev wavelet collocation and Adomian decomposition methods are proposed for solving Korteweg-de Vries equation. Integration of the Chebyshev wavelets operational matrices is derived. This problem is reduced to a system of non-linear algebraic equations by using their operational matrix. Thus, it becomes easier to solve KdV problem. The error estimation for the Chebyshev wavelet collocation method and ADM is investigated. The proposed method's validity and accuracy are demonstrated by numerical results. When the exact and approximate solutions are compared, for non-linear or linear partial differential equations, the Chebyshev wavelet collocation method is shown to be acceptable, efficient and accurate.

• Proceedings of the KSR Conference
• /
• 2000.05a
• /
• pp.78-85
• /
• 2000

Computer algorithms for the loadflow of the DC traction power supply system are examined. Algorithms to solve the nodal equation are reviewed and the two iterative methods to solve the nonlinear nature of the loadflow are analyzed and tested, which are so called conductance matrix method and current vector iterative mettled. The result of the analysis tells that the current vector iterative method makes faster convergency and needs less computing time, and it is verified by the test running of the programs based on each of the iterative methods.