• Computational Structural Engineering
• /
• v.11 no.1
• /
• pp.179-190
• /
• 1998

A geometrically nonlinear finite element formulation of spatial cable networks is presented using two cable elements. Firstly, derivation procedures of tangent stiffness and mass matrices for the space truss element and the elastic catenary cable element are summarized. The load incremental method based on Newton-Raphson iteration method and the dynamic relaxation method are presented in order to determine the initial static state of cable nets subjected to self-weights and support motions. Furthermore, static non-linear analysis of cable structures under additional live loads are performed based on the initial configuration. Challenging example problems are presented and discussed in order to demonstrate the feasibility of the present finite element method and investigate static nonlinear behaviors of cable nets.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.11 no.3
• /
• pp.1-10
• /
• 1991

The paper is concerned with the analysis of laminated composite shell structures using an improved degenerated shell element. In the formulation of the element stiffness, the combined use of three different techniques was made. They are; 1) an enhanced interpolation of transverse shear strains in the natural coordinate system to overcome the shear locking problem; 2) the reduced integration technique in in-plane strains to avoid the membrane locking behavior; and 3) selective addition of the nonconforming displacement modes to improve the element performances. This element is free of serious shear/membrane locking problems and undesirable compatible/commutable spurious kinematic deformation modes. An incremental total Lagrangian formulation is presented which allows the calculation of arbitrarily large displacements. The resulting non-linear equilibrium equations are solved by the Newton-Raphson method. The versatility and accuracy of this improved degenerated shell element are demonstrated by solving several numerical examples.

• Proceedings of the KIEE Conference
• /
• 1999.07a
• /
• pp.70-72
• /
• 1999

일반적으로 비선형 정자계 문제를 해석하기 위해서 뉴튼-�N슨(Newton -Raphson : NR)법이 이용된다. 하지만 뉴튼-�N슨법의 경우 각 반복계산 때마다 새로운 선형 시스템의 해를 구하기 위해서 LU-decomposition과 같은 과정을 매 반복계산 때마다 시행해야 하므로 절점(node)의 수가 증가할 경우 계산시간이 증가한다는 단점이 있다. 이러한 단점을 보완하기 위해서 최근 TLM (Transmission Line Modeling)법이 새로운 반복계산법으로 비선형 유한 요소 해석에 적용되었으며 뉴튼-�N슨법에 비해 훨씬 우수한 특성을 보여주었다. 하지만 지금까지의 TLM법은 2차원의 정식화만 이루어졌고 3차원에는 적용되지 못한 것이 사실이다. 본 논문에서는 3차원의 비선형 정자계 문제에 TLM법을 적용할 수 있는 수식을 최초로 제안하며 3차원 코어(core)모델에 대해 TLM법을 적용하여 그 타당성을 검증하기로 한다. 또한 3차원 비선형 TLM법을 이용한 해석 결과가 뉴튼-�N슨법에 의한 결과와 완전히 일치하며 수렴 속도에 있어서도 훨씬 향상된 결과를 나타냄을 보이도록 하겠다.

• Tribology and Lubricants
• /
• v.2 no.1
• /
• pp.46-52
• /
• 1986

얇은 유막위의 타이어 변형해석이나, 압착막댐퍼 설계등의 해석에 응용되는 변동 하중하의 탄성 경계상 유체 압착막의 해석을 Newton-Raphson 반복법을 사용하여 하였다. 유막두께의 계산은 등계수 요소를 사용하여 정확한 계산을 하였고 슬라이더의 궤적은 강쇄진동의 응답곡선과 유사함을 알 수 있었다. 특히 본 연구에서는 미끄럼 속도의 영향을 고려하였으며 유체가 베어링내에 생긴 포켓에 정체하고 있어서 미끄럼의 영향은 속도가 클때를 제외하고는 큰 영향을 미치지 못함을 알 수 있었다.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.3
• /
• pp.799-808
• /
• 1991

본 연구에서는 연속체 역학의 에너지 원리에서 출발하여, 동적 비선형 해석을 위한 유한요소 식들을 유도하고, 이를 이용하여 특이 경계조건을 갖는 고체의 대변위 동적 선형 현상과 비선형 현상에 관하여 연구하고자 한다.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.27 no.1
• /
• pp.193-200
• /
• 2021

In this study, an analysis technique using static condensation is proposed for an efficient local nonlinear static analysis. The static condensation method is a model reduction method based on the degrees of freedom, and the analysis model is divided into a target part and a condensed part to be omitted. In this study, the nonlinear and linear parts were designated to the target and the omitted parts, respectively, and both the stiffness matrix and load vector corresponding to the linear part were condensed into the nonlinear part. After model condensation, the reduced model comprising the stiffness matrix and the load vector for the nonlinear part is constructed, and only this reduced model was updated through the Newton-Raphson iteration for an efficient nonlinear analysis. Finally, the efficiency and reliability of the proposed analysis technique were presented by applying it to various numerical examples.

• Composites Research
• /
• v.25 no.6
• /
• pp.217-223
• /
• 2012

This study investigates a geometrical nonlinear dynamic behaviors of laminated skew plates made of advanced composite materials (ACM). Based on the first-order shear deformation plate theory (FSDT), the Newmark method and Newton-Raphson iteration are used for the nonlinear dynamic solution. The effects of cutout sizes, skew angles and lay up sequences on the nonlinear dynamic response for various parameters are studied using a nonlinear dynamic finite element program developed for this study. The several numerical results were in good agreement with those reported by other investigators for square composite plates with or without central cutouts, and the new results reported in this paper show the significant interactions between the cutout, skew angles and layup sequence in the laminate. Key observation points are discussed and a brief design guideline of skew laminates is given.

• Computational Structural Engineering
• /
• v.10 no.4
• /
• pp.217-228
• /
• 1997

The high precision analysis by the p-version of the finite element method are fairly well established as highly efficient method for linear elastic problems, especially in the presence of stress singularity. It has been noted that the merits of the p-version are accuracy, modeling simplicity, robustness, and savings in user's and CPU time. However, little has been done to exploit their benefits in elasto-plastic analysis. In this paper, the p-version finite element model is proposed for the materially nonlinear analysis that is based on the incremental theory of plasticity using the constitutive equation for work-hardening materials, and the associated flow rule. To obtain the solution of nonlinear equation, the Newton-Raphson method and initial stiffness method, etc are used. Several numerical examples are tested with the help of the square plates with cutout, the thick-walled cylinder under internal pressure, and the circular plate with uniformly distributed load. Those results are compared with the theoretical solutions and the numerical solutions of ADINA

• Journal of the Korean Institute of Telematics and Electronics
• /
• v.27 no.3
• /
• pp.101-107
• /
• 1990

In this work, one-dimensional simulation is carried out for PT-GaAs Schottky barrier diodes with finite difference method. Shockley's semiconductor governing equations: Poisson equation and current continuity equation are discertized, and linearized by Newton-Raphson method. The linear system of equation is solved by Gaussian elimination method until convergence is achieved. The boundary condition for this equation is taken from thermionic emission-diffusion theory. Simulation is done for PT-GaAs epitaxial-layer Schottky barrier diodes. The claculated results of electron and potential distribution are shown. Simulation results show exellent agreement with experiments.

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