Transient elastodynamic analysis with a combination of convolution quadrature method and pseudo-initial condition method

2018 ◽  
Vol 36 (1) ◽  
pp. 334-355
Author(s):  
Yuan Li ◽  
J. Zhang ◽  
Yudong Zhong ◽  
Xiaomin Shu ◽  
Yunqiao Dong
Keyword(s):  
Initial Conditions ◽  
Initial Condition ◽  
Quadrature Method ◽  
Memory Requirement ◽  
Content Type ◽  
Convolution Quadrature ◽  
Force Method ◽  
Cpu Time ◽  
Condition Method ◽  
Convolution Quadrature Method

Purpose The Convolution Quadrature Method (CQM) has been widely applied to solve transient elastodynamic problems because of its stability and generality. However, the CQM suffers from the problems of huge memory requirement in case of direct implementation in time domain or CPU time in case of its reformulation in Laplace domain. The purpose of this paper is to combine the CQM with the pseudo-initial condition method (PICM) to achieve a good balance between memory requirement and CPU time. Design/methodology/approach The combined methods first subdivide the whole analysis into a few sub-analyses, which is dealt with the PICM, namely, the results obtained by previous sub-analysis are used as the initial conditions for the next sub-analysis. In each sub-analysis, the time interval is further discretized into a number of sub-steps and dealt with the CQM. For non-zero initial conditions, the pseudo-force method is used to transform them into equivalent body forces. The boundary face method is employed in the numerical implementation. Three examples are analyzed. Results are compared with analytical solutions or FEM results and the results of reformulated CQM. Findings Results demonstrate that the computation time and the storage requirement can be reduced significantly as compared to the CQM, by using the combined approach. Originality/value The combined methods can be successfully applied to the problems of long-time dynamic response, which requires a large amount of computer memory when CQM is applied, while preserving the CQM stability. If the number of time steps is high, then the accuracy of the proposed approach can be deteriorated because of the pseudo-force method.

Boundary Element Method with Runge-Kutta Convolution Quadrature for Three-Dimensional Dynamic Poroelasticity

2014 ◽  
Vol 709 ◽  
pp. 101-104
Author(s):  
L.A. Igumnov ◽  
I.V. Vorobtsov ◽  
S.Yu. Litvinchuk
Keyword(s):  
Three Dimensional ◽  
Laplace Domain ◽  
Convolution Quadrature ◽  
Time Stepping ◽  
The Face ◽  
Highly Oscillatory ◽  
Runge Kutta ◽  
Oscillatory Quadrature ◽  
Convolution Quadrature Method ◽  
Laplace Transform Inversion

The paper contains a brief introduction to the state of the art in poroelasticity models, in BIE & BEM methods application to solve dynamic problems in Laplace domain. Convolution Quadrature Method is formulated, as well as Runge-Kutta convolution quadrature modification and scheme with a key based on the highly oscillatory quadrature principles. Several approaches to Laplace transform inversion, including based on traditional Euler stepping scheme and Runge-Kutta stepping schemes, are numerically compared. A BIE system of direct approach in Laplace domain is used together with the discretization technique based on the collocation method. The boundary is discretized with the quadrilateral 8-node biquadratic elements. Generalized boundary functions are approximated with the help of the Goldshteyn’s displacement-stress matched model. The time-stepping scheme can rely on the application of convolution theorem as well as integration theorem. By means of the developed software the following 3d poroelastodynamic problem were numerically treated: a Heaviside-shaped longitudinal load acting on the face of a column.

Sign in / Sign up

Export Citation Format

Share Document