Document Type : Regular Article
Authors
1 Maharashtra Institute of technology, Aurangabad (Maharashtra), India
2 Shreeyash College of Engineering, Aurangabad, India
Abstract
Highlights
Keywords
Main Subjects
Euler-Bernoulli hypothesis disregards the effects of the shear deformation and stress concentration which is in elementary theory of beam bending hence it is suitable for thin beams and is not suitable for deep beams since it is based on the assumption that the transverse normal to neutral axis remains so during bending and after bending, implying that the transverse shear strain is zero. Since theory neglects the transverse shear deformation. It underestimates deflections in case of thick beams where shear deformation effects are significant.
Timoshenko [1] showed that the effect of transverse vibration of prismatic bars. This theory is now widely referred to as Timoshenko beam theory or first order shear deformation theory (FSDT) in the literature. But in this theory transverse shear strain distribution is assumed to be constant through the thickness of beam and thus requires shear correction factor to appropriately represent the strain energy of deformation.
Cowper [2] has given refined expression for the shear correction factor for different cross-sections of beam. The accuracy of Timoshenko beam theory for transverse vibrations of simply supported beam in respect of the fundamental frequency is verified by Cowper [3] with a plane stress exact elasticity solution.
To remove the discrepancies in classical and first order shear deformation theories, higher order or refined shear deformation theories were developed and available in the open literature for static and vibration analysis of beam. Krishna Murthy [4], Baluch et al. [5], Bhimaraddi and Chandrashekhara [6] were presented parabolic shear deformation theories assuming a higher variation of axial displacement in terms of thickness coordinate. These theories satisfy shear stress free boundary conditions on top and bottom surfaces of beam and thus obviate the need of shear correction factor.
Kant and Gupta [7], and Heyliger and Reddy [8] presented finite element models based on higher order shear deformation uniform rectangular beams. However, these displacement based finite element models are not free from phenomenon of shear locking [9, 10].
Dahake and Ghugal [11] studied flexural analysis of thick simply supported beam using trigonometric shear deformation theory. Ghugal and Dahake [12, 13] given the flexural solution for the beam subjected to parabolic loading. Sawant and Dahake [14] developed the new hyperbolic shear deformation theory. Chavan and Dahake [15, 16] presented clamped-clamped beam using hyperbolic shear deformation theory. The displacement and stresses for thick beam given by Nimbalkar and Dahake [17].
Jadhav and Dahake [18] presented bending analysis of deep cantilever beam using steel as material. Manal et al [19] investigated the deep fixed beams using new displacement field. Patil and Dahake [20] carried out finite element analysis using 2D plane stress elements for thick beam. Dahake et al [21] studied flexural analysis of thick fixed beam subjected to cosine load. Tupe et al [22] compared various displacement fields for static analysis of thick isotropic beams.
In literature, most of the researchers have used steel as a beam material. As many parts of the spacecrafts, airplane structures are made up of aluminum due to its low weight density. In this research, an attempt has been made to analyze the aluminum deep cantilever beam subjected to cosine load.
The beam under consideration occupies in Cartesian coordinate system the region:
where x, y, z are Cartesian coordinates, L and b are the length and width of beam in the x and y directions respectively, and h is the thickness of the beam in the z-direction. The beam is made up of homogeneous, linearly elastic isotropic material.
The displacement field of the present beam theory is of the form as given below:
$\begin{array}{l} {u(x,z)=-z\frac{dw}{dx} +\frac{h}{\pi } \sin \frac{\pi z}{h} \phi (x)} \\ {w(x,z)=w(x)} \end{array}$ (1)
where is the axial displacement in x direction and w is the transverse displacement in z direction of the beam. The sinusoidal function is assigned according to the shear stress distribution through the thickness of the beam. The represents rotation of the beam at neutral axis, which is an unknown function to be determined.
Normal strain
$\varepsilon _{x} {\rm \; }{\rm \; \; }{\rm =}\frac{\partial u}{\partial x} =-z\frac{d^{2} w}{dx^{2} } +\frac{h}{\pi } \sin \frac{\pi z}{h} \frac{d\phi }{dx}$ (2)
Shear strain
(3)
Stress-Strain Relationships
(4)
Using the expressions for strains and stresses (2) through (4) and using the principle of virtual work, variationally consistent governing differential equations and boundary conditions for the beam under consideration can be obtained. The principle of virtual work when applied to the beam leads to:
(5)
where the symbol denotes the variational operator. Employing Green’s theorem in Eqn. (4) successively, we obtain the coupled Euler-Lagrange equations which are the governing differential equations and associated boundary conditions of the beam. The governing differential equations obtained are as follows:
(6)
(7)
The associated consistent natural boundary conditions obtained are of following form:
At the ends x = 0 and x = L
(8)
Thus the boundary value problem of the beam bending is given by the above variationally consistent governing differential equations and boundary conditions.
The general solution for transverse displacement w(x) and warping function(x) is obtained using Eqns. (6) and (7) using method of solution of linear differential equations with constant coefficients. Integrating and rearranging the first governing Eqn. (6), we obtain the following equation
(9)
where Q(x) is the generalized shear force for beam and it is given by .
Now second governing Eqn. (7) is rearranged in the following form:
(10)
A single equation in terms of is now obtained using Eqns (11) and (12) as:
(11)
where constants , and in Eqns. (10) and (11) are as follows
The general solution of Eqn. (11) is as follows:
(12)
The equation of transverse displacement w(x) is obtained by substituting the expression of (x) in Eqn. (12) and then integrating it thrice with respect to x. The general solution for w(x)is obtained as follows:
(13)
where, are arbitrary constants and can be obtained by imposing natural (forced) and / or geometric or kinematical boundary / end conditions of beam.
In order to prove the efficacy of the present theory, a numerical example is considered. For the static flexural analysis, a uniform beam of rectangular cross section, having span length ‘L’, width ‘b’ and thickness ‘h’ of homogeneous, elastic andisotropic material is considered. The following material properties for beam are used.
Table 1. Properties of Aluminum 6061-T6, 6061-T651 [13]
Physical Properties |
Quantity |
Density |
2700 kg/m3 |
Ultimate Tensile Strength |
310 MPa |
Modulus of Elasticity |
68.9 GPa |
Notched Tensile Strength |
|
Ultimate Bearing Strength |
|
Bearing Yield Strength |
|
Poisson's Ratio |
0.33 |
Fatigue Strength |
|
Shear Modulus |
|
Shear Strength |
The beam has its origin at left hand side fixed support at x = 0 and free at x = L. The beam is subjected to cosine load, on surface z = +h/2 acting in the downward z direction with maximum intensity of load.
L |
x, u |
q0 |
z, w |
Fig. 2. Cantilever beam with cosine load
Boundary conditions associated with this problem are as follows:
At Free end: at x = L and
At Fixed end: = 0 at x = 0
General expressions obtained for and are as follows:
(14)
(15)
The axial displacement and stresses obtained based on above solutions are as follows
(16)
(17)
(18)
(19)
In this paper, the results for inplane displacement, transverse displacement, inplane and transverse stresses are presented in the following non dimensional form for the purpose of presenting the results in this work.
For beam subjected to cosine load
The transverse shear stresses ( ) are obtained directly by constitutive relation and, alternatively, by integration of equilibrium equation of two dimensional elasticity and are denoted by ( ) and ( ) respectively. The transverse shear stress satisfies the stress free boundary conditions on the top and bottom surfaces of the beam when these stresses are obtained by both the above mentioned approaches.
Table 2: Non-Dimensional Axial Displacement ( ) at(x = L, z = h/2), Transverse Deflection ( ) at (x = L, z =0.0) Axial Stress ( ) at (x = 0, z = h/2)Maximum Transverse Shear Stresses (x=0.01L, z =0.0) and (x, z = 0.0) of the Cantilever Beam Subjected to Cosine Load for Aspect Ratio 4 and 10.
Source |
Aspect ratio |
Model |
|||||
Present |
4 |
TSDT |
-67.5989 |
6.1819 |
36.7529 |
1.8181 |
-2.7877 |
Sawant and Dahake [14] |
HPSDT |
-70.024 |
6.1928 |
39.8104 |
2.1609 |
-4.5581 |
|
Krishna Murty [4] |
HSDT |
-71.2291 |
6.1860 |
37.2887 |
1.9004 |
-2.8916 |
|
Timoshenko [1] |
FSDT |
23.1543 |
6.5444 |
22.2081 |
0.3076 |
3.7597 |
|
Bernoulli-Euler |
ETB |
23.1543 |
5.7541 |
22.2081 |
— |
3.7597 |
|
Present |
10 |
TSDT |
-1055.4548 |
5.8244 |
176.7877 |
7.7501 |
3.2052 |
Sawant and Dahake [14] |
HPSDT |
-1061.5175 |
5.8256 |
178.4137 |
8.3023 |
3.8042 |
|
Krishna Murty [4] |
HSDT |
-1064.5122 |
5.8255 |
172.1017 |
7.8208 |
3.7523 |
|
Timoshenko [1] |
FSDT |
361.7856 |
5.8805 |
138.8010 |
4.8073 |
9.3993 |
|
Bernoulli-Euler |
ETB |
361.7856 |
5.7541 |
138.8010 |
— |
9.3993 |
Fig. 2. Variation of axial displacement ( ) through the thickness of cantilever beam at (x = L, z) for aspect ratio 4.
Fig. 3. Variation of axial displacement ( ) through the thickness of cantilever beam at (x = L, z) for aspect ratio 10.
Fig. 4. Variation of maximum transverse displacement ( ) of beam at (x=L, z = 0) with aspect ratio S.
Fig. 5. Variation of axial stress ( ) through the thickness of beam at (x= 0, z) for aspect ratio 4.
Fig. 6. Variation of axial stress ( ) through the thickness of beam at ( x = 0, z) for aspect ratio 10.
Fig. 7. Variation of transverse shear stress ( ) through the thickness of beam at (x = 0.01L, z) obtain using CR for aspect ratio 4.
Fig. 8. Variation of transverse shear stress ( ) through the thickness of beam at (x = 0.01L, z) obtain using CR for aspect ratio 10.
Fig. 9. Variation of transverse shear stress ( ) through the thickness of beam at (x = 0.01L, z) obtain using EE for aspect ratio 4.
Fig. 10. Variation of transverse shear stress ( ) through the thickness of beam at (x = 0, z) obtain using EE for aspect ratio 10.
a) The axial displacement ( )
The present theory gives realistic results of this displacement component in commensurate with the other shear deformation theories. For cantilever beam with various loads, the result of present theory are nearly matching with those other higher order theory.
b) The transverse deflection ( )
For cantilever beam with cosine load, the transverse deflection given by present theory is in excellent agreement with that of other higher order shear deformation theories.
c) The axial stress ( )
The axial stress and its distribution across the thickness given by present theory is in excellent agreement with that of higher order shear deformation theories.
d) The transverse shear stresses and
For cantilever beam with cosine load, transverse shear stress and its distribution through the thickness of beam obtained from constitutive relation are in close agreement with that of other higher order refined theories; however, use of constitutive relation cannot predict the effect of stress concentration at the built-in end of the beam. The effect of stress concentration on variation of transverse shear stress is exactly predicted by the present theory with the use of equilibrium equation of two dimensional elasticity. The realistic variations of these stresses at the built-in end of various beams are presented. Hence the use of equilibrium equation is inevitable to predict the effect stress concentration in accordance with the higher / equivalent refined shear deformation theories.
In general, the use of present theory gives accurate results as seen from the numerical examples studied and it is capable of predicting the local effects in the vicinity of the built-in end of the cantilever beam. This validates the efficacy and credibility of trigonometric shear deformation theory.