Static and Free Vibration Analyses of Functionally Graded Plane Structures

: In recent years, the use of functionally graded materials has been the focus of several studies due to their intrinsic ability to be tailored according to the requirements of structures while minimising abrupt stress transitions commonly found in laminated composites. In most studies, the materials’ mixture gradient is established through a structural component, i.e., thickness, which is known to visibly enhance structural behaviour. However, depending on the type of structure, it is important to exploit the possibility of building a structure using other gradient directions. The innovative characteristic of this work, which aims to study plane truss and frame-type structures made of functionally graded materials, lies in the speciﬁcity that the materials’ mixture gradient occurs as a function of a geometric structure feature, i


Introduction
In engineering applications, conflicting property requirements necessitate combining materials, which can be in the molten state (or alloying), limited by thermodynamic equilibrium limits, and in the solid state, resulting in composite materials.Composite materials can be divided into laminates and functional gradients, with the latter being characterised by the minimisation of problems that occur in laminates, namely sharp interfaces between layers in which abrupt stress transitions occur and failures start [1,2].
Functionally graded materials are composite materials with properties that result from a combination of distinct materials in which the mixture proportions vary along one or more spatial coordinates.The improved properties of FGMs that have been achieved by using this tailoring flexibility in parallel with technological advances in manufacturing such as additive manufacturing, powder metallurgy, etc., (Kumar et al. [3]) justify their increasingly frequent use.FGMs have wide applicability in aircraft, spacecraft, and space structures, which space exploration depends on, as well as in car parts, electronic and medical equipment, thermal coatings for ceramic engines, the marine industry, gas turbines, and biomedicine, among other applications [4][5][6][7][8].Although functionally graded materials still have high production costs, this type of material is a topic that has been significantly discussed by numerous authors in recent years, with many publications on the subject being published, making it a very current and important topic that should be explored further [1,9,10].
To investigate the advantages of using FGMs, several authors have developed models through which static and free vibration analyses have been carried out to investigate the influence of the materials' characteristic parameters on their behaviour.In this sense, Chi and Chung [11] investigated elastic, rectangular, and simply supported FGM plates of medium thickness that were subjected to transverse loading.They considered constant Poisson's ratios and varying the elasticity modulus through the plate's thickness according to a power law, i.e., sigmoid, or exponential volume fraction function.The series solutions were obtained based on the classical plate theory and Fourier series expansion.The closedform solutions illustrated by the Fourier series expression were given in this paper, and the closed-form and finite element solutions were compared and discussed in the second part of the study (Chi and Chung [12]).
Reddy et al. [13] published a review paper about the governing equations and analytical solutions of the classical and shear deformation theories of FGM beams.More specifically, the classical, first-order, and third-order shear deformation theories accounted for the through thickness variation in a dual-phase FGM and modified couple stress (i.e., strain gradient) and von Kármán nonlinearity.This study also included analytical solutions for the bending of the linear theories to illustrate the influence of material variation, boundary conditions, and loads.One of the major issues when manufacturing FGMs is the presence of porosities.Considering this, Jouneghani et al. [14] investigated the vibrational behaviour of doubly curved shells made of FGMs, including porosities.The first-order shear deformation theory was taken as the theoretical framework, considering that the FGMs' mechanical properties would vary through the thickness direction according to two power laws of volume fraction distribution.The strain components were established in an orthogonal curvilinear coordinate system, and the governing equations were derived according to Hamilton's principle, and then Navier's solution method was used; the numerical results concerning three different types of shell structures were presented.Pham et al. [15] proposed the notion of combining the edge-based smoothed finite element method and mixed interpolation of tensorial components' triangular element to carry out a dynamic analysis of sandwich plates with a functionally graded porous core subjected to blast load.The configuration of the sandwich plates considered a homogeneous metal bottom layer, an all-ceramic top layer, and a FGP core layer with uneven porosity distribution.The proposed element aimed to enhance the accuracy and convergence of the MITC3 element as well as the traditional triangular element.The authors performed some studies to confirm the expected performance of the proposed method and other numerical studies to characterize the vibration behaviour of the sandwich plates.
Chakraborty et al. [16] verified that the presence of FGM layers in beams produced significantly different behaviours when compared with other beams constituted by only one of the constituent materials due to the combination of their properties.Koutoati et al. [17] studied the stresses and strains in the interlayer interfaces of a sandwich component, and they have found that the use of FGM allows for a reduction in the presence of residual stresses, thus minimizing the occurrence of delamination and microcracking at the interfaces.In a study by Soliman et al. [18], isotropic and orthotropic beams with functional gradient were studied, considering the variation in properties through the thickness.The results showed that the FGM beams with different power law exponents have greater resistance to distributed loads and concentrated loads when compared to isotropic and orthotropic beams for different boundary conditions and different lengths.Wang et al. [19] investigated the dynamic behaviour of a pinned-pinned spinning exponentially functionally graded shaft with unbalanced loads.For this, they used the Rayleigh beam model, considering rotary inertia and gyroscopic effects.The model validity was confirmed, and a numerical parametric study was carried out to assess the influence of the main parameters.The results achieved led to the conclusion that the vibration and instability of the spinning shaft strongly depended on the unbalanced load and material gradient.
In a study by Zohra et al. [20], the critical buckling temperatures and natural frequencies of beams with a thickness functional gradient were investigated, and it was observed that the effects of shear deformation on frequencies tend to be more significant when the beams become shorter (thicker).These effects were more evident for higher mode frequencies.In a study by Singh et al. [21], it was stated that most FGMs presented material properties variation through the thickness direction and also in the axial direction, although only to a minor extent.In addition, the authors noted that the selection of a suitable material for the intended application was an immediate and direct challenge for the future development of technology in this field of research.In a study by Cao et al. [22], there were studied beams with uniform and non-uniform cross-sections with an axially functionally graded material for different boundary conditions.In Maalawi's study [23], an optimisation model for improving the performance of different types of structural composites was presented, and the FGM concept was also taken into account.Several scenarios were considered to model the spatial variation in these materials' properties.The proposed optimisation strategies included maximising natural frequencies in thin composite beams, optimising transmission shafts against torsional buckling and rotational instability, and maximising the critical flight speed of subsonic aircraft wings.Wu et al. [24] developed a mixed finite element for the nonlinear free vibration analysis of FGM beams based on the Timoshenko beam theory under combinations of simply supported, free, and clamped edge conditions.Additionally, the finite element considered the von Kármán geometrical nonlinearity.The material properties of the beam varied through the thickness according to power-law volume fraction distributions, and the effective material properties were estimated via the rule of mixtures.A multilayer perceptron back propagation neural network was also developed to predict the nonlinear free vibration behaviour of the FGM beam.
In Katili et al.'s study [25], it was found that, in functionally graded beams with varying properties through the thickness, the evolution of the modulus of elasticity has the greatest influence on the beam's stress and displacement distributions.Thus, choosing the power law exponent adequately allowed the material properties to be modelled to comply with the requirements regarding the minimisation of stresses and displacements in beams.Alshorbagy et al. [26] analysed the free vibration characteristics and dynamic behaviour of beams with variation in properties along the thickness and axial directions using the finite element method.It was verified that the variation in the material distribution along the axial direction influenced the variation in the beam stiffness along its length, hence affecting the frequencies and the respective vibration modes.However, the vibration modes were not affected by the variation in properties in thickness.This study also found that the natural frequencies increased with increasing power law exponent when the ratio between the modulus of elasticity was less than 1 and decreased otherwise.
Murin et al. [27] presented a homogenised FGM beam finite element for modal analysis, considering a double symmetric cross-section.The element considered the shear deformation effect and the effect of longitudinally varying inertia and rotary inertia as well as a longitudinally varying Winkler elastic foundation and the influence of internal axial forces.The material effective quantities assumed a continuous variation and were determined using extended mixture rules and the multilayer method.The authors performed numerical simulations via conducting modal analyses of single FGM beams and spatial beam structures.Nguyen et al. [28] considered the first-order shear deformation theory to carry out static and free vibration analyses of axially loaded rectangular FGM beams.Improved transverse shear stiffness was derived from the in-plane stress and equilibrium equation; thus, the associated shear correction factor was obtained analytically.The authors presented analytical solutions for simply supported FGM beams, and the results were compared with existing solutions.The influence of the power-law exponent, the dissimilarity of the constituent materials, and the Poisson's ratio was investigated through different mechanical responses.A finite element model based on the first-order shear deformation theory for free vibration and the buckling of functionally graded beams was proposed by Kahya and Turan [29].The material properties varied continuously through the beam thickness according to the power-law expression, and the governing equations were derived via Lagrange's equations.Numerical parametric studies for the natural frequencies and buckling loads were performed for different boundary conditions, power-law exponents, and span/depth ratios.Banerjee and Ananthapuvirajah [30] studied the free vibration behaviour of FGM beams (FGBs) and frameworks containing FGBs by using the dynamic stiffness method.
In their study, it was considered that the material properties of the FGBs would vary continuously through the thickness according to a power law function.The numerical results were validated against alternative published results, and as mentioned by the authors, in the absence of published results for frameworks containing FGBs, consistency checks on the reliability of results were performed.
From our review of the literature, which has been briefly summarized in the above paragraphs, it was possible to understand that there are many studies about the behaviour of beams, plates, and even shells in FGMs; however, it was not possible to find published investigations on the behaviour of plane truss and/or frame-type structures, where the material gradient is associated with a structure dimensional characteristic, such as its height or its width.
One considers that this structure-dependent material gradient approach can be viewed as an additional design variable when dealing with such types of structures, so this was the main motivation to adopt a material gradient evolution that depends on the structure's vertical coordinate.
The first section of the present study considers some aspects related to the development of the proposed finite element model.This is followed by a section where the model is verified for some isotropic homogeneous and FGM beams.The results achieved are then compared against alternative solutions.Afterwards, several plane truss and frame structures, are analysed and discussed.
To the authors' knowledge, aside from the set of case studies made available for possible reproduction, there are no published works that consider this paper's proposed approach, and this constitutes the main and innovative contribution of the present work.

