• Structural Engineering and Mechanics
• /
• 제47권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
• /
• 제7권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
• /
• 제36권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.

• 한국재난정보학회 논문집
• /
• 제10권1호
• /
• pp.123-140
• /
• 2014

회전쉘의 정적문제 및 동적문제에 대해서 선점법에 의한 공간영역의 이산화에 대한 유효성을 비교 검토하였고, 정식화에 있어서 이산화 후의 연립방정식 및 시간에 관한 연립상미분방정식의 계수행렬에 대한 영향을 검토하였다. 수치해석법에는 동적특성의 문제를 간단하고 효과적인 방법을 위하여 선점법을 제안하였으며, 수치해석예로서 평균2승잔차법과의 비교검토로 고정도의 해석을 위한 정확성을 검증하였다.

• 한국기계가공학회지
• /
• 제21권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.

• 한국전산구조공학회논문집
• /
• 제27권4호
• /
• pp.297-303
• /
• 2014

페리다이나믹스 이론과 이진분해 기법의 병렬연산을 이용하여 동적 균열진전 문제에 대한 애조인 형상 설계민감도 해석법을 개발하였다. 페리다이나믹스에서는 균열의 연속적인 분기를 다룰 수 있으며, Explicit 시간적분법을 채택한다. 설계민감도 해석은 애조인 변수법은 경로의존성 문제에는 적합하지 않으나 여기서는 응답해석의 경로를 이미 알고 있으므로 채택하여 사용할 수 있었다. 얻어진 해석적 설계민감도는 유한차분과 비교하여 그 정확성을 검증하였다. 유한차분법은 설계섭동량에 민감하여 비선형성이 강한 페리다이나믹스 문제에서 부정확한 설계민감도를 제시할 수 있다. 정확한 설계민감도 해석을 위해서는 이산화과정에서 $C^1$ 연속성을 가지는 체적율이 필요함을 알 수 있었다.

• 한국해양공학회지
• /
• 제33권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.

• 대한토목학회논문집
• /
• 제10권4호
• /
• pp.141-150
• /
• 1990

하구나 연안에서 해수의 순환형상을 모사(模寫)하게 위해 천수방정식(淺水方程式)을 여러 가지 경계조건 하에서 수치해석하였다. 공간영역은 Galerkin방법으로 이산화(離散化)하였으며 시간영역에 대해서는 유한차분법(Crank-Nicolson방법)을 사용하였다. 네 가지 검정실험이 해석적인 해가 있는 일차원 수로에서 행하여졌으며, 해석해를 구할 수 없는 이차원 모형에도 적용되었다. 해석해가 있는 경우 수치모사 결과가 이와 잘 일치하였으며, 이차원 모형에서의 결과도 매우 합당함을 알 수 있었다. 또 일차원 문제에서 4점 bilinear요소와 삼각형 요소를 사용한 결과를 각각 비교하였으며 시간적분도 2단계 Lax-Wendroff방법을 사용하여 결과를 비교하였다. 음해법을 사용할 경우 비교적 정확한 결과를 얻을 수 있으나 요소의 갯수가 많아지면 구성되는 대수방정식(代數方程式)이 커지기 때문에 각 시간마다의 계산량이 엄청나게 늘어나게 되며 양해법을 사용할 때는 원하는 만큼의 정확한 결과를 얻기 위하여 시간간격이나 공간격자 간격을 선정하는데 각별히 유의하여야 할 것이다.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2005년도 추계 학술대회논문집
• /
• 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.

• 한국전산유체공학회지
• /
• 제10권4호통권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.