CHAMELEM_COMP.m

%--------------------------------------------------------------------------
% CHAMELEM.COMP
%--------------------------------------------------------------------------
% Creation of a CHAMELEM describing the stress field and the updated
% internal variables

%--------------------------------------------------------------------------
% Syntax
%--------------------------------------------------------------------------
CHAM4 = CHAMELEM.COMP(MO1,opt,CHAM1,CHAM2,CHAM3,istep)

%--------------------------------------------------------------------------
% Description
%--------------------------------------------------------------------------
% CHAMELEM.COMP allows to compute a stress field and to update the internal 3
% variables. In case of TRUSS elements, BEAM-like elements, T3G, etc., 
% the stresses are the generalized stresses.

%--------------------------------------------------------------------------
% Input arguments
%--------------------------------------------------------------------------
% MO1    : MODEL object
% opt    : integer, index to be set up according to the solving methodology
%          used to solve the global problem. opt can be equal to: 
%           - 1 : Compute the stress (standard procedure)
%           - 2 : Compute the stress for a given field of internal
%                 variables. That is to, the internal variables do not 
%                 evolve.
% CHAM1  : CHAMELEM object describing the material parameters
% CHAM2  : CHAMELEM object describing the strain field
% CHAM3  : CHAMELEM object describing the internal variable field
% istep  : integer, index to select the appropriate computation in case of
%          a non-local. istep can be equal to: 
%            - 0 : local computation
%            - 1 : computation of the averaged quantities
%            - 2 : computation of the stress, given the averaged quantities

%--------------------------------------------------------------------------
% Output arguments
%--------------------------------------------------------------------------
% CHAM4  : CHAMELEM object describing the stress field and the
%          corresponding internal variables

%--------------------------------------------------------------------------
% Example
%--------------------------------------------------------------------------
SOL

SOL = 

  SOLUTION with properties:

            problem: [1x1 PROBLEM]
               mesh: [1x1 MESH]
           topology: [1x1 TOPOLOGY]
       displacement: [1x100 CHPOINT]
           velocity: []
       acceleration: []
       cham_element: [1x100 struct]
           reaction: [1x100 CHPOINT]
    internal_forces: [1x100 CHPOINT]
     eigenfrequency: []
          modeshape: []
          
U = SOL.displacement(end)

U = 

  CHPOINT with properties:

    support: [1x602 NODE]
    loc_dof: [602x1 double]
    dof_val: [602x1 double]

CHAM1 = SOL.cham_element(end).cham

CHAM1 = 

  CHAMELEM with properties:

       model: [1x1 MODEL]
    cham_val: [1x2 struct]
    
CHAM_EPS = CHAMELEM.EPSI(MOT,U.dof_val,MAT)

CHAM_EPS = 

  CHAMELEM with properties:

       model: [1x1 MODEL]
    cham_val: [1x2 struct]
          
CHAM_HOOKE = CHAMELEM.HOOKE(MOT,MAT)

CHAM_HOOKE = 

  CHAMELEM with properties:

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

CHAM_STRESS = CHAMELEM.COMP(MOT,1,[MAT CHAM_HOOKE],CHAM_EPS,CHAM1,0)

CHAM_STRESS = 

  CHAMELEM with properties:

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