Next Article in Journal
Detection of Targets in Road Scene Images Enhanced Using Conditional GAN-Based Dehazing Model
Next Article in Special Issue
Analytical Thermal Analysis of Radially Functionally Graded Circular Plates with Coating or Undercoating under Transverse and Radial Temperature Distributions
Previous Article in Journal
Finite Element Simulation and Experimental Investigation of Nanostructuring Burnishing AISI 52100 Steel Using an Inclined Flat Cylindrical Tool
Previous Article in Special Issue
Thermo-Mechanical Buckling and Non-Linear Free Oscillation of Functionally Graded Fiber-Reinforced Composite Laminated (FG-FRCL) Beams
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimization on Elastoplasticity of Functionally Graded Shells of Revolution under Axisymmetric Loading

by
José Simões Moita
1,
Aurélio Lima Araújo
1,
Victor Franco Correia
2,*,
Cristóvão Mota Soares
1 and
José Herskovits
3
1
IDMEC, Instituto Superior Técnico, Universidade de Lisboa, Av. Rovisco Pais, 1049-001 Lisboa, Portugal
2
Escola Superior Náutica Infante D. Henrique, Av. Eng. Bonneville Franco, 2770-058 Paço de Arcos, Portugal
3
COPPE/UFRJ, Universidade Federal do Rio de Janeiro, Rio de Janeiro 21941-909, Brazil
*
Author to whom correspondence should be addressed.
Appl. Sci. 2023, 13(9), 5325; https://doi.org/10.3390/app13095325
Submission received: 10 March 2023 / Revised: 16 April 2023 / Accepted: 20 April 2023 / Published: 24 April 2023
(This article belongs to the Special Issue Latest Advances and Prospects of Functionally Graded Material)

Abstract

:
Depending on the load level, structures can experience a material nonlinearity known as elastoplasticity, which has an important role in the behaviour of structures. In order to avoid the elastoplastic behaviour, it is necessary to find the optimal thickness distribution, which corresponds to the minimum mass that provides an elastic behaviour for a certain load level. The elastoplasticity analysis of functionally graded axisymmetric shells under axisymmetric mechanical loading, and the subsequent optimization, was performed by using a simple conical frustum finite element model with two nodal circles; three degrees of freedom per node, which was based on Kirchhoff’s theory allowing for shear deformation; and using a reduced numerical integration procedure that is essential for its success when applied to thin shells. The formulation accounts for the calculation of the displacements and through-thickness stress distribution, including the effective stress. In this work, the thickness was the design variable in the optimization procedure and the mass was the objective function that needed to be minimized subject to a constraint imposed on the effective stress. The optimization solutions were obtained by using a feasible arc interior point gradient-based algorithm. Some illustrative examples were performed, and the corresponding results are presented and discussed.

1. Introduction

The FGM structures presented in this work were obtained by mixing two different component materials, in this case, a metal and a ceramic, in such a way that a continuous variation over the thickness direction of the material properties was verified, in contrast with the laminated composite material structures. Metal–ceramic FGM plates and shells are widely used in engineering in aircraft, space vehicles, and other engineering applications.
Depending on the load level, structures can experience material nonlinearity known as elastoplasticity, which should be avoided because of large deformations resulting from elastoplastic behaviour in isotropic and FGM structures.
In recent years, many research studies on FGM structures were done. Here, we focused our attention on axisymmetric-shell-related studies: Reddy et al. [1], Li et al. [2], Tran et al. [3] and Thai et al. [4] studied circular plates in linear and geometrically nonlinear behaviour using different theories and to provide different analyses: static bending, dynamic and buckling.
Regarding the structural optimization involving FGM plates and shells of revolution, there is a lack of publications specifically devoted to elastoplastic behaviour.
Elastoplasticity analysis has been less investigated in FGM structures. Zhang and Zhou [5] presented a model for FGM circular plates based on a physical neutral surface and higher-order shear deformation theory. Moita et al. [6] presented a formulation for linear and nonlinear static bending and elastoplasticity analysis of functionally graded axisymmetric plate/shell-type structures under mechanical loading. Moita et al. [7] presented structural and sensitivity analyses, allowing for the material distribution and sizing optimization of functionally graded material structures. The finite element model used in this work was the same as the one previously used by the authors in references [6,7]. Other studies that dealt with the optimization of FGM of different types of structures and objective functions are as follows: Cho et al. [8], Chen and Tong. [9], Xia and Wang [10], Ashjari and Khoshravan [11], Wang et al. [12], Shabana et al. [13] and Franco et al. [14]. In what concerns the elastoplastic behaviour of axisymmetric shells, the work of Popov and Sharifi [15] can be mentioned and was used here for comparison and validation purposes.
In this paper, we present the elastoplastic analysis and subsequent thickness optimization results of functionally graded axisymmetric structures subjected to axisymmetric loading. To avoid the elastoplastic behaviour of the structures, it was necessary to find the minimum thicknesses (which correspond to the minimum mass of the structure) that provide an elastic behaviour for a certain load level. This kind of optimization problem is not encountered in the literature, and therefore, it was the main objective of this work.
The solutions were obtained using a finite element model based on Zienckiewics et al. [16]. The analyses of some illustrative optimization applications were performed, and the results are presented and discussed.

2. Formulation of the P-FGM Model

