• Computational Structural Engineering
• /
• v.2 no.3
• /
• pp.97-103
• /
• 1989

In this study, a usefulness of the iterative solution procedures is reviewed for the elasto-plastic large deflection analysis of imperfect plates by finite element method. Three typical solution techniques such as simple incremental(SI) method, Newton-Raphson(NR) method and modified Newton-Raphson (mNR) method are compared. It is concluded that for thin plates which are given rise to the large deflection, iteration for the convergence of the unbalance force should be performed and in this case mNR method is more useful than NR method since the computing time of the former becomes to be a half of the latter, in which the accuracy of the result remains same. For thick plates or thin plates with large initial deflection, however, the use of SI method is quite better since the unbalance force may be negligible.

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

• Journal of Korea Water Resources Association
• /
• v.40 no.10
• /
• pp.801-810
• /
• 2007

In this study, validity and limitation of the estimation of roughness coefficient using the measured field data are investigated and the errors of the calculated roughness coefficient are analyzed. The assumption of uniform flow led to much difference of the computed results in low flow, and this is due to change of the cross-section informations such as flow area and hydraulic radius rather than the difference of velocity head. From the comparison between the estimations of average roughness coefficient in the reach which is relatively long, the calculation using the modified Newton-Raphson method is very efficient and accurate. In the measured roughness coefficient, the errors of measured flow and stage are included and the lower flow is, the larger the magnitude of error of measured roughness coefficient is. But the error of depth and velocity associated with uncertainty of roughness coefficient is less than about 5% in the both of low and high flow, and it shows the validity of measured roughness coefficient.

• Journal of Korea Water Resources Association
• /
• v.34 no.5
• /
• pp.543-549
• /
• 2001

This research estimated stream discharges indirectly by hydraulic channel routing. Only stage data from three stage stations and river cross section data were used to estimate Manning roughness coefficients and to compute stream discharges. When the discharges were estimated a stage-stage set of conditions was used for upstream-downstream boundary conditions. The research used the data from the upper Mississippi River. The hydraulic channel routings were performed by DWOPER (operational dynamic wave model). The global optimization program of SCE-UA was used to improve the roughness coefficient estimation module of the modified Newton-Raphson method in DWOPER. The results from SCE-US were better. For the case study of a flood, most estimated discharges except a few show errors within 10%.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.15 no.3
• /
• pp.481-490
• /
• 2002

The purpose of this paper is to evaluate the efficiency of various solution techniques and propose new efficient solution techniques for space trusses. Solution techniques used in this study are three load control methods (Newton-Raphson Method, modified Newton-Raphson Method, Secant-Newton Method), two load-displacement control methods(Arc-length Method, Work Increment Control Method) and three combined load-displacement control methods(Combined Arc-length Method I , Combined Arc-length MethodⅡ, Combined Work Increment Control Method). To evaluate the efficiency of these solution techniques, we must examine accuracy of their solutions, convergences and computing times of numerical examples. The combined load-displacement control methods are the most efficient in the geometric nonlinear solution techniques and in tracing post-buckling behavior of space truss. The combined work increment control method is the most efficient in tracing the buckling load of spate trusses with high degrees of freedom.

• Journal of IKEEE
• /
• v.18 no.2
• /
• pp.272-281
• /
• 2014

In electrical impedance tomography (EIT), modified Newton Raphson (mNR) method is widely used inverse algorithm for static image reconstruction due to its convergence speed and estimation accuracy. The unknown conductivity distribution is estimated iteratively by minimizing a cost functional such that the residual error namely the difference in measured and calculated voltages is reduced. Although, mNR method has good estimation performance, EIT inverse problem still suffers from ill-conditioned and ill-posedness nature. To mitigate the ill-posedness, generally, regularization methods are adopted. The inverse solution is highly dependent on the choice of regularization parameter. In most cases, the regularization parameter has a constant value and is chosen based on experience or trail and error approach. In situations, when the internal distribution changes or with high measurement noise, the solution does not get converged with the use of constant regularization parameter. Therefore, in this paper, in order to improve the image reconstruction performance, we propose a new scheme to determine the regularization parameter. The regularization parameter is computed based on residual error and updated every iteration. The proposed scheme is tested with numerical simulations and laboratory phantom experiments. The results show an improved reconstruction performance when using the proposed regularization scheme as compared to constant regularization scheme.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.30 no.7
• /
• pp.44-50
• /
• 2002

