#include <valve.h>

Inheritance diagram for Pressure_flow_control_valve:

Collaboration diagram for Pressure_flow_control_valve:

Public Member Functions
	Pressure_flow_control_valve (unsigned int uL, const DofOwner pD, HydraulicFluid hf, const PressureNode p1, const PressureNode p2, const PressureNode p3, const PressureNode p4, const PressureNode p5, const PressureNode p6, const DriveCaller *pDC, doublereal s0, doublereal s_mx, doublereal W, doublereal Loss_A, doublereal Valve_d, doublereal Valve_rho, doublereal cs, doublereal cv, doublereal ca, flag fOut)

	~Pressure_flow_control_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)

Public Member Functions inherited from DriveOwner
	DriveOwner (const DriveCaller *pDC=0)

	DriveOwner (const DriveOwner &drive)

virtual	~DriveOwner (void)

void	Set (const DriveCaller *pDC)

DriveCaller *	pGetDriveCaller (void) const

doublereal	dGet (const doublereal &dVar) const

doublereal	dGet (void) const

bool	bIsDifferentiable (void) const

doublereal	dGetP (const doublereal &dVar) const

doublereal	dGetP (void) const

Private Attributes
const PressureNode *	pNode1

const PressureNode *	pNode2

const PressureNode *	pNode3

const PressureNode *	pNode4

const PressureNode *	pNode5

const PressureNode *	pNode6

doublereal	start

doublereal	c_spost

doublereal	c_vel

doublereal	c_acc

doublereal	Force

doublereal	Cd

doublereal	width

doublereal	loss_area

doublereal	valve_diameter

doublereal	valve_density

doublereal	valve_area

doublereal	s_max

doublereal	A1

doublereal	A2

doublereal	A3

doublereal	A4

doublereal	Mass

doublereal	flow1

doublereal	flow2

doublereal	flow3

doublereal	flow4

doublereal	flow5

doublereal	flow6

doublereal	s

doublereal	v

doublereal	sp

doublereal	vp

doublereal	c1

doublereal	c2

doublereal	c3

doublereal	cf1

doublereal	cf2

doublereal	cf3

doublereal	deltaP

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

Protected Attributes inherited from DriveOwner
DriveCaller *	pDriveCaller

Detailed Description

Definition at line 324 of file valve.h.

Constructor & Destructor Documentation

Pressure_flow_control_valve::Pressure_flow_control_valve	(	unsigned int	uL,
		const DofOwner *	pD,
		HydraulicFluid *	hf,
		const PressureNode *	p1,
		const PressureNode *	p2,
		const PressureNode *	p3,
		const PressureNode *	p4,
		const PressureNode *	p5,
		const PressureNode *	p6,
		const DriveCaller *	pDC,
		doublereal	s0,
		doublereal	s_mx,
		doublereal	W,
		doublereal	Loss_A,
		doublereal	Valve_d,
		doublereal	Valve_rho,
		doublereal	cs,
		doublereal	cv,
		doublereal	ca,
		flag	fOut
	)

Definition at line 1072 of file valve.cc.

References ASSERT, Cd, PressureNode::GetNodeType(), Node::HYDRAULIC, M_PI, pNode1, pNode2, pNode3, pNode4, pNode5, pNode6, valve_area, valve_density, and valve_diameter.

 : Elem(uL, fOut),
 HydraulicElem(uL, pDO, hf, fOut), DriveOwner(pDC),
 pNode1(p1), pNode2(p2), pNode3(p3), pNode4(p4), pNode5(p5), pNode6(p6),
 start(s0), s_max(s_mx), width(W), loss_area(Loss_A), 
 valve_diameter(Valve_d), valve_density(Valve_rho),
 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(pNode4 != NULL);
    ASSERT(pNode4->GetNodeType() == Node::HYDRAULIC);
    ASSERT(pNode5 != NULL);
    ASSERT(pNode5->GetNodeType() == Node::HYDRAULIC);
    ASSERT(pNode6 != NULL);
    ASSERT(pNode6->GetNodeType() == Node::HYDRAULIC);
    ASSERT(s0 >= 0.);
    ASSERT(Valve_rho > std::numeric_limits<doublereal>::epsilon());
    ASSERT(Valve_d > std::numeric_limits<doublereal>::epsilon());
    ASSERT(W > std::numeric_limits<doublereal>::epsilon());
    ASSERT(Loss_A >= 0.);
    ASSERT(s_mx >= 0.);
     
    Cd = .611;                /* coefficiente di perdita */
    Mass = 1.175*valve_diameter*valve_diameter*valve_diameter*valve_density*M_PI; // massa della valvola 
    valve_area=0.25*valve_diameter*valve_diameter*M_PI;
 }

