Homec4science
Diffusion contact_internodes (master) Browse

rCONTACTINTERNODES/Code/Westergaard/setNumericalParameters.m
fccc43844265master

Last Change
Raw File

setNumericalParameters.m
View Options

function[Data]=setNumericalParameters(Mesh, Data, varargin)
% setNumericalParameters: Initializes the numerical parameters of the
% solver
% OPTIONAL INPUT:
% Contact:
% iterate: Parameter controlling whether iterations should
% take place (1) or not (0)
% Default: 1
% maxiter: Maximum number of iterations for solving the
% contact problem
% Default: 10
% rescaling: Rescaling parameter for the stiffness matrix
% Default: 1
% c: Number of nearest neighbors for computing the
% radius of support of the radial basis functions
% Default: 1
% radius: Intersection radius used to define the initial
% potential contact interface.
% Default: infinity
% gamma: Value of gamma for preventing quick loss of diagonal
% dominance of the matrices PhiMM
% Default: 0.5
% Gamma: Value of Gamma to avoid very small nonzero entries
% in the matrices PhiNM
% Default: 0.95
% d0: Measure of tolerance for contact detection.
% Default: 0.05
% tol: Tolerance for gap detection
% Default: 1e-1
% average: Parameter controlling whether a simple average (1)
% or a weighted average (0) of the normal vectors
% should be used.
% Default: 1
% Solver:
% name: Name of the solver to be used.
% Default: direct
% PostProcessing:
% quantity: Additional quantities to be computed (for example
% stresses, strains, fluxes,...)
% Default: none
% ampl_fact: Value of amplification factor to apply to
% displacements
% Default: auto
% eval_points: Number of evaluation points used for computing
% stresses or strains
% Default: 3 (P1 finite elements)
% 6 (P2 or higher order finite elements)
if nargin>2
if mod(length(varargin)-1, 2)~=0
error('Missing a name or value')
end
id=lower(varargin{1});
labels_in=varargin(2:2:end);
values_in=varargin(3:2:end);
end
% Set default Parameters
if ~isfield(Data, 'Contact')
% Boolean indicator for contact algorithm iterations
Data.Contact.iterate=true;
% Maximum number of iterations of the contact algorithm
Data.Contact.maxiter=10;
% Rescaling parameter for the INTERNODES matrix
Data.Contact.rescaling=1;
% Number of nearest neighbors
Data.Contact.c=1;
% Intersection radius
Data.Contact.radius=inf;
% gamma value
Data.Contact.gamma=0.5;
% Gamma value
Data.Contact.Gamma=0.95;
% Contact tolerance
Data.Contact.d0=0.05;
% Gap tolerance
Data.Contact.tol=1e-1;
% Average for normal vectors
Data.Contact.average=1;
end
if ~isfield(Data, 'Solver')
% Solver name
Data.Solver.name='direct';
end
if ~isfield(Data, 'PostProcessing')
% Additional quantities to be computed
Data.PostProcessing.quantity='none';
% Amplification factor
Data.PostProcessing.amplification='auto';
% Number of evaluation points per element for computing stresses/strains
eval_points={3,6};
if Mesh.order==1
Data.PostProcessing.eval_points=eval_points{1};
else
Data.PostProcessing.eval_points=eval_points{2};
end
end
% Override the defaults
if exist('id', 'var')
switch id
case 'contact'
%% Contact algorithm specific parameters
for i=1:length(labels_in)
arg=lower(labels_in{i});
val=lower(values_in{i});
switch arg
case {'iterate', 'average'}
if isa(val, 'double')
if val<0 || val>1
error([arg ' parameter must be a boolean'])
else
Data.Contact.(arg)=logical(val);
end
elseif isa(val, 'logical')
Data.Contact.(arg)=val;
else
error([arg ' parameter must be a boolean'])
end
case {'maxiter', 'c', 'radius', 'd0'}
if val<=0
error(['Parameter ' arg ' must be strictly greater than zero'])
else
Data.Contact.(arg)=val;
end
case 'rescaling'
if val==0
error(['Parameter ' arg ' must be different from zero'])
else
Data.Contact.(arg)=val;
end
case {'gamma', 'Gamma', 'tol'}
if val<=0 || val>=1
error([arg ' must be strictly greater than 0 and strictly smaller than 1'])
else
Data.Contact.(arg)=val;
end
otherwise
error('Unrecognized argument')
end
end
case 'solver'
%% Solver parameters
names={'direct'};
for i=1:length(labels_in)
arg=lower(labels_in{i});
val=values_in{i};
switch arg
case 'name'
if any(ismember(names, val))
Data.Solver.(arg)=val;
else
error('Unrecognized solver name')
end
otherwise
error('Unrecognized argument')
end
end
case 'postprocessing'
%% Post-processing parameters
quantities={'none', 'flux', 'error', 'stress', 'strain'};
eval_points={3,6};
for i=1:length(labels_in)
arg=lower(labels_in{i});
val=values_in{i};
switch arg
case 'quantity'
if any(ismember(quantities, val))
Data.PostProcessing.quantity=val;
else
error('Unrecognized quantity')
end
case 'eval_points'
if any(ismember(eval_points, val))
Data.PostProcessing.eval_points=val;
else
error('Invalid number of evaluation points')
end
case 'ampl_fact'
if isa(val, 'char')
if ~strcmp(val, 'auto')
error('Invalid amplification factor')
else
Data.PostProcessing.amplification=val;
end
elseif isa(val, 'double')
if val <= 0
error('Invalid amplification factor')
else
Data.PostProcessing.amplification=val;
end
else
error('Invalid amplification factor')
end
otherwise
error('Unrecognized argument')
end
end
otherwise
error('Unrecognized argument')
end
end
end
c4science · Help