Static Analysis of Laterized Concrete Cylindrical Shells under Food Grain Pressures by Initial Value Method

DOI : 10.17577/IJERTV2IS60026

Download Full-Text PDF Cite this Publication

Text Only Version

Static Analysis of Laterized Concrete Cylindrical Shells under Food Grain Pressures by Initial Value Method

T. U. Nwakonobi 1; N. N. Osadebe2 and C. C. Mbajiorgu3

1PhD, Senior Lecturer, Department of Agricultural and Environmental Engineering, University of Agriculture, P.M.B.2373, Makurdi, Benue State, Nigeria.

2Professor of Civil Engineering, Department of Civil Engineering, Faculty of Engineering, University of Nigeria, Nsukka, Enugu State, 410001, Nigeria.

3 Professor of Agricultural and Bioresources Engineering, Department of Agricultural and Bioresources Engineering, Faculty of Engineering,University of Nigeria, Nsukka, Enugu State, 410001, Nigeria.

Abstract: In this study, the performance of a laterized concrete cylindrical shell for storage of food grains was investigated. The structural characteristics values of laterized concrete at optimum mix proportion were applied in performing the static analysis of the cylindrical shell structure under the action of food grain pressures using Pasternacks equations formulated on the basis of theory of shell of revolution. The reduced equation of static equilibrium was solved by the initial- value method. A generalized solution of cylindrical shell made of laterized concrete was obtained and used to evaluate the performance of the structure through the determination of deflection and stresses. The results of the analysis indicate that the maximum stress developed due to soybean grain loading is greater than the material strength in silo thickness of 100mm, 20m in diameter and 8m in height. The wall thickness of 150mm was found adequate for silo of the same size.

