TY - JOUR
T1 - RAPID analysis of variable stiffness beams and plates
T2 - Legendre polynomial triple-product formulation
AU - O'Donnell, Matt
AU - Weaver, Paul
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 (LP) basis is employed for both trial functions and stiffness/load variations leading to a common form for the integrals encountered. The LP 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 (LP) basis is employed for both trial functions and stiffness/load variations leading to a common form for the integrals encountered. The LP 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 - Structures
KW - Composites
KW - Integration
U2 - 10.1002/nme.5528
DO - 10.1002/nme.5528
M3 - Article
VL - 112
SP - 86
EP - 100
JO - International Journal for Numerical Methods in Engineering
JF - International Journal for Numerical Methods in Engineering
SN - 0029-5981
IS - 1
ER -