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.9998;
00023
00024 ks = 18.0e12;
00025
00026 kt = 3.6e6;
00027
00028 kw = 2.0e9;
00029
00030
00031
00032 phi0 = 0.5;
00033
00034 kintr = 4.5e-13;
00035
00036 rhos = 2000.0;
00037
00038 rhow = 1000.0;
00039
00040 muw = 1.0e-3;
00041
00042
00043 }
00044
00045 sejtkrmatc::~sejtkrmatc()
00046 {}
00047
00048
00049
00050
00051
00052
00053
00054
00055
00056 void sejtkrmatc::read(XFILE *in)
00057 {
00058
00059 }
00060
00061
00062
00063
00064
00065
00066
00067
00068
00069
00070
00071 double sejtkrmatc::get_kuw(double pw)
00072 {
00073 double kuw;
00074
00075 kuw = -alpha;
00076
00077 return(kuw);
00078
00079 }
00080
00081
00082
00083
00084
00085
00086
00087
00088
00089
00090 double sejtkrmatc::get_kwu(double pw)
00091 {
00092 double kwu;
00093
00094 kwu = 0.0;
00095
00096 return(kwu);
00097
00098 }
00099
00100
00101
00102
00103
00104
00105
00106
00107
00108
00109 double sejtkrmatc::get_kww(double pw)
00110 {
00111 double kww;
00112
00113 kww = kintr/muw;
00114
00115 return(kww);
00116
00117 }
00118
00119
00120
00121
00122
00123
00124
00125
00126
00127
00128 double sejtkrmatc::get_capuw(double pw)
00129 {
00130 double capuw;
00131
00132 capuw = 0.0;
00133
00134 return(capuw);
00135
00136 }
00137
00138
00139
00140
00141
00142
00143
00144
00145
00146
00147 double sejtkrmatc::get_capwu(double pw)
00148 {
00149 double capwu;
00150
00151 capwu = alpha;
00152
00153 return(capwu);
00154
00155 }
00156
00157
00158
00159
00160
00161
00162
00163
00164
00165
00166 double sejtkrmatc::get_capww(double pw)
00167 {
00168 double n,capww;
00169
00170 n = phi;
00171
00172 capww = (alpha-n)/ks+n/kw;
00173
00174 return(capww);
00175
00176 }
00177
00178
00179
00180
00181
00182
00183
00184
00185
00186
00187 double sejtkrmatc::get_fw1(double pw)
00188 {
00189 double fw1;
00190
00191 fw1 = kintr/muw*rhow;
00192
00193 return(fw1);
00194
00195 }
00196
00197
00198
00199
00200
00201
00202
00203
00204
00205
00206
00207 double sejtkrmatc::get_fu1(double pw)
00208 {
00209 double fu1;
00210
00211 fu1 = rhos*(n-1) + n*rhow;
00212
00213 return(fu1);
00214
00215 }