Flow_valve Class Reference

#include <valve.h>

Inheritance diagram for Flow_valve:

Collaboration diagram for Flow_valve:

Public Member Functions
	Flow_valve (unsigned int uL, const DofOwner *pD, HydraulicFluid *hf, const PressureNode *p1, const PressureNode *p2, const PressureNode *p3, doublereal A_dia, doublereal mv, doublereal A_pipe, doublereal A_max, doublereal K, doublereal F0, doublereal w, doublereal s_mx, doublereal cs, doublereal cv, doublereal ca, flag fOut)
	~Flow_valve (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	SetValue (DataManager *pDM, VectorHandler &X, VectorHandler &XP, SimulationEntity::Hints *ph=0)
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 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 *	pNode1
const PressureNode *	pNode2
const PressureNode *	pNode3
doublereal	area_diaf
doublereal	mass
doublereal	area_pipe
doublereal	area_max
doublereal	Kappa
doublereal	force0
doublereal	width
doublereal	s_max
doublereal	c_spost
doublereal	c_vel
doublereal	c_acc
doublereal	c1
doublereal	c2
doublereal	c3
doublereal	c4
doublereal	cf1
doublereal	cf2
doublereal	cf3
doublereal	cf4
doublereal	h
doublereal	s
doublereal	sp
doublereal	v
doublereal	vp
doublereal	flow1
doublereal	flow2
doublereal	flow3

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 516 of file valve.h.

Constructor & Destructor Documentation

Flow_valve::Flow_valve	(	unsigned int	uL,
		const DofOwner *	pD,
		HydraulicFluid *	hf,
		const PressureNode *	p1,
		const PressureNode *	p2,
		const PressureNode *	p3,
		doublereal	A_dia,
		doublereal	mv,
		doublereal	A_pipe,
		doublereal	A_max,
		doublereal	K,
		doublereal	F0,
		doublereal	w,
		doublereal	s_mx,
		doublereal	cs,
		doublereal	cv,
		doublereal	ca,
		flag	fOut
	)

Definition at line 1840 of file valve.cc.

References ASSERT, PressureNode::GetNodeType(), h, Node::HYDRAULIC, pNode1, pNode2, pNode3, and s_max.

: Elem(uL, fOut),
HydraulicElem(uL, pDO, hf, fOut),
pNode1(p1), pNode2(p2), pNode3(p3),
area_diaf(A_dia), mass(mv), area_pipe(A_pipe), area_max(A_max), 
Kappa(K), force0(F0), width(w), s_max(s_mx), c_spost(cs), c_vel(cv), c_acc(ca) 
{
   ASSERT(pNode1 != NULL);
   ASSERT(pNode1->GetNodeType() == Node::HYDRAULIC);
   ASSERT(pNode2 != NULL);
   ASSERT(pNode2->GetNodeType() == Node::HYDRAULIC);
   ASSERT(pNode3 != NULL);
   ASSERT(pNode3->GetNodeType() == Node::HYDRAULIC);
   ASSERT(A_dia > std::numeric_limits<doublereal>::epsilon());
   ASSERT(mv > std::numeric_limits<doublereal>::epsilon()); 
   ASSERT(A_pipe > std::numeric_limits<doublereal>::epsilon());
   ASSERT(A_max > std::numeric_limits<doublereal>::epsilon());
   ASSERT(K > std::numeric_limits<doublereal>::epsilon() );
   ASSERT(F0 >= 0.);
   ASSERT(w > std::numeric_limits<doublereal>::epsilon()); 
   ASSERT(s_max >= 0.);
   
   h = .02; /* coefficiente di perdita di carico concentrata tra i nodi 1 e 2 (smorza il moto della valvola) */
}

Here is the call graph for this function:

Flow_valve::~Flow_valve

(

void

)

Definition at line 1875 of file valve.cc.

References NO_OP.

{
   NO_OP;
}

Member Function Documentation

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

virtual

Implements Elem.

Definition at line 1916 of file valve.cc.

References area_diaf, area_max, area_pipe, c1, c2, c3, cf1, cf2, cf3, copysign(), DEBUGCOUT, HydraulicFluid::dGetDensity(), ScalarAlgebraicNode::dGetX(), grad::fabs(), h, HydraulicElem::HF, Node::iGetFirstColIndex(), DofOwnerOwner::iGetFirstIndex(), Node::iGetFirstRowIndex(), Kappa, mass, pNode1, pNode2, pNode3, grad::pow(), FullSubMatrixHandler::PutCoef(), FullSubMatrixHandler::PutColIndex(), FullSubMatrixHandler::PutRowIndex(), FullSubMatrixHandler::Resize(), s, VariableSubMatrixHandler::SetFull(), grad::sqrt(), v, and width.

{
   DEBUGCOUT("Entering  Flow_valve::AssJac()" << std::endl);
   
   FullSubMatrixHandler& WM = WorkMat.SetFull();
   WM.Resize(5, 5);
   
   integer iNode1RowIndex = pNode1->iGetFirstRowIndex()+1;
   integer iNode2RowIndex = pNode2->iGetFirstRowIndex()+1;
   integer iNode3RowIndex = pNode3->iGetFirstRowIndex()+1;
   integer iNode1ColIndex = pNode1->iGetFirstColIndex()+1;
   integer iNode2ColIndex = pNode2->iGetFirstColIndex()+1;
   integer iNode3ColIndex = pNode3->iGetFirstColIndex()+1;
   integer iFirstIndex = iGetFirstIndex();
   
   WM.PutRowIndex(1, iNode1RowIndex);
   WM.PutRowIndex(2, iNode2RowIndex);
   WM.PutRowIndex(3, iNode3RowIndex);
   WM.PutColIndex(1, iNode1ColIndex);
   WM.PutColIndex(2, iNode2ColIndex);
   WM.PutColIndex(3, iNode3ColIndex);

   WM.PutRowIndex(4, iFirstIndex+1);
   WM.PutColIndex(4, iFirstIndex+1);
   WM.PutRowIndex(5, iFirstIndex+2);
   WM.PutColIndex(5, iFirstIndex+2);

   doublereal p1 = pNode1->dGetX();
   doublereal p2 = pNode2->dGetX();
   doublereal p3 = pNode3->dGetX();
   s = XCurr(iFirstIndex+1);  /* spostamento */
   v = XCurr(iFirstIndex+2);  /* velocita' */
   
   doublereal Cd = .6;        /* verificare */
   doublereal jumpPres12 = fabs(p1-p2);
   doublereal jumpPres23 = fabs(p2-p3);
   doublereal jumpPres13 = fabs(p1-p3);
   doublereal density = HF->dGetDensity();
    
   /* evito di dividere per un numero troppo piccolo */
   if (jumpPres12 < 1.e8*std::numeric_limits<doublereal>::epsilon()) {
      jumpPres12 = 1.e8*std::numeric_limits<doublereal>::epsilon();
   }
   if (jumpPres23 < 1.e8*std::numeric_limits<doublereal>::epsilon()) {
      jumpPres23 = 1.e8*std::numeric_limits<doublereal>::epsilon();
   }
   if (jumpPres13 < 1.e8*std::numeric_limits<doublereal>::epsilon()) {
      jumpPres13 = 1.e8*std::numeric_limits<doublereal>::epsilon();
   }
 
   doublereal Jac11 = 0.;
   doublereal Jac12 = 0.;
   doublereal Jac13 = 0.;
   doublereal Jac14 = 0.;
   doublereal Jac15 = 0.;
   doublereal Jac21 = 0.;
   doublereal Jac22 = 0.;
   doublereal Jac23 = 0.;
   doublereal Jac24 = 0.;
   doublereal Jac25 = 0.;
   doublereal Jac31 = 0.;
   doublereal Jac32 = 0.;
   doublereal Jac33 = 0.;
   doublereal Jac34 = 0.;
   doublereal Jac35 = 0.;
   doublereal Jac41 = 0.;
   doublereal Jac42 = 0.;
   doublereal Jac43 = 0.;
   doublereal Jac44 = 0.;
   doublereal Jac45 = 0.;
   doublereal Jac51 = 0.;
   doublereal Jac52 = 0.;
   doublereal Jac53 = 0.;
   doublereal Jac54 = 0.;
   doublereal Jac55 = 0.;

#ifdef HYDR_DEVEL
   DEBUGCOUT("Valore di p1: " << p1 << std::endl);
   DEBUGCOUT("Valore di p2: " << p2 << std::endl);
   DEBUGCOUT("Valore di p3: " << p3 << std::endl);
   DEBUGCOUT("Valore di s: " << s << std::endl);
#endif
   
   doublereal s12 = sqrt(2.*jumpPres12/density);
   doublereal s13 = sqrt(2.*jumpPres13/density);
   
   Jac11 = -width*s*Cd/s12-area_diaf*Cd/s13;
   Jac12 = width*s*Cd/s12;
   Jac13 = area_diaf*Cd/s13;
   Jac14 = -density*width*Cd*copysign(s12, p1-p2)*dCoef;
   Jac15 = 0.;
   
   Jac21 = width*s*Cd/s12;
   Jac22 = -width*s*Cd/s12;
   Jac23 = 0.;
   Jac24 = density*width*Cd*copysign(s12, p1-p2)*dCoef;
   Jac25 = 0.;
   
   Jac31 = area_diaf*Cd/s13;
   Jac32 = 0.;
   Jac33 = -area_diaf*Cd/s13;
   Jac34 = 0.;
   Jac35 = 0.;
   
   Jac41 = area_max;
   Jac42 = 0.;
   Jac43 =-area_max;
   Jac44 = -Kappa*dCoef-c1*dCoef-cf1*dCoef;
   
   
   doublereal Jac45old1;
   doublereal Jac45old2;
   Jac45old1 = -mass-c2*dCoef-c3-cf2*dCoef-cf3
     -density*area_max*pow(area_max/(area_diaf*Cd), 2.)*fabs(v)*dCoef;
   Jac45 = -mass-c2*dCoef-c3-cf2*dCoef-cf3
     -h*density*pow(area_max/area_pipe, 2.)*fabs(v)*dCoef;
   Jac45old2 = -mass-c2*dCoef-c3-cf2*dCoef-cf3-44.4*fabs(v)*dCoef;
   
   Jac51 = 0.;
   Jac52 = 0.;
   Jac53 = 0.;
   Jac54 = -1.;
   Jac55 = dCoef;
   
#ifdef HYDR_DEVEL
   DEBUGCOUT("Jac smorzatore " 
             << density*area_max*pow(area_max/(area_diaf*Cd), 2) << std::endl);
   DEBUGCOUT("Jac smorzatore " << h*density << std::endl);
   DEBUGCOUT("Jac smorzatore Jac45old1 " << Jac45old1 << std::endl);
   DEBUGCOUT("Jac smorzatore Jac45old2 " << Jac45old2 << std::endl);
   DEBUGCOUT("Jac smorzatore Jac45     " << Jac45 << std::endl);
   
   DEBUGCOUT("Jac11: " << Jac11 << std::endl);
   DEBUGCOUT("Jac12: " << Jac12 << std::endl);
   DEBUGCOUT("Jac13: " << Jac13 << std::endl);
   DEBUGCOUT("Jac14: " << Jac14 << std::endl);
   DEBUGCOUT("Jac15: " << Jac15 << std::endl);
   DEBUGCOUT("Jac21: " << Jac21 << std::endl);
   DEBUGCOUT("Jac22: " << Jac22 << std::endl);
   DEBUGCOUT("Jac23: " << Jac23 << std::endl);
   DEBUGCOUT("Jac24: " << Jac24 << std::endl);
   DEBUGCOUT("Jac25: " << Jac25 << std::endl);
   DEBUGCOUT("Jac31: " << Jac31 << std::endl);
   DEBUGCOUT("Jac32: " << Jac32 << std::endl);
   DEBUGCOUT("Jac33: " << Jac33 << std::endl);
   DEBUGCOUT("Jac34: " << Jac34 << std::endl);
   DEBUGCOUT("Jac35: " << Jac35 << std::endl);
   DEBUGCOUT("Jac41: " << Jac41 << std::endl);
   DEBUGCOUT("Jac42: " << Jac42 << std::endl);
   DEBUGCOUT("Jac43: " << Jac43 << std::endl);
   DEBUGCOUT("Jac44: " << Jac44 << std::endl);
   DEBUGCOUT("Jac45: " << Jac45 << std::endl);
   DEBUGCOUT("Jac51: " << Jac51 << std::endl);
   DEBUGCOUT("Jac52: " << Jac52 << std::endl);
   DEBUGCOUT("Jac53: " << Jac53 << std::endl);
   DEBUGCOUT("Jac54: " << Jac54 << std::endl);
   DEBUGCOUT("Jac55: " << Jac55 << std::endl);
#endif
   
   WM.PutCoef(1, 1, Jac11);
   WM.PutCoef(1, 2, Jac12);
   WM.PutCoef(1, 3, Jac13);
   WM.PutCoef(1, 4, Jac14);
   WM.PutCoef(1, 5, Jac15);
   WM.PutCoef(2, 1, Jac21);
   WM.PutCoef(2, 2, Jac22);
   WM.PutCoef(2, 3, Jac23);
   WM.PutCoef(2, 4, Jac24);
   WM.PutCoef(2, 5, Jac25);
   WM.PutCoef(3, 1, Jac31);
   WM.PutCoef(3, 2, Jac32);
   WM.PutCoef(3, 3, Jac33);
   WM.PutCoef(3, 4, Jac34);
   WM.PutCoef(3, 5, Jac35);
   WM.PutCoef(4, 1, Jac41);
   WM.PutCoef(4, 2, Jac42);
   WM.PutCoef(4, 3, Jac43);
   WM.PutCoef(4, 4, Jac44);
   WM.PutCoef(4, 5, Jac45);
   WM.PutCoef(5, 1, Jac51);
   WM.PutCoef(5, 2, Jac52);
   WM.PutCoef(5, 3, Jac53);
   WM.PutCoef(5, 4, Jac54);
   WM.PutCoef(5, 5, Jac55);
   
   return WorkMat;
}

Here is the call graph for this function:

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

virtual

Implements Elem.

Definition at line 2108 of file valve.cc.

References area_diaf, area_max, area_pipe, c1, c2, c3, c4, c_acc, c_spost, c_vel, cf1, cf2, cf3, cf4, copysign(), DEBUGCOUT, HydraulicFluid::dGetDensity(), ScalarAlgebraicNode::dGetX(), grad::fabs(), flow1, flow2, flow3, force0, h, HydraulicElem::HF, DofOwnerOwner::iGetFirstIndex(), Node::iGetFirstRowIndex(), Kappa, mass, pNode1, pNode2, pNode3, grad::pow(), SubVectorHandler::PutItem(), VectorHandler::Resize(), s, s_max, sp, grad::sqrt(), v, vp, and width.

{
   DEBUGCOUT("Entering Flow_valve::AssRes()" << std::endl);
   
   WorkVec.Resize(5);
   
   integer iNode1RowIndex = pNode1->iGetFirstRowIndex()+1;
   integer iNode2RowIndex = pNode2->iGetFirstRowIndex()+1;
   integer iNode3RowIndex = pNode3->iGetFirstRowIndex()+1;
   
   doublereal p1 = pNode1->dGetX();
   doublereal p2 = pNode2->dGetX();
   doublereal p3 = pNode3->dGetX();
   integer iFirstIndex = iGetFirstIndex();

   s = XCurr(iFirstIndex+1);        /* spostamento */
   v = XCurr(iFirstIndex+2);        /* velocita' */
   sp = XPrimeCurr(iFirstIndex+1);  /* velocita' */
   vp = XPrimeCurr(iFirstIndex+2);  /* accelerazione */
  
   doublereal Cd = .6;
   doublereal jumpPres12 = fabs(p1-p2);
   doublereal jumpPres23 = fabs(p2-p3);
   doublereal jumpPres13 = fabs(p1-p3);
   
   doublereal Res_1 = 0.;
   doublereal Res_2 = 0.;
   doublereal Res_3 = 0.;  
   doublereal Res_4 = 0.;
   doublereal Res_5 = 0.;
   doublereal density = HF->dGetDensity();
   
   doublereal x0 = ((p1-p3)*area_max-force0)/Kappa;
   
   if (s <= 0.) { 
      c1 = c_spost;

      if (sp < 0.) {
         c2 = c_vel; /* ho v negativa devo smorzare */
      } else {
         c2 = 0.;
      }

      if (vp < 0.) {
         c3 = c_acc; /* ho acc negativa devo smorzare */
      } else {
         c3 = 0.;
      }

      c4 = 0.; /* (force0-p1*area_max); */
   } else {
      c1 = 0.;
      c2 = 0.;
      c3 = 0.;
      c4 = 0.;
   }
   
   if (s >= s_max) { 
      cf1 = c_spost;
      
      if (sp > 0.) {
         cf2 = c_vel; /* ho v positiva devo smorzare */
      } else {
         cf2 = 0.;
      }

      if (vp > 0.) {
         cf3 = c_acc; /* ho acc positiva devo smorzare */
      } else {
         cf3 = 0.;
      }

      cf4 = 0.; /* (force0-p1*area_max); */
   } else {
      cf1 = 0.;
      cf2 = 0.;
      cf3 = 0.;
      cf4 = 0.;
   }

   doublereal s12 = sqrt(2.*jumpPres12/density);
   doublereal s13 = sqrt(2.*jumpPres13/density);
   
   Res_1 = density*width*s*Cd*copysign(s12, p1-p2)+
     density*area_diaf*Cd*copysign(s13, p1-p3);
   Res_2 = -density*width*s*Cd*copysign(s12, p1-p2);
   Res_3 = -density*area_diaf*Cd*copysign(s13, p1-p3);

   doublereal Res_4old;
   Res_4old = mass*vp-area_max*p1+force0+Kappa*s+area_max*p3+c1*s+c2*sp+c3*vp
     +c4+cf1*(s-s_max)+cf2*sp+cf3*vp+cf4
     +copysign(.5*density*area_max*pow(sp*area_max/(area_diaf*Cd), 2), sp);
   
   Res_4 = mass*vp-area_max*p1+force0+Kappa*s+area_max*p3+c1*s+c2*sp+c3*vp
     +c4+cf1*(s-s_max)+cf2*sp+cf3*vp+cf4
     +copysign(.5*h*density*pow(sp*area_max/area_pipe, 2), sp);
   
   Res_5 =sp-v;
   
   flow1=-Res_1;   /* per l'output */
   flow2=-Res_2;   /* per l'output */
   flow3=-Res_3;   /* per l'output */
   
#ifdef HYDR_DEVEL
   DEBUGCOUT("Res_4:         " <<  Res_4 << std::endl);
   DEBUGCOUT("Res_4old:      " <<  Res_4old << std::endl);
   DEBUGCOUT("smorazatore:   " 
             << copysign(.5*density*area_max*pow(sp*area_max/(area_diaf*Cd), 2), sp) << std::endl);
   DEBUGCOUT("smorzatoreold: " 
             << copysign(.5*h*density*pow(sp*area_max/area_pipe, 2), sp) << std::endl);
   
   DEBUGCOUT("width:         " << width << std::endl);
   DEBUGCOUT("Cd:            " << Cd << std::endl);
   DEBUGCOUT("jumpPres12:    " << jumpPres12  << std::endl);
   DEBUGCOUT("jumpPres23:    " << jumpPres23  << std::endl);
   DEBUGCOUT("jumpPres13:    " << jumpPres13  << std::endl);
   DEBUGCOUT("density:       " << density << std::endl);
   DEBUGCOUT("p1:            " << p1 << std::endl);
   DEBUGCOUT("p2:            " << p2 << std::endl);
   DEBUGCOUT("p3:            " << p3 << std::endl);
   DEBUGCOUT("s:             " << s << std::endl);
   DEBUGCOUT("sp:            " << sp << std::endl);
   DEBUGCOUT("v:             " << v << std::endl);
   DEBUGCOUT("vp:            " << vp << std::endl);
   DEBUGCOUT("area_max:      " << area_max << std::endl);
   DEBUGCOUT("Kappa:         " << Kappa << std::endl);
   DEBUGCOUT("force0:        " << force0 << std::endl);
   DEBUGCOUT("mass:          " << mass << std::endl);
   DEBUGCOUT("s_max:         " << s_max << std::endl);
   DEBUGCOUT("x0:            " << x0  << std::endl);
   DEBUGCOUT("c1:            " << c1 << std::endl);
   DEBUGCOUT("c2:            " << c2 << std::endl);
   DEBUGCOUT("c3:            " << c3 << std::endl);
   DEBUGCOUT("c4:            " << c4 << std::endl);
   DEBUGCOUT("cf1:           " << cf1 << std::endl);
   DEBUGCOUT("cf2:           " << cf2 << std::endl);
   DEBUGCOUT("cf3:           " << cf3 << std::endl);
   DEBUGCOUT("cf4:           " << cf4 << std::endl);
   DEBUGCOUT("PORTATE AI VARI NODI (positive se entranti)" << std::endl);
   DEBUGCOUT("-Res_1 (portata nodo1): " << -Res_1 << std::endl); 
   DEBUGCOUT("-Res_2 (portata nodo2): " << -Res_2 << std::endl); 
   DEBUGCOUT("-Res_3:(portata nodo3): " << -Res_3 << std::endl); 
   DEBUGCOUT("-Res_4:                 " << -Res_4 << std::endl); 
   DEBUGCOUT("-Res_5:                 " << -Res_5 << std::endl); 
#endif
   
   WorkVec.PutItem(1, iNode1RowIndex, Res_1);
   WorkVec.PutItem(2, iNode2RowIndex, Res_2);         
   WorkVec.PutItem(3, iNode3RowIndex, Res_3);
   WorkVec.PutItem(4, iFirstIndex+1, Res_4);    
   WorkVec.PutItem(5, iFirstIndex+2, Res_5);    
     
   return WorkVec;
}

Here is the call graph for this function:

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

inlinevirtual

Reimplemented from Elem.

Definition at line 590 of file valve.h.

References pNode1, pNode2, and pNode3.

                                                                                 {
     connectedNodes.resize(3);
     connectedNodes[0] = pNode1;
     connectedNodes[1] = pNode2;
     connectedNodes[2] = pNode3;
   };

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

virtual

Reimplemented from Elem.

Definition at line 1901 of file valve.cc.

References ASSERT, and DofOrder::DIFFERENTIAL.

{
   ASSERT(i >= 0 && i <= 1);
   return DofOrder::DIFFERENTIAL;
}

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

virtual

Implements HydraulicElem.

Definition at line 1882 of file valve.cc.

References HydraulicElem::FLOW_VALVE.

{
   return HydraulicElem:: FLOW_VALVE;
}

unsigned int Flow_valve::iGetNumDof ( void  ) const

virtual

Reimplemented from Elem.

Definition at line 1895 of file valve.cc.

{
   return 2;
}

void Flow_valve::Output ( OutputHandler &  OH ) const

virtual

Reimplemented from ToBeOutput.

Definition at line 2266 of file valve.cc.

References ToBeOutput::bToBeOutput(), flow1, flow2, flow3, WithLabel::GetLabel(), OutputHandler::Hydraulic(), s, v, and vp.

{
   if (bToBeOutput()) { 
      std::ostream& out = OH.Hydraulic();
      out << std::setw(8) << GetLabel()
        << " " << s  << " " << v  << " "<< vp  
        << " " << flow1  << " "<< flow2  << " "<< flow3  << std::endl;
   }  
}

Here is the call graph for this function:

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

virtual

Implements Elem.

Definition at line 1889 of file valve.cc.

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

void Flow_valve::SetValue ( DataManager *  pDM, VectorHandler &  X, VectorHandler &  XP, SimulationEntity::Hints *  ph = 0  )

virtual

Reimplemented from SimulationEntity.

Definition at line 2276 of file valve.cc.

References DofOwnerOwner::iGetFirstIndex(), and VectorHandler::PutCoef().

{   
   integer i = iGetFirstIndex();
   X.PutCoef(i+1, 0.);
   X.PutCoef(i+2, 0.);
   XP.PutCoef(i+1, 0.);
   XP.PutCoef(i+2, 0.);  
}

Here is the call graph for this function:

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

virtual

Implements Elem.

Definition at line 1909 of file valve.cc.

{
   *piNumRows = 5;
   *piNumCols = 5; 
}

Member Data Documentation

doublereal Flow_valve::area_diaf

private

Definition at line 521 of file valve.h.

Referenced by AssJac(), and AssRes().

doublereal Flow_valve::area_max

private

Definition at line 524 of file valve.h.

Referenced by AssJac(), and AssRes().

doublereal Flow_valve::area_pipe

private

Definition at line 523 of file valve.h.

Referenced by AssJac(), and AssRes().

doublereal Flow_valve::c1

private

Definition at line 532 of file valve.h.

Referenced by AssJac(), and AssRes().

doublereal Flow_valve::c2

private

Definition at line 533 of file valve.h.

Referenced by AssJac(), and AssRes().

doublereal Flow_valve::c3

private

Definition at line 534 of file valve.h.

Referenced by AssJac(), and AssRes().

doublereal Flow_valve::c4

private

Definition at line 535 of file valve.h.

Referenced by AssRes().

doublereal Flow_valve::c_acc

private

Definition at line 531 of file valve.h.

Referenced by AssRes().

doublereal Flow_valve::c_spost

private

Definition at line 529 of file valve.h.

Referenced by AssRes().

doublereal Flow_valve::c_vel

private

Definition at line 530 of file valve.h.

Referenced by AssRes().

doublereal Flow_valve::cf1

private

Definition at line 536 of file valve.h.

Referenced by AssJac(), and AssRes().

doublereal Flow_valve::cf2

private

Definition at line 537 of file valve.h.

Referenced by AssJac(), and AssRes().

doublereal Flow_valve::cf3

private

Definition at line 538 of file valve.h.

Referenced by AssJac(), and AssRes().

doublereal Flow_valve::cf4

private

Definition at line 539 of file valve.h.

Referenced by AssRes().

doublereal Flow_valve::flow1

private

Definition at line 546 of file valve.h.

Referenced by AssRes(), and Output().

doublereal Flow_valve::flow2

private

Definition at line 547 of file valve.h.

Referenced by AssRes(), and Output().

doublereal Flow_valve::flow3

private

Definition at line 548 of file valve.h.

Referenced by AssRes(), and Output().

doublereal Flow_valve::force0

private

Definition at line 526 of file valve.h.

Referenced by AssRes().

doublereal Flow_valve::h

private

Definition at line 541 of file valve.h.

Referenced by AssJac(), AssRes(), and Flow_valve().

doublereal Flow_valve::Kappa

private

Definition at line 525 of file valve.h.

Referenced by AssJac(), and AssRes().

doublereal Flow_valve::mass

private

Definition at line 522 of file valve.h.

Referenced by AssJac(), and AssRes().

const PressureNode* Flow_valve::pNode1

private

Definition at line 518 of file valve.h.

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

const PressureNode* Flow_valve::pNode2

private

Definition at line 519 of file valve.h.

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

const PressureNode* Flow_valve::pNode3

private

Definition at line 520 of file valve.h.

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

doublereal Flow_valve::s

private

Definition at line 542 of file valve.h.

Referenced by AssJac(), AssRes(), and Output().

doublereal Flow_valve::s_max

private

Definition at line 528 of file valve.h.

Referenced by AssRes(), and Flow_valve().

doublereal Flow_valve::sp

private

Definition at line 543 of file valve.h.

Referenced by AssRes().

doublereal Flow_valve::v

private

Definition at line 544 of file valve.h.

Referenced by AssJac(), AssRes(), and Output().

doublereal Flow_valve::vp

private

Definition at line 545 of file valve.h.

Referenced by AssRes(), and Output().

doublereal Flow_valve::width

private

Definition at line 527 of file valve.h.

Referenced by AssJac(), and AssRes().

---

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

• valve.h
• valve.cc