| 298 | |
| 299 | |
| 300 | int AlphaOSGeneralized_TP::formEleResidual(FE_Element *theEle) |
| 301 | { |
| 302 | theEle->zeroResidual(); |
| 303 | |
| 304 | // this does not work because for some elements damping is returned |
| 305 | // with the residual as well as the damping tangent |
| 306 | //theEle->addRtoResidual(alphaR); |
| 307 | //theEle->addD_Force(*Udot, -alphaD); |
| 308 | //theEle->addM_Force(*Udotdot, -alphaM); |
| 309 | |
| 310 | // instead use residual including the inertia terms and then correct |
| 311 | // the mass contribution (only works because alphaR = alphaD) |
| 312 | theEle->addRIncInertiaToResidual(alphaR); |
| 313 | theEle->addM_Force(*Udotdot, alphaR-alphaM); |
| 314 | |
| 315 | return 0; |
| 316 | } |
| 317 | |
| 318 | |
| 319 | int AlphaOSGeneralized_TP::formNodUnbalance(DOF_Group *theDof) |
nothing calls this directly
no test coverage detected