• Structural Engineering and Mechanics
• /
• v.38 no.2
• /
• pp.211-229
• /
• 2011

In this paper, we present a new explicit procedure using periodic cubic B-spline interpolation polynomials to solve linear and nonlinear dynamic equation of motion governing single degree of freedom (SDOF) systems. In the proposed approach, a straightforward formulation was derived from the approximation of displacement with B-spline basis in a fluent manner. In this way, there is no need to use a special pre-starting procedure to commence solving the problem. Actually, this method lies in the case of conditionally stable methods. A simple step-by-step algorithm is implemented and presented to calculate dynamic response of SDOF systems. The validity and effectiveness of the proposed method is demonstrated with four examples. The results were compared with those from the numerical methods such as Duhamel integration, Linear Acceleration and also Exact method. The comparison shows that the proposed method is a fast and simple procedure with trivial computational effort and acceptable accuracy exactly like the Linear Acceleration method. But its power point is that its time consumption is notably less than the Linear Acceleration method especially in the nonlinear analysis.

• Journal of the Korean Society for Precision Engineering
• /
• v.21 no.10
• /
• pp.115-126
• /
• 2004

In the present work a finite element formulation using dynamic-explicit time integration scheme is used for numerical analysis of auto-body panel stamping processes. The lumping scheme is employed for the diagonal mass matrix and dynamic explicit formulation. Analyzed auto-body panel stomping process correction of forming using software called Dynaform using dynamic extensive method. Further, the simulated results for the auto-body panel stamping processes are shown and discussed. Its application is being increased especially in the automotive industrial area for the cost reduction, weight saving, and improvement of strength.

• Structural Engineering and Mechanics
• /
• v.4 no.6
• /
• pp.589-600
• /
• 1996

An explicit 4th order time integration scheme for solving the convection-diffusion equation is discussed in this paper. A system of ordinary differential equations are derived first by discretizing the spatial derivatives of the relevant PDE using the finite difference method. The integration of the ODEs is then carried out using a 4th order scheme and a self-adaptive technique based on the spatial grid spacing. For a non-uniform spatial grid, different time step sizes are used for the integration of the ODEs defined at different spatial points, which improves the computational efficiency significantly. A numerical example is also discussed in the paper to demonstrate the implementation and effectiveness of the method.

• Journal of Industrial Technology
• /
• v.20 no.A
• /
• pp.211-218
• /
• 2000

It is desirable but difficult to predict springback quantitatively and accurately for successful tool and process design in sheet stamping operations. The result of springback analysis by the finite element method (FEM) is sensitively influenced by numerical factors such as blank element size, number of integration points, punch velocity, contact algorithm, etc. In the present work, a parametric study by Taguchi method is performed in order to evaluate the influence of numerical factors on the result of springback analysis quantitatively and to obtain the combination of numerical factors which gives the best approximation to experimental data. Since springback is determined by the residual stress after forming process, it is important to evaluate stress distribution accurately. The oscillation in the time history curve of stress obtained by the dynamic-explicit finite element method says that the stress solution at termination time is in very unstable state. Therefore, a variability study is also carried out in this study in order to assess the stability of implicit springback analysis starting from the stress solution by explicit forming simulation. The U-draw bending process, one of the NUMISHEET '93 benchmark problems, is adopted as an application model because it is most popular one for evaluating the springback characteristic.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2000.10a
• /
• pp.136-140
• /
• 2000

It is desirable to predict springback quantitatively and accurately for the tool and process design in sheet stamping operations, however, it is blown very difficult. The result of springback analysis by the finite element method is sensitively influenced by numerical factors such as blank element size, number of integration point, punch velocity, contact algorithm etc. In the present work, a parametric study by Taguchi method is performed in order to evaluate the influence of numerical factors on springback Quantitatively and to obtain the combination of numerical factors which yields the best approximation to experimental data. Since springback is determined by the residual stress after forming process, it is important to evaluate stress distribution accurately. The oscillation in the time history curve of stress obtained by explicit FEM says that the stress solution at termination time is in very unstable state. Therefore, a variability study is also carried out in this study in order to assess the stability of implicit springback analysis starting from the stress solution by explicit forming simulation. The 2D draw bending process, one of the NUMISHEET '93 benchmark problems, is adopted as an application model.

• Structural Engineering and Mechanics
• /
• v.60 no.1
• /
• pp.149-162
• /
• 2016

