Scribd
Upload
15 views

Variational Method 2

The document compares solutions to eigenvalue and finite element method problems using different orders of polynomials as trial solutions and different numbers of elements. Graphs are shown comparing the different solutions.

Uploaded by

Structural Spreadsheets
Copyright
© © All Rights Reserved
Available Formats
Download as XLSX, PDF, TXT or read online on Scribd
15 views

Variational Method 2

The document compares solutions to eigenvalue and finite element method problems using different orders of polynomials as trial solutions and different numbers of elements. Graphs are shown comparing the different solutions.

Uploaded by

Structural Spreadsheets
Copyright
© © All Rights Reserved
Available Formats
Download as XLSX, PDF, TXT or read online on Scribd
You are on page 1/ 2

COMPARISON GRAPHS:

EIGEN VALUE PROBLEM OF VARIATIONAL METHOD:


d2U/dx2 + Λu = 0 for 0<x<1, at x = 0, U = 0; at x = 1, U = 0

[a] Solution using 2nd order polynomial as approximate trial solution:

U1(x) = -x+x2 for 0<x<1 (1st Eigen Vector)

[b] Solution using 3rd order polynomial as approximate trial solution:

U1(x) = -x+x2 for 0<x<1 (1st Eigen Vector)

U2(x) = -x+x2 + (2/3)(-x+x3) for 0<x<1 (2nd Eigen Vector)

EIGEN VALUE PROBLEM


0.000
x U1(x) - [a] U1(x) - [b] U2(x) - [b] 0.00 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 1.00
0.00 0.000 0.000 0.000 -0.100
0.10 -0.090 -0.090 -0.156
-0.200
0.20 -0.160 -0.160 -0.288
0.30 -0.210 -0.210 -0.392 -0.300
0.40 -0.240 -0.240 -0.464
0.50 -0.250 -0.250 -0.500 -0.400
0.60 -0.240 -0.240 -0.496
0.70 -0.210 -0.210 -0.448 -0.500

0.80 -0.160 -0.160 -0.352


-0.600
0.90 -0.090 -0.090 -0.204
1.00 0.000 0.000 0.000 U1(x) - [a] U1(x) - [b] U2(x) - [b]

FINITE ELEMENT METHOD:


d2U/dx2 + x2 = 0 for 0<x<1, at x = 0, U = 1; at x = 1, dU/dx + 2U = 1

[a] Solution assuming one element:

U(x) = 1 - 0.25x for 0<x<1

[b] Solution assuming two elements:

U(x) = 1 - 0.177x for 0<x<1/2 U(x) = 1.073 - 0.323x for 1/2<x<1

[c] Solution using classical method as exact solution:

U(x) = 1 - x/6 - x4/12 for 0<x<1 FINITE ELEMENT METHOD


1.200
x U(x) - [a] U(x) - [b] U(x) - [c]
0.00 1.000 1.000 1.000 1.000
0.10 0.975 0.982 0.983 0.800
0.20 0.950 0.965 0.967
0.30 0.925 0.947 0.949 0.600
0.40 0.900 0.929 0.931
0.400
0.50 0.875 0.912 0.911
0.60 0.850 0.879 0.889 0.200
0.70 0.825 0.847 0.863
0.80 0.800 0.815 0.833 0.000
0.00 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 1.00
0.90 0.775 0.782 0.795
U(x) - [a] U(x) - [b] U(x) - [c]
0.200

0.000
0.00 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 1.00

1.00 0.750 0.750 0.750 U(x) - [a] U(x) - [b] U(x) - [c]

You might also like