Functionally Graded Materials
The prediction of the macroscopic properties of an FGM can be modelled according to different homogenisation models [31], perhaps the most used of which is the Rule of Mixtures.According to this rule, if one has a non-porous FGM composed by material 1 and material 2, a generic property P FGM , can be predicted by Equation (1), where V f corresponds to the volume fraction of material 1 (in the present work, a metallic or ceramic material).The volume fraction can be defined as a function of different space coordinates.Figure 1 shows the evolution of the volume fraction for different power law exponents.
As the focus of this work is the study of 2D structures with functional gradients relative to the vertical coordinate of a structure, a generic structure in which this variation in properties occurs is presented in Figure 3.In most published works that focus on beam structures, the volume fraction is expressed in terms of the thickness or the length of a beam.For illustrative purposes, Figure 2a and Equation (2) represent a through thickness mixture variation (direction z) that considers the typical power law.For this case, considering H as the maximum height of the structure, the volume fraction can be generically written as follows: Figure 2b and Equation (3) illustrate a material's mixture along the x direction that is associated with the length of a beam and its volume fraction law [32].

Displacement Field and Equilibrium Equations
As the focus of this work is the study of 2D structures with functional gradients relative to the vertical coordinate of a structure, a generic structure in which this variation in properties occurs is presented in Figure 3.
As the focus of this work is the study of 2D structures with functional gradients relative to the vertical coordinate of a structure, a generic structure in which this variation in properties occurs is presented in Figure 3.For this case, considering H as the maximum height of the structure, the volume fraction can be generically written as follows: For this case, considering H as the maximum height of the structure, the volume fraction can be generically written as follows:

Displacement Field and Equilibrium Equations
To carry out this study, we used the First-order Shear Deformation Theory (FSDT), and the displacement field for a beam under bending in the xy plane can be written as follows [32,33]: The beam schematic's representation and its relation to a local coordinate system, as well as its degrees of freedom, can be observed in Figure 4.
To carry out this study, we used the First-order Shear Deformation Theory (FSDT), and the displacement field for a beam under bending in the xy plane can be written as follows [32,33]: (, , ) (, ) =  (, )  (, ) +  θ (, ) 0 (5) The beam schematic's representation and its relation to a local coordinate system, as well as its degrees of freedom, can be observed in Figure 4.The degrees of freedom,  and  , are the mid-plane displacements along the x and z directions, which correspond to the length and thickness directions of the beam.θ is the rotation of the beam's mid-plane around the direction z.
The corresponding constitutive relations are in a compact form and can be written as follows [32,33]: where the generalized strain vectors, considering small deformations, are as follows: with the generalized displacements vector  consisting of the degrees of freedom  ,  , and θ , and the membrane-bending and shear elastic stiffness coefficients matrices are given as Considering the aim of the present work-focusing on the static and free vibration analyses of FGM 2D truss and frame-type structures-we considered the Hamilton principle to derive the equilibrium equations.Generically, this principle's formulation is given as shown in Equation (10) [34]: where, for a generic finite element (e), the respective values of the kinetic energy , the elastic strain energy , and the work performed by the external applied surface forces  are given as follows: The degrees of freedom, u 0 and v 0 , are the mid-plane displacements along the x and z directions, which correspond to the length and thickness directions of the beam.θ 0 z is the rotation of the beam's mid-plane around the direction z.
The corresponding constitutive relations are in a compact form and can be written as follows [32,33]: where the generalized strain vectors, considering small deformations, are as follows: with the generalized displacements vector {q} consisting of the degrees of freedom u 0 , v 0 , and θ z , and the membrane-bending and shear elastic stiffness coefficients matrices are given as Considering the aim of the present work-focusing on the static and free vibration analyses of FGM 2D truss and frame-type structures-we considered the Hamilton principle to derive the equilibrium equations.Generically, this principle's formulation is given as shown in Equation (10) [34]: where, for a generic finite element (e), the respective values of the kinetic energy T, the elastic strain energy U, and the work performed by the external applied surface forces W are given as follows: After some mathematical manipulation, one can yield the corresponding equilibrium equation at the element level, which, in a compact written format, can be written as follows: For the linear static analysis, the contribution of the kinetic energy is neglected; thus, the equation is simplified, leading to the following: Otherwise, to perform free vibration analysis, only the load vector will be null.Assuming free harmonic vibrations, the equilibrium equation will become Equation ( 14) [31,34,35]: where K e , M e , and {F e } represent the element stiffness, mass matrices, and element forces vector, respectively; q e and .. q e represent the element nodal degrees of freedom and corresponding accelerations; and ω denotes the j-th natural frequency.
To implement these equilibrium equations in the discretized model context, we used Lagrange quadratic beam elements with three degrees of freedom at a generic k-th node-the axial (u k ) and transverse (v k ) displacement and the rotation on the plane xy (θ zk ), k = 1,3-associated with the midplane surface of the beam [34].Figure 5 represents the one-dimensional quadratic element on its own reference in natural coordinates.

