• Journal of the Korean Society for Precision Engineering
• /
• v.12 no.5
• /
• pp.57-68
• /
• 1995

A general procedure for the numerical solution of coupled, nonlinear, differential two-point boundary-value problems, solutions of which are crucial to the controller design, has been developed and demonstrated. A fixed-end-points, free-terminal-time, optimal-control problem, which is derived from Pontryagin's Maximum Principle, is solved by an extension of Davidenko's method, a differential form of Newton's method, for algebraic root finding. By a discretization process like finite differences, the differential equations are converted to a nonlinear algebraic system. Davidenko's method reconverts this into a pseudo-time-dependent set of implicitly coupled ODEs suitable for solution by modern, high-performance solvers. Another important advantage of Davidenko's method related to the time-optimal problem is that the terminal time can be computed by treating this unkown as an additional variable and sup- plying the Hamiltonian at the terminal time as an additional equation. Davidenko's method uas used to produce optimal trajectories of a single-degree-of-freedom problem. This numerical method provides switching times for open-loop control, minimized terminal time and optimal input torque sequences. This numerical technique could easily be adapted to the multi-point boundary-value problems.

• Journal of the Korean Institute of Illuminating and Electrical Installation Engineers
• /
• v.15 no.6
• /
• pp.73-81
• /
• 2001

This study uses the distributed parameter systems resented by the spatial discretization technique. In this paper, the distributed parameter systems are converted into lumped parameter systems, End fast Walsh transform and the Picard's iteration method are allied to analysis the state of the systems. This thesis presents a new algorithm which usefully exercises the optimal contro1 in the distributed parameter systems. In exercising the optimal control of the distributed parameter systems, the excellent consequences are found without using the existing decentralized contro1 or hierarchical control method. This study can be applied to the linear time-varying systems and the non-linear systems. Farther researches are required to solve the problems of convergence in case of the numerous applicable intervals. The simulation proves the effectiveness of the proposed algorithm.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.11 no.4
• /
• pp.91-102
• /
• 1991

Traditionally, a dynamic programming(DP) technique has been used to the multireservoir control system. The algorithm has inherent problem to increase computational requirements exponentially due to discretization of variables and expanding the dimension of the system. To solve this problem, this paper describes transforming the optimal control system into a quadratic programming(QP), using Laguerre polynomials(LP) and its properties. The objective function of the proposed QP is independent of time variable. The solution of the QP is obtained by nonlinear programming(NLP) using augmented Lagrangian multiplier method. The numerical experiment shows that the water level of reservoirs is higher than Lee's and the evaluated benefit value is about the same as other researcher's.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.8 no.3
• /
• pp.1-10
• /
• 1988

A higher order plate bending finite element using cubic in-plane displacement profiles is proposed for geometrically nonlinear analysis of thin and thick plates. The higher order plate bending element has been derived from the three dimensional plate-like continuum by discretization of the equations of motion by Galerkin weighted residual method, together with enforcing higher order plate assumptions. Total Lagrangian formulation has been used for geometrically nonlinear analysis of plates and consistent linearization by Newton-Raphson method has been performed to solve the nonlinear equations. The element characteristics have been computed by, selective reduced integration technique using Gauss quadrature to avoid shear locking phenomenon in case of extremely thin plates. Several numerical examples were solved with FEAP macro program to demonstrate versatility and accuracy of the present higher order plate bending element.

• The Journal of Information Systems
• /
• v.17 no.3
• /
• pp.25-37
• /
• 2008

In a fault tree analysis, an uncertainty importance measure is often used to assess how much uncertainty of the top event probability (Q) is attributable to the uncertainty of a basic event probability ($q_i$), and thus, to identify those basic events whose uncertainties need to be reduced to effectively reduce the uncertainty of Q. For evaluating the measures suggested by many authors which assess a percentage change in the variance V of Q with respect to unit percentage change in the variance $v_i$ of $q_i$, V and ${\partial}V/{\partial}v_i$ need to be estimated analytically or by Monte Carlo simulation. However, it is very complicated to analytically compute V and ${\partial}V/{\partial}v_i$ for large-sized fault trees, and difficult to estimate them in a robust manner by Monte Carlo simulation. In this paper, we propose a method for evaluating the measure using discretization technique and Monte Carlo simulation. The proposed method provides a stable uncertainty importance of each basic event.

