Ball Poppet with Conical Seat with Moving Body (IL)
Path in the library:
/Physical Modeling/Fluids/Isothermal Liquid/Valves & Orifices/Spools & Poppets/Moving Body/Ball Poppet with Conical Seat with Moving Body (IL)
Description
Block Ball Poppet with Conical Seat (IL) It is a one-dimensional movement of a spherical valve with a conical seat.
The resulting force acting on the valve is due to the pressure force and external forces. It is assumed that the pressure in port B acts on the active region adjacent to the hole and tends to open the hole. The pressure in port A acts on the remaining area of the ball. These assumptions give the pressure force acting on the ball. This force can be adjusted using hydrodynamic force.
The displacement and speed of the piston are supplied to the port RA. There are no restrictions on the displacement value in the block, but restrictions can be provided by an attached block using end stops (Translational Hard Stop).
If the check box is selected Moving body, then the block is implemented Ball Poppet with Conical Seat with Moving Body (IL) and the body movement is simulated. In this case, the displacement and velocity of the case are transmitted to the CA port. There are no restrictions on the displacement value in the block, but restrictions can be provided by an attached block using end stops.
Rise () is a variable related to the movement of the piston and the movement of the body, if it is simulated. Naturally, the limits of this rise are associated with a restriction on the values of movements.
If the ball’s lift is more than 20% of the ball’s diameter, then the accuracy will be reduced.
The opening area should never exceed the opening area of the neck, determined by the diameter of the seat and the diameter of the stem (on the side of the seat). Nevertheless, it is sometimes useful to limit the hole area to a minimum and/or maximum value. The minimum area can be used to simulate a leak or a special flow hole, even when the ball is fully attached to the seat. The maximum area can be used to simulate the flow area adjacent to the opening when the valve is wide open.
Please note that the flow rate is calculated based on the movement of the ball.
The equations
If the check box is Moving body removed, and the body movement is not simulated, then the ball lift is calculated as:
where
— the rise corresponding to the zero offset, the value of the parameter Lift corresponding to zero displacement;
— movement of the piston, which is inserted into the port RA.
If the check box is Moving body is installed, and the body movement is simulated, then the ball lift is calculated as:
where — moving the case, which is inserted into the CA port.
The equations that the block uses depend on the model. Flow force model:
Simple — a simple model of hydrodynamic force;
Corrected by effective pressure area factor — a model of hydrodynamic force adjusted for the effective pressure area coefficient, for more information, see [1-2];
Advanced active area — a model of hydrodynamic force with a modification of the calculation of the active area upstream, for more information, see [3].
Equations for a simple model
The minimum flow area is determined by the curved surface of the truncated cone, as shown in the figure. It is assumed that this surface divides the area occupied by the liquid into two areas with different pressures. There is pressure on one of these areas. , and on the other — pressure . This assumption is reasonable if the lift of the ball is small compared to the diameter of the saddle. If the rise of the ball is large, it is obvious that at some point the smallest limitation will be the area of the neck.
The condition must be met:
In case of violation of this condition, the ball will not be able to rest on the conical seat.
If this condition is violated, then it is considered that the saddle is round as in the block. Ball Poppet with Sharp Edge Seat (IL). Make sure that the condition is not violated.
The area of the hole is defined as:
where
, where — half of the solution angle of the conical seat;
— the diameter of the ball.
The hydraulic diameter is calculated as:
The active diameter is calculated as:
Note that the value used for , limited between and the smaller of the values and , where — this is the lifting value at which the calculated area becomes equal to the area of the neck:
Meaning Always more .
Usually the value is it is zero, but it can be set higher to simulate the leakage rate. Meaning it is usually very high (for example, Inf), but a much lower value can be set to simulate an additional hole.
Liquid volume , the pressure in which is equal to the pressure , additional to the volume when the valve is closed, is calculated as:
The value of the additional volume it is important when calculating pressure dynamics (frequency analysis).
The derivative of the additional volume by calculated as:
If the check box Moving body removed, the volume of liquid that is discharged into port B is calculated as:
where — parameter value Volume at port B corresponding to zero lift.
And the volume of liquid that is discharged into port A is calculated as:
where — parameter value Volume at port A corresponding to zero lift.
If the check box Moving body installed:
where — the diameter of the conical seat on the side of the ball, the value of the parameter Seat cylinder diameter.
Flow coefficient calculated as:
where
— pressure difference between ports;
— hydraulic diameter;
— kinematic viscosity;
— the average density of the liquid.
The average density is calculated at an average pressure .
Expense ratio calculated as:
where
— maximum flow rate, parameter value Maximum flow coefficient;
— critical flow coefficient, parameter value Critical flow number.
For meaning practically does not change. For low meaning it changes linearly with the change .
A reasonable value the default value is `1000'. However, for holes with complex (rough) geometry, it may be less than `50'. For very smooth geometry, it can be set to `50,000'.
The average fluid velocity is:
The volume consumption is:
where
— the area of the passage hole;
— the density of a liquid at atmospheric pressure.
If the check box Moving body if withdrawn, the volume costs in ports B and A are calculated as:
where
— liquid density at port pressure B, ;
— liquid density at port pressure A, ;
— the speed of the rod in the port RA.
If the check box Moving body installed:
where
— the diameter of the conical seat on the side of the ball, the value of the parameter Seat cylinder diameter;
— the speed of the case in the port C_s.
The hydrodynamic force is determined by estimating the change in momentum. This force tends to close the valve. For a steady-state fluid flow, the hydrodynamic force is:
where — the angle of the jet:
Dependence of the hydrodynamic force from the rise It is defined as follows:
where — parameter value Lift corresponding to minimum area.
The power in port RA is calculated as:
where — the power that enters the RB port.
If the check box is selected Moving body, and the hull movement is modeled, then the force in the port CA is calculated as:
where
— limited from above by the value ;
— the power that enters the CB port.
Equations for the model adjusted for the effective pressure area coefficient
Model Corrected by effective pressure area factor differs from Simple the fact that the value of the active diameter and the hydrodynamic force is calculated taking into account the type of flow (laminar or turbulent).
This section contains the equations for calculating the adjusted parameters, the remaining parameters are calculated in the same way as for the model. Simple.
Active diameter defined as:
Active area defined as:
Real active area depends on the type of flow (laminar or turbulent):
If the flow is laminar:
If the current is turbulent:
In these formulas — parameter value Turbulent effective pressure area factor.
The hydrodynamic force also varies depending on the velocity of the liquid at the inlet and outlet of the hole.:
where
— the velocity of the liquid at the entrance to the hole:
— the velocity of the liquid at the outlet of the hole:
The power in port RA is calculated as:
where — the power that enters the RB port.
If the check box is selected Moving body, and the hull movement is modeled, then the force at the port CA is calculated as:
where
— limited from above by the value , calculated from the stem value[A_P]:
— the power that enters the CB port.
Equations with modification of calculation of active areas upstream
Model Advanced active area differs from Simple by calculating the active areas (minimum and maximum) upstream.
This section contains the equations for calculating the adjusted parameters, the remaining parameters are calculated in the same way as for the model. Simple.
The maximum active area upstream is calculated as:
The minimum active area upstream is calculated as follows:
The active area upstream is calculated as:
It should be noted that the calculation performed is valid only for the incoming jet and is not valid for the outgoing jet. Therefore, the updraft must be on the side of the saddle. Then the power in the port RA is calculated as:
where — the active diameter calculated from the value of the active area upstream .
If the check box is selected Moving body, and the hull movement is modeled, then the force at the port CA is calculated as:
A mechanical translational port corresponding to the body from the side opposite to the side of the seat.
Dependencies
To use this port, check the box Moving body.
Program usage name
case_flange_b
Parameters
Parameters
#Opening orientation —
the direction of movement of the spool corresponding to the opening of the hole
Positive relative rod displacement opens orifice | Negative relative rod displacement opens orifice
Details
The direction of movement of the element corresponding to the opening of the hole:
Positive orientation Positive relative rod displacement opens orifice this means that the positive movement of the spool opens the hole.
The negative direction Negative relative rod displacement opens orifice this means that the negative movement of the spool opens the hole.
Values
Positive relative rod displacement opens orifice | Negative relative rod displacement opens orifice
Check this box if you are modeling a movable enclosure.
If the flag is unchecked, it is assumed that the body is stationary.
Default value
—
Program usage name
moving_case
Evaluatable
No
#Flow force model —
the model of hydrodynamic force
Simple | Corrected by effective pressure area factor | Advanced active area
Details
The model of hydrodynamic force:
Simple — a simple model of hydrodynamic force;
Corrected by effective pressure area factor — a model of hydrodynamic force adjusted for the effective pressure area coefficient;
Advanced active area — a model of hydrodynamic force with a modification of the calculation of the active area upstream.
Values
Simple | Corrected by effective pressure area factor | Advanced active area
Default value
Simple
Program usage name
flow_force_model
Evaluatable
No
#Seat cylinder diameter —
diameter of the conical seat on the side of the ball
m | cm | ft | in | km | mi | mm | um | yd
Details
Diameter of the conical seat on the side of the ball, .
Dependencies
To use this option, check the box Moving body.
Units
m | cm | ft | in | km | mi | mm | um | yd
Default value
10.0 mm
Program usage name
seat_cylinder_diameter
Evaluatable
Yes
#Seat diameter (hole) —
diameter of the conical seat at the base
m | cm | ft | in | km | mi | mm | um | yd
Details
Seat diameter, .
Units
m | cm | ft | in | km | mi | mm | um | yd
Default value
5.0 mm
Program usage name
seat_diameter
Evaluatable
Yes
#Seat semi-angle (between 0 and 90) —
half of the solution angle of the conical seat
deg | rad | rev | mrad
Details
Half of the solution angle of the conical seat, .
Units
deg | rad | rev | mrad
Default value
45.0 deg
Program usage name
conical_seat_semi_angle
Evaluatable
Yes
#Ball diameter —
diameter of the ball
m | cm | ft | in | km | mi | mm | um | yd
Details
Diameter of the ball, .
Units
m | cm | ft | in | km | mi | mm | um | yd
Default value
10.0 mm
Program usage name
ball_diameter
Evaluatable
Yes
#Rod diameter (opposite to seat) —
the diameter of the stem on the side opposite to the seat
m | cm | ft | in | km | mi | mm | um | yd
Details
The diameter of the stem on the side opposite to the seat, .
Units
m | cm | ft | in | km | mi | mm | um | yd
Default value
0.0 mm
Program usage name
rod_diameter_at_seat_opposite_side
Evaluatable
Yes
#Rod diameter (seat side) —
diameter of the rod on the seat side
m | cm | ft | in | km | mi | mm | um | yd
Details
Diameter of the rod on the seat side, .
Units
m | cm | ft | in | km | mi | mm | um | yd
Default value
0.0 mm
Program usage name
rod_diameter_at_seat_side
Evaluatable
Yes
#Turbulent effective pressure area factor —
effective pressure area coefficient
Details
Effective pressure area coefficient, .
Dependencies
To use this parameter, set for the parameter Flow force model meaning Corrected by effective pressure area factor.
Default value
1.0
Program usage name
effective_pressure_area_factor
Evaluatable
Yes
#Lift corresponding to zero displacement —
the rise corresponding to the zero offset
m | cm | ft | in | km | mi | mm | um | yd
Details
The rise corresponding to the zero offset.
Units
m | cm | ft | in | km | mi | mm | um | yd
Default value
0.0 mm
Program usage name
lift_offset
Evaluatable
Yes
#Lift corresponding to minimum area —
lifting corresponding to the minimum area
m | cm | ft | in | km | mi | mm | um | yd
Details
Lift , corresponding to the minimum area of the passage hole.
Units
m | cm | ft | in | km | mi | mm | um | yd
Default value
0.0 mm
Program usage name
orifice_opening_at_min_area
Evaluatable
Yes
#Lift corresponding to maximum area —
lifting corresponding to the maximum area
m | cm | ft | in | km | mi | mm | um | yd
Details
Lift , corresponding to the maximum area of the passage hole.
Units
m | cm | ft | in | km | mi | mm | um | yd
Default value
Inf mm
Program usage name
orifice_opening_at_max_area
Evaluatable
Yes
#Volume at port A corresponding to zero lift —
the volume in port A corresponding to zero lift
l | gal | igal | m^3 | cm^3 | ft^3 | in^3 | km^3 | mi^3 | mm^3 | um^3 | yd^3 | N*m/Pa | N*m/bar | lbf*ft/psi | ft*lbf/psi
#Volume at port B corresponding to zero lift —
the volume in port B corresponding to zero lift
l | gal | igal | m^3 | cm^3 | ft^3 | in^3 | km^3 | mi^3 | mm^3 | um^3 | yd^3 | N*m/Pa | N*m/bar | lbf*ft/psi | ft*lbf/psi
#Jet force coefficient —
coefficient of hydrodynamic force
Details
The hydrodynamic force coefficient, which disables the hydrodynamic force when set to 0 (by default), and turns it on when set to 1. If there is experimental data for this coefficient, then you can adjust the model to fit this data.
The maximum flow rate affects the flow rate/pressure drop characteristics in the orifice. For most applications, this value can be left at the default.
The critical flow coefficient affects the flow rate/pressure drop characteristics in the orifice. For most applications, this value can be left at the default.
Default value
100.0
Program usage name
critical_flow_number
Evaluatable
Yes
Initial Conditions
#Initial rod displacement —
initial displacement of the rod
m | cm | ft | in | km | mi | mm | um | yd
Details
The initial displacement of the rod.
Units
m | cm | ft | in | km | mi | mm | um | yd
Default value
0.0 mm
Program usage name
rod_displacement_start
Evaluatable
Yes
#Initial case displacement —
initial displacement of the body
m | cm | ft | in | km | mi | mm | um | yd
Details
The initial displacement of the body.
Dependencies
To use this option, check the box Moving body.
Units
m | cm | ft | in | km | mi | mm | um | yd
Default value
0.0 mm
Program usage name
case_displacement_start
Evaluatable
Yes
Literature
N.Mittwollen, T.Michl, R.Breit "Parametric hydraulic valve model including transitional flow effects", 2nd MATHMOD Vienna 1997 (IMACS).
N.Mittwollen, "Hydraulic simulation of cavitation induced pressure fluctuations with peculiar periodicities in a fluid power unit", 8th Bath international fluid power workshop, September 1995
A. Clavier, M. Alirand, F. Vernat, B. Sagot, "Local approach to improve the global approach of hydraulic forces in ball poppet valves", 4 th Int. Symposium on Fluid Power, Wuhan, China, April 2003.