𝑊 = ( 𝑞 𝑓 ) 𝑑𝑆𝑒
After some mathematical manipulation, one can yield the corresponding equilibrium equation at the element level, which, in a compact written format, can be written as follows: +   =  (12) For the linear static analysis, the contribution of the kinetic energy is neglected; thus, the equation is simplified, leading to the following: Otherwise, to perform free vibration analysis, only the load vector will be null.Assuming free harmonic vibrations, the equilibrium equation will become Equation ( 14) [31,34,35]: where  ,  , and  represent the element stiffness, mass matrices, and element forces vector, respectively;  and  represent the element nodal degrees of freedom and corresponding accelerations; and  denotes the j-th natural frequency.
To implement these equilibrium equations in the discretized model context, we used Lagrange quadratic beam elements with three degrees of freedom at a generic k-th nodethe axial ( ) and transverse ( ) displacement and the rotation on the plane xy (θ ), k = 1,3-associated with the midplane surface of the beam [34].Figure 5 represents the onedimensional quadratic element on its own reference in natural coordinates.Considering that a generic degree of freedom is approximated as a linear combination of the corresponding nodal degrees of freedom and the shape functions, the previous element stiffness   and mass   matrices can be represented in a compact form as follows [34,36]: where   and   -the matrices that pair the deformations with the displacementsand -the matrix of the Lagrange quadratic shape functions  that are associated with the displacement field description-are, for the i-th node, constituted as follows: Considering that a generic degree of freedom is approximated as a linear combination of the corresponding nodal degrees of freedom and the shape functions, the previous element stiffness K e and mass M e matrices can be represented in a compact form as follows [34,36]: where B mb and B s -the matrices that pair the deformations with the displacements-and L-the matrix of the Lagrange quadratic shape functions N i that are associated with the displacement field description-are, for the i-th node, constituted as follows: Considering that the volume fraction expression depends on the structure height, it must be first rewritten within the domain of each element.Hence, the volume fraction expression can be rewritten as in Equation (18), where, in the structure referential, Y k,1 and Y k,3 correspond to the height coordinates of the initial and final nodes of the kth element, respectively, and ξ is the natural coordinate of the element local referential.
After the calculation of all elements' matrices and vectors, some may need to be transformed from their local coordinate system to the structure one, the global coordinate system, and when this happens, one should consider the usual transformation procedure [34]: Here, T is, in this case, the transformation matrix from the element to the structure coordinate system.The equilibrium equations' systems for the whole discretized domain are then obtained by assembling the element matrices and vector contributions, after which the boundary conditions are imposed, and the analyses can be performed.

Results
In the following sub-sections, we present some verification studies and, subsequently, a set of case studies to assess the influence of parameters that characterize the FGM plane structures.

Verification Studies
For each verification case presented in this section, the deviations in the obtained results compared to the ones from the literature were determined through using the following generic expression:

Static Analysis of a Simply Supported Functionally Graded Beam
In the first verification case, we considered the behaviour of a simply supported functionally graded beam with an aspect ratio of L/h = 20, where L represents beam length, and h represents section height (also studied by Nam et al. [37]).The beam is made of aluminium (E m = 70 GPa, ν = 0.3) and alumina (E c = 380 GPa, ν = 0.3), and the non-dimensional transverse displacement can be obtained using the following multiplier: For this purpose, a square cross-section with a unit edge and a twenty-fold higher length was assumed.The beam is subjected to a uniformly distributed load (q) of intensity of 100 N/m.The beam properties vary through the thickness according to the power law, and the rule of mixtures, which allows for the expression of this variation profile, is written as follows: From Table 1, good agreement between the results obtained by the present model and the ones in [37] can be observed for all of the power law exponents.The second verification case consists of a free vibration analysis of a simply supported functionally graded beam, considering a volume fraction distribution that varies both through the thickness and the length of the beam.Alshorbagy et al. [26] tested different aspect ratio relations and boundary conditions, as well as different power law exponent values.For the present work, we compared the present model's results with those of [26] for an aspect ratio of L/h = 20.
For the case of the functional gradient beam with variation in thickness, the rule of mixtures represents the beam being purely metallic on the bottom surface (L subscript) and purely ceramic at the top surface (U subscript), as given by the following: For the case of the functional gradient along the length, the rule of mixtures represents the beam being purely ceramic at the left extremity and purely metallic at the right extremity, according to the following expression: The results obtained, have been adimensionalised, through the following relations: The first fundamental frequencies for the functionally graded beam with variational properties through the thickness are given for the different elasticity moduli ratios in Table 2.The first fundamental frequencies for the functionally graded beams with length-based variational properties are presented in Table 3. From Tables 2 and 3, good agreement between the results obtained through the implemented model and the ones from [26] can be observed for each material distribution (both for thickness and length), modulus of elasticity ratio, and power law exponent.

Static and Free Vibration Analysis of an Isotropic Frame-Type Structure
The isotropic frame-type structure presented in Figure 6a is discretized into six beam elements in which each element has three nodes, as can be observed in Figure 6b.
For the static analysis of the structure, the elements were considered to be made of a metallic material with an elasticity modulus of 200 GPa and a square cross-section with a 0.02 m edge.The length L was set to 1 m, and the magnitudes of loads were set to Px = 100 N, Py = 1000 N, M0 = 250 N.m, and q = 90 N/m.
From Tables 2 and 3, good agreement between the results obtained through the implemented model and the ones from [26] can be observed for each material distribution (both for thickness and length), modulus of elasticity ratio, and power law exponent.

