1 | %FRICTIONWEERTMAN class definition
|
---|
2 | %
|
---|
3 | % Usage:
|
---|
4 | % frictionweertman=frictionhydro();
|
---|
5 |
|
---|
6 | classdef frictionhydro
|
---|
7 | properties (SetAccess=public)
|
---|
8 | q = NaN;
|
---|
9 | C = NaN;
|
---|
10 | As = NaN;
|
---|
11 | end
|
---|
12 | methods
|
---|
13 | function obj = frictionhydro(varargin) % {{{
|
---|
14 | switch nargin
|
---|
15 | case 0
|
---|
16 | obj=setdefaultparameters(obj);
|
---|
17 | otherwise
|
---|
18 | error('constructor not supported');
|
---|
19 | end
|
---|
20 | end % }}}
|
---|
21 | function obj = setdefaultparameters(obj) % {{{
|
---|
22 |
|
---|
23 | end % }}}
|
---|
24 | function md = checkconsistency(obj,md,solution,analyses) % {{{
|
---|
25 |
|
---|
26 | %Early return
|
---|
27 | if ~ismember(StressbalanceAnalysisEnum(),analyses) & ~ismember(ThermalAnalysisEnum(),analyses), return; end
|
---|
28 | md = checkfield(md,'fieldname','frictionhydro.q','NaN',1,'size',[md.mesh.numberofelements 1]);
|
---|
29 | md = checkfield(md,'fieldname','frictionhydro.C','NaN',1,'size',[md.mesh.numberofelements 1]);
|
---|
30 | md = checkfield(md,'fieldname','frictionhydro.As','NaN',1,'size',[md.mesh.numberofelements 1]);
|
---|
31 | end % }}}
|
---|
32 | function disp(obj) % {{{
|
---|
33 | disp(sprintf('Effective Pressure based friction law described in Gagliardini 2007')
|
---|
34 | fielddisplay(obj,'q','friction law exponent q>=1');
|
---|
35 | fielddisplay(obj,'C','friction law max value [SI]');
|
---|
36 | fielddisplay(obj,'As','Sliding Parameter without cavitation [m Pa^-n s^-1]');
|
---|
37 | end % }}}
|
---|
38 | function marshall(obj,md,fid) % {{{
|
---|
39 | WriteData(fid,'enum',FrictionLawEnum,'data',3,'format','Integer');
|
---|
40 | WriteData(fid,'class','friction','object',obj,'fieldname','q','format','DoubleMat','mattype',2);
|
---|
41 | WriteData(fid,'class','friction','object',obj,'fieldname','C','format','DoubleMat','mattype',2);
|
---|
42 | WriteData(fid,'class','friction','object',obj,'fieldname','As','format','DoubleMat','mattype',2);
|
---|
43 | end % }}}
|
---|
44 | end
|
---|
45 | end
|
---|