MBDyn-1.7.3
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups
Public Member Functions | Protected Attributes | List of all members
AbsoluteInternalForce Class Reference

#include <strforce_impl.h>

Inheritance diagram for AbsoluteInternalForce:
Inheritance graph
[legend]
Collaboration diagram for AbsoluteInternalForce:
Collaboration graph
[legend]

Public Member Functions

 AbsoluteInternalForce (unsigned int uL, const StructNode *pN1, const StructNode *pN2, const TplDriveCaller< Vec3 > *pDC, const Vec3 &TmpArm1, const Vec3 &TmpArm2, flag fOut)
 
 ~AbsoluteInternalForce (void)
 
virtual Force::Type GetForceType (void) const
 
virtual std::ostream & Restart (std::ostream &out) const
 
void WorkSpaceDim (integer *piNumRows, integer *piNumCols) const
 
virtual VariableSubMatrixHandler & AssJac (VariableSubMatrixHandler &WorkMat, doublereal dCoef, const VectorHandler &XCurr, const VectorHandler &XPrimeCurr)
 
virtual SubVectorHandler & AssRes (SubVectorHandler &WorkVec, doublereal dCoef, const VectorHandler &XCurr, const VectorHandler &XPrimeCurr)
 
virtual SubVectorHandler & AssRes (SubVectorHandler &WorkVec, const VectorHandler &, const VectorHandler &, const VectorHandler &, InverseDynamics::Order iOrder)
 
virtual void OutputPrepare (OutputHandler &OH)
 
virtual void Output (OutputHandler &OH) const
 
virtual void InitialWorkSpaceDim (integer *piNumRows, integer *piNumCols) const
 
virtual VariableSubMatrixHandler & InitialAssJac (VariableSubMatrixHandler &WorkMat, const VectorHandler &XCurr)
 
virtual SubVectorHandler & InitialAssRes (SubVectorHandler &WorkVec, const VectorHandler &XCurr)
 
- Public Member Functions inherited from Elem
 Elem (unsigned int uL, flag fOut)
 
virtual ~Elem (void)
 
virtual unsigned int iGetNumDof (void) const
 
virtual std::ostream & DescribeDof (std::ostream &out, const char *prefix="", bool bInitial=false) const
 
virtual void DescribeDof (std::vector< std::string > &desc, bool bInitial=false, int i=-1) const
 
virtual std::ostream & DescribeEq (std::ostream &out, const char *prefix="", bool bInitial=false) const
 
virtual void DescribeEq (std::vector< std::string > &desc, bool bInitial=false, int i=-1) const
 
virtual DofOrder::Order GetDofType (unsigned int) const
 
virtual void AssMats (VariableSubMatrixHandler &WorkMatA, VariableSubMatrixHandler &WorkMatB, const VectorHandler &XCurr, const VectorHandler &XPrimeCurr)
 
virtual bool bInverseDynamics (void) const
 
void SetInverseDynamicsFlags (unsigned uIDF)
 
unsigned GetInverseDynamicsFlags (void) const
 
bool bIsErgonomy (void) const
 
bool bIsRightHandSide (void) const
 
virtual VariableSubMatrixHandler & AssJac (VariableSubMatrixHandler &WorkMat, const VectorHandler &XCurr)
 
virtual int GetNumConnectedNodes (void) const
 
- Public Member Functions inherited from WithLabel
 WithLabel (unsigned int uL=0, const std::string &sN="")
 
virtual ~WithLabel (void)
 
void PutLabel (unsigned int uL)
 
void PutName (const std::string &sN)
 
unsigned int GetLabel (void) const
 
const std::string & GetName (void) const
 
- Public Member Functions inherited from SimulationEntity
 SimulationEntity (void)
 
virtual ~SimulationEntity (void)
 
virtual bool bIsValidIndex (unsigned int i) const
 
virtual DofOrder::Order GetEqType (unsigned int i) const
 
virtual void SetValue (DataManager *pDM, VectorHandler &X, VectorHandler &XP, SimulationEntity::Hints *h=0)
 
virtual Hint * ParseHint (DataManager *pDM, const char *s) const
 
virtual void BeforePredict (VectorHandler &, VectorHandler &, VectorHandler &, VectorHandler &) const
 
virtual void AfterPredict (VectorHandler &X, VectorHandler &XP)
 
virtual void Update (const VectorHandler &XCurr, const VectorHandler &XPrimeCurr)
 
virtual void DerivativesUpdate (const VectorHandler &XCurr, const VectorHandler &XPrimeCurr)
 
virtual void Update (const VectorHandler &XCurr, InverseDynamics::Order iOrder)
 
virtual void AfterConvergence (const VectorHandler &X, const VectorHandler &XP)
 
virtual void AfterConvergence (const VectorHandler &X, const VectorHandler &XP, const VectorHandler &XPP)
 
virtual unsigned int iGetNumPrivData (void) const
 
virtual unsigned int iGetPrivDataIdx (const char *s) const
 
virtual doublereal dGetPrivData (unsigned int i) const
 
virtual std::ostream & OutputAppend (std::ostream &out) const
 
virtual void ReadInitialState (MBDynParser &HP)
 
- Public Member Functions inherited from ToBeOutput
 ToBeOutput (flag fOut=fDefaultOut)
 
virtual ~ToBeOutput (void)
 
