Home Investigation of porosity effect on flexural analysis of doubly curved FGM conoids
Article Open Access

Investigation of porosity effect on flexural analysis of doubly curved FGM conoids

  • Md Irfan Ansari , Ajay Kumar EMAIL logo , Danuta Barnat-Hunek , Zbigniew Suchorab and Bartłomiej Kwiatkowski
Published/Copyright: October 1, 2019

Abstract

The flexural analysis of doubly curved functionally graded porous conoids was performed in the present paper. The porosities inside functionally graded materials (FGMs) can occur during the fabrication and lead to the occurrence of micro-voids in the materials. The mathematical model includes expansion of Taylor’s series up to the third degree in thickness coordinate and normal curvatures in in-plane displacement fields. Since there is a parabolic variation in transverse shear strain deformation across the thickness coordinate, the shear correction factor is not necessary. The condition of zero-transverse shear strain at upper and lower surface of conoidal shell is implemented in the present model. The improvement in the 2D mathematical model enables to solve problems of moderately thick FGM porous conoids. The distinguishing feature of the present shell from the other shells is that maximum transverse deflection does not occur at its centre. The improved mathematical model was implemented in finite element code written in FORTRAN. The obtained numerical results were compared with the results available in the literature. Once validated, the current model was employed to study the effect of porosity, boundary condition, volume fraction index, loading pattern and others geometric parameters.

1 Introduction

From last two decades, shell structures made of functionally graded materials (FGM) broadly used by civil, mechanical, aeronautical and marine engineer. The FGM is an inhomogeneous material, composed of two (or more) materials, organized with a view to having a smooth gradation in the desired direction. However, at the time of fabrication, porosities are infused in the FGM that create major issue. Apart from this, the FGM manufactured using sintering process retains porosities because of different solidification rate of material constituents. Koizumi [1] used FGM in advanced engineering structures experiencing elevated temperatures. Eslami et al. [2] developed a general solution for the 1-D steady-state mechanical and thermal stresses in an FGM hollow thick sphere. The static response of functionally graded cylindrical shells using the element-free kp-Ritz method was analysed by Zhao et al. [3]. An elastic solution for a sandwich panel with isotropic skins was presented by Kashtalyan and Menshykova [4]. Ameur et al. [5] developed a trigonometric shear deformation theory containing four unknowns for the bending analysis of an FGM plate. Natarajan and Manickam [6] considered a realistic variation in displacement along the thickness direction in higher order model for the bending and flexure vibration behaviour of an FGM sandwich plate. Ajay et al. [7] developed C0 finite element based on higher order shear deformation theory for the analysis of laminated composite rhombic shell. Bessaim et al. [8] developed five unknown-based shear deformation theory for the static analysis of FGM sandwich plate. The stresses in functionally graded doubly curved shells were calculated by incorporating the differential quadrature method in the first order shear deformation theory by Tornabene and Viola [9]. An analytical solution was developed by Sayyad and Ghugal [10] to account for the effect of transverse shear and transverse normal for the bending analysis of isotropic, laminated composite and sandwich plates. Viola et al. [11] used unconstrained third-order shear deformation theory for the static analysis of moderately thick functionally graded conical shells. Asemi et al. [12] used the classical theory for a linear thermoelastic analysis of thick truncated cone. Dai and Dai [13] developed

an analytical solution based upon the classical shell theory for FGM cylindrical shell under thermomechanical loading. Xiang and Liu [14] used a meshless global collocation method with nth-order shear deformation theory for the static analysis of FGM sandwich plate. Parihar et al. [15] pointed out the recent advancement in the manufacturing process of FGM. Alipour and Shariyat [16] conducted the free vibration analysis of sandwich plates with isotropic/orthotropic face sheet and different combinations of boundary conditions.

Since casting and fabrication of conoids is easy due to its singly ruled surface; hence, it is favoured in the construction industry. Conoids are structurally stiff, aesthetically appealing, and are used to cover the column-free large area in industrial structures, aircraft hangars and exhibition hall. The first study on simply supported and clamped conoids was carried out by Hadid [17] in 1964. A combined variational approach was used for bending analysis of elastic conoids for both types of boundary condition. Finite difference method was carried out on the conoidal shell by Das and Bandyopadhyay [18] for both experimental as well as theoretical investigation. Ghosh and Bandyopadhyay [19] engaged their own formulation to examine the influence of cutouts on the static analysis of conoidal shells. The bending analysis of stiffened conoids is studied by Das and Chakravorty [20] using three noded beam element. Kumari and Chakravorty [21] used FSDT for the study on the bending response of delaminated conoids. The FE method was used by Bakshi and Chakravorty [22] for the analysis of first ply failure occurs in laminated composite conoidal shell subjected to uniformly distributed loading. Malekzadeh Fard and Baghestani [23] explored the free vibration behavior of moderately thick doubly curved shell based on FSDT with classic boundary conditions. The bending and free vibration analysis of functionally graded porous plate was presented by Akbas [24]. He used FSDT model and Navier solution for solving the problem. An analytical model was developed by Al Rjoub and Hamad [25] to study the effect of porosity on free vibration analysis of porous beam. Eltaher et al. [26] developed a modified porosity model for the bending and vibration analysis of functionally graded nanobeams and mathematical model was solved using a finite element method. Kiran et al. [27] presented the effect of porosity on the static response of FGM magneto-electro-elastic skew plate. Gupta and Talha [28, 29, 30] explored the influence of porosities on the flexure and free vibration response of graded plate using higher deformation theory.

The literature review reveals that no results for the bending analysis of doubly curved FGM porous conoidal shell are available. Therefore, an attempt was in the present paper to study the bending behaviour of FGM porous conoidal shell with the help of improved mathematical model. The FE coding was done by using a C0 nine noded FE with seven nodal unknowns at each node for the present improved mathematical model developed by the authors. C1 continuity requirement associated with the present model is suitably circumvented. The present study facilitates the bending analysis by finite element (FE)modelling, keeping in mind the processing time in computer and simplicity of approach. This study is the first step towards enhancing our understanding of the bending problem of FGM porous conoidal shell.

