In a FG plate, mechanical and physical properties are varied continuously from one surface to another. You can approximate the variation using the elements of laminated (or multilayered) shell. In this case, you will have to specify the material properties of each layer according to the given law. This is similar to approximation of a continuous function by a piecewise constant function. The more layers you use the better the approximation. You can study convergence of the process increasing the number of layers. I guess, 10 or 20 layers would be sufficient.
For analysis of FGM problems by using a commercial finite element software, such as ABAQUS, assigning continuously variable properties (e.g. Young modulus, plasticity modulus, yield stress) across the wall thickness is not possible.
The virtual temperature technique that is mentioned in the manuscript is, indeed, a trick devised by the authors for modeling problem with continuously variable mechanical properties with excellent results.
In this technique, assigning continuously variable temperature across the thickness and then we create a one to one relationship between temperature and mechanical properties. So, Correlation of the distribution of temperature and mechanical properties of FGM material, allowed for the variation of the FGM mechanical properties in the thickness of cylinder.
Note that conductivity factor and other thermal parameters are zero, and temperature does not change in during the analysis time. indeed, the analysis is mechanical. In this technique, there is no need to subdivide the thick wall into several thin layers.
In this technique, assigning continuously variable temperature across the thickness and then we create a one to one relationship between temperature and mechanical properties. So, Correlation of the distribution of temperature and mechanical properties of FGM material, allowed for the variation of the FGM mechanical properties in the thickness of cylinder.
Note that conductivity factor and other thermal parameters are zero, and temperature does not change in during the analysis time. indeed, the analysis is mechanical. In this technique, there is no need to subdivide the thick wall into several thin layers.
You can divide your geometry to the tiny sections and attached the property that you want to each section. using an excel file to determine each section's property.