• Nuclear Engineering and Technology
• /
• 제55권1호
• /
• pp.310-323
• /
• 2023

The k-eigenvalue neutronics/thermal-hydraulics coupling calculation is a key issue for reactor design and analysis. Jacobian-free Newton-Krylov (JFNK) method, featured with super-linear convergence rate and high efficiency, has been attracting more and more attention to solve the multi-physics coupling problem. However, it may converge to the high-order eigenmode because of the multiple solutions nature of the k-eigenvalue form of multi-physics coupling issue. Based on our previous work, a modified JFNK with a line search method is proposed in this work, which can find the fundamental eigenmode together with thermal-hydraulics feedback in a wide range of initial values. In detail, the existing modified JFNK method is combined with the line search strategy, so that the intermediate iterative solution can avoid a sudden divergence and be adjusted into a convergence basin smoothly. Two simplified 2-D homogeneous reactor models, a PWR model, and an HTR model, are utilized to evaluate the performance of the newly proposed JFNK method. The results show that the performance of this proposed JFNK is more robust than the existing JFNK-based methods.

• 지구물리와물리탐사
• /
• 제7권3호
• /
• pp.207-212
• /
• 2004

자기지전류(MT) 자료의 3차원 역산에 대해 소개한다. MT 자료의 역산 문제는 기본적으로 악조건이므로 유일한 해가 존재하지 않는다. 이러한 비유일성을 줄이고 정확한 역산해를 구하기 위해서는 역산 시 사전정보를 추가하는 제약조건을 가해야 한다. 물리탐사 분야에서 비선형 역산에 사용되는 가장 일반적인 방법은 일련의 선형화된 역산문제를 푸는 Gauss-Newton법이다. 이 알고리듬은 수렴 시, 모델 공간에서 역산문제에 대한 목적함수를 최소화하는 최적해를 준다. 그러나 이러한 반복적 선형화기법은 3차원 MT 역산의 경우 Jacobian 행렬을 구하기 힘들기 때문에 그 유용성에 한계가 있다. 이러한 어려움은 CG법에 의해 완화할 수 있다. 선형 CG법은 Gauss-Newton 반복의 각 단계를 근사적으로 풀기 위해서 사용된다. 한편 비선형 CG법은 목적함수의 최소화에 직접적으로 적용된다. 이들 CG법은 Jacobian 행렬의 계산 및 대형 선형방정식의 해를 반복 당 세 번의 모델링으로 대치할 수 있어서 3차원 역산에 적합하다.

• Tribology and Lubricants
• /
• 제14권2호
• /
• pp.26-34
• /
• 1998

This paper presents a high-speed performance analysis of ball bearings with ceramic balls under thrust loads. The sliding velocity profiles between a ball and raceways were obtained by the 3-D quasi-dynamic equations of motion including both centrifugal force and gyroscopic moment derived by vector matrix algebra. The friction at the contact areas was obtained by the Bair-Winer's non-Newtonian rheological model and the Hamrock-Dowson's central film thickness in EHL analysis. The nonlinear equations were solved by the Newton-Raphson method and the underrelaxation iterative method. The friction torques and ball behaviors with various loads, ball materials, and contact angles were predicted by this model. It was shown that the friction torque was sensitive to thrust load and contact angle, and that the friction torque and the pitch angle of the bearing with ceramic balls are smaller than those of the bearing with steel balls.

• Steel and Composite Structures
• /
• 제21권5호
• /
• pp.1121-1144
• /
• 2016

This paper presents a displacement-based finite element procedure for second-order distributed plasticity analysis of planar steel frames with semi-rigid beam-to-column connections under static loadings. A partially strain-hardening elastic-plastic beam-column element, which directly takes into account geometric nonlinearity, gradual yielding of material, and flexibility of semi-rigid connections, is proposed. The second-order effects and distributed plasticity are considered by dividing the member into several sub-elements and meshing the cross-section into several fibers. A new nonlinear solution procedure based on the combination of the Newton-Raphson equilibrium iterative algorithm and the constant work method for adjusting the incremental load factor is proposed for solving nonlinear equilibrium equations. The nonlinear inelastic behavior predicted by the proposed program compares well with previous studies. Coupling effects of three primary sources of nonlinearity, geometric imperfections, and residual stress are investigated and discussed in this paper.

• Computers and Concrete
• /
• 제13권3호
• /
• pp.343-359
• /
• 2014

