• 제목/요약/키워드: Jacobi method

검색결과 84건 처리시간 0.028초

제한 3체문제에서의 작은 질량비에 대한 충돌궤도 (Collision Orbite for Small Mass Ratio in the Restricted Three Body Problem)

  • 조중현;박상영;이병선;최규홍
    • Journal of Astronomy and Space Sciences
    • /
    • 제5권1호
    • /
    • pp.19-30
    • /
    • 1988
  • 지구-달계의 질량비인 $\mu$-0.9878449 와 Jacobi 상수 2.9~3.4인 범위의 조건일 때 j종족 충돌궤도가 존재하는가를 조사하였다. 수치적인 방법으로는 Birckhoff의 정칙화 평면에서 Fehlberg(1968)의 5차 Runge-Kutta가변구간 수치적분방법을 사용하여 궤도의 주기성을 검증하였다. 그 결과 4개의 j종족 충돌궤도와 5개의 특이형궤도를 얻었고, 또한, 이 j 종족 충돌궤도가 Pinotsis와 Zikides(1984)가 제시한 주기와 궤도형태의 상관관계를 보였다.

  • PDF

THE RECURSIVE ALGOFITHM FOR OPTIMAL REGULATOR OF NONSTANCARD SINGULARLY PERTURVED SYSTEMS

  • Mukaidani, Hiroaki;Xu, Hau;Mizukami, Koichi
    • 제어로봇시스템학회:학술대회논문집
    • /
    • 제어로봇시스템학회 1995년도 Proceedings of the Korea Automation Control Conference, 10th (KACC); Seoul, Korea; 23-25 Oct. 1995
    • /
    • pp.10-13
    • /
    • 1995
  • This paper considers the linear-quadratic optimal regulator problem for nonstandard singularly perturbed systems making use of the recursive technique. We first derive a generalized Riccati differential equation by the Hamilton-Jacobi equation. In order to obtain the feedback gain, we must solve the generalized algebraic Riccati equation. Using the recursive technique, we show that the solution of the generalized algebraic Riccati equation converges with the rate of convergence of O(.epsilon.). The existence of a bounded solution of error term can be proved by the implicit function theorem. It is enough to show that the corresponding Jacobian matrix is nonsingular at .epsilon. = 0. As a result, the solution of optimal regulator problem for nonstandard singularly perturbed systems can be obtained with an accuracy of O(.epsilon.$^{k}$ ). The proposed technique represents a significant improvement since the existing method for the standard singularly perturbed systems can not be applied to the nonstandard singularly perturbed systems.

  • PDF

A receding contact problem of a layer resting on a half plane

  • Karabulut, Pembe Merve;Adiyaman, Gokhan;Birinci, Ahmet
    • Structural Engineering and Mechanics
    • /
    • 제64권4호
    • /
    • pp.505-513
    • /
    • 2017
  • In this paper, a receding contact problem for an elastic layer resting on a half plane is considered. The layer is pressed by two rectangular stamps placed symmetrically. 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 is 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 half contact length and contact pressures 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 length and the contact pressures are calculated under various stamp size, stamp position and material properties using both solutions. The analytic results are verified by comparison with finite element results.

Dynamic Robust Path-Following Using A Temporary Path Generator for Mobile Robots with Nonholonomic Constraints

  • Lee, Seunghee;Jongguk Yim;Park, Jong-Hyeon
    • 제어로봇시스템학회:학술대회논문집
    • /
    • 제어로봇시스템학회 2000년도 제15차 학술회의논문집
    • /
    • pp.515-515
    • /
    • 2000
  • The performance of dynamic path following of a wheeled mobile robot with nonholonomic constraints has some drawbacks such as the influence of the initial state. The drawbacks can be overcome by the temporary path generator and modified output. But with the previous input-output linearization method using them, it is difficult to tune the gains, and if there are some modeling errors, the low gain can make the system unstable. And if a high gain is used to overcome the model uncertainties, the control inputs are apt to be large so the system can be unstable. In this paper. an H$_{\infty}$ controller is designed to guarantee robustness to model parameter uncertainties and to consider the magnitude of control inputs. And the solution to Hamilton Jacobi (HJ) inequality, which is essential to H$_{\infty}$ control design, is obtained by nonlinear matrix inequality (NLMI).

  • PDF

