The code above calculates the response (curvature) to a moment. If you want to calculate the of a rectangular plate under uniform pressure $q_0$:
coordinate system of the plate using the transformation matrix based on the fiber angle SCIRP Open Access 4. Compute the ABD Matrix The ABD matrix relates the applied forces and moments to the mid-plane strains and curvatures SCIRP Open Access Composite Plate Bending Analysis With Matlab Code
) matrices are computed by integrating through the thickness. The code above calculates the response (curvature) to
% Bending stiffness matrix D = zeros(3,3); for i = 1:length(t) zi = z(i); D = D + (1/3) * (zi^3 - zi^2 * t(i)) * [Q11, Q12, Q16; Q12, Q22, Q26; Q16, Q26, Q66]; end Composite Plate Bending Analysis With Matlab Code