ThreeWayMinorLoss Class Reference

#include <hminor.h>

Inheritance diagram for ThreeWayMinorLoss:

Collaboration diagram for ThreeWayMinorLoss:

Public Member Functions
	ThreeWayMinorLoss (unsigned int uL, const DofOwner *pD, HydraulicFluid *hf, const PressureNode *p0, const PressureNode *p1, const PressureNode *p2, doublereal dK12, doublereal dK21, doublereal A1, doublereal A2, flag fOut)
	~ThreeWayMinorLoss (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_pNode0
const PressureNode *	m_pNode1
const PressureNode *	m_pNode2
const PressureNode *	m_pNodeN
doublereal	m_dKappa12
doublereal	m_dKappa21
doublereal	m_Area1
doublereal	m_Area2
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 104 of file hminor.h.

Constructor & Destructor Documentation

ThreeWayMinorLoss::ThreeWayMinorLoss	(	unsigned int	uL,
		const DofOwner *	pD,
		HydraulicFluid *	hf,
		const PressureNode *	p0,
		const PressureNode *	p1,
		const PressureNode *	p2,
		doublereal	dK12,
		doublereal	dK21,
		doublereal	A1,
		doublereal	A2,
		flag	fOut
	)

Definition at line 221 of file hminor.cc.

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

: Elem(uL, fOut),
HydraulicElem(uL, pDO, hf, fOut),
m_pNode0(p0), m_pNode1(p1), m_pNode2(p2), m_pNodeN(0),
m_dKappa12(dK12), m_dKappa21(dK21), m_Area1(A1), m_Area2(A2),
m_Area(0.),
flow(0.),
vel(0.),
m_dKappa(0.)
{
        ASSERT(m_pNode0 != NULL);
        ASSERT(m_pNode0->GetNodeType() == Node::HYDRAULIC);
        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(A1 > std::numeric_limits<doublereal>::epsilon());
        ASSERT(A2 > std::numeric_limits<doublereal>::epsilon());
}

Here is the call graph for this function:

ThreeWayMinorLoss::~ThreeWayMinorLoss

(

void

)

Definition at line 248 of file hminor.cc.

References NO_OP.

{
        NO_OP;
}

Member Function Documentation

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

virtual

Implements Elem.

Definition at line 288 of file hminor.cc.

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

{
        DEBUGCOUT("Entering ThreeWayMinorLoss::AssJac()" << std::endl);

        ASSERT(m_pNodeN != NULL);

        FullSubMatrixHandler& WM = WorkMat.SetFull();
        WM.Resize(2, 2);

        integer iNode0RowIndex = m_pNode0->iGetFirstRowIndex() + 1;
        integer iNodeNRowIndex = m_pNodeN->iGetFirstRowIndex() + 1;
        integer iNode0ColIndex = m_pNode0->iGetFirstColIndex() + 1;
        integer iNodeNColIndex = m_pNodeN->iGetFirstColIndex() + 1;

        WM.PutRowIndex(1, iNode0RowIndex);
        WM.PutRowIndex(2, iNodeNRowIndex);
        WM.PutColIndex(1, iNode0ColIndex);
        WM.PutColIndex(2, iNodeNColIndex);

        doublereal p0 = m_pNode0->dGetX();
        doublereal p = m_pNodeN->dGetX();

        doublereal jumpPres = fabs(p0-p);

        /* 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((p0 + p)/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 & ThreeWayMinorLoss::AssRes ( SubVectorHandler &  WorkVec, doublereal  dCoef, const VectorHandler &  XCurr, const VectorHandler &  XPrimeCurr  )

virtual

Implements Elem.

Definition at line 342 of file hminor.cc.

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

{
        DEBUGCOUT("Entering ThreeWayMinorLoss::AssRes()" << std::endl);

        WorkVec.Resize(2);

        doublereal p0 = m_pNode0->dGetX();
        doublereal p1 = m_pNode1->dGetX();
        doublereal p2 = m_pNode2->dGetX();
        doublereal p;

        m_pNodeN = NULL;

        if (p1 > p2) {
                m_pNodeN = m_pNode1;
                p = p1;
                m_Area = m_Area1;
        } else {
                m_pNodeN = m_pNode2;
                p = p2;
                m_Area = m_Area2;
        }

        doublereal jumpPres = fabs(p0-p);

        if (p0 > p) {
                m_dKappa = m_dKappa12;  /* flusso diretto da 0 a n */
        } else {
                m_dKappa = m_dKappa21;  /* flusso diretto da n a 0 */
        }

        doublereal density = HF->dGetDensity((p0 + p)/2.);
        flow = density*m_Area*sqrt(2./(m_dKappa*density))*copysign(sqrt(jumpPres), p0 - p);
        vel = flow/(density*m_Area);

#ifdef HYDR_DEVEL
        DEBUGCOUT("RES area :           " << m_Area << std::endl);
        DEBUGCOUT("RES flow:            " << flow << std::endl);
        DEBUGCOUT("RES p0:              " << p0 << std::endl);
        DEBUGCOUT("RES p:               " << p << 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 0 al nodo n " << std::endl);
        DEBUGCOUT("RES portata (nodo n):" << flow << std::endl);
        DEBUGCOUT("****************************************************" << std::endl);
#endif

        integer iNode0RowIndex = m_pNode0->iGetFirstRowIndex() + 1;
        integer iNodeNRowIndex = m_pNodeN->iGetFirstRowIndex() + 1;

        WorkVec.PutItem(1, iNode0RowIndex, flow);
        WorkVec.PutItem(2, iNodeNRowIndex, -flow);

        return WorkVec;
}

Here is the call graph for this function:

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

inlinevirtual

Reimplemented from Elem.

Definition at line 156 of file hminor.h.

References m_pNode0, m_pNode1, and m_pNode2.

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

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

virtual

Reimplemented from Elem.

Definition at line 274 of file hminor.cc.

References MBDYN_EXCEPT_ARGS.

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

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

virtual

Implements HydraulicElem.

Definition at line 255 of file hminor.cc.

References HydraulicElem::THREEWAYMINORLOSS.

{
        return HydraulicElem::THREEWAYMINORLOSS;
}

unsigned int ThreeWayMinorLoss::iGetNumDof ( void  ) const

virtual

Reimplemented from Elem.

Definition at line 268 of file hminor.cc.

{
        return 0;
}

void ThreeWayMinorLoss::Output ( OutputHandler &  OH ) const

virtual

Reimplemented from ToBeOutput.

Definition at line 403 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 & ThreeWayMinorLoss::Restart ( std::ostream &  out ) const

virtual

Implements Elem.

Definition at line 262 of file hminor.cc.

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

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

virtual

Implements Elem.

Definition at line 281 of file hminor.cc.

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

Member Data Documentation

doublereal ThreeWayMinorLoss::flow

private

Definition at line 117 of file hminor.h.

Referenced by AssRes(), and Output().

doublereal ThreeWayMinorLoss::m_Area

private

Definition at line 115 of file hminor.h.

Referenced by AssJac(), and AssRes().

doublereal ThreeWayMinorLoss::m_Area1

private

Definition at line 113 of file hminor.h.

Referenced by AssRes().

doublereal ThreeWayMinorLoss::m_Area2

private

Definition at line 114 of file hminor.h.

Referenced by AssRes().

doublereal ThreeWayMinorLoss::m_dKappa

private

Definition at line 119 of file hminor.h.

Referenced by AssJac(), and AssRes().

doublereal ThreeWayMinorLoss::m_dKappa12

private

Definition at line 111 of file hminor.h.

Referenced by AssRes(), and ThreeWayMinorLoss().

doublereal ThreeWayMinorLoss::m_dKappa21

private

Definition at line 112 of file hminor.h.

Referenced by AssRes(), and ThreeWayMinorLoss().

const PressureNode* ThreeWayMinorLoss::m_pNode0

private

Definition at line 106 of file hminor.h.

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

const PressureNode* ThreeWayMinorLoss::m_pNode1

private

Definition at line 107 of file hminor.h.

Referenced by AssRes(), GetConnectedNodes(), and ThreeWayMinorLoss().

const PressureNode* ThreeWayMinorLoss::m_pNode2

private

Definition at line 108 of file hminor.h.

Referenced by AssRes(), GetConnectedNodes(), and ThreeWayMinorLoss().

const PressureNode* ThreeWayMinorLoss::m_pNodeN

private

Definition at line 109 of file hminor.h.

Referenced by AssJac(), and AssRes().

doublereal ThreeWayMinorLoss::vel

private

Definition at line 118 of file hminor.h.

Referenced by AssRes(), and Output().

---

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

• hminor.h
• hminor.cc