virtual void Output (OutputHandler &OH, const VectorHandler &X, const VectorHandler &XP) const
 
virtual flag fToBeOutput (void) const
 
virtual bool bToBeOutput (void) const
 
virtual void SetOutputFlag (flag f=flag(1))
 
- Public Member Functions inherited from StructuralInternalForce
 StructuralInternalForce (unsigned int uL, const StructNode *pN1, const StructNode *pN2, const TplDriveCaller< Vec3 > *pDC, flag fOut)
 
virtual ~StructuralInternalForce (void)
 
virtual void GetConnectedNodes (std::vector< const Node * > &connectedNodes) const
 
- Public Member Functions inherited from Force
 Force (unsigned int uL, flag fOut)
 
virtual ~Force (void)
 
virtual Elem::Type GetElemType (void) const
 
virtual unsigned int iGetInitialNumDof (void) const
 
- Public Member Functions inherited from InitialAssemblyElem
 InitialAssemblyElem (unsigned int uL, flag fOut)
 
virtual ~InitialAssemblyElem (void)
 
- Public Member Functions inherited from SubjectToInitialAssembly
 SubjectToInitialAssembly (void)
 
virtual ~SubjectToInitialAssembly (void)
 

Protected Attributes

const Vec3 Arm1
 
const Vec3 Arm2
 
- Protected Attributes inherited from WithLabel
unsigned int uLabel
 
std::string sName
 
- Protected Attributes inherited from ToBeOutput
flag fOutput
 
- Protected Attributes inherited from StructuralInternalForce
TplDriveOwner< Vec3 > f
 
const StructNode * pNode1
 
const StructNode * pNode2
 

Additional Inherited Members

- Public Types inherited from Elem
enum  Type {
  UNKNOWN = -1, AIRPROPERTIES = 0, INDUCEDVELOCITY, AUTOMATICSTRUCTURAL,
  GRAVITY, BODY, JOINT, JOINT_REGULARIZATION,
  BEAM, PLATE, FORCE, INERTIA,
  ELECTRICBULK, ELECTRIC, THERMAL, HYDRAULIC,
  BULK, LOADABLE, DRIVEN, EXTERNAL,
  AEROMODAL, AERODYNAMIC, GENEL, SOCKETSTREAM_OUTPUT,
  RTAI_OUTPUT = SOCKETSTREAM_OUTPUT, LASTELEMTYPE
}
 
- Public Types inherited from SimulationEntity
typedef std::vector< Hint * > Hints
 
- Public Types inherited from ToBeOutput
enum  { OUTPUT = 0x1U, OUTPUT_MASK = 0xFU, OUTPUT_PRIVATE = 0x10U, OUTPUT_PRIVATE_MASK = ~OUTPUT_MASK }
 
- Public Types inherited from Force
enum  Type {
  UNKNOWN = -1, ABSTRACTFORCE = 0, ABSTRACTINTERNALFORCE, ABSOLUTEDISPFORCE,
  ABSOLUTEINTERNALDISPFORCE, ABSOLUTEFORCE, FOLLOWERFORCE, ABSOLUTECOUPLE,
  FOLLOWERCOUPLE, ABSOLUTEINTERNALFORCE, FOLLOWERINTERNALFORCE, ABSOLUTEINTERNALCOUPLE,
  FOLLOWERINTERNALCOUPLE, TOTALFORCE, TOTALINTERNALFORCE, EXTERNALSTRUCTURAL,
  MODALFORCE, EXTERNALMODAL, LASTFORCETYPE
}
 

Detailed Description

Definition at line 491 of file strforce_impl.h.

Constructor & Destructor Documentation

AbsoluteInternalForce::AbsoluteInternalForce ( unsigned int  uL,
const StructNode *  pN1,
const StructNode *  pN2,
const TplDriveCaller< Vec3 > *  pDC,
const Vec3 &  TmpArm1,
const Vec3 &  TmpArm2,
flag  fOut 
)

Definition at line 1498 of file strforce.cc.

References NO_OP.

1503 : Elem(uL, fOut),
1504 StructuralInternalForce(uL, pN1, pN2, pDC, fOut),
1505 Arm1(TmpArm1), Arm2(TmpArm2)
1506 #ifdef USE_NETCDF
1507 ,
1508 Var_A1(0),
1509 Var_A2(0)
1510 #endif // USE_NETCDF
1511 {
1512  NO_OP;
1513 }
NO_OP
#define NO_OP
Definition: myassert.h:74
StructuralInternalForce::StructuralInternalForce
StructuralInternalForce(unsigned int uL, const StructNode *pN1, const StructNode *pN2, const TplDriveCaller< Vec3 > *pDC, flag fOut)
Definition: strforce.cc:1457
Elem::Elem
Elem(unsigned int uL, flag fOut)
Definition: elem.cc:41
AbsoluteInternalForce::~AbsoluteInternalForce ( void  )

Definition at line 1516 of file strforce.cc.

References NO_OP.

1517 {
1518  NO_OP;
1519 }
NO_OP
#define NO_OP
Definition: myassert.h:74

Member Function Documentation

VariableSubMatrixHandler & AbsoluteInternalForce::AssJac ( VariableSubMatrixHandler &  WorkMat,
doublereal  dCoef,
const VectorHandler &  XCurr,
const VectorHandler &  XPrimeCurr 
)
virtual

Implements Elem.

Definition at line 1556 of file strforce.cc.

