[27478] | 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) % {{{
|
---|
[28030] | 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:');
|
---|
[27478] | 44 | disp(' C^2 |u|^(1/m) ');
|
---|
| 45 | disp(' tau_b = - ____________________________');
|
---|
[28030] | 46 | disp(' (|u|/u0 + 1)^(1/m) ');
|
---|
[27478] | 47 | disp(' ');
|
---|
| 48 | fielddisplay(self,'C','friction coefficient [SI]');
|
---|
| 49 | fielddisplay(self,'m','m exponent (set to m=3 in original paper)');
|
---|
[28030] | 50 | fielddisplay(self,'u0','velocity controlling plastic limit');
|
---|
[27478] | 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');
|
---|
[28030] | 56 | WriteData(fid,prefix,'class','friction','object',self,'fieldname','C','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1,'yts',yts);
|
---|
[27478] | 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
|
---|