• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.4 no.2
• /
• pp.99-110
• /
• 2000

This paper is devoted to the parallelism of a numerical matrix eigenvalue problem. The eigenproblem arises in a variety of applications, including engineering, statistics, and economics. Especially we try to approach the industrial techniques from mathematical modeling. This paper has developed a parallel algorithm to find all eigenvalues. It is contributed to solve a specific practical problem, a vibration problem in the industry. Also we compare the runtime between the serial algorithm and the parallel algorithm for the given problems.

• Proceedings of the KIEE Conference
• /
• 2002.07d
• /
• pp.2721-2723
• /
• 2002

Electrical impedance tomograpy(EIT) determines the resistivity distribution inside an inhomogeneous target by means of voltage and current measurements conducted at the target boundary. In this paper, a simultaneous perturbation stochastic approximation(SPSA) approach is proposed for the solution of the EIT image reconstruction. Results of numerical experiments of EIT solved by the SPSA approach are presented and compared to that obtained by the modified Newton-Raphson(mNR) method.

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

• Advances in nano research
• /
• v.13 no.4
• /
• pp.393-406
• /
• 2022

In the present article, functionally graded small-scaled plates based on modified strain gradient theory (MSGT) are studied for analyzing the nonlinear bending and post-buckling responses. Von-Karman's assumptions are applied to incorporate geometric nonlinearity and the first-order shear deformation theory is used to model the plates. Modified strain gradient theory includes three length scale parameters and is reduced to the modified couple stress theory (MCST) and the classical theory (CT) if two or all three length scale parameters become zero, respectively. The Ritz method with Legendre polynomials are used to approximate the unknown displacement fields. The solution is found by the minimization of the total potential energy and the well-known Newton-Raphson technique is used to solve the nonlinear system of equations. In addition, numerical results for the functionally graded small-scaled plates are obtained and the effects of different boundary conditions, material gradient index, thickness to length scale parameter and length to thickness ratio of the plates on nonlinear bending and post-buckling responses are investigated and discussed.

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

• Structural Engineering and Mechanics
• /
• v.2 no.4
• /
• pp.395-402
• /
• 1994

An improved finite element-transfer matrix method is applied to the transient analysis of plates with large displacement under various excitations. In the present method, the transfer of state vectors from left to right in a combined finite element-transfer matrix method is changed into the transfer of generally incremental stiffness equations of every section from left to right. Furthermore, in this method, the propagation of round-off errors occurring in recursive multiplications of transfer and point matrices is avoided. The Newmark-${\beta}$ method is employed for time integration and the modified Newton-Raphson method for equilibrium iteration in each time step. An ITNONDL-W program based on this method using the IBM-PC/AT microcomputer is developed. Finally numerical examples are presented to demonstrate the accuracy as well as the potential of the proposed method for dynamic large deflection analysis of plates with random boundaries under various excitations.

• Structural Engineering and Mechanics
• /
• v.8 no.6
• /
• pp.561-578
• /
• 1999

The response of a reinforced concrete structure to loading is both immediate and time-dependent. Under a sustained load, the deflections caused by creep and shrinkage may be several times their instantaneous values. The paper describes a general finite element procedure, based on the so-called layered model, to analyse reinforced concrete members, and shows in particular how the simple Step by Step Method may be incorporated into this procedure. By invoking the Modified Newton Raphson Method as a solution procedure, the accuracy of the finite element method is verified against independent test results, and then applied to a variety of problems in order to demonstrate its efficacy. The method forms a general method for analysing highly indeterminate concrete structures in the time domain.

• 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.15 no.4
• /
• pp.395-414
• /
• 2003

In this paper we recall briefly the constitutive equations for solids subjected to thermal strain taking in account the bounded tensile stress of the material. In view to solve the equilibrium problem via the finite element method using the Newton Raphson procedure, we show that the tangent elasticity tensor is semi-definite positive. Therefore, in order to obtain a convergent numerical method, the constitutive equation needs to be modified. Specifically, the dependency of the stress by the anelastic deformation is made explicit by means of a parameter ${\delta}$, varying from 0 to 1, that factorizes the elastic tensor. This parameterization, for ${\delta}$ near to 0, assures the positiveness of the tangent elasticity tensor and enforces the convergence of the numerical method. Some numerical examples are illustrated.

• Journal of Biomedical Engineering Research
• /
• v.18 no.3
• /
• pp.211-216
• /
• 1997

An inverse method with regularization has been developed to determine the viscoelastic properties of trabecular bone. A series of stress relaxation experiments were performed under the condition of uniaxial compression stress state. Optimization has been formulated within the framework of nonlinear least-squares and a modified Gauss-Newton method with a zeroth-order regularization technique. The stress relaxation behavior of trabecular bone was analyzed using a standard viscoelastic model. The present study clearly shows that trabecular bone exhibits typical viscoelastic stress relaxation behavior, and the obtained material parameters well represent the viscoelastic behavior of trabecular bone.