References FullSubMatrixHandler::Add(), Arm1, Arm2, DEBUGCOUT, StructuralInternalForce::f, TplDriveOwner< T >::Get(), StructNode::GetRRef(), StructDispNode::iGetFirstMomentumIndex(), StructDispNode::iGetFirstPositionIndex(), MatCrossCross, StructuralInternalForce::pNode1, StructuralInternalForce::pNode2, FullSubMatrixHandler::PutColIndex(), FullSubMatrixHandler::PutRowIndex(), FullSubMatrixHandler::ResizeReset(), VariableSubMatrixHandler::SetFull(), and FullSubMatrixHandler::Sub().

1560 {
1561  DEBUGCOUT("Entering AbsoluteInternalForce::AssJac()" << std::endl);
1562 
1563  FullSubMatrixHandler& WM = WorkMat.SetFull();
1564 
1565  /* Dimensiona e resetta la matrice di lavoro */
1566  WM.ResizeReset(6, 6);
1567 
1568  integer iFirstPositionIndex1 = pNode1->iGetFirstPositionIndex() + 3;
1569  integer iFirstMomentumIndex1 = pNode1->iGetFirstMomentumIndex() + 3;
1570 
1571  integer iFirstPositionIndex2 = pNode2->iGetFirstPositionIndex() + 3;
1572  integer iFirstMomentumIndex2 = pNode2->iGetFirstMomentumIndex() + 3;
1573 
1574  for (integer iCnt = 1; iCnt <= 3; iCnt++) {
1575  WM.PutRowIndex(iCnt, iFirstMomentumIndex1 + iCnt);
1576  WM.PutColIndex(iCnt, iFirstPositionIndex1 + iCnt);
1577 
1578  WM.PutRowIndex(3 + iCnt, iFirstMomentumIndex2 + iCnt);
1579  WM.PutColIndex(3 + iCnt, iFirstPositionIndex2 + iCnt);
1580  }
1581 
1582  /* Dati */
1583  Vec3 TmpArm1(pNode1->GetRRef()*Arm1);
1584  Vec3 TmpArm2(pNode2->GetRRef()*Arm2);
1585  Vec3 TmpDir = f.Get()*dCoef;
1586 
1587  /* | F/\ | | F |
1588  * | | Delta_g = | |
1589  * | (d/\F)/\ | | d/\F |
1590  */
1591 
1592  WM.Sub(1, 1, Mat3x3(MatCrossCross, TmpDir, TmpArm1));
1593  WM.Add(4, 4, Mat3x3(MatCrossCross, TmpDir, TmpArm2));
1594 
1595  return WorkMat;
1596 }
StructuralInternalForce::f
TplDriveOwner< Vec3 > f
Definition: strforce_impl.h:461
FullSubMatrixHandler::PutColIndex
void PutColIndex(integer iSubCol, integer iCol)
Definition: submat.h:325
StructNode::GetRRef
virtual const Mat3x3 & GetRRef(void) const
Definition: strnode.h:1006
Vec3
Definition: matvec3.h:98
VariableSubMatrixHandler::SetFull
FullSubMatrixHandler & SetFull(void)
Definition: submat.h:1168
StructuralInternalForce::pNode1
const StructNode * pNode1
Definition: strforce_impl.h:462
FullSubMatrixHandler::Add
void Add(integer iRow, integer iCol, const Vec3 &v)
Definition: submat.cc:209
DEBUGCOUT
#define DEBUGCOUT(msg)
Definition: myassert.h:232
TplDriveOwner::Get
T Get(const doublereal &dVar) const
Definition: tpldrive.h:109
StructDispNode::iGetFirstMomentumIndex
virtual integer iGetFirstMomentumIndex(void) const =0
StructDispNode::iGetFirstPositionIndex
virtual integer iGetFirstPositionIndex(void) const
Definition: strnode.h:452
FullSubMatrixHandler::ResizeReset
virtual void ResizeReset(integer, integer)
Definition: submat.cc:182
MatCrossCross
const MatCrossCross_Manip MatCrossCross
Definition: matvec3.cc:640
FullSubMatrixHandler::PutRowIndex
void PutRowIndex(integer iSubRow, integer iRow)
Definition: submat.h:311
StructuralInternalForce::pNode2
const StructNode * pNode2
Definition: strforce_impl.h:463
FullSubMatrixHandler::Sub
void Sub(integer iRow, integer iCol, const Vec3 &v)
Definition: submat.cc:215
integer
long int integer
Definition: colamd.c:51

Here is the call graph for this function:

SubVectorHandler & AbsoluteInternalForce::AssRes ( SubVectorHandler &  WorkVec,
doublereal  dCoef,
const VectorHandler &  XCurr,
const VectorHandler &  XPrimeCurr 
)
virtual

Implements Elem.

Definition at line 1601 of file strforce.cc.

References VectorHandler::Add(), Arm1, grad::Cross(), DEBUGCOUT, StructuralInternalForce::f, TplDriveOwner< T >::Get(), StructNode::GetRCurr(), StructDispNode::iGetFirstMomentumIndex(), StructuralInternalForce::pNode1, StructuralInternalForce::pNode2, SubVectorHandler::PutRowIndex(), VectorHandler::ResizeReset(), VectorHandler::Sub(), and WorkSpaceDim().

