Modeling the Thermoelastic Bending of Ferric Oxide (Fe2O3) Nanoparticles-Enhanced RC Slabs

Nanoparticles, by virtue of their amorphous nature and high specific surface area, exhibit ideal pozzolanic activity which leads to the formation of additional C-S-H gel by reacting with calcium hydroxide, resulting in a denser matrix. The proportions of ferric oxide (Fe2O3), silicon dioxide (SiO2), and aluminum oxide (Al2O3) in the clay, which interact chemically with the calcium oxide (CaO) during the clinkering reactions, influence the final properties of the cement and, therefore, of the concrete. Through the phases of this article, a refined trigonometric shear deformation theory (RTSDT), taking into account transverse shear deformation effects, is presented for the thermoelastic bending analysis of concrete slabs reinforced with ferric oxide (Fe2O3) nanoparticles. Thermoelastic properties are generated using Eshelby’s model in order to determine the equivalent Young’s modulus and thermal expansion of the nano-reinforced concrete slab. For an extended use of this study, the concrete plate is subjected to various mechanical and thermal loads. The governing equations of equilibrium are obtained using the principle of virtual work and solved using Navier’s technique for simply supported plates. Numerical results are presented considering the effect of different variations such as volume percent of Fe2O3 nanoparticles, mechanical loads, thermal loads, and geometrical parameters on the thermoelastic bending of the plate. According to the results, the transverse displacement of concrete slabs subjected to mechanical loading and containing 30% nano-Fe2O3 was almost 45% lower than that of a slab without reinforcement, while the transverse displacement under thermal loadings increased by 10%.


