420
VI=VI+VIorg;
421
}
422
423
if (/*(LAYERPROB)||*/(UI<UIorg))
424
{
425
VI = VI - 2.0*VIorg;
426
VIorg=VIorg/2.0;
427
VI = VI + VIorg;
428
//
LAYERPROB = 0;
429
}
430
if (VI<0.0)
431
{
432
VI = V(0.5/mmpdTomps);
433
}
434
435
return 1;
436
437
}
438
439
int Trigidice::slope(){
440
if(GRADPROB){
441
orgresolution=orgresolution/2;
442
LCMAX=LCMAX*2;
443
return 1;
444
}
445
else{
446
return 0;
447
}
448
}
449
450
451
void Trigidice::quasi(){
452
453
VB = VBorg;
454
GRADTB = GRADTBorg;
455
VI = V(0.5/mmpdTomps);
456
VIorg = V(10/mmpdTomps);
457
resolution=orgresolution;
458
do{
459
do{
460
once();
461
}while(slope());
462
} while (pressdiff());
463
}
464
465
double Trigidice::heave(){
466
467
quasi();
468
return (VI * gammaIA*lamda/(eta*nu)* mmpdTomps);
469
}
470
471
double Trigidice::pressure(){
472
once();
473
return ((UI*gammaIA/lamda)/kP2P);
26