Here is the call graph for this function:

Pressure_flow_control_valve::~Pressure_flow_control_valve ( void )

Definition at line 1117 of file valve.cc.

References NO_OP.

 {
    NO_OP;
 }

Member Function Documentation

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

virtual

Implements Elem.

Definition at line 1154 of file valve.cc.

References A1, A2, A3, A4, c1, c2, c3, Cd, cf1, cf2, cf3, copysign(), DEBUGCOUT, deltaP, HydraulicFluid::dGetDensity(), ScalarAlgebraicNode::dGetX(), grad::fabs(), HydraulicElem::HF, Node::iGetFirstColIndex(), DofOwnerOwner::iGetFirstIndex(), Node::iGetFirstRowIndex(), pNode1, pNode2, pNode3, pNode4, pNode5, pNode6, FullSubMatrixHandler::PutCoef(), FullSubMatrixHandler::PutColIndex(), FullSubMatrixHandler::PutRowIndex(), FullSubMatrixHandler::ResizeReset(), s, VariableSubMatrixHandler::SetFull(), sp, grad::sqrt(), v, valve_area, valve_diameter, and width.

 {
    DEBUGCOUT("Entering Control_valve::AssJac()" << std::endl);
    
    FullSubMatrixHandler& WM = WorkMat.SetFull();
    WM.ResizeReset(8, 8);
    
    integer iNode1RowIndex = pNode1->iGetFirstRowIndex()+1;
    integer iNode2RowIndex = pNode2->iGetFirstRowIndex()+1;
    integer iNode3RowIndex = pNode3->iGetFirstRowIndex()+1;
    integer iNode4RowIndex = pNode4->iGetFirstRowIndex()+1;
    integer iNode5RowIndex = pNode5->iGetFirstRowIndex()+1;
    integer iNode6RowIndex = pNode6->iGetFirstRowIndex()+1;
    integer iNode1ColIndex = pNode1->iGetFirstColIndex()+1;
    integer iNode2ColIndex = pNode2->iGetFirstColIndex()+1;
    integer iNode3ColIndex = pNode3->iGetFirstColIndex()+1;
    integer iNode4ColIndex = pNode4->iGetFirstColIndex()+1;
    integer iNode5ColIndex = pNode5->iGetFirstColIndex()+1;
    integer iNode6ColIndex = pNode6->iGetFirstColIndex()+1;
    integer iFirstIndex = iGetFirstIndex();
    
    WM.PutRowIndex(1, iNode1RowIndex);
    WM.PutRowIndex(2, iNode2RowIndex);
    WM.PutRowIndex(3, iNode3RowIndex);
    WM.PutRowIndex(4, iNode4RowIndex);
    WM.PutRowIndex(5, iNode5RowIndex);
    WM.PutRowIndex(6, iNode6RowIndex);
    WM.PutColIndex(1, iNode1ColIndex);
    WM.PutColIndex(2, iNode2ColIndex);
    WM.PutColIndex(3, iNode3ColIndex);
    WM.PutColIndex(4, iNode4ColIndex);
    WM.PutColIndex(5, iNode5ColIndex);
    WM.PutColIndex(6, iNode6ColIndex);
  
    WM.PutRowIndex(7, iFirstIndex+1);
    WM.PutColIndex(7, iFirstIndex+1);
    WM.PutRowIndex(8, iFirstIndex+2);
    WM.PutColIndex(8, iFirstIndex+2);
    
    doublereal p1 = pNode1->dGetX();
    doublereal p2 = pNode2->dGetX();
    doublereal p3 = pNode3->dGetX();
    doublereal p4 = pNode4->dGetX();
    doublereal p5 = pNode5->dGetX();
    doublereal p6 = pNode6->dGetX();
    doublereal density = HF->dGetDensity();
 
    s = XCurr(iFirstIndex+1); /* spostamento */
    v = XCurr(iFirstIndex+2); /* velocita' */
    
    doublereal jumpPres12 = fabs(p1-p2); /* salto di pressione nodo1 & nodo2 */
    doublereal jumpPres13 = fabs(p1-p3); /* salto di pressione nodo1 & nodo3 */
    doublereal jumpPres24 = fabs(p2-p4); /* salto di pressione nodo2 & nodo4 */
    doublereal jumpPres34 = fabs(p3-p4); /* salto di pressione nodo3 & nodo4 */
    
    if (jumpPres12 < 1.e8*std::numeric_limits<doublereal>::epsilon()) {
       jumpPres12 = 1.e8*std::numeric_limits<doublereal>::epsilon();
    }
    if (jumpPres13 < 1.e8*std::numeric_limits<doublereal>::epsilon()) {
       jumpPres13 = 1.e8*std::numeric_limits<doublereal>::epsilon();
    }
    if (jumpPres24 < 1.e8*std::numeric_limits<doublereal>::epsilon()) {
       jumpPres24 = 1.e8*std::numeric_limits<doublereal>::epsilon();
    }
    if (jumpPres34 < 1.e8*std::numeric_limits<doublereal>::epsilon()) {
       jumpPres34 = 1.e8*std::numeric_limits<doublereal>::epsilon();
    }
    doublereal Jac17;
    doublereal Jac27;
    doublereal Jac37;
    doublereal Jac47;
     
    doublereal primo   = Cd*A1/sqrt(2.*jumpPres12/density);
    doublereal secondo = Cd*A2/sqrt(2.*jumpPres13/density);
    doublereal terzo   = Cd*A4/sqrt(2.*jumpPres24/density);
    doublereal quarto  = Cd*A3/sqrt(2.*jumpPres34/density);
 
    doublereal Jac11 = -primo-secondo;
    doublereal Jac12 = primo;
    doublereal Jac13 = secondo;
    doublereal Jac14 = 0.;
    doublereal Jac21 = primo;
    doublereal Jac22 = -primo-terzo;
    doublereal Jac23 = 0.;
    doublereal Jac24 = terzo;
    doublereal Jac31 = secondo;
    doublereal Jac32 = 0.;
    doublereal Jac33 = -secondo-quarto;
    doublereal Jac34 = quarto;
    doublereal Jac41 = 0.;
    doublereal Jac42 = terzo;
    doublereal Jac43 = quarto;
    doublereal Jac44 = -terzo-quarto;
   
    doublereal Jac57 = -valve_area;
    doublereal Jac67 = +valve_area;
      
    doublereal costante = density*Cd*width*valve_diameter;
    
    doublereal Jac71  = -.2*costante/sqrt(density*deltaP)*sp-.43*width*s; 
    doublereal Jac72  = 0.;
    doublereal Jac73  = 0.;
    doublereal Jac74  = 0.;
    doublereal Jac77  = -.4*valve_diameter*Cd*width*sqrt(density*deltaP)-dCoef*.43*width*deltaP-dCoef*c1-c2-dCoef*cf1-cf2;
    doublereal Jac78  = -Mass-c3-cf3;
 
    if (s > 0.) {  /* collegamento diretto 1->2   3->4 */
       Jac17 = dCoef*density*Cd*width*copysign(sqrt(2.*jumpPres12/density), p1-p2);
       Jac27 = -Jac17;
       Jac37 = dCoef*density*Cd*width*copysign(sqrt(2.*jumpPres34/density), p3-p4);
       Jac47 = -Jac37;
    } else {       /* collegamento inverso 1->3     2->4 */
       Jac17 = -dCoef*density*Cd*width*copysign(sqrt(2.*jumpPres13/density), p1-p3);
       Jac27 = -dCoef*density*Cd*width*copysign(sqrt(2.*jumpPres24/density), p2-p4);
       Jac37 = -Jac17;
       Jac47 = -Jac27;
    }
    
    doublereal Jac87  = -1;
    doublereal Jac88  = dCoef;  
 
    WM.PutCoef(1, 1, Jac11);
    WM.PutCoef(1, 2, Jac12);
    WM.PutCoef(1, 3, Jac13);
    WM.PutCoef(1, 4, Jac14);
    WM.PutCoef(1, 7, Jac17);
    WM.PutCoef(2, 1, Jac21);
    WM.PutCoef(2, 2, Jac22);
    WM.PutCoef(2, 3, Jac23);
    WM.PutCoef(2, 4, Jac24);
    WM.PutCoef(2, 7, Jac27);
    WM.PutCoef(3, 1, Jac31);
    WM.PutCoef(3, 2, Jac32);
    WM.PutCoef(3, 3, Jac33);
    WM.PutCoef(3, 4, Jac34);
    WM.PutCoef(3, 7, Jac37);
    WM.PutCoef(4, 1, Jac41);
    WM.PutCoef(4, 2, Jac42);
    WM.PutCoef(4, 3, Jac43);
    WM.PutCoef(4, 4, Jac44);
    WM.PutCoef(4, 7, Jac47);
    WM.PutCoef(5, 7, Jac67);
    WM.PutCoef(6, 7, Jac67);
    WM.PutCoef(7, 1, Jac71);
    WM.PutCoef(7, 2, Jac72);
    WM.PutCoef(7, 3, Jac73);
    WM.PutCoef(7, 7, Jac77);
    WM.PutCoef(7, 8, Jac78);
    WM.PutCoef(8, 7, Jac87);
    WM.PutCoef(8, 8, Jac88);
  
    return WorkMat;
 }

