OptFEM2DP1 Toolbox  V1.2
Matlab/Octave Optimized P1-Lagrange Finite Element Method in 2D
base/ElemStiffElasMatBb2DP1.m File Reference

Computation of the element stiffness elasticity matrix for $P_1$-Lagrange method. The method for numbering the degrees of freedom is local bloc numbering (classical method) More...

Go to the source code of this file.

Functions

function Ke = ElemStiffElasMatBb2DP1 (ql, area, C)
 Computation of the element stiffness elasticity matrix for $P_1$-Lagrange method. The method for numbering the degrees of freedom is local bloc numbering (classical method)

Detailed Description

Computation of the element stiffness elasticity matrix for $P_1$-Lagrange method. The method for numbering the degrees of freedom is local bloc numbering (classical method)

Definition in file ElemStiffElasMatBb2DP1.m.


Function Documentation

function Ke = ElemStiffElasMatBb2DP1 (   ql,
  area,
  C 
)

Computation of the element stiffness elasticity matrix for $P_1$-Lagrange method. The method for numbering the degrees of freedom is local bloc numbering (classical method)

Hooke's matrix
C=[L + 2*M L 0] [ L L + 2*M 0] [ 0 0 M]
Numbering of local points in reference element is
P=[(0, 0), (1, 0), (0, 1)]

Example:

        ql=[0 1 0;0 0 1];
        area=1/2.;
        lambda=1.; mu=1.;
        C=[lambda+2*mu, lambda, 0;lambda, lambda + 2*mu, 0;0, 0, mu];
        Elem=ElemStiffElasMatBb2DP1(ql,area,C);
      
Copyright
See License issues
Parameters:
qlarray of coordinates of the vertices of the triangle, 2-by-3 matrix (double)
areatriangle area (double)
CHooke's matrix (3-by-3 double)
Return values:
Elemelement stiffness elasticity matrix, 6-by-6 matrix (double)

Definition at line 17 of file ElemStiffElasMatBb2DP1.m.

 All Files Functions