#include <valve.h>

Inheritance diagram for Dynamic_control_valve:

Collaboration diagram for Dynamic_control_valve:

Public Member Functions
	Dynamic_control_valve (unsigned int uL, const DofOwner pD, HydraulicFluid hf, const PressureNode p1, const PressureNode p2, const PressureNode p3, const PressureNode p4, 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)

	~Dynamic_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

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	s_max

doublereal	A1

doublereal	A2

doublereal	A3

doublereal	A4

doublereal	Mass

doublereal	flow1

doublereal	flow2

doublereal	flow3

doublereal	flow4

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

Constructor & Destructor Documentation

Dynamic_control_valve::Dynamic_control_valve	(	unsigned int	uL,
		const DofOwner *	pD,
		HydraulicFluid *	hf,
		const PressureNode *	p1,
		const PressureNode *	p2,
		const PressureNode *	p3,
		const PressureNode *	p4,
		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 628 of file valve.cc.

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

 : Elem(uL, fOut),
 HydraulicElem(uL, pDO, hf, fOut), DriveOwner(pDC),
 pNode1(p1), pNode2(p2), pNode3(p3), pNode4(p4),
 start(s0), c_spost(cs), c_vel(cv), c_acc(ca),
 width(W), loss_area(Loss_A),
 valve_diameter(Valve_d), valve_density(Valve_rho), s_max(s_mx)
 {
    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(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 
    
    // W = .005;                 /* larghezza del condotto (m): A=x*W 0.005; */
    // diameter = .01;           /* diametro della valvola (m) */
    // densityvalve = 7900.;     /* densita' del corpo della valvola (acciaio) (kg/m^3) */
    // Mass = 1.175*diameter*diameter*diameter*densityvalve*M_PI; /* massa della valvola (Kg) */
    // s_max = 2.;               /* corsa massima della valvola (m) */
 }

Here is the call graph for this function:

Dynamic_control_valve::~Dynamic_control_valve ( void )

Definition at line 675 of file valve.cc.

References NO_OP.

 {
    NO_OP;
 }

Member Function Documentation

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

virtual

Implements Elem.

Definition at line 712 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, FullSubMatrixHandler::PutCoef(), FullSubMatrixHandler::PutColIndex(), FullSubMatrixHandler::PutRowIndex(), FullSubMatrixHandler::ResizeReset(), s, VariableSubMatrixHandler::SetFull(), sp, grad::sqrt(), v, valve_diameter, and width.

 {
    DEBUGCOUT("Entering Control_valve::AssJac()" << std::endl);
    
    FullSubMatrixHandler& WM = WorkMat.SetFull();
    WM.ResizeReset(6, 6);
    
    integer iNode1RowIndex = pNode1->iGetFirstRowIndex()+1;
    integer iNode2RowIndex = pNode2->iGetFirstRowIndex()+1;
    integer iNode3RowIndex = pNode3->iGetFirstRowIndex()+1;
    integer iNode4RowIndex = pNode4->iGetFirstRowIndex()+1;
    integer iNode1ColIndex = pNode1->iGetFirstColIndex()+1;
    integer iNode2ColIndex = pNode2->iGetFirstColIndex()+1;
    integer iNode3ColIndex = pNode3->iGetFirstColIndex()+1;
    integer iNode4ColIndex = pNode4->iGetFirstColIndex()+1;
    integer iFirstIndex = iGetFirstIndex();
    
    WM.PutRowIndex(1, iNode1RowIndex);
    WM.PutRowIndex(2, iNode2RowIndex);
    WM.PutRowIndex(3, iNode3RowIndex);
    WM.PutRowIndex(4, iNode4RowIndex);
    WM.PutColIndex(1, iNode1ColIndex);
    WM.PutColIndex(2, iNode2ColIndex);
    WM.PutColIndex(3, iNode3ColIndex);
    WM.PutColIndex(4, iNode4ColIndex);
    
    WM.PutRowIndex(5, iFirstIndex+1);
    WM.PutColIndex(5, iFirstIndex+1);
    WM.PutRowIndex(6, iFirstIndex+2);
    WM.PutColIndex(6, iFirstIndex+2);
    
    doublereal p1 = pNode1->dGetX();
    doublereal p2 = pNode2->dGetX();
    doublereal p3 = pNode3->dGetX();
    doublereal p4 = pNode4->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 Jac15;
    doublereal Jac25;
    doublereal Jac35;
    doublereal Jac45;
     
    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 costante = density*Cd*width*valve_diameter;
 
 #if 0
    doublereal Jac51  = -.2*costante*sp/sqrt(density*deltaP)-.43*width*s; 
    doublereal Jac52  = -Jac51;
    doublereal Jac53  = Jac51;
    doublereal Jac54  = -Jac51;
    doublereal Jac55  = -.4*valve_diameter*Cd*width*sqrt(density*deltaP)-dCoef*.43*width*deltaP-dCoef*c1-c2-dCoef*cf1-cf2;
    doublereal Jac56  = -Mass-c3-cf3;
 #endif
    
    doublereal Jac51  = -.2*costante*sp/sqrt(density*deltaP)-.43*width*s; 
    doublereal Jac52  = 0.;
    doublereal Jac53  = 0.;
    doublereal Jac54  = 0.;
    doublereal Jac55  = -.4*valve_diameter*Cd*width*sqrt(density*deltaP)-dCoef*.43*width*deltaP-dCoef*c1-c2-dCoef*cf1-cf2;
    doublereal Jac56  = -Mass-c3-cf3;
    
    if (s > 0.) {  /* collegamento diretto 1->2   3->4 */
       Jac15 = dCoef*density*Cd*width*copysign(sqrt(2.*jumpPres12/density), p1-p2);
       Jac25 = -Jac15;
       Jac35 = dCoef*density*Cd*width*copysign(sqrt(2.*jumpPres34/density), p3-p4);
       Jac45 = -Jac35;
    } else {       /* collegamento inverso 1->3     2->4 */
       Jac15 = -dCoef*density*Cd*width*copysign(sqrt(2.*jumpPres13/density), p1-p3);
       Jac25 = -dCoef*density*Cd*width*copysign(sqrt(2.*jumpPres24/density), p2-p4);
       Jac35 = -Jac15;
       Jac45 = -Jac25;
    }
    
    doublereal Jac65  = -1;
    doublereal Jac66  = dCoef;  
 
    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, 5, Jac55);
    WM.PutCoef(5, 6, Jac56);
    WM.PutCoef(6, 5, Jac65);
    WM.PutCoef(6, 6, Jac66);
  
    return WorkMat;
 }

Here is the call graph for this function:

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

virtual

Implements Elem.

Definition at line 864 of file valve.cc.

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

 {
    DEBUGCOUT("Entering Dynamic_control_valve::AssRes()" << std::endl);
    
    WorkVec.Resize(6);
    
    integer iNode1RowIndex = pNode1->iGetFirstRowIndex()+1;
    integer iNode2RowIndex = pNode2->iGetFirstRowIndex()+1;
    integer iNode3RowIndex = pNode3->iGetFirstRowIndex()+1;
    integer iNode4RowIndex = pNode4->iGetFirstRowIndex()+1;
    integer iFirstIndex = iGetFirstIndex();
    
    doublereal p1 = pNode1->dGetX();
    doublereal p2 = pNode2->dGetX();
    doublereal p3 = pNode3->dGetX();
    doublereal p4 = pNode4->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 area_max = width*s_max;         // Area massima
    doublereal area_min = area_max*loss_area;  // Area minimo 
    doublereal deltaA = s*width;
 
    if (s > 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;
    }
 
     if (s <= -s_max) { 
       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.;
       }
 
    } else {
        c1 = 0.;
        c2 = 0.;
        c3 = 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.;
       }
       
    } else {
       cf1 = 0.;
       cf2 = 0.;
       cf3 = 0.;
    }
       
    
    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;
    
    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_5 = -Force+Mass*vp+C*sp+K*s+c1*(s+s_max)+c2*sp+c3*vp+cf1*(s-s_max)+cf2*sp+cf3*vp;
    doublereal Res_6 = sp-v;
 
    flow1 = -Res_1;
    flow2 = -Res_2;
    flow3 = -Res_3;
    flow4 = -Res_4;
 
 #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("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  eq.dinamica   : " << -Res_5 << std::endl); 
    DEBUGCOUT("-Res_6 sp-v           : " << -Res_6 << 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, iFirstIndex+1, Res_5);
    WorkVec.PutItem(6, iFirstIndex+2, Res_6);    
     
    return WorkVec;
 }

