• Computers and Concrete
• /
• v.6 no.6
• /
• pp.505-521
• /
• 2009

The analysis of prismatic members made of reinforced concrete under inclined bending, especially the computation of ultimate loads, is a pronounced non-linear problem which is frequently solved by discretizing the stress distribution in the cross-section using interpolation functions. In the approach described in the present contribution the exact analytical stress distribution is used instead. The obtained expressions are integrated by means of a symbolic manipulation package and automatically converted to optimized Fortran code. The direct problem-computation of ultimate internal forces given the position of the neutral axis-is first described. Subsequently, two kinds of inverse problem are treated: the computation of rupture envelops and the dimensioning of reinforcement, given design internal forces. An iterative Newton-Raphson procedure is used. Examples are presented.

• Computers and Concrete
• /
• v.13 no.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.

• Interaction and multiscale mechanics
• /
• v.1 no.1
• /
• pp.83-101
• /
• 2008

A derivation is given of two-scale models that are able to describe deformation and flow in a fluid-saturated and progressively fracturing porous medium. From the micromechanics of the flow in the cavity, identities are derived that couple the local momentum and the mass balances to the governing equations for a fluid-saturated porous medium, which are assumed to hold on the macroscopic scale. By exploiting the partition-of-unity property of the finite element shape functions, the position and direction of the fractures are independent from the underlying discretization. The finite element equations are derived for this two-scale approach and integrated over time. The resulting discrete equations are nonlinear due to the cohesive crack model and the nonlinearity of the coupling terms. A consistent linearization is given for use within a Newton-Raphson iterative procedure. Finally, examples are given to show the versatility and the efficiency of the approach.

• Structural Engineering and Mechanics
• /
• v.12 no.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.06a
• /
• 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.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 2000.11a
• /
• pp.445-450
• /
• 2000

In this paper, a numerical method is proposed for cam synthesis. kinematics of closed loop system with cam and follower is presented using relative coordinates. The system is transformed into an open loop system by cutting fictitiously higher-pair contact of cam and follower and envelope constraint equations are derived. Follower constraint equations are derived from the motion of the follower ends. The joint variables and follower profile parameters are calculated from the envelope constraint equations and follower constraint equations by using the Newton - Raphson iterative method. Algorithms for cam synthesis are presented and simulations are done to verify the effectiveness of the proposed method.

• Proceedings of the KIEE Conference
• /
• 2004.07d
• /
• pp.2242-2244
• /
• 2004

This paper deals with the design of a reduced-order stabilizing controller for the linear system. The coupled lineal matrix inequality (LMI) problem subject to a rank condition is solved by a sequential semidefinite programming (SDP) approach. The nonconvex rank constraint is incorporated into a strictly linear penalty function, and the computation of the gradient and Hessian function for the Newton method is not required. The penalty factor and related term are updated iteratively. Therefore the overall procedure leads to a successive LMI relaxation method. Extensive numerical experiments illustrate the proposed algorithm.

• Journal of applied mathematics & informatics
• /
• v.6 no.3
• /
• pp.697-726
• /
• 1999

In this paper we develop a now procedure to control stepsize for linear multistep methods applied to semi-explicit index 1 differential-algebraic equations. in contrast to the standard approach the error control mechanism presented here is based on monitoring and contolling both the local and global errors of multistep formulas. As a result such methods with the local-global stepsize control solve differential-algebraic equation with any prescribed accuracy (up to round-off errors). For implicit multistep methods we give the minimum number of both full and modified Newton iterations allowing the iterative approxima-tions to be correctly used in the procedure of the local-global stepsize control. We also discuss validity of simple iterations for high accuracy solving differential-algebraic equations. Numerical tests support the the-oretical results of the paper.

• Journal of applied mathematics & informatics
• /
• v.8 no.3
• /
• pp.935-942
• /
• 2001

A finite element approximation of a fourth-order nonlinear boundary value problem is given. In the direct implementation, a nonlinear system will be obtained and also a full size matrix will be introduced when Newton’s method is adopted to solve the system. To avoid this difficulty we introduce an iterative scheme which can be shown to converge the positive solution of the system lying between 0 and $sin{\pi}x$.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.5
• /
• pp.1640-1648
• /
• 1991

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