%--------------------------------------------------------------------------
% 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]