Introduction
Nanotechnology has become one of the most important and exciting areas of many fields such as physics, chemistry, biology, and, very recently, civil engineering. It has also attracted a great deal of attention in the field of material science and technology. To contribute in the development of concrete physical and mechanical conventional properties, researchers intent to replace long-used pozzolanic powder additives in concrete composition, such as micro silica fume (MS), blast furnace slag, and fly ash, with new materials to achieve a more efficient and durable concrete while considering the economy aspects [1,2].
Intrinsically, the specific objective of researchers lately was to find the preference for how to use nanoparticles with amorphous form to add even more capabilities to a concrete matrix. The required improvement of concrete properties reinforced with nanoparticles is mostly to be achieved due to the unique properties of those nano-sized elements, such as their high strength, elevated modulus of elasticity, high specific surface area, exorbitant electrical conductivity, enhancement of calcium silicate gel (C-S-H), and specific chemical activities (Bahari et al. [3], Bartos et al. [4], Bittnar et al. [5], Sikora et al. [6]).
On the experimental level, several studies have been carried out over the last decades to test the mechanical performance of structures enhanced by incorporating nanoparticles in concrete matrix, namely silicon dioxide SiO 2 (Senff et al.) [7], titanium dioxide TiO 2 (Feng et al. [8]), zinc dioxide ZnO 2 (Nazari et al.) [9], and iron oxide Fe 2 O 3 (Nazari et al.) nanoparticles [10]. For example, Rong et al. [11] demonstrated that adding nano-SiO 2 to concrete matrix resulted in higher compressive and better tensile and flexural properties; due to their size in the range of 1-500 nm, they produced shorter setting time and water permeability, as well as strong resistance to chemical attacks. Because of this, they can form part of many concrete structures we see nowadays. Moreover, other considerable experimental studies have shown that the addition of silica nanoparticles SiO 2 improves not only the mechanical characteristics but also the physio-chemical characteristics of the fresh and hardened concrete [12]. Using an experimental approach, Mondal and his colleagues [13] confirmed that silica nanoparticle reinforcement is not only better for the environment but also yields better results. It can increase the durability and performance of concrete as well as save cement and lower the project's overall cost and environmental impact. However, the exorbitant price of nanoparticles may limit their utilization. Nazari et al. [10,14] revealed the chemical effect of incorporating oxidized nanoparticles in indoors fresh cement by increasing the degree of hydration, which leads to a greater consistency and workability of the cement paste. Ghannam et al. [15] conducted an experimental study to explore the possibility of using granite powder and iron powder to partially replace sand in concrete. The percentages of granite powder and iron powder added to replace the sand were 5%, 10%, 15%, and 20% of the sand by weight. It has been observed through experiments that the substitution of up to 20% of sand by weight with iron powder in concrete has resulted in an increase in the compressive and flexural strength of concrete [16]. Lee et al. [17] studied the effect of iron oxide additives on concrete interlocking blocs (CIBs). His study concluded a definite relationship between the flexural strength and the absorption rate of pigment-stained blocks. Kaïkea et al. [18] conducted an experimental study on six different sets of high performance fiber-reinforced concrete specimens by adding different types of mineral additives and checked their influence on the overall behavior of the material. The obtained results showed that by adding high volume content of fibers and iron oxide powder in the mix, the HPFRC material achieved very good mechanical performance.
Theoretically speaking, only a minority of research has been exposed using analytical modeling for the reinforcement of concrete structures with nanoparticles. Harrat et al. [19] conducted an analytical study of concrete beams incorporated with nano-silica particles using the Voigt model to describe the agglomeration effect of silica nanoparticles in the matrix and then studied the static behavior of nano-SiO 2 reinforced concrete beams resting on the Winkler-Pasternak foundation. When Chatbi et al. [20] also investigated the bending behavior of SiO 2 -reinforced concrete slabs resting on Kerr's elastic foundation, the most findings of those studies indicate that the use of an optimum amount of SiO 2 nanoparticles on concretes increases better mechanical behavior. In addition, they concluded that the elastic foundation has a significant impact on the bending of concrete structures. Amnieh et al. [21] compared experimentally and theoretically the dynamic analysis of SiO 2 nanoparticle mixed concrete blocks subjected to blast loading using Mindlin's plate deformation theory (FSDT) and proved that the experimental model and the theoretical one were very close. Shokravi [22] analytically investigated the nonlinear vibration of nano-silica reinforced concrete columns based on Timoshenko's beam (CBT) using the Mori-Tanaka homogenization model and then used the Differential Quadrature Method (DQM) to predict the frequency of the system. Many authors have been looking for ways to promote higher-order shear deformation theories (HSDT) to satisfy equilibrium conditions at the top and bottom of the plate without using shear correction factors (Reddy [23], Touratier [24], Soldatos [25], Mantari et al. [26]) and stay to this day very efficient. However, more importantly, on the accuracy of predicting the plate response, a recent series of articles developed a new refined and robust plate theory for the mechanical behavior of a simply supported plate with only four unknown variables. Bachir Bouiadjra et al. [27], Bourada et al. [28], and Kettaf et al. [29] used the refined theory for studying the thermal buckling of functionally graded sandwich plates. Zenkour et al. [30] studied the thermoelastic bending analysis of functionally graded ceramic-metal sandwich plates, using several kinds of sandwich plates and taking into account the symmetry of the plate and the thickness of layers. Tounsi et al. [31] used a refined trigonometric shear deformation theory to analyze the thermoelastic bending of functionally graded sandwich plates. Sayyad et al. [32] have also investigated the thermoelastic bending but on laminated composite plates according to various shear deformation theories. Radwan [33] presented a hyperbolic displacement model to investigate the nonlinear hygrothermal and mechanical buckling responses of FG sandwich plates resting on two-parameter elastic foundations while considering a non-linear Fourier temperature distribution as a special case. Hellal et al. [34] proposed a new simple "four-variable shear deformation" plate model to demonstrate the hygrothermal environment effects on dynamic and buckling of functionally graded material sandwich plates supported by Winkler-Pasternak elastic foundations.
In order to increase the strength and improve the mechanical properties of concrete, first the industrial iron wastes that pollute the environment were recovered and evaluated; then their chemical composition and their use in different forms in the concrete mix were analyzed.
Regarding the reinforcement of concrete by iron oxide nanoparticles, analytical studies are not found in the literature contrary to the experimental researches [35,36], and, therefore, there is a great lack in this subject despite the attractiveness of the concept of using Fe 2 O 3 nanoparticles in enhancing concrete. For this purpose, the authors aim through this paper to present a mathematical model simulating an equivalent system. In this study, thermoelastic bending analysis is conducted to analyze the behavior of concrete slabs reinforced with ferric oxide and subjected to different varying thermal loads. A recently developed refined trigonometric shear plate theory (RTSDT) is adopted in this work to simulate plate deformations. The theory contains only four unknown functions, against five in the case of other shear deformation theories. The transverse shear stresses vary parabolically through the thickness, which satisfies the free surface conditions without shear stress. The equilibrium equations are obtained based on the principle of virtual work. Using Navier's analytical solutions for thermoelastic bending, the effects of different parameters such as volume percent of nanoparticles, thermal loading type, and geometry parameters on the stress and displacement analysis of the plate are determined and presented.

