• Structural Engineering and Mechanics
• /
• v.55 no.5
• /
• pp.1055-1086
• /
• 2015

The dynamic response of two-dimensional unbounded domain on the rigid bedrock in the time domain is numerically obtained. It is realized by the modified scaled boundary finite element method (SBFEM) in which the original scaling center is replaced by a scaling line. The formulation bases on expanding dynamic stiffness by using the continued fraction approach. The solution converges rapidly over the whole time range along with the order of the continued fraction increases. In addition, the method is suitable for large scale systems. The numerical method is employed which is a combination of the time domain SBFEM for far field and the finite element method used for near field. By using the continued fraction solution and introducing auxiliary variables, the equation of motion of unbounded domain is built. Applying the spectral shifting technique, the virtual modes of motion equation are eliminated. Standard procedure in structural dynamic is directly applicable for time domain problem. Since the coefficient matrixes of equation are banded and symmetric, the equation can be solved efficiently by using the direct time domain integration method. Numerical examples demonstrate the increased robustness, accuracy and superiority of the proposed method. The suitability of proposed method for time domain simulations of complex systems is also demonstrated.

• Structural Engineering and Mechanics
• /
• v.21 no.4
• /
• pp.437-450
• /
• 2005

In this study, a complete analysis of soil-structure interaction problems is presented which includes a modelling of the near surrounding of the building (near-field) and a special description of the wave propagation process in larger distances (far-field). In order to reduce the computational effort which can be very high for time domain analysis of wave propagation problems, a special approach based on similarity transformation of the infinite domain on the near-field/far-field interface is applied for the wave radiation of the far-field. The near-field is discretised with standard Finite Elements, which also allows to introduce non-linear material behaviour. In this paper, a new approach to calculate the involved convolution integrals is presented. This approximation in time leads to a dramatically reduced computational effort for long simulation times, while the accuracy of the method is not affected. Finally, some benchmark examples are presented, which are compared to a coupled Finite Element/Boundary Element approach. The results are in excellent agreement with those of the coupled Finite Element/Boundary Element procedure, while the accuracy is not reduced. Furthermore, the presented approach is easy to incorporate in any Finite Element code, so the practical relevance is high.

• Journal of The Korean Society of Agricultural Engineers
• /
• v.48 no.6
• /
• pp.45-56
• /
• 2006

The finite element method is a practical tool to compute the response of the irregularly layered soil deposit to the base-rock motions. The method is useful not only in estimating the interaction between the structure and the surrounding soil as a whole and the local behavior of the contacting area in detail, but also in predicting the resulting behavior of the superstructure affected by such soil-structure interactions. However, the computation of finite element analysis is marched in the time domain (TD), while the site response analysis has been carried out mostly in the frequency domain (FD) with equivalent linear analysis. This study is intended to compare the results of the TD and FD analysis with focus on the peak response accelerations and the predominant frequencies, and thus to evaluate the applicability and the validity of the finite element analysis in the site response analysis. The comparison shows that one can obtain the results very close to that of FD analysis, from the finite element analysis by including sufficiently large width of foundation in the model and further by applying partial mode superposition. The finite element analysis turned out to be well agreeing with FD analysis in their computed results of the peak acceleration and the acceleration response spectra, especially at the surface layer.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.25 no.4
• /
• pp.574-581
• /
• 2001

Transient linear elastodynamic problems are numerically analyzed in a time-domain by the Finite Element Method, for which the variational formulation based upon the equations of motion in convolution integral is newly derived. This formulation is implicit and does not include the time derivative terms so that the computation procedure is simple and less assumptions are required comparing to the conventional time-domain dynamic numerical algorithms, being able to get the improved numerical accuracy and stability. That formulation is expanded using the semi-discrete approximation to obtain the finite element equations. In the temporal approximation, the time axis is divided equally and constant and linear time variations are assumed in those intervals. It is found that unconditionally stable numerical results are obtained in case of the constant time variation. Some numerical examples are given to show the versatility of the presented formulation.

• Earthquakes and Structures
• /
• v.9 no.1
• /
• pp.173-194
• /
• 2015