The FGM structures considered in the present work were made by mixing two distinct isotropic material phases, in this case, a ceramic and a metal. Using the volume fraction of the constituent materials, as proposed by the power law function in Bao and Wang [17], the through-thickness material properties of an FGM plate–shell structure can be obtained, as in Moita et al. [6,7]. The volume fraction of the ceramic and metal phases for each virtual layer k is defined according to the power law
V c k = ( 0.5 + z ¯ h )   p ;   V m k = 1.0 V c k
where z ¯ is the thickness coordinate of the mid-surface of each virtual layer k.
Knowing the volume fractions V c k and V m k and using the rule of mixtures, the material properties, such as Young’s modulus E or mass density ρ, of each virtual layer of an FGM can be determined:
E k = V c k   E c + V m k   E m
These equations show that when p = 0, the plate is fully ceramic and when p = ∞, the plate is fully metal.

3. Displacement and Stress Fields

Figure 1 shows an axisymmetric shell that is subjected to axisymmetric loading. Considering Kirchhoff’s theory, the displacement field is given by
u   ( s , z ) = u 0 z   ( d w / ds )   ;   w   ( s , z ) = w
where u and w are displacements of a generic point in the middle plane with respect to the local axes.
Imposing a constraint C   ( w , β ) = dw / ds β = 0 , hence γ sz = dw / ds β , where β is the rotation of the normal to the middle plane, which can vary independently, one can obtain the transverse shear strain.
For the conical frustum finite element with straight elements, as in Zienckiewics et al. [16], applied to axisymmetric shells under axisymmetric loading, we have R = ,   u θ = 0 and d / d θ = 0. The linear strain components in the local curvilinear system are given as follows:
{ ε s ε θ κ s κ θ γ sz } = { du / ds ( u cos ϕ w sin ϕ ) / r d β / ds ( β cos ϕ ) / r dw / ds β }
The relations between the local (s, z’) and global (r, z) coordinates, as well as the local–global displacements, can be found in Moita et al. [7], and, after this transformation, the strain–displacement relations written in terms of global displacements are
{ ε s ε θ κ s κ θ γ sz } = { ( u r / s ) cos ϕ + ( u z / ds ) sin ϕ u r / r d β / ds ( β cos ϕ ) / r ( du r / ds )   sin ϕ + ( du z / ds )   cos ϕ β }
Equation (3) is then written as
{ u ( s , z ) w ( s , z ) } = [ 1 0 z 0 1 0 ]   { u 0 w β } ;   u = Z   d ,
and the stress–strain relations for each virtual layer k are given by
σ k = Q k   ε k  
where σ k = {   σ s   m   σ θ m   σ s b   σ θ b   τ sz s   } T   and ε k = {   ε s   m   ε θ m   ε s b   ε θ b   γ sz   } T are the stress and strain vectors, respectively, and Q k   is the elasticity matrix for a virtual layer k:
σ k = { σ s m σ θ m σ s b σ θ b τ sz } k = [ Q 1 0 0   0 Q 1 0   0 0 Q 2 ] k   { ε s ε θ z ¯ k   κ s z ¯ k   κ θ γ sz } k   ;   Q 1 k = [ E k 1 ν 2 ν E k 1 ν 2 ν E k 1 ν 2 E k 1 ν 2 ]   ;   Q 2 k = 5 6 E k 2 ( 1 + ν k )
where [ Q 1 ] and [0] are (2 × 2) sub-matrices and z ¯ k   is the middle surface z coordinate of a generic virtual layer k.
The constitutive equation for a virtual layer k is then given by
σ ^ k = D ^ k   ε k
where σ ^ k are the resultant forces and moments and D ^ k is the constitutive matrix.
The functionally graded materials have nonsymmetric properties through the thickness, and for this reason, the bending–stretching coupling exists, and thus, Equation (9) takes the form given in Moita et al. [6]:
{ N M Q } = [ A B 0 B D 0 0 0 A s ]   { ε m ε b ε s }
( A , B , D ) k = Q k   z k - 1 z k ( 1 , z , z 2 )   dz ;   A s k = Q k   z k 1 z k dz

4. Elastoplastic Formulation for FGM Structures

As described in Moita et al. [6], the present study used an extended Tamura–Tomota–Ozawa (TTO) model to describe the elastic–plastic behaviour of ceramic/metal FGM. This extended TTO model uses the stress–strain transfer parameter q * , which depends on the constituent material properties and the microstructural interaction in the FG material, and is given by Jin et al. [18]:
q * = σ c σ m ε c ε m
Using this parameter, the through-thickness variation in Young’s modulus and the yield stress σ Y may be obtained as follows:
E k = ( V m k E m q * + E c q * + E m + V c k E c ) ( V m k q * + E c q * + E m + V c k )
σ Y k = σ Y m   ( V m k + ( q * + E m q * + E c ) E c E m V c k )
With these equations, the variations in Young´s modulus and the yield stress for the case of constituents Si3N4 and SUS304 are shown in Figure 2a,b, respectively.

4.1. Elastoplastic Constitutive Relation

