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