2 Formulation

The FGM porous conoidal shell of sides a, b and thickness h are depicted in Fig. 1 (a). The upper layer (z = +h/2) of the porous shell surface is ceramic rich, while the bottom portion (z = −h/2) of the porous shell surface is metal rich, as shown in Fig. 1(b)., with a gradation zone having a smooth variation of material properties in between the two surfaces. Using the power law, the final properties of the FGM porous conoidal shell at any point in thickness coordinate (z) can be stated as

Figure 1 (a) Geometry of conoidal shell and (b) porosities models for even porosity in thickness coordinate
Figure 1

(a) Geometry of conoidal shell and (b) porosities models for even porosity in thickness coordinate

(1)Pz=PcVcz+PmVmzPc+Pmβ2

where,

(2)Vcz=12+zhn,0n,Vcz+Vmz=1

where β (β << 1) is the volume fraction of porosities. P(z) implies material properties like the Young’s modulus of elasticity (E), material density (ρ) and Poisson’s ratio (ν) of the FGM porous conoidal shell.

Once the volume fraction index of material constituents is known, the Young’s modulus at any point ‘z’ can be calculated as per the rule of mixture.

(3)E(z)=EcVc(z)+EmVm(z)Ec+Emβ2

where Ec, Vc and Em, Vm are the material properties and volume fraction of ceramic and metal, respectively, and n is the volume-fraction index.

2.1 Theory and formulation

2.1.1 Displacement field and strains

In order to develope the mathematical model, the displacement fields for FGM shell are considered on the basis of the third order shell theory and are as follows:

(4)u(x,y,z)=1+zRxu0(x,y)+zθx(x,y)+z2ξx(x,y)+z3ζx(x,y)v(x,y,z)=1+zRyv0(x,y)+zθy(x,y)+z2ξy(x,y)+z3ζy(x,y)w(x,y,z)=w0(x,y)

where the middle section is taken as a reference for material coordinates (x, y, z). u, v, w denoted the displacements of along the (x, y, z) coordinates, u0, v0, w0 are corresponding displacements on the mid-plane. θx,θy are the bending rotations normal to the mid-plane about the y axis and x-axes, respectively. The functions ξx, ξy , ζx and ζy are higher order terms of Taylor’s series expansion at the midplane of the conoidal shell. By imposing the boundary condition (zero transverse shear strain at top and bottom) in Eq. (4), the function ξx, ξy , ζx and ζy will be calculated as

(5)ξx=ξy=0,ζx=43h2θx+wxv0Rxy,ζy=43h2θy+wyu0Rxy

Replacing the unknown in Eq. (5) into Eq. (4), we obtain

(6)u(x,y,z)=1+zRxu0+θxz4z33h2w0x4z33h2+v04z33h2Rxyv(x,y,z)=1+zRyv0+θyz4z33h2w0y4z33h2+u04z33h2Rxyw(x,y,z)=w0

For omitting C1 continuity problem associated with TSDT, the out of plane derivatives are replaced by the following relations

(7)ψx=w0x,ψy=w0y

The final form of displacement filed expression owning C0 continuity is written as

(8)u(x,y,z)=1+zRxu0+θxz4z33h2ψx4z33h2+v04z33h2Rxyv(x,y,z)=1+zRyv0+θyz4z33h2ψy4z33h2+u04z33h2Rxyw(x,y,z)=w0

Hence, the field variables per node taken in the present investigation are u0, v0, w0, θx, θy, Ψ x and Ψ y. Mathematically, it may be expressed as

(9)δ=u0,v0,w0,θx,θy,ψx,ψyT

where {δ} is termed as a displacement vector.

The strain vector may be expressed as

(10)ε=εx,εy,εz,γxy,γxz,γyzT

Further, the strain vector {ε} can be associated with global displacement vector {X} by the subsequent relationship

(11)ε=BX

Here strain-displacement matrix [B] contains the derivatives of shape functions.

The in-plane and transverse shear strains are

(12)εx=ux+wRxεy=vy+wRyγxy=vx+uy+2wRxyγxz=uz+wxu0Rxv0Rxyγyz=vz+wyv0Ryu0Rxy

The strain relationship can be written as:

(13)εxεyγxy=εx0εy0γxy0+zk11k21k614z33h2k13k23k63γxzγyz=γxz0γyz04z2h2k42k52

where

(14)εx0=u0x+w0RX,εy0=v0y+w0Ry,γxy0=v0x+u0y+2w0Rxyγxz0=w0x+θxv0Rxy,γyz0=w0y+θyu0Rxyk11=θxx+u0x1Rx,k21=θyy+v0y1Ryk61=θyx+θxy+u0y1Rx+v0x1Ryk13=θxx+ψxxv0x1Rxy,k23=θyy+ψyyu0y1Rxyk63=θyx+ψyx+θxy+ψxy+u0x+v0y1Rxyk42=θx+ψxv0Rxy,k52=θy+ψyu0Rxy

2.1.2 Constitutive relationship

The linear stress-strain constitutive relationship for the FGM conoids are

(15)σ=Qε

where the constitutive matrix

(16)Q=Q11Q12000Q21Q2200000Q6600000Q4400000Q55

where the term Qij can be obtained with the help of the Young’s modulus (E) and Poisson’s ratio (ν), in which E is the function of thickness coordinate.

(17)Q11=Q22=Ez1ν2,Q12=Q21=νEz1ν2,Q66=Q44=Q55=Ez21+ν

3 Finite element modelling

3.1 Introduction

For the present 5C0 finite element (FE) model, nine noded isoparametric Lagrangian element with seven degrees of freedom at each node is utilized in the present investigation. The shape function (interpolation function) is used to express the generalized displacement vector and element geometry at any point within an element as

(18)δ=i=19Niξ,ηδix=i=19Niξ,ηxiy=i=19Niξ,ηyi

The shape functions Ni of nine noded isoprametric Lagrangian element are depicted below.

For corner nodes:

N1=14ξ1η1ξη,
N3=14ξ+1η1ξη,
N7=14ξ1η+1ξη
N9=14ξ+1η+1ξη

