#include <hminor.h>

Inheritance diagram for MinorLoss:

Collaboration diagram for MinorLoss:

Public Member Functions
	MinorLoss (unsigned int uL, const DofOwner pD, HydraulicFluid hf, const PressureNode p1, const PressureNode p2, doublereal dK12, doublereal dK21, doublereal A, flag fOut)

	~MinorLoss (void)

virtual HydraulicElem::Type	GetHydraulicType (void) const

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

virtual unsigned int	iGetNumDof (void) const

virtual DofOrder::Order	GetDofType (unsigned int i) const

virtual void	WorkSpaceDim (integer piNumRows, integer piNumCols) const

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

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

virtual void	Output (OutputHandler &OH) const

virtual void	GetConnectedNodes (std::vector< const Node * > &connectedNodes) const

Public Member Functions inherited from Elem
	Elem (unsigned int uL, flag fOut)

virtual	~Elem (void)

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 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 SubVectorHandler &	AssRes (SubVectorHandler &WorkVec, const VectorHandler &XCurr, const VectorHandler &XPrimeCurr, const VectorHandler &XPrimePrimeCurr, InverseDynamics::Order iOrder=InverseDynamics::INVERSE_DYNAMICS)

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	OutputPrepare (OutputHandler &OH)

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 HydraulicElem
	HydraulicElem (unsigned int uL, const DofOwner pDO, HydraulicFluid hf, flag fOut)

virtual	~HydraulicElem (void)

virtual Elem::Type	GetElemType (void) const

Public Member Functions inherited from ElemWithDofs
	ElemWithDofs (unsigned int uL, const DofOwner *pDO, flag fOut)

virtual	~ElemWithDofs (void)

Public Member Functions inherited from DofOwnerOwner
	DofOwnerOwner (const DofOwner *pDO)

virtual	~DofOwnerOwner ()

virtual const DofOwner *	pGetDofOwner (void) const

virtual integer	iGetFirstIndex (void) const

virtual void	SetInitialValue (VectorHandler &X)

Private Attributes
const PressureNode *	m_pNode1

const PressureNode *	m_pNode2

doublereal	m_dKappa12

doublereal	m_dKappa21

doublereal	m_Area

doublereal	flow

doublereal	vel

doublereal	m_dKappa

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 HydraulicElem
enum	Type { UNKNOWN = -1, MINOR_LOSS = 0, THREEWAYMINORLOSS, CONTROL_VALVE, DYNAMIC_CONTROL_VALVE, PRESSURE_FLOW_CONTROL_VALVE, PRESSURE_VALVE, FLOW_VALVE, ORIFICE, ACCUMULATOR, TANK, PIPE, DYNAMIC_PIPE, HYDRAULIC_ACTUATOR, ACTUATOR, LASTHYDRAULICTYPE }

Protected Attributes inherited from WithLabel
unsigned int	uLabel

std::string	sName

Protected Attributes inherited from ToBeOutput
flag	fOutput

Protected Attributes inherited from HydraulicElem
HydraulicFluid *	HF

Detailed Description

Definition at line 44 of file hminor.h.

Constructor & Destructor Documentation

MinorLoss::MinorLoss	(	unsigned int	uL,
		const DofOwner *	pD,
		HydraulicFluid *	hf,
		const PressureNode *	p1,
		const PressureNode *	p2,
		doublereal	dK12,
		doublereal	dK21,
		doublereal	A,
		flag	fOut
	)

Definition at line 46 of file hminor.cc.

References ASSERT, PressureNode::GetNodeType(), Node::HYDRAULIC, m_dKappa12, m_dKappa21, m_pNode1, and m_pNode2.

 : Elem(uL, fOut),
 HydraulicElem(uL, pDO, hf, fOut),
 m_pNode1(p1), m_pNode2(p2),
 m_dKappa12(dK12), m_dKappa21(dK21), m_Area(A),
 flow(0.),
 vel(0.),
 m_dKappa(0.)
 {
         ASSERT(m_pNode1 != NULL);
         ASSERT(m_pNode1->GetNodeType() == Node::HYDRAULIC);
         ASSERT(m_pNode2 != NULL);
         ASSERT(m_pNode2->GetNodeType() == Node::HYDRAULIC);
         ASSERT(m_dKappa12 >= 0.);
         ASSERT(m_dKappa21 >= 0.);
         ASSERT(A > std::numeric_limits<doublereal>::epsilon());
 }