1605 {
1606  DEBUGCOUT("Entering AbsoluteInternalForce::AssRes()" << std::endl);
1607 
1608  integer iNumRows;
1609  integer iNumCols;
1610  WorkSpaceDim(&iNumRows, &iNumCols);
1611  WorkVec.ResizeReset(iNumRows);
1612 
1613  /* Indici delle incognite del nodo */
1614  integer iFirstMomentumIndex1 = pNode1->iGetFirstMomentumIndex();
1615  integer iFirstMomentumIndex2 = pNode2->iGetFirstMomentumIndex();
1616  for (integer iCnt = 1; iCnt <= 6; iCnt++) {
1617  WorkVec.PutRowIndex(iCnt, iFirstMomentumIndex1 + iCnt);
1618  WorkVec.PutRowIndex(6 + iCnt, iFirstMomentumIndex2 + iCnt);
1619  }
1620 
1621  /* Dati */
1622  Vec3 F(f.Get());
1623  Vec3 M1((pNode1->GetRCurr()*Arm1).Cross(F));
1624  Vec3 M2(F.Cross(pNode2->GetRCurr()*Arm1)); /* - x2 /\ F */
1625 
1626  WorkVec.Add(1, F);
1627  WorkVec.Add(4, M1);
1628  WorkVec.Sub(7, F);
1629  WorkVec.Sub(10, M2);
1630 
1631  return WorkVec;
1632 }
StructuralInternalForce::f
TplDriveOwner< Vec3 > f
Definition: strforce_impl.h:461
Vec3
Definition: matvec3.h:98
VectorHandler::ResizeReset
virtual void ResizeReset(integer)
Definition: vh.cc:55
StructNode::GetRCurr
virtual const Mat3x3 & GetRCurr(void) const
Definition: strnode.h:1012
StructuralInternalForce::pNode1
const StructNode * pNode1
Definition: strforce_impl.h:462
VectorHandler::Sub
virtual void Sub(integer iRow, const Vec3 &v)
Definition: vh.cc:78
SubVectorHandler::PutRowIndex
virtual void PutRowIndex(integer iSubRow, integer iRow)=0
AbsoluteInternalForce::WorkSpaceDim
void WorkSpaceDim(integer *piNumRows, integer *piNumCols) const
Definition: strforce.cc:1523
DEBUGCOUT
#define DEBUGCOUT(msg)
Definition: myassert.h:232
TplDriveOwner::Get
T Get(const doublereal &dVar) const
Definition: tpldrive.h:109
StructDispNode::iGetFirstMomentumIndex
virtual integer iGetFirstMomentumIndex(void) const =0
grad::Cross
VectorExpression< VectorCrossExpr< VectorLhsExpr, VectorRhsExpr >, 3 > Cross(const VectorExpression< VectorLhsExpr, 3 > &u, const VectorExpression< VectorRhsExpr, 3 > &v)
Definition: matvec.h:3248
VectorHandler::Add
virtual void Add(integer iRow, const Vec3 &v)
Definition: vh.cc:63
StructuralInternalForce::pNode2
const StructNode * pNode2
Definition: strforce_impl.h:463
integer
long int integer
Definition: colamd.c:51

Here is the call graph for this function:

SubVectorHandler & AbsoluteInternalForce::AssRes ( SubVectorHandler &  WorkVec,
const VectorHandler &  ,
const VectorHandler &  ,
const VectorHandler &  ,
InverseDynamics::Order  iOrder 
)
virtual

Reimplemented from Elem.

Definition at line 1636 of file strforce.cc.

References VectorHandler::Add(), Arm1, grad::Cross(), DEBUGCOUT, StructuralInternalForce::f, TplDriveOwner< T >::Get(), StructNode::GetRCurr(), StructDispNode::iGetFirstPositionIndex(), StructuralInternalForce::pNode1, StructuralInternalForce::pNode2, SubVectorHandler::PutRowIndex(), VectorHandler::ResizeReset(), VectorHandler::Sub(), and WorkSpaceDim().

