• Structural Engineering and Mechanics
• /
• v.6 no.3
• /
• pp.339-345
• /
• 1998

The computation size and accuracy in the boundary element method are mutually coupled and strongly influenced by the formulations in boundary discretization and integration. This aspect is studied numerically for two-dimensional elastodynamic problems in the frequency-domain. The localized nature of error is observed in the computed results. A boundary discretization criterion is examined. The number of integration points in the boundary integration is studied to find the optimum number for accuracy. Useful information is obtained concerning the optimization in boundary discretization and integration.

• Structural Engineering and Mechanics
• /
• v.42 no.3
• /
• pp.353-362
• /
• 2012

In this paper decay and mapped elastodynamic infinite elements, based on modified Bessel shape functions and appropriate for Soil-Structure Interaction problems are described and discussed. These elements can be treated as a new form of the recently proposed Elastodynamic Infinite Elements with United Shape Functions (EIEUSF) infinite elements. The formulation of 2D horizontal type infinite elements (HIE) is demonstrated, but by similar techniques 2D vertical (VIE) and 2D corner (CIE) infinite elements can also be formulated. It is demonstrated that the application of the elastodynamical infinite elements is the easier and appropriate way to achieve an adequate simulation including basic aspects of Soil-Structure Interaction. Continuity along the artificial boundary (the line between finite and infinite elements) is discussed as well and the application of the proposed elastodynamical infinite elements in the Finite Element Method is explained in brief. Finally, a numerical example shows the computational efficiency of the proposed infinite elements.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2009.04a
• /
• pp.534-535
• /
• 2009

This paper deals with BEM analysis of transient elastodynamic problems using domain decomposition method and particular integrals. The particular method is used to approximate the acceleration term in the governing equation. The domain decomposition method is examined to consider multi-region problems. The domain of the original problem is subdivided into sub-regions, which are modeled by the particular integral BEM. The iterative coupling employing Schwarz algorithm is used for the successive update of the interface boundary conditions until convergence is achieved. The numerical results, compared with those by ABAQUS, demonstrate the validity of the present formulation.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.5
• /
• pp.883-890
• /
• 1992

Analysis of dynamic or impact problems is very important in engineering fields such as airplanes and automobiles. In the present study, two-dimensional elastodynamic BEM program with Laplace transformation is developed to analyze dynamic or impact problems. Accuracy and efficiency of the BEM program are tested by making the comparision of impact analysis of some models with other's published results. The BEM developed is applied to the impact crack problem and the dynamic stress intensity factors of some impact cracks is obtained by the displacement extrapolation method. It is confirmed to be possible to analyze impact problems accurately with only a little elements in simple models. And also it is found to be careful to use the singular element usually using in static crack problems because that the elastodynamic fundamental solution usually using in static crack problems because that the elastodynamic fundamental solution has more sensitive singularity than the static fundamental solution and to determine the boundary conditions in dynamic problems.

• Journal of the Korean Society for Nondestructive Testing
• /
• v.20 no.2
• /
• pp.130-137
• /
• 2000

Ultrasonic technique which is one of the most common and reliable nondestructive evaluation techniques has been applied to evaluate the integrity of structures by analyzing the characteristics of signal scattered from internal defects. Therefore, the numerical analysis of the ultrasonic scattered field is absolutely necessary for the accurate and quantitative estimation of internal defects. Various modeling techniques now play an important role in nondestructive evaluation and have been employed to solve elastic wave scattering problems. Because the elastodynamic boundary element method is useful to analyze the scattered field in infinite media. it has been used to calculate the ultrasonic wavefields scattered from internal defects. In this study, a review of the boundary element method used for elastic wave scattering problems is presented and, as examples of the boundary element method, the scattered fields due to a circular cavity subjected to incident SH-wave and due to a surface-breaking crack subjected to incident Rayleigh wave are illustrated.

• Structural Engineering and Mechanics
• /
• v.19 no.3
• /
• pp.281-295
• /
• 2005