Here is the call graph for this function:

MinorLoss::~MinorLoss ( void )

Definition at line 68 of file hminor.cc.

References NO_OP.

 {
         NO_OP;
 }

Member Function Documentation

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

virtual

Implements Elem.

Definition at line 108 of file hminor.cc.

References DEBUGCOUT, HydraulicFluid::dGetDensity(), ScalarAlgebraicNode::dGetX(), grad::fabs(), HydraulicElem::HF, Node::iGetFirstColIndex(), Node::iGetFirstRowIndex(), m_Area, m_dKappa, m_pNode1, m_pNode2, FullSubMatrixHandler::PutCoef(), FullSubMatrixHandler::PutColIndex(), FullSubMatrixHandler::PutRowIndex(), FullSubMatrixHandler::Resize(), VariableSubMatrixHandler::SetFull(), and grad::sqrt().

 {
         DEBUGCOUT("Entering MinorLoss::AssJac()" << std::endl);
 
         FullSubMatrixHandler& WM = WorkMat.SetFull();
         WM.Resize(2, 2);
 
         integer iNode1RowIndex = m_pNode1->iGetFirstRowIndex() + 1;
         integer iNode2RowIndex = m_pNode2->iGetFirstRowIndex() + 1;
         integer iNode1ColIndex = m_pNode1->iGetFirstColIndex() + 1;
         integer iNode2ColIndex = m_pNode2->iGetFirstColIndex() + 1;
 
         WM.PutRowIndex(1, iNode1RowIndex);
         WM.PutRowIndex(2, iNode2RowIndex);
         WM.PutColIndex(1, iNode1ColIndex);
         WM.PutColIndex(2, iNode2ColIndex);
 
         doublereal p1 = m_pNode1->dGetX();
         doublereal p2 = m_pNode2->dGetX();
 
         doublereal jumpPres = fabs(p1-p2);
 
         /* evito di dividere per un numero troppo piccolo */
         if (jumpPres < 1.e8*std::numeric_limits<doublereal>::epsilon()) {
                 jumpPres = 1.e8*std::numeric_limits<doublereal>::epsilon();
         }
 
         /*
          * se voglio usare un fluido comprimibile, metto la pressione
          * media nel condotto:
          */
 
         doublereal density = HF->dGetDensity((p1 + p2)/2.);
 
         /* altrimenti lascio la densita' di riferimento
          * doublereal density = HF->dGetDensity();
          */
 
         doublereal Jac = -density*.5*m_Area*sqrt(2./(m_dKappa*density*jumpPres));
 
         WM.PutCoef(1, 1, Jac);
         WM.PutCoef(1, 2, -Jac);
         WM.PutCoef(2, 1, -Jac);
         WM.PutCoef(2, 2, Jac);
 
         return WorkMat;
 }

Here is the call graph for this function:

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

virtual

Implements Elem.

Definition at line 160 of file hminor.cc.

References copysign(), DEBUGCOUT, HydraulicFluid::dGetDensity(), ScalarAlgebraicNode::dGetX(), grad::fabs(), flow, HydraulicElem::HF, Node::iGetFirstRowIndex(), m_Area, m_dKappa, m_dKappa12, m_dKappa21, m_pNode1, m_pNode2, SubVectorHandler::PutItem(), VectorHandler::Resize(), grad::sqrt(), and vel.

 {
         DEBUGCOUT("Entering MinorLoss::AssRes()" << std::endl);
 
         WorkVec.Resize(2);
 
         integer iNode1RowIndex = m_pNode1->iGetFirstRowIndex() + 1;
         integer iNode2RowIndex = m_pNode2->iGetFirstRowIndex() + 1;
 
         doublereal p1 = m_pNode1->dGetX();
         doublereal p2 = m_pNode2->dGetX();
 
         doublereal jumpPres = fabs(p1-p2);
 
         if (p1 > p2) {
                 m_dKappa = m_dKappa12;  /* flusso diretto da 1 a 2 */
         } else {
                 m_dKappa = m_dKappa21;  /* flusso diretto da 2 a 1 */
         }
 
         doublereal density = HF->dGetDensity((p1 + p2)/2.);
         flow = density*m_Area*sqrt(2./(m_dKappa*density))*copysign(sqrt(jumpPres), p1 - p2);
         vel = flow/(density*m_Area);
 
 #ifdef HYDR_DEVEL
         DEBUGCOUT("RES area :           " << m_Area << std::endl);
         DEBUGCOUT("RES flow:            " << flow << std::endl);
         DEBUGCOUT("RES p1:              " << p1 << std::endl);
         DEBUGCOUT("RES p2:              " << p2 << std::endl);
         DEBUGCOUT("RES dKappa:          " << m_dKappa << std::endl);
         DEBUGCOUT("****************************************************" << std::endl);
         DEBUGCOUT("RES velocita':       " << vel << std::endl);
         DEBUGCOUT("    se positiva il fluido va dal nodo 1 al nodo 2 " << std::endl);
         DEBUGCOUT("RES portata (nodo2): " << flow << std::endl);
         DEBUGCOUT("****************************************************" << std::endl);
 #endif
 
         WorkVec.PutItem(1, iNode1RowIndex, flow);
         WorkVec.PutItem(2, iNode2RowIndex, -flow);
 
         return WorkVec;
 }

