• 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.

• Structural Engineering and Mechanics
• /
• v.14 no.1
• /
• pp.71-84
• /
• 2002

A finite element formulation for the time-domain analysis of linear transient elastodynamic problems is presented based on the weak form obtained by applying the Galerkin's method to the integro-differential equations which contain the initial conditions implicitly and does not include the inertia terms. The weak form is extended temporally under the assumptions of the constant and linear time variations of field variables, since the time-stepping algorithms such as the Newmark method and the Wilson ${\theta}$-method are not necessary, obtaining two kinds of implicit finite element equations which are tested for numerical accuracy and convergency. Three classical examples having finite and infinite domains are solved and numerical results are compared with the other analytical and numerical solutions to show the versatility and accuracy of the presented 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.102-109
• /
• 2000

This paper describes two different theories used to model the scattering of ultrasound by a volumetric flaw and a crack-like flaw. The elastodynamic Kirchhoff approximation (EKA) and the geometrical theory of diffraction (GTD) are applied respectively to a cylindrical cavity and a semi-infinite crack. These methods are known as high frequency approximations. The 2-D elastodynamic scattering problems of a plane wave incident on these model defects are considered and the scattered fields are expressed in terms of the reflection and diffraction coefficients. The ratio of the scattered far field amplitude to the incident wave amplitude is computed as a function of the angular location and compared with the boundary element solutions.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1994.04a
• /
• pp.1-8
• /
• 1994

An automated procedure which allows adaptation of spatial and time discretization simultaneously in finite element analysis of linear elastodynamic problems is presented. For dynamic problems having responses dominated by high frequency modes, such as those with impact, explosive, traveling and earthquake loads high gradient stress regions change their locations from time to time. And the time step size may need to vary in order to deal with whole process ranging from transient phase to steady state phase. As the sizes of elements in space vary in different regions, the procedure also permits different time stepping. In such a way, the best performance attainable by the finite element method can be achieved. In this study, we estimate both of the kinetic energy error and stran energy error induced by spatial and time discretization in a consistent manner. Numerical examples are used to demonstrate the performance of the procedure.

• Journal of the Korean Society of Industry Convergence
• /
• v.12 no.3
• /
• pp.149-155
• /
• 2009

A new finite element equation is derived by applying quadratic and cubic time integration scheme to the variational formulation in time-integral for the analysis of the transient elastodynamic problems to increase the numerical accuracy and stability. Emphasis is focused on methodology for cubic time integration scheme procedure which are never presented before. In this semidiscrete approximations of the field variables, the time axis is divided equally and quadratic and cubic time variation is assumed in those intervals, and space is approximated by the usual finite element discretization technique. It is found that unconditionally stable numerical results are obtained in case of the cubic time variation. Some numerical examples are given to show the versatility of the presented formulation.

• Journal of the Korean Society for Precision Engineering
• /
• v.9 no.2
• /
• pp.81-91
• /
• 1992

A control methodology for suppressing the elastodynamic responses of high-speed flexible linkage mechanisms is presented by adopting the concept of smart structures featuring piezoelectric films. The dynamic modeling of the proposed mechanism is accomplished by employing a finite element formulation which accounts for dynamic motion in both inertial and elastic coordinates. The dynamics of piezoelectric actuators and sensors bonded on the original flexible structure are developed for one-dimensional beam in conjunction with the modal analysis. The linear optimal controller which consists of a feedback control law and a Luenberger observer is employed. Numerical simulation is performed to evaluate the improvement of elastodynamic responses.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.23 no.11 s.170
• /
• pp.1993-2006
• /
• 1999

A Volume Integral Equation Method (VIEM) is applied for the effective analysis of elastic wave scattering problems and plane elastostatic problems in unbounded solids containing general anisotropic inclusions. It should be noted that this newly developed numerical method does not require the Green's function for anisotropic inclusions to solve this class of problems since only Green's function for the unbounded isotropic matrix is involved in their formulation for the analysis. This new method can also be applied to general two-dimensional elastodynamic and elastostatic problems with arbitrary shapes and number of anisotropic inclusions and voids. Through the analysis of plane elastodynamic and elastostatic problems in unbounded isotropic matrix with orthotropic inclusions and voids, it will be established that this new method is very accurate and effective for solving plane elastic problems in unbounded solids containing general anisotropic inclusions and voids.

• 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.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.12
• /
• pp.3109-3117
• /
• 1994

To operate a flexible mechanism in high speed its weight must be reduced as far as the structural strength does not decrease too much, but a light-weighted mechanism causes undesirable elastodynamic responses deteriorating the system performance. Besides, clearance within the connections of mechanisms causes rapid wear, increased noise and vibration. Even if the problems described above must be considered in the initial design stage, there has been no effective design process which takes account of the correlation between dynamic characteristics of flexible mechanism and the clearance effect at the joint. In this study, the generalized elastodynamic governing equations which include dynamic characteristics and boundary conditions of flexible mechanism are derived by variational calculus and solved by using FFM theory. To take the clearance effect at joint into account a new dynamic model is presented and also the method of modified stiffness/damping matrix is proposed to activate the dynamic clearance model, which cooperates with the developed governing equation very easily. As the results of this study, the proposed method(modified stiffness/damping matrix) to calculate clearance effect was proved to be superior to the existing one(force reaction method) in solution convergency and calculation performance. Besides this method can be easily adopted to the complex shape joint without calculation of reaction force direction.