Abstract
For linear elasticity problems the finite element method is an extremely successful method to model complicated structures. The successful implementation requires the solution of very large, sparse, positive definite linear systems of algebraic equations. A new technique for solving these systems using the preconditioned conjugate gradient method is proposed. Using ideas from both additive Schwarz methods and iterative substructuring methods, it is proven that the condition number of the resulting system does not grow as the substructures are made smaller and the mesh is refined. This result holds for two and three dimensions. Numerical experiments have been performed to demonstrate the power of this method. For linear elasticity problems in two dimensions the condition numbers are observed numerically to be less than four when using a regular mesh.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
More From: SIAM Journal on Scientific and Statistical Computing
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.