%--------------------------------------------------------------------------
% MAT_ELEM
%--------------------------------------------------------------------------
% Creation of a MAT_ELEM object.
%--------------------------------------------------------------------------
% Syntax
%--------------------------------------------------------------------------
MA1 = MAT_ELEM('STIFF',MO1,IR,D,cara,EL1);
MA1 = MAT_ELEM('MASS',MO1,IR,D,cara,EL1);
MA1 = MAT_ELEM('DIRI',cur_node,cur_node_loc);
MA1 = MAT_ELEM('MADD',cur_node,cur_node_loc,val1);
MA1 = MAT_ELEM('CONF_RELA',cur_node,cur_node_loc,elem_node_cur,val2);
MA1 = MAT_ELEM('POINT_RELA',cur_node,cur_node_loc,val2);
%--------------------------------------------------------------------------
% Description
%--------------------------------------------------------------------------
% MAT_ELEM creates a MAT_ELEM object. This type of object can be used to
% describe an elementary matrix. This matrix can describe different situations,
% according to the needs. This object has the following properties:
% - type : STRING - type of the desired matrix.
% - mat_val : matrix - values of the elementary matrix stored in
% sparse mode
% - support : NODE object - set of NODE supporting the matrix
% - loc_dof : scalar, vector - set of indices describing the local
% type of degrees of freedom
% - num_rela : scalar - number of kinematic relationships prescribed
% - mat_full : matrix - values of the uncondensed elementary matrix
% stored in sparse mode (for FCQ elements)
%--------------------------------------------------------------------------
% Input arguments
%--------------------------------------------------------------------------
% type : STRING - type can be equal to:
% - STIFF : elementary stiffness matrix
% - MASS : elementary mass matrix
% - DIRI : elementary Dirichlet BC matrix (Lagrange
% multuplier technique)
% - MADD : elementary lumped mass matrix
% - CONF_RELA : elementary conformity kinematic relationship
% matrix
% - POINT_RELA : kinematic relationship matrix between two
% points
% MO1 : MODEL object
% IR : INT_RULE object
% D : scalar, matrix - local modulus (Hooke matrix, secant matrix,
% tangent matrix)
% cara : struct object - material properties
% EL1 : ELEMENT object
% cur_node : NODE object - set of nodes
% cur_node_loc : scalar, vector - set of indices to specify the local type
% of the degrees of freedom
% val1 : scalar - value of the lumped mass
% elem_node_cur: NODE object - set of nodes
% val2 : scalar - right side of the kinematic relationship
%--------------------------------------------------------------------------
% Output arguments
%--------------------------------------------------------------------------
% MA1 : MAT_ELEM object
%--------------------------------------------------------------------------
% Example
%--------------------------------------------------------------------------
M_ELEM
M_ELEM =
MAT_ELEM with properties:
type: 'STIFF'
mat_val: [6x6 double]
support: [1x3 NODE]
loc_dof: [1 2]
num_rela: []
mat_full: []
M_ELEM.mat_val
ans =
1.0e+09 *
(1,1) 2.5894
(3,1) -2.5804
(4,1) -2.5000
(5,1) -0.0091
(6,1) 2.5000
(2,2) 6.4736
(3,2) -1.2500
(4,2) -6.4509
(5,2) 1.2500
. .
. .
. .