• Structural Engineering and Mechanics
• /
• 제3권2호
• /
• pp.185-192
• /
• 1995

A numerical algorithm for plane stress and shell elasto-plasticity is presented in this paper. The proposed strain decomposition (SD) algorithm is an elastic predictor/plastic corrector algorithm, and in the context of operator splitting, is a return mapping algorithm. However, it differs significantly from other return mapping algorithms in that only the necessary response functions are used without invoking their gradients, and the stress increment is updated only at the end of the time step. This makes the proposed SD algorithm more suitable for materials with complex yield surfaces and will guard against error accumulation during the time step. Comparative analyses of structural systems using the proposed strain decomposition (SD) algorithm and the iterative radial return (IRR) algorithm are presented. The results demonstrate the accuracy and usefulness of the proposed algorithm.

• 대한기계학회논문집B
• /
• 제24권12호
• /
• pp.1570-1579
• /
• 2000

Finite element analysis of fluid flow with moving free surface has been carried out in two and tree dimensions. The new VOF-based numerical algorithm that has been proposed by the present authors was applied to several 2-D and 3-D free surface flow problems. The proposed free surface tracking scheme is based on two numerical tools that have been newly introduced by the present authots; the orientation vector to represent the free surface orientation in each cell and the baby-cell to determine the fluid volume flux at each cell boundary. The proposed numerical algorithm has been applied to 2-D and 3-D cavity filling and sloshing problems, which demonstrated versatility and effectiveness of the new free surface tracking scheme as well as the overall solution procedure. The proposed numerical algorithm resolved successfully the interacting free surface with each other. The simulated results demonstrated the applicability of proposed numerical algorithm to the practical problems of large free surface motion. Also, it has been demonstrated that the proposed free surface tracking scheme can be easily implemented in any irregular non-uniform grid systems and can be extended to the 3-D free surface flow problem without additional efforts.

• 정보처리학회논문지A
• /
• 제14A권1호
• /
• pp.73-76
• /
• 2007

단위원의 내부로부터 Jordan 영역으로의 등각사상을 구하는 것은 일반적으로 비선형방정식인 Theodorsen 방정식을 푸는 것으로 귀결된다. 저자는 이 비선형 방정식의 수치적 해법 중 가장 효율적인 방법으로 알려진 Wegmann의 해법에 저주파 필터를 적용하여 개선하고 새로운 산법의 수렴성을 이론적으로 증명한 바 있다[1, 2]. 또한 이 해법에 있어 참값을 모르더라도 오차평가가 가능한 방법을 제안하였다[3]. 본 논문에서는 참값을 모르더라도 오차평가가 가능한 연구결과를 이용하여 주어진 문제영역과 허용오차에 따라 자동으로 수치등각사상이 결정되는 알고리즘을 제안한다. 이 알고리즘에서는 지금까지 경험에 의존했었던 표본수와 저주파 필터 파라메터가 주어진 문제영역에 따라 자동으로 결정된다. 이것은 문제의 난이도가 문제영역의 변형에 의존한다는 전제로 문제영역의 모양을 결정하는 함수를 Fourier 급수로 전개, 분석하여 얻을 수 있다. 수치실험을 통해 그 유효성을 입증한다.

• 한국군사과학기술학회지
• /
• 제11권6호
• /
• pp.164-171
• /
• 2008

A modified generalized particle algorithm, MGPA, was suggested to improve the computational efficiency of standard SPH method in numerical analysis of high speed impact behavior. This method uses a numerical failure mechanism than material failure models to describe the target penetration. MGPA algorithm was more effective to describe the impact phenomena and new boundaries produced during the calculation process were well recognized and treated in the target penetration problem of a bullet. When bullet perforation problems were analyzed by this method, MGPA algorithm calculation gives the stable numerical solution and stress oscillation or particle penetration phenomena were not shown. The error range in ballistic velocity limit is less than $2{\sim}13%$ for various target thickness.

• Journal of information and communication convergence engineering
• /
• 제8권3호
• /
• pp.312-316
• /
• 2010

Conformal mapping has been a familiar tool of science and engineering for generations. The methods of numerical mapping are usually classified into those which construct the map from standard domain such as the unit disk onto the 'problem domain', and those which construct the map in the reverse direction. We treat numerical conformal mapping from the unit disk onto the Jordan regions as the problem domain in this paper. The traditional standard methods of this type are based on Theodorsen integral equation. Wegmann's method is well known as a Newton-like efficient one for solving Theodorsen equation. An improved method for convergence by applying low frequency pass filter to the Wegmann's method was proposed. In this paper we propose an effective algorithm for numerical conformal mapping based on the improved method. This algorithm is able to determine the discrete numbers and initial values automatically in accordance with the given region and the required accuracy. This results come from analyzing the shape of given domain as seen in the Fourier Transform.

