00001
00002
00003
00004
00005
00006
00007 #include <stdio.h>
00008 #include <stdlib.h>
00009 #include <math.h>
00010 #include "aliast.h"
00011 #include "onemediumc.h"
00012 #include "coupmatu.h"
00013 #include "coupmatl.h"
00014 #include "globalc.h"
00015 #include "global.h"
00016 #include "globalt.h"
00017 #include "isotrmatc.h"
00018 #include "intpoints.h"
00019 #include "consol_awf1.h"
00020
00021 medc1::medc1()
00022 {
00023
00024
00025
00026
00027 scale_t = 1.0;
00028 scale_u = 1.0;
00029 }
00030 medc1::~medc1()
00031 {}
00032
00033
00034
00035
00036
00037
00038
00039
00040
00041
00042
00043
00044 void medc1::matcond_u (matrix &d,long ri,long ci,long ipp)
00045 {
00046 long i,m;
00047 double k;
00048 m = d.m;
00049
00050 i = Cmu->ip[ipp].idm;
00051
00052 switch (Cp->mednam){
00053 case mech_heat:{
00054
00055 switch (Cmu->ip[ipp].tm){
00056 case isotransmatc:{
00057 matrix sd(d.m,d.n);
00058
00059 k = Cmu->itrmc[i].get_alpha();
00060 k = -1.0*k;
00061
00062 switch (m){
00063 case 1:{
00064 fillm(0.0,sd);
00065
00066 sd[0][0] = k;
00067 break;
00068 }
00069 case 3:{
00070 fillm(0.0,sd);
00071
00072 sd[0][0] = k;
00073 sd[1][0] = k;
00074 sd[2][0] = 0.0;
00075 break;
00076 }
00077 case 6:{
00078 fillm(0.0,sd);
00079
00080 sd[0][0]=k;
00081 sd[1][0]=k;
00082 sd[2][0]=k;
00083 sd[3][0]=0.0;
00084 sd[4][0]=0.0;
00085 sd[5][0]=0.0;
00086 break;
00087 }
00088 default:{
00089 print_err("\n Unknown number of components of stress tensor is required in ", __FILE__, __LINE__, __func__);
00090 }
00091 }
00092 matrix s(d.m,d.m);
00093
00094 Cmu->matstiff (s,ipp);
00095 mxm(s,sd,d);
00096
00097 destrm (s); destrm (sd);
00098
00099 break;
00100 }
00101 default:{
00102 print_err("\n Unknown material type is required in ", __FILE__, __LINE__, __func__);
00103 }
00104 }
00105 break;
00106 }
00107 case mech_water:{
00108 switch (Cmu->ip[ipp].tm){
00109 case sejtkrc:{
00110 double pw;
00111
00112 pw = Cmu->ip[ipp].av[0];
00113 k = Cmu->sejtkrmc[i].get_kuw(pw);
00114 break;
00115 }
00116 case consolawf1c:{
00117 double pw;
00118
00119 pw = Cmu->ip[ipp].av[0];
00120 k = Cmu->consol_awf1c[i].get_kuw(pw);
00121 break;
00122 }
00123 default:{
00124 print_err("\n Unknown material type is required in ", __FILE__, __LINE__, __func__);
00125 }
00126 }
00127
00128 switch (m){
00129 case 1:{
00130 d[0][0] = k;
00131 break;
00132 }
00133 case 3:{
00134 d[0][0] = k;
00135 d[1][0] = k;
00136 d[2][0] = 0.0;
00137 break;
00138 }
00139 case 6:{
00140 d[0][0]=k;
00141 d[1][0]=k;
00142 d[2][0]=k;
00143 d[3][0]=0.0;
00144 d[4][0]=0.0;
00145 d[5][0]=0.0;
00146 break;
00147 }
00148 default:{
00149 print_err("\n Unknown number of components of stress tensor is required in ", __FILE__, __LINE__, __func__);
00150 }
00151 }
00152 break;
00153 }
00154 default:{
00155 print_err("\n Unknown media name is required in ", __FILE__, __LINE__, __func__);
00156 }
00157 }
00158 }
00159
00160
00161
00162
00163
00164
00165
00166
00167
00168
00169
00170
00171
00172
00173 void medc1::matcap_u (matrix &d,long ri,long ci,long ipp)
00174 {
00175 long i,m;
00176 double c;
00177 m = d.m;
00178
00179 i = Cmu->ip[ipp].idm;
00180
00181 switch (Cp->mednam){
00182 case mech_water:
00183 case mech_heat:{
00184
00185 switch (Cmu->ip[ipp].tm){
00186 case isotransmatc:{
00187
00188
00189 c = 0.0;
00190
00191 break;
00192 }
00193 case sejtkrc:{
00194 double pw;
00195
00196 pw = Cmu->ip[ipp].av[0];
00197 c = Cmu->sejtkrmc[i].get_capuw(pw);
00198 break;
00199 }
00200
00201 case consolawf1c:{
00202 double pw;
00203
00204 pw = Cmu->ip[ipp].av[0];
00205 c = Cmu->consol_awf1c[i].get_capuw(pw);
00206 break;
00207 }
00208
00209 default:{
00210 print_err("\n Unknown material type is required in ", __FILE__, __LINE__, __func__);
00211 }
00212 }
00213
00214 switch (m){
00215 case 1:{
00216 fillm(0.0,d);
00217
00218 d[0][0] = c;
00219 break;
00220 }
00221 case 3:{
00222 fillm(0.0,d);
00223
00224 d[0][0] = c;
00225 d[1][0] = c;
00226 d[2][0] = 0.0;
00227 break;
00228 }
00229 case 6:{
00230 fillm(0.0,d);
00231
00232 d[0][0]=c;
00233 d[1][0]=c;
00234 d[2][0]=c;
00235 d[3][0]=0.0;
00236 d[4][0]=0.0;
00237 d[5][0]=0.0;
00238 break;
00239 }
00240 default:{
00241 print_err("\n Unknown number of components of stress tensor is required in ", __FILE__, __LINE__, __func__);
00242 }
00243 }
00244
00245 break;
00246 }
00247 default:{
00248 print_err("\n Unknown media name is required in ", __FILE__, __LINE__, __func__);
00249 }
00250 }
00251 }
00252
00253
00254
00255
00256
00257
00258
00259
00260
00261
00262
00263
00264 void medc1::matcond_l (matrix &d,long ri,long ci,long ipp)
00265 {
00266 long i,n;
00267 double k;
00268 n = d.n;
00269
00270 i = Cml->ip[ipp].idm;
00271
00272 switch (Cp->mednam){
00273 case mech_water:
00274 case mech_heat:{
00275
00276 switch (Cml->ip[ipp].tm){
00277 case isotransmatc:{
00278
00279 k = 0.0;
00280 break;
00281 }
00282 case sejtkrc:{
00283 double pw;
00284
00285 pw = Cmu->ip[ipp].av[0];
00286 k = Cmu->sejtkrmc[i].get_kwu(pw);
00287 break;
00288 }
00289 case consolawf1c:{
00290 double pw;
00291
00292 pw = Cmu->ip[ipp].av[0];
00293 k = Cmu->consol_awf1c[i].get_kwu(pw);
00294 break;
00295 }
00296 default:{
00297 print_err("\n Unknown material type is required in ", __FILE__, __LINE__, __func__);
00298 }
00299 }
00300
00301 switch (n){
00302 case 1:{
00303 fillm(0.0,d);
00304
00305 d[0][0] = k;
00306 break;
00307 }
00308 case 3:{
00309 fillm(0.0,d);
00310
00311 d[0][0] = k;
00312 d[0][1] = k;
00313 d[0][2] = 0.0;
00314 break;
00315 }
00316 case 6:{
00317 fillm(0.0,d);
00318
00319 d[0][0]=k;
00320 d[0][1]=k;
00321 d[0][2]=k;
00322 d[0][3]=0.0;
00323 d[0][4]=0.0;
00324 d[0][5]=0.0;
00325 break;
00326 }
00327 default:{
00328 print_err("\n Unknown number of components of stress tensor is required in ", __FILE__, __LINE__, __func__);
00329 }
00330 }
00331
00332 break;
00333 }
00334 default:{
00335 print_err("\n Unknown media name is required in ", __FILE__, __LINE__, __func__);
00336 }
00337 }
00338 }
00339
00340
00341
00342
00343
00344
00345
00346
00347
00348
00349
00350
00351
00352
00353 void medc1::matcap_l (matrix &d,long ri,long ci,long ipp)
00354 {
00355 long i,n;
00356 double c;
00357 n = d.n;
00358
00359 i = Cml->ip[ipp].idm;
00360
00361 switch (Cp->mednam){
00362 case mech_water:
00363 case mech_heat:{
00364
00365 switch (Cml->ip[ipp].tm){
00366 case isotransmatc:{
00367
00368
00369 c = 0.0;
00370 break;
00371 }
00372 case sejtkrc:{
00373 double pw;
00374
00375 pw = Cmu->ip[ipp].av[0];
00376 c = Cmu->sejtkrmc[i].get_capwu(pw);
00377 break;
00378 }
00379 case consolawf1c:{
00380 double pw;
00381
00382 pw = Cmu->ip[ipp].av[0];
00383 c = Cmu->consol_awf1c[i].get_capwu(pw);
00384 break;
00385 }
00386 default:{
00387 print_err("\n Unknown material type is required in ", __FILE__, __LINE__, __func__);
00388 }
00389 }
00390
00391 switch (n){
00392 case 1:{
00393 fillm(0.0,d);
00394
00395 d[0][0] = c;
00396 break;
00397 }
00398 case 3:{
00399 fillm(0.0,d);
00400
00401 d[0][0] = c;
00402 d[0][1] = c;
00403 d[0][2] = 0.0;
00404 break;
00405 }
00406 case 6:{
00407 fillm(0.0,d);
00408
00409 d[0][0]=c;
00410 d[0][1]=c;
00411 d[0][2]=c;
00412 d[0][3]=0.0;
00413 d[0][4]=0.0;
00414 d[0][5]=0.0;
00415 break;
00416 }
00417 default:{
00418 print_err("\n Unknown number of components of stress tensor is required in ", __FILE__, __LINE__, __func__);
00419 }
00420 }
00421 break;
00422 }
00423 default:{
00424 print_err("\n Unknown media name is required in ", __FILE__, __LINE__, __func__);
00425 }
00426 }
00427 }
00428
00429
00430
00431
00432
00433
00434
00435
00436
00437
00438
00439
00440
00441 void medc1::rhs_u (matrix &d,long ri,long ci,long ipp)
00442 {
00443 long i,m;
00444 double k;
00445 m = d.m;
00446
00447 i = Cmu->ip[ipp].idm;
00448
00449 switch (Cp->mednam){
00450 case mech_water:{
00451 switch (Cmu->ip[ipp].tm){
00452 case consolawf1c:{
00453 Cmu->consol_awf1c[i].rhs_volume(d,ri,ci,ipp);
00454 break;
00455 }
00456 case sejtkrc:{
00457 break;
00458 }
00459 default:{
00460 print_err("\n Unknown material type is required in ", __FILE__, __LINE__, __func__);
00461 }
00462 }
00463 break;
00464 }
00465 case mech_heat:{
00466
00467 switch (Cmu->ip[ipp].tm){
00468 case isotransmatc:{
00469 matrix sd(d.m,d.n);
00470
00471
00472 k = Cmu->itrmc[i].get_alpha();
00473 k = -1.0*k;
00474
00475 switch (m){
00476 case 1:{
00477 fillm(0.0,sd);
00478
00479 sd[0][0] = k;
00480 break;
00481 }
00482 case 3:{
00483 fillm(0.0,sd);
00484
00485 sd[0][0] = k;
00486 sd[1][0] = k;
00487 sd[2][0] = 0.0;
00488 break;
00489 }
00490 case 6:{
00491 fillm(0.0,sd);
00492
00493 sd[0][0]=k;
00494 sd[1][0]=k;
00495 sd[2][0]=k;
00496 sd[3][0]=0.0;
00497 sd[4][0]=0.0;
00498 sd[5][0]=0.0;
00499 break;
00500 }
00501 default:{
00502 print_err("\n Unknown number of components of stress tensor is required in ", __FILE__, __LINE__, __func__);
00503 }
00504 }
00505 matrix s(d.m,d.m);
00506
00507 Cmu->matstiff (s,ipp);
00508 mxm(s,sd,d);
00509
00510 destrm (s); destrm (sd);
00511
00512 break;
00513 }
00514 default:{
00515 print_err("\n Unknown material type is required in ", __FILE__, __LINE__, __func__);
00516 }
00517 }
00518 break;
00519 }
00520 default:{
00521 print_err("\n Unknown media name is required in ", __FILE__, __LINE__, __func__);
00522 }
00523 }
00524 }