## 3 Mixed FE-GDQ Numerical Formulation
\begin{table}
\begin{tabular}{c c c c c c c c c c} \hline \hline \multirow{2}{*}{**Aspect Ratio**} & \multirow{2}{*}{**Buckling Load (10 GPa)**} & \multicolumn{3}{c}{**PNND**} & \multicolumn{3}{c}{**PND**} & \multicolumn{3}{c}{**PUD**} \\ \cline{3-11}  & & \(\beta=0.0\) & \(\beta=0.6\) & \(\beta=0.9\) & \(\beta=0.0\) & \(\beta=0.6\) & \(\beta=0.9\) & \(\beta=0.0\) & \(\beta=0.6\) & \(\beta=0.9\) \\ \hline \multirow{4}{*}{\(a/b=1\)} & \(\lambda_{1}\) & 0.288 & 0.303 & 0.325 & 0.361 & 0.362 & 0.389 & 0.198 & 0.216 & 0.246 \\  & \(\lambda_{2}\) & 0.683 & 0.723 & 0.781 & 0.825 & 0.830 & 0.846 & 0.472 & 0.523 & 0.603 \\  & \(\lambda_{3}\) & 0.683 & 0.723 & 0.781 & 0.825 & 0.830 & 0.846 & 0.472 & 0.523 & 0.603 \\  & \(\lambda_{4}\) & 0.103 & 0.108 & 0.146 & 0.123 & 0.127 & 0.131 & 0.711 & 0.763 & 0.843 \\  & \(\lambda_{1}\) & 0.171 & 0.186 & 0.207 & 0.213 & 0.221 & 0.229 & 0.119 & 0.138 & 0.168 \\  & \(\lambda_{2}\) & 0.319 & 0.369 & 0.362 & 0.397 & 0.400 & 0.405 & 0.219 & 0.239 & 0.273 \\  & \(\lambda_{3}\) & 0.513 & 0.564 & 0.635 & 0.616 & 0.634 & 0.638 & 0.362 & 0.426 & 0.480 \\  & \(\lambda_{4}\) & 0.573 & 0.603 & 0.644 & 0.697 & 0.702 & 0.713 & 0.393 & 0.427 & 0.526 \\ \hline \hline \end{tabular}
\end{table} 
Table 10: First four biaxial buckling loads for a SSSS FG-saturated porous plate with different aspect ratios, \(a/b\), porosity distributions and Skempton coefficients, \(\beta\), while keeping \(e_{0}=0.6\) fixed. \begin{table}
\begin{tabular}{c c c c c c c c c c} \hline \hline \multirow{2}{*}{**Aspect Ratio**} & \multirow{2}{*}{**Buckling Load (10 GPa)**} & \multicolumn{3}{c}{**PNND**} & \multicolumn{3}{c}{**PNDD**} & \multicolumn{3}{c}{**PUD**} \\ \cline{3-11}  & & \(\beta=0.0\) & \(\beta=0.6\) & \(\beta=0.9\) & \(\beta=0.0\) & \(\beta=0.6\) & \(\beta=0.9\) & \(\beta=0.0\) & \(\beta=0.6\) & \(\beta=0.6\) & \(\beta=0.9\) \\ \hline \multirow{4}{*}{\(a/b=1\)} & \(\lambda_{1}\) & 0.848 & 0.897 & 0.955 & 0.969 & 0.978 & 0.989 & 0.583 & 0.636 & 0.709 \\  & \(\lambda_{2}\) & 0.853 & 0.900 & 0.967 & 0.978 & 0.983 & 0.993 & 0.583 & 0.645 & 0.728 \\  & \(\lambda_{3}\) & 1.231 & 1.267 & 1.299 & 1.329 & 1.337 & 1.346 & 0.831 & 0.865 & 0.878 \\  & \(\lambda_{4}\) & 1.240 & 1.275 & 1.306 & 1.338 & 1.343 & 1.371 & 0.836 & 0.868 & 0.879 \\  & \(\lambda_{1}\) & 0.903 & 0.954 & 1.021 & 1.031 & 1.039 & 1.042 & 0.618 & 0.680 & 0.765 \\  & \(\lambda_{2}\) & 0.904 & 0.955 & 1.026 & 1.032 & 1.040 & 1.068 & 0.619 & 0.680 & 0.765 \\  & \(\lambda_{3}\) & 1.318 & 1.360 & 1.399 & 1.423 & 1.402 & 1.416 & 0.890 & 0.931 & 0.955 \\ \hline \hline \end{tabular}
\end{table} 
Table 8: First four shear buckling loads for a FCFC FG-saturated porous plate with different aspect ratios, \(a/b\), porosity distributions and Skempton coefficients, \(\beta\), while keeping \(e_{0}=0.6\) fixed. \begin{table}
\begin{tabular}{c c c c c c c c c c} \hline \hline \multirow{2}{*}{**Aspect Ratio**} & \multirow{2}{*}{**Buckling Load (10 GPa)**} & \multicolumn{3}{c}{**PNND**} & \multicolumn{3}{c}{**PNDD**} & \multicolumn{3}{c}{**PUD**} \\ \cline{3-11}  & & \(\beta=0.0\) & \(\beta=0.6\) & \(\beta=0.9\) & \(\beta=0.0\) & \(\beta=0.6\) & \(\beta=0.9\) & \(\beta=0.0\) & \(\beta=0.6\) & \(\beta=0.9\) \\ \hline \multirow{4}{*}{\(a/b=1\)} & \(\lambda_{1}\) & 0.848 & 0.897 & 0.955 & 0.969 & 0.978 & 0.989 & 0.583 & 0.636 & 0.709 \\  & \(\lambda_{2}\) & 0.853 & 0.900 & 0.967 & 0.978 & 0.983 & 0.993 & 0.583 & 0.645 & 0.728 \\  
Figure 5: First four buckling mode shapes of a FG-saturated porous plate subjected to a shear load (FCFC, \(a=2\) m, \(b=1\) m). Figure 6: _Cont._  
Figure 6: First four buckling mode shapes of a FG-saturated porous plate subjected to a shear load (SSSS, \(a=2\) m, \(b=1\) m). Figure 7: First four buckling mode shapes of a FG-saturated porous plate subjected to uniaxial load (SSSS, \(a=2\) m, \(b=1\) m).