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 "sejtkrmat.h"
00013 #include "globalt.h"
00014
00015 sejtkrmat::sejtkrmat()
00016 {
00017
00018 emod = 1.3e6;
00019
00020 nu = 0.4;
00021
00022 alpha = 0.0001;
00023
00024 ks = 1.0e12;
00025
00026 kw = 2.0e9;
00027
00028
00029
00030 phi0 = 0.2975;
00031
00032 kintr = 4.5e-13;
00033
00034 rhos = 2000.0;
00035
00036 rhow = 1000.0;
00037
00038 muw = 1.0e-3;
00039
00040
00041 }
00042
00043 sejtkrmat::~sejtkrmat()
00044 {}
00045
00046
00047
00048
00049
00050
00051
00052
00053
00054 void sejtkrmat::read(XFILE *in)
00055 {
00056
00057 }
00058
00059
00060
00061
00062
00063
00064
00065
00066
00067
00068 double sejtkrmat::get_sw(double pw)
00069 {
00070
00071 double sw,pc;
00072
00073 pc = -pw;
00074
00075 sw = 1.0;
00076
00077 return(sw);
00078 }
00079
00080
00081
00082
00083
00084
00085
00086
00087
00088
00089 double sejtkrmat::get_dsw_dpw(double pw)
00090 {
00091
00092 double dsw_dpc,pc;
00093
00094 pc = -pw;
00095
00096 dsw_dpc = 0.0;
00097
00098 return(dsw_dpc);
00099 }
00100
00101
00102
00103
00104
00105
00106
00107
00108
00109
00110 double sejtkrmat::get_krw(double pw)
00111 {
00112 double krw,s;
00113
00114 s = get_sw(pw);
00115 krw = pow(s,10.0);
00116
00117 return(krw);
00118 }
00119
00120
00121
00122
00123
00124
00125
00126
00127
00128 double sejtkrmat::get_phi()
00129 {
00130 double phi;
00131
00132 phi = phi0;
00133
00134 return(phi);
00135 }
00136
00137
00138
00139
00140
00141
00142
00143
00144
00145 double sejtkrmat::get_kintr()
00146 {
00147 return(kintr);
00148 }
00149
00150
00151
00152
00153
00154
00155
00156
00157
00158 double sejtkrmat::get_rhos()
00159 {
00160 return(rhos);
00161 }
00162
00163
00164
00165
00166
00167
00168
00169
00170
00171 double sejtkrmat::get_alpha()
00172 {
00173 return(alpha);
00174 }
00175
00176
00177
00178
00179
00180
00181
00182
00183
00184 double sejtkrmat::get_emod()
00185 {
00186 return(emod);
00187 }
00188
00189
00190
00191
00192
00193
00194
00195
00196 double sejtkrmat::get_nu()
00197 {
00198 return(nu);
00199 }
00200
00201
00202
00203
00204
00205
00206
00207
00208
00209 double sejtkrmat::get_ks()
00210 {
00211 return(ks);
00212 }
00213
00214
00215
00216
00217
00218
00219
00220
00221
00222 double sejtkrmat::get_kw()
00223 {
00224 return(kw);
00225 }
00226
00227
00228
00229
00230
00231
00232
00233
00234
00235 double sejtkrmat::get_muw()
00236 {
00237 return(muw);
00238 }
00239
00240
00241
00242
00243
00244
00245
00246
00247
00248
00249 double sejtkrmat::get_kuw(double pw)
00250 {
00251 double alpha,s,kuw;
00252
00253 alpha = get_alpha();
00254 s = get_sw(pw);
00255
00256 kuw = -alpha*s;
00257
00258 return(kuw);
00259
00260 }
00261
00262
00263
00264
00265
00266
00267
00268
00269
00270
00271 double sejtkrmat::get_kwu(double pw)
00272 {
00273 double kwu;
00274
00275 kwu = 0.0;
00276
00277 return(kwu);
00278
00279 }
00280
00281
00282
00283
00284
00285
00286
00287
00288
00289
00290 double sejtkrmat::get_kww(double pw)
00291 {
00292 double krw,muw,kww;
00293
00294 krw = get_krw(pw);
00295 kintr = get_kintr();
00296 muw = get_muw();
00297
00298 kww = krw*kintr/muw;
00299
00300 return(kww);
00301
00302 }
00303
00304
00305
00306
00307
00308
00309
00310
00311
00312
00313 double sejtkrmat::get_capuw(double pw)
00314 {
00315 double capuw;
00316
00317 capuw = 0.0;
00318
00319 return(capuw);
00320
00321 }
00322
00323
00324
00325
00326
00327
00328
00329
00330
00331
00332 double sejtkrmat::get_capwu(double pw)
00333 {
00334 double alpha,s,capwu;
00335
00336 alpha = get_alpha();
00337 s = get_sw(pw);
00338
00339 capwu = alpha*s;
00340
00341 return(capwu);
00342
00343 }
00344
00345
00346
00347
00348
00349
00350
00351
00352
00353
00354 double sejtkrmat::get_capww(double pw)
00355 {
00356 double alpha,sw,ks,dsw_dpw,kw,n,capww;
00357
00358 alpha = get_alpha();
00359 sw = get_sw(pw);
00360 dsw_dpw = get_dsw_dpw(pw);
00361 n = get_phi();
00362 ks = get_ks();
00363 kw = get_kw();
00364
00365 capww = (alpha-n)/ks*sw*(sw+dsw_dpw/n*pw)+n*sw/kw + dsw_dpw;
00366
00367 return(capww);
00368
00369 }
00370
00371
00372
00373
00374
00375
00376
00377
00378
00379
00380 double sejtkrmat::get_fw1(double pw)
00381 {
00382 double krw,muw,fw1;
00383
00384 krw = get_krw(pw);
00385 kintr = get_kintr();
00386 muw = get_muw();
00387 rhow = 1000.0;
00388
00389 fw1 = krw*kintr/muw*rhow;
00390
00391 return(fw1);
00392
00393 }
00394
00395
00396
00397
00398
00399
00400
00401
00402
00403
00404
00405 double sejtkrmat::get_fu1(double pw)
00406 {
00407 double sw,n,rhos,fu1;
00408
00409 sw = get_sw(pw);
00410 n = get_phi();
00411 rhos = get_rhos();
00412 rhow = 1000.0;
00413
00414
00415 fu1 = rhos*(n-1) + sw*n*rhow;
00416
00417 return(fu1);
00418
00419 }