• Structural Engineering and Mechanics
• /
• v.47 no.5
• /
• pp.599-622
• /
• 2013

Seismic response of two dimensional liquid tanks is numerically simulated using fully nonlinear velocity potential theory. Galerkin-weighted-residual based finite element method is used for solving the governing Laplace equation with fully nonlinear free surface boundary conditions and also for velocity recovery. Based on mixed Eulerian-Lagrangian (MEL) method, fourth order explicit Runge-Kutta scheme is used for time integration of free surface boundary conditions. A cubic-spline fitted regridding technique is used at every time step to eliminate possible numerical instabilities on account of Lagrangian node induced mesh distortion. An artificial surface damping term is used which mimics the viscosity induced damping and brings in numerical stability. Four earthquake motions have been suitably selected to study the effect of frequency content on the dynamic response of tank-liquid system. The nonlinear seismic response vis-a-vis linear response of rectangular liquid tank has been studied. The impulsive and convective components of hydrodynamic forces, e.g., base shear, overturning base moment and pressure distribution on tank-wall are quantified. It is observed that the convective response of tank-liquid system is very much sensitive to the frequency content of the ground motion. Such sensitivity is more pronounced in shallow tanks.

• Geomechanics and Engineering
• /
• v.7 no.6
• /
• pp.613-631
• /
• 2014

To complement the method of field-scale seismic ground motion simulations by buried blast techniques, the application and evaluation of the capability of a numerical modeling platform to simulate buried explosion-induced ground motion at a real soil site is presented in this paper. Upon a layout of the experimental setup at a level site wherein multiple charges that were buried over a large-diameter circle and detonated in a planned sequence, the formulation of a numerical model of the soil and the explosives using the finite element code LS-DYNA is developed for the evaluation of the resulting ground motion and surface subsidence. With a compact elastoplastic cap model calibrated for the loess soils on the basis of the site and laboratory test program, numerical solutions are obtained by explicit time integration for various dynamic aspects and their relation with the field blast experiment. Quantitative comparison of the computed ground acceleration time histories at different locations and induced spatial subsidence on the surface afterwards is given for further engineering insights in regard to the capabilities and limitations of both the numerical and experimental approaches.

• Structural Engineering and Mechanics
• /
• v.36 no.3
• /
• pp.301-319
• /
• 2010

The paper seeks to explore some aspects of the current state of knowledge on progressive collapse in the technical literature covering blast loads and structural analysis procedure applicable to reinforced concrete (RC) buildings. The paper describes the progressive collapse analysis of a commercial RC building located in the city of Riyadh and subjected to different blast scenarios. A 3-D finite element model of the structure was created using LS-DYNA, which uses explicit time integration algorithms for solution. Blast loads were treated as dynamic pressure-time history curves applied to the exterior elements. The inherent shortcomings of notional member removal have been taken care of in the present paper by simulating the damage of structural elements through the use of solid elements with the provision of element erosion. Effects of erosion and cratering are studied for different scenarios of the blast.

• Journal of the Society of Disaster Information
• /
• v.10 no.1
• /
• pp.123-140
• /
• 2014

This paper proposes a simple and effective method for determining the dynamic characteristics of revolution shells. This is a weighted residual method in which the collocation points are taken at the roots of orthogonal polynomial. In this paper the collocation method is employed to replace a partical differential eqations by a system of ordinary differential equations in time, and the resulting equations are solved by two different numerical methods of time integration : an implicit method and an explicit method. The proposed approach is formulated in some detail. The versatility and accuracy are illustrated through several numerical examples. The method appears to be relatively easy to set up and gives satisfactory results.

• Journal of the Korean Society of Manufacturing Process Engineers
• /
• v.21 no.8
• /
• pp.93-98
• /
• 2022

Adiabatic blanking is a method to improve productivity through an autocatalytic cycle that occurs repeatedly through plastic deformation and thermal softening caused by impact energy. In this study, an axisymmetric analysis model comprising a punch, die, holder, and specimen was developed to confirm the temperature and deformation characteristics caused by an impact load. Through this, the impact energy, diameter of the punch, gap between the punch and die, and the effect of the fillet were analyzed. Because this process occurs in a very short time, adiabatic analysis can be performed using the explicit time-integration method. The analysis, confirmed that it is necessary to design a structure capable of increasing the local temperature and plastic deformation by controlling the impact energy, working area, gap, and the fillet.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.4
• /
• pp.297-303
• /
• 2014

