## 3 Mixed FE-GDQ Numerical Formulation
\begin{table}
\begin{tabular}{c c c c c c c c c c c} \hline \multirow{2}{*}{**Aspect Ratio**} & \multirow{2}{*}{**Buckling Load (10 GPa)**} & \multicolumn{3}{c}{**PNND**} & \multicolumn{3}{c}{**PNND**} & \multicolumn{3}{c}{**PUD**} \\ \cline{3-11}  & & \(\mathbf{e_{0}=0.3}\) & \(\mathbf{e_{0}=0.6}\) & \(\mathbf{e_{0}=0.9}\) & \(\mathbf{e_{0}=0.3}\) & \(\mathbf{e_{0}=0.6}\) & \(\mathbf{e_{0}=0.9}\) & \(\mathbf{e_{0}=0.3}\) & \(\mathbf{e_{0}=0.6}\) & \(\mathbf{e_{0}=0.9}\) \\ \hline \multirow{4}{*}{_a/b = 1_} & \(\lambda_{1}\) & 2.666 & 1.935 & 1.226 & 2.674 & 1.961 & 1.263 & 2.269 & 1.212 & 0.275 \\  & \(\lambda_{2}\) & 2.666 & 1.936 & 1.226 & 2.674 & 1.961 & 1.263 & 2.269 & 1.212 & 0.275 \\  & \(\lambda_{3}\) & 2.958 & 2.255 & 1.442 & 2.998 & 2.265 & 1.390 & 2.595 & 1.482 & 0.354 \\  & \(\lambda_{4}\) & 2.984 & 2.261 & 1.447 & 3.016 & 2.266 & 1.390 & 2.607 & 1.483 & 0.354 \\  & \(\lambda_{1}\) & 2.281 & 1.863 & 1.162 & 2.392 & 2.005 & 1.269 & 2.083 & 1.281 & 0.2952 \\  & \(\lambda_{2}\) & 2.339 & 1.865 & 1.162 & 2.438 & 2.029 & 1.271 & 2.105 & 1.282 & 0.2954 \\  & \(\lambda_{3}\) & 2.558 & 1.968 & 1.237 & 2.627 & 2.068 & 1.379 & 2.248 & 1.334 & 0.3464 \\  & \(\lambda_{4}\) & 2.719 & 2.044 & 1.307 & 2.761 & 2.094 & 1.388 & 2.350 & 1.356 & 0.3467 \\ \hline \end{tabular}
\end{table} 
Table 2: First four shear buckling loads for a SSSS FG-saturated porous plate with different aspect ratios, _a/b_, porosity distributions and porosity coefficients, \(e_{0}\), keeping \(\beta=0.6\) fixed. Figure 2: Comparative evaluation of the first three buckling mode shapes for a SSSS homogenous plate under a shear loading condition. \begin{table}
\begin{tabular}{c c c c c c c c c c c} \hline \hline \multirow{2}{*}{**Aspect Ratio**} & \multirow{2}{*}{**Buckling Load (10 GPa)**} & \multicolumn{4}{c}{**PNND**} & \multicolumn{4}{c}{**PNND**} & \multicolumn{4}{c}{**PUD**} \\ \cline{3-11}  & & \(\mathbf{c_{0}=0.3}\) & \(\mathbf{c_{0}=0.6}\) & \(\mathbf{c_{0}=0.9}\) & \(\mathbf{c_{0}=0.3}\) & \(\mathbf{c_{0}=0.6}\) & \(\mathbf{c_{0}=0.9}\) & \(\mathbf{c_{0}=0.3}\) & \(\mathbf{c_{0}=0.6}\) & \(\mathbf{c_{0}=0.9}\) \\ \hline \multirow{6}{*}{\(a\)/\(b\) = 1} & \(\lambda_{1}\) & 1.143 & 0.897 & 0.536 & 1.199 & 1.006 & 0.692 & 1.020 & 0.636 & 0.177 \\  & \(\lambda_{2}\) & 1.150 & 0.900 & 0.543 & 1.208 & 1.007 & 0.693 & 1.020 & 0.645 & 0.182 \\  & \(\lambda_{3}\) & 1.653 & 1.267 & 0.773 & 1.705 & 1.434 & 0.845 & 1.455 & 0.865 & 0.219 \\  & \(\lambda_{4}\) & 1.664 & 1.275 & 0.778 & 1.717 & 1.350 & 0.849 & 1.464 & 0.868 & 0.219 \\  & \(\lambda_{1}\) & 0.310 & 0.255 & 0.153 & 0.320 & 0.292 & 0.241 & 0.287 & 0.186 & 0.055 \\  & \(\lambda_{2}\) & 0.497 & 0.393 & 0.229 & 0.538 & 0.470 & 0.375 & 0.449 & 0.281 & 0.078 \\ \(a\)/\(b\) = 2 & \(\lambda_{3}\) & 0.598 & 0.466 & 0.281 & 0.610 & 0.520 & 0.405 & 0.528 & 0.325 & 0.096 \\  & \(\lambda_{4}\) & 0.805 & 0.622 & 0.374 & 0.817 & 0.682 & 0.504 & 0.708 & 0.430 & 0.116 \\ \hline \hline \end{tabular}
\end{table} 
Table 4: First four shear buckling loads for a CFCF FG-saturated porous plate with different aspect ratios, \(a/b\), porosity distributions and porosity coefficients, \(e_{0}\), keeping \(\beta=0.6\) fixed. Figure 3: Comparative evaluation of the first three buckling mode shapes for a SSSS homogenous plate under a biaxial compression.