% Local coordinates x = xi * a_elem; y = eta * b_elem;
% Stress = Q_bar * Strain stress_bot = Q_bar_store(:,:,k) * strain_bot; stress_top = Q_bar_store(:,:,k) * strain_top; Composite Plate Bending Analysis With Matlab Code
if norm(B,'fro') > 1e-6 warning('Unsymmetrical laminate: using bending stiffness only (approximate)'); end % Local coordinates x = xi * a_elem;
[ A_11 \frac\partial u_0\partial x + \dots - B_11 \frac\partial^2 w\partial x^2 - \dots = 0 ] [ D_11 \frac\partial^4 w\partial x^4 + \dots - B_11 \frac\partial^3 u_0\partial x^3 - \dots = q ] y = eta * b_elem