Using the bond-based peridynamics and the parallel computation with binary decomposition, an adjoint shape design sensitivity analysis(DSA) method is developed for the dynamic crack propagation problems. The peridynamics includes the successive branching of cracks and employs the explicit scheme of time integration. The adjoint variable method is generally not suitable for path-dependent problems but employed since the path of response analysis is readily available. The accuracy of analytical design sensitivity is verified by comparing it with the finite difference one. The finite difference method is susceptible to the amount of design perturbations and could result in inaccurate design sensitivity for highly nonlinear peridynamics problems with respect to the design. It turns out that $C^1$-continuous volume fraction is necessary for the accurate evaluation of shape design sensitivity in peridynamic discretization.

• Journal of Ocean Engineering and Technology
• /
• v.33 no.5
• /
• pp.400-410
• /
• 2019

In this study, the numerical code for the 3D nonlinear dynamic analysis of an SLWR (Steel Lazy Wave Riser) was developed using the lumped mass line model in a FORTRAN environment. Because the lumped mass line model is an explicit method, there is no matrix operation. Thus, the numerical algorithm is simple and fast. In the lumped mass line model, the equations of motion for the riser were derived by applying the various forces acting on each node of the line. The applied forces at the node of the riser consisted of the tension, shear force due to the bending moment, gravitational force, buoyancy force, riser/ground contact force, and hydrodynamic force based on the Morison equation. Time integration was carried out using a Runge-Kutta fourth-order method, which is known to be stable and accurate. To validate the accuracy of the developed numerical code, simulations using the commercial software OrcaFlex were carried out simultaneously and compared with the results of the developed numerical code. To understand the nonlinear dynamic characteristics of an SLWR, dynamic simulations of SLWRs excited at the hang-off point and of SLWRs in regular waves were carried out. From the results of these dynamic simulations, the displacements at the maximum bending moments at important points of the design, like the hang-off point, sagging point, hogging points, and touch-down point, were observed and analyzed.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.10 no.4
• /
• pp.141-150
• /
• 1990

Numerical experiments of sballow water equations are performed under various boundary conditions by finite element method to simulate the circulation in estuaries and coastal areas. Galerkin method is employed to discretize spatial domain, and for time integration, finite difference method (Crank-Nicolson scheme) is used. This method is tested in five problems, in which first four cases have analytic solutions. The computed values are well in agreement with the analytic solutions in four experiments and the result of the last 2-dimensional ease is resonable. Implicit and two step Lax-Wendroff schemes in time domain are compared, and the results when using four node bilinear and triangular elements are presented. Consequently it takes very long time for complex problems requiring many elements to integrate all the time steps using the implicit schemes. And the explicit scheme requires careful consideration in selecting the time step and the grid size to obtain the desired accuracy.

• 한국전산유체공학회:학술대회논문집
• /
• 2005.10a
• /
• pp.235-238
• /
• 2005

The unsteady supersonic flow over two- and three-Dimensional cavities has been analyzed by the integration of unsteady Reynolds-Averaged Navier-Stokes(RANS) with the k - w turbulence model. The unsteady flow is characterized by the periodicity due to the mutual relation between the shear layer and the internal flow in cavities. Numerical method is upwind TVD scheme based on the flux vector split with the Van Leer limiters, and time accuracy is used explicit 4th stage Runge-Kutta scheme. Cavity flows are Comparison of two- and three-dimensional. The cavity has a L/D ratio of 3 for two-dimensional case. and same L/D and W/D ratio is 1 for three-dimensional case. The Mach and Reynolds numbers are held constant at 1.5 and 450000 respectively. For the three-dimensional case, the flow field is observed to oscillate in the 'shear layer mode' with a feedback mechanism that follow Rossiter's formula. On the other hand, the self-sustained oscillating flow transitions to a 'wake mode' for the two-dimensional simulation, with more violent fluctuations inside the cavity.

• Journal of computational fluids engineering
• /
• v.10 no.4 s.31
• /
• pp.51-58
• /
• 2005

The unsteady supersonic flow over two- and three-dimensional cavities has been analyzed by the integration of unsteady Reynolds-Averaged Navier-Stokes(RANS) with the k-$\omega$ turbulence model. The unsteady flow is characterized by the periodicity due to the mutual relation between the shear layer and the internal flow in the cavity. An explicit 4th order Runge-Kutta scheme and an upwind TVD scheme based on the flux vector split with the van Leer limiters are used for time and space discritizations, respectively. The cavity has a L/D ratio of 3 for two-dimensional case, and same L/D and W/D ratio of I for three-dimensional case. The Mach and Reynolds numbers are 1.5 and 450000 respectively. In the three-dimensional flow, the field is observed to oscillate in the 'shear layer mode' with a feedback mechanism that follows Rossiter's formula. In the two-dimensional simulation, the self-sustained oscillating flow has more violent fluctuation inside the cavity. The primary fluctuating frequencies of two- and three- dimensional flow agree very well with the 2nd mode of Rossiter's frequency. In the three-dimensional flow, the 1st mode of frequency could be seen.