1 | %FRICTIONREGCOULOMB class definition
|
---|
2 | %
|
---|
3 | % Usage:
|
---|
4 | % frictionregcoulomb=frictionregcoulomb();
|
---|
5 |
|
---|
6 | classdef frictionregcoulomb
|
---|
7 | properties (SetAccess=public)
|
---|
8 | C = NaN;
|
---|
9 | u0 = 0.;
|
---|
10 | m = NaN;
|
---|
11 | end
|
---|
12 | methods
|
---|
13 | function self = frictionregcoulomb(varargin) % {{{
|
---|
14 | switch nargin
|
---|
15 | case 0
|
---|
16 | self=setdefaultparameters(self);
|
---|
17 | case 1
|
---|
18 | self=structtoobj(frictionregcoulomb(),varargin{1});
|
---|
19 | otherwise
|
---|
20 | error('constructor not supported');
|
---|
21 | end
|
---|
22 | end % }}}
|
---|
23 | function self = extrude(self,md) % {{{
|
---|
24 | self.C = project3d(md,'vector',self.C,'type','node');
|
---|
25 | self.m = project3d(md,'vector',self.m,'type','element');
|
---|
26 | end % }}}
|
---|
27 | function self = setdefaultparameters(self) % {{{
|
---|
28 |
|
---|
29 | self.u0 = 1000;
|
---|
30 |
|
---|
31 | end % }}}
|
---|
32 | function md = checkconsistency(self,md,solution,analyses) % {{{
|
---|
33 |
|
---|
34 | %Early return
|
---|
35 | if ~ismember('StressbalanceAnalysis',analyses) & ~ismember('ThermalAnalysis',analyses), return; end
|
---|
36 | md = checkfield(md,'fieldname','friction.C','timeseries',1,'NaN',1,'Inf',1,'>=',0.);
|
---|
37 | md = checkfield(md,'fieldname','friction.u0','NaN',1,'Inf',1,'>',0.,'numel',1);
|
---|
38 | md = checkfield(md,'fieldname','friction.m','NaN',1,'Inf',1,'>',0.,'size',[md.mesh.numberofelements,1]);
|
---|
39 | end % }}}
|
---|
40 | function disp(self) % {{{
|
---|
41 | %See Joughin et al. 2019 (equivalent form by Matt Trevers, poster at AGU 2022) https://agupubs.onlinelåibrary.wiley.com/doi/full/10.1029/2019GL082526
|
---|
42 | disp('Regularized Coulomb friction law (Joughin et al., 2019) parameters:');
|
---|
43 | disp(' Regularized Coulomb friction law reads:');
|
---|
44 | disp(' C^2 |u|^(1/m) ');
|
---|
45 | disp(' tau_b = - ____________________________');
|
---|
46 | disp(' (|u|/u0 + 1)^(1/m) ');
|
---|
47 | disp(' ');
|
---|
48 | fielddisplay(self,'C','friction coefficient [SI]');
|
---|
49 | fielddisplay(self,'m','m exponent (set to m=3 in original paper)');
|
---|
50 | fielddisplay(self,'u0','velocity controlling plastic limit');
|
---|
51 | end % }}}
|
---|
52 | function marshall(self,prefix,md,fid) % {{{
|
---|
53 | yts=md.constants.yts;
|
---|
54 |
|
---|
55 | WriteData(fid,prefix,'name','md.friction.law','data',14,'format','Integer');
|
---|
56 | WriteData(fid,prefix,'class','friction','object',self,'fieldname','C','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1,'yts',yts);
|
---|
57 | WriteData(fid,prefix,'class','friction','object',self,'fieldname','u0','format','Double','scale',1/yts);
|
---|
58 | WriteData(fid,prefix,'class','friction','object',self,'fieldname','m','format','DoubleMat','mattype',2);
|
---|
59 | end % }}}
|
---|
60 | end
|
---|
61 | end
|
---|