We use a (size 2a×2b in local coordinates). Each node has 3 DOF: w, θx = ∂w/∂y, θy = -∂w/∂x.
fprintf('\nMaximum deflection: %.4f mm\n', max(W)*1000); fprintf('Minimum deflection: %.4f mm\n', min(W)*1000); Composite Plate Bending Analysis With Matlab Code
Running the code for a 0.5m×0.5m, 5mm thick cross-ply [0/90]s plate under 1 kPa pressure yields: We use a (size 2a×2b in local coordinates)
The code below solves a under uniform pressure. We compare center deflection with analytical series solution. θx = ∂w/∂y