Static and Free Vibration Analysis of an Isotropic Frame-Type Structure
The isotropic frame-type structure presented in Figure 6a is discretized into six beam elements in which each element has three nodes, as can be observed in Figure 6b.For the static analysis of the structure, the elements were considered to be made of a metallic material with an elasticity modulus of 200 GPa and a square cross-section with a 0.02 m edge.The length L was set to 1 m, and the magnitudes of loads were set to Px = 100 N, Py = 1000 N, M0 = 250 N.m, and q = 90 N/m.Different discretisations were considered by dividing the members (I, II, III, IV, V, VI) of the structure in such a way that the lengths of the elements were close to each other.The support reactions obtained are presented in Tables 4 and 5.  Different discretisations were considered by dividing the members (I, II, III, IV, V, VI) of the structure in such a way that the lengths of the elements were close to each other.The support reactions obtained are presented in Tables 4 and 5. To verify the results obtained using the present model, we considered the analytical solution provided through Castigliano's second theorem and a numerical solution obtained by using an open-source finite element application named Ftool [38] (Table 6).As it is possible to conclude from Tables 4-6, the results obtained by the present model are reliable, based on the convergence of the results for different discretisation, where variations are very small.

Case Studies
The present section examines a set of case studies that focus on truss and frame-type structures made of functionally graded materials whose constituent materials' mixture varies as a function of the height of the structure.It was considered that the thickness of the different members would be much less than their length and that the volume fraction would be associated with the beam's mid-surface.In the latter assumption, this means that when the thickness coincides with the direction of the height of the structure, the volume fraction remains constant.
As the objective of this study is to focus on the linear static and free vibration behaviour of these structures, buckling was not considered at this stage, although we do plan to consider it in future studies.The studies were restricted to aspect ratios within the application domain of the first order shear deformation theory; thus, we restricted the aspect ratio so that it varied between 10 and 20.The first case study that we will examine is a simple truss consisting of three members and subject to a concentrated load P at the free node, represented in Figure 7, for which the nodal displacements and the first five natural frequencies of the structure are determined for different values of the power law exponents.These exponents determine the variation profile of the quantity of a given material (metallic or ceramic) within the mixture along the vertical coordinate (Y) of the structure.In this case, we considered lengths of 0.2 m for the horizontal and vertical members, with a square cross-section of 20 mm edge, so that the aspect ratio varied between 10 and 15.The coordinate system used can be observed in Figure 8; based on this, the expression of the metal volume fraction was defined as follows:  The first parametric study performed involved an FGM structure (the volume fraction power law of which is described in Equation ( 29)) that was associated with the variation of the metallic material as a function of the structure's vertical coordinate.The nodal displacements are presented in Table 7, and the first five natural frequencies are shown in Table 8.
The properties' evolution (Figure 9a) along the height can be described as follows:   The first parametric study performed involved an FGM structure (the volume fraction power law of which is described in Equation ( 29)) that was associated with the variation of the metallic material as a function of the structure's vertical coordinate.The nodal displacements are presented in Table 7, and the first five natural frequencies are shown in Table 8.The properties' evolution (Figure 9a) along the height can be described as follows: The first parametric study performed involved an FGM structure (the volume fraction power law of which is described in Equation ( 29)) that was associated with the variation of the metallic material as a function of the structure's vertical coordinate.The nodal displacements are presented in Table 7, and the first five natural frequencies are shown in Table 8.
The properties' evolution (Figure 9a) along the height can be described as follows:  The results obtained for the displacements of the node under study allowed us to verify that there was a decrease in the vertical displacement with the increase in the power law exponent due to the increasing quantity of material with a higher modulus of elasticity (in this case, ceramic).
On the other hand, if the power law is considered to characterize the distribution of the ceramic material as a function of the vertical coordinate of the structure, a behavioural trend in the opposite sense is verified, as expected.These results are presented in Tables 9 and 10 for the nodal displacements and natural frequencies, respectively.The properties' generic distribution (Figure 9b) in this case is now characterized by the following: In this case, the nodal displacements presented for a situation in which the power law rules the distribution of the ceramic material along the structure's height exhibits an antithetical behaviour compared to the previous situation (power law relative to the distribution of the metallic material as a function of the vertical coordinate), meaning that the smaller the amount of ceramic material, the larger are the nodal displacements obtained.From Tables 7-10, several trends were observed, namely, the displacements more directly affected by the load to which Truss 1 is subjected to and a decrease with an increase in the ceramic phase with a higher modulus of elasticity.
As for the displacement along the X direction, it can be observed that when the rule of mixtures is related to the distribution of aluminium (Table 7), it remains constant with the displacement verified to the totally metallic structure, except for the cases wherein the structure is totally ceramic.In the opposite case, when the rule of mixtures is related to the distribution of zirconia (Table 9), it remains constant with the displacement verified for the totally ceramic structure regardless of the power law exponents, except when the structure is totally metallic.
From Table 8, in which the first five natural frequencies of the structure are presented, it is possible to observe that when the volume fraction is associated with the metallic material, the first natural frequency increases from ey = 0 (totally metallic) to ey = 1 and decreases from ey = 1 to ey = 10.For the remaining frequencies, this increase occurs from ey = 0 to ey = 10.
In the case of the volume fraction being directly related to the distribution of the ceramic material (see Table 10 for the natural frequencies), one can observe a decrease in the fundamental frequency from ey = 0 (all-ceramic structure) to ey = 1, followed by an increase from ey = 1 to ey = 10.For the higher-order frequencies, a decrease occurs from ey = 0 to ey = 10.As expected, the lowest results are verified when the structure is totally metallic, and the highest ones are verified when Truss 1 is totally ceramic.