Here is the call graph for this function:

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

inlinevirtual

Reimplemented from Elem.

Definition at line 310 of file valve.h.

References pNode1, pNode2, pNode3, and pNode4.

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

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

virtual

Reimplemented from Elem.

Definition at line 697 of file valve.cc.

References ASSERT, and DofOrder::DIFFERENTIAL.

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

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

virtual

Implements HydraulicElem.

Definition at line 681 of file valve.cc.

References HydraulicElem::DYNAMIC_CONTROL_VALVE.

 {
    return HydraulicElem::DYNAMIC_CONTROL_VALVE;
 }

unsigned int Dynamic_control_valve::iGetNumDof ( void ) const

virtual

Reimplemented from Elem.

Definition at line 692 of file valve.cc.

 {
    return 2;
 }

void Dynamic_control_valve::Output ( OutputHandler & OH ) const

virtual

Reimplemented from ToBeOutput.

Definition at line 1042 of file valve.cc.

References A1, A2, A3, A4, ToBeOutput::bToBeOutput(), flow1, flow2, flow3, flow4, 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
         << " " << A1 << " "  << A2 << " "  << A3 << " " << A4 << " " << std::endl;
    }
 }

Here is the call graph for this function:

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

virtual

Implements Elem.

Definition at line 687 of file valve.cc.

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

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