Keywords: static analysis; Food grain; storage; laterized concrete; cylindrical shell; Initial-value method

  1. Introduction

    The post harvest handling of food grains and seeds has constituted a matter of great concern to agricultural engineers and other agricultural experts. In developed countries, efficient storage and processing of food have been recognized as major factors in the solution of food problems (Aboaba, 1996).

    In Nigeria, about 70% of farming population stores their farm produce using traditional storage structures as the use of modern structures has become unaffordable. The use of these structures are associated with numerous problems which, according to several writers (Olumeko 1996; Adejum and Raji 2007; Nwakonobi and Mbajiougu, 2009); they are not durable, not watertight, not rodent proof, susceptible to insert attack, low viability of grains, change in grain texture, appearance and flavour and of low capacity. These defects and shortcoming of the traditional structures usually result in tremendous losses estimated to be up to 30 40% (Igbeka and Olumeko, 1991). To overcome these problems efforts are now being directed towards evolving low-cost, effective storage structures made of locally available materials for use by low to medium scale farmers. Laterized concrete has been identified by (Lasisi and Osunade, 1990) as a possible low-cost construction material.

    Laterized concrete is concrete in which some or the entire fine aggregate is from lateritic soil. Sand (fine aggregate) which is the second largest of the concrete volume is not always available. According to Osunade (1997), sand is the second most costly item per unit volume of concrete produced. Thus, replacing sand with lateritic soil can reduce construction cost since the later is abundant and readily available at construction sites. Several investigators have studied the applications of lateritic soils as a component in concrete (Osunade et al 1990; Osunade, 1993; [Osunade, 1994 and Osadebe and Nwakonobi, 2007, Nwakonobi, et al, 2012) and some of the research findings were very encouraging. However a lot still remain to be done in order to establish the behaviour of laterized concrete cylindrical shells under different types of loadings.

    The primary objective of this study was to establish the behaviour of laterized concrete cylindrical shells under food grain loadings. The specific objectives are to:

    • Analyze the cylindrical shell structure made of laterized concrete using pasternacks semi- moment equation which was formulated on the basis of cylindrical theory of shell of revolution.

    • Develop a general solution for the design of a shell cylindrical silo built of laterized concrete.

    • Evaluate the performance of the structure through the determination of deflection and stresses

1.1 Background

A known problem in the study of shells is how to determine whether analysis based on a particular theorem is adequate, satisfactory or successful. Because the analysis of a shell structure is very tedious and cumbersome, theories have been developed to simplify the design of a shell structures. Some of the shell theories includes; Membrane theory, variational theory, Bending or moment theory and Half (or semi) moment theory.

According to Zingoni et al, (2000) membrane theory is applicable to shell with finite bending rigidity but in which the moments that are developed are so small as to be negligible. The theory assumes that silo walls are subjected to tensile forces only. However, there are some limitations associated with the application of membrane theory in the analysis of the shell structure. This theory according to Timoshenko and Woinowsky-Krieger (1959) fails to represent the true stress in those portions of the shell close to the edges, since the edge conditions usually cannot be completely satisfied by considering only membrane stresses. Zingoni et al, (2000) also noted similar conditions under which membrane theory is not valid which include conditions at certain locations in the interior of the shell. There exist at these locations transverse shearing force and moments to satisfy the condition of continuity of the shell. Any component of the forces reaction at the support in the direction normal to the shells middle surface would induce transverse shear in the shell, which the, membrane theory does not admit.

The variational theory is deficient as bending term are not adequately accounted for due to approximation in the strain energy function Muhammad (2000). For the bending or moment theory, it can be used in obtaining reduced (finite) set of equations which involve no more than stress resultants and stress couples for the description of stress distribution, together with the approximate boundary conditions. Bending theory is the most accurate of all the theories but the analysis is very cumbersome as the analyst are involved in the generally very difficult task of searching for the actual particular solution corresponding to a given surface loading. But in this theory, transverse shearing force and bending terms are accounted for in the stress analysis of shell structure. However, half or semi moment theory which is

also known as pasternacks theory tried to simplify the moment theory and this theory was adopted in this study based on its strength.

Several researchers using various methods solve reduced equation of static equilibrium of circular cylindrical shell structure and Ren and Kuan-Chan (2001) tried to develop series solutions to the differential equations for a complete cylindrical shell under various point loads. Instead of expanding the solution series in both axial and circumferential coordinates and each terms in the series represent an exact solution to the set of basic equations. An exact three- dimensional elasticity solution is obtained from an infinite long, thick transversely isotropic circular cylindricalshell panel, simply supported along the longitudinal edges and subjected to a radial patch load (Chandrashekhera and Rao, 1996). The boundary value problem is reduced to Bessels differential equation using a set of three displacement functions. In another study, Chandrashekhera and Rao,(1996) presented an approximate three dimensional elastic solution for an infinitely cross-ply laminated circular cylindrical shell panel with simply supported boundary conditions, subjected to an arbitrary discontinuous transverse loading. The solution is based on the principal assumption that the ratio of the thickness of the lamina to its middle surface radius is negligible compared to unity.

Mohammad (2000) carried out analysis of the cylindrical shell subject to local and continuous axi- symmetric loads using pasternacks equations. The reduced differential equation was solved using classical and initial value methods. The results obtained in each of the methods are the same indicating the applicability of either of the two methods in solving equation of static equilibrium of cylindrical shell. Initial value method has been applied in the analysis of circular cylindrical shell under hydrostatic pressure and ring force (Amodou, 2002). Reinforced concrete dome structures designed with reinforcement ratios were analyzed to check the applicability of the finite element analysis techniques (Song et al, 2002). Their results showed that the technique could be applied effectively to failure analysis of various types of reinforced concrete shell structures.




+ +

Fig. 1: Direct and shear forces in an element of a shell




Fig. 2: Bending and twisting moments in an element of a shell

  1. Stress Equilibrium Equations

    The, and are taken at a given point on the middle surface to be the direction to the axis of

    the cylinder, the tangent to the circumference, and the normal to the middle surface of the shell respectively. Assuming a transverse distributed load whose components in , and directions are ,

    and respectively, R is the radius of the cylinder, and are the sizes of the element in and

    directions. The direct and shear forces as well as bending and twisting moments in an element of a shell are shown in Figs.1 and 2, respectively.

    The equilibrium of forces , and -direction are given as follows:

    + + X = 0 (1)

    + + + y = 0 (2)

    Z + + – = 0 (3)

    The moment equilibrium with respect to x and y- axes in:

    + – = 0 (4)

    + – = 0 (5)

    In semi- moment (or pasternacks) theory as applied to cylindrical shell, twisting moments are usually neglected. This result in:

    = = 0 (6)

    Based on the law of reciprocity of shearing stresses, we have:

    = = S (7)

    Substitution of equation (6) into equations (4) and (5) gives:

    = (8)

    = (9)

    Substitution of equations (7), (8) and (9) into equations (1), (2) and (3) to obtain:

    + + = 0 (10)

    + + + = 0 (11)

    + + + = 0 (12)

    Equations (10) through (12) represent the equation of the static equilibrium of a cylindrical shell under

    distributed forces of , Y and Z in and directions only. The systems are statically indeterminate since there are five unknowns for only three equations. In order to solve the systems, the

    following stress strain relations were needed to facilitate the determination of the unknowns namely and :

    = ( + ) (13)

    = ( + ) (14)

    = (15)

    = (16)

    = (17)


    = (18)

    = + (19)

    = (20)

    = ( (21)

    = (22)

    = (23)

    Substituting equations (18) through (23) into equations (13) through (17) gives the stress resultants and stress couples in terms of the displacements and.

    = [ + + )] (24)

    = [ + +] (25)

    = [ + ] (26)

    = [ – ( – )] (27)

    = [ ( ) ( )] (28)

    Substitution of above expressions into equations (10), (11) and (12) gives:


    ) + . )- – ] + . (30)



    From the Pasternaks theory on cylindrical shell, the following assumptions were made:

    1. The middle plane of the shell is inextensible in y direction such .

    2. It is assumed that the stress generated by imposed load from stored grain do not depend on the y

      but only on the length

    3. The normal force acting on the transverse section of the shell is neglected.

      Consequently, all derivatives with the respect to y are set to be zero.

      In case of grain pressure, the external forces acting on the cylindrical shell (see Fig.3) are such that y and where stands for the weight per unit volume of stored grain material (soybean).

      W(x) L


      D= 2R

      Fig.3: A section of circular cylindrical shell

      Applying the above assumptions to equations (24) through (28) to get.




      = M = – (35)

      = – (36)


        1. Differential Equation of Equilibrium

          Making use of the previous results, it can be seen that equations (29) and (30) are satisfied. Substituting these same results into equation (31) to get:

          + 4

          Which can be written as + = (37)




          Equation (37) is due to pasternack and is differential equation for cylindrical shell.

        2. General Solution of Static Equilibrium by Initial Value Method

          In case of food grain pressure, the distribution of forces is triangular (see Fig. 4)

          The static differential equation is given by:

          + = (38)

          Where; and are as defined in section 3.1 and




          L R= radius

          grain pressure

          h= thickness


          Fg.4: Grain pressure distribution.

          Where; is the lateral- to- vertical pressure ratio given by Rankine as: = [- ], angle of internal friction of grain

          = unit weight of grain Equation (38) can be written as;

          + = (39)

          From the initial value method; Let Z =

          Differentiating above with respect to gives:

          = (40)


          = ( = (41)

          Substitution of equations (40) and (41) into equation (39) gives:

          = + 4 =

          Dividing all through by gives:

          = + 4 = (42)

          Substitute for and as defined in equation (37), to get:

          + 4 = z (43)

        3. The Initial Value Homogenous Solution

      The Homogenous equation is given by:

      + 4 = 0 (44)

      This solution is found by using the result of the classical homogenous solution obtained as: = + + + + + + + (44)

      By taking and ; equation (44) is put as:

      = + + + (45)

      Where: are arbitrary constants +

      cosz +

      Equation (45) can now be put as:


      Where: = (47)

      The constants , and ( are 16 in number and are obtained by differentiating

      equation (47) with respect to z till the fourth order and making use of the supplementary initial conditions which are given as:

      (0) = 0 (0) = 0 (0) = 0

      (0) = 1 (0) = 0 (0) = 0

      (0) = 0 (0) = 1 (0) = 0

      (0) = 0 (0) = 0 (0) = 1

      The constants are obtained as:

      = 1 = 0 = 0 = 0

      = 0 = = 0 =

      = 0 = = 0 = –

      = 0 = 0 = = 0

      These constants are substituted into equation (47) to get:



      (sin h z sin z) (50)


      , are four independent expressions. Therefore, the homogenous solution is written as:

      (z) + (z) + (z) .. (52)

      The arbitrary constants and are determined using the initial conditions.

      The relationships between the deflection and the rotation moment , shear force Q(z) and hoop tension N(z) are derived as follows:

      Recalling equation (46)

      Using it together in the following expression = to obtain:


      = (54)

      = (55)

      = (56)

      Taking and as initial value o and respectively, the

      arbitrary constants are are found by making use of the initial supplementary conditions given above. The constants are obtained as follows:


      Substituting for and into equations (46) and (53) through (56) to obtain





      3.4 The Initial Value Particular Solution






      Q A Q A




      dQ B

      dx dQ B



      Fig. 5: Elemental Force Acting on side of Grain Silo


      Considering the above Figure 5, the particular solution is found by assuming an origin transformation from A to B;



      The distributed load at new origin B is given by


      Therefore the elemental force is expressed as:

      = (66)

      To get the particular solution, use is made of the homogeneous solution in which only parameters similar

      to are considered. Therefore, the coefficients of in the homogeneous solution are used with the appropriate sign to formulate the particular solutions:

      Integrating the above expressions gives the actual particular solutions;

      [Z – (z)] (72)






      [1 – (z)] (73)

      (z) (74)

      (z) (75)

      [z – (z)] (76)

      3.5 The Initial Value General Solution

      This is obtained by adding homogeneous and particular solutions:


      (z)(z) (77)

      (z) (z) (z)

      (z)(z) (78)


      (z)(z) (79)


      (z)(z)(z) (80)


      (z)] (81)

      Boundary conditions are applied in equations (76) through (80) At = = 0

      At (82)

      Substitution of and equation (81) into equations (76) and (77) which gives:



      Solving equations (82) and (83), to obtain

      = (85)

      = (85)

      Substituting , , , and in equations (76) through (80) to get:

      kR2 Y (L) Y (L)Y (L) Y (L)Y (L) LY (L)







      W (z)

      1 4 3

      Y2 (z)

      Eh Y2 (

      L)Y3 (

      L) Y1 (

      L)Y4 ( L)

      kR2 Y2 (L) LY1(L)






      EhY (L)Y (L) Y (L)Y

      (L) 4

      • Eh

      z Y2 (z)

      2 3 1 4

      kR2 Y4 (L) Y2 (L)Y3 (L) Y1(L)Y4 (L) LY3 (L)

      (z) EhY (L)Y (L) Y (L)Y (L)


      2 3 1 4

      kR2 Y2 (L) LY1(L)



      EhY (L)Y (L) Y (L)Y (L)Y3 (z)

      Eh 1 Y1(z)

      2 3 1 4

      M (z) k Y4 (L) Y2 (L)Y3 (L) Y1 (L)Y4 (L) LY3 (L)Y (z)

      3 Y (L)Y (L) Y (L)Y (L) 4

      2 3 1 4

      kY2 (L) LY1(L) Y


      • k Y



      43Y (L)Y (L) Y (L)Y (L) 2

      3 4

      2 3 1 4

      Q(z) kY4 (L) Y2 (L)Y3 (L) Y1(L)Y4 (L) LY3 (L)Y (z)

      2 Y (L)Y (L) Y (L)Y (L) 3

      2 3 1 4

      kY2 (L) LY1(L)

      Y (z)

      k Y (z)


      42Y (L)Y (L) Y (L)Y

      (L) 1

      2 3

      2 3 1 4

      N (z) kRY4 (L) Y2 (L)Y3 (L) Y1(L)Y4 (L) LY3 (L)Y


      Y (L)Y (L) Y (L)Y (L) 2

      2 3 1 4

      kRY2 (L) LY1(L) Y


      • kR z Y



      Y (L)Y (L) Y (L)Y

      (L) 4 2

      2 3 1 4

      The initial value general solutions given in equations (86) through (90) are used to evaluate the

      performance of the cylindrical shell under food grain pressures.

  2. Materials and Method

    The values of laterized concrete properties at optimum mix proportion as determined by Osadebe and Nwakonobi (2007) were adopted in performing the static analysis of the circular cylindrical shell under the action of food grain pressure using pasternack,s equations formulated on the basis of cylindrical theory of shell of revolution

      1. Characteristics of the Silo Shell

        The silo size 20m in diameter (D), 8m in height (L) and wall thickness (H) of 0.1m and 0.150m were selected. The value of laterized concrete properties such as young modules, EL and poissons ratio, are

        18,888.9kN/ and 0.26 respectively (Osadebe and Nwakonobi, 2007). The characteristics of the soybean grain adopted from Nwakonobi and Idike (2002) are given as follows:

        Unit weight of soybean, = .

        Angle of internal friction for soybean grain, = .

      2. Static Analysis of Laterized Cylindrical Shell The general solutions of the initial value method given in equations (86) through (90) were used to evaluate the performance of a cylindrical shell under soybean grain pressures. The values for different parameters

        given in section 4.1 were substituted into equations (86) through (90) and the values of deflection bending , slope shear force and hoop tension computed for each of the

        ten segment of the silo height for the silo wall thicknesses of 100mm and 150mm as presented in Tables 1 and 2, respectively.

        Table 1: Results of analysis by Initial Value method for a shell thickness of 100mm















































































        Table 2: Results of analysis by Initial Value method for a shell thickness of 150mm









































































        The maximum stress developed due to the bending moment and hoop tension for both wall thicknesses are as presented in Table 3.

        Table 3: Maximum stress developed by two different wall thicknesses

        Wall thickness mm

        Bending stress N/mm2

        Hoop tension stress N/mm2







          1. Results and Discussion

            Tables 1 and 2 show the results of stresses and deflection analysis for a grain silo shell of thickness

            100 and 150, respectively.

            • Direct deflection

    From Table 2, it is observed that the deflection is zero at both end of silo. The maximum value of deflection is and occur at a height of 5.6m from the top of silo which is high and

    150 thick.

    • Rotation

    From Table 2 the rotation is 0.0001425 radians at the top and is zero at the bottom end. The maximum rotation is negative of 0.0006294 radians and occurs at a height 7.2 from the top.

        • Bending moment

          The bending moment is zero at the top of the grain silo (see Table 2). The maximum sagging moment is 2.9506 kN-m and occurs at height 6.4 from the top while the maximum hogging

          moment is 12.8405 kN-m and occurs at the bottom.

        • Shear force

          Table 2 showed the shear force to be 0.007588 kN at the top and is -28.93818 kN at the bottom end of the silo. The maximum positive shear is 1.96735 kN and occurs at height 5.6 from the

          top of the silo.

        • Hoop tension

          The maximum hoop tension is and occurs at height 5.6 from the top of the silo (see Table 2).

        • Silo Shell thickness

    It is seen from Table 3 that the maximum stress developed due to soybean grain loading is 5.26N/mm2 in wall thickness of , in diameter and in height. While the maximum

    stress developed within wall thickness of was 3.42N/mm2. These stresses developed were compared with the allowable stress of the laterized concrete material which is 4.12 (Nwakonobi and Osadebe, 2007). The wall thickness of 150 was within safe limit while that of 100mm was not.

  3. Conclusion

From the cylindrical shell analysis using the initial-value method, a wall thickness of was found inadequate to retain soybean grain to the silo height of and in diameter without reinforcement. A wall thickness of was found adequate. It is however recommended that minimum reinforcement may be provided in a laterized concrete silo to prevent crack due to thermal effect or temperature

variations. Cement plaster should be applied on both inner and outer surface for more environmental protection and longevity. But in order to validate this analytical study it is recommended that the laterized concrete shell be constructed to measure the deflection that will develop during the storage of food grains. This study has however demonstrated that the theory of shell structure can therefore be applied in the analysis and design of structure built of laterized concrete.


Aboaba, F.O., (1996). Storage and processing of major crops in Nigeria, AFP symposium proceeding, 131-135.

Amodou, A., (2002). Application of the initial value method in the analysis of circular cylindrical shell under hydrostatic pressure and ring forces. M.Eng thesis submitted to the Department of civil engineering, university of Nigeria, Nsukka.

Adejumo, B.A.; Raji, A.O., (2007). Technical appraisal of storage systems in the nigerian sudan savanna.Agricultural Engineering International: the CIGR Ejournal, IX (11).

Chandrashekhera, K. Rao, K.S.N., (1996). Analysis of a long thick orthotropic circular cylindrical shell panel. J. Eng. Mech., 122(6) 575 579.

Igbeka, J.C.; Olumeko, D.O.,(1991). An appraisal of village level Grain Practices in Nigeria. Agricultural Mechanization in Asia, Africa and Latin America. 27( 2) 29-33.

Lasisi, F.; Osunade, J.A., (1990). Laterized concrete masonry as an alternative in building and rural infrastructures. in proc. of international conference on Structural Engineering, seam 2, Kumasi, Ghana.

Muhammad, A.K., (2000). Analysis of thin cylindrical shell subjected to local and continuous axi- symmetric loads. unpublished project report submitted to the Department of Civil Engineering, KNUST Kumasi-Ghana.

Nwakonobi, T.U; Idike, F.I., (2002). Physical properties of soybean. Nigeria J. Tech, 21(1) 90 91.

Nwakonobi, T.U.; Mbajiorgu, C.C., (2009). Review of agricultural storage structures for rice in Nigeria. Paper presented at south-east regional conference of the Nigerian society of Agricultural Engineers held at Abakiliki, Ebonyi State, 16 18 of September.

Nwakonobi, T.U.; Osadebe, N.N.; Mbajiorgu, C.C., (2012). Application of initial value method in the static analysis of laterised concrete cylindrical shells under hydrostatic loading. Int. J. of Structural Engineering, 3(4) 229-249).