For mid nodes:

N2=121ξ2η2η,
N4=12ξ2ξ1η2,
N6=12ξ2+ξ1η2
N8=121ξ2η2+η

For centre node:

N5=1ξ21η2

3.2 Governing equation

The strain energy may be expressed as

(19)U=12{ϵ}T{σ}dxdydz

By using the Eq. (15), the above expression can be represented as

(20)U=12{ϵ}T[D]{ϵ}dxdy

where, [D] = ∫ [H]T [Q] [H] dx3 in which [H] is the matrix that contains the terms involving z and h.

By utilizing Eq. (11) the stiffness matrix [K] is written as

(21)K=[B]TDBdxdy

4 Results and discussion

In this section, the bending analysis of FGM porous conoidal shells was analysed under various type of mechanical loading. Parameters like thickness ratio, aspect ratio, hl/hh ratio, porosities and volume fraction index are also accomplished for numerical results. Unless stated otherwise, the following non-dimensional factors are utilized in the upcoming examples.

Non-dimensional quantities used are

w¯=100wh3Ecq0a4,σ¯x=σxa2,b2,h2100h2q0a2
(22)σ¯y=σya2,b2,h2100h2q0a2,τ¯xy=τxy0,0,h2100h2q0a2
τ¯xz=τxz0,a2,0100h2q0a2,τ¯yz=τyza2,0,0100h2q0a2

The loading patterns used in the present analysis are written below:

(23)q=q0,q=q0sinπxasinπyb,q=q0cosπxasinπybq=q0cosπxacosπyb

The details of some of the boundary condition (BC) used are as follows:

  1. Simply supported (SSSS):

    At x = 0, a v = w = θy = ψ y = 0

    At y = 0, b u = w = θx = ψ x = 0

  2. Clamped (CCCC):

    At x = 0, a and y = 0, b u = v = w = θx = θy = ψ x = ψ y = 0

  3. Clamped and simply supported (CCSS):

    At x = 0, a u = v = w = θx = θy = ψ x = ψ y = 0

    At y = 0, b u = w = θx = ψ x = 0

4.1 Comparison and convergence study

Three appropriate examples have been solved to check the consistency and the stability of present FE results. No results for FGM porous conoidal shellwere found in the literature; hence, the effectiveness of the present formulation has been tested by relating dimensionless deflection for FGM porous plate and isotropic conoidal shell. The material properties of the FGM components used in the present study are as follows:

FGM-1 (Al/ZrO2): Ec = 151 GPa, Em = 70 GPa, νc = νm = 0.3,

ρc = 3000 kg/m3, ρm = 2707 kg/m3

FGM-2 (Al/Al2O3): Ec = 380GPa, Em = 70 GPa, νc = νm = 0.3,

ρc = 3800 kg/m3, ρm = 2707 kg/m3

Example 1. The convergence study of the dimensionless maximum deflection of FGM-1 (Al/ZrO2) plate was conducted for four power law indices (n = 0, 0.5, 1,∞) and presented in Table 1. For this particular example, we used hl = 0 and hh = 0 in the present FE code of the FGM-1 conoidal shell for converting it into FGM plate. It was found that at 16 × 16 mesh size, the results converge for the present nine noded isoparametric elements. From Table 1, it was noted that our numerical results are consistent with Ferreira et al. [31].

Table 1

Comparison of non-dimensional maximum deflection at the centroid of simply supported FGM plate.

Volume fraction index
MeshCeramic0.51Metal
6x60.024800.033290.036800.05350
8x80.024810.033310.036850.05352
12x120.024820.033330.036880.05353
16x160.024820.033330.036880.05353
Ferreira et al. [31]0.02480.03300.03680.0536

Example 2. Fig. 2 represents the comparison of dimensionless maximum transverse displacement of simply supported FGM porous plate. The bottom layer of FGM plate is metallic and the top layer is ceramic. For various value of volume fraction indices, the results were compared with Akbaş [24] and again it showed a decent agreement.

Figure 2 Comparison of non-dimensional deflection of FGM-1 porous plate subjected to uniform loading.
Figure 2

Comparison of non-dimensional deflection of FGM-1 porous plate subjected to uniform loading.

Table 2

Comparison of deflection (x10−2) of simply supported isotropic conoid under uniformly distributed load along y/b = 0.5.

x1/aBakshi and Chakravorty [22]Hadid [17]Present
0.102.36802.52312.5621
0.405.15574.63334.7609
0.604.05103.54483.8575
0.703.41293.31493.2970
0.802.48092.55442.5949

Example 3. No results for FGM porous conoidal shell were available in literature; thus, the presented FE result of FGM porous conoidal shell were compared with numerical result of isotropic conoidal shell. The Poisson’s ratio ν = 0.15, hl/hh = 0.5 and side-to-thickness ratio a/h = 19 was used to validate present result. Table 4 shows the validation of present FE formulation with Hadid [17] and Bakshi and Chakravorty [22]. The validation study confirms the improvement of the present FE result over Bakshi and Chakravorty [22] as the current results are closer to the elasticity results presented by Hadid [17].

Table 3

Variation of non-dimensional deflection of FGM porous conoidal shell subjected to uniform loading.

