A numerical algorithm based on finite element method (FEM) and boundary element method (BEM) is proposed for the analyses of acoustic-structure coupled problems. By this algorithm, the structural domain and the acoustic domain are modelled by FEM and BEM, respectively, which are coupled with each other through the consideration of the appropriate compatibility and equilibrium conditions on the interface of the two domains. To improve the computational efficiency, the adaptive cross approximation (ACA) approach is incorporated into the proposed algorithm to deal with the nonsymmetric and fully populated matrices resulting from the coupling of the FEM and BEM. The validity and the high efficiency of the present approach are demonstrated by two examples.
Purpose
Based on the error analysis, the authors proposed a new kind of high accuracy boundary element method (BEM) (HABEM), and for the large-scale problems, the fast algorithm, such as adaptive cross approximation (ACA) with generalized minimal residual (GMRES) is introduced to develop the high performance BEM (HPBEM). It is found that for slender beams, the stress analysis using iterative solver GMRES will difficult to converge. For the analysis of slender beams and thin structures, to enhance the efficiency of GMRES solver becomes a key problem in the development of the HPBEM. The purpose of this paper is study on the preconditioning method to solve this convergence problem, and it is started from the 2D BE analysis of slender beams.
Design/methodology/approach
The conventional sparse approximate inverse (SAI) based on adjacent nodes is modified to that based on adjacent nodes along the boundary line. In addition, the authors proposed a dual node variable merging (DNVM) preprocessing for slender thin-plate beams. As benchmark problems, the pure bending of thin-plate beam and the local stress analysis (LSA) of real thin-plate cantilever beam are applied to verify the effect of these two preconditioning method.
Findings
For the LSA of real thin-plate cantilever beams, as GMRES (m) without preconditioning applied, it is difficult to converge provided the length to height ratio greater than 50. Even with the preconditioner SAI or DNVM, it is also difficult to obtain the converged results. For the slender real beams, the iteration of GMRES (m) with SAI or DNVM stopped at wrong deformation state, and the computation failed. By changing zero initial solution to the analytical displacement solution of conventional beam theory, GMRES (m) with SAI or DNVM will not be stopped at wrong deformation state, but the stress error is still difficult to converge. However, by GMRES (m) combined with both SAI and DNVM preconditioning, the computation efficiency enhanced significantly.
Originality/value
This paper presents two preconditioners: DNVM and a modified SAI based on adjacent nodes along the boundary line of slender thin-plate beam. In the LSA, by using GMRES (m) combined with both DNVM and SAI, the computation efficiency enhanced significantly. It provides a reference for the further development of the 3D HPBEM in the LSA of real beam, plate and shell structures.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.