87 tilde_f.
Write(out,
", ") <<
", orientation, reference, node, 1, ",
161 size_t l =
STRLENOF(
"constitutiveLaw.");
162 if (strncmp(s,
"constitutiveLaw.", l) == 0) {
209 return tilde_Theta(i);
220 return tilde_dPrime(i - 3);
230 return tilde_Omega(i - 6);
239 return GetF()(i - 9);
266 for (
int iCnt = 1; iCnt <= 6; iCnt++) {
298 for (
int iCnt = 1; iCnt <= 6; iCnt++) {
326 for (
int iCnt = 1; iCnt <= 6; iCnt++) {
327 WorkVec.
PutRowIndex(iCnt, iNodeFirstMomIndex + iCnt);
364 for (
int iCnt = 1; iCnt <= 6; iCnt++) {
365 WorkVec.
PutRowIndex(iCnt, iNodeFirstMomIndex + iCnt);
388 for (
int iCnt = 1; iCnt <= 6; iCnt++) {
389 WorkVec.
PutRowIndex(iCnt, iNodeFirstPosIndex + iCnt);
411 WMB.
Add(1, 1, F_dPrime);
413 WMB.
Add(4, 1, MTmp2);
419 WMA.
Add(1, 4, F_dPrime*MTmp);
420 WMA.
Add(4, 4, MTmp2*MTmp);
495 for (
int iCnt = 1; iCnt <= 6; iCnt++) {
Mat3x3 MultRMRt(const Mat3x3 &m, const Mat3x3 &R)
void PutColIndex(integer iSubCol, integer iCol)
Vec3 Cross(const Vec3 &v) const
void AssMats(FullSubMatrixHandler &WMA, FullSubMatrixHandler &WMB, doublereal dCoef)
virtual void AfterConvergence(const VectorHandler &X, const VectorHandler &XP)
std::ostream & Write(std::ostream &out, const char *sFill=" ") const
Vec3 MultRV(const Vec3 &v, const Mat3x3 &R)
virtual const Mat3x3 & GetRRef(void) const
virtual bool bToBeOutput(void) const
MatrixExpression< TransposedMatrix< MatrixDirectExpr< Matrix< T, N_rows, N_cols > > >, N_cols, N_rows > Transpose(const Matrix< T, N_rows, N_cols > &A)
ConstitutiveLaw< T, Tder > * pGetConstLaw(void) const
virtual void ResizeReset(integer)
const Vec3 & GetVec2(void) const
const MatCross_Manip MatCross
FullSubMatrixHandler & SetFull(void)
virtual const Mat3x3 & GetRCurr(void) const
virtual Node::Type GetNodeType(void) const
void Add(integer iRow, integer iCol, const Vec3 &v)
virtual unsigned int iGetPrivDataIdx(const char *s) const
virtual Hint * ParseHint(DataManager *pDM, const char *s) const
std::ostream & Write(std::ostream &out, const FullMatrixHandler &m, const char *s, const char *s2)
virtual void Sub(integer iRow, const Vec3 &v)
virtual SubVectorHandler & AssRes(SubVectorHandler &WorkVec, doublereal dCoef, const VectorHandler &XCurr, const VectorHandler &XPrimeCurr)
virtual void Output(OutputHandler &OH) const
virtual unsigned int iGetNumPrivData(void) const
virtual Hint * ParseHint(DataManager *pDM, const char *s) const
std::vector< Hint * > Hints
virtual ~ViscousBody(void)
virtual void WorkSpaceDim(integer *piNumRows, integer *piNumCols) const
virtual VariableSubMatrixHandler & InitialAssJac(VariableSubMatrixHandler &WorkMat, const VectorHandler &XCurr)
Vec3 GetVec(unsigned short int i) const
virtual doublereal dGetPrivData(unsigned int i) const
Vec3 VecRot(const Mat3x3 &Phi)
Vec3 MulTV(const Vec3 &v) const
virtual VariableSubMatrixHandler & AssJac(VariableSubMatrixHandler &WorkMat, doublereal dCoef, const VectorHandler &XCurr, const VectorHandler &XPrimeCurr)
virtual void PutRowIndex(integer iSubRow, integer iRow)=0
virtual std::ostream & Restart(std::ostream &out) const
const Vec3 & GetVec1(void) const
const Vec6 Zero6(0., 0., 0., 0., 0., 0.)
virtual integer iGetFirstMomentumIndex(void) const =0
virtual integer iGetFirstPositionIndex(void) const
virtual const Vec3 & GetWCurr(void) const
virtual std::ostream & Restart(std::ostream &out) const
virtual void AfterPredict(VectorHandler &X, VectorHandler &XP)
std::ostream & Joints(void) const
virtual unsigned int iGetNumPrivData(void) const
void AfterConvergence(const T &Eps, const T &EpsPrime=mb_zero< T >())
#define ASSERT(expression)
virtual void ResizeReset(integer, integer)
virtual SubVectorHandler & InitialAssRes(SubVectorHandler &WorkVec, const VectorHandler &XCurr)
virtual doublereal dGetPrivData(unsigned int i) const
Mat3x3 Transpose(void) const
virtual unsigned int iGetPrivDataIdx(const char *s) const
void PutRowIndex(integer iSubRow, integer iRow)
const T & GetF(void) const
virtual const Vec3 & GetVCurr(void) const
void Sub(integer iRow, integer iCol, const Vec3 &v)
ViscousBody(unsigned int uL, const DofOwner *pDO, const ConstitutiveLaw6D *pCL, const StructNode *pN, const Vec3 &tilde_f, const Mat3x3 &tilde_Rh, const OrientationDescription &od, flag fOut)
void AssVec(SubVectorHandler &WorkVec)
virtual void SetValue(DataManager *pDM, VectorHandler &X, VectorHandler &XP, SimulationEntity::Hints *h=0)
virtual bool bInverseDynamics(void) const
void SetValue(DataManager *pDM, VectorHandler &X, VectorHandler &XP, SimulationEntity::Hints *ph=0)
const Tder & GetFDEPrime(void) const
std::ostream & Output(std::ostream &out, const char *sJointName, unsigned int uLabel, const Vec3 &FLocal, const Vec3 &MLocal, const Vec3 &FGlobal, const Vec3 &MGlobal) const
void Update(const T &Eps, const T &EpsPrime=mb_zero< T >())
unsigned int GetLabel(void) const
virtual void InitialWorkSpaceDim(integer *piNumRows, integer *piNumCols) const