Site response analysis is an important topic in earthquake engineering. A time-domain numerical method called as one-dimensional (1D) finite element artificial boundary method is proposed to simulate the homogeneous plane elastic wave propagation in a layered half space subjected to the obliquely incident plane body wave. In this method, an exact artificial boundary condition combining the absorbing boundary condition with the inputting boundary condition is developed to model the wave absorption and input effects of the truncated half space under layer system. The spatially two-dimensional (2D) problem consisting of the layer system with the artificial boundary condition is transformed equivalently into a 1D one along the vertical direction according to Snell's law. The resulting 1D problem is solved by the finite element method with a new explicit time integration algorithm. The 1D finite element artificial boundary method is verified by analyzing two engineering sites in time domain and by comparing with the frequency-domain transfer matrix method with fast Fourier transform.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2005.11b
• /
• pp.85-89
• /
• 2005

For the analysis of a vibrating two dimensional structure such as the simply supported rectangular plate, Spectral Finite Element Method (SFEM) has been studied. Under the condition that two parallel edges are simply supported at least and the other two edges can be arbitrary, Spectral Finite Element has been developed. Using this element SFEM is applied to the vibrating rectangular plate which all edges are simply supported, and obtain the frequency response function in frequency domain and the dynamic response in time domain. To evaluate these results normal mode method and finite element method (FEM) are also accomplished and compared. It is seen that SFEM is more powerful analysis tool than FEM in high frequency range.

• Bulletin of the Korean Mathematical Society
• /
• v.35 no.2
• /
• pp.345-362
• /
• 1998

We introduce and anlyze a naturally parallelizable frequency-domain method for parabolic problems with a Neumann boundary condition. After taking the Fourier transformation of given equations in the space-time domain into the space-frequency domain, we solve an indefinite, complex elliptic problem for each frequency. Fourier inversion will then recover the solution of the original problem in the space-time domain. Existence and uniqueness of a solution of the transformed problem corresponding to each frequency is established. Fourier invertibility of the solution in the frequency-domain is also examined. Error estimates for a finite element approximation to solutions fo transformed problems and full error estimates for solving the given problem using a discrete Fourier inverse transform are given.

• Ocean Systems Engineering
• /
• v.1 no.4
• /
• pp.317-336
• /
• 2011

A time domain finite element based method is employed to analyze wave radiation by multiple cylinders. The nonlinear free surface and body surface boundary conditions are satisfied based on the perturbation method up to the second order. The first- and second-order velocity potential problems at each time step are solved through a finite element method (FEM). The matrix equation of the FEM is solved through an iteration and the initial solution is obtained from the result at the previous time step. The three-dimensional (3D) mesh required is generated based on a two-dimensional (2D) hybrid mesh on a horizontal plane and its extension in the vertical direction. The hybrid mesh is generated by combining an unstructured grid away from cylinders and two structured grids near the cylinder and the artificial boundary, respectively. The fluid velocity on the free surface and the cylinder surface are calculated by using a differential method. Results for various configurations including two-cylinder and four-cylinder cases are provided to show the mutual influence due to cylinders on the first and second waves and forces.

• Journal of the Korean Society of Safety
• /
• v.8 no.4
• /
• pp.201-206
• /
• 1993

A space-time finite element method was presented for time dependent problem. The method which treat both the space and time unformly were proposed and numerically tested. The weighted residual process was used to formulate a finite element method in a space-time domain based upon continuous Galerkin method. This method leads to a conditional stabie high-order accurate solver.

• Journal of Ocean Engineering and Technology
• /
• v.23 no.1
• /
• pp.16-23
• /
• 2009

The three-dimensional ship motion with forward speed was solved by a finite element method in the time domain. A boundary value problem was described in the frame of a fixed-body reference, and the problem was formulated according to Double-Body and Neumann-Kelvin linearizations. Laplace's equation with boundary conditions was solved by a classical finite element method based on the weak formulation. Chebyshev filtering was used to get rid of an unwanted saw-tooth wave and a wave damping zone was adopted to impose a numerical radiation condition. The time marching of the free surface was performed by the 4th order Adams-Bashforth-Moulton method. Wigley I and Wigely III models were considered for numerical validation. The hydrodynamic coefficients and wave exciting forces were validated by a comparison with experimental data and the numerical results of the Wigley I. The effects of the linearization are also discussed. The motion RAO was also checked with a Wigley III model through mono-chromatic and multi-chromatic regular waves.