OptFEM3DP1 Toolbox  V1.0
Matlab/Octave Optimized P1-Lagrange Finite Element Method in 3D
 All Files Functions Variables Pages
valid_MECA.m
Go to the documentation of this file.
1 function Test=valid_MECA()
2 % Copyright:
3 % See \ref license
4 i=1;
5 Test(i).u={@(x,y,z) x - 2.*y, @(x,y,z) x + y - z, @(x,y,z) 3.*x + 2.*z};
6 Test(i).cu={'x - 2*y','x + y - z','3*x + 2*z'};
7 Test(i).v={@(x,y,z) x + 2.*y + 4.*z, @(x,y,z) 2.*x - y + 4.*z, @(x,y,z) 3.*x - 2.*y};
8 Test(i).cv={'x + 2*y + 4*z','2*x - y + 4*z','3*x - 2*y'};
9 Test(i).Mass=1.000000000000000;
10 Test(i).Stiff=30.00000000000000;
11 Test(i).lambda=1;
12 Test(i).mu=2;
13 Test(i).degree=1;
14 i=2;
15 Test(i).u={@(x,y,z) x.^2 - 2.*x.*y + x.*z, @(x,y,z) y.^2 - y.*z + z.^2 + x, @(x,y,z) x.^2 - x.*z - y.*z - z.^2};
16 Test(i).cu={'x^2 - 2*x*y + x*z','y^2 - y*z + z^2 + x','x^2 - x*z - y*z - z^2'};
17 Test(i).v={@(x,y,z) x.^2 + 2.*y.^2 - x.*z, @(x,y,z) 2.*x.^2 - x.*y + y.*z, @(x,y,z) x.*y - y.*z + z.^2};
18 Test(i).cv={'x^2 + 2*y^2 - x*z','2*x^2 - x*y + y*z','x*y - y*z + z^2'};
19 Test(i).Mass=0.4166666666666667;
20 Test(i).Stiff=-6.666666666666667;
21 Test(i).lambda=1;
22 Test(i).mu=2;
23 Test(i).degree=2;
24 i=3;
25 Test(i).u={@(x,y,z) x.^2 - 2.*x.*y + x.*z, @(x,y,z) x.^3 + y.^2 - y.*z + z.^2, @(x,y,z) -x.^2.*z - x.*y.*z + x.^2 - z.^2};
26 Test(i).cu={'x^2 - 2*x*y + x*z','x^3 + y^2 - y*z + z^2','-x^2*z - x*y*z + x^2 - z^2'};
27 Test(i).v={@(x,y,z) -x.*z.^2 + x.^2 + 2.*y.^2, @(x,y,z) 2.*x.^2 - x.*y + y.*z, @(x,y,z) x.*y - y.*z + z.^2};
28 Test(i).cv={'-x*z^2 + x^2 + 2*y^2','2*x^2 - x*y + y*z','x*y - y*z + z^2'};
29 Test(i).Mass=0.3375000000000000;
30 Test(i).Stiff=-2.958333333333333;
31 Test(i).lambda=1;
32 Test(i).mu=2;
33 Test(i).degree=3;
34 i=4;
35 Test(i).u={@(x,y,z) x.^4 - 2.*x.*y + x.*z, @(x,y,z) y.^2 - y.*z + z.^2 + x, @(x,y,z) -x.*y.*z + x.^2 - y.*z - z.^2};
36 Test(i).cu={'x^4 - 2*x*y + x*z','y^2 - y*z + z^2 + x','-x*y*z + x^2 - y*z - z^2'};
37 Test(i).v={@(x,y,z) x.^2 + 2.*y.^2 - x.*z, @(x,y,z) y.*z.^3 + 2.*x.^2 - x.*y, @(x,y,z) z.^4 + x.*y - y.*z};
38 Test(i).cv={'x^2 + 2*y^2 - x*z','y*z^3 + 2*x^2 - x*y','z^4 + x*y - y*z'};
39 Test(i).Mass=0.3277777777777778;
40 Test(i).Stiff=-6.270833333333333;
41 Test(i).lambda=1;
42 Test(i).mu=2;
43 Test(i).degree=4;
44 i=5;
45 Test(i).u={@(x,y,z) -2.*x.^3.*y.^2 + x.^2 + x.*z, @(x,y,z) y.^2 - y.*z + z.^2 + x, @(x,y,z) -x.^3.*z.^2 + x.^2 - y.*z - z.^2};
46 Test(i).cu={'-2*x^3*y^2 + x^2 + x*z','y^2 - y*z + z^2 + x','-x^3*z^2 + x^2 - y*z - z^2'};
47 Test(i).v={@(x,y,z) x.^5 + 2.*y.^2 - x.*z, @(x,y,z) 2.*x.^2 - x.*y + y.*z, @(x,y,z) x.*y.^4 - y.*z + z.^2};
48 Test(i).cv={'x^5 + 2*y^2 - x*z','2*x^2 - x*y + y*z','x*y^4 - y*z + z^2'};
49 Test(i).Mass=0.7736044973544974;
50 Test(i).Stiff=-5.226190476190476;
51 Test(i).lambda=1;
52 Test(i).mu=2;
53 Test(i).degree=5;
54 i=6;
55 Test(i).u={@(x,y,z) x.^6 - 2.*x.^3.*y.^2 + x.*z, @(x,y,z) y.^2 - y.*z + z.^2 + x, @(x,y,z) -x.^3.*z.^3 + x.^2 - y.*z - z.^2 - 1};
56 Test(i).cu={'x^6 - 2*x^3*y^2 + x*z','y^2 - y*z + z^2 + x','-x^3*z^3 + x^2 - y*z - z^2 - 1'};
57 Test(i).v={@(x,y,z) x.^5.*z + 2.*y.^2 - x.*z + 2, @(x,y,z) 2.*x.^2.*y.*z.^2 - x.*y + y.*z, @(x,y,z) x.*y.^4 - y.*z + z.^2};
58 Test(i).cv={'x^5*z + 2*y^2 - x*z + 2','2*x^2*y*z^2 - x*y + y*z','x*y^4 - y*z + z^2'};
59 Test(i).Mass=0.2857010582010582;
60 Test(i).Stiff=-3.547123015873016;
61 Test(i).lambda=1;
62 Test(i).mu=2;
63 Test(i).degree=6;
64 i=7;
65 Test(i).u={@(x,y,z) cos(1./2.*x + y - 1./2.*z), @(x,y,z) x.^2 - 1./2.*y + 1./2.*z - 1, @(x,y,z) cos(-1./3.*x + 2./3.*y + 1./3.*z)};
66 Test(i).cu={'cos(1/2*x + y - 1/2*z)','x^2 - 1/2*y + 1/2*z - 1','cos(-1/3*x + 2/3*y + 1/3*z)'};
67 Test(i).v={@(x,y,z) sin(1./2.*x + 1./2.*y - 1./4.*z), @(x,y,z) sin(4./3.*x + 1./3.*y - z + 1./3), @(x,y,z) 2./3.*y.^2 - 2./3.*x - 1./3.*z};
68 Test(i).cv={'sin(1/2*x + 1/2*y - 1/4*z)','sin(4/3*x + 1/3*y - z + 1/3)','2/3*y^2 - 2/3*x - 1/3*z'};
69 Test(i).Mass=-0.3329851257339365;
70 Test(i).Stiff=-0.3284478441464897;
71 Test(i).lambda=1;
72 Test(i).mu=2;
73 Test(i).degree=-1;