Scribd
Upload
55 views

National Taiwan University Finite Element Method: Lab Grading Policy: Basic 100%, Bonus 20%

This document provides instructions for a finite element method lab assignment. Students are asked to modify provided MATLAB code to include an interactive preprocessor for discretizing a domain using bilinear quadrilateral elements. They must generate input data files and output displacement results. For full credit, they are to add a convergence study plotting error versus mesh density. The document provides sample output and specifies requirements for submitting assignment files.

Uploaded by

agam9
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
55 views

National Taiwan University Finite Element Method: Lab Grading Policy: Basic 100%, Bonus 20%

This document provides instructions for a finite element method lab assignment. Students are asked to modify provided MATLAB code to include an interactive preprocessor for discretizing a domain using bilinear quadrilateral elements. They must generate input data files and output displacement results. For full credit, they are to add a convergence study plotting error versus mesh density. The document provides sample output and specifies requirements for submitting assignment files.

Uploaded by

agam9
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 3

National Taiwan University Finite Element Method

Department of Civil Engineering Instructor: T.-Y. Wu

Lab Assignment 7, 05/28/2020, 1800 -- 2000


Lab Grading Policy: Basic 100%, Bonus 20%

You are expected to complete the basic part during the Lab. In case you have difficulty in finishing
the basic part on time, you should upload them before 2100 on Saturday and a penalty of 20% discount
will be applied on your score. You are encouraged to complete the bonus part (no penalty applied).
Basic and/or bonus parts should be submitted by 2100 on Saturday. No late submission.

Download Sec4.2+2DQ4Taper.zip from the course website and unzip it. You will find a
folder containing a problem2dQ4Taper.m file with two data file and seven functions.

1. (100%) Consider the same trapezoidal panel given in the lecture note. The left edge is fixed. The
bottom and the right edges are traction free. Traction 𝑡𝑦 = −20 N/m is applied along the top
edge. Material and sectional properties are Young’s modulus 𝐸 = 3 × 107 Pa, Poisson’s ratio
𝑣 = 0.3 and thickness = 1 m. Assume plane stress conditions.

Modify the MATLAB code to include an interactive preprocessor to discretize the domain using
MxN bilinear Q4 quadrilateral elements. The modified code can output the input data file and also
the displacements. Name the input data file as Lab7+input.mat . Submit the modified code
as a zip file named Lab7+1.zip . TAs will run your code directly with arbitrary M and N to
see if your code can generate the corresponding data file and correct displacement results for
grading. Below is a sample run:

MxN = 2x2
Enter the discretization M and N: 2 2
Displacements
Node UX UY
1 0.0000e+00 0.0000e+00
2 -2.7321e -06 -6.7654e -06
3 -1.2329e -06 -1.8611e -05
4 0.0000e+00 0.0000e+00
5 4.2819e-07 -6.6965e-06
6 1.9583e-06 -1.8676e-05

1
7 0.0000e+00 0.0000e+00
8 3.9956e-06 -7.1385e-06
9 5.0961e-06 -1.8767e-05

MxN = 6x8
Enter the discretization M and N: 6 8
Displacements
Node UX UY
1 0.0000e+00 0.0000e+00
2 -1.5059e -06 -1.5078e -06
3 -2.7139e -06 -3.5074e -06
4 -3.5266e -06 -6.3540e -06
5 -3.8586e -06 -9.8898e -06
6 -3.7282e -06 -1.3937e -05
7 -3.1931e -06 -1.8308e -05
8 -2.3350e -06 -2.2810e -05
9 -1.2876e -06 -2.7261e -05
10 0.0000e+00 0.0000e+00
11 -1.0115e-06 -1.2789e -06
12 -1.7865e-06 -3.3574e -06
13 -2.1936e-06 -6.2442e -06
14 -2.2521e-06 -9.8186e -06
15 -1.9892e-06 -1.3896e -05
16 -1.4504e-06 -1.8291e -05
17 -7.0034e-07 -2.2808e -05
18 1.7661e-07 -2.7263e-05
19 0.0000e+00 0.0000e+00
20 -5.3390e-07 -1.1651e -06
21 -8.4436e-07 -3.2900e -06
22 -8.8144e-07 -6.2222e -06
23 -6.8454e-07 -9.8154e -06
24 -2.9706e-07 -1.3904e -05
25 2.4671e-07 -1.8305e-05
26 9.0417e-07 -2.2825e-05
27 1.6263e-06 -2.7276e-05
28 0.0000e+00 0.0000e+00
29 -4.0174e-08 -1.1386e -06
30 1.2035e-07 -3.3139e-06
31 4.5058e-07 -6.2753e-06
32 8.8374e-07 -9.8700e-06
33 1.3831e-06 -1.3952e-05
34 1.9281e-06 -1.8346e-05
35 2.4963e-06 -2.2858e-05
36 3.0696e-06 -2.7303e-05
37 0.0000e+00 0.0000e+00
38 4.9698e-07 -1.2001e-06
39 1.1564e-06 -3.4352e-06
40 1.8461e-06 -6.3971e-06
41 2.4914e-06 -9.9731e-06
42 3.0859e-06 -1.4036e-05
43 3.6205e-06 -1.8411e-05
44 4.0929e-06 -2.2905e-05
45 4.5156e-06 -2.7344e-05
46 0.0000e+00 0.0000e+00
47 1.1375e-06 -1.3870e-06
48 2.3453e-06 -3.6598e-06

2
49 3.3440e-06 -6.5774e-06
50 4.1786e-06 -1.0119e-05
51 4.8462e-06 -1.4148e-05
52 5.3499e-06 -1.8495e-05
53 5.7087e-06 -2.2965e-05
54 5.9709e-06 -2.7396e-05
55 0.0000e+00 0.0000e+00
56 2.1521e-06 -1.8027e-06
57 3.7152e-06 -3.9475e-06
58 4.9889e-06 -6.8171e-06
59 5.9858e-06 -1.0303e-05
60 6.7007e-06 -1.4284e-05
61 7.1437e-06 -1.8591e-05
62 7.3556e-06 -2.3033e-05
63 7.4366e-06 -2.7452e-05

2. (Bonus, 20%) You should finish Problem 1 first. Conduct a convergence study on Problem 1 by
setting M = N and change M from 1, 2, 4, 16, 32, and 64. Monitor the deflection along the y
direction on node 4. Use the answer from N = 64 as your reference solution and plot the difference
in % in the y-axis and M in the x-axis.

Submit the modified code as a zip file named Lab7+2.zip . TAs will run your code directly
to see if your code can generate the requested plot for grading. Below is a sample plot:

You might also like