| 406 | |
| 407 | |
| 408 | void Foam::functionObjects::forces::writeForces() |
| 409 | { |
| 410 | Log << type() << " " << name() << " write:" << nl |
| 411 | << " sum of forces:" << nl |
| 412 | << " pressure : " << sum(force_[0]) << nl |
| 413 | << " viscous : " << sum(force_[1]) << nl |
| 414 | << " porous : " << sum(force_[2]) << nl |
| 415 | << " sum of moments:" << nl |
| 416 | << " pressure : " << sum(moment_[0]) << nl |
| 417 | << " viscous : " << sum(moment_[1]) << nl |
| 418 | << " porous : " << sum(moment_[2]) |
| 419 | << endl; |
| 420 | |
| 421 | writeTime(file(MAIN_FILE)); |
| 422 | file(MAIN_FILE) << tab << setw(1) << '(' |
| 423 | << sum(force_[0]) << setw(1) << ' ' |
| 424 | << sum(force_[1]) << setw(1) << ' ' |
| 425 | << sum(force_[2]) << setw(3) << ") (" |
| 426 | << sum(moment_[0]) << setw(1) << ' ' |
| 427 | << sum(moment_[1]) << setw(1) << ' ' |
| 428 | << sum(moment_[2]) << setw(1) << ')' |
| 429 | << endl; |
| 430 | |
| 431 | if (localSystem_) |
| 432 | { |
| 433 | vectorField localForceN(coordSys_.localVector(force_[0])); |
| 434 | vectorField localForceT(coordSys_.localVector(force_[1])); |
| 435 | vectorField localForceP(coordSys_.localVector(force_[2])); |
| 436 | vectorField localMomentN(coordSys_.localVector(moment_[0])); |
| 437 | vectorField localMomentT(coordSys_.localVector(moment_[1])); |
| 438 | vectorField localMomentP(coordSys_.localVector(moment_[2])); |
| 439 | |
| 440 | writeTime(file(MAIN_FILE)); |
| 441 | file(MAIN_FILE) << tab << setw(1) << '(' |
| 442 | << sum(localForceN) << setw(1) << ' ' |
| 443 | << sum(localForceT) << setw(1) << ' ' |
| 444 | << sum(localForceP) << setw(3) << ") (" |
| 445 | << sum(localMomentN) << setw(1) << ' ' |
| 446 | << sum(localMomentT) << setw(1) << ' ' |
| 447 | << sum(localMomentP) << setw(1) << ')' |
| 448 | << endl; |
| 449 | } |
| 450 | } |
| 451 | |
| 452 | |
| 453 | void Foam::functionObjects::forces::writeBins() |