Numerical Computation
In this study, results for radial stress, tangential stress and displacement are obtained numerically by using the FEM based solver COMSOL Multiphysics
. The governing differential equation as given by Equation (
3), under the prescribed boundary conditions given by Equation (
8), is solved by using a solid mechanics module and axisymmetric geometric conditions. Using the global analytical function, Young’s modulus is defined over the domain of an axisymmetric functionally graded cylinder. As depicted in
Figure 3, the domain of the cylinder is discretized with triangular elements into an extremely fine mesh by considering 5034 domain elements and 240 boundary elements. The model is then solved using a linear direct PARDISO solver and results for radial stress, tangential stress and displacement are obtained for a functionally graded cylinder with a relative tolerance of
. Computation of stresses and displacement is also carried out for non-FGM material using a parametric sweep. The numerical solution obtained is then compared with an iterative solution for both FGM and non-FGM materials. Radial displacement and stresses have been solved up to two iterations using the iterative technique. After two iterations, a good agreement can be observed between the results obtained from iterative technique and FEM, which can further be improved by considering a greater number of iterations.
In
Table 1, a comparison of the radial displacement in the FGM cylinder by FEM and iterative methods is shown, where the values in both methods are decreasing from the inner to the outer radius. Minimum and maximum absolute percentage errors of
and
occur at internal and external radii respectively. In
Table 2, a comparison of radial stress in the FGM cylinder obtained by FEM and iterative methods is shown for different radial points. This table shows that radial stress decreases along the radius of the cylinder. At the internal radius, the absolute error is
and at the outer radius, an absolute error of
is observed; overall, % error along the radius of the cylinder is less than 1. Maximum error can be seen at radius
which is
. Derived values of tangential stress for the FGM cylinder by FEM and iterative methods are shown for different radial points in
Table 3. Tangential stress is decreasing from the inner to the outer radius of the FGM cylinder. From
Table 3, we can observe that the propagation of error has a decreasing trend from inner to outer radial points. Maximum absolute error is obtained at the internal radius as
and the minimum error at the external radius is
.
Figure 4,
Figure 5,
Figure 6,
Figure 7,
Figure 8 and
Figure 9 present stresses and displacement in FGM and non-FGM cylinders under internal/external pressure conditions obtained using the iterative technique.
Figure 4 shows radial stress under internal pressure,
MPa. Under the effect of internal pressure, radial stress in the FGM cylinder is higher in magnitude as compared to non-FGM. It can be observed due to higher internal pressure and increasing elasticity modulus from inner to outer radii of the FGM cylinder. It can be noted from
Figure 4 that radial stress is on the higher side at the internal radial points of the cylinder due to low elasticity modulus at the inner radial points as compared to the outer radial points. Iterative radial stress in the cylinder under external pressure
MPa is depicted in
Figure 5. It can be observed that the magnitude of radial stress is lower in the FGM cylinder as compared to the non-FGM cylinder.
Figure 6 presents tangential stress in the cylinder under the influence of internal pressure
MPa. As seen from
Figure 6, tangential stress under internal pressure is tensile throughout the radius of the cylinder whereas from
Figure 7, tangential stress under the effect of external pressure is found to be compressive throughout the radius of the cylinder. Moreover, under the internal pressure, the magnitude of tangential stress is found to decrease from inner to outer radii of the cylinder.
Tangential stress under internal pressure is on the higher side in the non-FGM cylinder but at outer radial points it is found to be higher for the FGM cylinder. Under the effect of external pressure, tangential stress in the FGM cylinder increases in magnitude towards the outer radius of the cylinder whereas in the non-FGM cylinder, tangential stress becomes less compressive as it moves along the outer radius of the cylinder.
As seen from
Figure 8, radial displacement under the influence of internal pressure has tensile values for both FGM and non-FGM cylinders whereas radial displacement under the influence of external pressure exhibits the compressive behaviour, as shown in and
Figure 9. Additionally, the magnitude of radial displacement in the case of the FGM cylinder is lower than that of the non-FGM cylinder under both internal and external pressure conditions as shown in
Figure 8 and
Figure 9. This behaviour can be attributed to the higher values of Young’s modulus for the FGM cylinder and the lower values for the non-FGM cylinder.
In
Figure 10, it can be observed that under the effect of high internal pressure and low external pressure, the magnitude of radial stress decreases from inner to outer radial points, for both FGM and non-FGM cylinders. However, for cylinders made up of High Carbon Steel–Magnesium Oxide functionally graded material, the magnitude of radial stress is found to be on the higher side as compared to the non-FGM (High Carbon Steel) cylinder. This mechanical response of the cylinder under high internal pressure and low external pressure is due to its material property, namely, Young’s modulus, which increases from inner to outer radii, causing a high resisting force due to which the radial stress generated is of a higher magnitude in comparison to the non-FGM cylinder.
As seen from
Figure 11, under high external pressure and low internal pressure, radial stress for the FG cylinder is less as compared to the non-FGM cylinder, from internal to external radii of the cylinder. Hence, the amount of pressure required to move points from its original location is less for the FG cylinder. As the slope of radial stress decreases from inner to outer radii, this leads to producing a lower resistive force towards the outer radius of the functionally graded cylinder. This happens as Young’s modulus increases towards the outer radius; then the cylinder requires a lower resistive force to counter the external pressure. Hence, radial stress in the FG cylinder is lower in magnitude as compared to the non-FGM cylinder due to homogeneous Young’s modulus along the radial points.
As seen from
Figure 12, tangential stress from internal to external radii of the cylinder is found to be decreasing under high internal and low external pressure, but for the FGM cylinder, the magnitude of tangential stress is on the higher side at outer radial points as compared to the non-FGM cylinder. Tangential stresses under this case are found to be elastic throughout the radius of FGM and non-FGM cylinders.
From
Figure 13, under low internal and high external pressures, the magnitude of tangential stress for the FGM cylinder increases from inner to outer radii whereas in the non-FGM cylinder tangential stress is higher at inner radial points and decreases towards the outer radius of the cylinder.
It can be observed that the resistance force required by material points to move out from its external surface is less in the FGM cylinder as compared to the non-FGM cylinder and, hence, it saves a lot of stress.
From
Figure 14, it can be observed that under high internal and low external pressure conditions, displacement decreases throughout the cylinder and is found to be lowest for the FGM cylinder. The high internal pressure causes compressive strains with displacement values decreasing along the radius of the cylinder. Due to graded elasticity in the FGM cylinder, it has a lower magnitude of stress as compared to the non-FGM cylinder. Moreover, displacement in both FGM and non-FGM cylinders is found to be compressive towards the outer radial points.
Under low internal and high external pressure, as shown in
Figure 15, displacement in the non-FGM cylinder has a higher magnitude than in the FGM cylinder and is also found to be increasing along the radius of the cylinder whereas it decreases in the case of the FGM cylinder.
Figure 16 and
Figure 17 show a good agreement between iterative and finite element methods for the results of radial stress and tangential stress, respectively, in the case of the non-FGM cylinder.
From
Figure 18, it is observed that finite element method solution for radial displacement in non-FGM cylinder varies with iterative solution only at fifth position after decimal.
Table 4 represents radial displacement in tabular form for different iteration stages a
,
and
. In this table,
,
, and
show initial, first and second iterations for radial displacement, respectively. From this table, it can be observed that relative error
is approximately less than and equal to
as we iterate from
to
. The minimum relative error of
occurs at internal surface
of the cylinder and the maximum relative error of
occurs at external surface
of the cylinder.
Table 5 shows radial stress values in tabular form for the considered number of iterations in the study. It can be noted that relative error from the first iteration to the second iteration is much less (<0.5
. It can also be seen that radial stress has zero error at the inner and outer radii of the cylinder. Tabular values of tangential stress for different iterations are presented in
Table 6. From the table, it is found that relative error
from first iteration to second iteration is 1.5%We can see from
Table 4,
Table 5 and
Table 6 that we have an overall relative error of less than
including error for radial displacement, radial stress and tangential stress, thus the obtained results of the iterative method are acceptable and are in good agreement with the solution obtained using the element method.
Figure 19,
Figure 20 and
Figure 21 present the behaviour of stresses and the displacement in the FGM cylinder under internal pressure (
MPa)–external pressure (
) MPa for
,
and
iteration steps. Examining
Figure 19 for the radial stresses obtained at different iterations, it can be noted that solution values of radial stress are refined at the second iteration,
. Moreover, in the case of tangential stress,
Figure 20 can be observed to understand the efficacy of the iterative technique as tangential stress values obtained at first and second iterations have an average relative error of
. A similar outcome of the solution refinement can also be observed for displacement solution values from
Figure 4, in which first and second iterations generate displacement values with
average relative error between them.