Parallel Robust $H_{\infty}$ Control for Weakly Coupled Bilinear Systems with Parameter Uncertainties Using Successive Galerkin Approximation

  • Kim, Young-Joong;Lim, Myo-Taeg
    • International Journal of Control, Automation, and Systems
    • /
    • 제4권6호
    • /
    • pp.689-696
    • /
    • 2006
  • This paper presents a new algorithm for the closed-loop $H_{\infty}$ composite control of weakly coupled bilinear systems with time-varying parameter uncertainties and exogenous disturbance using the successive Galerkin approximation(SGA). By using weak coupling theory, the robust $H_{\infty}$ control can be obtained from two reduced-order robust $H_{\infty}$ control problems in parallel. The $H_{\infty}$ control theory guarantees robust closed-loop performance but the resulting problem is difficult to solve for uncertain bilinear systems. In order to overcome the difficulties inherent in the $H_{\infty}$ control problem, two $H_{\infty}$ control laws are constructed in terms of the approximated solution to two independent Hamilton-Jacobi-Isaac equations using the SGA method. One of the purposes of this paper is to design a closed-loop parallel robust $H_{\infty}$ control law for the weakly coupled bilinear systems with parameter uncertainties using the SGA method. The other is to reduce the computational complexity when the SGA method is applied to the high order systems.

유한요소법을 이용한 코오드-고무 복합판의 동적특성에 관한 연구 (A study on the dynamic characteristics of the cord-rubber laminates rectangular plate by finite element method)

  • 김두만;김항욱
    • 오토저널
    • /
    • 제8권2호
    • /
    • pp.51-64
    • /
    • 1986
  • There has been considerable interest over the last twenty years in the subject of the elastic properties of the cord-rubber laminate. This has been due to the rather intensive study of the composites materials characteristics brought about by the increased use of rigid composites materials characteristics brought about by the increased use of rigid composites in many structural applications. The object of this study is to obtain the natural frequencies and modes of the simply supported cord-rubber laminate plates prior to the study on the analysis of the dynamic properties of the pneumatic tire. To obtain these natural frequencies and modes, the 12 degrees of freedom orthotropic rectangular plate finite elements are developed. By using classical lamination theory, the stress-strain relations are represented. The governing equation for the finite element is derived by energy method. To find the natural frequencies and modes, he eigenvalues and corresponding eigenvectors are computed by the well known Jacobi power method. In order to verify the capability of this present finite element, the results of the specially orthotropic plate and the angle-ply laminate plate are compared with the analytical solution. The analytical and numberical results are in good agreement. The following problems of the simply supported plate are analyzed by the present finite element. a) the natural frequencies and mode shapes of the cord-rubber laminate plate for various aspect ratio. b) The natural frequencies and mode shapes of the orthotropic plate with the rectangular hole in its center.

  • PDF

Investigation of nonlinear free vibration of FG-CNTRC cylindrical panels resting on elastic foundation

  • J.R. Cho
    • Structural Engineering and Mechanics
    • /
    • 제88권5호
    • /
    • pp.439-449
    • /
    • 2023
  • Non-linear vibration characteristics of functionally graded CNT-reinforced composite (FG-CNTRC) cylindrical shell panel on elastic foundation have not been sufficiently examined. In this situation, this study aims at the profound numerical investigation of the non-linear vibration response of FG-CNTRC cylindrical panels on Winkler-Pasternak foundation by introducing an accurate and effective 2-D meshfree-based non-linear numerical method. The large-amplitude free vibration problem is formulated according to the first-order shear deformation theory (FSDT) with the von Karman non-linearity, and it is approximated by Laplace interpolation functions in 2-D natural element method (NEM) and a non-linear partial derivative operator HNL. The complex and painstaking numerical derivation on the curved surface and the crucial shear locking are overcome by adopting the geometry transformation and the MITC3+ shell elements. The derived nonlinear modal equations are iteratively solved by introducing a three-step iterative solving technique which is combined with Lanczos transformation and Jacobi iteration. The developed non-linear numerical method is estimated through the benchmark test, and the effects of foundation stiffness, CNT volume fraction and functionally graded pattern, panel dimensions and boundary condition on the non-linear vibration of FG-CNTRC cylindrical panels on elastic foundation are parametrically investigated.