An analytical method is presented to solve the elastodynamic problem of finitely long hollow cylinder subjected to torsional impact often occurs in engineering mechanics. The analytical solution is composed of a solution of quasi-static equation satisfied with the non-homogeneous boundary condition and a solution of dynamic equation satisfied with homogeneous boundary condition. The quasi-static solution is obtained directly by solving the quasi-static equation satisfied with the non-homogeneous boundary condition. The solution of the non-homogeneous dynamic equation is obtained by means of finite Hankel transform on the radial variable, r, Laplace transform on time variable, t, and finite Fourier transform on axial variable, z. Thus, the solution for finitely long, hollow cylinder subjected to torsion impact is obtained. In the calculating examples, the response histories and distributions of shear stress in the finitely long hollow cylinder subjected to an exponential decay torsion load are obtained, and the results have been analyzed and discussed. Finally, a dynamic finite element for the same problem is carried out by using ABAQUS finite element analysis. Comparing the analytical solution with the finite element solution, it can be found that two kinds of results obtained by means of two different methods agree well. Therefore, it is further concluded that the analytical method and computing process presented in the paper are effective and accurate.

• Coupled systems mechanics
• /
• v.4 no.3
• /
• pp.263-277
• /
• 2015

This paper presents a coupled FEM-BEM strategy for the numerical analysis of elastodynamic problems where infinite-domain models and complex heterogeneous media are involved, rendering a configuration in which neither the Finite Element Method (FEM) nor the Boundary Element Method (BEM) is most appropriate for the numerical analysis. In this case, the coupling of these methodologies is recommended, allowing exploring their respective advantages. Here, frequency domain analyses are focused and an iterative FEM-BEM coupling technique is considered. In this iterative coupling, each sub-domain of the model is solved separately, and the variables at the common interfaces are iteratively updated, until convergence is achieved. A relaxation parameter is introduced into the coupling algorithm and an expression for its optimal value is deduced. The iterative FEM-BEM coupling technique allows independent discretizations to be efficiently employed for both finite and boundary element methods, without any requirement of matching nodes at the common interfaces. In addition, it leads to smaller and better-conditioned systems of equations (different solvers, suitable for each sub-domain, may be employed), which do not need to be treated (inverted, triangularized etc.) at each iterative step, providing an accurate and efficient methodology.

• Structural Engineering and Mechanics
• /
• v.17 no.6
• /
• pp.735-749
• /
• 2004

This work presents a direct integration scheme, based on a fourth order finite difference approach, for elastodynamics. The proposed scheme was chosen as an alternative for attenuating the errors due to the use of the central difference method, mainly when the time-step length approaches the critical time-step. In addition to eliminating the spurious numerical oscillations, the fourth order finite difference scheme keeps the advantages of the central difference method: reduced computer storage and no requirement of factorisation of the effective stiffness matrix in the step-by-step solution. A study concerning the stability of the fourth order finite difference scheme is presented. The Finite Element Method and the Boundary Element Method are employed to solve elastodynamic problems. In order to verify the accuracy of the proposed scheme, two examples are presented and discussed at the end of this work.

• The Journal of the Acoustical Society of Korea
• /
• v.18 no.1E
• /
• pp.37-43
• /
• 1999

This paper deals with a hybrid finite element method for wave scattering problems in infinite domains. Scattering of waves involving complex geometries, in conjunction with infinite domains is modeled by introducing a mathematical boundary within which a finite element representation is employed. On the mathematical boundary, the finite element representation is matched with a known analytical solution in the infinite domain in terms of fields and their derivatives. The derivative continuity is implemented by using a slope constraint. Drilling degrees of freedom at each node of the finite element model are introduced to make the numerical model more sensitive to the transverse component of the elastodynamic field. To verify the effects of drilling degrees freedom and slope constraints individually, reflection of normally incident P and SV waves on a traction free half spaces is considered. For the P-wave incidence, the results indicate that the use of slope constraint is more effective because it suppresses artificial reflection at the mathematical boundary. For the SV-wave case, the use of drilling degrees freedom is more effective by reducing numerical error at irregular frequencies.

• Structural Engineering and Mechanics
• /
• v.67 no.4
• /
• pp.317-326
• /
• 2018

The iterative decomposition coupling formulation of the precise integration finite element method (FEM) and the time domain boundary element method (TD-BEM) is presented for elstodynamic problems. In the formulation, the FEM node and the BEM node are not required to be coincident on the common interface between FEM and BEM sub-domains, therefore, the FEM and BEM are independently discretized. The force and displacement converting matrices are used to transfer data between FEM and BEM nodes on the common interface between the FEM and BEM sub-domains, to renew the nodal variables in the process of the iterations for the un-coincident FEM node and BEM node. The iterative coupling formulation for elastodynamics in current paper is of high modeling accuracy, due to the semi-analytical solution incorporated in the precise integration finite element method. The decomposition coupling formulation for elastodynamics is verified by examples of a cantilever bar under a Heaviside-type force and a harmonic load.