nβBoundary condition
SSSSCCCCCCSSCSCSCCFFCFCF
001.29810.88571.00191.02072.606329.1371
(0.40,0.50)(0.45,0.50)(0.45,0.50)(0.50,0.55)(0.475,0.0)(1,1)
0.11.42180.96841.09531.11682.795331.2321
(0.40,0.50)(0.45,0.50)(0.45,0.50)(0.50,0.55)(0.475,0.0)(1,1)
0.21.56421.06361.20281.22743.014633.6526
(0.40,0.50)(0.45,0.50)(0.45,0.50)(0.50,0.55)(0.475,0.0)(1,1)
0.201.42061.01831.11851.14513.029634.3131
(0.40,0.50)(0.45,0.50)(0.45,0.50)(0.475,0.55)(0.475,0.0)(1,1)
0.11.57051.12791.23531.26663.294337.3592
(0.40,0.50)(0.45,0.50)(0.45,0.50)(0.475,0.55)(0.475,0.0)(1,1)
0.21.74601.25751.37261.40953.612640.9732
(0.40,0.50)(0.45,0.50)(0.45,0.50)(0.475,0.55)(0.475,0.0)(1,1)
0.501.53101.12321.21741.25003.343337.9973
(0.40,0.50)(0.45,0.50)(0.45,0.50)(0.475,0.55)(0.475,0.0)(1,1)
0.11.70641.25741.35671.39573.676241.8726
(0.40,0.50)(0.45,0.50)(0.45,0.50)(0.475,0.55)(0.475,0.0)(1,1)
0.21.91641.42031.52441.57144.088546.6175
(0.40,0.50)(0.45,0.50)(0.45,0.50)(0.475,0.55)(0.475,0.0)(1,1)
101.66781.24461.33651.37553.694342.0038
(0.40,0.50)(0.45,0.50)(0.45,0.50)(0.475,0.55)(0.475,0.0)(1,1)
0.11.87811.41101.50621.55374.115846.9387
(0.40,0.50)(0.45,0.50)(0.45,0.50)(0.475,0.55)(0.475,0.0)(1,1)
0.22.13681.62001.71671.77514.657053.2180
(0.40,0.50)(0.45,0.50)(0.45,0.50)(0.475,0.55)(0.475,0.0)(1,1)
1002.25171.61421.77561.82754.653851.6012
(0.40,0.50)(0.475,0.50)(0.45,0.50)(0.50,0.55)(0.475,0.0)(1,1)
0.12.62801.89282.07252.13915.340159.2021
(0.40,0.50)(0.475,0.50)(0.45,0.50)(0.50,0.55)(0.475,0.0)(1,1)
0.23.13142.27222.47272.56046.285169.5130
(0.40,0.50)(0.475,0.50)(0.45,0.50)(0.50,0.55)(0.475,0.0)(1,1)
02.80021.91072.16132.20185.622162.8531
(0.40,0.50)(0.45,0.50)(0.45,0.50)(0.50,0.55)(0.475,0.0)(1,1)
0.13.38672.30822.60922.66026.676375.0419
(0.40,0.50)(0.45,0.50)(0.45,0.50)(0.50,0.55)(0.475,0.0)(1,1)
0.24.20942.86223.23673.30308.112490.5591
(0.40,0.50)(0.45,0.50)(0.45,0.50)(0.50,0.55)(0.475,0.0)(1,1)
Table 4

Variation of non-dimensional deflection of FGM porous conoidal shell subjected to sin-sin loading.

nβBoundary condition
SSSSCCCCCCSSCSCSCCFFCFCF
000.88030.66280.72510.73281.34249.8687
(0.425,0.50)(0.475,0.50)(0.475,0.50)(0.50,0.525)(0.475,0.50)(1,1)
0.10.96350.72440.79240.80141.460210.5537
(0.425,0.50)(0.475,0.50)(0.475,0.50)(0.50,0.525)(0.475,0.50)(1,1)
0.21.05930.79540.86980.88031.596111.3454
(0.425,0.50)(0.475,0.50)(0.475,0.50)(0.50,0.525)(0.475,0.50)(1,1)
0.20
0.97320.76310.81600.82901.552711.5935
(0.425,0.50)(0.475,0.50)(0.475,0.50)(0.50,0.525)(0.475,0.50)(1,1)
0.11.07600.84510.90160.91711.712012.5893
(0.425,0.50)(0.475,0.50)(0.475,0.50)(0.50,0.525)(0.475,0.50)(1,1)
0.21.19670.94211.00251.02111.900813.7711
(0.425,0.50)(0.475,0.50)(0.475,0.50)(0.50,0.525)(0.475,0.50)(1,1)
0.501.05430.84190.89130.90811.712712.8317
(0.425,0.50)(0.475,0.50)(0.475,0.50)(0.50,0.525)(0.475,0.50)(1,1)
0.11.17610.94240.99431.01471.909214.1007
(0.425,0.50)(0.475,0.50)(0.475,0.50)(0.50,0.525)(0.475,0.50)(1,1)
0.21.32251.06461.11871.14372.148615.6543
(0.425,0.50)(0.475,0.50)(0.475,0.50)(0.475,0.525)(0.475,0.50)(1,1)
101.15320.93280.98081.00161.894314.1854
(0.425,0.50)(0.475,0.50)(0.475,0.50)(0.50,0.525)(0.475,0.50)(1,1)
0.11.30061.05751.10701.13282.139215.8053
(0.425,0.50)(0.475,0.50)(0.475,0.50)(0.50,0.525)(0.475,0.50)(1,1)
0.21.48311.21431.26441.29662.448217.8653
(0.425,0.50)(0.475,0.50)(0.475,0.50)(0.475,0.525)(0.475,0.50)(1,1)
1001.54411.20571.29181.31852.415617.5095
(0.45,0.50)(0.475,0.50)(0.475,0.50)(0.50,0.525)(0.475,0.50)(1,1)
0.11.80381.41311.50871.54372.815120.0428
(0.45,0.50)(0.475,0.50)(0.475,0.50)(0.50,0.525)(0.475,0.50)(1,1)
0.22.15291.69561.80191.84913.358723.4823
(0.45,0.50)(0.475,0.50)(0.475,0.50)(0.50,0.525)(0.475,0.50)(1,1)
Metal01.89891.42971.56411.58082.895821.2883
(0.425,0.50)(0.475,0.50)(0.475,0.50)(0.50,0.525)(0.475,0.50)(1,1)
0.12.29561.72701.88811.90963.490225.3677
(0.425,0.50)(0.475,0.50)(0.475,0.50)(0.50,0.525)(0.475,0.50)(1,1)
0.22.85062.14042.34072.36894.295230.5302
(0.425,0.50)(0.475,0.50)(0.475,0.50)(0.50,0.525)(0.475,0.50)(1,1)

4.2 Results and discussion

In order to analyze FGM conoidal shells under various type of transverse loading, different composition of material constituents, different combination of boundary condition, numerous value of volume fraction index, side-to-thickness ratio, aspect ratio and hl/hh ratio are considered.

