CHAMELEM_HOOKE.m

%--------------------------------------------------------------------------
% CHAMELEM.HOOK
%--------------------------------------------------------------------------
% Creation of a CHAMELEM describing the local stiffness matrix (Hooke,
% secant or tangent)

%--------------------------------------------------------------------------
% Syntax
%--------------------------------------------------------------------------

% 1. Computation of the Hooke matrix field
CHAM2 = CHAMELEM.HOOKE(MO1,CHAM1)

% 2. Computation of either the Hooke, secant or tangent matrix field
%    according to the value of label
CHAM2 = CHAMELEM.HOOKE(MO1,CHAM1,label,CHAM3)

%--------------------------------------------------------------------------
% Description
%--------------------------------------------------------------------------
% CHAMELEM.HOOKE allows to compute a field of local moduli. The local
% moduli can be the Hooke matrix, the secant matrix or the tangent matrix.
% This field is then used when computing the elemental stiffness matrix. 

%--------------------------------------------------------------------------
% Input arguments
%--------------------------------------------------------------------------
% MO1    : MODEL object
% CHAM1  : CHAMELEM object describing the material parameters
% label  : STRING object allowing to select the type of local stiffness to
%          be computed. label can be:  
%             - LINEAR   : Hooke matrix
%             - SECANT   : Secant matrix
%             - TANGENT  : Tangent matrix
%
% CHAM3  : CHAMELEM object describing the field of internal variables. This
%          field is mandatory when computing the secant and tangent matrix.

%--------------------------------------------------------------------------
% Output arguments
%--------------------------------------------------------------------------
% CHAM2   : CHAMELEM object

%--------------------------------------------------------------------------
% Example
%--------------------------------------------------------------------------
HO1 = CHAMELEM.HOOKE(MO1,MA1)

HO1 = 

  CHAMELEM with properties:

       model: [1x1 MODEL]
    cham_val: [1x2 struct]

HO1.cham_val(1).val_elem(1)

ans = 

    maho: {[3x3 double]}

HO1.cham_val(1).val_elem(1).maho{1}

ans =

   1.0e+10 *

    3.7500    0.7500         0
    0.7500    3.7500         0
         0         0    1.5000
    
HO1.cham_val(2).val_elem(1)

ans = 

    maho: {[3x3 double]  [3x3 double]  [3x3 double]  [3x3 double]}