To conduct the elastoplastic analysis, the material was assumed to obey the von Mises yielding criterion, and the corresponding yield surface was assumed. The yield condition can be expressed as given in Nayak and Zienkiewicz [19]:
F ( σ , κ ) = f ( σ ) σ Y ( κ ) = 0
where the yield level σ Y can be a function of the hardening parameter κ .
For the case of the virtual layer approach, for each layer of an FGM structure, and considering the present type of shell and type of loading, the effective stress f (σ, κ ) is given by
f ( σ k ) = σ ¯ k = [ σ s k 2 + σ θ k 2 σ s k σ θ k + 3   τ sn k 2 ]   1 / 2
Furthermore, the effect of the shear stress can be neglected. In this case, we have
f ( σ ) = σ ¯ 2 = ( σ s ) 2 + ( σ θ ) 2 ( σ s   σ θ ) .
Thus, the yield surface for each virtual layer k is given by
F k = f ( σ k ) = σ ¯ k = [ σ s k 2 + σ θ k 2 σ s k σ θ k ]   1 / 2 = σ 0 k
where f ( σ k ) = σ ¯ k and σ 0 k are the effective stress and the constant uniaxial yield stress of a virtual layer, respectively. It should be emphasized that for the FGM, the mechanical properties vary through the thickness because of the changes in V c k (or V m k ) through the thickness.
The previous equations were then modified by separating the membrane and bending stresses to consider the simplest yield condition proposed by Ilyushin [20]:
F k = σ ¯ k = [ ( σ s k m ) 2 + ( σ θ k m ) 2 ( σ s k m ) ( σ θ k m ) + ( σ s k b ) 2 + ( σ θ k b ) 2 ( σ s k b ) ( σ θ k b ) ]     1 / 2     = σ 0 k

4.2. Flow Rule

Following the development of Nayak and Zienckiewics [19], the plastic strain increment is defined as being proportional to the stress gradient of a plastic potential Q, which is taken to be equal to the yield surface condition for an associated flow rule. At the beginning of plasticity, the total increment of strains in elastoplasticity can be calculated by summing up the elastic and the plastic strain components:
d   ε = d   ε e + d   ε p
Introducing the flow vector a = d Q / d σ dF / d σ , the following can be written:
a 1 = F k σ s k m = 2   σ s k m σ θ k m 2 σ ¯ k ;   a 2 = F k σ θ k m = 2   σ θ k m σ s k m 2 σ ¯ k a 3 = F k σ s k b = 2   σ s k b σ θ k b 2 σ ¯ k ;   a 4 = F k σ θ k b = 2   σ θ k b σ s k b 2 σ ¯ k .
For elastic behaviour given in Equation (8), the incremental constitutive elastoplastic relation is then given by
d σ k = Q k ep   d ε k
and the elastoplastic constitutive matrix is given by
Q k ep = Q k Q k   a k   a k T   Q k A + a k T   Q k   a k
where the hardening parameter A = d σ k / d ε p can be obtained from uniaxial testing.
σ k = { σ s m σ θ m σ s b σ θ b τ sz } k = [ Q 1 ep 0 0   0 Q 1 ep 0   0 0 Q 2 ] k   { ε s ε θ z ¯ k   κ s z ¯ k   κ θ γ sz } k
In this study, ideal plasticity (i.e., with no hardening) was considered, leading to A = 0. To satisfy the yield criterion, the state of stress cannot move outside the yield surface. Following the procedure of Owen and Hinton [21], this discrepancy can be practically eliminated by ensuring that after the beginning of elastoplasticity, the load increments considered in the solution are sufficiently small. Moreover, the state of stress can be reduced to the yield surface by scaling the strain vector.

5. Finite Element Approach

As mentioned before, for a conical frustum, which is a straight finite element, as shown in Figure 3, with two nodes and three degrees of freedom per node, the displacements and rotation uri, uri and β i were used in this work.
The simplest interpolations were given by Zienkiewics et al. [8]: N i = ( 1 + ξ i ξ ) / 2 , with ξ i = 1 and ξ = 2   s / L . Thus, N 1 = ( 1 ξ ) / 2 and N 2 = ( 1 + ξ ) / 2 . For numerical integration, reduced integration with one Gauss point was used at the coordinate ξ = 0 , leading to accurate results for thin and thick structures, as validated elsewhere. From previous expressions, N 1 = N 2 = 1 / 2 , dN 1 / d ξ = 1 / 2 , dN 2 / d ξ = 1 / 2 and d ξ / ds = 2 / L .
The displacement field can be represented in matrix form as
u = Z   ( i = 1 2 N i   d i ) = Z   N   a e ;   d = i = 1 2 N i   d i = N   a e
where ae is the element displacement vector, d i = { u r i u z i β i } is the displacement vector of node i and N i are the shape functions. Thus, we have
{ u r u z β } = [ N i 0 0 0 N i 0 0 0 N i ]   { u r i u z i β i }
From this displacement field, the linear strains are as given as follows:
ε = i = 1 2 B i   d i
B i = [ dN i ds cos ϕ dN i ds sin ϕ 0 N i r 0 0 0 0 dN i ds 0 0 N i cos ϕ r dN i ds sin ϕ dN i ds cos ϕ N i ] ;   dN i ds = dN i d ξ d ξ ds = 1 L
where ϕ is the angle represented in Figure 3.

Elastoplastic Analysis

