| 381 | |
| 382 | |
| 383 | int KRAlphaExplicit_TP::formEleResidual(FE_Element *theEle) |
| 384 | { |
| 385 | theEle->zeroResidual(); |
| 386 | |
| 387 | // this does not work because for some elements damping is returned |
| 388 | // with the residual as well as the damping tangent |
| 389 | //theEle->addRtoResidual(alphaR); |
| 390 | //theEle->addD_Force(*Udot, -alphaD); |
| 391 | //theEle->addM_Force(*Udotdot, -alphaM); |
| 392 | |
| 393 | // instead use residual including the inertia terms and then correct |
| 394 | // the mass contribution (only works because alphaR = alphaD) |
| 395 | theEle->addRIncInertiaToResidual(alphaR); |
| 396 | theEle->addM_Force(*Udotdot, alphaR-alphaM); |
| 397 | |
| 398 | return 0; |
| 399 | } |
| 400 | |
| 401 | |
| 402 | int KRAlphaExplicit_TP::formNodUnbalance(DOF_Group *theDof) |
nothing calls this directly
no test coverage detected