Tables 3 and 4 represent the dimensionless maximum deflection and their location subjected to uniform and sin-sin loading for the FGM-1 porous conoidal shell, respectively.

The numerical results were computed for a/h = 10, a/b = 1, hl/hh = 0.25 (hl = 0.05, hh = 0.2), and SSSS, CCCC, CCSS, CSCS, CCFF and CFCF boundary condition. It can be observed from that the maximum deflection of the conoidal shell occurs at a different location, unlike plate structure where maximum deflection occurs at the midpoint. These numerical values revealed that the dimensionless deflection of porous conoidal shell increases along with porosity percentage because of the decrease in the equivalent modulus of elasticity. Apart from that as we move from n = 0 to n = ∞, the dimensionless deflection of the porous conoidal shell is increased and it may be attributed to the higher volume fraction index which leads to lesser ceramic content, thus reducing its stiffness. Intrestingly, around 200% increase in the maximum dimensionless deflection was noticed as the volume fraction index changed from 0 to ∞ and almost 30% growths were displayed when the volume fraction index changed from 0 to 1 for all boundary conditions. Remarkably, CFCF boundary condition yields the highest deflection while CCCC retains the lowest deflection value among all considered boundary conditions. The influence of thickness ratio on dimensionless deflection of FGM-1 porous conoidal shell is tabulated in Tables 5 and 6 subjected to uniform and sin-sin loading, respectively. These numerical results highlighted that non-dimensional deflection increases along with the thickness of porous conoidal shell. The correlation between hl/hh ratio and dimensionless deflection of FGM-1 porous conoids subjected to uniform loading is tested in Table 7. Interestingly, for higher value of hl/hh ratio, lower value of dimensionless deflection is noticed. The decrease in hl/hh ratio reduces the curvature of the lower end (hl) of the conoidal shell, due to this the stiffness of shell reduces, therefore deflection increased. Fig. 3 depicts the variation of dimensionless deflection along the centre line subjected to uniform loading. The numerical results of dimensionless deflection of FGM-1 porous conoids subjected to various type of transverse load were plotted in Fig. 4. a/h = 10, a/b = 1 and hl/hh = 0.25 was used. The maximum and minimum non-dimensional deflections were found under uniform loading and cos-cos loading respectively. In Fig. 5, the maximum dimensionless deflection of FGM-1 and FGM-2 porous conoidal shellwas plotted against the porosity coefficient. With the rise in the porosity parameter, the transverse deflection of the porous conoidal shell increases for both types of FGM. The dimensionless value of stresses increases along with the thickness ratio of porous conoidal shell. Fig. 6 shows the dimensionless value of axial stress σ¯xfor FGM porous conoidal shell having various types of boundary conditions. It is remarkable to notice that for all types of BC except CFCF, the nature of graph along thickness coordinate is the same except for CFCF type boundary condition, in which the nature of graph is opposite to the other considered boundary conditions. The non-dimensional numerical value of axial stress decreases with an increase in constraints at end support. The dimensionless values of axial stress σ¯yand shear stress τ¯xyare shown in Fig. 7 and Fig. 8 respectively.

Figure 3 Variation of non-dimensional deflection of FGM-1 porous conoidal shell under uniform loading.
Figure 3

Variation of non-dimensional deflection of FGM-1 porous conoidal shell under uniform loading.

Figure 4 Variation of non-dimensional deflection of FGM-1 porous conoidal shell under uniform loading for β = 0.1.
Figure 4

Variation of non-dimensional deflection of FGM-1 porous conoidal shell under uniform loading for β = 0.1.

Figure 5 Variation of non-dimensional deflection of FGM porous conoidal shell under uniform loading.
Figure 5

Variation of non-dimensional deflection of FGM porous conoidal shell under uniform loading.

Figure 6 Effect of porosity on dimensionless axial stress of FGM-1 porous conoidal shell under uniform loading.
Figure 6

Effect of porosity on dimensionless axial stress of FGM-1 porous conoidal shell under uniform loading.

Figure 7 Effect of porosity on dimensionless axial stress of FGM-1 porous conoidal shell under uniform loading.
Figure 7

Effect of porosity on dimensionless axial stress of FGM-1 porous conoidal shell under uniform loading.

Figure 8 Effect of porosity on dimensionless shear stress of FGM-1 porous conoidal shell under uniform loading.
Figure 8

Effect of porosity on dimensionless shear stress of FGM-1 porous conoidal shell under uniform loading.

Table 5

Variation of non-dimensional deflection of simply supported FGM porous conoidal shell subjected to uniform loading.

nβa/h
5102050100
002.99821.29810.47000.11280.0380
(0.45,0.50)(0.40,0.50)(0.325,0.50)(0.25,0.50)(0.20,0.50)
0.13.27701.42180.51480.12350.0416
(0.45,0.50)(0.40,0.50)(0.325,0.50)(0.225,0.50)(0.20,0.50)
0.23.59781.56420.56630.13580.0457
(0.45,0.50)(0.40,0.50)(0.325,0.50)(0.225,0.50)(0.20,0.50)
0.203.28371.42060.52220.12810.0435
(0.45,0.50)(0.40,0.50)(0.325,0.50)(0.25,0.50)(0.20,0.50)
0.13.62381.57050.57780.14190.0481
(0.45,0.50)(0.40,0.50)(0.325,0.50)(0.25,0.50)(0.20,0.50)
0.24.02141.74600.64320.15830.0537
(0.45,0.50)(0.40,0.50)(0.325,0.50)(0.225,0.50)(0.20,0.50)
0.503.52201.53100.56930.14160.0482
(0.45,0.50)(0.40,0.50)(0.325,0.50)(0.25,0.50)(0.20,0.50)
0.13.91691.70640.63580.15860.0539
(0.45,0.50)(0.40,0.50)(0.325,0.50)(0.25,0.50)(0.20,0.50)
0.24.38821.91640.71610.17920.0610
(0.45,0.50)(0.40,0.50)(0.325,0.50)(0.25,0.50)(0.20,0.50)
103.81541.66780.62690.15770.0538
(0.45,0.50)(0.40,0.50)(0.325,0.50)(0.25,0.50)(0.20,0.50)
0.14.28371.87810.70830.17900.0611
(0.45,0.50)(0.40,0.50)(0.325,0.50)(0.25,0.50)(0.20,0.50)
0.24.85762.13680.80960.20570.0702
(0.45,0.50)(0.40,0.50)(0.325,0.50)(0.25,0.50)(0.20,0.50)
1005.10822.25170.84540.21080.0718
(0.45,0.50)(0.40,0.50)(0.325,0.50)(0.25,0.50)(0.20,0.50)
0.15.93532.62800.99130.24840.0847
(0.45,0.50)(0.40,0.50)(0.325,0.50)(0.25,0.50)(0.20,0.50)
0.27.03623.13141.18900.30010.1025
(0.45,0.50)(0.40,0.50)(0.325,0.50)(0.25,0.50)(0.20,0.50)
Metal06.46752.80021.01390.24330.0820
(0.45,0.50)(0.40,0.50)(0.325,0.50)(0.25,0.50)(0.20,0.50)
0.17.82973.38671.22450.29340.0988
(0.45,0.50)(0.40,0.50)(0.325,0.50)(0.25,0.50)(0.20,0.50)
0.29.68184.20941.52380.36540.1231
(0.45,0.50)(0.40,0.50)(0.325,0.50)(0.225,0.50)(0.20,0.50)
Table 6