Olumeko, D.O., (1996). Assessment of On-farm storage structures in Nigeria. In a Paper presented at NCAM national training workshops, Ilorin, Kwara State, Nigeria.Osunade, J.A.,(1997). Laterized concrete masonry as an alternative in building and rural infrastructures. International Conference Proc. on Structural Engineering, 151 161.

Osunade, J.A.,(1990). Adeyefa PO, F. Lasisi F. Shear and Tensile Properties of Unreinforced Laterized Concrete. Ife J. of Tech, 2 (1) 8-17.

Osunade, J.A., ( 1993). The compressive strength of laterized concrete. the effect of types and sizes of coarse aggregates. Masonry International J. Bulletin, ( 1) 104.

Osunade, J.A., (1994). Effect of grain size ranges of laterite fine aggregate on the shear and tensile strengths of laterized concrete. International J. for Housing Sci. and its Application, USA, 18 (2) 91 104.

Osadebe, N.N.; Nwakonobi, T.U., (2007). Structural characteristics of laterized concrete at optimum mix proportion. Nigeria J. of Tech, 26(1)12 17.

Ren, D.; Kuan-Chen, F., (2001). Solution of complete circular cylindrical shell under concentrated loads. J. Eng. Mech., 127(3) 248 253.

Song, H.; Shim.; S.; Byun, K., (2002). Failure analysis of reinforced concrete shell structures using layered shell element with pressure mode. J. Struct. Eng. 2002; 128(5) 655 664.

Timoshenko, S.; Woinosky-krieyer, S., (1959). Theory of plate and shell. McGraw Hill Book. Inc. New York.

Zingoni, A.J.; Nwakal, I.A; Salahuddin, A., (2000). Theory and analysis of structures: trusses, beams, frames, plates, and shells. UNESCO Publ. Paris, France, Chap 5, 165 214.


b = Breadth of section d = Depth of section

D = Cylindrical rigidity

Ec = Static modulus of elasticity of laterized concrete Fc = Compressive strength of Laterized- concrete

Ftb, = Tensile strenth of laterized concrete h = Thickness of shell

k = Lateral to vertical pressure ratio L = Height of cylindrical silo

M = Bending moment N = Hoop tension

Pw = Hydrostatic pressure Pg = Grain pressure

Q = Shear force

R = Radius of cylindrical shell W =,Deflection

= A constant

= Unit weight of water

g = Unit weight of grain

= Rotation

= Angle of internal friction of grain

= Stress

= Poisson's ratio

= displacement component along the – axis

displacement component along the – axis

displacement component along the axi

thickness of the shell

= is the flexural rigidity of the shel Poissons ratio

Young modulus of the shell material

Leave a Reply