| 275 | |
| 276 | |
| 277 | int HHTGeneralizedExplicit_TP::formEleResidual(FE_Element *theEle) |
| 278 | { |
| 279 | theEle->zeroResidual(); |
| 280 | |
| 281 | // this does not work because for some elements damping is returned |
| 282 | // with the residual as well as the damping tangent |
| 283 | //theEle->addRtoResidual(alphaR); |
| 284 | //theEle->addD_Force(*Udot, -alphaD); |
| 285 | //theEle->addM_Force(*Udotdot, -alphaM); |
| 286 | |
| 287 | // instead use residual including the inertia terms and then correct |
| 288 | // the mass contribution (only works because alphaR = alphaD) |
| 289 | theEle->addRIncInertiaToResidual(alphaR); |
| 290 | theEle->addM_Force(*Udotdot, alphaR-alphaM); |
| 291 | |
| 292 | return 0; |
| 293 | } |
| 294 | |
| 295 | |
| 296 | int HHTGeneralizedExplicit_TP::formNodUnbalance(DOF_Group *theDof) |
nothing calls this directly
no test coverage detected