• 대한전기학회:학술대회논문집
• /
• 대한전기학회 2000년도 하계학술대회 논문집 B
• /
• pp.831-833
• /
• 2000

In this paper, a numerical analysis algorithm of eddy current testing(ECT) for heat exchanger tube with axi-symmetric defects using finite element method(FEM) is presented. In the ECT FEM analysis, we used trianglular and rectangular elements for exact signal of ECT for variable shape of defects. This paper presents a systematic and efficient numerical analysis algorithm for ECT. We employ the LU decomposition and Cholesky method for solving the system matrix. This numerical analysis algorithm is effectively applied to heat exchanger tube with defects.

• Structural Engineering and Mechanics
• /
• 제70권6호
• /
• pp.649-659
• /
• 2019

Vibration-based structural damage detection through optimization algorithms and minimization of objective function has recently become an interesting research topic. Application of various objective functions as well as optimization algorithms may affect damage diagnosis quality. This paper proposes a new damage identification method using Moth-Flame Optimization (MFO). MFO is a nature-inspired algorithm based on moth's ability to navigate in dark. Objective function consists of a term with modal assurance criterion flexibility and natural frequency. To show the performance of the said method, two numerical examples including truss and shear frame have been studied. Furthermore, Los Alamos National Laboratory test structure was used for validation purposes. Finite element model for both experimental and numerical examples was created by MATLAB software to extract modal properties of the structure. Mode shapes and natural frequencies were contaminated with noise in above mentioned numerical examples. In the meantime, one of the classical optimization algorithms called particle swarm optimization was compared with MFO. In short, results obtained from numerical and experimental examples showed that the presented method is efficient in damage identification.

• 대한기계학회논문집B
• /
• 제24권12호
• /
• pp.1555-1569
• /
• 2000

Numerical simulation of fluid flow with moving free surface has been carried out. For the free surface flow, a VOF(Volume of Fluid)-based algorithm utilizing a fixed grid system has been investigated. In order to reduce numerical smearing at the free surface represented on a fixed grid system, a new free surface tracking algorithm based on the donor-acceptor scheme has been presented. Novel features of the proposed algorithm are characterized as two numerical tools; the orientation vector to represent the free surface orientation in each cell and the baby-cell to determine the fluid volume flux at each cell boundary. The proposed algorithm can be easily implemented in any irregular non-uniform grid systems that are usual in finite element method (FEM). Moreover, the proposed algorithm can be extended and applied to the 3-D free surface flow problem without additional efforts. For computation of unsteady incompressible flow, a finite element approximation based on the explicit fractional step method has been adopted. In addition, the SUPG(streamline upwind/Petrov-Galerkin) method has been implemented to deal with convection dominated flows. Combination of the proposed free surface tracking scheme and explicit fractional step formulation resulted in an efficient solution algorithm. Validity of the present solution algorithm was demonstrated from its application to the broken dam and the solitary wave propagation problems.

• 대기
• /
• 제30권4호
• /
• pp.335-346
• /
• 2020

Three free parameters included in a cumulus parameterization are optimized by using micro-genetic algorithm for three precipitation cases occurred in the Korea Peninsula during the summer season in order to reduce biases in a regional model associated with the uncertainties of the parameters and thus to improve the predictability of precipitation. The first parameter is the one that determines the threshold in convective trigger condition. The second parameter is the one that determines boundary layer forcing in convective closure. Finally, the third parameter is the one used in calculating conversion parameter determining the fraction of condensate converted to convective precipitation. Optimized parameters reduce the occurrence of convections by suppressing the trigger of convection. The reduced convection occurrence decreases light precipitation but increases heavy precipitation. The sensitivity experiments are conducted to examine the effects of the optimized parameters on the predictability of precipitation. The predictability of precipitation is the best when the three optimized parameters are applied to the parameterization at the same time. The first parameter most dominantly affects the predictability of precipitation. Short-range forecasts for July 2018 are also conducted to statistically assess the precipitation predictability. It is found that the predictability of precipitation is consistently improved with the optimized parameters.

• Journal of applied mathematics & informatics
• /
• 제26권1_2호
• /
• pp.15-27
• /
• 2008

We present and discuss an algorithm for the numerical solution of initial value problems of the form $D_*^\alpha$y(t) = f(t, y(t)), y(0) = y0, where $D_*^\alpha$y is the derivative of y of order $\alpha$ in the sense of Caputo and 0<${\alpha}{\leq}1$. The algorithm is based on the fractional Euler's method which can be seen as a generalization of the classical Euler's method. Numerical examples are given and the results show that the present algorithm is very effective and convenient.