Structure 2
In this second case study, we analyse the truss shown in Figure 10 (which consists of nineteen members and is subjected to a concentrated load P at its upper vertex).The displacement of this point and the first five natural frequencies of the structure are presented for different values of the volume fraction law exponent.
In this case, the distance between the supports was considered to be 1.4 m, and the structure's height was considered to be half of that distance; the elements have a square section with a 20 mm side, and the aspect ratio varies between 10 and 20.The origin of the structure coordinate system is located at the left support, as depicted in Figure 11, so the expression of the metal volume fraction is defined by the following: Figure 11 schematically presents the minimum discretisation of the structure with 19 quadratic elements and the evolution of the volume fraction as a function of the vertical Y coordinate.
structure coordinate system is located at the left support, as depicted in Figure 11, so the expression of the metal volume fraction is defined by the following: Figure 11 schematically presents the minimum discretisation of the structure with 19 quadratic elements and the evolution of the volume fraction as a function of the vertical Y coordinate.The first parametric study performed considers that the volume fraction described in Equation ( 32) is directly associated with the variation of the metallic material in the structure as a function of the vertical coordinate of the structure (Figure 12a).The nodal displacements are presented in Table 11, and the natural frequencies are shown in Table 12.
(a) (b) . Schematic representation of Truss 2 with load P applied to its highest point.
structure coordinate system is located at the left support, as depicted in Figure 11, so expression of the metal volume fraction is defined by the following:   The first parametric study performed considers that the volume fraction describe Equation ( 32) is directly associated with the variation of the metallic material in structure as a function of the vertical coordinate of the structure (Figure 12a).The no displacements are presented in Table 11, and the natural frequencies are shown in T 12.The first parametric study performed considers that the volume fraction described in Equation ( 32) is directly associated with the variation of the metallic material in the structure as a function of the vertical coordinate of the structure (Figure 12a).The nodal displacements are presented in Table 11, and the natural frequencies are shown in Table 12.
structure coordinate system is located at the left support, as depicted in Figure 11, so the expression of the metal volume fraction is defined by the following:   The first parametric study performed considers that the volume fraction described in Equation ( 32) is directly associated with the variation of the metallic material in the structure as a function of the vertical coordinate of the structure (Figure 12a).The nodal displacements are presented in Table 11, and the natural frequencies are shown in Table 12.The results obtained for the nodal displacements allow us to verify that there is a decrease in the vertical displacement of the selected highest node with the increase in the power law exponent due to the associated increasing quantity of material with a higher modulus of elasticity (in this case, the ceramic phase).On the other hand, if it is considered that the power law characterizes the distribution of the ceramic material as a function of the vertical coordinate (Figure 12b), an inverse behaviour to the one considered in the previous study is verified.The results of nodal displacements and natural frequencies are presented in Tables 13 and 14, respectively, for this configuration.As can be observed in Table 13, where the volume fraction is directly related to the distribution of the ceramic material along the structure' height, one can conclude that, as expected, the smaller the amount of ceramic material whose modulus of elasticity is higher than that of the metallic material, the larger are the nodal displacements obtained.
From Tables 11-14, several trends can be observed.In concrete, one can observe that the displacements more visibly affected by the load to which Truss 2 is subjected to (displacement along y direction) decrease with the increase in the ceramic phase.As for the displacement along the X direction and when the volume fraction is directly related to the distribution of aluminium (Table 11), this displacement remains equal to the one presented from the totally ceramic structure, regardless of the power law exponent.In the opposite case, when the volume fraction is related to the distribution of zirconia (Table 13), the displacement remains equal to the one presented by the metallic structure despite the power law exponent values.As for the rotation, there were no clear trends to highlight.
When the volume fraction characterizes the distribution of the metallic material phase, it is also possible to conclude from Table 12 that the first natural frequency decreases from ey = 0.1 to ey = 10; for the fourth frequency, the maximum value is obtained when ey = 5, and for the remaining ones, the maximum value is obtained when ey = 10.The highest value for the first frequency is obtained when this truss has a power law exponent that is equal to 0.1, and for the remaining frequencies, the highest values are obtained when truss 2 is totally ceramic.As expected, this structure presents the lowest fundamental frequency when it is fully metallic.
In the opposite situation, when the power law volume fraction characterises the distribution of the ceramic material, the results of the natural frequencies depicted in Table 14 show an increasing trend from ey = 0.1 to ey = 1 and a decreasing one from ey = 1 to ey = 10.For the remaining ones, there is a decrease from ey = 0.1 to ey = ∞.Also, as expected, the lowest frequencies are obtained when the structure is completely metallic.The highest values for the first, third, and fifth frequencies are obtained when truss 2 is totally ceramic, but for the remaining ones, the highest values are obtained when ey = 0.1.

Structure 3
The third case study is a simple frame consisting of three members that is subjected to a distributed load in the transverse direction of the upper horizontal member, which is additionally subjected to a concentrated load P along its length, as presented in Figure 13.
The upper nodal displacements and the first five natural frequencies of the structure are presented in Tables 15-18, as for the previous case studies.The upper nodal displacements and the first five natural frequencies of the structure are presented in Tables 15-18, as for the previous case studies.The characteristics of the cross-section are equal to those of the previous case; the distance between the supports is 0.2 m, and the height of the structure was set so that it could be double that distance.The origin of the coordinate system coincides with the lefthand side support.The metal volume fraction is in this case defined as follows:  The characteristics of the cross-section are equal to those of the previous case; the distance between the supports is 0.2 m, and the height of the structure was set so that it could be double that distance.The origin of the coordinate system coincides with the left-hand side support.The metal volume fraction is in this case defined as follows: Figure 14 schematically presents the minimum discretisation of the structure with 19 quadratic beam elements and the evolution of the volume fraction as a function of the vertical Y coordinate.The upper nodal displacements and the first five natural frequencies of the structure are presented in Tables 15-18, as for the previous case studies.The characteristics of the cross-section are equal to those of the previous case; the distance between the supports is 0.2 m, and the height of the structure was set so that it could be double that distance.The origin of the coordinate system coincides with the lefthand side support.The metal volume fraction is in this case defined as follows: In the first parametric study, we considered the volume fraction from Equation (33) to be directly associated with the variation in the metallic material in the structure as a function of the vertical coordinate of the structure (Figure 15a).The nodal displacements achieved are presented in Table 15, and the natural frequencies can be observed in Table 16.In the first parametric study, we considered the volume fraction from Equation (33) to be directly associated with the variation in the metallic material in the structure as a function of the vertical coordinate of the structure (Figure 15a).The nodal displacements achieved are presented in Table 15, and the natural frequencies can be observed in Table 16.The results obtained for the nodal displacements allowed us to conclude that there is a decrease in the vertical displacement of the nodes under study with an increase in the power law exponent due to the increasing quantity of the material with a higher modulus of elasticity (in this case, ceramic) (Figure 15b).The results obtained for the nodal displacements allowed us to conclude that there is a decrease in the vertical displacement of the nodes under study with an increase in the power law exponent due to the increasing quantity of the material with a higher modulus of elasticity (in this case, ceramic) (Figure 15b).
Again, if the volume fraction characterizes the distribution of the ceramic material as a function of the vertical coordinate of the structure (Figure 15b), an expected opposite behaviour is shown (Tables 18 and 19).From Table 17, it is possible to see that the smaller the amount of ceramic material, whose modulus of elasticity, in the present case, is higher than that of the metallic material, the larger are the nodal displacements obtained.
From Tables 15-18, several trends can be observed, namely, that displacements decrease with the increase in the ceramic phase characterized by a higher modulus of elasticity when compared to the that of the metallic phase.
For the Frame 1 structure, the first five natural frequencies of which are presented in Table 16, one can conclude that when the power law volume fraction is associated with the metallic material, the first natural frequency increases from ey = 0 to ey = 2 and decreases from ey = 2 to ey = ∞.For the second frequency, the values increase from ey = 0 to ey = 1 and decrease from ey = 1 to ey = 10, being higher when ey = ∞ due to the complete ceramic character of the structure.For the remaining frequencies, one can observe that there is an increasing trend from ey = 0 to ey = ∞.
In the case of the power law characterising the distribution of the ceramic material, the results presented in Table 18 show that, for the first and second frequencies, there is a decreasing trend (from ey = 0 to ey = 10).The third frequency decreases from ey = 0 to ey = 1 and increases from ey = 1 to ey = 5 before decreasing.The remaining higher-order frequencies decrease from ey = 0 to ey = ∞.