Here is the call graph for this function:

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

virtual

Implements Elem.

Definition at line 1312 of file valve.cc.

References A1, A2, A3, A4, c1, c2, c3, Cd, cf1, cf2, cf3, copysign(), DEBUGCOUT, deltaP, DriveCaller::dGet(), HydraulicFluid::dGetDensity(), ScalarAlgebraicNode::dGetX(), grad::fabs(), flow1, flow2, flow3, flow4, flow5, flow6, HydraulicElem::HF, DofOwnerOwner::iGetFirstIndex(), Node::iGetFirstRowIndex(), loss_area, DriveOwner::pGetDriveCaller(), pNode1, pNode2, pNode3, pNode4, pNode5, pNode6, SubVectorHandler::PutItem(), VectorHandler::Resize(), s, s_max, sp, grad::sqrt(), v, valve_area, valve_density, valve_diameter, vp, and width.

 {
    DEBUGCOUT("Entering Pressure_flow_control_valve::AssRes()" << std::endl);
    
    WorkVec.Resize(8);
    
    integer iNode1RowIndex = pNode1->iGetFirstRowIndex()+1;
    integer iNode2RowIndex = pNode2->iGetFirstRowIndex()+1;
    integer iNode3RowIndex = pNode3->iGetFirstRowIndex()+1;
    integer iNode4RowIndex = pNode4->iGetFirstRowIndex()+1;
    integer iNode5RowIndex = pNode5->iGetFirstRowIndex()+1;
    integer iNode6RowIndex = pNode6->iGetFirstRowIndex()+1;
    integer iFirstIndex = iGetFirstIndex();
    
    doublereal p1 = pNode1->dGetX();
    doublereal p2 = pNode2->dGetX();
    doublereal p3 = pNode3->dGetX();
    doublereal p4 = pNode4->dGetX();
    doublereal p5 = pNode5->dGetX();
    doublereal p6 = pNode6->dGetX();
    doublereal density = HF->dGetDensity();
 
    Force = pGetDriveCaller()->dGet();
    
    s = XCurr(iFirstIndex+1);       /* spostamento */
    v = XCurr(iFirstIndex+2);       /* velocita' */
    sp = XPrimeCurr(iFirstIndex+1); /* velocita' */
    vp = XPrimeCurr(iFirstIndex+2); /* accelerazione */
    
    doublereal jumpPres12 = fabs(p1-p2); /* salto di pressione nodo1 & nodo2 */
    doublereal jumpPres13 = fabs(p1-p3); /* salto di pressione nodo1 & nodo3 */
    doublereal jumpPres24 = fabs(p2-p4); /* salto di pressione nodo2 & nodo4 */
    doublereal jumpPres34 = fabs(p3-p4); /* salto di pressione nodo3 & nodo4 */
 
    /* qui decido se sono a fondo od ad inizio corsa e di conseguenza calcolo
     * i giusti coefficienti */
 
    doublereal x = s;  /* usato per calcolare le aree */
    
    doublereal area_max = width*s_max;         // Area massima
    doublereal area_min = area_max*loss_area;  // Area minimo 
    doublereal deltaA = x*width;
 
    if (x > 0.) { 
       A1 = area_min+deltaA;
       A2 = area_min;
       A3 = area_min+deltaA;
       A4 = area_min;
    } else { /* ho deltaA negativo */
       A1 = area_min;
       A2 = area_min-deltaA;
       A3 = area_min;
       A4 = area_min-deltaA;
    }
    
    doublereal Q12 = density*Cd*A1*copysign(sqrt(2.*jumpPres12/density), p1-p2);
    doublereal Q13 = density*Cd*A2*copysign(sqrt(2.*jumpPres13/density), p1-p3);
    doublereal Q24 = density*Cd*A4*copysign(sqrt(2.*jumpPres24/density), p2-p4);
    doublereal Q34 = density*Cd*A3*copysign(sqrt(2.*jumpPres34/density), p3-p4);
 
    doublereal Res_1 = Q12+Q13;
    doublereal Res_2 = -Q12+Q24;
    doublereal Res_3 = -Q13+Q34;
    doublereal Res_4 = -Q34-Q24;
    doublereal Res_5 = valve_area*sp;
 #warning "????????????? Res_6 = -Res_6 ?"
    doublereal Res_6 = -Res_6;
  
    deltaP = p1;   
    if (deltaP == 0.) {
       deltaP = 10.; /* evito di dividere per zero nello jacobiano */
    }
    doublereal C = .4*valve_diameter*Cd*width*sqrt(density*deltaP);
    doublereal K = .43*width*deltaP;
    doublereal Res_7 = -Force+Mass*vp+C*sp+K*s+c1*s+c2*sp+c3*vp+cf1*(s-s_max)+cf2*sp+cf3*vp;
    doublereal Res_8 = sp-v;
 
    flow1 = -Res_1;
    flow2 = -Res_2;
    flow3 = -Res_3;
    flow4 = -Res_4;
    flow5 = -Res_5;
    flow6 = -Res_6;
 
 #ifdef HYDR_DEVEL
    DEBUGCOUT("Force:         " << Force << std::endl);
    DEBUGCOUT("X equilibrio:  " << Force/K << std::endl);
    DEBUGCOUT("A1:            " << A1 << std::endl);
    DEBUGCOUT("A2:            " << A2  << std::endl);
    DEBUGCOUT("A3:            " << A3 << std::endl);
    DEBUGCOUT("A4:            " << A4 << std::endl);
    DEBUGCOUT("p1:            " << p1 << std::endl);
    DEBUGCOUT("p2:            " << p2 << std::endl);
    DEBUGCOUT("p3:            " << p3 << std::endl);
    DEBUGCOUT("p4:            " << p4 << std::endl);
    DEBUGCOUT("Cd:            " << Cd << std::endl);
    DEBUGCOUT("s_max :        " << s_max << std::endl);
    DEBUGCOUT("x :            " << x << std::endl);
    DEBUGCOUT("s :            " << s << std::endl);
    DEBUGCOUT("sp:            " << sp << std::endl);
    DEBUGCOUT("v :            " << v << std::endl);
    DEBUGCOUT("vp:            " << vp << std::endl);
    DEBUGCOUT("Valve_diameter:" << valve_diameter << std::endl);
    DEBUGCOUT("massa:         " << Mass << std::endl);
    DEBUGCOUT("smorzatore:    " << C << std::endl);
    DEBUGCOUT("molla:         " << K << std::endl);
    DEBUGCOUT("density:       " << density << std::endl);
    DEBUGCOUT("Valve_density: " << valve_density << std::endl);
    DEBUGCOUT("Area_max:      " << area_max << std::endl);
    DEBUGCOUT("Width:         " << width << std::endl);
    DEBUGCOUT("Loss_area:     " << loss_area << std::endl);
    DEBUGCOUT("c1:            " << c1 << std::endl);
    DEBUGCOUT("c2:            " << c2 << std::endl);
    DEBUGCOUT("c3:            " << c3 << std::endl);
    DEBUGCOUT("cf1:           " << cf1 << std::endl);
    DEBUGCOUT("cf2:           " << cf2 << std::endl);
    DEBUGCOUT("cf3:           " << cf3 << std::endl);
    DEBUGCOUT("Q12:           " << Q12 << std::endl);
    DEBUGCOUT("Q13:           " << Q13 << std::endl);
    DEBUGCOUT("Q24:           " << Q24 << std::endl);
    DEBUGCOUT("Q34:           " << Q34 << 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 (portata nodo4): " << -Res_4 << std::endl); 
    DEBUGCOUT("-Res_5 (portata nodo3): " << -Res_5 << std::endl); 
    DEBUGCOUT("-Res_6 (portata nodo4): " << -Res_6 << std::endl); 
    DEBUGCOUT("-Res_6   eq dinamica  : " << -Res_7 << std::endl); 
    DEBUGCOUT("-Res_7 sp-v           : " << -Res_8 << std::endl); 
 #endif
    
    WorkVec.PutItem(1, iNode1RowIndex, Res_1);
    WorkVec.PutItem(2, iNode2RowIndex, Res_2);         
    WorkVec.PutItem(3, iNode3RowIndex, Res_3);
    WorkVec.PutItem(4, iNode4RowIndex, Res_4);    
    WorkVec.PutItem(5, iNode5RowIndex, Res_5);
    WorkVec.PutItem(6, iNode6RowIndex, Res_6);    
    WorkVec.PutItem(7, iFirstIndex+1, Res_7);
    WorkVec.PutItem(8, iFirstIndex+2, Res_8);    
     
    return WorkVec;
 }

Here is the call graph for this function:

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

inlinevirtual

Reimplemented from Elem.

Definition at line 414 of file valve.h.

References pNode1, pNode2, pNode3, pNode4, pNode5, and pNode6.

                                                                                  {
      connectedNodes.resize(6);
      connectedNodes[0] = pNode1;
      connectedNodes[1] = pNode2;
      connectedNodes[2] = pNode3;
      connectedNodes[3] = pNode4;
      connectedNodes[4] = pNode5;
      connectedNodes[5] = pNode6;
   };

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

virtual

Reimplemented from Elem.

Definition at line 1139 of file valve.cc.

References ASSERT, and DofOrder::DIFFERENTIAL.

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

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

virtual

Implements HydraulicElem.

Definition at line 1123 of file valve.cc.

References HydraulicElem::PRESSURE_FLOW_CONTROL_VALVE.

 {
    return HydraulicElem::PRESSURE_FLOW_CONTROL_VALVE;
 }

unsigned int Pressure_flow_control_valve::iGetNumDof ( void ) const

virtual

Reimplemented from Elem.

Definition at line 1134 of file valve.cc.

 {
    return 2;
 }

void Pressure_flow_control_valve::Output ( OutputHandler & OH ) const

virtual

Reimplemented from ToBeOutput.

Definition at line 1460 of file valve.cc.

References A1, A2, A3, A4, ToBeOutput::bToBeOutput(), flow1, flow2, flow3, flow4, flow5, flow6, WithLabel::GetLabel(), OutputHandler::Hydraulic(), s, sp, and vp.

 {
    if (bToBeOutput()) { 
       std::ostream& out = OH.Hydraulic();
       out << std::setw(8) << GetLabel()
         << " " << s << " "  << sp << " "  << vp 
         << " " << flow1 << " " << flow2 << " " << flow3 << " " << flow4
         << " " << flow5 << " " << flow6
         << " " << A1 << " "  << A2 << " "  << A3 << " " << A4 << " " << std::endl;
    }
 }

Here is the call graph for this function:

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

virtual

Implements Elem.

Definition at line 1129 of file valve.cc.

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

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

virtual

Reimplemented from SimulationEntity.

Definition at line 1474 of file valve.cc.

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

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

Here is the call graph for this function:

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

virtual

Implements Elem.

Definition at line 1147 of file valve.cc.

 {
    *piNumRows = 8; 
    *piNumCols = 8; 
 }

Member Data Documentation

doublereal Pressure_flow_control_valve::A1

private

Definition at line 346 of file valve.h.

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

doublereal Pressure_flow_control_valve::A2

private

Definition at line 347 of file valve.h.

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

doublereal Pressure_flow_control_valve::A3

private

Definition at line 348 of file valve.h.

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

doublereal Pressure_flow_control_valve::A4

private

Definition at line 349 of file valve.h.

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

doublereal Pressure_flow_control_valve::c1

private

Definition at line 363 of file valve.h.

Referenced by AssJac(), and AssRes().

doublereal Pressure_flow_control_valve::c2

private

Definition at line 364 of file valve.h.

Referenced by AssJac(), and AssRes().

doublereal Pressure_flow_control_valve::c3

private

Definition at line 365 of file valve.h.

Referenced by AssJac(), and AssRes().

doublereal Pressure_flow_control_valve::c_acc

private

Definition at line 335 of file valve.h.

doublereal Pressure_flow_control_valve::c_spost

private

Definition at line 333 of file valve.h.

doublereal Pressure_flow_control_valve::c_vel

private

Definition at line 334 of file valve.h.

doublereal Pressure_flow_control_valve::Cd

private

Definition at line 338 of file valve.h.

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

doublereal Pressure_flow_control_valve::cf1

private

Definition at line 366 of file valve.h.

Referenced by AssJac(), and AssRes().

doublereal Pressure_flow_control_valve::cf2

private

Definition at line 367 of file valve.h.

Referenced by AssJac(), and AssRes().

doublereal Pressure_flow_control_valve::cf3

private

Definition at line 368 of file valve.h.

Referenced by AssJac(), and AssRes().

doublereal Pressure_flow_control_valve::deltaP

private

Definition at line 369 of file valve.h.

Referenced by AssJac(), and AssRes().

doublereal Pressure_flow_control_valve::flow1

private

Definition at line 352 of file valve.h.

Referenced by AssRes(), and Output().

doublereal Pressure_flow_control_valve::flow2

private

Definition at line 353 of file valve.h.

Referenced by AssRes(), and Output().

doublereal Pressure_flow_control_valve::flow3

private

Definition at line 354 of file valve.h.

Referenced by AssRes(), and Output().

doublereal Pressure_flow_control_valve::flow4

private

Definition at line 355 of file valve.h.

Referenced by AssRes(), and Output().

doublereal Pressure_flow_control_valve::flow5

private

Definition at line 356 of file valve.h.

Referenced by AssRes(), and Output().

doublereal Pressure_flow_control_valve::flow6

private

Definition at line 357 of file valve.h.

Referenced by AssRes(), and Output().

doublereal Pressure_flow_control_valve::Force

private

Definition at line 336 of file valve.h.

doublereal Pressure_flow_control_valve::loss_area

private

Definition at line 340 of file valve.h.

Referenced by AssRes().

doublereal Pressure_flow_control_valve::Mass

private

Definition at line 351 of file valve.h.

const PressureNode* Pressure_flow_control_valve::pNode1

private

Definition at line 326 of file valve.h.

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

const PressureNode* Pressure_flow_control_valve::pNode2

private

Definition at line 327 of file valve.h.

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

const PressureNode* Pressure_flow_control_valve::pNode3

private

Definition at line 328 of file valve.h.

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

const PressureNode* Pressure_flow_control_valve::pNode4

private

Definition at line 329 of file valve.h.

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

const PressureNode* Pressure_flow_control_valve::pNode5

private

Definition at line 330 of file valve.h.

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

const PressureNode* Pressure_flow_control_valve::pNode6

private

Definition at line 331 of file valve.h.

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

doublereal Pressure_flow_control_valve::s

private

Definition at line 359 of file valve.h.

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

doublereal Pressure_flow_control_valve::s_max

private

Definition at line 344 of file valve.h.

Referenced by AssRes().

doublereal Pressure_flow_control_valve::sp

private

Definition at line 361 of file valve.h.

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

doublereal Pressure_flow_control_valve::start

private

Definition at line 332 of file valve.h.

Referenced by SetValue().

doublereal Pressure_flow_control_valve::v

private

Definition at line 360 of file valve.h.

Referenced by AssJac(), and AssRes().

doublereal Pressure_flow_control_valve::valve_area

private

Definition at line 343 of file valve.h.

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

doublereal Pressure_flow_control_valve::valve_density

private

Definition at line 342 of file valve.h.

Referenced by AssRes(), and Pressure_flow_control_valve().

doublereal Pressure_flow_control_valve::valve_diameter

private

Definition at line 341 of file valve.h.

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

doublereal Pressure_flow_control_valve::vp

private

Definition at line 362 of file valve.h.

Referenced by AssRes(), and Output().

doublereal Pressure_flow_control_valve::width

private

Definition at line 339 of file valve.h.

Referenced by AssJac(), and AssRes().

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