Documentation |
Variable-displacement reversible hydraulic machine with regime-dependable efficiency
The Variable-Displacement Hydraulic Machine block represents a variable-displacement hydraulic machine of any type as a data-sheet-based model. The model accounts for the power flow direction and simulates the machine in both the motor and pump mode. The efficiency of the machine is variable, and you can set it in accordance with experimental data provided in the catalog or data sheet.
The machine displacement is controlled by the signal provided through the physical signal port C. The machine efficiency is simulated by implementing regime-dependable leakage and friction torque based on the experimentally established correlations between the machine efficiencies and pressure, angular velocity, and displacement.
With respect to the relationship between the control signal and the displacement, two block parameterization options are available:
By the maximum displacement and stroke — The displacement is assumed to be linearly dependent on the control member position.
By table-specified relationship between the control member position and the machine displacement — The displacement is determined by one-dimensional table lookup based on the control member position. You have a choice of three interpolation methods and two extrapolation methods.
The variable-displacement machine is represented with the following equations:
$$q=D\xb7\omega -{k}_{m}\xb7{q}_{L}$$
$$T=D\xb7p+{k}_{m}\xb7{T}_{fr}$$
$$D=\{\begin{array}{l}\frac{{D}_{\mathrm{max}}}{{x}_{\mathrm{max}}}\cdot x\hfill \\ D(x)\text{}\hfill \end{array}$$
$$p={p}_{A}-{p}_{B}$$
where
q | Machine flow rate |
p | Pressure differential across the machine |
p_{A,}p_{B} | Gauge pressures at the block terminals |
D | Machine instantaneous displacement |
D_{max} | Machine maximum displacement |
x | Control member displacement |
x_{max} | Control member maximum stroke |
T | Torque at the machine shaft |
ω | Machine shaft angular velocity |
q_{L} | Leakage flow |
T_{fr} | Friction torque |
k_{m} | Machine type coefficient. k_{m} = 1 for the pump, k_{m} = –1 for the motor. |
The key parameters that determine machine efficiency are its leakage and friction on the shaft. In the block, these parameters are specified with experimentally-based correlations similar to [1]
$${q}_{L}=D\xb7\omega \xb7{k}_{L1}{\left(\frac{p}{{p}_{nom}}\right)}^{{k}_{LP}}{\left(\frac{D}{{D}_{\mathrm{max}}}\right)}^{{k}_{LD}}{\left(\frac{\omega}{{\omega}_{nom}}\right)}^{{k}_{L\omega}}$$
$${T}_{fr}=D\xb7p\xb7{k}_{F1}{\left(\frac{p}{{p}_{nom}}\right)}^{{k}_{FP}}{\left(\frac{D}{{D}_{\mathrm{max}}}\right)}^{{k}_{FD}}{\left(\frac{\omega}{{\omega}_{nom}}\right)}^{{k}_{F\omega}}$$
where
p_{nom} | Nominal pressure |
ω_{nom} | Nominal angular velocity |
k_{L1} | Leakage proportionality coefficient |
k_{F1} | Friction proportionality coefficient |
k_{LP}, k_{LD}, k_{Lω}, k_{FP}, k_{FD}, k_{Fω} | Approximating coefficients |
The approximating coefficients are determined from the efficiency plots, usually provided by the machine manufacturer. With the leakage known, the pump volumetric efficiency can be expressed as
$${\eta}_{vp}=\frac{D\cdot \omega -{q}_{L}}{D\cdot \omega}=1-{k}_{L1}{\left(\frac{p}{{p}_{nom}}\right)}^{{k}_{LP}}{\left(\frac{D}{{D}_{\mathrm{max}}}\right)}^{{k}_{LD}}{\left(\frac{\omega}{{\omega}_{nom}}\right)}^{{k}_{L\omega}}$$
For a motor, the expression looks like the following
$${\eta}_{vm}=\frac{D\cdot \omega}{D\cdot \omega +{q}_{L}}=\frac{1}{1+{k}_{L1}{\left(\frac{p}{{p}_{nom}}\right)}^{{k}_{LP}}{\left(\frac{D}{{D}_{\mathrm{max}}}\right)}^{{k}_{LD}}{\left(\frac{\omega}{{\omega}_{nom}}\right)}^{{k}_{L\omega}}}$$
The mechanical efficiency is based on the known friction torque
$${\eta}_{mp}=\frac{D\cdot p}{D\cdot p+{T}_{fr}}=\frac{1}{1+{k}_{F1}{\left(\frac{p}{{p}_{nom}}\right)}^{{k}_{FP}}{\left(\frac{D}{{D}_{\mathrm{max}}}\right)}^{{k}_{FD}}{\left(\frac{\omega}{{\omega}_{nom}}\right)}^{{k}_{F\omega}}}$$
$${\eta}_{mm}=\frac{D\cdot p-{T}_{fr}}{D\cdot p}=1-{k}_{F1}{\left(\frac{p}{{p}_{nom}}\right)}^{{k}_{FP}}{\left(\frac{D}{{D}_{\mathrm{max}}}\right)}^{{k}_{FD}}{\left(\frac{\omega}{{\omega}_{nom}}\right)}^{{k}_{F\omega}}$$
The curve-fitting procedure is based on the comparison of the efficiency, determined with one of the above expressions, and the experimental data $${\eta}_{\mathrm{exp}}=f(p,D,\omega )$$, an example of which is shown in the following plot.
The procedure can be performed with the Optimization Toolbox™ software. For instance, the pump volumetric efficiency approximating coefficients can be found by solving the following problem:
$$\underset{x}{\mathrm{min}}F(x)$$
$$x=\left[{k}_{L1},{k}_{LP},{k}_{LD},{k}_{L\omega}\right]$$
$$F(x)={{\displaystyle \sum _{i}{\displaystyle \sum _{j}{\displaystyle \sum _{k}\left({\eta}_{\mathrm{exp}}\left({p}_{i},{D}_{j},{\omega}_{k}\right)-\left(1-{k}_{L1}{\left(\frac{{p}_{i}}{{p}_{nom}}\right)}^{{k}_{LP}}{\left(\frac{{D}_{j}}{{D}_{\mathrm{max}}}\right)}^{{k}_{LD}}{\left(\frac{{\omega}_{k}}{{\omega}_{nom}}\right)}^{{k}_{L\omega}}\right)\right)}}}}^{2}$$
where
i | Number of experimental pressure points, from 1 to n |
j | Number of experimental displacement points, from 1 to m |
k | Number of experimental angular velocity points, from 1 to l |
Connections A and B are hydraulic conserving ports associated with the machine inlet and outlet, respectively. Connection S is a mechanical rotational conserving port associated with the machine shaft. Connection C is a physical signal port that controls machine displacement. The flow rate from port A to port B causes the shaft to rotate in positive direction, provided positive signal is applied to port C.
Fluid compressibility is neglected.
No inertia on the machine shaft is considered.
The model is applicable only for fluid and fluid temperature at which the approximating coefficients have been determined.
Exercise extreme caution to not exceed the limits within which the approximating coefficients have been determined. The extrapolation could result in large errors.
Select one of the following block parameterization options:
By maximum displacement and control member stroke — Provide values for maximum machine displacement and maximum stroke. The displacement is assumed to be linearly dependent on the control member position. This is the default method.
By displacement vs. control member position table — Provide tabulated data of machine displacements and control member positions. The displacement is determined by one-dimensional table lookup. You have a choice of three interpolation methods and two extrapolation methods.
Machine maximum displacement. The default value is 5e-6 m^3/rad.
Maximum control member stroke. The default value is 0.005 m. This parameter is used if displacement is specified as By maximum displacement and control member stroke.
Specify the vector of machine displacements as a one-dimensional array. The machine displacements vector must be of the same size as the control member positions vector. The default values, in m^3/rad, are [-5e-06 -3e-06 0 3e-06 5e-06]. This parameter is used if displacement is specified as By displacement vs. control member position table.
Specify the vector of input values for control member position as a one-dimensional array. The input values vector must be strictly increasing. The values can be nonuniformly spaced. The minimum number of values depends on the interpolation method: you must provide at least two values for linear interpolation, at least three values for cubic or spline interpolation. The default values, in meters, are [-0.0075 -0.0025 0 0.0025 0.0075]. This parameter is used if displacement is specified as By displacement vs. control member position table.
Select one of the following interpolation methods for approximating the output value when the input value is between two consecutive grid points:
Linear — Uses a linear interpolation function.
Cubic — Uses the Piecewise Cubic Hermite Interpolation Polinomial (PCHIP).
Spline — Uses the cubic spline interpolation algorithm.
For more information on interpolation algorithms, see the PS Lookup Table (1D) block reference page. This parameter is used if displacement is specified as By displacement vs. control member position table.
Select one of the following extrapolation methods for determining the output value when the input value is outside the range specified in the argument list:
From last 2 points — Extrapolates using the linear method (regardless of the interpolation method specified), based on the last two output values at the appropriate end of the range. That is, the block uses the first and second specified output values if the input value is below the specified range, and the two last specified output values if the input value is above the specified range.
From last point — Uses the last specified output value at the appropriate end of the range. That is, the block uses the last specified output value for all input values greater than the last specified input argument, and the first specified output value for all input values less than the first specified input argument.
For more information on extrapolation algorithms, see the PS Lookup Table (1D) block reference page. This parameter is used if displacement is specified as By displacement vs. control member position table.
Nominal pressure differential across the machine. The default value is 1e7 Pa.
Nominal angular velocity of the output shaft. The default value is 188 rad/s.
The friction torque on the machine shaft ideally should be introduced as T_{fr}sign(ω). To avoid discontinuity at ω –> 0, the friction is defined as T_{fr}tanh(4ω /ω_{max} ), where ω_{max} is a small velocity, representing the shaft velocity at peak friction, at which tanh(4ω /ω_{max} ) is equal to 0.999. The default value of ω_{max} is 0.01 rad/s.
Approximating coefficient k_{L1} in the block description preceding. The default value is 0.05.
Approximating coefficient k_{LP} in the block description preceding. The default value is 0.65.
Approximating coefficient k_{Lω} in the block description preceding. The default value is -0.2.
Approximating coefficient k_{LD} in the block description preceding. The default value is -0.8.
Approximating coefficient k_{F1} in the block description preceding. The default value is 0.06.
Approximating coefficient k_{FP} in the block description preceding. The default value is -0.65.
Approximating coefficient k_{Fω} in the block description preceding. The default value is 0.2.
Approximating coefficient k_{FD} in the block description preceding. The default value is -0.75.