Structure-dependent integration methods seem promising for structural dynamics applications since they can integrate unconditional stability and explicit formulation together, which can enable the integration methods to save many computational efforts when compared to an implicit method. A newly developed structure-dependent integration method can inherit such numerical properties. However, an unusual overshooting behavior might be experienced as it is used to compute a forced vibration response. The root cause of this inaccuracy is thoroughly explored herein. In addition, a scheme is proposed to modify this family method to overcome this unusual overshooting behavior. In fact, two improved formulations are proposed by adjusting the difference equations. As a result, it is verified that the two improved formulations of the integration methods can effectively overcome the difficulty arising from the inaccurate integration of the steady-state response of a high frequency mode.

• Structural Engineering and Mechanics
• /
• v.65 no.1
• /
• pp.121-128
• /
• 2018

Three structure-dependent integration methods with no numerical dissipation have been successfully developed for time integration. Although these three integration methods generally have the same numerical properties, such as unconditional stability, second-order accuracy, explicit formulation, no overshoot and no numerical damping, there still exist some different numerical properties. It is found that TLM can only have unconditional stability for linear elastic and stiffness softening systems for zero viscous damping while for nonzero viscous damping it only has unconditional stability for linear elastic systems. Whereas, both CEM and CRM can have unconditional stability for linear elastic and stiffness softening systems for both zero and nonzero viscous damping. However, the most significantly different property among the three integration methods is a weak instability. In fact, both CRM and TLM have a weak instability, which will lead to an adverse overshoot or even a numerical instability in the high frequency responses to nonzero initial conditions. Whereas, CEM possesses no such an adverse weak instability. As a result, the performance of CEM is much better than for CRM and TLM. Notice that a weak instability property of CRM and TLM might severely limit its practical applications.

• Journal of computational fluids engineering
• /
• v.21 no.1
• /
• pp.10-18
• /
• 2016

Numerical boundary conditions are as important as the governing equations when analyzing the fluid flows numerically. An explicit boundary condition method updates the solutions at the boundaries with extrapolation from the interior of the computational domain, while the implicit boundary condition method in conjunction with an implicit time integration method solves the solutions of the entire computational domain including the boundaries simultaneously. The implicit boundary condition method, therefore, is more robust than the explicit boundary condition method. In this paper, steady compressible 2-Dimensional Navier-Stokes solver is developed. We present the implicit boundary condition method coupled with LU-SGS(Lower Upper Symmetric Gauss Seidel) method. Also, the explicit boundary condition method is implemented for comparison. The preconditioning Navier-Stokes equations are solved on unstructured meshes. The numerical computations for a number of flows show that the implicit boundary condition method can give accurate solutions.

• Journal of Ocean Engineering and Technology
• /
• v.12 no.1
• /
• pp.10-22
• /
• 1998

In the present work the elastic-plastic FE formulations using dynamic explicit time integration schemes are used for numerical analysis of a large auto-body panel stamping processes. For analyses of more complex cases with larger and more refined meshes, the explicit method is more time effective than implicit method, and has no convergency problem and has the robust nature of contact and friction algorithms while implicit method is widely used because of excellent accuracy and reliability. The elastic-plastic scheme is more reliable and rigorous while the rigid-plastic scheme require small computation time. In finite element simulation of auto-body panel stamping processes, the roobustness and stability of computation are important requirements since the computation time and convergency become major points of consideration besides the solution accuracy due to the complexity of geometry conditions. The performnce of the dynamic explicit algorithms are investigated by comparing the simulation results of formaing of complicate shaped autobody parts, such as a fuel tank and a rear hinge, with the experimental results. It has been shown that the proposed dynamic explicit elastic-plastic finite element method enables an effective computation for complicated auto-body panel stamping processes.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1996.06a
• /
• pp.97-109
• /
• 1996

In the present work computations are carried out for analysis of complicated sheet metal forming process such as forming of a rear hinge. Finite element formulation using dynamic explicit time integration scheme and step-wise combined Implicit/Explicit scheme are introduced for numerical analysis of sheet metal forming process. The rigid-plastic finite element method based on membrane elements has long been employed as a useful numerical technique for the analysis of sheet metal forming because of its time effectiveness. The explicit scheme in general use is based on the elastic-plastic modelling of material requiring large computation time. In finite element simulation of sheet metal forming processes, the robustness and stability of computation are important requirements since the computation time and convergency become major points of consideration besides the solution accuracy due to the complexity of geometry and boundary conditions. The implicit scheme employs a more reliable and rigorous scheme in considering the equilibrium at each step of deformation, while in the explicit scheme the problem of convergency is eliminated at the cost of solution accuracy. The explicit approach and the implicit approach have merits and demerits, respectively. In order to combine the merits of these two methods a step-wise combined implicit/explicit scheme has been developed.