• Computational Structural Engineering
• /
• v.11 no.1
• /
• pp.205-216
• /
• 1998

A solution method is presented to solve the eigenvalue problem arising in the dynamic analysis of nonclassicary damped structural systems with multiple eigenvalues. The proposed method is obtained by applying the modified Newton-Raphson technique and the orthonormal condition of the eigenvectors to the linear eigenproblem through matrix augmentation of the quadratic eigenvalue problem. In the iteration methods such as the inverse iteration method and the subspace iteration method, singularity may be occurred during the factorizing process when the shift value is close to an eigenvalue of the system. However, even though the shift value is an eigenvalue of the system, the proposed method provides nonsingularity, and that is analytically proved. Since the modified Newton-Raphson technique is adopted to the proposed method, initial values are need. Because the Lanczos method effectively produces better initial values than other methods, the results of the Lanczos method are taken as the initial values of the proposed method. Two numerical examples are presented to demonstrate the effectiveness of the proposed method and the results are compared with those of the well-known subspace iteration method and the Lanczos method.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 2004.10a
• /
• pp.347-350
• /
• 2004

In this paper, the acceleration is studied for the rigid-plastic FEM of metal forming simulation. In the FEM, the direct iteration and Newton-Raphson iteration are applied to obtain the initial solution and accurate solution respectively. In general, the acceleration scheme for the direct iteration is not used. In this paper, an Aitken accelerator is applied to the direct iteration. In the modified Newton-Raphson iteration, the step length or the deceleration coefficient is used for the fast and robust convergence. The step length can be determined by using the accelerator. The numerical experiments have been performed for the comparisons. The faster convergence is obtained with the acceleration in the direct and Newton-Raphson iterations.

• The Korean Journal of Applied Statistics
• /
• v.6 no.2
• /
• pp.319-328
• /
• 1993

The actual convergence rate of Newton-Raphson iteration method at each step is studied under the regularity conditions for the limiting distribution: The convergence rate of it is accelerated with good starting values. Hence we can decide a number of iterations according to our purposes.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.18 no.4 s.70
• /
• pp.385-393
• /
• 2005

This paper established the dynamic model of a flexible Timoshenko beam capable of geometrical nonlinearities subject to large overall motions by using the finite element method. Equations of motion are derived by using Hamilton principle and are formulated in terms of finite elements in which the nonlinear constraint equations are adjoined to the system using Lagrange multipliers. The Newmark direct integration method and the Newton-Raphson iteration are employed here for the numerical study which is to demonstrate the efficiency of the proposed formulation.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2008.05a
• /
• pp.640-644
• /
• 2008

하상경사가 커서 동수역학적 부정류 계산모형을 안정적으로 적용하기 어렵고, 홍수파의 감쇄효과가 적은 중소하천에 적합한 준정상류 계산모형을 개발하였다. 수립된 모형은 매 시각 유량에 대하여 1차원 하천 부등류 지배방정식인 단면 평균된 1차원 에너지 방정식을 풀도록 구성되어 있으며, 수치해법으로는 Newton-Raphson 방법을 적용한 표준축차법을 사용하였다. Newton-Raphson 방법을 적용하기 위해서는 통수면적, 하폭, 윤변, 동수반경 및 수위에 대한 윤변의 변화율 등의 변수들이 필요하다. 이와 같은 변수들은 각 계산점에서 수위를 계산하기에 앞서 단면자료를 사용하여 0.1 m 간격으로 모든 수위에 대하여 그 값들을 미리 구한 후, 반복 계산 단계에서 사용되는 수위에 대하여 필요한 변수들을 앞서 계산된 변수들과 선형 보간하여 사용하도록 하였다. 하천 구간내에 보가 존재하는 경우에는 보가 위치한 상 하류 간의 지배방정식으로 에너지 방정식 대신에 월류 유량 관계식을 사용하였으며, 이때의 수치해법 역시 Newton-Raphson 방법을 사용하였다. 수립된 모형을 한탄강 하류 구간에 적용하여 HEC-RAS 모형과 모의 결과를 비교한 결과, 두 모형의 계산결과가 잘 일치하는 것으로 나타났다. 에너지 경사항의 근사 방법에 따른 민감도 분석을 실시하였다.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.12 no.4
• /
• pp.642-651
• /
• 1988

The present paper develops finite element procedures to calculate displacements, strains and stresses in initially stressed elastic solids subjected to static or time-dependent loading conditions. As a point of departure, we employ Hamilton's principle to obtain nonlinear equations of motion characterizing the displacement in a solid. The equations of motion reduce to linear equations of motion if incremental stresses are assumed to be infinitesimal. In the case of linear problem, finite element solutions are obtained by Newmark's direct integration method and by modal analysis. An analytic solution is referred to compare with the linear finite element solution. In the case of nonlinear problem, finite element solutions are obtained by Newton-Raphson iteration method and compared with the linear solution. Finally, the effect of the order of Gauss-Legendre numerical integration on the nonlinear finite element solution, has been investigated.

• The Transactions of the Korean Institute of Electrical Engineers A
• /
• v.51 no.3
• /
• pp.127-133
• /
• 2002

This paper presents an efficient improvement of the iterative eigenvalue calculation method of the AESOPS algorithm. The intuitively and heuristically approximated iterative eigenvalue calculation method of the AESOPS algorithm is transformed to the Second Order Newton Raphson Method which is generally used in numerical analysis. The equations of second order partial differentiation of external torque, terminal and internal voltages are derived from the original AESOPS algorithm. Therefore only a few calculation steps are added to transform the intuitively and heuristically approximated AESOPS algorithm to the Second Order Newton Raphson Method, while the merits of original algorithm are still preserved.

• Journal of the Society of Naval Architects of Korea
• /
• v.28 no.1
• /
• pp.182-188
• /
• 1991

The modified arc-length algorithms for the automatic incremental solution of nonlinear finite element equations proposed by Riks are presented, which comprise the cylindrical arc-length method and the normal arc-length method. These methods are developed to trace the nonlinear path of large displacement problems such as a pre and post bucking/collapse response of general structures. These methods are applied to analyse the nonlinear behavior of arch and shell problems in parallel with the standard and modified Newton-Raphson method.

• Economic and Environmental Geology
• /
• v.22 no.3
• /
• pp.267-276
• /
• 1989

This paper presents results of interpretaton of gravity data by iterative nonlinear inversion methods. The gravity data are obtained by a theoretical formula for two-dimensional 2-layer structure. Depths to the basement of the structure are determined from the gravity data by four interative inversion methods. The four inversion methods used here are the Gradient, Gauss-Newton, Newton-Raphson, and Full Newton methods. Inversions are performed by using different initial guesses of depth for the over-determined, even-determined, and under-determined cases. This study shows that the depth can be determined well by all of the methods and most efficiently by the Newton-Raphson method.

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