Homogenization Model
Homogenization approaches are widely used today to determine the effective properties of heterogeneous mediums. They introduce two distinct scales: the microscopic scale, where the heterogeneity of the materials is obvious, and the macroscopic scale (engineering scale), where the materials are considered homogeneous ( Figure 1). Moreover, the scales must be strictly separated. This results in homogenized properties that do not particularly depend on the size of the microstructural elements. The effective properties depend only on the structural morphology and phase properties. In this paper, the analytical approach aims to determine the thermoelastic properties of reinforced concrete on the Representative Volume Element (R.V.E) scale using Eshelby's homogenization model [37], which is limited to the average properties of composite materials with simple microstructures.  [37], which is limited to the average properties of composite materials with simple microstructures. Eshelby's homogenization model is basically the most important analytical model for predicting the properties of nano-composite reinforced matrix. Strictly speaking, it only applies to the placement of an ellipsoidal inclusion in an infinite matrix; therefore, in this study, the iron nanoparticles are considered of a spherical form. The stiffness tensor CT for the nano-composite is given by Equation (1): in which I is the identity matrix, and Cm and Cr are the stiffness tensors for the concrete matrix and the nano-iron reinforcement, respectively. Meanwhile, Vm and Vr are the volume fraction of the matrix and reinforcement, and S is the Eshelby tensor which is related to the Poisson ratios of nanoparticles. For both isotropic materials, the stiffnesses Cm and Cr are expressed in Equation (2): In which E is the Young's modulus of either the concrete matrix or the iron nanoparticle reinforcement, and υ is Poisson's ratios. The indexes 1, 2, and 3 conform to the x, y, z directions of the composite Cartesian co-ordinate system, respectively.
Reinforcement with spherical form S is given in Equation (3) (Clyne and Withers [38]): where the Eshelby's matrix components are given in Equation (4): Eshelby's homogenization model is basically the most important analytical model for predicting the properties of nano-composite reinforced matrix. Strictly speaking, it only applies to the placement of an ellipsoidal inclusion in an infinite matrix; therefore, in this study, the iron nanoparticles are considered of a spherical form. The stiffness tensor CT for the nano-composite is given by Equation (1): in which I is the identity matrix, and C m and C r are the stiffness tensors for the concrete matrix and the nano-iron reinforcement, respectively. Meanwhile, V m and V r are the volume fraction of the matrix and reinforcement, and S is the Eshelby tensor which is related to the Poisson ratios of nanoparticles. For both isotropic materials, the stiffnesses C m and C r are expressed in Equation (2): In which E is the Young's modulus of either the concrete matrix or the iron nanoparticle reinforcement, and υ is Poisson's ratios. The indexes 1, 2, and 3 conform to the x, y, z directions of the composite Cartesian co-ordinate system, respectively.
Reinforcement with spherical form S is given in Equation ( where the Eshelby's matrix components are given in Equation (4): In relation (4), υ r denotes the Poisson's ratio of nanoparticle reinforcements. For particulate composites, Eshelby's approach for the prediction of the thermal expansion α T is presented in Equation (5): In Equation (5), α r and α m are the reinforcement and matrix expansion tensors, respectively. C m and C r are the same stiffness tensors, but S tensor is now appropriate for heat properties, in which For all others S ij = 0 (6)

Mathematical Modeling of the Plate
A simply supported flat concrete slab having a length 'a', width 'b', and total thickness 'h' reinforced with iron oxide particles is considered in this study. The Fe 2 O 3 nanoparticles are supposed to be randomly distributed in the concrete matrix as illustrated in Figure 2. The chosen coordinate system (x, y, z) is also shown in Figure 2, in which Materials 2023, 16 In relation (4), υr denotes the Poisson's ratio of nanoparticle reinforcements. For particulate composites, Eshelby's approach for the prediction of the thermal expansion αT is presented in Equation (5): In Equation (5), αr and αm are the reinforcement and matrix expansion tensors, respectively. Cm and Cr are the same stiffness tensors, but S tensor is now appropriate for heat properties, in which

Mathematical Modeling of the Plate
A simply supported flat concrete slab having a length 'a', width 'b', and total thickness 'h' reinforced with iron oxide particles is considered in this study. The Fe2O3 nanoparticles are supposed to be randomly distributed in the concrete matrix as illustrated in Figure 2. The chosen coordinate system (x, y, z) is also shown in Figure 2, in which 0 ≤ ≤ ; 0 ≤ ≤ ; −ℎ/2 ≤ ≤ ℎ/2. For comparison purposes, different higher order shear deformation theories (HSDT's) were used in our analysis. The HSDT displacement field of a material point located at (x, y, z) in the plate can be described in Equation where U, V, and W are displacements in x, y, z directions; u0, v0, and w0 are mid-plane translations, and θx and θy are rotations of normal to the mid-plane around y-axis and xaxis, respectively. f(z) represents the shape function that determines the distribution of For comparison purposes, different higher order shear deformation theories (HSDT's) were used in our analysis. The HSDT displacement field of a material point located at (x, y, z) in the plate can be described in Equation (8): where U, V, and W are displacements in x, y, z directions; u 0 , v 0 , and w 0 are mid-plane translations, and θx and θy are rotations of normal to the mid-plane around y-axis and x-axis, respectively. f (z) represents the shape function that determines the distribution of transverse shear strain and stress through the thickness. The displacement of parabolic shear deformation plate theory (PSDPT) is obtained by setting (Reddy [23]) Additionally, the sinusoidal shear deformation plate theory (SSDPT) is obtained by setting (Touratier et al. [24]) In addition, the displacement of exponential shear deformation plate theory (ESDPT) is obtained by setting (Karama et al. [39])

Assumptions of the Trigonometric Refined Plate Theory (RTSDT)
Unlike other theories, the number of unknown functions included in the refined triangular shear strain theory (RTSDT) is only four, compared to five in other shear strain theories (Reddy [23], Touratier et al. [24], Karama et al. [39]). RTSDT does not require a shear correction factor, such as in the case of the first order deformation theory (Mindlin 1951), and introduces a parabolic variation of transverse shear stress through the thickness of the plate.
The main assumptions of the refined trigonometric shear deformation theory are recalled in the following:

•
The displacements are small besides the plate total thickness 'h'; therefore, strains involved are considered infinitesimal.

•
The transverse displacement 'w' is divided into two components, including 'w b ' and 'w s ', accounting for the bending and shear effects; these components are the function of the spatial coordinates x, y and of time t only.
As a result of the previous assumptions, the displacements 'U' in 'x' direction and 'V' in 'y' direction can be expressed by means of the translation, bending, and shear components: The bending components 'u b ' and 'v b ' are assumed to be identical to the displacements given by the classical plate theory. Therefore, the expression for 'u b ' and 'v b ' can be given by Equation (14): The shear components 'u s ' and 'v s ' are defined in this theory to express the shear displacement effect 'w s ' to different variations of shear strains γ xz and γ yz and, hence, to shear stresses σ xz and σ yz across the thickness of the plate in such a way that shear stresses σ xz and σ yz are zero at the top and bottom faces of the plate. Consequently, 'u s ' and 'v s ' can be put as in Equation (15): Here, 'f (z)' represents the shape function determining the transverse shear distribution through the plate thickness. The transverse normal stress 'σ z ' is assumed to be negligible in comparison with the in-plane stresses 'σ x ' and 'σ y '. The body force of the plate is also neglected in this analysis.

Kinematics
Based on Equations (12)-(15), the displacement field in the trigonometric refined theory can be written as: where 'u 0 ' and 'v 0 are the mid-plane displacements of the plate along the x and y direction. 'w b ' and 'w s ' are the bending and shear components of transverse displacement in z direction, respectively. In this analysis, the trigonometric shape function is presented as It should be noted that, unlike the first-order shear deformation theory, refined plate theories do not require shear correction factors.
The linear strains components can be derived from the displacement field in Equation (16) and are given in Equation (18): The constitutive stress-strain relations of the nano-composite can be defined as where C ij T are the reduced elastic constants of the concrete slab reinforced with Fe 2 O 3 , obtained using Eshelby's homogenization model.

Principle of Virtual Work
The virtual work principle is applied to develop the equations of motion in the present case yields: where δΨ and δΦ are the virtual variation of the internal strain energy and the virtual work performed by external forces, respectively. The expression of the virtual strain energy is By substituting Equation (18) into Equation (22), one finds where stress resultants can be expressed in Equation (24): For the concrete plates subjected to bending loads "q", the virtual work performed by external loadings is The equilibrium equations associated with the trigonometric refined shear deformation theory can be derived from Equation (21) by integrating parts and collecting the coefficients of δu 0 , δv 0 , δw b , and δw s and are given in Equation (26): By substituting Equation (20) into Equation (24), one obtains the stress resultants in form of material stiffness and displacement components: where A ij , B ij , D ij , Bs ij , Ds ij , Hs ij , and As ij are the plate stiffness, defined in Equation (28): The normal stress and moment resultants N x T = N It is assumed that the variation of the temperature field through the thickness is expressed as Here, T 1 , T 2 and T 3 are thermal loadings, and for the case of the present refined trigonometric theory (RTSDT), Substituting from Equation (27) into Equation (26), the following equations are obtained: The generalized force vector is described as and d i jlm are differential operators expressed as the following: The components of the generalized force vector {f } are

Navier's Technique
To formulate the closed-form solutions for bending problems of simply supported rectangular plates for the refined trigonometric shear deformation theory, Navier's method is employed: where U mn , V mn , W bmn , and W smn are the arbitrary parameters to be determined, α = mπ/a, and ϕ = nπ/b. The transverse load 'q' for the mechanical bending is expressed by the double-Fourier sine series by The various load patterns are illustrated in Figure 3, where the coefficients Q mn are presented below:   Navier presented the transverse temperature loads T 1 , T 2 , and T 3 in the form of a double trigonometric series as Substituting Equation (35) into Equation (26), one obtains the closed-form solutions in matrix form: where The components of the elastic stiffness matrix [K ij ] are as follows: where the elements of the stiffness matrix are given in Equation (42): The components of the generalized thermal force vector {f } = {f 1 , f 2 , f 3 , f 4 } T are given by

Results and Discussion
To analyze the static behavior of simply supported nano-iron-reinforced concrete slabs, various numerical examples are presented and discussed to predict the mechanical and thermoelastic bending of the slab by using the refined trigonometric shear deformation theory (RTSDT).
Herein is presented a concrete slab with Young's modulus of E m = 20 GPa and thermal expansion α m = 13.5 (10 −6 /K). The concrete matrix is incorporated with iron nanoparticles for Young's modulus of E r = 200 GPa and thermal expansion α r = 12 (10 −6 /K). Poisson's ratios are υ m = 0.3 and υ r = 0.291 for concrete and iron nanoparticles, respectively. Dimensionless displacements and stresses generated by the mechanical external loads are presented according to the following definitions: Dimensionless displacements and stresses generated by the temperature loadings are expressed by the following definitions: To demonstrate the preceding thermal-structural analysis, different sample problems are considered. For the sake of brevity, only linearly varying (across the thickness) temperature distribution T = zT 2 , non-linearly varying (across the thickness) temperature distribution T = Ψ(z)T 3 , and a combination of both T = zT 2 + Ψ(z)T 3 are treated, where the reference values are E 0 = 1 GPa and α 0 = 10 −6 /K. Note that in most of the literature, the thermal stress problems are considered under a steady state temperature distribution that is linear with respect to the thickness direction.

Validation
First, it is decisive to verify the adequacy of the present adopted (RTSDT) mathematical model since analytical results for the mechanical bending analysis of concrete slabs reinforced with iron nanoparticles are not available in the literature. Considering the material and geometric parameters similar to Thai and Choi [40], the results in terms of transverse displacements (ω) as well as normal stresses (σ x , σ xy ) of functionally graded (FG) slabs while varying the power index 'P' are embraced for comparison with the present refined trigonometric shear deformation theory as presented in Table 1. As shown in Table 1, the results are generally in fair agreement. Using different shear deformation theories or shape functions gives similar results of displacements and normal stresses. However, the refined trigonometric shear deformation theory (RTSDT) predicts almost the same results as Zenkour's refined sinusoidal shear deformation plate theory (RSSDPT).
In order to validate the elastic properties estimated by the Eshelby analytical model of a concrete matrix reinforced with ferric oxide (Fe 2 O 3 ) nanoparticles, and in the absence of studies in the literature dealing with a similar analytical model, a comparison was made between the elastic stiffnesses (C ij ) obtained from a concrete slab reinforced with iron nanoparticles (using the Eshelby homogenization model) and the elastic stiffnesses (C ij ) obtained from a concrete slab reinforced with silica nanoparticles using the Voigt homogenization model without taking into account the agglomeration effect of (SiO 2 ) nanoparticles in the concrete matrix. Figure 4 shows that both types of reinforcement (Fe 2 O 3 and SiO 2 ) have the same effect on the concrete, in which the reduced elastic stiffnesses (C ij ) of a concrete slab increase with the concentration of reinforcements (V r ). Yet, it is especially noticeable in Figure 4 that, for iron nanoparticles (Fe 2 O 3 ) reinforcements, the stiffnesses are further more improved compared to nano-silica (SiO 2 ) reinforcements, particularly in the case of (C 12 ) which denotes the elastic stiffness of the plate in (x, y) plan. This improvement is mainly due to the nanoparticle's high mechanical properties such as Young's modulus.

Part I: Bending Analysis
To begin with, a comparison between the non-dimensional transverse displacement values obtained using the refined trigonometric theory (RTSDT) is shown in Figure 5 as a function of reinforcement volumes. One of the curves expresses the development of the non-dimensional deflection (ω) of a concrete slab reinforced with iron nanoparticles Fe 2 O 3 . These material properties are determined by Eshelby's homogenization model, and the other one expresses the variation of (ω) as well but this time of a concrete slab reinforced by silica nanoparticles, on which these elastic properties are determined by the Voigt homogenization law without considering the agglomeration effect of the SiO 2 nanoparticles (Harrat et al., 2021 [19]).

Part I: Bending Analysis
To begin with, a comparison between the non-dimensional transverse displacemen values obtained using the refined trigonometric theory (RTSDT) is shown in Figure 5 as function of reinforcement volumes. One of the curves expresses the development of th non-dimensional deflection (ω) of a concrete slab reinforced with iron nanoparticles Fe2O These material properties are determined by Eshelby's homogenization model, and th other one expresses the variation of (ω) as well but this time of a concrete slab reinforce by silica nanoparticles, on which these elastic properties are determined by the Voigt h mogenization law without considering the agglomeration effect of the SiO2 nanoparticle (Harrat et al., 2021 [19]).
It is clear from Figure 5 that the reinforcement by iron nanoparticles is more effectiv than by nano-silica, and this is probably due to the high physical properties of the nan iron (Young's elastic modulus, precisely), which leads to a higher stiffness of the concre nano-composite.
Thus, we should retain from Figure 5 that when the plate is reinforced by nano-iro volumes, it gives us a lower deflection than a plate reinforced by SiO2. It is also wort noting that the higher the percentage of particles is, the lower the deflection is.   Figure 6 shows the variation of the non-dimensional transverse displacement of concrete plate reinforced with different volumes (Vr) of iron nanoparticles. Several theo ries are used to determine the accuracy of the theory used in our analysis (the refine trigonometric shear deformation theory) in predicting the mechanical behaviors of con crete nano-composite.
From Figure 6, it should be noted that all the theories have the same behavior, bu the (RTSDT) theory slightly overestimates the deflection compared to the others, and th is due to the used shape function that expresses the transverse shear stresses that evolv through the thickness of the plate. It is clear from Figure 5 that the reinforcement by iron nanoparticles is more effective than by nano-silica, and this is probably due to the high physical properties of the nanoiron (Young's elastic modulus, precisely), which leads to a higher stiffness of the concrete nano-composite.
Thus, we should retain from Figure 5 that when the plate is reinforced by nano-iron volumes, it gives us a lower deflection than a plate reinforced by SiO 2 . It is also worth noting that the higher the percentage of particles is, the lower the deflection is. Figure 6 shows the variation of the non-dimensional transverse displacement of a concrete plate reinforced with different volumes (V r ) of iron nanoparticles. Several theories are used to determine the accuracy of the theory used in our analysis (the refined trigonometric shear deformation theory) in predicting the mechanical behaviors of concrete nano-composite. Figure 6 shows the variation of the non-dimensional transverse displacement of a concrete plate reinforced with different volumes (Vr) of iron nanoparticles. Several theories are used to determine the accuracy of the theory used in our analysis (the refined trigonometric shear deformation theory) in predicting the mechanical behaviors of concrete nano-composite.
From Figure 6, it should be noted that all the theories have the same behavior, but the (RTSDT) theory slightly overestimates the deflection compared to the others, and this is due to the used shape function that expresses the transverse shear stresses that evolve through the thickness of the plate.  From Figure 6, it should be noted that all the theories have the same behavior, but the (RTSDT) theory slightly overestimates the deflection compared to the others, and this is due to the used shape function that expresses the transverse shear stresses that evolve through the thickness of the plate.
Different subjected loading patterns on a rectangular (a = b), concrete slab reinforced with different concentrations of iron nanoparticles (Fe 2 O 3 ) are shown in Figure 7.  Figure 7.
The slab is considered simply supported ( Figure 2). As it was noticed in the previous figures, the reinforcement of concrete slabs with iron nanoparticles has a strengthening role on the structure, but from Figure 7, we can say that this is actual regardless of the type of loading. It should also be noted that uniform loading has a greater effect on the bending of reinforced plates than other types of loading, while concentrated loading has the least effect.   The slab is considered simply supported ( Figure 2). As it was noticed in the previous figures, the reinforcement of concrete slabs with iron nanoparticles has a strengthening role on the structure, but from Figure 7, we can say that this is actual regardless of the type of loading. It should also be noted that uniform loading has a greater effect on the bending of reinforced plates than other types of loading, while concentrated loading has the least effect. Figure 8 is used to discuss the non-dimensional transverse displacement (ω) of a simply supported square concrete plate reinforced with iron nanoparticles over the entire length (x/a) of the plate by varying the reinforcement volume from 0% (non-reinforced concrete plate) to V r = 30%. The plate is supposed to be subjected to sinusoidal loading.  Figure 9 shows the variation of the non-dimensional axial displacement (U) of a non reinforced concrete plate (Figure 9a) and of a concrete plate reinforced by iron nanoparti cles Vr = 30% (Figure 9b). The plates are assumed square and simply supported on thei line edges. Different loading patterns (single sine, uniformly, concentrated, and linearl distributed loads) are applied on the plates to determine the loading effect on the axia displacement of concrete plates. It is to be noted that the uniform load has a much highe effect on the axial displacement of concrete plates reinforced by iron nanoparticles com pared to the other types of loadings, which is the same as the case of transverse displace ments mentioned before. Markedly, iron nanoparticles reinforcement in a concrete matri reduces the non-dimensional axial displacement (U) of the composite. Table 2 shows the non-dimensional normal ( , ) and shear stresses ( ) of simply supported rectangular (a = b) concrete slab under sinusoidal loading. The concret slab is reinforced with different volumes of ferric oxide nanoparticles (Vr = 0%, 10%, 20% 30%), while varying the geometric ratio (a/h). The constraints became nearly unchangeabl when a/h exceeded 30.   Figure 8, shows that the higher deflections are obtained in the mid-length of the plate for all nano-Fe 2 O 3 concentrations. It also confirms the stiffening effect caused by iron nanoparticles by decreasing the deflection of the nano-composite plate, knowing that the deflection of a reinforced plate (V r = 30%) is way lower (almost by half) than the non-reinforced concrete plate (V r = 0%). Figure 9 shows the variation of the non-dimensional axial displacement (U) of a nonreinforced concrete plate (Figure 9a) and of a concrete plate reinforced by iron nanoparticles Vr = 30% (Figure 9b). The plates are assumed square and simply supported on their line edges. Different loading patterns (single sine, uniformly, concentrated, and linearly distributed loads) are applied on the plates to determine the loading effect on the axial displacement of concrete plates. It is to be noted that the uniform load has a much higher effect on the axial displacement of concrete plates reinforced by iron nanoparticles compared to the other types of loadings, which is the same as the case of transverse displacements mentioned before. Markedly, iron nanoparticles reinforcement in a concrete matrix reduces the non-dimensional axial displacement (U) of the composite.
reduces the non-dimensional axial displacement (U) of the composite. Table 2 shows the non-dimensional normal ( , ) and shear stresses ( ) of a simply supported rectangular (a = b) concrete slab under sinusoidal loading. The concrete slab is reinforced with different volumes of ferric oxide nanoparticles (Vr = 0%, 10%, 20%, 30%), while varying the geometric ratio (a/h). The constraints became nearly unchangeable when a/h exceeded 30.  Table 2 shows the non-dimensional normal (σ x , σ xy ) and shear stresses (τ xz ) of a simply supported rectangular (a = b) concrete slab under sinusoidal loading. The concrete slab is reinforced with different volumes of ferric oxide nanoparticles (V r = 0%, 10%, 20%, 30%), while varying the geometric ratio (a/h). The constraints became nearly unchangeable when a/h exceeded 30. Table 2. Variation of the dimensionless shear stresses σ x , σ xy , and τ xz of reinforced concrete square plates using the refined trigonometric shear deformation theory (RTSDT). In Table 2, it can be seen also that the reinforcement concentrations of iron nanoparticles have a very slight effect on the shear stresses (τ xz ) which evolve through the thickness of the composite slab. It can be taken that the shear stresses (τ xz ) are almost stable when varying reinforcement volumes (V r ). In reverse, the reinforcement concentration (V r ) interestingly decreases the normal stress (σ x ) and increases (σ xy ) of the nano-reinforced plate.

Part II: Thermoelastic Bending
As with mechanical bending, mathematical modeling using the refined trigonometric shear deformation theory (RTSDT) must be verified analytically for predicting the thermoelastic bending of slabs. Since there are no models of concrete slabs reinforced with iron nanoparticles in the literature, it is quite appropriate to check the accuracy of the present mathematical model (RTSDT) using the functionally graded (FG) sandwich-structured plate of Zenkour (2008) [30].
By setting 'P = 0' and considering the material and geometrical parameters similar to those of Zenkour (2008) [30], the results for the transverse displacements (ω T ) as well as the normal and shear stresses (σ x , τ xZ ) of the sandwich plates subjected to a linearly varying thermal load (T 2 = 100) through the thickness of the plate are listed for comparison in Table 3. Table 3. Validation of present theory against other published theories in predicting the thermomechanical bending of FG plate, (a/h = 10, a = b, T 2 = 100 , T 3 = T 1 = 0). As observed, the different plate theories are very consistent in predicting thermoelastic bending. Similar displacement, normal, and shear stresses results can be obtained using different shear deformation theories or shape functions. Figure 10 shows the non-dimensional transverse displacement (ω T ) predicted by different plate theories as a function of the concentration of iron nanoparticles (Vr) in concrete matrix. The concrete slab is assumed to be simply supported and subjected to a thermal load that varies linearly through the plate thickness (T 2 = 100, T 3 = T 1 = 0). The parabolic shear deformation theory (PSDPT) seems to have a slight overestimation over the other plate theories, while the exponential (ESDPT) and sinusoidal (SSDPT) shear theories are corresponding. Otherwise, the plate theories are in accordance when predicting the transverse deflection of the nano-composite concrete slab under thermal load. In contrast to the case of the mechanical bending, here the reinforcement by nanoiron (Fe2O3) has a weakening effect on the slab behavior when considering the thermal load. It can be clearly seen that by increasing the volume of nano-iron (Vr), the deflection increases, and this is due to the thermoelastic properties of the nano-iron reinforcements (the high thermal expansion of the iron nanoparticles). The effect of geometric ratios (a/h and a/b) on the non-dimensional transverse displacement ( ) of a simply supported concrete plate reinforced with various proportions of iron oxide (Fe2O3) nanoparticles is illustrated in Figure 11.

Shape Function
In Figure 11a, the effect of the length/thickness (a/h) ratio of a plate subjected to a linearly varying thermal load is shown with regard to the different nanometric ferric oxide incorporation volumes (Vr), whereas Figure 11b shows us the effect of the length/width (a/b) ratio of a concrete plate subjected to a combination of linear and non-linear thermal field ( = = 100). In contrast to the case of the mechanical bending, here the reinforcement by nano-iron (Fe 2 O 3 ) has a weakening effect on the slab behavior when considering the thermal load. It can be clearly seen that by increasing the volume of nano-iron (V r ), the deflection increases, and this is due to the thermoelastic properties of the nano-iron reinforcements (the high thermal expansion of the iron nanoparticles).
The effect of geometric ratios (a/h and a/b) on the non-dimensional transverse displacement (ω T ) of a simply supported concrete plate reinforced with various proportions of iron oxide (Fe 2 O 3 ) nanoparticles is illustrated in Figure 11. A comparison between the thermoelastic deflection ( ) along the length (x/a) of an unreinforced concrete slab and a concrete slab incorporated with iron nanoparticles of volume Vr = 30% is shown in Figure 12. Both plates are considered simply supported and subjected to linearly varying thermal load ( = 100). It is too clear by seeing Figure 12 that when Vr is 30%, the plate deflects even more in comparison to the non-reinforced concrete plate.  Figure 13 illustrates the variation of the non-dimensional transverse displacement of a reinforced concrete slab (Vr = 30%) subjected to different thermal loads, as a function of different geometric ratios.
In Figure 13, it can be seen that the length-to-thickness (a/h) ratio has an increasing effect on the deflection of the plate. However, when the ratio (a/h) exceeds 10, the deflection starts to be stable. On the other hand, the length-to-width ratio (a/b) has a sliding In Figure 11a, the effect of the length/thickness (a/h) ratio of a plate subjected to a linearly varying thermal load is shown with regard to the different nanometric ferric oxide incorporation volumes (V r ), whereas Figure 11b shows us the effect of the length/width (a/b) ratio of a concrete plate subjected to a combination of linear and non-linear thermal field (T 2 = T 3 = 100).
As mentioned before, these results also confirm nano-iron's weakening effect on the plate, and this is mainly due to the thermoelastic properties (high thermal expansion) of the nano-iron incorporated in the concrete matrix. Therefore, the higher the percentage of particles present in the concrete matrix, the higher the deflection became, regardless of the geometric ratios.
A comparison between the thermoelastic deflection (ω T ) along the length (x/a) of an unreinforced concrete slab and a concrete slab incorporated with iron nanoparticles of volume V r = 30% is shown in Figure 12. Both plates are considered simply supported and subjected to linearly varying thermal load (T 2 = 100). It is too clear by seeing Figure 12 that when Vr is 30%, the plate deflects even more in comparison to the non-reinforced concrete plate.
In Figure 13, it can be seen that the length-to-thickness (a/h) ratio has an increasing effect on the deflection of the plate. However, when the ratio (a/h) exceeds 10, the deflection starts to be stable. On the other hand, the length-to-width ratio (a/b) has a sliding effect on the thermoelastic deflection of the plate. It can be noticed also from Figure 13 that the application of combined linear and non-linear thermal loads (T 3 = 100) has a greater effect on the behavior of the plate in comparison to the other loads. In general, these results reveal that the variation of transverse displacement is very sensitive to the variation of the thermal load T 3 value. Figure 13 illustrates the variation of the non-dimensional transverse displacement of a reinforced concrete slab (V r = 30%) subjected to different thermal loads, as a function of different geometric ratios.
A comparison between the thermoelastic deflection ( ) along the length (x/a) of an unreinforced concrete slab and a concrete slab incorporated with iron nanoparticles of volume Vr = 30% is shown in Figure 12. Both plates are considered simply supported and subjected to linearly varying thermal load ( = 100). It is too clear by seeing Figure 12 that when Vr is 30%, the plate deflects even more in comparison to the non-reinforced concrete plate.  Figure 13 illustrates the variation of the non-dimensional transverse displacement of a reinforced concrete slab (Vr = 30%) subjected to different thermal loads, as a function of different geometric ratios.
In Figure 13, it can be seen that the length-to-thickness (a/h) ratio has an increasing effect on the deflection of the plate. However, when the ratio (a/h) exceeds 10, the deflection starts to be stable. On the other hand, the length-to-width ratio (a/b) has a sliding effect on the thermoelastic deflection of the plate. It can be noticed also from Figure 13 that the application of combined linear and non-linear thermal loads ( = 100) has a greater effect on the behavior of the plate in comparison to the other loads. In general, these results reveal that the variation of transverse displacement is very sensitive to the variation of the thermal load value. Figure 13. Effect of geometric ratios on the transverse displacement of a concrete slab under different thermal loads. Figure 14 contains the axial stress ( ) plot through the thickness (z/h) of an unreinforced concrete slab (Vr = 0%) and a ferric oxide Fe2O3 reinforced concrete slab (Vr = 30%). It is evident that the stresses are tensile below the mid-plane (z/h = 0) and compressive above the mid-plane for both plates (Figure 14a,b). The axial stress is continuous through the thickness of the slab. The results also show a nonlinear variation of the axial stress through the thickness of the plate for both cases. The maximum axial stresses are those generated by the thermal load ( = 100).   Figure 14 contains the axial stress (σ x ) plot through the thickness (z/h) of an unreinforced concrete slab (V r = 0%) and a ferric oxide Fe 2 O 3 reinforced concrete slab (V r = 30%). It is evident that the stresses are tensile below the mid-plane (z/h = 0) and compressive above the mid-plane for both plates (Figure 14a,b). The axial stress is continuous through the thickness of the slab. The results also show a nonlinear variation of the axial stress through the thickness of the plate for both cases. The maximum axial stresses are those generated by the thermal load (T 3 = 100). Figure 15 shows the through-thickness distributions of transverse shear stresses τ xz in the unreinforced and reinforced concrete plates. It is observed that the maximum value occurs at the center of the plate as in the case of the homogeneous plate. It is also noticeable that the maximum stresses are marked when the reinforcement volume is 30%.
forced concrete slab (Vr = 0%) and a ferric oxide Fe2O3 reinforced concrete slab (Vr = 30%). It is evident that the stresses are tensile below the mid-plane (z/h = 0) and compressive above the mid-plane for both plates (Figure 14a,b). The axial stress is continuous through the thickness of the slab. The results also show a nonlinear variation of the axial stress through the thickness of the plate for both cases. The maximum axial stresses are those generated by the thermal load ( = 100).  Figure 15 shows the through-thickness distributions of transverse shear stresses τxz in the unreinforced and reinforced concrete plates. It is observed that the maximum value occurs at the center of the plate as in the case of the homogeneous plate. It is also noticeable that the maximum stresses are marked when the reinforcement volume is 30%.

Conclusions
To design a nano-concrete, it is imperative to perform an analytical model to make preliminary study on the effect of different added nanostructures on the physical charac teristics of concrete.
In this paper, the static behavior of concrete slabs reinforced with iron nanoparticle was analyzed. A refined trigonometric plate theory (RTSDT) was used to simulate th plate, which is assumed simply supported and subjected to different thermal and mechan ical loads. The Eshelby homogenization model was used to determine the thermoelasti properties of the nano-composite and virtual work principle to determine the governin equations.
Based on this study, the findings are the following: • The use of iron nanoparticles as reinforcement in a concrete matrix can improve it

Conclusions
To design a nano-concrete, it is imperative to perform an analytical model to make a preliminary study on the effect of different added nanostructures on the physical characteristics of concrete.
In this paper, the static behavior of concrete slabs reinforced with iron nanoparticles was analyzed. A refined trigonometric plate theory (RTSDT) was used to simulate the plate, which is assumed simply supported and subjected to different thermal and mechanical loads. The Eshelby homogenization model was used to determine the thermoelastic properties of the nano-composite and virtual work principle to determine the governing equations.