%--------------------------------------------------------------------------
% ELEMENT.FINT_ELEM
%--------------------------------------------------------------------------
% Computation of the elementary internal forces.
%--------------------------------------------------------------------------
% Syntax
%--------------------------------------------------------------------------
[out, out2] = FINT_ELEM(EL1,MO1,cara,sigf,IR)
[out, out2] = ELEMENT.FINT_ELEM(EL1,MO1,cara,sigf,IR)
%--------------------------------------------------------------------------
% Description
%--------------------------------------------------------------------------
% The FINT_ELEM method of the ELEMENT class is used to compute the
% elementary internal forces for a unique ELEMENT.
%--------------------------------------------------------------------------
% Input arguments
%--------------------------------------------------------------------------
% EL1 : ELEMENT object
% MO1 : MODEL object related to the ELEMENT object EL1
% cara : struct - description of the material parameters related to
% the ELEMENT object EL1.
% sigf : cell - description of the stress at each integration points.
% sigf{i} contains the stress evaluated at the ith integration
% point.
% IR : INT_RULE object describing the integration rule
%--------------------------------------------------------------------------
% Output arguments
%--------------------------------------------------------------------------
% out/out2 : scalar, vector - elementary internal force to be assembled.
%
%--------------------------------------------------------------------------
% Example
%--------------------------------------------------------------------------
% Case of a QUA4 ELEMENT.
cara =
rho: {[0] [0] [0] [0]}
alpha: {1x4 cell}
visq: {[1] [1] [1] [1]}
epai: {[1] [1] [1] [1]}
youn: {1x4 cell}
nu: {[0.2000] [0.2000] [0.2000] [0.2000]}
out = FINT_ELEM(EL1,MO1,cara,sigf,IR)
out =
-0.4893
-24.1518
5.6237
28.0728
-1.7027
20.1457
-3.4317
-24.0667