The virtual work principle applied to elastoplasticity analysis is given by
k = 1 N { A e k 1 k k δ   ε k L   Q k e p   ε k L   d z   d A e   } = δ   a T F e x t e k = 1 N   A e k 1 k k δ   ε k L   σ ^ k L   d z   d A e
Integrating through the thickness, Equation (26) can be written in the following form:
A e δ   a T B T D ^ e p   Β   a   d A e = δ   a T F e x t e δ   a T F i n t e
where
K T e = A e B T   D ^ e p   B   dA e ;   F e x t e = A e N T   f   dA e + F c ;   F i n t e = A e B T   σ ^ e   dA e
A e = 2   π   r m   ds = π   r m   L   d ξ = 2   π   r m   L .
Considering all elements in the domain, a system of equilibrium equations is obtained, and by using the Newton–Raphson incremental/iterative (I;i) method, the incremental equilibrium path can be obtained as follows:
K T   Δ q = F e x t F i n t   = ψ
At any increment iteration i, we have
δ q i = ( K T ) 1 ψ i 1 .
At the end of each increment I, Δ q I = Δ q I 1   + δ q i and the residual force vector is very small, i.e.,
ψ <   F ext   .   Δ
where Δ is a predefined tolerance.

6. Optimization

In the present work, the objective of the optimization problem was to find the design variable h (shell thickness) that minimized an objective function Ω (mass of the structure), subjected to a set of constraints (here, the material yield stresses) for a specific applied load level and a specific gradient index p:
m i n   ( h )
subjected to σ ¯ k / σ Y k 1 ,  p = constant and the mass of the structure calculated using
= 1 NE ( k = 1 N m k   ) = 1 NE ( k = 1 N ρ k   h k   2 π r m L )
For each virtual layer k, we have h k = ( h / V L ) at any optimization iteration, where VL = 20 is the number of virtual layers used in this work. More details justifying the use of this value can be found in previous works from the first author.
The sensitivities of the objective function and constraint functions can be evaluated analytically, as was done in the work of Moita et al. [7], or evaluated by using global finite differences. The optimization problem was solved by using a gradient-based algorithm, namely, the feasible arc interior point algorithm developed by Herskovits [22] and Herskovits et al. [23]. In this work, the sensitivities were obtained mainly using the global finite difference method, except otherwise stated.

7. Applications

From the work of Moita et al. [6], where the analysis of elastoplasticity of shells of revolution was previously performed, it was observed that the results obtained using the present model when compared with results from other sources were in excellent agreement and, therefore, we assumed that the validation of the finite element model was previously done. Nevertheless, the accuracy comparison is presented here first for linear analysis in Section 7.1 and later for elastoplasticity in Section 7.4. The remaining subsections show several optimization problems conducted for some illustrative cases, which can be used for benchmarking purposes.

7.1. Linear Analysis of a Clamped Circular FGM Plate under a Uniform Pressure Load