다중회귀모형으로 추정된 모수에 의한 최적단위유량도의 유도에 관한 연구 (A Study on the Derivation of the Unit Hydrograph using Multiple Regression Model)

  • 이종남;김채원;황창현
    • 물과 미래
    • /
    • 제25권1호
    • /
    • pp.93-100
    • /
    • 1992
  • Abstract A study on the Derivation of the Unit Hydrograph using Multiple Regression Moe이. The purpose of this study is to deriver an optimal unit hydrograph suing the multiple regression model, particularly when only small amount of data is available. The presence of multicollinearity among the input data can cause serious oscillations in the derivation of the unit hydrograph. In this case, the oscillations in the unit hydrograph ordinate are eliminated by combining the data. The data used in this study are based upon the collection and arrangement of rainfall-runoff data(1977-1989) at the Soyang-river Dam site. When the matrix X is the rainfall series, the condition number and the reciprocal of the minimum eigenvalue of XTX are calculated by the Jacobi an method, and are compared with the oscillation in the unit hydrograph. The optimal unit hydrograph is derived by combining the numerous rainfall-runoff data. The conclusions are as follows; 1)The oscillations in the derived unit hydrograph are reduced by combining the data from each flood event. 2) The reciprocals of the minimum eigen\value of XTX, 1/k and the condition number CN are increased when the oscillations are active in the derived unit hydrograph. 3)The parameter estimates are validated by extending the model to the Soyang river Dam site with elimination of the autocorrelation in the disturbances. Finally, this paper illustrates the application of the multiple regression model to drive an optimal unit hydrograph dealing with the multicollinearity and the autocorrelation which cause some problems.

  • PDF

AN ASSESSMENT OF PARALLEL PRECONDITIONERS FOR THE INTERIOR SPARSE GENERALIZED EIGENVALUE PROBLEMS BY CG-TYPE METHODS ON AN IBM REGATTA MACHINE

  • Ma, Sang-Back;Jang, Ho-Jong
    • Journal of applied mathematics & informatics
    • /
    • 제25권1_2호
    • /
    • pp.435-443
    • /
    • 2007
  • Computing the interior spectrum of large sparse generalized eigenvalue problems $Ax\;=\;{\lambda}Bx$, where A and b are large sparse and SPD(Symmetric Positive Definite), is often required in areas such as structural mechanics and quantum chemistry, to name a few. Recently, CG-type methods have been found useful and hence, very amenable to parallel computation for very large problems. Also, as in the case of linear systems proper choice of preconditioning is known to accelerate the rate of convergence. After the smallest eigenpair is found we use the orthogonal deflation technique to find the next m-1 eigenvalues, which is also suitable for parallelization. This offers advantages over Jacobi-Davidson methods with partial shifts, which requires re-computation of preconditioner matrx with new shifts. We consider as preconditioners Incomplete LU(ILU)(0) in two variants, ever-relaxation(SOR), and Point-symmetric SOR(SSOR). We set m to be 5. We conducted our experiments on matrices from discretizations of partial differential equations by finite difference method. The generated matrices has dimensions up to 4 million and total number of processors are 32. MPI(Message Passing Interface) library was used for interprocessor communications. Our results show that in general the Multi-Color ILU(0) gives the best performance.

정적 연소실내 난류 예혼합화염 전파의 시뮬레이션 (Simulation of Turbulent Premixed Flame Propagation in a Closed Vessel)

  • 권세진
    • 대한기계학회논문집
    • /
    • 제19권6호
    • /
    • pp.1510-1517
    • /
    • 1995
  • A theoretical method is described to simulate the propagation of turbulent premixed flames in a closed vessel. The objective is to develop and test an efficient technique to predict the propagation speed of flame as well as the geometric structure of the flame surfaces. Flame is advected by the statistically generated turbulent flow field and propagates as a wave by solving twodimensional Hamilton-Jacobi equation. In the simulation of the unburned gas flow field, following turbulence properties were satisfied: mean velocity field, turbulence intensities, spatial and temporal correlations of velocity fluctuations. It is assumed that these properties are not affected by the expansion of the burned gas region. Predictions were compared with existing experimental data for flames propagating in a closed vessel charged with hydrogen/air mixture with various turbulence intensities and Reynolds numbers. Comparisons were made in flame radius growth rate, rms flame radius fluctuations, and average perimeter and fractal dimensions of the flame boundaries. Two dimensional time dependent simulation resulted in correct trends of the measured flame data. The reasonable behavior and high efficiency proves the usefulness of this method in difficult problems of flame propagation such as in internal combustion engines.