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

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2009년도 정기 학술대회
• /
• 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.

• 대한기계학회논문집
• /
• 제16권5호
• /
• pp.883-890
• /
• 1992

본 연구에서는 Laplace 변환법에 의한 2차원 동적 문제의 경계요소 프로그램 을 작성하고, 간단한 모델을 해석하여 프로르램의 정도 및 그 유용성을 검토하고, 응 용문제로 동적하중을 받는 균열문제의 몇 가지 모델에 대하여 변위 외삽법으로 균열의 파괴 역학적 파라미터인 동적응력확대계수(dynamic stress intensity factor)를 구하 여 결과를 검토하여 보고자 한다.

• 비파괴검사학회지
• /
• 제20권2호
• /
• pp.130-137
• /
• 2000

대표적인 비파괴 평가 기술들 중의 하나인 초음파응용 기술은 각종 구조물에 존재하는 내부결함에 의한 산란신호를 통해 건전성을 평가하는 기법이므로 결함의 신뢰성 높은 정량적 평가를 위해서는 결함으로부터의 초음파 산란신호특성에 대한 기본적 이해가 필수적이며 따라서 이를 위한 모델링 수치해석 연구가 요구된다. 모델링 기법들은 비파괴 평가기술에 있어서 중요한 역할을 하고 있고 많은 모델링 기법들이 결함의 탄성파 산란문제를 해석하기 위하여 사용되어 오고 있다. 본 연구에서는 다양한 수치기법들 중 탄성파 산란문제 해석에 있어 많은 장점을 가지고 있는 동탄성 경계요소법에 대하여 자세히 소개하고, 응용 예로서 경계요소법을 이용한 기공결함의 수평횡파 산란장 해석과 표면균열의 표면파 산란장 해석을 소개한다.

• Structural Engineering and Mechanics
• /
• 제19권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
• /
• 제4권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
• /
• 제17권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
• /
• 제18권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
• /
• 제67권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.