The accuracy of the present finite element model (PM) is illustrated in this numerical application. An FGM clamped circular plate, with thickness h and radius R and made of titanium/zirconium was studied ( E m = 110.25 GPa, ν m = 0.288, Ec = 278.41 GPa, ν c = 0.288). Based on the rule of mixtures, the effective modulus was estimated following Reddy et al. [1]: E ( z ) = ( E m E c )   V m +   E c   with V m k = ( 0.5 z ¯ / h )   p ; V c k = 1 V m k . Table 1 shows the normalized maximum deflections (at the centre of the plate) when the plate was subjected to a uniform pressure q, calculated by w ¯ c = w c ( 64 D c / q R 4 , D c = E c h 3 / 12 ( 1 ν 2 ) . An excellent agreement was observed between the results obtained with the present model and the results obtained using Mindlin theory (Reddy et al. [1]), elasticity theory (Li et al. [2]) and HSDT theory (Tran et al. [3]). The through-thickness distribution of the nondimensional radial stress is presented in Figure 4. A comparison with the solution obtained by Li et al. [2] revealed very good agreement.

7.2. Elastoplastic Optimization of a Circular FGM Plate under a Uniform Pressure Load

A clamped P-FGM circular plate with a radius-to-thickness ratio R/h = 40 was considered, and its geometry was modelled using 10 finite elements. The constituents of the FGM were Si3N4 as the ceramic and SUS304 as the metal, with the following material properties: E c = 322.27 GPa, E m = 207.79 GPa, ν c = 0.24, ν m = 0.318, ρc = 2370 kg/m3 and ρm = 8160 kg/m3. The yield stress for SUS304 was considered to be σ Y m = 400 MPa. The radius of the circular plate was R = 0.5 m, and the initial thickness was h = 12.5 mm. In this numerical example, the gradient index was specified as p = 2.0, and the stress–strain transfer parameter was q* = 60 × 109 Pa. The non-dimensional loads and displacements were given by Q = p 0   R 4 / E m h 4 and wc/h, respectively.
Solving the optimization problem consisted of the minimization of the mass of the plate, subjected to a constraint of the effective stress: min  ϕ =     ( h ) subjected to σ ¯ k σ Y k . The solution obtained for each different index p is shown in Table 2, and Figure 5 gives the load–displacement paths obtained for the case of p = 2. Since the final design corresponds to higher thicknesses, this means that the elastic behaviour condition was not met in the initial design.

7.3. Elastoplastic Optimization of a Toroidal Shell under Internal Radial Pressure

A numerical example with an FGM toroidal shell, as represented in Figure 6, was considered. The toroidal shell was modelled using 24 finite elements in a sequence, starting at the inner point and moving in a clockwise direction (where node 1 was located at the inner point of the circular section of the right side of Figure 6 such that they were numbered in sequence from 1 to 24 in the clockwise direction). The geometry was defined by r = 254 mm, R = 381 mm and initial thickness h = 12.7 mm corresponding to a mass of m = 335.263 kg. The constituents were zirconia and stainless steel, i.e., ZrO2/SUS304, with the following properties at 300 K: E m = 207.79 GPa, ν m = 0.3176, E c = 168.0 GPa, ν c = 0.2978, ρ m = 8160 kg/m3, ρ c = 5700 kg/m3 and q* = 90 × 109 Pa. The inner surface was ceramic and the outer surface was metal. The gradient index p = 1 was considered, and the shell was subjected to an inner radial pressure load. Figure 7 shows the original and deformed shapes in which the displacements were enlarged 100 times. In this analysis, a finite element mesh with 48 elements was used.
The optimization of the shell was performed with the objective of minimizing the mass of the vessel, with a constraint on the effective stress. It was required for the load level of 11.16 MPa that the structure should not display elastoplastic behaviour. The structure could be modelled using one or more groups of finite elements to allow for different thicknesses. For the optimization problem, a finite element mesh with 24 elements was used. Considering a uniform thickness, the results obtained were as follows: minimum thickness h = 13.8 mm, which corresponded to a mass of m = 364.301 kg, and a maximum effective stress equal to 395 MPa, which equalized the yield stress (395 MPa). The deformed shapes (enlarged 100 times), considering the initial thickness h = 12.7 mm and the optimized h = 13.8 mm, are shown in Figure 8, with the plastic zone in yellow, and the load–displacement paths are presented in Figure 9.
Considering two groups of elements, with one group formed by four elements, namely, the ones that plasticize (elements 1, 2, 23 and 24), and another group with 20 elements for the remaining structure, the results obtained were as follows: h1 = 13.92 mm and h2 =11.21 mm, resulting in a mass of m = 300.396 kg. Considering four groups of elements, each one with six elements, as described in Table 3 (where, as described before, node 1 was located at the inner point of the circular section of the right side of Figure 6), different results in the final mass were obtained, as shown in Table 4.

7.4. Elastoplastic Analysis of a Pressure Vessel with a Torispherical End

To compare the capability of the present model, the elastoplastic analysis was carried out by considering the proposed application of Popov and Sharifi [15]: a pressure vessel with torispherical ends made of isotropic material, subjected to a radial pressure load. The dimensions were as follows: radius of curvature of the spherical crown R = 100 in, ϕ0 = 27.89°, cylinder diameter D = 100 in, radius of the toroidal part r = 6 in and thickness t = 0.8 in. The material properties were as follows: E = 30 × 106 lb/in2 and σY = 30 × 103 lb/in2. In Figure 10, the results obtained with the present model and the results obtained by the authors mentioned before are shown. A very good agreement was observed.
The static elastoplastic analysis of an FGM cylindrical pressure vessel with a torispherical end was performed. The constituents were zirconia and stainless steel, i.e., ZrO2/SUS304, at a temperature of 300 K. This structure is represented in Figure 11, where R = Di = 135 mm, t = 1.27 mm, r1 = 10.1 mm, and is subjected to an inner pressure load.
The material properties of the constituents were as follows: E m = 207.7877 GPa,   ν m = 0. 3175, ρ m = 8160 kg/m3, E c = 168.0629 GPa,   ν c = 0.2978 and ρ c = 5700 kg/m3. The yield stress for SUS304 was σ Y m = 400 MPa and the stress–strain transfer parameter was q* = 90 × 109 Pa. The gradient index considered was p = 0.5. The initial mass was m = 0.379 kg. The deformed shape (with uz enlarged 10 times) of the toroidal shell, considering an internal pressure of 2 MPa, is presented in Figure 12.
The optimization was performed once more by considering the objective function, i.e., the mass of the vessel, and the constraint, i.e., the effective stress. It was required that for the load level of 3.75 MPa, the structure should not display elastoplastic behaviour. The results obtained, as presented in Table 5, were as follows: minimum thickness h = 2.0 mm, which corresponded to a mass m = 0.597 kg, and a maximum effective stress equal to 395.8 MPa and a yield stress equal to 395.8 MPa. For index p = 0.5, the deformed shapes (enlarged 10 times), considering h = 1.27 mm and optimized h = 2.0 mm, are presented in Figure 13, with the plastic zone in yellow. The respective load–displacement curves for the apex are shown in Figure 14.

7.5. Elastoplastic Analysis of an FGM Spherical Cap under Radial Pressure

A clamped spherical cap (Figure 15) with a radius of curvature R = 2.286 m, thickness h = 0.0762 m, half-opening angle φ 0 = 35 ° , made of constituents zirconia and aluminium (Em = 70.0 × 109 N/m2, ν m = 0.3, ρm = 3000 kg/m3, Ec = 151.0 × 109 N/m2, ν c = 0.3, ρc = 2700 kg/m3) and subjected to an outer radial pressure load was considered. The yield stress for aluminium was σ Y m = 250 MPa and q* = 80 × 109. First, p = 1.0 was considered and the initial mass was m = 1288.34 kg. Next, the optimization was performed for the load level of 20 MPa, and the respective load–displacement curves are shown in Figure 16. Afterwards, different gradient indexes were considered. For the spherical caps with exponent indexes p = 2.0 and 5.0, this load level was not reached because they collapsed before the uniform radial external pressure of 20 MPa could be achieved, as summarized in Table 6. The spherical cap for p = 2.0 collapsed at 18.90 MPa and uz = 30.2 mm, and for p = 5.0, it collapsed at 17.20 MPa and uz = 30.5 mm.
In addition, we considered the case of the spherical cap divided into two groups of elements to allow for different thicknesses in each group. Group 1 had the elements 1, 2 and 3, starting at the apex, and group 2 had the elements 4, 5, 6 and 7. The gradient index was taken as p = 1.0. Table 7 shows the results obtained, where it is observed that by allowing for different thicknesses in the structure, a better optimal design was reached.

8. Conclusions

The elastoplastic analysis and optimization of axisymmetric shells made of functionally graded material under axisymmetric loading were performed in this work.
The simple finite element model for axisymmetric shells used in this work required only a reduced number of finite elements to model even complex axisymmetric structures, and the discrete approach of virtual layers used to model the continuous variation of the mechanical properties through the thickness was shown to be sufficiently accurate and efficient. These are important requisites of a finite element model to be integrated into a gradient-based optimization procedure. Therefore, we concluded that the present numerical optimization package can be successfully applied for the elastoplastic design of FGM axisymmetric shell structures.
The simple numerical optimization problems for the elastoplastic design of axisymmetric shell structures can be used by researchers for comparison purposes, considering the scarcity of these studies in the literature.
The optimization examples shown in this paper showed that the gradient index had a great influence on the optimal solutions, which was due to the dependency of all the mechanical properties of this parameter (see Table 2 and Table 6). The presented optimization scheme demonstrated its ability to efficiently find a suitable design that avoided the elastoplastic behaviour of the FGM structure and calculate the appropriate thickness for a given p-index value and loading condition.

Author Contributions

Conceptualization, J.S.M.; methodology, J.S.M.; software, J.S.M.; optimization software, J.H., validation, J.S.M. and A.L.A.; formal analysis, all authors; investigation, all authors; resources, all authors; data curation, J.S.M.; writing—original draft preparation, J.S.M.; writing—review and editing, V.F.C.; supervision, C.M.S.; project administration, C.M.S. All authors read and agreed to the published version of the manuscript.

Funding

This research was funded by FCT, Fundação para a Ciência e Tecnologia, through IDMEC under LAETA, project UIDB/50022/2020, and CNPq, CAPES AND FAPERJ, from Brazil.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data supporting the reported results can be obtained from the first author.

Acknowledgments

The authors dedicate this paper to the memory of José Herskovits Norman (1945–2022), who sadly passed away on 20 June 2022. He was a world-leading academic and distinguished expert in the fields of multidisciplinary optimization, non-linear programming, solid mechanics, finite element techniques, advanced composite materials, computational methods and characterization of material identification properties. We extend our condolences to his wife and family; his colleagues at COPPE-UFRJ Federal University of Rio de Janeiro, Brazil; and to all his friends and colleagues all over the world. He was a friend of the authors and influenced their professional careers through his research collaboration and advice.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Reddy, J.N.; Wang, C.M.; Kitiporncha, S. Axisymmetric bending of functionally graded circular and annular plates. Eur. J. Mech. A/Solids 1999, 18, 185–199. [Google Scholar] [CrossRef]
  2. Li, X.Y.; Ding, H.J.; Chen, W.Q. Elasticity solutions for a transversely isotropic functionally graded circular plate subject to an axisymmetric transverse load qrk. Int. J. Solids Struct. 2008, 45, 191–210. [Google Scholar] [CrossRef]
  3. Tran, L.V.; Ferreira, A.J.M.; Nguyen-Xuan, H. Isogeometric analysis of functionally graded plates using higher-order shear deformation theory. Compos. Part B 2013, 51, 368–383. [Google Scholar] [CrossRef]
  4. Thai, C.H.; Kulasegaram, S.; Tran, L.V.; Nguyen-Xuan, H. Generalized shear deformation theory for functionally graded isotropic and sandwich plates based on isogeometric approach. Comput. Struct. 2014, 141, 94–112. [Google Scholar] [CrossRef]
  5. Zhang, D.G.; Zhou, H.M. Nonlinear bending analysis of FGM circular plates based on physical neutral surface and higher-order shear deformation theory. Aerosp. Sci. Technol. 2015, 41, 90–98. [Google Scholar] [CrossRef]
  6. Moita, J.S.; Mota Soares, C.M.; Mota Soares, C.A.; Ferreira, A.M. Elastoplastic and nonlinear analysis of Functionally Graded axisymmetric shell structures under thermal environment, using a conical frustum finite element model. Compos. Struct. 2019, 226, 111186. [Google Scholar] [CrossRef]
  7. Moita, J.S.; Araújo, A.L.; Correia, V.M.; Mota Soares, C.M.; Herskovits, J. Material Distribution and Sizing Optimization of Functionally Graded Plate-Shell Structures. Compos. Part B 2018, 142, 263–272. [Google Scholar] [CrossRef]
  8. Cho, J.R.; Shin, S.W. Material composition optimization for heat-resisting FGMs by artificial neural network. Compos. Appl. Sci. Manuf. 2004, 35, 585–594. [Google Scholar] [CrossRef]
  9. Chen, B.; Tong, L. Thermomechanically coupled sensitivity analysis and design optimization of functionally graded materials. Comput. Methods Appl. Mech. Eng. 2005, 194, 1891–1911. [Google Scholar] [CrossRef]
  10. Xia, Q.; Wang, M.Y. Simultaneous optimization of the material properties and the topology of functionally graded structures. Comput.-Aided Des. 2008, 40, 660–675. [Google Scholar] [CrossRef]
  11. Ashjari, M.; Khoshravan, M.R. Mass optimization of functionally graded plate for mechanical loading in the presence of deflection and stress constraints. Compos. Struct. 2014, 110, 118–132. [Google Scholar] [CrossRef]
  12. Wang, Z.W.; Zhang, Q.; Xia, L.Z.; Wu, J.T.; Liu, P.Q. Stress analysis and parameter optimization of an FGM pressure vessel subjected to thermo-mechanical loading. Proc. Eng. 2015, 130, 374–389. [Google Scholar] [CrossRef]
  13. Shabana, Y.M.; Elsawaf, A.; Khalaf, H.; Khalil, Y. Stresses minimization in functionally graded cylinders using particle swarm optimization technique. Int. J. Press. Vessel. Pip. 2017, 154, 1–10. [Google Scholar] [CrossRef]
  14. Franco Correia, V.M.; Aguilar Madeira, J.F.; Araújo, A.L.; Mota Soares, C.M. Multiobjective optimization of ceramic-metal functionally graded plates using a higher order model. Compos. Struct. 2018, 183, 146–160. [Google Scholar] [CrossRef]
  15. Popov, E.P.; Sharifi, P. A refined curved element for thin shells of revolution. Int. J. Numer. Methods Eng. 1971, 3, 495–508. [Google Scholar] [CrossRef]
  16. Zienkiewics, O.C.; Bauer, J.; Morgan, K.; Onate, E. A simple and efficient element for axisymmetric shells. Int. J. Numer. Methods Eng. 1977, 11, 1545–1558. [Google Scholar] [CrossRef]
  17. Bao, G.; Wang, L. Multiple cracking in functionally graded ceramic/metal coatings. Int. J. Solids Struct. 1995, 32, 2853–2871. [Google Scholar] [CrossRef]
  18. Jin, Z.H.; Paulino, G.H.; Dodds, R.H., Jr. Cohesive fracture modelling of elastic–plastic crack growth in functionally graded materials. Eng. Fract. Mech. 2003, 70, 1885–1912. [Google Scholar] [CrossRef]
  19. Nayak, G.C.; Zienckiewics, O.C. Elastoplastic stress analysis. A generalization for various constitutive relations including strain softening. Int. J. Numer. Methods Eng. 1972, 5, 113–135. [Google Scholar] [CrossRef]
  20. Ilyushin, A.A. Plasticité-Deformations Elasto-Plastic Stress Analysis; Eyrolles: Paris, France, 1956. [Google Scholar]
  21. Owen, D.R.J.; Hinton, E. Finite Element in Plasticity; Pineridge Press Ltd.: Swansea, Wales, 1980. [Google Scholar]
  22. Herskovits, J. Feasible Direction Interior-Point Technique for Nonlinear Optimization. J. Optim. Theory Appl. 1998, 99, 121–146. [Google Scholar] [CrossRef]
  23. Herskovits, J.; Mappa, P.; Goulart, E.; Mota Soares, C.M. Mathematical programming models, and algorithms for engineering design optmization. Comput. Methods Appl. Mech. Eng. 2005, 194, 3244–3268. [Google Scholar] [CrossRef]
Figure 1. Axisymmetric shell.
Figure 1. Axisymmetric shell.
Applsci 13 05325 g001
Figure 2. (a) Variation in Young’s modulus through the thickness; (b) variation of the yield stress through the thickness.
Figure 2. (a) Variation in Young’s modulus through the thickness; (b) variation of the yield stress through the thickness.
Applsci 13 05325 g002
Figure 3. Conical frustum finite element.
Figure 3. Conical frustum finite element.
Applsci 13 05325 g003
Figure 4. Through-thickness distribution of the nondimensional stress for different p-indexes. Li et al. [2].
Figure 4. Through-thickness distribution of the nondimensional stress for different p-indexes. Li et al. [2].
Applsci 13 05325 g004
Figure 5. Load–displacement paths for the original and optimized clamped FGM plate.
Figure 5. Load–displacement paths for the original and optimized clamped FGM plate.
Applsci 13 05325 g005
Figure 6. Toroidal shell with thickness h subjected to internal pressure p.
Figure 6. Toroidal shell with thickness h subjected to internal pressure p.
Applsci 13 05325 g006
Figure 7. Deformed shapes of an FGM toroidal shell for index p = 1.
Figure 7. Deformed shapes of an FGM toroidal shell for index p = 1.
Applsci 13 05325 g007
Figure 8. Deformed and optimized shapes of an FGM toroidal shell (index p = 1).
Figure 8. Deformed and optimized shapes of an FGM toroidal shell (index p = 1).
Applsci 13 05325 g008
Figure 9. Load–displacement paths for the original and optimized toroidal shell.
Figure 9. Load–displacement paths for the original and optimized toroidal shell.
Applsci 13 05325 g009
Figure 10. Load–displacement paths of an isotropic pressure vessel.
Figure 10. Load–displacement paths of an isotropic pressure vessel.
Applsci 13 05325 g010
Figure 11. Pressure vessel with a torispherical end.
Figure 11. Pressure vessel with a torispherical end.
Applsci 13 05325 g011
Figure 12. Initial and deformed shapes of the vessel.
Figure 12. Initial and deformed shapes of the vessel.
Applsci 13 05325 g012
Figure 13. Initial and optimized deformed shapes of the vessel.
Figure 13. Initial and optimized deformed shapes of the vessel.
Applsci 13 05325 g013
Figure 14. Load–displacement paths for the original and optimized vessel.
Figure 14. Load–displacement paths for the original and optimized vessel.
Applsci 13 05325 g014
Figure 15. Spherical cap under radial pressure.
Figure 15. Spherical cap under radial pressure.
Applsci 13 05325 g015
Figure 16. Load–displacement paths for the original and optimized cap.
Figure 16. Load–displacement paths for the original and optimized cap.
Applsci 13 05325 g016
Table 1. Dimensionless central deflection w ¯ c of an FGM clamped circular plate.
Table 1. Dimensionless central deflection w ¯ c of an FGM clamped circular plate.
h/RSourcePower Index p
0241050100Metal
0.05Li et al. [2]2.5611.4051.2841.1571.0491.0321.015
Reddy et al. [1]2.5541.4021.2821.1551.0461.0291.011
Tran et al. [3]2.5481.3991.2791.1521.0441.0271.009
PM2.5601.4041.2831.1551.0431.0241.015
0.1Li et al. [2]2.6671.4561.3291.2011.0911.0741.057
Reddy et al. [1]2.6391.4441.3201.1901.0801.0631.045
Tran et al. [3]2.6301.4391.3141.1861.0761.0591.042
PM2.6431.4461.3211.1891.0761.0581.046
Table 2. Optimization of a clamped plate.
Table 2. Optimization of a clamped plate.
Initial DesignFinal Design
pThickness (mm)Mass (kg)Thickness (mm)Mass (kg)
5.0h = 12.570.175h = 17.1996.504
2.0h = 12.560.782h = 17.0983.101
2.0h = 12.560.782h = 17.08 *83.052
1.0h = 12.551.394h = 16.8369.198
0.5h = 12.541.983h = 16.4555.250
0.2h = 12.532.574h = 15.9341.513
*—calculated by using analytical sensitivities.
Table 3. Four groups of elements.
Table 3. Four groups of elements.
Group NumberElements in the Group
11 2 3 22 23 24
24 5 6 19 20 21
37 8 9 16 17 18
410 11 12 13 14 15
Table 4. Optimization of a toroidal shell.
Table 4. Optimization of a toroidal shell.
Initial DesignFinal Design
Number of
Groups
p Index and
Thickness (mm)
Mass (kg)Thickness (mm)Mass (kg)
1p = 1; h = 12.70335.263h = 13.80364.301
1p = 1; h = 12.70335.263h = 13.88 *366.402
2p = 1; h = 12.70335.263h1 = 13.92
h2 = 11.21
300.396
4p = 1; h = 12.70335.263h1 = 13.86
h2 = 9.37
h3 = 7.51
h4 = 7.32
222.418
*—calculated by using analytical sensitivities.
Table 5. Optimization of a pressure vessel with a torispherical end.
Table 5. Optimization of a pressure vessel with a torispherical end.
Initial DesignFinal Design
pThickness (mm)Mass(kg)Thickness (mm)Mass (kg)
0.5h = 1.270.379h = 2.00.597
Table 6. Optimization of the spherical cap.
Table 6. Optimization of the spherical cap.
Initial DesignFinal Design
p Index and
Thickness (mm)
Mass(kg)Thickness (mm)Mass (kg)
p = 1.0; h = 76.21288.340h = 110.01859.808
p = 1.0; h = 76.21288.340h = 110.7 *1871.643
p = 0.5; h = 76.21311.028h = 97.61679.216
p = 2.0; h = 76.21265.709h = 1262092.905
p = 5.0; h = 76.21243.065h = 1402283.845
*—calculated by using analytical sensitivities.
Table 7. Optimization of a spherical cap (index p = 1.0).
Table 7. Optimization of a spherical cap (index p = 1.0).
Initial DesignFinal Design
No. of
Groups
p Index and
Thickness (mm)
Mass(kg)Thickness (mm)Mass (kg)
1p = 1.0; h = 76.21288.340h = 1101859.808
2p = 1.0; h = 76.21288.340h1 = 110.6
h2 = 97.4
1688.825
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Simões Moita, J.; Araújo, A.L.; Franco Correia, V.; Soares, C.M.; Herskovits, J. Optimization on Elastoplasticity of Functionally Graded Shells of Revolution under Axisymmetric Loading. Appl. Sci. 2023, 13, 5325. https://doi.org/10.3390/app13095325

AMA Style

Simões Moita J, Araújo AL, Franco Correia V, Soares CM, Herskovits J. Optimization on Elastoplasticity of Functionally Graded Shells of Revolution under Axisymmetric Loading. Applied Sciences. 2023; 13(9):5325. https://doi.org/10.3390/app13095325

Chicago/Turabian Style

Simões Moita, José, Aurélio Lima Araújo, Victor Franco Correia, Cristóvão Mota Soares, and José Herskovits. 2023. "Optimization on Elastoplasticity of Functionally Graded Shells of Revolution under Axisymmetric Loading" Applied Sciences 13, no. 9: 5325. https://doi.org/10.3390/app13095325

APA Style

Simões Moita, J., Araújo, A. L., Franco Correia, V., Soares, C. M., & Herskovits, J. (2023). Optimization on Elastoplasticity of Functionally Graded Shells of Revolution under Axisymmetric Loading. Applied Sciences, 13(9), 5325. https://doi.org/10.3390/app13095325

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop