00001 #include <errno.h>
00002 #include "dampermeability.h"
00003 #include "globalt.h"
00004
00005 dampermeability::dampermeability (void)
00006 {
00007
00008 kinvis = new gfunct (tab,12);
00009
00010 kinvis->tabf->x[0]=0.0;
00011 kinvis->tabf->y[0]=1.787e-6;
00012
00013 kinvis->tabf->x[1]=5.0;
00014 kinvis->tabf->y[1]=1.519e-6;
00015
00016 kinvis->tabf->x[2]=10.0;
00017 kinvis->tabf->y[2]=1.307e-6;
00018
00019 kinvis->tabf->x[3]=20.0;
00020 kinvis->tabf->y[3]=1.004e-6;
00021
00022 kinvis->tabf->x[4]=30.0;
00023 kinvis->tabf->y[4]=0.801e-6;
00024
00025 kinvis->tabf->x[5]=40.0;
00026 kinvis->tabf->y[5]=0.658e-6;
00027
00028 kinvis->tabf->x[6]=50.0;
00029 kinvis->tabf->y[6]=0.553e-6;
00030
00031 kinvis->tabf->x[7]=60.0;
00032 kinvis->tabf->y[7]=0.475e-6;
00033
00034 kinvis->tabf->x[8]=70.0;
00035 kinvis->tabf->y[8]=0.413e-6;
00036
00037 kinvis->tabf->x[9]=80.0;
00038 kinvis->tabf->y[9]=0.365e-6;
00039
00040 kinvis->tabf->x[10]=90.0;
00041 kinvis->tabf->y[10]=0.326e-6;
00042
00043 kinvis->tabf->x[11]=100.0;
00044 kinvis->tabf->y[11]=0.294e-6;
00045
00046
00047
00048
00049
00050
00051
00052
00053
00054
00055
00056
00057
00058
00059
00060
00061
00062
00063
00064
00065
00066
00067
00068
00069
00070
00071
00072
00073
00074
00075
00076
00077
00078
00079
00080
00081
00082
00083
00084
00085
00086
00087
00088
00089
00090
00091
00092
00093
00094
00095
00096
00097
00098
00099
00100
00101
00102
00103
00104
00105
00106
00107
00108
00109
00110
00111
00112
00113
00114
00115
00116
00117
00118
00119
00120
00121
00122
00123
00124
00125
00126
00127
00128
00129
00130
00131
00132 }
00133
00134 dampermeability::~dampermeability (void)
00135 {
00136 delete kinvis;
00137 }
00138
00139
00140
00141
00142
00143
00144
00145 void dampermeability::read (XFILE *in)
00146 {
00147 }
00148
00149
00150
00151
00152
00153
00154
00155 void dampermeability::print (FILE *out)
00156 {
00157 }
00158
00159
00160
00161
00162
00163
00164
00165
00166
00167
00168 void dampermeability::matcond (matrix &d,long ipp)
00169 {
00170 long i,n;
00171 double omega,l,wc,alpha,beta,k,kd,kl,km,t,logk;
00172 double tmp;
00173
00174
00175 alpha = 11.3;
00176 beta = 1.64;
00177
00178
00179 n = d.n;
00180
00181
00182 k = d[0][0];
00183
00184
00185
00186 t = 300.0;
00187
00188
00189 omega = Tm->givenontransq (scal_iso_damage,ipp);
00190
00191 if (omega == 0.0)
00192 return;
00193
00194
00195 l = Tm->givenontransq (proc_zone_length,ipp);
00196
00197 wc = Tm->givenontransq (crack_width,ipp);
00198
00199 if (omega<0.15){
00200
00201 errno = 0;
00202 km = k*(1.0+pow((alpha*omega),beta)+pow((alpha*omega),2.0*beta)/2.0);
00203 check_math_errel(Tm->elip[ipp]);
00204 }
00205
00206
00207
00208
00209
00210
00211
00212
00213
00214
00215
00216
00217
00218
00219
00220
00221
00222
00223
00224
00225
00226
00227
00228
00229
00230
00231
00232
00233
00234
00235
00236
00237
00238
00239
00240
00241 if (0.15 <= omega){
00242
00243
00244 tmp = omega;
00245 if (omega > 0.9)
00246 omega = 0.9;
00247 errno = 0;
00248 kd = k*(1.0+pow((alpha*omega),beta)+pow((alpha*omega),2.0*beta)/2.0);
00249 check_math_errel(Tm->elip[ipp]);
00250 errno = 0;
00251 kl = wc*wc*wc/12.0/l/kinvis->getval (t-273.15);
00252 check_math_errel(Tm->elip[ipp]);
00253
00254 if (kd>kl)
00255 km = kd+kl;
00256 else
00257 {
00258 logk = (1.0-omega)*log (kd) + omega *log (kl);
00259 errno = 0;
00260 if (test_math_err())
00261 km = kd+kl;
00262
00263 errno = 0;
00264 km = exp (logk);
00265 if (test_math_err())
00266 km = kd+kl;
00267 }
00268 omega = tmp;
00269 }
00270 if (omega > 0.9)
00271 {
00272 kd = wc*wc*wc/12.0/l/kinvis->getval (t-273.15);
00273 if (kd > km)
00274 km = kd;
00275 }
00276 for (i=0;i<n;i++){
00277 d[i][i]=km;
00278 }
00279 }
00280
00281
00282
00283
00284
00285
00286
00287
00288
00289
00290
00291
00292
00293
00294 void dampermeability::give_reqntq(long *antq)
00295 {
00296
00297 antq[scal_iso_damage-1] = 1;
00298 }
00299
00300
00301