virtual

Reimplemented from SimulationEntity.

Definition at line 1055 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 Dynamic_control_valve::WorkSpaceDim	(	integer *	piNumRows,
		integer *	piNumCols
	)		const

virtual

Implements Elem.

Definition at line 705 of file valve.cc.

 {
    *piNumRows = 6; 
    *piNumCols = 6; 
 }

Member Data Documentation

doublereal Dynamic_control_valve::A1

private

Definition at line 246 of file valve.h.

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

doublereal Dynamic_control_valve::A2

private

Definition at line 247 of file valve.h.

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

doublereal Dynamic_control_valve::A3

private

Definition at line 248 of file valve.h.

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

doublereal Dynamic_control_valve::A4

private

Definition at line 249 of file valve.h.

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

doublereal Dynamic_control_valve::c1

private

Definition at line 260 of file valve.h.

Referenced by AssJac(), and AssRes().

doublereal Dynamic_control_valve::c2

private

Definition at line 261 of file valve.h.

Referenced by AssJac(), and AssRes().

doublereal Dynamic_control_valve::c3

private

Definition at line 262 of file valve.h.

Referenced by AssJac(), and AssRes().

doublereal Dynamic_control_valve::c_acc

private

Definition at line 235 of file valve.h.

Referenced by AssRes().

doublereal Dynamic_control_valve::c_spost

private

Definition at line 233 of file valve.h.

Referenced by AssRes().

doublereal Dynamic_control_valve::c_vel

private

Definition at line 234 of file valve.h.

Referenced by AssRes().

doublereal Dynamic_control_valve::Cd

private

Definition at line 238 of file valve.h.

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

doublereal Dynamic_control_valve::cf1

private

Definition at line 263 of file valve.h.

Referenced by AssJac(), and AssRes().

doublereal Dynamic_control_valve::cf2

private

Definition at line 264 of file valve.h.

Referenced by AssJac(), and AssRes().

doublereal Dynamic_control_valve::cf3

private

Definition at line 265 of file valve.h.

Referenced by AssJac(), and AssRes().

doublereal Dynamic_control_valve::deltaP

private

Definition at line 266 of file valve.h.

Referenced by AssJac(), and AssRes().

doublereal Dynamic_control_valve::flow1

private

Definition at line 252 of file valve.h.

Referenced by AssRes(), and Output().

doublereal Dynamic_control_valve::flow2

private

Definition at line 253 of file valve.h.

Referenced by AssRes(), and Output().

doublereal Dynamic_control_valve::flow3

private

Definition at line 254 of file valve.h.

Referenced by AssRes(), and Output().

doublereal Dynamic_control_valve::flow4

private

Definition at line 255 of file valve.h.

Referenced by AssRes(), and Output().

doublereal Dynamic_control_valve::Force

private

Definition at line 236 of file valve.h.

doublereal Dynamic_control_valve::loss_area

private

Definition at line 240 of file valve.h.

Referenced by AssRes().

doublereal Dynamic_control_valve::Mass

private

Definition at line 251 of file valve.h.

const PressureNode* Dynamic_control_valve::pNode1

private

Definition at line 228 of file valve.h.

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

const PressureNode* Dynamic_control_valve::pNode2

private

Definition at line 229 of file valve.h.

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

const PressureNode* Dynamic_control_valve::pNode3

private

Definition at line 230 of file valve.h.

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

const PressureNode* Dynamic_control_valve::pNode4

private

Definition at line 231 of file valve.h.

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

doublereal Dynamic_control_valve::s

private

Definition at line 256 of file valve.h.

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

doublereal Dynamic_control_valve::s_max

private

Definition at line 244 of file valve.h.

Referenced by AssRes().

doublereal Dynamic_control_valve::sp

private

Definition at line 258 of file valve.h.

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

doublereal Dynamic_control_valve::start

private

Definition at line 232 of file valve.h.

Referenced by SetValue().

doublereal Dynamic_control_valve::v

private

Definition at line 257 of file valve.h.

Referenced by AssJac(), and AssRes().

doublereal Dynamic_control_valve::valve_density

private

Definition at line 242 of file valve.h.

Referenced by AssRes(), and Dynamic_control_valve().

doublereal Dynamic_control_valve::valve_diameter

private

Definition at line 241 of file valve.h.

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

doublereal Dynamic_control_valve::vp

private

Definition at line 259 of file valve.h.

Referenced by AssRes(), and Output().

doublereal Dynamic_control_valve::width

private

Definition at line 239 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