![IJESM Volume 2, Issue 2 ISSN: 2320-0294
_________________________________________________________
A Quarterly Double-Blind Peer Reviewed Refereed Open Access International e-Journal - Included in the International Serial Directories
Indexed & Listed at: Ulrich's Periodicals Directory ©, U.S.A., Open J-Gage, India as well as in Cabell’s Directories of Publishing Opportunities, U.S.A.
International Journal of Engineering, Science and Mathematics
http://www.ijmra.us
72
June
2013
1. Introduction
A beam is a structural element or member designed to support loads applied at various points
along the element. Beams are one of the components used in structural engineering and can be in
1-dimension, 2-dimension or 3-dimension. They can be horizontal, vertical and also at angles.
We have analyzed a uniform elastic beam, simply supported length L and subjected to vertical
forces acting in the principal plane of a symmetrical cross-section as shown below.
Figure 1: A uniform elastic beam.
q – Centre concentrated load.
0q
– Uniform distributed load intensity.
L – The length AB.
)t,x(u – deflection at (x,t). where x is a one dimension spatial variable point at time t.
The flexture of the uniform elastic length L whose ends are simply supported can be modeled by
the equation:
txf
x
u
c
t
u
,4
4
2
2
2
(1.1)
subject to boundary conditions:
0,,0 tLutu (1.1a)
0,,0 tLutu xxxx (1.1b)
According to [12] ,finite element method has been used to approximate the deflection u(x) of a
simply supported beam. ZaferAhsan[17] applied Laplace transform method, Singh[12] used
Rayleigh method and Osongo[15] used direct integration method to the beam equation.](https://image.slidesharecdn.com/7ijmra-esm3219-170815140153/85/Finite-Element-Approach-to-the-Solution-of-Fourth-order-Beam-Equation-2-320.jpg)
![IJESM Volume 2, Issue 2 ISSN: 2320-0294
_________________________________________________________
A Quarterly Double-Blind Peer Reviewed Refereed Open Access International e-Journal - Included in the International Serial Directories
Indexed & Listed at: Ulrich's Periodicals Directory ©, U.S.A., Open J-Gage, India as well as in Cabell’s Directories of Publishing Opportunities, U.S.A.
International Journal of Engineering, Science and Mathematics
http://www.ijmra.us
73
June
2013
W.T.Thomson[14] employed mode-summation method which gave results in series form with no
immediate practical evaluation.
Since all the above mentioned methods presented their solutions in non closed form or single
degree of freedom or which lacked immediate practical evaluation, there was a need to look for
an alternative numerical approach to the two dimensional one degree of freedom beam equation
),(2
txfucu xxxxtt .
2 preliminaries
Basing on [6,7,10,11,12], Finite element method involves discretization, development of local
mass, stiffness and force matrices and then assembly to global load, stiffness and force matrix
using finite element method thereafter approximate resultant ordinary differential equations.
Consider the BVP
txf
x
u
EI
xt
u
,2
2
2
2
2
2
(2.1)
subject to boundary conditions,
(i) u(0,t)=u(L,t)=0
(ii) uxx (0,t)=uxx(L,t)=0 .
The function of approximation is given by:
xt)t,x(u i
i
i
5
1
(2.2)
where
xi is the shape function or basis function
and )(ti
the Fourier coefficients
then , solution is assumed to be in form;
L
xi
ttxu
i
i sin,
5
1
. (2.3)
and
0,, 2
2
2
2
2
2
txf
x
u
EI
xt
u
txrresidual (2.4)](https://image.slidesharecdn.com/7ijmra-esm3219-170815140153/85/Finite-Element-Approach-to-the-Solution-of-Fourth-order-Beam-Equation-3-320.jpg)
![IJESM Volume 2, Issue 2 ISSN: 2320-0294
_________________________________________________________
A Quarterly Double-Blind Peer Reviewed Refereed Open Access International e-Journal - Included in the International Serial Directories
Indexed & Listed at: Ulrich's Periodicals Directory ©, U.S.A., Open J-Gage, India as well as in Cabell’s Directories of Publishing Opportunities, U.S.A.
International Journal of Engineering, Science and Mathematics
http://www.ijmra.us
77
June
2013
+
11
10
9
8
7
6
5
4
3
2
1
3
3
2.39135612880221416035.10783.6599.3520.15427.46927.5
356106.2639209015541.10707.6663.3621.16099.48239.6
28806.26393126.184814029.9801.6199.3393.15156.46951.5
221420906.18485.24311787.8331.5317.2934.13158.402032.5
16031554140211781639.6544205.2339.10448.32165.4
5.10781.10709.9807.8339.6544.634.3033.16927.7666.230384.3
3.6597.6666191.5314204.30346.219.10959.4941.159825.1
9.3523.3629.3397.2935.2333.1699.10981.18849.27666.81152.1
0.1541.16025.1514.1319.10427.7659.4985.2784.12924.35049.0
72.4699.4856.4658.4048.3266.2341.15666.8924.3804.31572.0
927.5239.6951.52032.5165.40384.39825.11152.15049.01572.02939.0
2L
EI
=
0865.0
1.0
1057.0
1011.0
0840.0
0515.0
0
7725.0
1959.0
4045.0
9511.0
),( Ltxf
. (2.16)
From [5, 6], stability depends only on mass and stiffness of the system and is independent of the
number of subdivisions of the beam.
3 NUMERICAL SOLUTIONS
The approximate solution of Equation (2.15) is given by
)(
)(
)(
)(
)(
5.628285.04.251857.05.1120348.067.330067.0167.40033.0
4.251257.14366.157.1169495.067.421482.0033.60875.0
5.1121919.037.11695495.05.133232.154.300303.15.41515.0
67.330067.067.421482.054.300303.13199.1885.19428.0
167.40437.0033.60875.05.41515.0885.19428.05.0471.1
5
4
3
2
1
t
t
t
t
t
nmnmnmnmnm
nmmnmnmnm
nmnmnmnmnm
nmnmnmmnm
nmnmnmnmnm
=
5
),(7071.0
4
),(
3
),(7071.0
0
),(7071.0
txf
txf
txf
txf
(3.1)
4984.27,99.762,01.34373 nm
See([6],[12]).
Substituting the values above in thematrix (3.1) with boundary conditions (1.1a),](https://image.slidesharecdn.com/7ijmra-esm3219-170815140153/85/Finite-Element-Approach-to-the-Solution-of-Fourth-order-Beam-Equation-7-320.jpg)
![IJESM Volume 2, Issue 2 ISSN: 2320-0294
_________________________________________________________
A Quarterly Double-Blind Peer Reviewed Refereed Open Access International e-Journal - Included in the International Serial Directories
Indexed & Listed at: Ulrich's Periodicals Directory ©, U.S.A., Open J-Gage, India as well as in Cabell’s Directories of Publishing Opportunities, U.S.A.
International Journal of Engineering, Science and Mathematics
http://www.ijmra.us
81
June
2013
REFERENCES
[1] J.S.Rao(1992) ADVANCED THEORY OF VIBRATION Wiley Eastern Limited.
[2] Madhu B. Kanchi,(1993) Matrix Methods of Structural Analysis Wiley Eastern Limited. Pp
350-367.
[3] MajidMirzael, May 9,(2008) Finite Element Methods 1 lecture notes,
http://www.modares.ac.ir/eng/mmirzaei/finite.htm . Pp 1-40.
[4] P.W. Hemker,(2004) Discretisation of PDEs, Finite Element Method OWI
Amsterdam.Pp17-90.
[5] S.S. Sastry, (1988) Introduction Methods of Numerical Analysis 3rd
edition Prentice Hall
of India Private Ltd Pp 325-360 .
[6] Singh Sadhu,(2002) Theory of Elasticity by Khana publishers.Pp499-510,Pp555-
569,Pp609-640.
[7] W.T.Thomson, (1988) Theory of Vibration with Applications 3rd
edition CBS Publisher and
distributors Pp 233-240,Pp100-110.
[8] W.M. Onsongo,(1993) Analysis of Stastically Determinate Structure University of Nairobi
Press Pp 213.
[9] ZaferAhsan(2004) Differential Equations and their Applications Prentice-Hall of India.](https://image.slidesharecdn.com/7ijmra-esm3219-170815140153/85/Finite-Element-Approach-to-the-Solution-of-Fourth-order-Beam-Equation-11-320.jpg)
Finite Element Approach to the Solution of Fourth order Beam Equation:
- 1. IJESM Volume 2, Issue 2 ISSN: 2320-0294 _________________________________________________________ A Quarterly Double-Blind Peer Reviewed Refereed Open Access International e-Journal - Included in the International Serial Directories Indexed & Listed at: Ulrich's Periodicals Directory ©, U.S.A., Open J-Gage, India as well as in Cabell’s Directories of Publishing Opportunities, U.S.A. International Journal of Engineering, Science and Mathematics http://www.ijmra.us 71 June 2013 Finite Element Approach to the Solution of Fourth order Beam Equation: 2 ,tt xxxxu c u f x t Oduor Michael E. Okoya T. J. O. Aminer* Hagai Amakobe James Nthiiri Joyce Kagendo Abstract Finite element method is a class of mathematical tool which approximates solutions to initial and boundary value problems. Finite element, basic functions, stiffness matrices,systems of ordinary differential equations and hence approximate solutions of partial differential equations which involves rendering the partial differential equation into system of ordinary differential equations. The ordinary differential equations are then numerically integrated. We present a finite element approach in solving fourth order linear beam equation: 2 ,tt xxxxu c u f x t , which arises in model studies of building structures wave theory. In physical application of waves in building structures, coefficient 2 c , has the meaning of flexural rigidity per linear mass density and ,f x t external forcing term. In this paper, we give a solution to the beam equation with 2 139c and , 100.f x t Keywords: beam equation, finite element,approximation functions, stiffness matrix Department of Mathematics and Actuarial Science, JaramogiOgingaOdinga University of Science and Technology,Bondo, Kenya ACK Milimani Girls Secondary School, Musanda-Kenya Department of Mathematics, MasindeMuliro University of Science and Technology, Kakamega, Kenya
- 2. IJESM Volume 2, Issue 2 ISSN: 2320-0294 _________________________________________________________ A Quarterly Double-Blind Peer Reviewed Refereed Open Access International e-Journal - Included in the International Serial Directories Indexed & Listed at: Ulrich's Periodicals Directory ©, U.S.A., Open J-Gage, India as well as in Cabell’s Directories of Publishing Opportunities, U.S.A. International Journal of Engineering, Science and Mathematics http://www.ijmra.us 72 June 2013 1. Introduction A beam is a structural element or member designed to support loads applied at various points along the element. Beams are one of the components used in structural engineering and can be in 1-dimension, 2-dimension or 3-dimension. They can be horizontal, vertical and also at angles. We have analyzed a uniform elastic beam, simply supported length L and subjected to vertical forces acting in the principal plane of a symmetrical cross-section as shown below. Figure 1: A uniform elastic beam. q – Centre concentrated load. 0q – Uniform distributed load intensity. L – The length AB. )t,x(u – deflection at (x,t). where x is a one dimension spatial variable point at time t. The flexture of the uniform elastic length L whose ends are simply supported can be modeled by the equation: txf x u c t u ,4 4 2 2 2 (1.1) subject to boundary conditions: 0,,0 tLutu (1.1a) 0,,0 tLutu xxxx (1.1b) According to [12] ,finite element method has been used to approximate the deflection u(x) of a simply supported beam. ZaferAhsan[17] applied Laplace transform method, Singh[12] used Rayleigh method and Osongo[15] used direct integration method to the beam equation.
- 3. IJESM Volume 2, Issue 2 ISSN: 2320-0294 _________________________________________________________ A Quarterly Double-Blind Peer Reviewed Refereed Open Access International e-Journal - Included in the International Serial Directories Indexed & Listed at: Ulrich's Periodicals Directory ©, U.S.A., Open J-Gage, India as well as in Cabell’s Directories of Publishing Opportunities, U.S.A. International Journal of Engineering, Science and Mathematics http://www.ijmra.us 73 June 2013 W.T.Thomson[14] employed mode-summation method which gave results in series form with no immediate practical evaluation. Since all the above mentioned methods presented their solutions in non closed form or single degree of freedom or which lacked immediate practical evaluation, there was a need to look for an alternative numerical approach to the two dimensional one degree of freedom beam equation ),(2 txfucu xxxxtt . 2 preliminaries Basing on [6,7,10,11,12], Finite element method involves discretization, development of local mass, stiffness and force matrices and then assembly to global load, stiffness and force matrix using finite element method thereafter approximate resultant ordinary differential equations. Consider the BVP txf x u EI xt u ,2 2 2 2 2 2 (2.1) subject to boundary conditions, (i) u(0,t)=u(L,t)=0 (ii) uxx (0,t)=uxx(L,t)=0 . The function of approximation is given by: xt)t,x(u i i i 5 1 (2.2) where xi is the shape function or basis function and )(ti the Fourier coefficients then , solution is assumed to be in form; L xi ttxu i i sin, 5 1 . (2.3) and 0,, 2 2 2 2 2 2 txf x u EI xt u txrresidual (2.4)
- 4. IJESM Volume 2, Issue 2 ISSN: 2320-0294 _________________________________________________________ A Quarterly Double-Blind Peer Reviewed Refereed Open Access International e-Journal - Included in the International Serial Directories Indexed & Listed at: Ulrich's Periodicals Directory ©, U.S.A., Open J-Gage, India as well as in Cabell’s Directories of Publishing Opportunities, U.S.A. International Journal of Engineering, Science and Mathematics http://www.ijmra.us 74 June 2013 Multiplying r(x,t) by test function v(x) and varying the integral to weak form and integrate the product of r(x,t) and v(x) and equate to zero; 0),( 1 2 2 2 2 2 2 dxtxfv x u EI x v t u v e e x x (2.5) If assumed approximate solution is; 5 1 ),( i ii txtxu (2.6) then; 5 1 ),( i ii txtx t u (2.7) 5 1 2 2 i ii tx t u (2.8) 5 1 2 2 2 2 ˆ i i IV i tx x u x . (2.9) Substituting equations (2.8) and (2.9) into (2.5) we obtain: dxtxfvtxvEItxv e e x x i i IV i i ii 1 ),( 5 1 5 1 =0 (2.10) Equation (2.10) can be divided into three major parts (I,II and III): 5 1 )()()( 1 i ii x x dxtxxvvI e e (2.11) this gives mass matrix e e x x ji i iij dxxxtC 1 )()()( 5 1 in matrix form we have: 5 4 3 2 1 1.06285.05.03367.01667.0 6285.008081.06667.03771.0 5.08081.01667.08485.05.0 3367.06667.08485.004714.0 1667.03771.05.04714.05.0 2 L Cij
- 5. IJESM Volume 2, Issue 2 ISSN: 2320-0294 _________________________________________________________ A Quarterly Double-Blind Peer Reviewed Refereed Open Access International e-Journal - Included in the International Serial Directories Indexed & Listed at: Ulrich's Periodicals Directory ©, U.S.A., Open J-Gage, India as well as in Cabell’s Directories of Publishing Opportunities, U.S.A. International Journal of Engineering, Science and Mathematics http://www.ijmra.us 75 June 2013 e e x x i i iv i dxtxvII 1 5 1 )()( (2.12) gives the stiffness matrix dx dx ud )t( dx d )x(vEIK i i X X ij e e 2 25 1 2 2 1 . In matrix form the stiffness matrix is; 5 4 3 2 1 3 3 5.6241.2515.11267.33167.4 41.251037.11667.42033.6 5.11237.1165.1354.305.4 67.3367.4254.300885.1 167.4033.65.4885.15.0 2L EI Kij e e X X dx)t,x(vfIII 1 (2.13) we obtain the nodal force vector dxtxfxvF e e x x j ),()( 1 In matrix form wehave ),(14142.0 ),(25.0 ),(2357.0 0 ),(.7071.0 txf txf txf txf L F (2.14) The equilibrium equation is;
- 6. IJESM Volume 2, Issue 2 ISSN: 2320-0294 _________________________________________________________ A Quarterly Double-Blind Peer Reviewed Refereed Open Access International e-Journal - Included in the International Serial Directories Indexed & Listed at: Ulrich's Periodicals Directory ©, U.S.A., Open J-Gage, India as well as in Cabell’s Directories of Publishing Opportunities, U.S.A. International Journal of Engineering, Science and Mathematics http://www.ijmra.us 76 June 2013 5 4 3 2 1 1.06285.05.03367.01667.0 6285.008081.06667.03771.0 5.08081.01667.08485.05.0 3367.06667.08485.004714.0 1667.03771.05.04714.05.0 2 L + 5 4 3 2 1 3 3 5.6241.2515.11267.33167.4 41.251037.11667.42033.6 5.11237.1165.1354.305.4 67.3367.4254.300885.1 167.4033.65.4885.15.0 2L EI = ),(14142.0 ),(25.0 ),(2357.0 0 ),(.7071.0 txf txf txf txf L . (2.15) A similar equation with ten subdivisions(case b) with have an equilibrium equation of the form 11 10 9 8 7 6 5 4 3 2 1 0267.02943.02939.02860.02705.02476.02179.01823.01414.00965.00490.0 2943.003253.03266.03173.02972.02667.02264.01778.01225.00624.0 2939.03253.03417.03567.03533.03364.03057.02622.02075.01437.00735.0 2860.03266.03567.00594.03757.03618.03319.02868.02281.01585.00813.0 2705.03173.03533.03757.00679.03713.03429.02978.02378.01657.0085.0 2476.02972.03364.03618.03713.00490.03371.02939.02354.01643.00844.0 2179.02667.03057.03319.03429.03371.002747.02204.01541.00793.0 1823.02264.02622.02868.02978.02939.02747.00735.01934.01354.00697.0 1414.01778.02075.02281.02378.02354.02204.01934.01585.01090.00561.0 0965.01225.01437.01585.01657.01643.01541.01354.01090.02378.00393.0 0490.00624.00735.00813.0085.00844.00793.00697.00561.00393.02939.0 2 L
- 7. IJESM Volume 2, Issue 2 ISSN: 2320-0294 _________________________________________________________ A Quarterly Double-Blind Peer Reviewed Refereed Open Access International e-Journal - Included in the International Serial Directories Indexed & Listed at: Ulrich's Periodicals Directory ©, U.S.A., Open J-Gage, India as well as in Cabell’s Directories of Publishing Opportunities, U.S.A. International Journal of Engineering, Science and Mathematics http://www.ijmra.us 77 June 2013 + 11 10 9 8 7 6 5 4 3 2 1 3 3 2.39135612880221416035.10783.6599.3520.15427.46927.5 356106.2639209015541.10707.6663.3621.16099.48239.6 28806.26393126.184814029.9801.6199.3393.15156.46951.5 221420906.18485.24311787.8331.5317.2934.13158.402032.5 16031554140211781639.6544205.2339.10448.32165.4 5.10781.10709.9807.8339.6544.634.3033.16927.7666.230384.3 3.6597.6666191.5314204.30346.219.10959.4941.159825.1 9.3523.3629.3397.2935.2333.1699.10981.18849.27666.81152.1 0.1541.16025.1514.1319.10427.7659.4985.2784.12924.35049.0 72.4699.4856.4658.4048.3266.2341.15666.8924.3804.31572.0 927.5239.6951.52032.5165.40384.39825.11152.15049.01572.02939.0 2L EI = 0865.0 1.0 1057.0 1011.0 0840.0 0515.0 0 7725.0 1959.0 4045.0 9511.0 ),( Ltxf . (2.16) From [5, 6], stability depends only on mass and stiffness of the system and is independent of the number of subdivisions of the beam. 3 NUMERICAL SOLUTIONS The approximate solution of Equation (2.15) is given by )( )( )( )( )( 5.628285.04.251857.05.1120348.067.330067.0167.40033.0 4.251257.14366.157.1169495.067.421482.0033.60875.0 5.1121919.037.11695495.05.133232.154.300303.15.41515.0 67.330067.067.421482.054.300303.13199.1885.19428.0 167.40437.0033.60875.05.41515.0885.19428.05.0471.1 5 4 3 2 1 t t t t t nmnmnmnmnm nmmnmnmnm nmnmnmnmnm nmnmnmmnm nmnmnmnmnm = 5 ),(7071.0 4 ),( 3 ),(7071.0 0 ),(7071.0 txf txf txf txf (3.1) 4984.27,99.762,01.34373 nm See([6],[12]). Substituting the values above in thematrix (3.1) with boundary conditions (1.1a),
- 8. IJESM Volume 2, Issue 2 ISSN: 2320-0294 _________________________________________________________ A Quarterly Double-Blind Peer Reviewed Refereed Open Access International e-Journal - Included in the International Serial Directories Indexed & Listed at: Ulrich's Periodicals Directory ©, U.S.A., Open J-Gage, India as well as in Cabell’s Directories of Publishing Opportunities, U.S.A. International Journal of Engineering, Science and Mathematics http://www.ijmra.us 78 June 2013 0)()( 51 tt we obtain; 46.687 14.648 0 )( )( )( 27.4938024.5611770.27462 64.5596473.5578279.12112 70.2746279.1211294.45368 4 3 2 t t t Using mat lab we obtain solutions as: .0021783.0)(,0097102.0)(,0012739.0)( 432 ststst (3.2) From (2.6), we obtain the deflection equation as : xxxtxu sin0021783.075.0sin0097102.05.0sin0012739.0),( (3.3) The approximate solution of Equation (2.16) in matrix form is given by nmnmnmnmnmnmnmnmnmnmnm nmmnmnmnmnmnmnmnmnmnm nmnmnmnmnmnmnmnmnmnmnm nmnmnmnmnmnmnmnmnmnmnm nmnmnmnmnmnmnmnmnmnmnm nmnmnmnmnmnmnmnmnmnmn nmnmnmnmnmnmnmnmnmnmnm nmnmnmnmnmnmnmnmnmnmnm nmnmnmnmnmnmnmnmnmnmnm nmnmnmnmnmnmnmnmnmnmnm nmnmnmnmnmnmnmnmnmmnnm 3912.03561268.02880008.02214008.01603007.01079007.0659006.0353005.0154004.07.46003.093.5001.0 2880589.0620.02640324.0209901.0155501.0107001.0667009.0362008.0160007.049005.024.6002.0 2880263.02040015.0312001.0184602.01402014.0981014.0619013.034001.0151009.06.40006.06003.0 2214245.0209001.01849327.02446.0117833.0834016.0531015.0294014.0131011.06.40008.02.5004.0 1603224.0155501.01402014.0117833.01636.0655332.0420017.0234015.0105012.05.32009.017.4004.0 1079198.0107001.0981014.0834016.065511.05.63252.030307.0169005.03.76013.07.32009.004.3005.0 659169.066701.0619013.0531015.0420017.0303331.022612.0110329.06.49012.04.15009.098.1005.0 353138.0362373.0340011.0294014.0233015.0169015.0110329.019615.09.27325.066.8008.012.1004.0 154105.0160067.0151009.0131011.0105012.03.76013.05.49012.09.27325.08.1262.092.3320.0505.0003.0 7.46071.049005.06.40006.06.40008.05.32009.07.23009.04.15009.067.8008.092.3320.08.3624.0157.0316.0 93.5036.024.6002.095.5003.020.5004.017.4004.004.3005.098.1004.012.1054.0505.0003.0271.0157.0294.0549.0 0865.0).,( 1.0).,(. 1057.0).,(. 1011.0).,(. 0840.0).,(. 0515.0).,(. 0 7725.0).,(. 1959.0).,(. 4045.0).,(. 9511.0).,(. 11 10 9 8 7 6 5 4 3 2 1 txf txf txf txf txf txf txf txf txf txf
- 9. IJESM Volume 2, Issue 2 ISSN: 2320-0294 _________________________________________________________ A Quarterly Double-Blind Peer Reviewed Refereed Open Access International e-Journal - Included in the International Serial Directories Indexed & Listed at: Ulrich's Periodicals Directory ©, U.S.A., Open J-Gage, India as well as in Cabell’s Directories of Publishing Opportunities, U.S.A. International Journal of Engineering, Science and Mathematics http://www.ijmra.us 79 June 2013 4984.27,99.762,01.34373 nm and sssssss 010725.0,010774.0,012808.0,013978.0,420568.0,228128.1,699109.0 8765432 ss 013529.0,012832.0 109 From (2.6) deflection equation is given by: xxxxxtxu 5.1sin012808.025.1sin0139778.0sin420568.075.0sin228128.15.0sin699109.0),( xxxx 5.2sin013529.025.2sin012832.02sin010725.075.1sin010774.0 (3.4) Graphical output of Equation (3.3) Figure 2: Deflections ))((,),( mxlengthagainsttxu Graphical output of Equation (3.4) - 0 . 0 1 - 0 . 0 0 5 0 0 . 0 0 5 0 . 0 1 0 . 0 1 5 0 0 . 5 1 1 . 5 2 2 . 5 3 3 . 5 4
- 10. IJESM Volume 2, Issue 2 ISSN: 2320-0294 _________________________________________________________ A Quarterly Double-Blind Peer Reviewed Refereed Open Access International e-Journal - Included in the International Serial Directories Indexed & Listed at: Ulrich's Periodicals Directory ©, U.S.A., Open J-Gage, India as well as in Cabell’s Directories of Publishing Opportunities, U.S.A. International Journal of Engineering, Science and Mathematics http://www.ijmra.us 80 June 2013 Figure 3: Deflections ),( txu against length(x)(m) . 4 CONCLUSIONS We have solved the beam equation usingthefinite element method and the graphical outputs (fig 2) and ( fig 3) show that the solution obtained satisfy boundary conditions hence the solutions within Ltxu ),(0 are assumed consistent with exact values. Deflections in (fig 2) and fig (3) portray beam deflection of a similar pattern. The deflections in the graphical output smoothens as the number of subdivisions increases. - 1 . 5 - 1 - 0 . 5 0 0 . 5 1 1 . 5 2 0 0 . 5 1 1 . 5 2 2 . 5 3 3 . 5 4
- 11. IJESM Volume 2, Issue 2 ISSN: 2320-0294 _________________________________________________________ A Quarterly Double-Blind Peer Reviewed Refereed Open Access International e-Journal - Included in the International Serial Directories Indexed & Listed at: Ulrich's Periodicals Directory ©, U.S.A., Open J-Gage, India as well as in Cabell’s Directories of Publishing Opportunities, U.S.A. International Journal of Engineering, Science and Mathematics http://www.ijmra.us 81 June 2013 REFERENCES [1] J.S.Rao(1992) ADVANCED THEORY OF VIBRATION Wiley Eastern Limited. [2] Madhu B. Kanchi,(1993) Matrix Methods of Structural Analysis Wiley Eastern Limited. Pp 350-367. [3] MajidMirzael, May 9,(2008) Finite Element Methods 1 lecture notes, http://www.modares.ac.ir/eng/mmirzaei/finite.htm . Pp 1-40. [4] P.W. Hemker,(2004) Discretisation of PDEs, Finite Element Method OWI Amsterdam.Pp17-90. [5] S.S. Sastry, (1988) Introduction Methods of Numerical Analysis 3rd edition Prentice Hall of India Private Ltd Pp 325-360 . [6] Singh Sadhu,(2002) Theory of Elasticity by Khana publishers.Pp499-510,Pp555- 569,Pp609-640. [7] W.T.Thomson, (1988) Theory of Vibration with Applications 3rd edition CBS Publisher and distributors Pp 233-240,Pp100-110. [8] W.M. Onsongo,(1993) Analysis of Stastically Determinate Structure University of Nairobi Press Pp 213. [9] ZaferAhsan(2004) Differential Equations and their Applications Prentice-Hall of India.