%--------------------------------------------------------------------------
% CHAMELEM.CONNEC
%--------------------------------------------------------------------------
% Computation of a connectivity field described by a CHAMELEM.
%--------------------------------------------------------------------------
% Syntax
%--------------------------------------------------------------------------
CHAM2 = CHAMELEM.CONNEC(MO1,lc)
%--------------------------------------------------------------------------
% Description
%--------------------------------------------------------------------------
% CHAMELEM.CONNEC allows to compute a connectivity field used when the
% non local approach is activated.
%--------------------------------------------------------------------------
% Input arguments
%--------------------------------------------------------------------------
% MO1 : MODEL object
% lc : scalar - internal length
%--------------------------------------------------------------------------
% Output arguments
%--------------------------------------------------------------------------
% CHAM2 : CHAMELEM object describing the connectivity field. CHAM2 has 5
% components:
% - coor_pint : coordinates of the integration points in the
% physical domain
% - conn_pint : connectivity between the current integration point
% and the ones which are located inside the averaging
% domain; weights are also stored here
% - detJ_pint : determinant of the Jacobian matrix
% - weight_pint : weigth when computing the non local integral
% - conn_pint_elt: connectivity between the element in which the
% current integration is and the ones which contain
% the integration points which are connected
%--------------------------------------------------------------------------
% Example
%--------------------------------------------------------------------------
CHAM2 = CHAMELEM.CONNEC(MO1,0.5)
CHAM2 =
CHAMELEM with properties:
model: [1x1 MODEL]
cham_val: [1x1 struct]
CHAM2.cham_val.val_elem(1)
ans =
coor_pint: {[0.0211 0.0789] [0.0211 0.0211] [0.0789 0.0789] [0.0789 0.0211]}
conn_pint: {[100x3 double] [89x3 double] [111x3 double] [100x3 double]}
detJ_pint: {[0.0025] [0.0025] [0.0025] [0.0025]}
weight_pint: {[1] [1] [1] [1]}
conn_pint_elt: {[100x4 double] [89x4 double] [111x4 double] [100x4 double]}
CHAM2.cham_val.val_elem(1).conn_pint{1}
ans =
1.0000 1.0000 0.0351
1.0000 2.0000 0.0332
1.0000 3.0000 0.0332
1.0000 4.0000 0.0315
1.0000 5.0000 0.0299
1.0000 6.0000 0.0283
1.0000 7.0000 0.0235
1.0000 8.0000 0.0223
% This table means that the current integration point is connected to the
% integration points 1 to 8, belonging to the cluser 1 and with weigths of
% 0.0351 to 0.0223
CHAM2.cham_val.val_elem(1).conn_pint_elt{1}
ans =
1 1 1 1
1 1 1 2
1 1 1 3
1 1 1 4
1 1 2 1
1 1 2 2
1 1 2 3
% This table is similar to the previous one but is traduced in terms of
% MODEL, sub-domain, element and integration point number.