• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1998.10a
• /
• pp.431-438
• /
• 1998

An efficient solution method is presented to solve the eigenvalue problem arising in tile dynamic analysis of non-proportionally damped structural systems with multiple or close eigenvalues. The proposed method is obtained by applying the modified Newton-Raphson technique and the orthonormal condition of the eigenvectors to the quadratic eigenvalue problem. Even if the shift value is an eigenvalue of the system, the proposed method guarantees nonsingularity, which is analytically proved. The initial values of the proposed method can be taken as the intermediate results of iteration methods or results of approximate methods. Two numerical examples are also 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.

• Structural Engineering and Mechanics
• /
• v.36 no.5
• /
• pp.529-544
• /
• 2010

Nonlinear equations of structures are generally solved numerically by the iterative solution of linear equations. However, this iterative procedure diverges when the tangent stiffness is ill-conditioned which occurs near limit points. In other words, a major challenge with simple iterative methods is failure caused by a singular or near singular Jacobian matrix. In this paper, using the Newton-Raphson algorithm based on Davidenko's equations, the iterations can traverse the limit point without difficulty. It is argued that the propose algorithm may be both more computationally efficient and more robust compared to the other algorithm when tracing path through severe nonlinearities such as those associated with structural collapse. Two frames are analyzed using the proposed algorithm and the results are compared with the previous methods. The ability of the proposed method, particularly for tracing the limit points, is demonstrated by those numerical examples.

• Smart Structures and Systems
• /
• v.24 no.6
• /
• pp.683-692
• /
• 2019

Traditional damage detection methods for nonlinear structures are often based on simplified models, such as the mass-spring-damper and shear-building models, which are insufficient for predicting the vibration responses of a real structure. Conventional global nonlinear finite element model updating methods are computationally intensive and time consuming. Thus, they cannot be applied to practical structures. A decentralized approach for identifying the nonlinear material parameters is proposed in this study. With this technique, a structure is divided into several small zones on the basis of its structural configuration. The unknown material parameters and measured vibration responses are then divided into several subsets accordingly. The structural parameters of each subset are then updated using the vibration responses of the subset with the Newton-successive-over-relaxation (SOR) method. A reinforced concrete and steel frame structure subjected to earthquake loading is used to verify the effectiveness and accuracy of the proposed method. The parameters in the material constitutive model, such as compressive strength, initial tangent stiffness and yielding stress, are identified accurately and efficiently compared with the global nonlinear model updating approach.

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

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.10 no.1
• /
• pp.57-69
• /
• 1990

The static nonlinear analysis of offshore pipeline is carried out by the finite element method. The governing equilibrium equation are derived by the principle of minimum potential energy and the modified Newton-Raphson procedure is used to solve the system of nonlinear algebraic equation. Geometrically nonlinear beam elements and spring elements are utilized to model the pipeline, stinger, pipe supports and seabed simultaneously. The beam element developed can be used to model redundant structures. It provides for both the torsional deformation and elongation of pipeline, and permits the use of different physical properties in each principal direction. The validity of this method is investigated by comparing the results with these obtained by other methods.

• The Journal of Churna Manual Medicine for Spine and Nerves
• /
• v.7 no.2
• /
• pp.39-52
• /
• 2012

Objectives : The purpose of this study was to investigate the correlation between lumbar lordotic angle and the power of trunk flexors, extensors in normal adults Methods : 34 normal participants participated in this study. Their lumbar lordotic angle(L1-S1 Cobb's angle and L1-L5 cobb's angle) was measured by x-ray taken on lateral direction, erect cross-arm position. And muscle power of trunk flexors and extensors of each participant measured using Cybex HUMAC NORM. Results : 1. The average of L1-S1 Cobb's angle was $47.21{\pm}8.88^{\circ}$ and the average of L1-L5 Cobb's angle was $36.32{\pm}9.62^{\circ}$(Table IV). 2. The average ratio of trunk flexors/extensors was $6.44{\pm}19.31%$(Table V). The average power of the trunk flexors was $165.18{\pm}55.05$(Newton-Meter/kg), and the power of trunk extensors was $257.18{\pm}85.53$ (Newton-Meter/kg)(Table VI). 3. Lumbar lordotic angle has no relation to the ratio of trunk flexors/extensors(Table VII, Fig. 4). 4. Lumbar lordotic angle has no relation to both the power of the trunk flexors and extensors(Table VIII, Fig. 5, Fig. 6). Conclusions : These results suggest that the lumbar lordotic angle measured by radiograph could not evaluate the power and ratio of trunk flexors, extensors.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.6 no.2
• /
• pp.1-15
• /
• 1986

