TY - JOUR
T1 - RAPID analysis of variable stiffness beams and plates
T2 - Legendre polynomial triple-product formulation
AU - O'Donnell, Matthew P.
AU - Weaver, Paul M.
N1 - Publisher Copyright:
Copyright © 2017 John Wiley & Sons, Ltd.
PY - 2017/10/5
Y1 - 2017/10/5
N2 - Numerical integration techniques are commonly employed to formulate the system matrices encountered in the analysis of variable stiffness beams and plates using a Ritz based approach. Computing these integrals accurately is often computationally costly. Herein, a novel alternative is presented, the Recursive Analytical Polynomial Integral Definition (RAPID) formulation. The RAPID formulation offers a significant improvement in the speed of analysis, achieved by reducing the number of numerical integrations that are performed by an order of magnitude. A common Legendre Polynomial basis is employed for both trial functions and stiffness/load variations leading to a common form for the integrals encountered. The Legendre Polynomial basis possesses algebraic recursion relations that allow these integrals to be reformulated as triple-products with known analytical solutions, defined compactly using the Wigner (3j) coefficient. The satisfaction of boundary conditions, calculation of derivatives and transformation to other bases is achieved through combinations of matrix multiplication, with each matrix representing a unique boundary condition or physical effect, therefore permitting application of the RAPID approach to a variety of problems. Indicative performance studies demonstrate the advantage of the RAPID formulation when compared to direct analysis using Matlab's ‘integral’ and ‘integral2’.
AB - Numerical integration techniques are commonly employed to formulate the system matrices encountered in the analysis of variable stiffness beams and plates using a Ritz based approach. Computing these integrals accurately is often computationally costly. Herein, a novel alternative is presented, the Recursive Analytical Polynomial Integral Definition (RAPID) formulation. The RAPID formulation offers a significant improvement in the speed of analysis, achieved by reducing the number of numerical integrations that are performed by an order of magnitude. A common Legendre Polynomial basis is employed for both trial functions and stiffness/load variations leading to a common form for the integrals encountered. The Legendre Polynomial basis possesses algebraic recursion relations that allow these integrals to be reformulated as triple-products with known analytical solutions, defined compactly using the Wigner (3j) coefficient. The satisfaction of boundary conditions, calculation of derivatives and transformation to other bases is achieved through combinations of matrix multiplication, with each matrix representing a unique boundary condition or physical effect, therefore permitting application of the RAPID approach to a variety of problems. Indicative performance studies demonstrate the advantage of the RAPID formulation when compared to direct analysis using Matlab's ‘integral’ and ‘integral2’.
KW - composites
KW - integration
KW - structures
UR - http://www.scopus.com/inward/record.url?scp=85017336994&partnerID=8YFLogxK
U2 - 10.1002/nme.5528
DO - 10.1002/nme.5528
M3 - Article
AN - SCOPUS:85017336994
SN - 0029-5981
VL - 112
SP - 86
EP - 100
JO - International Journal for Numerical Methods in Engineering
JF - International Journal for Numerical Methods in Engineering
IS - 1
ER -