OptFEM2DP1 Toolbox  V1.2b3
Matlab/Octave Optimized P1-Lagrange Finite Element Method in 2D
 All Files Functions Pages
ElemStiffElasMatP1BaOptV0.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 alternate numbering (classical method) More...

Go to the source code of this file.

Functions

function Ke = ElemStiffElasMatP1BaOptV0 (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 alternate 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 alternate numbering (classical method)

Definition in file ElemStiffElasMatP1BaOptV0.m.

Function Documentation

function Ke = ElemStiffElasMatP1BaOptV0 (   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 alternate 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=ElemStiffElasMatP1BaOptV0(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 ElemStiffElasMatP1BaOptV0.m.