In this work a numerical method to simulate the response of reinforced concrete structures subject to cyclically imposed displacements is proposed. The method consists of a combination of a displacement and load controlled version of the Newton-Raphson iterative technique, used for the loading and the unloading part of the cycles respectively. The whole procedure is combined with a relatively simple concrete model whose only material parameter is its uniaxial compressive strength. The proposed methodology may realistically simulate, in an easy way, the physical process of any experimentally tested RC structure under imposed displacements cycles. The efficiency of the approach is demonstrated through a series of analyses of experimentally tested specimens reported in the literature.

• Structural Engineering and Mechanics
• /
• 제12권1호
• /
• pp.35-50
• /
• 2001

A numerical methodology is presented in this paper for the geometrically non-linear analysis of slender uni-dimensional structural elements under unilateral contact constraints. The finite element method together with an updated Lagrangian formulation is used to study the structural system. The unilateral constraints are imposed by tensionless supports or foundations. At each load step, in order to obtain the contact regions, the equilibrium equations are linearized and the contact problem is treated directly as a minimisation problem with inequality constraints, resulting in a linear complementarity problem (LCP). After the resulting LCP is solved by Lemke's pivoting algorithm, the contact regions are identified and the Newton-Raphson method is used together with path following methods to obtain the new contact forces and equilibrium configurations. The proposed methodology is illustrated by two examples and the results are compared with numerical and experimental results found in literature.

• 제어로봇시스템학회:학술대회논문집
• /
• 제어로봇시스템학회 2005년도 ICCAS
• /
• pp.906-909
• /
• 2005

In this study, the dynamic flying characteristics of the worn head sliders are investigated theoretically due to the change in head geometry caused by head and disk contact. The film shapes can be approximated as taper- truncated cycloidal-flat film. Two-dimensional time dependent modified Reynolds equation included molecular slip effect are formulated with neglected the roughness effect. The motion of head slider was assumed to have two degree of freedom in this work. Finite difference approximation with Newton Raphson iterative technique and the fourth order Runge-Kutta method were implemented to obtain the transient response of the slider head with various change in head geometry numerically and compared with the transient response of the IBM3380 type head slider. The simulation results show the film shape has affects significantly on the static and dynamic characteristic of slider head in magnetic storage systems.

• 대한기계학회논문집
• /
• 제15권5호
• /
• pp.1640-1648
• /
• 1991

본 연구에서는 이러한 편사 함수 최소화의 방법을 적용함에 있어 보다 안정된 수렴성과 계산 시간을 단축시키기 위하여 근접 위치 방법(near position method)을 개 발하여 적용하였다. 근접 위치 방법이란 이론적 해석법으로 풀기가 불가능한 기구학 을 갖는 6관절 로봇을 반복적 해석법을 사용한다는 것을 전제로 하여, 초기 위치를 목 표 위치에 가능한 근접하게 잡아서 반복 계산을 수행하는 방법으로써 로봇의 기구학적 자세에 따른 수렴의 불안정성을 방지하고, 계산 시간을 단축하는데 그 목적이 있다.

• Computers and Concrete
• /
• 제1권2호
• /
• pp.131-150
• /
• 2004

We present in this work a coupled phenomenological chemo-mechanical model that represents the degradation of concrete-like materials. The chemical behaviour is described by the nowadays well known simplified calcium leaching approach. And the mechanical damage behaviour is described by a continuum damage model which involves the gradient of the damage quantity. The coupled nonlinear problem at hand is addressed within the context of the finite element method. For the equation governing the calcium dissolution-diffusion part of the problem, special care is taken to treat the highly nonlinear calcium conductivity and solid calcium functions. The algorithmic design is based on a Newton-type iterative scheme where use is made of a recently proposed relaxed linearization procedure. And for the equation governing the damage part of the problem, an augmented Lagrangian formulation is used to take into account the damage irreversibility constraint. Finally, numerical simulations are compared with experimental results on cement paste.

• Structural Engineering and Mechanics
• /
• 제43권5호
• /
• pp.583-605
• /
• 2012

This paper presents a spatial catenary cable element for the nonlinear analysis of cable-supported structures. An incremental-iterative solution based on the Newton-Raphson method is adopted for solving the equilibrium equation. As a result, the element stiffness matrix and nodal forces are determined, wherein the effect of self-weight and pretension are taken into account. In the case of the initial cable tension is given, an algorithm for form-finding of cable-supported structures is proposed to determine precisely the unstressed length of the cables. Several classical numerical examples are solved and compared with the other available numerical methods or experiment tests showing the accuracy and efficiency of the present elements.