• Computational Structural Engineering
• /
• v.9 no.3
• /
• pp.187-197
• /
• 1996

An explicit direct time integration method based solution algorithm is presented to predict dynamic buckling response of Euler column. On the basis of large deflection beam theory, a plane frame finite element is formulated and implemented into the solution algorithm. The element formulation takes into account geometrical nonlinearity and overall buckling of steel structural frames. The solution algorithm employs the central difference method. Using the computer program developed by the author, dynamic instability behavior of Euler column under impact loading is investigated by considering the time variation of load, load magnitude, and load duration. The free vibration of Euler column caused by a short duration impact load is also studied. The validity and efficiency of the present formulation and solution algorithm are verified through illustrative numerical examples.

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

• Transactions of the Korean Society of Automotive Engineers
• /
• v.4 no.5
• /
• pp.16-25
• /
• 1996

In the present work a rigid-plastic finite element formulation using dynamic explicit time integration scheme is proposed for numerical analysis of auto-body panel stamping processes. 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. A damping scheme is proposed in order to achieve a stable solution procedure in dynamic sheet forming problems. In order to improve the drawbacks of the conventional membrane elements, BEAM(abbreviated from Bending Energy Augmented Membrane) elements are employed. Rotational damping and spring about the drilling direction are introduced to prevent a zero energy mode. The lumping scheme is employed for the diagonal mass matrix and linearizing dynamic formulation. A contact scheme is developed by combining the skew boundary condition and the direct trial-and-error method. Computations are carried out for analysis of complicated auto-body panel stamping processes such as forming of an oilpan, a fuel tank and a front fender. The numerical results of explicit analysis are compared with the implicit results with good agreements and it is shown that the explicit scheme requires much shorter computational time, especially when the problem becomes more complicated. It is thus shown that the proposed dynamic explicit rigid-plastic finite element method enables an effective computation for complicated autobody panel stamping processes.

• Journal of the Korea institute for structural maintenance and inspection
• /
• v.6 no.3
• /
• pp.139-149
• /
• 2002

An explicit direct time integration method based solution algorithm is proposed to predict dynamic postbuckling response of thin plates. Based on the von Karman's plate equations and Marquerre's shallow shell theory, a rectangular plate finite element is formulated and utilized in this study. The element formulation takes into account geometrical nonlinearity and initial deflection of plates. The solution algorithm employs the central difference method. Using the computer program developed by the authors, dynamic postbuckling behavior of elastic thin plates under impact loading is investigated by considering the time variation of load and load duration. The efficiency of the proposed solution algorithm is examined through illustrative numerical examples.

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

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

• Earthquakes and Structures
• /
• v.7 no.2
• /
• pp.119-139
• /
• 2014

A new method is proposed for random vibration anaylsis of hysteretic systems subjected to non-stationary random excitations. With the Bouc-Wen model, motion equations of hysteretic systems are first transformed into quasi-linear equations by applying the concept of equivalent excitations and decoupling of the real and hysteretic displacements, and the derived equation system can be solved by either the precise time integration or the Newmark-${\beta}$ integration method. Combining the numerical solution of the auxiliary differential equation for hysteretic displacements, an explicit iteration algorithm is then developed for the dynamic response analysis of hysteretic systems. Because the computational cost for a large number of deterministic analyses of hysteretic systems can be significantly reduced, Monte-Carlo simulation using the explicit iteration algorithm is now viable, and statistical characteristics of the non-stationary random responses of a hysteretic system can be obtained. Numerical examples are presented to show the accuracy and efficiency of the present approach.

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

• Structural Engineering and Mechanics
• /
• v.22 no.1
• /
• pp.107-130
• /
• 2006

The paper deals with the numerical solution of the dynamic problem of masonry structures. Masonry is modelled as a non-linear elastic material with zero tensile strength and infinite compressive strength. Due to the non-linearity of the adopted constitutive equation, the equations of the motion must be integrated directly. In particular, we apply the Newmark or the Hilber-Hughes-Taylor methods implemented in code NOSA to perform the time integration of the system of ordinary differential equations obtained from discretising the structure into finite elements. Moreover, with the aim of evaluating the effectiveness of these two methods, some dynamic problems, whose explicit solutions are known, have been solved numerically. Comparisons between the exact solutions and the corresponding approximate solutions obtained via the Newmark and Hilber-Hughes-Taylor methods show that in the cases under consideration both numerical methods yield satisfactory results.

• Journal of Ocean Engineering and Technology
• /
• v.14 no.2
• /
• pp.28-35
• /
• 2000

In this paper a complicated structural behavior in collision and its effect of energy transmission to the collision bulkhead was examined through a methodology of the numerical simulation to obtain a ideal bow construction and a location of collision bulkhead against heat on collision. At present the bow structure is normally designed in consideration of its specific structural arrangement and internal and external loads in these areas such as hydrostatic and dynamic pressure wave impact and bottom slamming in accordance with the Classification rules and the specific location of collision bulkhead by SOLAS requirement. By these studies the behavior of the bow collapse due to collision was synthetically evaluated for the different size of tankers and its operational speed limits and by the results of these simulation it provides the optimal design concept for the bow construction to prevent the subsequent plastic deformation onto or near to the collision bulkhead boundary and to determine the rational location of collision bulkhead.