1641 {
1642  DEBUGCOUT("Entering AbsoluteInternalForce::AssRes()" << std::endl);
1643 
1644  integer iNumRows;
1645  integer iNumCols;
1646  WorkSpaceDim(&iNumRows, &iNumCols);
1647  WorkVec.ResizeReset(iNumRows);
1648 
1649  /* Indici delle incognite del nodo */
1650  integer iFirstMomentumIndex1 = pNode1->iGetFirstPositionIndex();
1651  integer iFirstMomentumIndex2 = pNode2->iGetFirstPositionIndex();
1652  for (integer iCnt = 1; iCnt <= 6; iCnt++) {
1653  WorkVec.PutRowIndex(iCnt, iFirstMomentumIndex1 + iCnt);
1654  WorkVec.PutRowIndex(6 + iCnt, iFirstMomentumIndex2 + iCnt);
1655  }
1656 
1657  /* Dati */
1658  Vec3 F(f.Get());
1659  Vec3 M1((pNode1->GetRCurr()*Arm1).Cross(F));
1660  Vec3 M2(F.Cross(pNode2->GetRCurr()*Arm1)); /* - x2 /\ F */
1661 
1662  WorkVec.Add(1, F);
1663  WorkVec.Add(4, M1);
1664  WorkVec.Sub(7, F);
1665  WorkVec.Sub(10, M2);
1666 
1667  return WorkVec;
1668 }
StructuralInternalForce::f
TplDriveOwner< Vec3 > f
Definition: strforce_impl.h:461
Vec3
Definition: matvec3.h:98
VectorHandler::ResizeReset
virtual void ResizeReset(integer)
Definition: vh.cc:55
StructNode::GetRCurr
virtual const Mat3x3 & GetRCurr(void) const
Definition: strnode.h:1012
StructuralInternalForce::pNode1
const StructNode * pNode1
Definition: strforce_impl.h:462
VectorHandler::Sub
virtual void Sub(integer iRow, const Vec3 &v)
Definition: vh.cc:78
SubVectorHandler::PutRowIndex
virtual void PutRowIndex(integer iSubRow, integer iRow)=0
AbsoluteInternalForce::WorkSpaceDim
void WorkSpaceDim(integer *piNumRows, integer *piNumCols) const
Definition: strforce.cc:1523
DEBUGCOUT
#define DEBUGCOUT(msg)
Definition: myassert.h:232
TplDriveOwner::Get
T Get(const doublereal &dVar) const
Definition: tpldrive.h:109
StructDispNode::iGetFirstPositionIndex
virtual integer iGetFirstPositionIndex(void) const
Definition: strnode.h:452
grad::Cross
VectorExpression< VectorCrossExpr< VectorLhsExpr, VectorRhsExpr >, 3 > Cross(const VectorExpression< VectorLhsExpr, 3 > &u, const VectorExpression< VectorRhsExpr, 3 > &v)
Definition: matvec.h:3248
VectorHandler::Add
virtual void Add(integer iRow, const Vec3 &v)
Definition: vh.cc:63
StructuralInternalForce::pNode2
const StructNode * pNode2
Definition: strforce_impl.h:463
integer
long int integer
Definition: colamd.c:51

Here is the call graph for this function:

virtual Force::Type AbsoluteInternalForce::GetForceType ( void  ) const
inlinevirtual

Implements Force.

Definition at line 513 of file strforce_impl.h.

References Force::ABSOLUTEINTERNALFORCE.

513  {
514  return Force::ABSOLUTEINTERNALFORCE;
515  };
VariableSubMatrixHandler & AbsoluteInternalForce::InitialAssJac ( VariableSubMatrixHandler &  WorkMat,
const VectorHandler &  XCurr 
)
virtual

Reimplemented from Force.

Definition at line 1729 of file strforce.cc.

References FullSubMatrixHandler::Add(), Arm1, Arm2, DEBUGCOUT, StructuralInternalForce::f, TplDriveOwner< T >::Get(), StructNode::GetRRef(), StructNode::GetWRef(), StructDispNode::iGetFirstPositionIndex(), MatCross, MatCrossCross, StructuralInternalForce::pNode1, StructuralInternalForce::pNode2, FullSubMatrixHandler::PutColIndex(), FullSubMatrixHandler::PutRowIndex(), FullSubMatrixHandler::ResizeReset(), VariableSubMatrixHandler::SetFull(), and FullSubMatrixHandler::Sub().

