• Computational Structural Engineering
• /
• v.7 no.2
• /
• pp.69-75
• /
• 1994

The procedure for the stress and displacement analysis of realistic viscoelastic materials by time domain boundary element method(BEM) has been discussed. The fundamental solutions and stress kernels have been obtained using the elastic-viscoelastic correspondence principle. The relaxation function is expanded in a sum of exponentials and the transformed fundamental solutions and stress kernels are inverted numerically into real time space. The proposed procedure requires a small computational effort and it is applicable in time domain boundary element analysis of realistic viscoelastic problems. Numerical results of example problem show the effectiveness and applicability of the proposed method.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2016.05a
• /
• pp.470-474
• /
• 2016

Madsen et al. (2002)이 제안한 일차원 고차 Boussinesq 방정식에 대하여 불연속갤러킨 유한요소법(Discontinuous Galerkin Finite Element Method)을 적용하였다. 연속적인 Boussinesq 방정식에서 각 요소경계에 불연속을 허용할 수 있도록 공간차분하고, 시간방향으로 4차 Runge-Kutta 시간적분법, 각 요소사이에는 Lax-Friedrichs 수치흐름률을 사용하였다. 계산영역의 양쪽에 불필요한 파랑의 반사를 억제하도록 흡수층을 설치하였으며, 영역 내부에서 조파할 수 있도록 하였다. Luth et al.(1994)의 수중잠제 실험에 적용하여 관측값과 잘 일치함을 확인하였다.

• Journal of Power System Engineering
• /
• v.16 no.5
• /
• pp.40-46
• /
• 2012

For free surface flow problem, a high-order spectral/boundary element method is adapted as an efficient numerical tool. This method is one of the most efficient numerical methods by which the nonlinear gravity waves can be simulated and hydrodynamic forces also can be calculated in time domain. In this method, the velocity potential is expressed as the sum of surface potential and body potential. Then, surface potential is solved by using the high-order spectral method and body potential is solved by using the high-order boundary element method. Using the combination of these two methods, the free surface flow problems of a submerged moving body are solved in time domain. In the present study, lifting surface theory is added to the former work to include effects of lift force. Therefore, a new formulation for the basic mathematical theory is introduced to contain the lift body in calculation.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.26 no.4
• /
• pp.213-221
• /
• 2013

This paper introduces a mesh continuation scheme for a one-dimensional inverse medium problem to reconstruct the spatial distribution of elastic wave velocities in heterogeneous semi-infinite solid domains. To formulate the inverse problem, perfectly-matched-layers(PMLs) are introduced as wave-absorbing boundaries that surround the finite computational domain truncated from the originally semi-infinite extent. To tackle the inverse problem in the PML-truncated domain, a partial-differential-equations(PDE)-constrained optimization approach is utilized, where a least-squares misfit between calculated and measured surface responses is minimized under the constraint of PML-endowed wave equations. The optimization problem iteratively solves for the unknown wave velocities with their updates calculated by Fletcher-Reeves conjugate gradient algorithms. The optimization is performed using a mesh continuation scheme through which the wave velocity profile is reconstructed in successively denser mesh conditions. Numerical results showed the robust performance of the mesh continuation scheme in reconstructing target wave velocity profile in a layered heterogeneous solid domain.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.3
• /
• pp.25-34
• /
• 1993

The long-term behavior, such as in excavation problems of weak medium, can be dealt with by the elasto-viscoplasticity models. In this paper, a combined formulation of elasto-viscoplasticity using boundary elements and finite elements without using internal cells is presented. The domain integral introduced due to the viscoplastic stresses is transformed into a boundary integral applying direct integration in cylindrical coordinates. The results of the developed boundary element analysis are compared with those from the explicit solution and from the finite element analysis. It is observed that the boundary element analysis without internal cells results in some error because of its deficiency in handling the nonlinearity in local stress concentration. Therefore, a coupled analysis of boundary elements and finite elements, in which finite elements are used in the area of stress concentration, is developed. The coupled method is applied to a time dependent inelastic problem with semi-infinite boundaries. It results in reasonable solution compared with other methods where relatively higher degree of freedoms are employed. Thus, it is concluded that the combined analysis may be used for such problems in the effective manner.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.31 no.2
• /
• pp.87-95
• /
• 2018