Variation of non-dimensional deflection of simply supported FGM porous conoidal shell subjected to sin-sin loading.

nβa/h
5102050100
001.99140.88030.33550.09720.0419
(0.475,0.50)(0.425,0.50)(0.325,0.50)(0.25,0.50)(0.20,0.50)
0.12.17420.96350.36720.10640.0457
(0.475,0.50)(0.425,0.50)(0.325,0.50)(0.225,0.50)(0.20,0.50)
0.22.38471.05930.40360.11690.0501
(0.475,0.50)(0.425,0.50)(0.325,0.50)(0.225,0.50)(0.20,0.50)
0.202.19450.97320.37880.11320.0487
(0.475,0.50)(0.425,0.50)(0.325,0.50)(0.25,0.50)(0.20,0.50)
0.12.42031.07600.41930.12560.0539
(0.475,0.50)(0.425,0.50)(0.325,0.50)(0.25,0.50)(0.20,0.50)
0.22.68481.19670.46720.14030.0601
(0.475,0.50)(0.425,0.50)(0.325,0.50)(0.25,0.50)(0.20,0.50)
0.502.36301.05430.41620.12620.0544
(0.475,0.50)(0.425,0.50)(0.325,0.50)(0.25,0.50)(0.20,0.50)
0.12.62741.17610.46550.14170.0609
(0.475,0.50)(0.425,0.50)(0.325,0.50)(0.25,0.50)(0.20,0.50)
0.22.94411.32250.52540.16070.0688
(0.475,0.50)(0.425,0.50)(0.325,0.50)(0.25,0.50)(0.20,0.50)
102.56901.15320.46090.14110.0609
(0.475,0.50)(0.425,0.50)(0.325,0.50)(0.25,0.50)(0.20,0.50)
0.12.88531.30060.52210.16080.0691
(0.475,0.50)(0.425,0.50)(0.325,0.50)(0.25,0.50)(0.20,0.50)
0.23.27491.48310.59880.18570.0796
(0.475,0.50)(0.425,0.50)(0.325,0.50)(0.25,0.50)(0.20,0.50)
1003.43131.54410.61210.18250.0790
(0.475,0.50)(0.45,0.50)(0.325,0.50)(0.25,0.50)(0.20,0.50)
0.13.98871.80380.71870.21510.0930
(0.475,0.50)(0.45,0.50)(0.325,0.50)(0.25,0.50)(0.20,0.50)
0.24.73422.15290.86380.26000.1123
(0.475,0.50)(0.45,0.50)(0.325,0.50)(0.25,0.50)(0.20,0.50)
Metal04.29571.89890.72380.20980.0903
(0.475,0.50)(0.425,0.50)(0.325,0.50)(0.25,0.50)(0.20,0.50)
0.15.19682.29560.87350.25310.1089
(0.475,0.50)(0.425,0.50)(0.325,0.50)(0.225,0.50)(0.20,0.50)
0.26.41712.85061.08620.31470.1349
(0.475,0.50)(0.425,0.50)(0.325,0.50)(0.225,0.50)(0.20,0.50)
Table 7

Variation of non-dimensional deflection of simply supported FGM porous conoidal shell subjected to uniform loading (a/h = 10).

hl/hhβn
Ceramic0.20.5110Metal
0.2501.29811.42061.5311.66782.25172.8002
(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)
0.11.42181.57051.70641.87812.6283.3867
(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)
0.21.56421.7461.91642.13683.13144.2094
(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)
0.2001.37561.50601.62201.76542.37832.9674
(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)
0.11.50661.66481.80771.98772.77443.5895
(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)
0.21.65741.85072.03012.26133.30364.4601
(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)
0.1501.45631.59511.71711.86742.51003.1416
(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)
0.11.59491.76331.91372.10242.92663.8006
(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)
0.21.75431.96022.14892.39153.48274.7209
(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)
0.1001.53981.68771.81601.97352.64613.3215
(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)
0.11.68611.86552.02382.22163.08394.0189
(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)
0.21.85442.07382.27252.52683.66764.9902
(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)
0.0501.62531.78291.91792.08302.78583.5059
(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)
0.11.77941.97072.13732.34463.24514.2425
(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)
0.21.95682.19072.39982.66653.85725.2658
(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)
0.0001.71191.88002.02212.19492.92783.6928
(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)
0.11.87402.07792.25332.47043.40914.4690
(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)
0.22.06052.30982.53002.80944.04995.5448
(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)(0.40,0.50)

5 Conclusion

