• Structural Engineering and Mechanics
• /
• v.14 no.1
• /
• pp.71-84
• /
• 2002

A finite element formulation for the time-domain analysis of linear transient elastodynamic problems is presented based on the weak form obtained by applying the Galerkin's method to the integro-differential equations which contain the initial conditions implicitly and does not include the inertia terms. The weak form is extended temporally under the assumptions of the constant and linear time variations of field variables, since the time-stepping algorithms such as the Newmark method and the Wilson ${\theta}$-method are not necessary, obtaining two kinds of implicit finite element equations which are tested for numerical accuracy and convergency. Three classical examples having finite and infinite domains are solved and numerical results are compared with the other analytical and numerical solutions to show the versatility and accuracy of the presented formulation.

• Structural Engineering and Mechanics
• /
• v.3 no.2
• /
• pp.145-162
• /
• 1995

A damped trapezoidal rule method for numerical time-integration is presented, and its application in analyses of dynamic response of damped structures is discussed. It is shown that the damped trapezoidal rule method has features that make it an attractive approach for applications in dynamic analyses of structures. Accuracy and stability analyses are developed for the damped single-degree-of-freedom systems. Error analyses are also performed for the Newmark beta method and compared with the damped trapezoidal rule method as a basis for discussion of the relative merits of the proposed method. The procedure is fully explicit and easy to implement. However, since the method is an explicit method, it is conditionally stable. The methodology is applied to several example problems to illustrate its strengths, limitations and inherent simplicity.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 2003.06a
• /
• pp.1545-1548
• /
• 2003

In rigid-plastic finite element method, there is a heavy computation time and convergence problem. In this study. static-explicit rigid-plastic finite element method will be introduced. This method is the way that restrict the convergence interval. In result, convergence problem and computation time due to large non-linearity in the existing numerical analysis method were no longer a critical problem. Also, we investigated the effect of punch stroke and the strain increment this method. It is expected that various results from the numerical analysis will give very useful information for the design of tools in sheet metal forming process.

• Journal of the Institute of Electronics Engineers of Korea TC
• /
• v.39 no.10
• /
• pp.38-45
• /
• 2002

This paper presents the new Compact 2D ADI-FDTD(Alternating-Direction Implicit Finite-Difference Time-Domain) method, where the time step is no longer restricted by the numerical stability condition. This method is an accelerating algorithm for the conventional Compact 2D FDTD method. To validate this algorithm, we have analyzed the dispersion characteristics of the hollow rectangular waveguide and the shielded microstrip line. The results of the proposed method are very well agreed with those of both the conventional analytic method and the Compact 2D FDTD method. The CPU time for analysis of this method is very much reduced compared with the conventional Compact 2D FDTD method. The proposed method is valuable as a fast algorithm in the research of dispersion characteristics of the waveguide structures.

• Journal of applied mathematics & informatics
• /
• v.33 no.5_6
• /
• pp.485-502
• /
• 2015

A computational method for solving singularly perturbed boundary value problem of differential equation with shift arguments of mixed type is presented. When shift arguments are sufficiently small (o(ε)), most of the existing method in the literature used Taylor's expansion to approximate the shift term. This procedure may lead to a bad approximation when the delay argument is of O(ε). The main idea for this work is to deal with constant shift arguments, which are independent of ε. In the present method, we construct the formally asymptotic solution of the problem using the method of composite expansion. The reduced problem is solved numerically by using operator compact implicit method, and the second problem is solved analytically. Error estimate is derived by using the maximum norm. Numerical examples are provided to support the theoretical results and to show the efficiency of the proposed method.

• Structural Engineering and Mechanics
• /
• v.85 no.3
• /
• pp.369-379
• /
• 2023

