Functionally graded square and rectangular plates of different thicknesses placed on the elastic foundation modeled according to the Winkler-Pasternak theory have been studied. The thermal and mechanical characteristics, apart from Poisson’s ratio, are considered to continuously differ through the thickness of the studied material as stated in a power-law distribution. A mathematical model of functionally graded plate which include interaction with elastic foundation is defined. The equilibrium and stability equations are derived using high order shear deformation theory that comprises various kinds of shape function and the von Karman nonlinearity. A new analytically integrable shape function has been introduced. Hamilton’s principle has been applied with the purpose of acquiring the equations of motion. An analytical method for identifying both natural frequencies and critical buckling temperature for cases of linear and nonlinear temperature change through the plate thickness has been established. In order to verify the derived theoretical results on numerical examples, an original program code has been implemented within software MATLAB. Critical buckling temperature and natural frequencies findings are shown below. Previous scientific research and papers confirms that presented both the theoretical formulation and the numerical results are accurate. The comparison has been made between newly established findings based on introduced shape function and the old findings that include 13 different shape functions available in previously published articles. The final part of the research provides analysis and conclusions related to the impact of the power-law index, foundation stiffness, and temperature gradient on critical buckling temperature and natural frequencies of the functionally graded plates.

The bending analysis of thick and moderately thick functionally graded square and rectangular plates as well as plates on Winkler–Pasternak elastic foundation subjected to sinusoidal transverse load is presented in this paper. The plates are assumed to have isotropic, two-constituent material distribution through the thickness, and the modulus of elasticity of the plate is assumed to vary according to a power-law distribution in terms of the volume fractions of the constituents. This paper presents the methodology of the application of the high order shear deformation theory based on the shape functions. A new shape function has been developed and the obtained results are compared to the results obtained with 13 different shape functions presented in the literature. Also, the validity and accuracy of the developed theory was verified by comparing those results with the results obtained using the third order shear deformation theory and 3D theories. In order to determine the procedure for the analysis and the prediction of behavior of functionally graded plates, the new program code in the software package MATLAB has been developed based on the theories studied in this paper. The effects of transversal shear deformation, side-to-thickness ratio, and volume fraction distributions are studied and appropriate conclusions are given.

This study investigates the free vibration behaviors of functionally graded (FGM) plates with a porous structure, resting on a Kerr-type elastic foundation, while accounting for thermal effects and complex material property distributions. Within the framework of higher-order shear deformation theory (HSDT), two novel shape functions are introduced to accurately model transverse shear deformation across the plate thickness without employing shear correction factors. These functions are constructed to satisfy shear stress boundary conditions and capture nonlinear effects induced by material gradation and porosity. A variational formulation is developed to describe the dynamic response of FGM plates in a thermo-mechanical environment, incorporating temperature-dependent material properties and three porosity distributions: uniform, linear, and trigonometric. Numerical solutions are obtained using in-house MATLAB codes, allowing complete control over the formulation and interpretation of the results. The model is validated through detailed comparisons with existing literature, demonstrating high accuracy. The findings reveal that the porosity distribution pattern and gradient intensity significantly influence natural frequencies and mode shapes. The trigonometric porosity distribution exhibits favorable dynamic performance due to preserved stiffness in the surface regions. Additionally, the Kerr-type elastic foundation enables fine tuning of the dynamic response, depending on its specific parameters. The proposed approach provides a reliable and efficient tool for analyzing FGM structures under complex loading conditions and lays the groundwork for future extensions involving nonlinear, time-dependent, and multiphysics analyses.

A thermal buckling analysis of functionally graded thick rectangular plates accord?ing to von Karman non-linear theory is presented. The material properties of the functionally graded plate, except for the Poisson?s ratio, were assumed to be graded in the thickness direction, according to a power-law distribution, in terms of the volume fractions of the metal and ceramic constituents. Formulations of equilibrium and stability equations are derived using the high order shear deformation theory based on different types of shape functions. Analytical method for determination of the critical buckling temperature for uniform increase of temperature, linear and non-linear change of temperature across thickness of a plate is developed. Numeri?cal results were obtained in ?ATLAB software using combinations of symbolic and numeric values. The paper presents comparative results of critical buckling tempera?ture for different types of shape functions. The accuracy of the formulation presented is verified by comparing to results available from the literature.

The influence of various process parameters on the deep drawing process is a current research topic in sheet metal forming technology. Starting from the application of the previously constructed original testing device, an original tribological model was developed based on the process of sheet metal strip sliding between flat contact surfaces under variable pressures. A complex experiment was executed using an Al alloy sheet, tool contact surfaces of different roughness, two types of lubricants and variable contact pressures. The procedure included analytically pre-defined contact pressure functions based on which, for each of the mentioned conditions, the dependencies of the drawing forces and friction coefficients were obtained. The pressure in function P1 constantly decreased from a high initial value until the minimum, while in function P3 the pressure increased until the minimum value at the halfway point of the stroke, after which it increased up to the initial value. On the other hand, the pressure in function P2 constantly increased from the initial minimum value until the maximum value, while in function P4 the pressure increased until reaching the maximum value at the halfway point of the stroke, after which it decreased to the minimum value. This enabled the determination of the influence of tribological factors on the process parameters of intensity of traction (deformation force) and coefficient of friction. The pressure functions starting with decreasing trends produced higher values for the traction forces and the friction coefficient. In addition, it was established that the roughness of the contact surfaces of the tool, especially those with titanium nitride coating, has a significant influence on the process parameters. For surfaces of lower roughness (polished), a tendency of the Al thin sheet to form a glued-on layer was noticed. This was especially prominent for lubrication with MoS₂-based grease under conditions of high contact pressure (functions P1 and P4 at the beginning of the contact).