Structure 4
The fourth case study focused on the frame-type structure represented in Figure 16, which is consists of 19 members and is subjected to a horizontal concentrated load P at the upper-left node and to a transverse uniformly distributed load in the upper-horizontal member.ceramic character of the structure.For the remaining frequencies, one can observe that there is an increasing trend from ey = 0 to ey = ∞.
In the case of the power law characterising the distribution of the ceramic material, the results presented in Table 18 show that, for the first and second frequencies, there is a decreasing trend (from ey = 0 to ey = 10).The third frequency decreases from ey = 0 to ey = 1 and increases from ey = 1 to ey = 5 before decreasing.The remaining higher-order frequencies decrease from ey = 0 to ey = ∞.

Structure 4
The fourth case study focused on the frame-type structure represented in Figure 16, which is consists of 19 members and is subjected to a horizontal concentrated load P at the upper-left node and to a transverse uniformly distributed load in the upper-horizontal member.
A parametric study concerning the influence of the power law exponent was performed to assess the upper nodal displacements and the first five natural frequencies of the structure.The distance between the structure's supports and the cross-section is equal to the previous case study.The structure's height is twice that distance.The metal volume fraction for this structure, considering the coordinate system origin shown in Figure 17, is defined by the following: Figure 17 schematically presents the minimum discretisation of the structure with 19 quadratic elements and the evolution of the volume fraction as a function of the vertical Y structure coordinate.
In Figure 18, one can observe a schematic representation of the graded material evolution, considering the writing of the rule of mixtures homogenisation scheme.
The results obtained for the nodal displacements (presented in Table 19) denote a decrease in the vertical displacement of the upper nodes with the increase in the volume fraction law exponent due to the increasing quantity of ceramic material (with a higher modulus of elasticity).In Table 20, the first five fundamental frequencies of Frame 2 are presented, with the volume fraction expression related to the metallic material phase.A parametric study concerning the influence of the power law exponent was performed to assess the upper nodal displacements and the first five natural frequencies of the structure.
The distance between the structure's supports and the cross-section is equal to the previous case study.The structure's height is twice that distance.The metal volume fraction for this structure, considering the coordinate system origin shown in Figure 17, is defined by the following: Figure 17 schematically presents the minimum discretisation of the structure with 19 quadratic elements and the evolution of the volume fraction as a function of the vertical Y structure coordinate.
In Figure 18, one can observe a schematic representation of the graded material evolution, considering the writing of the rule of mixtures homogenisation scheme.The results obtained for the nodal displacements (presented in Table 19) denote a decrease in the vertical displacement of the upper nodes with the increase in the volume fraction law exponent due to the increasing quantity of ceramic material (with a higher modulus of elasticity).In Table 20, the first five fundamental frequencies of Frame 2 are presented, with the volume fraction expression related to the metallic material phase.
Again, if one considers the homogenisation scheme written in the opposite direction, that is, the volume fraction law characterising the distribution of the ceramic material as a function of the vertical coordinate of the structure, an inverted behaviour is found as expected and as observed in Tables 21 and 22.
From Tables 19 and 21, the same trends that occurred for Frame 1 can be observed.It is clear that the existence of the intermediate horizontal beam on Frame 2 yields a significant reduction in the nodal displacements and a remarkable increase in the fundamental frequencies.From Table 20, the results obtained for the first five natural frequencies of the structure when the power law is relative to the metallic material show that, for the first, second, and fifth natural frequencies, they increase from ey = 0 to ey = 2 and decrease from ey = 2 to ey = 10.For the remaining frequencies, there is an increasing trend from ey = 0 to ey = 5 before they decrease from ey = 5 to ey = 10.When ey = ∞, which means the structure is fully ceramic, the highest frequencies are observed.
When the power law volume fraction characterises the distribution of the ceramic material (for which Table 22 shows all of the first five frequencies), except the for the fourth frequency, a decreasing trend can be observed from ey = 0 to ey = 5, and an increase from ey = 5 to ey = ∞ can also be observed.For the fourth frequency, there is a decrease from ey = 0 to ey = ∞.

Differently Oriented Cross-Section Areas with the Same Value
This last study was conducted to compare the structural response of the first frametype structure when different aspect ratio cross-sections are considered.For this study, we maintained the value of the cross-sectional area but changed the relationship between the sides.Initially, a square cross-section was considered, but we subsequently considered that one of the sides would have four times the dimension of the other.Figure 19 schematically represents the three sections analysed.This study enabled us to assess the influence that the second area moment has on the behaviour of the FGM structure.In this case, when the section is changed from A to B, the second moment increases four times; when the section is changed from A to C, the second moment decreases four times.
structure when the power law is relative to the metallic material show that, for the first, second, and fifth natural frequencies, they increase from ey = 0 to ey = 2 and decrease from ey = 2 to ey = 10.For the remaining frequencies, there is an increasing trend from ey = 0 to ey = 5 before they decrease from ey = 5 to ey = 10.When ey = ∞, which means the structure is fully ceramic, the highest frequencies are observed.
When the power law volume fraction characterises the distribution of the ceramic material (for which Table 22 shows all of the first five frequencies), except the for the fourth frequency, a decreasing trend can be observed from ey = 0 to ey = 5, and an increase from ey = 5 to ey = ∞ can also be observed.For the fourth frequency, there is a decrease from ey = 0 to ey = ∞.