Low-velocity impact on composite sandwich panel has been investigated. Contact force is computed from a proposed modified Hertzian contact law. The Hertzian contact law is constructed by adjusting numerical value of the exponent and reducing the through-the- thickness elastic constant of honeycomb core. The equivalent transverse elastic constant is calculated from the rule of mixture. Nonlinear equation to calculate the contact force is solved by the Newton-Raphson method and time integration is done by the Newmark-beta method. A finite element program for the low-velocity impact analysis is coded by implementing these techniques and an 18-node assumed strain solid element. Behaviors of composite sandwich panels subjected to low-velocity impact are analyzed for various cases with different geometry and lay-ups. It has been found that the present code with the proposed contact law can predict measured contact forces and contact times for most cases within reasonable error bounds.

• Journal of Korean Society of Steel Construction
• /
• v.15 no.4 s.65
• /
• pp.389-401
• /
• 2003

The purpose of this study is to develop an algorithm, which can be designed the optimal sections of the composite framed structures constituted with the beams and the columns consisted of H type of steel section and concrete considering the reliability index. The optimized problem or the composite framed structures is formulated with the objective function and the constraints taking the section sizes as the design variables. The objective functions are constituted by the total costs of constructions. Also, the constraints are derived by considering the reliability index of section stress and allowable stress. The algorithm optimized the section of the composite framed structures utilizes the SUMT method using the modified Newton-Raphson direction method. The optimizing algorithm developed in this study is applied to the numerical examples with respecting a one-bay, one-story composite framed structure and a one-bay five-story one for the practical utilization of design on the composite framed structures using the reliability indices$({\beta})$ three and zero. In addition, their numerical results are compared and analyzed to examine the possibility of optimization the applicability, and the convergence this algorithm.

• Journal of the Society of Disaster Information
• /
• v.17 no.4
• /
• pp.839-847
• /
• 2021

Purpose: This study investigates mechanical behavior of functionally graded (FG) carbon nanotube-reinforced composite (CNTRC) plate in flexure. Isogeometric analysis (IGA) method coupled with shear deformable theory of higher-order (HSDT) to analyze the nonlinear bending response is presented. Method: Shear deformable plate theory into which a polynomial shear shape function and the von Karman type geometric nonlinearity are incorporated is used to derive the nonlinear equations of equilibrium for FG-CNTRC plate in bending. The modified Newton-Raphson iteration is adopted to solve the system equations. Result: The dispersion pattern of carbon nanotubes, plate geometric parameter and boundary condition have significant effects on the nonlinear flexural behavior of FG-CNTRC plate. Conclusion: The proposed IGA method coupled with the HSDT can successfully predict the flexural behavior of FG-CNTRC plate.

• Proceedings of the KIEE Conference
• /
• 1997.07c
• /
• pp.1094-1096
• /
• 1997

이제까지 전력계통의 상태를 알기 위한 조류계산은 전력방정식에 대한 직접적인 풀이로써 이루어졌다. 그런데 만일 자코비안 행열이 singular나 유사 singular가 된다면 그 전력 방정식의 해를 구할 수가 없게 된다. singular나 유사 singular가 되는 자코비안 행열을 가지는 전력방정식을 풀기 위하여 보조 방정식을 추가하여 미분불능인 점을 미분가능으로 변환하는 continuous method가 있다. 그런데 continuous method에서 보조 방정식은 원함수의 성질에 따라 달라지므로 어떤 일반적인 형태를 가지지 않는다. 따라서 본 논문에서는 전력방정식을 제약조건으로 가지고 부하 모선의 전력이 수용가의 요구량과 거의 일치되도록 하는 것을 목적 함수로 하는 최적화 문제로 조류계산을 변환하여 풀이하였다. 이러한 최적화 기법을 이용함으로써 continuous method에서 보조방정식이 원함수에 따라 계속적으로 변하는 문제를 개선하였다.