Reliability-Based Design Optimization (RBDO) is an appropriate framework for obtaining optimal designs by taking uncertainties into account. Large-scale problems with implicit limit state functions and problems with discrete design variables are two significant challenges to traditional RBDO methods. To overcome these challenges, this paper proposes a hybrid method to perform RBDO of structures that links Firefly Algorithm (FA) as an optimization tool to advanced (finite element) reliability methods. Furthermore, the Genetic Algorithm (GA) and the FA are compared based on the design cost (objective function) they achieve. In the proposed method, Weighted Simulation Method (WSM) is utilized to assess reliability constraints in the RBDO problems with explicit limit state functions. WSM is selected to reduce computational costs. To performing RBDO of structures with finite element modeling and implicit limit state functions, a First-Order Reliability Method (FORM) based on the Direct Differentiation Method (DDM) is utilized. Four numerical examples are considered to assess the effectiveness of the proposed method. The findings illustrate that the proposed RBDO method is applicable and efficient for RBDO problems with discrete and continuous design variables and finite element modeling.

• Journal of the Korean Society of Manufacturing Process Engineers
• /
• v.3 no.3
• /
• pp.39-45
• /
• 2004

In rigid-plastic finite element method, there is a heavy computation time and convergence problem. In this study, static-explicit rigid-plastic finite element method will be introduced. This method is the way that restrict the convergence interval. In result, convergence problem and computation time due to large non-linearity in the existing numerical analysis method were no longer a critical problem. It is expected that various results from the numerical analysis will give very useful information for the design of tools in sheet metal forming process.

• Journal of the Korean Society of Manufacturing Process Engineers
• /
• v.2 no.2
• /
• pp.91-97
• /
• 2003

In rigid-plastic finite element method, there is a heavy computation time and convergence problem. In this study, revised rigid-plastic finite element method Will be introduced. This method is the way that restrict the convergence interval. In result, convergence problem and computation time due to large non-linearity in the existing numerical analysis method were no longer a critical problem. It is expected that various results from the numerical analysis will give very useful information for the design of tools in sheet metal forming process.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.6
• /
• pp.675-682
• /
• 2010

The flood impact pressure acting on a vertical wall resulting from a dam-breaking problem is simulated using a navier-Stokes(N-S) solver. The N-S solver uses Eulerian Finite Volume Method(FVM) along with Volume Of Fluid(VOF) method for 2-D incompressible free surface flows. A Split Lagrangian Advection(SLA) scheme for VOF method is implemented in this paper. The SLA scheme is developed based on an algorithm of Piecewise Linear Interface Calculation(PLIC). The coupling between the continuity and momentum equations is affected by using a well-known Semi-Implicit Method for Pressure-Linked Equations (SIMPLE) algorithm. Several two-dimensional numerical simulations of the dam-breaking problem are presented to validate the accuracy and demonstrate the capability of the present algorithm. The significance of the time step and grid resolution are also discussed. The computational results are compared with experimental data and with computations by other numerical methods. The results showed a favorable agreement of water impact pressure as well as the global fluid motion.

• International Union of Geodesy and Geophysics Korean Journal of Geophysical Research
• /
• v.26 no.1
• /
• pp.1-14
• /
• 1998

Various numerical methods for the two dimensional shallow water equations have been applied to the problems of flood routing, tidal circulation, storm surges, and atmospheric circulation. These methods are often based on the Alternating Direction Implicity(ADI) method. However, the ADI method results in inaccuracies for large time steps when dealing with a complex geometry or bathymetry. Since this method reduces the performance considerably, a fully implicit method developed by Wilders et al. (1998) is used to improve the accuracy for a large time step. Finite Difference Methods are defined on a rectangular grid. Two drawbacks of this type of grid are that grid refinement is not possibile locally and that the physical boundary is sometimes poorly represented by the numerical model boundary. Because of the second deficiency several purely numerical boundary effects can be involved. A boundary fitted curvilinear coordinate transformation is used to reduce these difficulties. It the curvilinear coordinate transformation is used to reduce these difficulties. If the coordinate transformation is orthogonal then the transformed shallow water equations are similar to the original equations. Therefore, an orthogonal coorinate transformation is used for defining coordinate system. A multigrid (MG) method is widely used to accelerate the convergence in the numerical methods. In this study, a technique using a MG method is proposed to reduce the computing time and to improve the accuracy for the orthogonal to reduce the computing time and to improve the accuracy for the orthogonal grid generation and the solutions of the shallow water equations.