TY - JOUR
T1 - Schur decomposition in the scaled boundary finite element method in elastostatics
AU - Li, M.
AU - Song, H.
AU - Guan, H.
AU - Zhang, H.
PY - 2010
Y1 - 2010
N2 - The scaled boundary finite element method (SBFEM) was originally proposed for modelling elastodynamics in bounded and unbounded media. The method has demonstrated its superiority to the finite element method and the boundary element method when dealing with problems involving unbounded computational domains, or with difficulties of irregular frequencies and sharp corners. The SBFEM transforms the governing equations from partial differential equations to ordinary differential equations (ODEs). In addition, only the boundary of the study domain needs to be discretised which significantly reduces the computational cost. In the existing solution procedure, an eigenvalue problem of the Hamiltonian matrix, formulated from the coefficient matrices of ODEs, needs to be solved. Subsequently, the eigenvectors of the Hamiltonian matrix are arranged in a matrix form for the stiffness matrix in the nodal force-displacement relationship. However, the matrix formulated by the eigenvectors is close to singular when multiple eigenvalues with parallel eigenvectors exist, which leads to an inaccurate solution. In the present study, this problem is eliminated by using the Schur decomposition instead of the eigenvalue decomposition. A three-dimensional study of a cylindrical pile subjected to uniformly distributed load, is carried out. The performance and efficiency of the Schur decomposition are discussed in some detail for achieving more accurate solutions in using the SBFEM.
AB - The scaled boundary finite element method (SBFEM) was originally proposed for modelling elastodynamics in bounded and unbounded media. The method has demonstrated its superiority to the finite element method and the boundary element method when dealing with problems involving unbounded computational domains, or with difficulties of irregular frequencies and sharp corners. The SBFEM transforms the governing equations from partial differential equations to ordinary differential equations (ODEs). In addition, only the boundary of the study domain needs to be discretised which significantly reduces the computational cost. In the existing solution procedure, an eigenvalue problem of the Hamiltonian matrix, formulated from the coefficient matrices of ODEs, needs to be solved. Subsequently, the eigenvectors of the Hamiltonian matrix are arranged in a matrix form for the stiffness matrix in the nodal force-displacement relationship. However, the matrix formulated by the eigenvectors is close to singular when multiple eigenvalues with parallel eigenvectors exist, which leads to an inaccurate solution. In the present study, this problem is eliminated by using the Schur decomposition instead of the eigenvalue decomposition. A three-dimensional study of a cylindrical pile subjected to uniformly distributed load, is carried out. The performance and efficiency of the Schur decomposition are discussed in some detail for achieving more accurate solutions in using the SBFEM.
UR - http://www.scopus.com/inward/record.url?scp=84907694803&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=84907694803&partnerID=8YFLogxK
U2 - 10.1088/1757-899X/10/1/012243
DO - 10.1088/1757-899X/10/1/012243
M3 - Article
AN - SCOPUS:84907694803
SN - 1757-8981
VL - 10
SP - 1
EP - 9
JO - IOP Conference Series: Materials Science and Engineering
JF - IOP Conference Series: Materials Science and Engineering
IS - 1
M1 - 012243
ER -