1731 {
1732  DEBUGCOUT("Entering AbsoluteInternalForce::InitialAssJac()"
1733  << std::endl);
1734 
1735  FullSubMatrixHandler& WM = WorkMat.SetFull();
1736 
1737  /* Dimensiona e resetta la matrice di lavoro */
1738  WM.ResizeReset(12, 12);
1739 
1740  integer iFirstPositionIndex1 = pNode1->iGetFirstPositionIndex() + 3;
1741  integer iFirstVelocityIndex1 = iFirstPositionIndex1 + 6;
1742 
1743  integer iFirstPositionIndex2 = pNode2->iGetFirstPositionIndex() + 3;
1744  integer iFirstVelocityIndex2 = iFirstPositionIndex2 + 6;
1745 
1746  for (integer iCnt = 1; iCnt <= 3; iCnt++) {
1747  WM.PutRowIndex(iCnt, iFirstPositionIndex1 + iCnt);
1748  WM.PutRowIndex(3 + iCnt, iFirstVelocityIndex1 + iCnt);
1749 
1750  WM.PutColIndex(iCnt, iFirstPositionIndex1 + iCnt);
1751  WM.PutColIndex(3 + iCnt, iFirstVelocityIndex1 + iCnt);
1752 
1753  WM.PutRowIndex(6 + iCnt, iFirstPositionIndex2 + iCnt);
1754  WM.PutRowIndex(9 + iCnt, iFirstVelocityIndex2 + iCnt);
1755 
1756  WM.PutColIndex(6 + iCnt, iFirstPositionIndex2 + iCnt);
1757  WM.PutColIndex(9 + iCnt, iFirstVelocityIndex2 + iCnt);
1758  }
1759 
1760  /* Dati */
1761  Vec3 TmpArm1(pNode1->GetRRef()*Arm1);
1762  Vec3 TmpArm2(pNode2->GetRRef()*Arm2);
1763  Vec3 TmpDir = f.Get();
1764  Vec3 Omega1(pNode1->GetWRef());
1765  Vec3 Omega2(pNode2->GetWRef());
1766 
1767  /* | F/\ | | F |
1768  * | | Delta_g = | |
1769  * | (d/\F)/\ | | d/\F |
1770  */
1771 
1772  Mat3x3 MTmp(MatCrossCross, TmpDir, TmpArm1);
1773  WM.Sub(1, 1, MTmp);
1774  WM.Sub(4, 1, Mat3x3(MatCrossCross, TmpDir, Omega1)*Mat3x3(MatCross, TmpArm1));
1775  WM.Sub(4, 4, MTmp);
1776 
1777  MTmp = Mat3x3(MatCrossCross, TmpDir, TmpArm2);
1778  WM.Add(7, 7, MTmp);
1779  WM.Add(10, 7, Mat3x3(MatCrossCross, TmpDir, Omega2)*Mat3x3(MatCross, TmpArm2));
1780  WM.Add(10, 10, MTmp);
1781 
1782  return WorkMat;
1783 }
StructuralInternalForce::f
TplDriveOwner< Vec3 > f
Definition: strforce_impl.h:461
FullSubMatrixHandler::PutColIndex
void PutColIndex(integer iSubCol, integer iCol)
Definition: submat.h:325
StructNode::GetRRef
virtual const Mat3x3 & GetRRef(void) const
Definition: strnode.h:1006
Vec3
Definition: matvec3.h:98
MatCross
const MatCross_Manip MatCross
Definition: matvec3.cc:639
VariableSubMatrixHandler::SetFull
FullSubMatrixHandler & SetFull(void)
Definition: submat.h:1168
StructuralInternalForce::pNode1
const StructNode * pNode1
Definition: strforce_impl.h:462
FullSubMatrixHandler::Add
void Add(integer iRow, integer iCol, const Vec3 &v)
Definition: submat.cc:209
StructNode::GetWRef
virtual const Vec3 & GetWRef(void) const
Definition: strnode.h:1024
DEBUGCOUT
#define DEBUGCOUT(msg)
Definition: myassert.h:232
TplDriveOwner::Get
T Get(const doublereal &dVar) const
Definition: tpldrive.h:109
StructDispNode::iGetFirstPositionIndex
virtual integer iGetFirstPositionIndex(void) const
Definition: strnode.h:452
FullSubMatrixHandler::ResizeReset
virtual void ResizeReset(integer, integer)
Definition: submat.cc:182
MatCrossCross
const MatCrossCross_Manip MatCrossCross
Definition: matvec3.cc:640
FullSubMatrixHandler::PutRowIndex
void PutRowIndex(integer iSubRow, integer iRow)
Definition: submat.h:311
StructuralInternalForce::pNode2
const StructNode * pNode2
Definition: strforce_impl.h:463
FullSubMatrixHandler::Sub
void Sub(integer iRow, integer iCol, const Vec3 &v)
Definition: submat.cc:215
integer
long int integer
Definition: colamd.c:51

Here is the call graph for this function:

SubVectorHandler & AbsoluteInternalForce::InitialAssRes ( SubVectorHandler &  WorkVec,
const VectorHandler &  XCurr 
)
virtual

Implements SubjectToInitialAssembly.

Definition at line 1788 of file strforce.cc.

References VectorHandler::Add(), Arm1, Arm2, grad::Cross(), DEBUGCOUT, StructuralInternalForce::f, TplDriveOwner< T >::Get(), StructNode::GetRCurr(), StructNode::GetWCurr(), StructDispNode::iGetFirstPositionIndex(), InitialWorkSpaceDim(), StructuralInternalForce::pNode1, StructuralInternalForce::pNode2, SubVectorHandler::PutRowIndex(), VectorHandler::ResizeReset(), and VectorHandler::Sub().

