M3 Test of FEM

M3 - Test of FEM module

The form functions

The frist thing that we have to take care is about the base element geometry, and it's contained in the file geometrie.py

For that, the page Base Element's Geometry has the definition of all the base elements geometries and show the functions and its derivatives.

And the points of integration are show in the page Gauss Integration Points.

Runing the code

This test is to verify if we calculate the values of F and E correctly. For that, we say that we know the value of u for each element and its value is given by

def fct(x):
    """
    Function that translates the position X into the vector U displacement
    Field:
        U(X1,X2) = X1*X2(a*e1 + b*e2)
    # for a=b=1 (Second and Third part of the exercise)
    """
    a = b = 1
    f0 = a * x[0] * x[1]
    f1 = b * x[0] * x[1]
    return [f0, f1]

As result, we have the figures given by

Result of F	Result of E