00001
00002
00003
00004
00005
00006
00007
00008 #include <stdio.h>
00009 #include <stdlib.h>
00010 #include <math.h>
00011 #include "constrel.h"
00012 #include "sejtkrmatc.h"
00013 #include "globalt.h"
00014
00015 sejtkrmatc::sejtkrmatc()
00016 {
00017
00018 emod = 10.0e6;
00019
00020 nu = 0.4;
00021
00022 alpha = 0.5;
00023
00024 kz = 18.0e12;
00025
00026 ks = 3.6e6;
00027
00028 kk = 2.0e9;
00029
00030 phi0 = 0.5;
00031
00032 k = 1.13e-10;
00033
00034 rhok = 998.0;
00035
00036 g = 9.806;
00037
00038 lambda = 0.0;
00039
00040 c = 4.1667e-8;
00041 }
00042
00043 sejtkrmatc::~sejtkrmatc()
00044 {}
00045
00046
00047
00048
00049
00050
00051
00052
00053
00054 void sejtkrmatc::read(XFILE *in)
00055 {
00056
00057 alpha = 1.0 - ks/kz;
00058 }
00059
00060
00061
00062
00063
00064
00065
00066
00067
00068
00069
00070 void sejtkrmatc::rhs_volume (matrix &d,long ri,long ci,long ipp)
00071 {
00072 long m;
00073 m = d.m;
00074
00075 switch (m){
00076 case 1:{
00077 rhs1d1 (d,ri,ci,ipp);
00078 break;
00079 }
00080 case 2:{
00081 rhs2d1 (d,ri,ci,ipp);
00082 break;
00083 }
00084 case 3:{
00085 rhs3d1 (d,ri,ci,ipp);
00086 break;
00087 }
00088 default:{
00089 fprintf (stderr,"\n unknown number of components of stress tensor is required");
00090 fprintf (stderr,"\n in function (%s, line %d).\n",__FILE__,__LINE__);
00091 }
00092 }
00093 }
00094
00095
00096
00097
00098
00099
00100
00101
00102
00103 void sejtkrmatc::rhs1d1 (matrix &d,long ri,long ci,long ipp)
00104 {
00105 double f,pw;
00106 double g;
00107
00108 pw = Tm->ip[ipp].av[0];
00109
00110 f = get_fw1(pw);
00111 g = Tp->gr1;
00112
00113 fillm(0.0,d);
00114 d[0][0] = f*g;
00115 }
00116
00117
00118
00119
00120
00121
00122
00123
00124 void sejtkrmatc::rhs2d1 (matrix &d,long ri,long ci,long ipp)
00125 {
00126 double f,pw;
00127 double *g;
00128
00129 pw = Tm->ip[ipp].av[0];
00130
00131 g = new double [2];
00132 f = get_fw1(pw);
00133
00134 g[0] = Tp->gr1;
00135 g[1] = Tp->gr2;
00136
00137 fillm(0.0,d);
00138 d[0][0] = f*g[0];
00139 d[1][0] = f*g[1];
00140 }
00141
00142
00143
00144
00145
00146
00147
00148
00149 void sejtkrmatc::rhs3d1 (matrix &d,long ri,long ci,long ipp)
00150 {
00151 double f,pw;
00152 double *g;
00153
00154 pw = Tm->ip[ipp].av[0];
00155
00156 g = new double [3];
00157 f = get_fw1(pw);
00158
00159 g[0] = Tp->gr1;
00160 g[1] = Tp->gr2;
00161 g[2] = Tp->gr3;
00162
00163 fillm(0.0,d);
00164 d[0][0] = f*g[0];
00165 d[1][0] = f*g[1];
00166 d[2][0] = f*g[2];
00167 }
00168
00169
00170
00171
00172
00173
00174
00175
00176
00177
00178
00179 double sejtkrmatc::get_kuw(double pw)
00180 {
00181 double kuw;
00182
00183 kuw = -alpha;
00184
00185 return(kuw);
00186
00187 }
00188
00189
00190
00191
00192
00193
00194
00195
00196
00197
00198 double sejtkrmatc::get_kwu(double pw)
00199 {
00200 double kwu;
00201
00202 kwu = 0.0;
00203
00204 return(kwu);
00205
00206 }
00207
00208
00209
00210
00211
00212
00213
00214
00215
00216
00217 double sejtkrmatc::get_kww(double pw)
00218 {
00219 double kww;
00220
00221 kww = c/ks;
00222
00223 return(kww);
00224
00225 }
00226
00227
00228
00229
00230
00231
00232
00233
00234
00235
00236 double sejtkrmatc::get_capuw(double pw)
00237 {
00238 double capuw;
00239
00240 capuw = 0.0;
00241
00242 return(capuw);
00243
00244 }
00245
00246
00247
00248
00249
00250
00251
00252
00253
00254
00255 double sejtkrmatc::get_capwu(double pw)
00256 {
00257 double capwu;
00258
00259 capwu = alpha;
00260
00261 return(capwu);
00262
00263 }
00264
00265
00266
00267
00268
00269
00270
00271
00272
00273
00274 double sejtkrmatc::get_capww(double pw)
00275 {
00276 double capww;
00277
00278 lambda = phi0/kk + (1-phi0)/kz - (1 - alpha)/kz;
00279
00280 capww = lambda;
00281
00282 return(capww);
00283
00284 }
00285
00286
00287
00288
00289
00290
00291
00292
00293
00294
00295 double sejtkrmatc::get_fw1(double pw)
00296 {
00297 double fw1;
00298
00299 fw1 = c/ks;
00300
00301 return(fw1);
00302
00303 }
00304
00305
00306
00307
00308
00309
00310
00311
00312
00313
00314
00315 double sejtkrmatc::get_fu1(double pw)
00316 {
00317 double fu1;
00318
00319 fu1 = 0.0;
00320
00321 return(fu1);
00322
00323 }