• Advances in aircraft and spacecraft science
• /
• v.5 no.3
• /
• pp.349-362
• /
• 2018

A non-linear numerical simulation technique for predicting the unsteady performances of an airbreathing engine is developed. The study focuses on the simulation of integrated propulsion systems, where a closer coupling is needed between the airframe and the engine dynamics. In fact, the solution of the fully unsteady flow governing equations, rather than a lumped volume gas dynamics discretization, is essential for modeling the coupling between aero-servoelastic modes and engine dynamics in highly integrated propulsion systems. This consideration holds for any propulsion system when a full separation between the fluid dynamic time-scale and engine transient cannot be appreciated, as in the case of flow instabilities (e.g., rotating stall, surge, inlet unstart), or in case of sudden external perturbations (e.g., gas ingestion). Simulations of the coupling between external and internal flow are performed. The flow around the nacelle and inside the engine ducts (i.e., air intakes, nozzles) is solved by CFD computations, whereas the flow evolution through compressor and turbine bladings is simulated by actuator disks. Shaft work balance and rotor dynamics are deduced from the estimated torque on each turbine/compressor blade row.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2000.06a
• /
• pp.309-313
• /
• 2000

A structural-acoustic coupling problem involving fluid in a cavity divided with flexible walls and porous materials is investigated in this paper. In many practical problems, the use of finite elements to discretize the fluid region leads to large stiffness and mass matrices. But, since the acoustic boundary element discretization requires to put elements only on the surface of structure, the size of matrices is reduced considerably. Here, we developed a numerical analysis program for the structural-acoustic coupling problems of the multi-region cavity, using boundary elements for the fluid regions and finite elements for the structure. By considering sound transmission through layered systems placed in a cavity, the accuracy of the coupled acoustical-structural finite element model has been verified by comparing its transmission loss predictions with analytical sloutions. Example problems are included to investigate the characteristics of the multi-region structural-acoustic coupling system with porous material.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.12 no.4_1
• /
• pp.1-8
• /
• 1992

This paper presents a p-version finite element approach for modeling the stress distribution around a circular hole in a finite strip subjected to membrane and flexural behaviors. Also, same problem with a crack emanating from a perforated tension strip was solved by virtual crack extension method. The p-version of the finite element method based on integrals of Legendre polynomials is shown to perform very well for modeling geometries with very steep stress gradients in the vicinity of a circular cutout. Here, the transfinite mapping technique for circular boundaries was used to avoid the discretization errors. The numerical results from the proposed scheme have a good comparison with those by Nisida, Howland, Newman etc. and the conventional finite element approach.

• The KIPS Transactions:PartD
• /
• v.10D no.3
• /
• pp.501-512
• /
• 2003

This paper proposes an indexing technique for fast retrieval of similar image subsequences using the multi-dimensional time warping distance. The time warping distance is a more suitable similarity measure than Lp distance in many applications where sequences may be of different lengths and/or different sampling rates. Our indexing scheme employs a disk-based suffix tree as an index structure and uses a lower-bound distance function to filter out dissimilar subsequences without false dismissals. It applies the normaliration for an easier control of relative weighting of feature dimensions and the discretization to compress the index tree. Experiments on medical and synthetic image sequences verify that the proposed method significantly outperforms the naive method and scales well in a large volume of image sequence databases.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.8
• /
• pp.8-17
• /
• 2005

It is well known that preconditioning methods are efficient for convergence acceleration in the compressible low Mach number flows. In this study, the original Euler equations and three differently nondimensionalized preconditioning methods are implemented in two dimensional inviscid bump flows using the 3rd order MUSCL and DADI schemes as numerical flux discretization and time integration, respectively. The multigrid and local time stepping methods are also used to accelerate the convergence. The test case indicates that a properly modified local preconditioning technique involving concepts of a global preconditioning allows Mach number independent convergence. Besides, an asymptotic analysis for properties of preconditioning methods is added.