EVOL_REAC.m

%--------------------------------------------------------------------------
% EVOL.REAC
%--------------------------------------------------------------------------
% Computation of the reaction curve in case of an incremental iterative 
% analysis

%--------------------------------------------------------------------------
% Syntax
%--------------------------------------------------------------------------
EV2 = EVOL.REAC(SOL1,BLO1,EV2,idx)
EV2 = EVOL.REAC(SOL1,BLO1,EV2,idx,label)

%--------------------------------------------------------------------------
% Description
%--------------------------------------------------------------------------
% The method REAC of the EVOL class allows for computing the reaction curve 
% in case of an incremental iterative analysis

%--------------------------------------------------------------------------
% Input arguments
%--------------------------------------------------------------------------
% SOL1    : SOLUTION object coming fom the SOLVE method
% BLO1    : MATRICE object describing the Dirichlet boundary conditions
% EV2     : EVOL object discribing the time evolution of the loading
% idx     : INTEGER object describing the local type of the degrees of 
%           freedom. It is used to to select the appropriate abscissa axis.
%           idx can be : 
%              - 1 : Direction 1 (displacement)
%              - 2 : Direction 2 (displacement)
%              - 3 : Direction 3 (displacement)
%              - 4 : Direction 1 (rotation)
%              - 5 : Direction 2 (rotation)
%              - 6 : Direction 3 (rotation)
%
% label   : STING object describing the label of the node whose the 
%           displacement should be plotted on the abscissa axis

%--------------------------------------------------------------------------
% Output arguments
%--------------------------------------------------------------------------
% EV2     : EVOL object describing the reaction curve.

%--------------------------------------------------------------------------
% Example
%--------------------------------------------------------------------------
PB1 = PROBLEM('model',MO1,'mater',MA1,'diric',CLT,'loadt',CHT,'comp_time',0:0.1:1);

SOL = SOLVE(PB1);

SOL = 

  SOLUTION with properties:

            problem: [1x1 PROBLEM]
               mesh: [1x1 MESH]
           topology: [1x1 TOPOLOGY]
       displacement: [1x10 CHPOINT]
           velocity: []
       acceleration: []
       cham_element: [1x10 struct]
           reaction: [1x10 CHPOINT]
    internal_forces: [1x10 CHPOINT]
     eigenfrequency: []
          modeshape: []

EV2 = EVOL.REAC(SOL,CL1,EV1,1,'P3')

EV2 = 

  EVOL with properties:

       absc: [1x11 double]
       ordo: [1x11 double]
    x_label: 'Displacement (m)'
    y_label: 'Reaction (N)'