source: issm/trunk-jpl/src/m/classes/frictionregcoulomb.m@ 28030

Last change on this file since 28030 was 28030, checked in by jdquinn, 15 months ago

CHG: MATLAB -> Python; cleanup

File size: 2.4 KB
Line 
1%FRICTIONREGCOULOMB class definition
2%
3% Usage:
4% frictionregcoulomb=frictionregcoulomb();
5
6classdef 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
61end
Note: See TracBrowser for help on using the repository browser.