1790 {
1791  DEBUGCOUT("Entering AbsoluteInternalForce::InitialAssRes()"
1792  << std::endl);
1793 
1794  integer iNumRows;
1795  integer iNumCols;
1796  InitialWorkSpaceDim(&iNumRows, &iNumCols);
1797  WorkVec.ResizeReset(iNumRows);
1798 
1799  /* Indici delle incognite del nodo */
1800  integer iFirstPositionIndex1 = pNode1->iGetFirstPositionIndex();
1801  integer iFirstVelocityIndex1 = iFirstPositionIndex1 + 6;
1802 
1803  integer iFirstPositionIndex2 = pNode2->iGetFirstPositionIndex();
1804  integer iFirstVelocityIndex2 = iFirstPositionIndex2 + 6;
1805 
1806  for (integer iCnt = 1; iCnt <= 6; iCnt++) {
1807  WorkVec.PutRowIndex(iCnt, iFirstPositionIndex1 + iCnt);
1808  WorkVec.PutRowIndex(6 + iCnt, iFirstVelocityIndex1 + iCnt);
1809 
1810  WorkVec.PutRowIndex(12 + iCnt, iFirstPositionIndex2 + iCnt);
1811  WorkVec.PutRowIndex(18 + iCnt, iFirstVelocityIndex2 + iCnt);
1812  }
1813 
1814  /* Dati */
1815  Vec3 TmpDir(f.Get());
1816  Vec3 TmpArm1(pNode1->GetRCurr()*Arm1);
1817  Vec3 TmpArm2(pNode2->GetRCurr()*Arm2);
1818  const Vec3& Omega1(pNode1->GetWCurr());
1819  const Vec3& Omega2(pNode2->GetWCurr());
1820 
1821  WorkVec.Add(1, TmpDir);
1822  WorkVec.Add(4, TmpArm1.Cross(TmpDir));
1823 
1824  /* In 7 non c'e' nulla */
1825  WorkVec.Add(10, (Omega1.Cross(TmpArm1)).Cross(TmpDir));
1826 
1827  WorkVec.Sub(7, TmpDir);
1828  WorkVec.Sub(10, TmpArm2.Cross(TmpDir));
1829 
1830  /* In 7 non c'e' nulla */
1831  WorkVec.Sub(16, (Omega2.Cross(TmpArm2)).Cross(TmpDir));
1832 
1833  return WorkVec;
1834 }
StructuralInternalForce::f
TplDriveOwner< Vec3 > f
Definition: strforce_impl.h:461
Vec3
Definition: matvec3.h:98
VectorHandler::ResizeReset
virtual void ResizeReset(integer)
Definition: vh.cc:55
StructNode::GetRCurr
virtual const Mat3x3 & GetRCurr(void) const
Definition: strnode.h:1012
StructuralInternalForce::pNode1
const StructNode * pNode1
Definition: strforce_impl.h:462
VectorHandler::Sub
virtual void Sub(integer iRow, const Vec3 &v)
Definition: vh.cc:78
SubVectorHandler::PutRowIndex
virtual void PutRowIndex(integer iSubRow, integer iRow)=0
DEBUGCOUT
#define DEBUGCOUT(msg)
Definition: myassert.h:232
TplDriveOwner::Get
T Get(const doublereal &dVar) const
Definition: tpldrive.h:109
StructDispNode::iGetFirstPositionIndex
virtual integer iGetFirstPositionIndex(void) const
Definition: strnode.h:452
StructNode::GetWCurr
virtual const Vec3 & GetWCurr(void) const
Definition: strnode.h:1030
AbsoluteInternalForce::InitialWorkSpaceDim
virtual void InitialWorkSpaceDim(integer *piNumRows, integer *piNumCols) const
Definition: strforce.cc:1533
grad::Cross
VectorExpression< VectorCrossExpr< VectorLhsExpr, VectorRhsExpr >, 3 > Cross(const VectorExpression< VectorLhsExpr, 3 > &u, const VectorExpression< VectorRhsExpr, 3 > &v)
Definition: matvec.h:3248
VectorHandler::Add
virtual void Add(integer iRow, const Vec3 &v)
Definition: vh.cc:63
StructuralInternalForce::pNode2
const StructNode * pNode2
Definition: strforce_impl.h:463
integer
long int integer
Definition: colamd.c:51

Here is the call graph for this function:

void AbsoluteInternalForce::InitialWorkSpaceDim ( integer *  piNumRows,
integer *  piNumCols 
) const
virtual

Implements SubjectToInitialAssembly.

Definition at line 1533 of file strforce.cc.

Referenced by InitialAssRes().

1536 {
1537  *piNumRows = 24;
1538  *piNumCols = 12;
1539 }
void AbsoluteInternalForce::Output ( OutputHandler &  OH) const
virtual

Reimplemented from ToBeOutput.

Definition at line 1698 of file strforce.cc.

References Arm1, Arm2, ToBeOutput::bToBeOutput(), StructuralInternalForce::f, OutputHandler::FORCES, OutputHandler::Forces(), TplDriveOwner< T >::Get(), OutputHandler::GetCurrentStep(), WithLabel::GetLabel(), StructNode::GetRCurr(), StructDispNode::GetXCurr(), Vec3::pGetVec(), StructuralInternalForce::pNode1, StructuralInternalForce::pNode2, OutputHandler::UseNetCDF(), and OutputHandler::UseText().

1699 {
1700  if (bToBeOutput()) {
1701 #ifdef USE_NETCDF
1702  if (OH.UseNetCDF(OutputHandler::FORCES)) {
1703  Var_F->put_rec(f.Get().pGetVec(), OH.GetCurrentStep());
1704  Var_A1->put_rec((pNode1->GetXCurr() + pNode1->GetRCurr()*Arm1).pGetVec(), OH.GetCurrentStep());
1705  Var_A2->put_rec((pNode2->GetXCurr() + pNode2->GetRCurr()*Arm2).pGetVec(), OH.GetCurrentStep());
1706  }
1707 #endif // USE_NETCDF
1708 
1709  if (OH.UseText(OutputHandler::FORCES)) {
1710  Vec3 F(f.Get());
1711  OH.Forces()
1712  << GetLabel()
1713  << " " << pNode1->GetLabel()
1714  << " " << F
1715  << " " << pNode1->GetXCurr() + pNode1->GetRCurr()*Arm1
1716  << " " << pNode2->GetLabel()
1717  << " " << -F
1718  << " " << pNode2->GetXCurr() + pNode2->GetRCurr()*Arm2
1719  << std::endl;
1720  }
1721 
1722  /* TODO: NetCDF */
1723  }
1724 }
StructuralInternalForce::f
TplDriveOwner< Vec3 > f
Definition: strforce_impl.h:461
ToBeOutput::bToBeOutput
virtual bool bToBeOutput(void) const
Definition: output.cc:890
Vec3
Definition: matvec3.h:98
OutputHandler::UseNetCDF
bool UseNetCDF(int out) const
Definition: output.cc:491
StructNode::GetRCurr
virtual const Mat3x3 & GetRCurr(void) const
Definition: strnode.h:1012
StructuralInternalForce::pNode1
const StructNode * pNode1
Definition: strforce_impl.h:462
OutputHandler::GetCurrentStep
long GetCurrentStep(void) const
Definition: output.h:116
TplDriveOwner::Get
T Get(const doublereal &dVar) const
Definition: tpldrive.h:109
StructDispNode::GetXCurr
virtual const Vec3 & GetXCurr(void) const
Definition: strnode.h:310
Vec3::pGetVec
const doublereal * pGetVec(void) const
Definition: matvec3.h:192
StructuralInternalForce::pNode2
const StructNode * pNode2
Definition: strforce_impl.h:463
WithLabel::GetLabel
unsigned int GetLabel(void) const
Definition: withlab.cc:62
OutputHandler::Forces
std::ostream & Forces(void) const
Definition: output.h:450
OutputHandler::UseText
bool UseText(int out) const
Definition: output.cc:446