Differently Oriented Cross-Section Areas with the Same Value
This last study was conducted to compare the structural response of the first frametype structure when different aspect ratio cross-sections are considered.For this study, we maintained the value of the cross-sectional area but changed the relationship between the sides.Initially, a square cross-section was considered, but we subsequently considered that one of the sides would have four times the dimension of the other.Figure 19 schematically represents the three sections analysed.This study enabled us to assess the influence that the second area moment has on the behaviour of the FGM structure.In this case, when the section is changed from A to B, the second moment increases four times; when the section is changed from A to C, the second moment decreases four times.A comparison between the most significant nodal displacements is presented in Table 23, where the displacement in the X direction and the rotation of the upper-left node are displayed for each section.It is important to note that, in this case, the aspect ratio for the section B situation varies between 5 and 10; for the section C situation, it varies between 20 and 40.
For each section considered, the displacements visibly vary depending on the exponent of the power law, confirming a significant influence of the second area moment.A comparison between the most significant nodal displacements is presented in Table 23, where the displacement in the X direction and the rotation of the upper-left node are displayed for each section.It is important to note that, in this case, the aspect ratio for the section B situation varies between 5 and 10; for the section C situation, it varies between 20 and 40.For each section considered, the displacements visibly vary depending on the exponent of the power law, confirming a significant influence of the second area moment.
From Table 23, it is clear that when the section is changed from A to B, which means that the section's height doubles and the width halves, the displacements evaluated decrease by a factor between 3.7 and 3.8.On the other hand, if the section is changed from A to C, which means that section's height halves and the width is doubled, the displacements evaluated increase by a factor of between 3.9 and 4.
The fundamental frequencies of this structure for each cross-section are next presented in Table 24.From Table 24, it can be observed that when the section is changed from A to B, the fundamental frequencies double; on the other hand, if the section is changed from A to C, they are halved.We verified the expected trends of both situations.Finally, the nodal generalized displacements of the central node of the horizontal element are presented in Tables 25-27, according to each section considered.From Table 25, one can observe the same behaviour already observed in Tables 26 and 27.When the section is altered from A to B, the displacements and rotation evaluated decrease significantly; if section A is altered to the C configuration, the displacements and rotation evaluated experience a significant increase.From Table 25, it is clear that the evaluated central node displacement decreases by a factor of approximately 3.8 if the section is altered from A to B, while if the section is altered from A to C, the evaluated displacement increases by a factor of approximately 3.95.From Table 26, one can observe that if the section is altered from the A to the B configuration, the evaluated displacement decreases by a factor between 1.5 and 1.86; if the section is altered to the C configuration, the displacement increases by a factor between 2.44 and 2.84.The rotation of the central node presented in Table 27 shows a decreasing trend by a factor between 4.6 and 5.4 when the section is altered from A to B. If section A is altered to C, the rotation increases a factor between 4.1 and 4.3.

Conclusions
From the results achieved, one can conclude that the base and the added material phases play an important role, particularly when the materials' dissimilarity is relevant, as is the case of the materials considered in the present study.In fact, if one considers the metallic phase as the base material and the ceramic as the added one, such inclusion results in improved structure performance.On the contrary, the mixture properties will worsen and degrade the structure's mechanical behaviour.
Overall, the nodal displacements show a decreasing trend with increasing ceramic phase quantity, which possesses a higher modulus of elasticity when compared to the metallic phase.Similarly, the natural frequencies present an opposite trend, although there is not a linear relation between frequency evolution and power law exponent progress.Specifically, in terms of the free vibration dynamic behaviour, if it is intended to maximize the fundamental frequency, it is visible that the best option for Truss 1, in terms of a functionally graded structure, is a power law exponent of 0.1, and the volume fraction is directly associated with the ceramic material.However, in this case, the all-ceramic structure enables a higher fundamental frequency.In the remaining cases, as symmetrical structures, it was noted that constructing the structures with homogeneous material is not the option that provides the best dynamic behaviour in the free vibration regime.For Truss 2, the best solution is a power law exponent of 0.1, but in this case, the volume fraction law describes the distribution of metallic material.For both Frame 1 and Frame 2, the best solution is a power law exponent of 2 and a volume fraction that describes the metallic material distribution.
The orientation of the cross-section significantly impacts the static and free vibration performance of the structures, which confirms our expectations.
Globally, the possibility of considering the variation in the materials' mixture ruled by structure-dependent volume fraction distributions can benefit the mechanical responses of structures and contribute to tailoring them toward specific operating conditions.These additional alternative or complementary solutions may enable improved free vibration and static performances, as shown by the case studies presented in this work.

Figure 1 .Figure 1 .
Figure 1.Volume fraction of material 1 as a function of a structure's height.

Figure 1 .Figure 2 .
Figure 1.Volume fraction of material 1 as a function of a structure's height.

Figure 3 .
Figure 3. Structure with variation in its properties as a function of its height.

Figure 2 .
Figure 2. Structure with variation in its properties (a) through the thickness and (b) along the length.

Figure 1 .
Figure 1.Volume fraction of material 1 as a function of a structure's height.

Figure 2 .
Figure 2. Structure with variation in its properties (a) through the thickness and (b) along the length.

Figure 3 .
Figure 3. Structure with variation in its properties as a function of its height.

Figure 3 .
Figure 3. Structure with variation in its properties as a function of its height.

Figure 4 .
Figure 4. Beam representation and local coordinate system.

Figure 4 .
Figure 4. Beam representation and local coordinate system.

Figure 5 .
Figure 5. One-dimensional quadratic element on its natural coordinate system.

2 Figure 5 .
Figure 5. One-dimensional quadratic element on its natural coordinate system.

Figure 6 .
Figure 6.Schematic representation of (a) an isotropic frame-type structure and the loadings considered in the verification study and (b) the structure's minimum discretisation.

Figure 6 .
Figure 6.Schematic representation of (a) an isotropic frame-type structure and the loadings considered in the verification study and (b) the structure's minimum discretisation.

Figure 7 .
Figure 7. Schematic representation of a three-member truss with tip vertical load.Figure 7. Schematic representation of a three-member truss with tip vertical load.

Figure 7 .
Figure 7. Schematic representation of a three-member truss with tip vertical load.Figure 7. Schematic representation of a three-member truss with tip vertical load.

Figure 7 .
Figure 7. Schematic representation of a three-member truss with tip vertical load.

Figure 8 .
Figure 8. Volume fraction evolution through the vertical coordinate of the truss.

Figure 9 .
Figure 9. Schematic Representation of Truss 1 (constituted by an aluminium-zirconia-graded material along the structure height).(a) Volume fraction associated with aluminium; (b) volume fraction associated with zirconia.

Figure 8 .
Figure 8. Volume fraction evolution through the vertical coordinate of the truss.

Figure 8
Figure 8 schematically presents the minimum discretisation of the structure with three quadratic elements and the evolution of the volume fraction as a function of the vertical Y coordinate.The first parametric study performed involved an FGM structure (the volume fraction power law of which is described in Equation (29)) that was associated with the variation of the metallic material as a function of the structure's vertical coordinate.The nodal displacements are presented in Table7, and the first five natural frequencies are shown in Table8.

Figure 8 .
Figure 8. Volume fraction evolution through the vertical coordinate of the truss.

Figure 9 .
Figure 9. Schematic Representation of Truss 1 (constituted by an aluminium-zirconia-graded material along the structure height).(a) Volume fraction associated with aluminium; (b) volume fraction associated with zirconia.

Figure 9 .
Figure 9. Schematic Representation of Truss 1 (constituted by an aluminium-zirconia-graded material along the structure height).(a) Volume fraction associated with aluminium; (b) volume fraction associated with zirconia.

Figure 10 .
Figure 10.Schematic representation of Truss 2 with load P applied to its highest point.

Figure 11 .
Figure 11.Representation of the variation in the volume fraction through the vertical coordinate of Truss 2.

