//get residual with inertia terms
const Vector&
NineNodeMixedQuad::getResistingForceIncInertia( )
{
int tang_flag = 0 ; //don't get the tangent
static Vector res(18);
//do tangent and residual here
formResidAndTangent( tang_flag ) ;
//inertia terms
formInertiaTerms( tang_flag ) ;
res = resid;
// add the damping forces if rayleigh damping
if (alphaM != 0.0 || betaK != 0.0 || betaK0 != 0.0 || betaKc != 0.0)
res += this->getRayleighDampingForces();
// subtract external loads
if (load != 0)
res -= *load;
return res;
}
const Matrix &
ZeroLengthContact2D::getInitialStiff(void)
{
//opserr<<this->getTag()<< " ZeroLengthContact2D::getInitialStiff()" <<endln;
int tang_flag = 1 ; //get the tangent
//do tangent and residual here
formResidAndTangent( tang_flag ) ;
return stiff ;
}
const Vector &
ZeroLengthContact3D::getResistingForce()
{
int tang_flag = 0 ; //don't get the tangent
formResidAndTangent( tang_flag ) ;
return resid ;
}
//return stiffness matrix
const Matrix& ShellMITC9::getTangentStiff( )
{
int tang_flag = 1 ; //get the tangent
//do tangent and residual here
formResidAndTangent( tang_flag ) ;
return stiff ;
}
//return stiffness matrix
const Matrix& ConstantPressureVolumeQuad :: getTangentStiff( )
{
int tang_flag = 1 ; //get the tangent
//do tangent and residual here
formResidAndTangent( tang_flag ) ;
return stiff ;
}
const Vector &
ZeroLengthContact3D::getResistingForceIncInertia()
{
// there is no Inertia
int tang_flag = 0 ; //don't get the tangent
formResidAndTangent( tang_flag ) ;
return resid ;
}
const Matrix &
ZeroLengthContact3D::getInitialStiff(void)
{
int tang_flag = 1 ; //get the tangent
//do tangent and residual here
formResidAndTangent( tang_flag ) ;
return stiff ;
}
//return stiffness matrix
const Matrix&
NineNodeMixedQuad::getTangentStiff( )
{
int tang_flag = 1 ; //get the tangent
//do tangent and residual here
formResidAndTangent( tang_flag ) ;
return stiff ;
}
//get residual
const Vector& BbarBrick::getResistingForce( )
{
int tang_flag = 0 ; //don't get the tangent
formResidAndTangent( tang_flag ) ;
if (load != 0)
resid -= *load;
return resid ;
}
//get residual
const Vector& ShellMITC9::getResistingForce( )
{
int tang_flag = 0 ; //don't get the tangent
formResidAndTangent( tang_flag ) ;
// subtract external loads
if (load != 0)
resid -= *load;
return resid ;
}
//get residual
const Vector& ConstantPressureVolumeQuad :: getResistingForce( )
{
int tang_flag = 0 ; //don't get the tangent
formResidAndTangent( tang_flag ) ;
// subtract external loads
if (load != 0)
resid -= *load;
return resid ;
}
const Vector &
ZeroLengthContact2D::getResistingForceIncInertia()
{
// there is no Inertia
//opserr<< this->getTag()<< " ZeroLengthContact2D::getResistingForceIncInertia()" <<endln;
int tang_flag = 0 ; //don't get the tangent
formResidAndTangent( tang_flag ) ;
//opserr<< resid;
return resid ;
}
//get residual with inertia terms
const Vector& BbarBrick::getResistingForceIncInertia( )
{
static Vector res(24);
int tang_flag = 0 ; //don't get the tangent
//do tangent and residual here
formResidAndTangent( tang_flag ) ;
formInertiaTerms( tang_flag ) ;
res = resid;
// add the damping forces if rayleigh damping
if (alphaM != 0.0 || betaK != 0.0 || betaK0 != 0.0 || betaKc != 0.0)
res += this->getRayleighDampingForces();
if (load != 0)
res -= *load;
return res ;
}