Here is the call graph for this function:

void AbsoluteInternalForce::OutputPrepare ( OutputHandler &  OH)
virtual

Reimplemented from ToBeOutput.

Definition at line 1671 of file strforce.cc.

References ASSERT, ToBeOutput::bToBeOutput(), OutputHandler::FORCES, WithLabel::GetLabel(), OutputHandler::IsOpen(), OutputHandler::NETCDF, and OutputHandler::UseNetCDF().

1672 {
1673  if (bToBeOutput()) {
1674 #ifdef USE_NETCDF
1675  if (OH.UseNetCDF(OutputHandler::FORCES)) {
1676  ASSERT(OH.IsOpen(OutputHandler::NETCDF));
1677 
1678  std::ostringstream os;
1679  os << "elem.force." << GetLabel();
1680  (void)OH.CreateVar(os.str(), "internal absolute");
1681 
1682  // joint sub-data
1683  os << '.';
1684  Var_F = OH.CreateVar<Vec3>(os.str() + "F", "N",
1685  "global force components (x, y, z)");
1686 
1687  Var_A1 = OH.CreateVar<Vec3>(os.str() + "Arm1", "m",
1688  "node 1 arm in global frame (x, y, z)");
1689 
1690  Var_A2 = OH.CreateVar<Vec3>(os.str() + "Arm2", "m",
1691  "node 2 arm in global frame (x, y, z)");
1692  }
1693 #endif // USE_NETCDF
1694  }
1695 }
ToBeOutput::bToBeOutput
virtual bool bToBeOutput(void) const
Definition: output.cc:890
Vec3
Definition: matvec3.h:98
OutputHandler::UseNetCDF
bool UseNetCDF(int out) const
Definition: output.cc:491
OutputHandler::IsOpen
bool IsOpen(int out) const
Definition: output.cc:395
ASSERT
#define ASSERT(expression)
Definition: colamd.c:977
WithLabel::GetLabel
unsigned int GetLabel(void) const
Definition: withlab.cc:62

Here is the call graph for this function:

std::ostream & AbsoluteInternalForce::Restart ( std::ostream &  out) const
virtual

Implements Elem.

Definition at line 1544 of file strforce.cc.

References Arm1, Arm2, StructuralInternalForce::f, WithLabel::GetLabel(), TplDriveOwner< T >::pGetDriveCaller(), StructuralInternalForce::pNode1, StructuralInternalForce::pNode2, TplDriveCaller< T >::Restart(), Force::Restart(), and Vec3::Write().

1545 {
1546  Force::Restart(out) << ", absolute internal, "
1547  << pNode1->GetLabel() << ", position, reference, node, ",
1548  Arm1.Write(out, ", ") << ", "
1549  << pNode2->GetLabel() << ", position, reference, node, ",
1550  Arm2.Write(out, ", ") << ", ";
1551  return f.pGetDriveCaller()->Restart(out) << ';' << std::endl;
1552 }
StructuralInternalForce::f
TplDriveOwner< Vec3 > f
Definition: strforce_impl.h:461
Vec3::Write
std::ostream & Write(std::ostream &out, const char *sFill=" ") const
Definition: matvec3.cc:738
StructuralInternalForce::pNode1
const StructNode * pNode1
Definition: strforce_impl.h:462
Force::Restart
virtual std::ostream & Restart(std::ostream &out) const
Definition: force.cc:52
TplDriveOwner::pGetDriveCaller
TplDriveCaller< T > * pGetDriveCaller(void) const
Definition: tpldrive.h:105
TplDriveCaller::Restart
virtual std::ostream & Restart(std::ostream &out) const =0
StructuralInternalForce::pNode2
const StructNode * pNode2
Definition: strforce_impl.h:463
WithLabel::GetLabel
unsigned int GetLabel(void) const
Definition: withlab.cc:62

Here is the call graph for this function:

void AbsoluteInternalForce::WorkSpaceDim ( integer *  piNumRows,
integer *  piNumCols 
) const
virtual

Implements Elem.

Definition at line 1523 of file strforce.cc.

Referenced by AssRes().

1526 {
1527  *piNumRows = 12;
1528  *piNumCols = 6;
1529 }

Member Data Documentation

const Vec3 AbsoluteInternalForce::Arm1
protected

Definition at line 494 of file strforce_impl.h.

Referenced by AssJac(), AssRes(), InitialAssJac(), InitialAssRes(), Output(), and Restart().

const Vec3 AbsoluteInternalForce::Arm2
protected

Definition at line 495 of file strforce_impl.h.

Referenced by AssJac(), InitialAssJac(), InitialAssRes(), Output(), and Restart().

The documentation for this class was generated from the following files: