Method formNodUnbalance

SRC/analysis/integrator/PFEMIntegrator.cpp:503–534

}    

int
PFEMIntegrator::formNodUnbalance(DOF_Group *theDof)
{
    if (sensitivityFlag == 0) {  // NO SENSITIVITY ANALYSIS

        this->TransientIntegrator::formNodUnbalance(theDof);

    } else {  // ASSEMBLE ALL TERMS

        theDof->zeroUnbalance();

        // The term -M*(-1/dt*dvn)
        theDof->addM_Force(dVn, c3);

        // The term -dM/dh*acc
        theDof->addM_ForceSensitivity(*Udotdot, -1.0);

        // The term -C*(dvn)
        //theDof->addD_Force(dVn ,-1.0);

	
        // The term -dC/dh*vel
        theDof->addD_ForceSensitivity(*Udot,-1.0);


        // In case of random loads (have already been formed by 'applyLoadSensitivity')
        theDof->addPtoUnbalance();

    }

    return 0;
}    


int PFEMIntegrator::domainChanged()