This paper introduces a new explicit spectral element method for the simulation of SH-waves in semi-infinite domains. To simulate the wave motion in unbounded domains, it is necessary to reduce the infinite extent to a finite computational domain of interest. To prevent the wave reflection from the trunctated boundaries, perfectly matched layer(PML) wave-absorbing boundary is introduced. The forward problem for simulating SH-waves in PML-truncated domains can be formulated as second-order PDEs. The second-order semi-discrete form of the governing PDEs is constructed by using a mixed spectral elements with Legendre-gauss-Lobatto quadrature method, which results in a diagonalized mass matrix. Then the second-order semi-discrete form is transformed to a first-order, whose solutions are calculated by the fourth-order Runge-Kutta method. Numerical examples showed that solutions of SH-wave in the two-dimensional analysis domain resulted in stable and accurate, and reflections from truncated boundaries could be reduced by using PML boundaries. Elastic wave propagation analysis using explicit time integration method may be apt for solving larger domain problems such as three-dimensional elastic wave problem more efficiently.

• Computational Structural Engineering
• /
• v.9 no.1
• /
• pp.85-91
• /
• 1996

The focus of the present work is on the computation of the stress intensity factor for the crack at the elastic-viscoelastic bimaterial interface. First, the stress intensity factor for an interface crack in dissimilar elastic and viscoelastic materials is dervied by applying the correspondence principle to associated elastic expression. Then the time-domain boundary element analysis is performed to calculate the stress intensity factor. Numerical results show that the proposed method is very useful for the analysis of the interface crack in elastic and viscoelastic materials.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.28 no.2
• /
• pp.123-130
• /
• 2015

This paper presents a new formulation for transient simulations of microwave propagation in heterogeneous unbounded domains. In particular, perfectly-matched-layers(PMLs) are introduced to allow for wave absorption at artificial boundaries used to truncate the infinite extent of the physical domains. The development of the electromagnetic PML targets the application to engineering mechanics problems such as structural health monitoring and inverse medium problems. To formulate the PML for plane electromagnetic waves, a complex coordinate transformation is introduced to Maxwell's equations in the frequency-domain. Then the PML-endowed partial differential equations(PDEs) for transient electromagnetic waves are recovered by the application of the inverse Fourier transform to the frequency-domain equations. A mixed finite element method is employed to solve the time-domain PDEs for electric and magnetic fields in the PML-truncated domain. Numerical results are presented for plane microwaves propagating through concrete structures, and the accuracy of solutions is investigated by a series of error analyses.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.26 no.4
• /
• pp.223-231
• /
• 2013

This paper presents a time-domain Gauss-Newton full-waveform inversion method for the material profile reconstruction in heterogeneous semi-infinite solid media. To implement the inverse problem in a finite computational domain, perfectly-matchedlayers( PMLs) are introduced as wave-absorbing boundaries within which the domain's wave velocity profile is to be reconstructed. The inverse problem is formulated in a partial-differential-equations(PDE)-constrained optimization framework, where a least-squares misfit between measured and calculated surface responses is minimized under the constraint of PML-endowed wave equations. A Gauss-Newton-Krylov optimization algorithm is utilized to iteratively update the unknown wave velocity profile with the aid of a specialized regularization scheme. Through a series of one-dimensional examples, the solution of the Gauss-Newton inversion was close enough to the target profile, and showed superior convergence behavior with reduced wall-clock time of implementation compared to a conventional inversion using Fletcher-Reeves optimization algorithm.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.9
• /
• pp.56-65
• /
• 2005

A dynamic analysis method that freezes a time domain by discretization and solves the spatial propagation equation has a unique feature that provides a degree of freedom on spatial domain compared with the space discretization or space-time discretization finite element method. Using this feature, the time finite element analysis can be effectively applied to optimize the spatial characteristics of distributed type actuators. In this research, the time domain finite element method was used to discretize the model. A state variable vector was used in the discretization to include arbitrary initial conditions. A performance index was proposed on spatial domain to consider both potential and vibrational energy, so that the resulting shape of the distributed actuator was optimized for dynamic control of the structure. It is assumed that the structure satisfies the final rest condition using the realizable control scheme although the initial disturbance can affect the system response. Both equations on states and costates were derived based on the selected performance index and structural model. Ricatti matrix differential equations on state and costate variables were derived by the reconfiguration of the sub-matrices and application of time/space boundary conditions, and finally optimal actuator distribution was obtained. Numerical simulation results validated the proposed actuator shape optimization scheme.