This paper focuses on the influence of porosity on the bending analysis of FGM porous conoidal shell based on TSDT subjected to the various types of loads using an efficient C0 FE model. The subsequent outcomes of the present study were written below for the volume fraction indices, skew angles, thickness ratios, hl/hh ratios and various types of end support.

  1. The presence of porosity increases the deflection in comparison with perfect conoidal shell.

  2. The dimensionless deflection of the FGM porous conoidal shell increases along with the porosity volume fraction and volume fraction index.

  3. The porosity shows an insignificant effect on the dimensionless stresses of FGM porous conoidal shell except CFCF boundary condition.

  4. The non-dimensional deflection and dimensionless stresses increase with increase in thickness of shell.

  5. The dimensionless deflection and normal stresses at the top of the rhombic shell increase with a reduction in the hl/hh ratio subjected to uniform and sin-sin loading.

References

[1] Koizumi M. Compos. Part B 1997, 28B, 1–4.10.1016/S1359-8368(96)00016-9Search in Google Scholar

[2] Eslami MR, Babaei MH, Poultangari R. Int. J. Press. Vessel. Pip. 2005, 82, 522–527.10.1016/j.ijpvp.2005.01.002Search in Google Scholar

[3] Zhao X, Lee YYÃ, Liew KM. Int. J. Mech. Sci. 2009, 51, 694–707.10.1016/j.ijmecsci.2009.08.001Search in Google Scholar

[4] Kashtalyan M, Menshykova M. Compos. Struct. 2009, 87, 36–43.10.1016/j.compstruct.2007.12.003Search in Google Scholar

[5] Ameur M, Tounsi A, Mechab I, El Bedia AA. KSCE J. Civ. Eng. 2011, 15, 1405-1414.10.1007/s12205-011-1361-zSearch in Google Scholar

[6] Natarajan S, Manickam G. Finite Elem. Anal. Des. 2012, 57, 32–42.10.1016/j.finel.2012.03.006Search in Google Scholar

[7] Kumar A, Chakrabarti A, Ketkar M. Lat. Am. j. solids struct. 2013, 10, 891-919.10.1590/S1679-78252013000600002Search in Google Scholar

[8] Bessaim A, Houari MS, Tounsi A, Mahmoud S, Bedia EAA. J. Sandw. Struct. Mater. 2013, 15, 671–703.10.1177/1099636213498888Search in Google Scholar

[9] Tornabene F, Viola E. Meccanica 2013, 48, 901–930.10.1007/s11012-012-9643-1Search in Google Scholar

[10] Sayyad AS, Ghugal YM. Int. J. Mech. Mater. Des. 2014, 10, 247–267.10.1007/s10999-014-9244-3Search in Google Scholar

[11] Viola E, Rossetti L, Fantuzzi N, Tornabene F. Compos. Struct. 2014, 112, 44–65.10.1016/j.compstruct.2014.01.039Search in Google Scholar

[12] Asemi K, Salehi M, Akhlaghi M. Int. J. Press. Vessel Pip. 2014, 119, 52–61.10.1016/j.ijpvp.2014.03.002Search in Google Scholar

[13] Dai H, Dai T. Meccanica 2014, 49, 1069–1081.10.1007/s11012-013-9853-1Search in Google Scholar

[14] Xiang S, Liu Y-Q. J. Sandw. Struct. Mater. 2016, 18, 579–596.10.1177/1099636216647928Search in Google Scholar

[15] Parihar RS, Setti SG, Sahu RK. Sci. Eng. Compos. Mater. 2016, 25, 309–336.10.1515/secm-2015-0395Search in Google Scholar

[16] Alipour MM, Shariyat M. Int. J. Mech. Mater. Des. 2017, 13, 125–157.10.1007/s10999-015-9321-2Search in Google Scholar

[17] Hadid HA. PhD Dissertation, University of Southampton, 1964.Search in Google Scholar

[18] Das AK, Bandyopadhyay JN. Comput. Struct. 1993, 49, 531–536.10.1016/0045-7949(93)90054-HSearch in Google Scholar

[19] Ghosh B, Bandyopadhyay JN. Comput. Struct. 1994, 53, 9–18.10.1016/0045-7949(94)90124-4Search in Google Scholar

[20] Das HS, Chakravorty D. Adv. Vib. Eng. 2009, 8, 321–328.Search in Google Scholar

[21] Kumari S, Chakravorty D. Int. J. Eng. Sci. Technol. 2010, 2, 54–70.10.4314/ijest.v2i4.59199Search in Google Scholar

[22] Bakshi K, Chakravorty D. Thin-Walled Struct. 2014, 76, 1–7.10.1016/j.tws.2013.10.021Search in Google Scholar

[23] Malekzadeh Fard K, Baghestani AM. Aerosp. Sci. Technol. 2017, 69, 136–48.10.1016/j.ast.2017.06.021Search in Google Scholar

[24] Akbaş ŞD. J. Appl. Comput. Mech. 2017, 3, 199–207.Search in Google Scholar

[25] Al Rjoub YS, Hamad AG. KSCE J. Civ. Eng. 2017, 21, 792–806.10.1007/s12205-016-0149-6Search in Google Scholar

[26] Eltaher MA, Fouda N, El-midany T, Sadoun AM. J. Brazilian Soc. Mech. Sci. Eng. 2018, 40, 141.10.1007/s40430-018-1065-0Search in Google Scholar

[27] Kiran MC, Kattimani SC, Vinyas M. Compos. Struct. 2018, 191, 36–77.10.1016/j.compstruct.2018.02.023Search in Google Scholar

[28] Gupta A, Talha M. Int. J. Struct. Stab. Dy. 2018, 18. doi:10.1142/S021945541850013X.10.1142/S021945541850013XSearch in Google Scholar

[29] Gupta A, Talha M. Int. J. Mech. Mater. Des. 2017, 14. doi:10.1007/s10999-017-9369-2.10.1007/s10999-017-9369-2Search in Google Scholar

[30] Gupta A, Talha M. Arab. J. Sci. Eng. 2018, 43, 4931. doi:10.1007/s13369-018-3240-0.10.1007/s13369-018-3240-0Search in Google Scholar

