%--------------------------------------------------------------------------
% 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)'