Bahir Dar University
Bahir Dar Institute of Technology (BiT)
Faculty of Civil and Water Resource Engineering
Post Graduate Program in Structural Engineering
Finite Element Method in Structures
Assignment-3
Prepared by: Natinael Wolde
ID: BDU1022155
Submitted to: Instructor Seyfe N.
Submission Date: June 23, 2020
Finite Element Method in Structures Assignment -3
#Q-1.Solve the following equation using the Galerkin’s method (use at least a
two -parameter solution)
( (
Solution
Step-1:Assume a trial solution ,I observe that it is not easy to select a trial function that
satiesfies all the boundary condition.for now let us choose polynomial function
From boundary condition
Then
(
Step-2:Find the domain of residual
( then substitute to differential equation
( ( (
( (
( (
( (
Step-2:The minimize the residual ,the weighting function to be used in the Galerkin
technique are the trial function themselves i.e
( ( , ( (
Thus the weighted residual statements
∫ ( (
Master of Structural Engineering Page 1
Finite Element Method in Structures Assignment -3
Galerkin’s method with respect to (
∫ ( (
∫( ( (
∫ ( (
Integrate each by each
∫ ( ( |
∫ ( ( |
∫ ( |
∫ ( |
Then substitute
Galerkin’s method with respect to (
∫ ( (
∫( ( (
∫ ( (
∫ ( ( |
Master of Structural Engineering Page 2
Finite Element Method in Structures Assignment -3
∫ ( ( |
∫ ( |
∫ ( |
∫ ( |
Then substitute
Then substitute Equation -1 to Equation-2
From Equation-1 to substitute the value of C3
#Q-2:Develop the weak form and find the solution
v(x)=0, ( , ( (
Master of Structural Engineering Page 3
Finite Element Method in Structures Assignment -3
Solution:
Step-1: let as assume a trial solution
From boundary condition
(
(
( (
(
Step-2: Compute the weighted residual statement by weak form .The governing differential
equation is given by
∫ ( (
∫ ( ( ∫ ( ∫ ( ( ∫ (
Recall the standard formula for integration by parts (for a definite integral)
∫ ( ⌊ ∫ ( ⌊ ( ⌊ ∫
In our case ( integrating the first term by parts,we get
* ( + ∫( ∫ (
Integrating by parts once again,we obtain
Master of Structural Engineering Page 4
Finite Element Method in Structures Assignment -3
(
* ( + * + ∫( ∫ (
Substituing for the prescribed natural boundary condition ( ( ,we
have the resulting weak form:
(
( [ ] ∫( ∫ (
let W(x) be the weighting function.
From (
Let us assume weighting function ( and (
( and (
( and (
( , ( and ( , (
( , ( ( , (
( , ( and ( , (
Weak form with respect to (
(
( [ ] ∫( ∫ (
(
( ( [ ] [ ] ∫( ∫ (
* + ∫ ( ∫
Weak form with respect to (
Master of Structural Engineering Page 5
Finite Element Method in Structures Assignment -3
(
( [ ] ∫( ∫ (
(
( ( [ ] [ ] ∫( ∫ (
* + ∫ ( ∫
(
Then substitute Equation-1 to Equation-2
then collect like term
Then substitute Equation-3 to Equation-1
( (
(
(
Therefore our approximate solution to the given problem is
( ( ( ( ( )
( ( (
( *( ( +
Master of Structural Engineering Page 6
Finite Element Method in Structures Assignment -3
#Q-3: Find the solution for a deformation v(x) for a simply supported beam
under uniform load of qO as shown in Figure us ing Rayleigh‐Ritz (R‐R)
method.
Solution:
Step-1:Setting an approximation function either polynominal or trigonometric function,for
now let us take trigonometric function
Step-2: compute Strain Energy(U) , potential energy due to external force and total
a) compute Strain Energy(U)
∫ ( ) dx
∫( ) ∫( ) ∫( )
then substitute
(∫ ∫ )= ( | )|
b) Potential Energy due to External Force(V)
Master of Structural Engineering Page 7
Finite Element Method in Structures Assignment -3
∫ ∫ | (
c) Evaluate Total Potential Energy(
=Internal strain energy in the structure (U)+ Potential Energy due to External Force(V)
Step-3: Set up and solve the system of equations ,from the principle of stationary
total potential (PSTP),we have
Thus, The final solution is
Master of Structural Engineering Page 8
