• 대한수학회지
• /
• 제34권1호
• /
• pp.87-118
• /
• 1997

This paper is concerned with a mixed Lagrangian formulation of the wiscous, stationary, incompressible Navier-Stokes equations $$ (1.1) -\nu\Delta u + (u \cdot \nabla)u + \nabla_p = f in \Omega $$ and $$ (1.2) \nubla \cdot u = 0 in \Omega $$ along with inhomogeneous Dirichlet boundary conditions on a portion of the boundary $$ (1.3) u = ^{0 on \Gamma_0 _{g on \Gamma_g, $$ where $\Omega$ is a bounded open domain in $R^d, d = 2 or 3$, or with a boundary $\Gamma = \partial\Omega$, which is composed of two disjoint parts $\Gamma_0$ and $\Gamma_g$.

• Coupled systems mechanics
• /
• 제7권5호
• /
• pp.583-610
• /
• 2018

This paper discusses the issues associated with modeling frictional contact between solid bodies undergoing large deformations. The most common model for friction on contact interfaces in solid mechanics is the Coulomb friction model, in which two distinct responses are possible: stick and slip. Handling the transition between these two phases computationally has been a source of algorithmic instability, lack of convergence and non-unique solutions, particularly in the presence of large deformations. Most computational models for frictional contact have used penalty or updated Lagrangian approaches to enforce frictional contact conditions. These two approaches, however, present some computational challenges due to conditioning issues in penalty-type implementations and the iterative nature of the updated Lagrangian formulation, which, particularly in large simulations, may lead to relatively slow convergence. Alternatively, a plasticity-inspired implementation of frictional contact has been shown to handle the stick-slip conditions in a local, algorithmically efficient manner that substantially reduces computational cost and successfully avoids the issues of instability and lack of convergence often reported with other methods (Laursen and Simo 1993). The formulation of this approach, however, has been limited to the small deformations realm, a fact that severely limited its application to contact problems where large deformations are expected. In this paper, we present an algorithmically consistent formulation of this method that preserves its key advantages, while extending its application to the realm of large-deformation contact problems. We show that the method produces results similar to the augmented Lagrangian formulation at a reduced computational cost.

• International Journal of Aeronautical and Space Sciences
• /
• 제16권2호
• /
• pp.177-189
• /
• 2015

This paper presents an all-direct domain decomposition approach for large-scale structural analysis. The proposed approach achieves computational robustness and efficiency by enforcing the compatibility of the displacement field across the sub-domain boundaries via local Lagrange multipliers and augmented Lagrangian formulation (ALF). The proposed domain decomposition approach was compared to the existing FETI approach in terms of the computational time and memory usage. The parallel implementation of the proposed algorithm was described in detail. Finally, a preliminary validation was attempted for the proposed approach, and the numerical results of two- and three-dimensional problems were compared to those obtained through a dual-primal FETI approach. The results indicate an improvement in the performance as a result of the implementing the proposed approach.

• 한국군사과학기술학회지
• /
• 제12권3호
• /
• pp.272-281
• /
• 2009

In this paper, three-dimensional optimal evasive maneuver patterns for air-to-surface attack missiles against proportionally navigated anti-air defense missiles were investigated. An interception error of the defense missile is produced by an evasive maneuver of the attack missile. It is assumed that the defense missiles are continuously launched during the flight of attack missile. The performance index to be minimized is then defined as the negative square integral of the interception errors. The direct parameter optimization technique based on SQP and a co-evolution method based on the augmented Lagrangian formulation are adopted to get the attack missile's optimal evasive maneuver patterns. The overall shape of the resultant optimal evasive maneuver is represented as a deformed barrel-roll.

• 한국전산구조공학회논문집
• /
• 제25권6호
• /
• pp.469-476
• /
• 2012

본 논문에서는 대규모 구조해석을 위하여 국부(local) 및 전역-국부 혼합(mixed) Lagrange 승수(Lagrange multiplier)를 이용한 새로운 유한요소 영역분할 기법을 제시한다. 제시되는 FETI 알고리즘은 계산 효율성을 향상시키기 위하여 기존의 FETI 기법들에서 사용되어 온 전통적인 Lagrange 승수법과는 달리, 국부 및 전역-국부 혼합 Lagrange 승수를 도입하고 ALF(Augmented Lagrangian Formulation)과의 결합을 유도하여 공유면 문제(interface problem)의 해의 수렴성을 향상 시켰다. 추가적으로, 몇 가지 수치예제 계산을 통해 기존의 FETI-DP 기법과 비교하여 유연도 행렬의 조건수, 계산 시간 그리고 메모리 사용량에 대한 계산결과를 제시하였다.

• Journal of Magnetics
• /
• 제16권1호
• /
• pp.77-82
• /
• 2011

This paper presents a new self-adjoint material sensitivity formulation for optimal designs and inverse problems in the high frequency domain. The proposed method is based on the continuum approach using the augmented Lagrangian method. Using the self-adjoint formulation, there is no need to solve the adjoint system additionally when the goal function is a function of the S-parameter. In addition, the algorithm is more general than most previous approaches because it is independent of specific analysis methods or gridding techniques, thereby enabling the use of commercial EM simulators and various custom solvers. For verification, the method was applied to the several numerical examples of dielectric material reconstruction problems in the high frequency domain, and the results were compared with those calculated using the conventional method.

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

• 대한전기학회:학술대회논문집
• /
• 대한전기학회 2007년도 추계학술대회 논문집 전력기술부문
• /
• pp.226-228
• /
• 2007

This paper proposes Nonlinear Interior Point Method (NIPM) including discrete control variables optimal power flow formulations. The algorithm utilizes the robustness in terms of starting point and fast convergence for large scale power system of NIPM and an introduction of rounding penalty function which is augmented in the Lagrangian function to handle discrete control variables. The derived formulation shows a simplified approach to deal with discrete control problems which is implementable in real large scale systems.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제18권1호
• /
• pp.17-26
• /
• 2014

This paper presents development of an improved domain decomposition method for large scale structural problem that aims to provide high computational efficiency. In the previous researches, we developed the domain decomposition algorithm based on augmented Lagrangian formulation and proved numerical efficiency under both serial and parallel computing environment. In this paper, new computational analysis by the proposed domain decomposition method is performed. For this purpose, reduction in computational time achieved by the proposed algorithm is compared with that obtained by the dual-primal FETI method under serial computing condition. It is found that the proposed methods significantly accelerate the computational speed for a linear structural problem.

• Journal of Mechanical Science and Technology
• /
• 제19권spc1호
• /
• pp.388-394
• /
• 2005

Recently, the authors have developed a method for real-time dynamics of multibody systems, which combines a semi-recursive formulation to derive the equations of motion in dependent relative coordinates, along with an augmented Lagrangian technique to impose the loop closure conditions. The following numerical integration procedures, which can be grouped into the so-called structural integrators, were tested : trapezoidal rule, Newmark dissipative schemes, HHT rule, and the Generalized-${\alpha}$ family. It was shown that, for large multi body systems, Newmark dissipative was the best election since, provided that the adequate parameters were chosen, excellent behavior was achieved in terms of efficiency and robustness with acceptable levels of accuracy. In the present paper, the performance of the described method in combination with another group of integrators, the Implicit Runge-Kutta family (IRK), is analyzed. The purpose is to clarify which kind of IRK algorithms can be more suitable for real-time applications, and to see whether they can be competitive with the already tested structural family of integrators. The final objective of the work is to provide some practical criteria for those interested in achieving real-time performance for large and complex multibody systems.