| OptFEM2DP1 Toolbox 1.2b4 Matlab/Octave Optimized P1-Lagrange Finite Element Method in 2D

- Main Page
- Related Pages
- Files

- File List
- File Members

Functions
Opt/ElemStiffElasMatP1BaOptV0.m File Reference

Computation of the element stiffness elasticity matrix for -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 -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 -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 -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:

  ql | array of coordinates of the vertices of the triangle, 2-by-3 matrix (double)
  area | triangle area (double)
  C | Hooke's matrix (3-by-3 double)

Return values:

  Elem | element stiffness elasticity matrix, 6-by-6 matrix (double)

Definition at line 17 of file ElemStiffElasMatP1BaOptV0.m.

All Files Functions