Here is the call graph for this function:

virtual void MinorLoss::GetConnectedNodes ( std::vector< const Node * > & connectedNodes ) const

inlinevirtual

Reimplemented from Elem.

Definition at line 91 of file hminor.h.

References m_pNode1, and m_pNode2.

                                                                                       {
                 connectedNodes.resize(2);
                 connectedNodes[0] = m_pNode1;
                 connectedNodes[1] = m_pNode2;
         };

DofOrder::Order MinorLoss::GetDofType ( unsigned int i ) const

virtual

Reimplemented from Elem.

Definition at line 94 of file hminor.cc.

References MBDYN_EXCEPT_ARGS.

 {
         silent_cerr("MinorLoss has no dofs!" << std::endl);
         throw ErrGeneric(MBDYN_EXCEPT_ARGS);
 }

HydraulicElem::Type MinorLoss::GetHydraulicType ( void ) const

virtual

Implements HydraulicElem.

Definition at line 75 of file hminor.cc.

References HydraulicElem::MINOR_LOSS.

 {
         return HydraulicElem::MINOR_LOSS;
 }

unsigned int MinorLoss::iGetNumDof ( void ) const

virtual

Reimplemented from Elem.

Definition at line 88 of file hminor.cc.

 {
         return 0;
 }

void MinorLoss::Output ( OutputHandler & OH ) const

virtual

Reimplemented from ToBeOutput.

Definition at line 207 of file hminor.cc.

References ToBeOutput::bToBeOutput(), flow, WithLabel::GetLabel(), OutputHandler::Hydraulic(), and vel.

 {
         if (bToBeOutput()) {
                 std::ostream& out = OH.Hydraulic();
                 out << std::setw(8) << GetLabel()
                         << " " << vel  << " " << flow << std::endl;
         }
 }

Here is the call graph for this function:

std::ostream & MinorLoss::Restart ( std::ostream & out ) const

virtual

Implements Elem.

Definition at line 82 of file hminor.cc.

 {
         return out << "MinorLoss not implemented yet!" << std::endl;
 }

void MinorLoss::WorkSpaceDim	(	integer *	piNumRows,
		integer *	piNumCols
	)		const

virtual

Implements Elem.

Definition at line 101 of file hminor.cc.

 {
         *piNumRows = 2;
         *piNumCols = 2;
 }

Member Data Documentation

doublereal MinorLoss::flow

private

Definition at line 53 of file hminor.h.

Referenced by AssRes(), and Output().

doublereal MinorLoss::m_Area

private

Definition at line 51 of file hminor.h.

Referenced by AssJac(), and AssRes().

doublereal MinorLoss::m_dKappa

private

Definition at line 55 of file hminor.h.

Referenced by AssJac(), and AssRes().

doublereal MinorLoss::m_dKappa12

private

Definition at line 49 of file hminor.h.

Referenced by AssRes(), and MinorLoss().

doublereal MinorLoss::m_dKappa21

private

Definition at line 50 of file hminor.h.

Referenced by AssRes(), and MinorLoss().

const PressureNode* MinorLoss::m_pNode1

private

Definition at line 46 of file hminor.h.

Referenced by AssJac(), AssRes(), GetConnectedNodes(), and MinorLoss().

const PressureNode* MinorLoss::m_pNode2

private

Definition at line 47 of file hminor.h.

Referenced by AssJac(), AssRes(), GetConnectedNodes(), and MinorLoss().

doublereal MinorLoss::vel

private

Definition at line 54 of file hminor.h.

Referenced by AssRes(), and Output().

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

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation