In the present paper a fast solver for dual boundary element analysis of 3D anisotropic crack problems is formulated, implemented and tested. The fast solver is based on the use of hierarchical matrices for the representation of the collocation matrix. The admissible low rank blocks are computed by adaptive cross approximation (ACA). The performance of ACA against the accuracy of the adopted computational scheme for the evaluation of the anisotropic kernels is investigated, focusing on the balance between the kernel representation accuracy and the accuracy required for ACA. The system solution is computed by a preconditioned GMRES and the preconditioner is built exploiting the hierarchical arithmetic and taking full advantage of the hierarchical format. The effectiveness of the proposed technique for anisotropic crack problems has been numerically demonstrated, highlighting the accuracy as well as the significant reduction in memory storage and analysis time. In particular, it has been numerically shown that the computational cost grows almost linearly with the number of degrees of freedom, obtaining up to solution speedups of order 10 for systems of order 104. Moreover, the sensitivity of the performance of the numerical scheme to materials with different degrees of anisotropy has been assessed.

Benedetti, I., Milazzo, A., Aliabadi, M.H. (2009). A fast dual boundary element method for 3D anisotropic crack problems. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, 80, 1356-1378 [10.1002/nme.2666].

A fast dual boundary element method for 3D anisotropic crack problems

BENEDETTI, Ivano;MILAZZO, Alberto;
2009-01-01

Abstract

In the present paper a fast solver for dual boundary element analysis of 3D anisotropic crack problems is formulated, implemented and tested. The fast solver is based on the use of hierarchical matrices for the representation of the collocation matrix. The admissible low rank blocks are computed by adaptive cross approximation (ACA). The performance of ACA against the accuracy of the adopted computational scheme for the evaluation of the anisotropic kernels is investigated, focusing on the balance between the kernel representation accuracy and the accuracy required for ACA. The system solution is computed by a preconditioned GMRES and the preconditioner is built exploiting the hierarchical arithmetic and taking full advantage of the hierarchical format. The effectiveness of the proposed technique for anisotropic crack problems has been numerically demonstrated, highlighting the accuracy as well as the significant reduction in memory storage and analysis time. In particular, it has been numerically shown that the computational cost grows almost linearly with the number of degrees of freedom, obtaining up to solution speedups of order 10 for systems of order 104. Moreover, the sensitivity of the performance of the numerical scheme to materials with different degrees of anisotropy has been assessed.
2009
Settore ING-IND/04 - Costruzioni E Strutture Aerospaziali
Benedetti, I., Milazzo, A., Aliabadi, M.H. (2009). A fast dual boundary element method for 3D anisotropic crack problems. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, 80, 1356-1378 [10.1002/nme.2666].
File in questo prodotto:
File Dimensione Formato  
A-fast-dual-boundary-element-method-for-3D-anisotropic-crack-problems_2009_International-Journal-for-Numerical-Methods-in-Engineering.pdf

Solo gestori archvio

Dimensione 703.29 kB
Formato Adobe PDF
703.29 kB Adobe PDF   Visualizza/Apri   Richiedi una copia

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/10447/42004
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 49
  • ???jsp.display-item.citation.isi??? 44
social impact