The algorithm Proposed utilizes the two-levels technique. In the first level which consists of two phases, the cross-sectional area of the truss member is optimized by transforming the nonlinear problem into SUMT, and solving SUMT utilizing the modified Newton-Rahson method. In the second level, the geometric shape is optimized utilizing the unindirectional search technique of the Powell method which make it possible to minimize only the objective function. The algorithm Proposed in this study is numerically tested for several truss structures with various shapes, loading conditions and design criteria, and compared with the results of the other algorithms to examine its applicability and stability. The numerical comparisons show that the two-Levels algorithm Proposed in this study is safely applicable to any design criteria, and the convergency rate is relathely fast and stable compared with other iteration methods for the geometric optimization of truss structures.

• Journal of the Korean GEO-environmental Society
• /
• v.15 no.3
• /
• pp.41-48
• /
• 2014

Experimental study of sedimentation and self-weight consolidation has been primary research area in dredged soil. However, good quality of the dredged soil and minimum water pollution caused by the pumping of reclaimed soil require intensive study of the flow characteristics of dredged material due to dumping. In this study, continuity and the equilibrium equations for mass flow assuming single phase was derived to simulate mass flow in dredged containment area. To optimize computation and modeling time for three dimensional geometry and boundary conditions, depth integration is applied to governing equations to consider three dimensional topography of the site. Petrov-Galerkin formulation is applied in spatial discretization of governing equations. Generalized trapezoidal rule is used for time integration, and Newton iteration process approximated the solution. DG and CDG technique were used for weighting matrix in discontinuous test function in dredged flow analysis, and numerical stability was evaluated by performed a square slump simulation. A comparative analysis for numerical methods showed that DG method applied to SU / PG formulation gives minimal pseudo oscillation and reliable numerical results.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.15 no.4
• /
• pp.661-674
• /
• 2002

The objective of this study is to develop efficient numerical method to enable solution of optimal control problems of Navier-Stokes flows and to apply these technique to the problem of viscous drag minimization on a bluff body by controlling boundary velocities on the surface of the body. In addition to the industrial importance of the drag reduction problem, it serves as a model for other more complex flow optimization settings, and allows us to study, modify, and improve the behavior of the optimal control methods proposed here. The control is affected by the suction or injection of fluid on portions of the boundary, and the objective function represents the rate at which energy is dissipated in the fluid. This study shows how reduced Hessian successive quadratic programming method, which avoid converging the flow equations at each iteration, can be tailored to these problems.

• Journal of Korean Society of Steel Construction
• /
• v.19 no.4
• /
• pp.345-355
• /
• 2007

During winter, the CWR (continuous welded rail) may be broken when a temperature drop below the neutral level changes the axial force, causing tensile fracture and creating a rail gap. The passage of a train on a rail with an open gap may lead to very costly derailments. In this paper, the use of a track-and-train-coupled model whose rail has an open gap is proposed for dynamic interaction analysis. Linear track and train systems were coupled in this study by a nonlinear Herzian contact spring, and the complete system matrices of the total track-train system were constructed. Moreover, the interaction phenomenon considering the presence of an open gap in the rail was toughly defined by assigning the irregularity functions between the two sides of the gap. Time history analysis, which has an iteration scheme such as the Newmark-$\beta$ method (based on the Modified Newton-Raphson methods), was conducted to solve the nonlinear equation. .Finally, numerical studies were conducted to assess the effect of the various parameters of the system when applied to various speeds, open-gap sizes, and support stiffnesses of the rail.