[31] Ferreira AJM, Roque CMC, Jorge RMN, Fasshauer GE, Batra RC. Mech. Adv. Mater. Struct. 2007, 14, 577–87.10.1080/15376490701672732Search in Google Scholar

Received: 2018-09-25
Accepted: 2019-05-27
Published Online: 2019-10-01
Published in Print: 2019-01-28

© 2019 Md Irfan Ansari et al., published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 Public License.

Articles in the same Issue

  1. Analysis on the impact response of fiber-reinforced composite laminates: an emphasis on the FEM simulation
  2. Artificial neural network for predicting the flexural bond strength of FRP bars in concrete
  3. Cyclic behavior of GFRP strengthened infilled RC frames with low and normal strength concrete
  4. Durability of basalt fiber-reinforced polymer bars in wet-dry cycles alkali-salt corrosion
  5. Effect of B4C particle size on the mechanical properties of B4C reinforced aluminum matrix layered composite
  6. Enhanced dielectric properties of BaTiO3 ceramics with cerium doping, manganese doping and Ce-Mn co-doping
  7. Free and forced vibration analysis of rectangular/circular/annular plates made of carbon fiber-carbon nanotube-polymer hybrid composites
  8. Influence of nano-SiO2 on the bonding strength and wear resistance properties of polyurethane coating
  9. Investigation of wear behavior of nanoalumina and marble dust-reinforced dental composites
  10. Negative effect of clay fillers on the polyvinyl alcohol biodegradation: technical note
  11. Photocatalytic activity of Cu2O/ZnO nanocomposite for the decomposition of methyl orange under visible light irradiation
  12. Sub-surface mechanical properties and sub-surface creep behavior of wood-plastic composites reinforced by organoclay
  13. Surface integrity in wire-EDM tangential turning of in situ hybrid metal matrix composite A359/B4C/Al2O3
  14. The influence of the WC-Co composite microstructure model on stress field heterogeneity at the microstructure level: FEM based study
  15. Vibration-damping characterization of the basalt/epoxy composite laminates containing graphene nanopellets
  16. A review on nanocomposite hydrogels and their biomedical applications
  17. Optimization and simulation analysis of structure parameters of OPCM ultrasonic longitudinal wave actuating element
  18. Research Article
  19. Preparation of POSS-triol/wollastonite composite particles by Liquid phase mechanochemical method and its application in UV curable coatings
  20. Research on preload relaxation for composite pre-tightened tooth connections
  21. Dough moulding compound reinforced silicone rubber insulating composites using polymerized styrene butadiene rubber as a compatibilizer
  22. Hydration And Microstructure Of Astm Type I Cement Paste
  23. Effects of NiO content on the microstructure and mechanical properties of AgSnO2NiO composites
  24. Overall buckling behaviour of laminated CFRP tubes with off-axis ply orientation in axial compression
  25. UV sensing optode for composite materials environment monitoring
  26. On crushing characteristics of hybrid sandwich aluminum-cardboard panels reinforced with glass fiber composite rods
  27. Preparation and characterization of Ni-Cu composite nanoparticles for conductive paints
  28. A research on the preparation of oil-adsorbing hydrophobic porous resins by high internal phase emulsions (HIPEs) template
  29. Material characteristics of random glass-mat-reinforced thermoplastic under cryogenic thermal cycles
  30. Differentiation of non-black fillers in rubber composites using linear discriminant analysis of principal components
  31. Research Article
  32. Efficiency of TiO2 catalyst supported by modified waste fly ash during photodegradation of RR45 dye
  33. Synthesis and performance of polyurethane/silicon oxide nano-composite coatings
  34. Study on preparation of magnesium-rich composite coating and performance enhancement by graft modification of epoxy resin
  35. Research Article
  36. Mechanical and wear properties of polyetheretherketone composites filled with basalt fibres
  37. Mechanical Properties of Al 25 wt.% Cu Functionally Graded Material
  38. Research Article
  39. Weight reduction of a carbon fibre composite wheel
  40. Synthesis, electrical properties, and kinetic thermal analysis of polyaniline/ polyvinyl alcohol - magnetite nanocomposites film
  41. Seismic Behaviour of TRC-Strengthened RC Columns under Different Constraint Conditions
  42. Characterization of neat and modified asphalt binders and mixtures in relation to permanent deformation
  43. Microstructures, interface structure and room temperature tensile properties of magnesium materials reinforced by high content submicron SiCp
  44. Research Article
  45. Effect of Cutting Temperature on Bending Properties of Carbon Fibre Reinforced Plastics
  46. Mechanical and tribological properties of B-C-N coatings sliding against different wood balls
  47. Thermal conductivity of unidirectional composites consisting of randomly dispersed glass fibers and temperature-dependent polyethylene matrix
  48. Effects of Waste Eggshells addition on Microstructures, Mechanical and Tribological Properties of Green Metal Matrix Composite
  49. Investigation of porosity effect on flexural analysis of doubly curved FGM conoids
  50. Review Article
  51. Utilization of tailings in cement and concrete: A review
  52. Research Article
  53. Equivalent stiffness prediction and global buckling analysis using refined analytical model of composite laminated box beam
  54. Mechanochemical synthesis of zincite doped with cadmium in various amounts
  55. Size-dependent vibration analysis of graphene-PMMA lamina based on non-classical continuum theory
  56. Automated, Quality Assured and High Volume Oriented Production of Fiber Metal Laminates (FML) for the Next Generation of Passenger Aircraft Fuselage Shells
  57. Research Article
  58. An investigation of the stitching effect on single lap shear joints in laminated composites
  59. The low-velocity impact and compression after impact (CAI) behavior of foam core sandwich panels with shape memory alloy hybrid face-sheets
  60. Effect of granulometric distribution on electromagnetic shielding effectiveness for polymeric composite based on natural graphite
  61. The enhancement of filament winding in marine launching rubber gasbag
  62. Research on ELID Grinding Mechanism and Process Parameter Optimization of Aluminum-Based Diamond Composites for Electronic Packaging
Downloaded on 2.11.2025 from https://www.degruyterbrill.com/document/doi/10.1515/secm-2019-0026/html?lang=en
Scroll to top button