Figure 11 schematically
Figure 11 schematically presents the minimum discretisation of the structure wit quadratic elements and the evolution of the volume fraction as a function of the vert Y coordinate.

Figure 10 .
Figure 10.Schematic representation of Truss 2 with load P applied to its highest point.

Figure 11 .
Figure 11.Representation of the variation in the volume fraction through the vertical coordina Truss 2.

Figure 11 .
Figure 11.Representation of the variation in the volume fraction through the vertical coordinate of Truss 2.

Figure 11 schematically
Figure 11 schematically presents the minimum discretisation of the structure with 19 quadratic elements and the evolution of the volume fraction as a function of the vertical Y coordinate.

Figure 10 .
Figure 10.Schematic representation of Truss 2 with load P applied to its highest point.

Figure 11 .
Figure 11.Representation of the variation in the volume fraction through the vertical coordinate of Truss 2.

Figure 12 .
Figure 12.Schematic representation of Truss 2. (a) Volume fraction characterizing the distribution of metallic material (in this case, aluminium) through the structure's height; (b) volume fraction characterizing the distribution of ceramic material (in this case, zirconia) through the structure's height.

Figure 13 .
Figure 13.Schematic representation of Frame 1 and the loading considered.

33 )Figure 14 Figure 13 .
Figure 14 schematically presents the minimum discretisation of the structure with 19 quadratic beam elements and the evolution of the volume fraction as a function of the vertical Y coordinate.

Figure 13 .
Figure 13.Schematic representation of Frame 1 and the loading considered.

Figure 14 Figure 14 .
Figure 14 schematically presents the minimum discretisation of the structure with 19 quadratic beam elements and the evolution of the volume fraction as a function of the vertical Y coordinate.

Figure 14 .
Figure 14.Representation of the variation in the volume fraction through the vertical coordinate of Frame1.

Figure 15 .
Figure 15.Schematic representation of Frame 1.(a) Volume fraction characterizes the distribution of metallic material (aluminium) through the structure's height; (b) volume fraction characterizes the distribution of ceramic material (zirconia) through the structure's height.

Figure 15 .
Figure 15.Schematic representation of Frame 1.(a) Volume fraction characterizes the distribution of metallic material (aluminium) through the structure's height; (b) volume fraction characterizes the distribution of ceramic material (zirconia) through the structure's height.

Figure 16 .
Figure 16.Schematic representation of Frame 2 and the loading considered.

Figure 16 .
Figure 16.Schematic representation of Frame 2 and the loading considered.

Figure 17 .Figure 18 .
Figure 17.Representation of the variation in the volume fraction through the vertical coordinate of the truss considered.

Figure 17 .
Figure 17.Representation of the variation in the volume fraction through the vertical coordinate of the truss considered.

Figure 17 .Figure 18 .
Figure 17.Representation of the variation in the volume fraction through the vertical coordinate of the truss considered.

Figure 18 .
Figure 18.Schematic representation of Frame 2 considering that the rule of mixtures is such that (a) the volume fraction characterizes the distribution of metallic material (in this case, aluminium) through the structure's height or (b) the volume fraction characterizes the distribution of ceramic material (in this case, zirconia) through the structure's height.

Figure 19 .
Figure 19.Cross-sections considered with the same cross-section area.

Figure 19 .
Figure 19.Cross-sections considered with the same cross-section area.
3.1.2.Free Vibration Analysis of a Simply Supported Functionally Graded Beam

Table 2 .
Non-dimensional fundamental frequencies for functionally graded beams with variational properties through the thickness.

Table 3 .
Non-dimensional fundamental frequencies for functionally graded beams with variation in properties along their length.

Table 4 .
Convergence study on the reactions in the clamped support at point A.

Table 5 .
Convergence study on the reactions in the clamped support at point B.

Table 4 .
Convergence study on the reactions in the clamped support at point A.

Table 5 .
Convergence study on the reactions in the clamped support at point B.

Table 6 .
Reactions in the boundary conditions of the frame-type structure obtained through using Castigliano's second theorem and the free-to-use Ftool application.

Table 7 .
Nodal displacements of the free node of the truss when the volume fraction characterizes the distribution of aluminium.

Table 8 .
First five fundamental frequencies of Truss 1 when the volume fraction characterizes the distribution of aluminium.

Table 9 .
Nodal displacements of the free node of the truss when the volume fraction characterizes the distribution of zirconia.

Table 10 .
First five fundamental frequencies of Truss 1 when the volume fraction characterizes the distribution of zirconia.

Table 11 .
Nodal displacements of Truss 2's highest node, where load P is Applied, when the volume fraction characterizes the distribution of aluminium.

Table 12 .
First five fundamental frequencies of Truss 2 when the volume fraction characterizes the distribution of aluminium.

Table 13 .
Nodal displacements of Truss 2's highest node, where load P is applied, when the volume fraction characterizes the distribution of zirconia.

Table 14 .
First five Fundamental frequencies of Truss 2 when the volume fraction characterizes the distribution of zirconia.

Table 15 .
Nodal displacements of Frame 1's upper nodes when the volume fraction characterizes the distribution of aluminium.

Table 16 .
First five fundamental frequencies of Frame 1 when the volume fraction characterizes the distribution of aluminium.

Table 17 .
Nodal displacements of Frame 1's upper nodes when the volume fraction characterizes the distribution of zirconia.

Table 18 .
First five fundamental frequencies of Frame 1 when the volume fraction characterizes the distribution of zirconia.

Table 15 .
Nodal displacements of Frame 1's upper nodes when the volume fraction characterizes the distribution of aluminium.

Table 16 .
First five fundamental frequencies of Frame 1 when the volume fraction characterizes the distribution of aluminium.

Table 19 .
Nodal displacements of Frame 2's upper nodes when the volume fraction characterizes the distribution of aluminium.

Table 19 .
Nodal displacements of Frame 2's upper nodes when the volume fraction characterizes the distribution of aluminium.

Table 19 .
Nodal displacements of Frame 2's upper nodes when the volume fraction characterizes the distribution of aluminium.

Table 20 .
First five fundamental frequencies of Frame 2 when the volume fraction characterizes the distribution of aluminium.

Table 21 .
Nodal displacements of Frame 2's upper nodes when the volume fraction characterizes the distribution of zirconia.

Table 22 .
First five fundamental frequencies of Frame 2 when the volume fraction characterizes the distribution of zirconia.

Table 23 .
Nodal displacements of Frame 1's upper-left node when the volume fraction characterizes the distribution of aluminium for each cross-section.

Table 24 .
Fundamental Frequencies for Frame 1 structures for each cross-section.

Table 25 .
Displacement according to the global x direction of the horizontal member's central node.

Table 26 .
Displacement according to the global y direction of the horizontal member's central node.

Table 27 .
Rotation of horizontal member's central node.