SRI LANKA INSTITUE OF INFORMATION

TECHNOLOGY

Faculty of Engineering Department of civil

Finite Element Methos in Structural Engineering
CE 4011

EN20408510 – Bandara J.M.T.S

1
 Assessment Cover Sheet
Faculty of Engineering

Module Details

Module Code CE 4011 Module Title Finite Element Methos in Structural Engineering

Program: SLIIT Course: BSc

Stream: Civil

Assessment details

Title Finite Element Methos in Structural Group assignment NO

Engineering computer lab 03
If yes, Group No.

Lecturer/ Instructor Date of Performance

Due date 17/04/2023 Date submitted 17/04/2023

Student statement and signature

By this declaration, I/we confirm my/our understanding and acceptance that the work reported in this report is my/our own work. I/we also understand the
consequences of engaging in plagiarism or copying others work without proper citation. Any material used in this work (whether from published sources, the
internet or elsewhere) have been fully acknowledged and referenced and are without fabrication or falsification of data.

[Copying or plagiarism will result in a “0” mark for the continuous assessment and “F” for the module after an investigation on academic misconduct;
All academic misconduct is considered seriously and defined as dishonest and in direct opposition to the values of a learning community.
Misconduct may result in penalties from failure to exclusion from the campus.
Further help and guidance on how to avoid academic misconduct can be obtained from your academic advisor/tutor]

By this declaration, I/we confirm my understanding and acceptance that-

• I/we have adhered to relevant ethical guidelines and procedures in the completion of the assignment.
• I/we have not allowed another student to have access to or copy from this work. This work has not been submitted previously.
[The Institute may request an electronic copy of this work for submission to the Plagiarism detection facility (TURNITIN). You must make sure that an electronic
copy of your work is available in these circumstances]

Details of the student/s submitting the assignment Signature

ID Number Name (As per the institute records )

EN20408510 BANDARA J.M.T.S

OFFICE USE ONLY

Receiving Officer Specific comments about the work (including overall comments and guidelines
(seal, signature, date) for improvement)

Tutor: Signature: Date:

Marks: [ All marks are subject to external moderation and approval of board of
examinations]

2
 CE 4011- Finite Element Methods in Structural Engineering
Week – 9 - Computer Lab Tutorial
PROBLEM 01

The coordinates of gridlines are,

angle in angle in
Point radius X(R cos) Y( R sin)
degree radians
1 0 0 0.02 0.0200 0.0000
2 22.5 0.392699082 0.02 0.0185 0.0077
3 45 0.785398163 0.02 0.0141 0.0141
4 67.5 1.178097245 0.02 0.0077 0.0185 0-90 up
5 90 1.570796327 0.02 0.0000 0.0200 rhs
6 0 0 0.02 0.0200 0.0000
7 -22.5 -0.392699082 0.02 0.0185 -0.0077
8 -45 -0.785398163 0.02 0.0141 -0.0141
9 -67.5 -1.178097245 0.02 0.0077 -0.0185 0-90 down
10 -90 -1.570796327 0.02 0.0000 -0.0200 rhs
11 90 1.570796327 0.02 0.0000 0.0200
11 112.5 1.963495408 0.02 -0.0077 0.0185
12 135 2.35619449 0.02 -0.0141 0.0141
13 157.5 2.748893572 0.02 -0.0185 0.0077 0-90 up
14 180 3.141592654 0.02 -0.0200 0.0000 lhs
1.1. Model the plate with a mesh of about 10 x 40 quad elements.

1.1.a. Plot the variation of normal stress σx along the vertical axis of symmetry.
Y Coordinates σx (KN/m2)
0.1 48715.352
0.08 52352.86
0.06 53855.202
0.04 58557.274
0.02 120782.26
0 0
-0.02 120782.26
-0.04 58557.274
-0.06 53855.202
-0.08 52352.86
-0.1 48715.352
 σx (KN/m2)
140000
120782.259 120782.259
120000

100000
σx (KN/m2)

80000
58557.274 58557.274
53855.202
52352.86 60000 53855.202
52352.86
48715.352 48715.352

40000

20000
0
0
-0.15 -0.1 -0.05 0 0.05 0.1 0.15
y Coordinate

1.1.b. What is the largest value of σx and where does it occur?

Largest value - σx is 120782.26KNm-2
Location - at (0,0.020) and (0,-0.020)
1.2. Using symmetry, analyze a quarter of the plate with the same mesh as in (1.1) above.
Repeat (a) and (b) of (1.1) Compare the results of (1.1) and (1.2). Are they the same?

Values obtained from quarter of the plate.

Y Coordinates σx (KN/m2)

0.1 48737.222
0.08 52262.062
0.06 54225.138
0.04 58749.285
0.02 113205.34
0 0
 σx (KN/m2)
120000 113205.34

100000

80000

58749.285
60000 54225.138 52262.062
48737.222

40000

20000

0
0
0 0.02 0.04 0.06 0.08 0.1 0.12

Largest value - σx is 120532.19 KNm-2

Location at - (0,0.02)

Comparison of the results of the above two.

Y Coordinates σx (KN/m2) σx (KN/m2)

(Full plate) (Quarter plate)

0.1 48715.352 48737.222
0.08 52352.86 52262.062
0.06 53855.202 54225.138
0.04 58557.274 58749.285
0.02 120782.259 113205.34
0 0 (HOLE) 0 (HOLE)
 COMPARISON OF STRESS
140000

120782.259
σx (KN/m2) (Full plate) σx (KN/m2) (Quarter plate)

113205.34
120000

100000

80000

58749.285
58557.274

54225.138
53855.202

52262.062
52352.86
60000

48737.222
48715.352
40000

20000

0
0

0 0.02 0.04 0.06 0.08 0.1 0.12

When comparing the results of two cases, it is evident that they are nearly identical, which can
be attributed to the symmetry of the structure. Therefore, when analyzing a symmetrical
structure, it is not necessary to model the entire structure. By modeling only, the symmetrical
component and conducting an analysis on it, we can obtain results that can be extrapolated to
the entire structure. Since the structure being analyzed is symmetrical, it can be assumed that
analyzing only the symmetrical component is sufficient for the analysis of the entire structure.
1.3. Refine the mesh of the quarter plate till the analysis results converge. You may use
the value of the largest stress for checking convergence. Show the results (number of
elements and the value of largest stress) in a table. For the final FEM model, plot the
variation of σx along the vertical axis of symmetry. Estimate the stress concentration
factor for the hole based on the analysis results.

80 x 20 Quad Elements
When using a mesh with 80 x 20 Quad Elements it gives 175492.27 KN/m2 as the maximum
stress.