Cylinder Components

Cylinder Component General Description

The pressure created by a fed working fluid is the basis for the operation of a hydraulic cylinder. Flow Simulator provides two types of piston cylinders. In a single-acting hydraulic cylinder, the working fluid comes into the cylinder from one hydraulic tank and presses down on the piston from one side, creating a return stroke of the rod. The return spring actuates a forward stroke.

In a double-acting cylinder, the return spring is absent, its role is performed by the second hydraulic tank, creating pressure on the piston from the other side and, as a result, providing reciprocating motions of the mechanism.

Quick Guide for Cylinder Component Creation in the GUI

There are two Cylinder subtypes under “Compressible Gas Elements” and “Incompressible Liquid Elements” sections in Flow Simulator.

Cylinder Component Inputs

Table of the inputs for the Cylinder Component.

Cylinder Component Input Variables
Index UI Name (. flo label) Description
1 Type (CPTYPE) Type of component. See common component inputs table above for a full list.
2 Subtype (SUBTYPE)

Subtype of piston cylinder component;

  1. Single: having a control rod on one side
  2. Double: having working fluid on both sides
3 Cylinder Diameter (CYL_DIAM) The diameter of the cylinder base face
4 Side 1 Rod Diameter (ROD_DIAM1) Rod Diameter of 1st slide
6 Maximum Rod Displacement (MAX_ROD_TRAVEL) The maximum displacement that the rod is allowed to travel by
8 Piston and Rod Mass (PISTON_ROD_MASS) Total mass of the piston and rod
9 Leakage Coefficient (LEAK_COEF)

The coefficient which describes the level of leakage between the sides of the cylinder base and the walls of the element. The leakage flow is proportional to the pressure difference across the piston. For more advanced leakage modeling, you can use a separate flow element that connects the one side of the cylinder and the other.

Note 1: If there are different materials on either side of the pistons, the leakage coefficient should be 0, because for example, oil should not leak into gas. For mixtures, please use a Custom Material definition.

Note 2: Only for Double Acting Cylinder piston

10 Reference Viscosity (MU_REF) Viscosity value to calculate viscous force and viscous force is added in total drag force
6 Columb Friction Coefficient (COULOMB_FRIC) Columb type of friction coefficient to be added in total drag force
7 (STICTION_FORCE) Stiction force value to be added in total drag force
8 (STICTION_VEL)
10 Spring #1 Force at Position=0 (F0_SPRING1) The initial, preloaded force on the spring at position 0
11 Spring #1 Constant (K_SPRING1) Spring constant, k of 1st spring to be used in total spring force
12 Spring #2 Force at Position=0 (F0_SPRING2) The initial, preloaded force on the spring at position 0
13 Spring #2 Constant (K_SPRING2) Spring constant, k of 2nd spring to be used in total spring force
14 Initial Position (ROD_POSITION0) The starting position of the cylinder base face
15 Slide #1 External Pressure on Rod (PS_EXTERNAL1) The external pressure amount on the rod
16 (F_EXTERNAL1) The initial, preloaded force on the spring at position L (Maximum Rod Displacement)
17 Slide #2 External Pressure on Rod (PS_EXTERNAL2) External pressure on the rod
18 (F_EXTERNAL2) The initial, preloaded force on the spring at position L (Maximum Rod Displacement)
19 Initial Velocity (ROD_VELOCITY0) The initial velocity of the cylinder rod at the beginning of the transient run
20 Slide #1 Fluid Volume at Position =0 (VOL1_AT_POS0) Fluid volume when the piston is at initial position
21 Slide #1 Fluid Volume at Position =0 (VOL2_AT_POS0) Fluid volume when the piston is at initial position
22 Slide #2 Rod Diameter (ROD_DIAM2) Slide 2nd rod diameter
23 Spring #1 Constant (DISCRETIZATION) Spring constant, k of 1st spring

Cylinder Component Theory Manual

Cylinder Governing Equations

Piston-Rod Force Variable Definitions

Force Balance Equation (Newton’s 2nd Law)

The Force Balance Equation can be re-written as a 2nd Order linear ODE

There are two complications.

  1. The Coulomb (Dynamic) Friction term is either positive or negative depending on direction of piston motion. This means the governing ODE of piston motion needs to be split into two parts, one for left-to-right motion and the other for right-to-left motion. In both cases, the initial position and velocity are the initial conditions.

  2. The total pressure force and external force () is changing within each timestep. Floinhance assumes linear variation within each timestep, which is 2nd order accurate, and helps when you want to take relatively large fluid timesteps. However, it should be noted that you can always get achieve greater accuracy by taking smaller fluid timesteps.

The governing ODE can be split into four distinct ODE regimes by introducing and parameters. (LR signifies left-to-right motion and RL is right-to-left motion.)

Depending on the ODE coefficients of the two non-trivial ODEs, the ODE solution can follow one of eight different equation forms.

Under-damped when

Solution Parameters based on ODE Coefficients and Initial Conditions:

Solution Form:

Over-damped when

Solution Parameters based on ODE Coefficients and Initial Conditions:

Critically Damped when

Solution Parameters based on ODE Coefficients and Initial Conditions:

Spring-less when K=0

Solution Parameters based on ODE Coefficients and Initial Conditions:

Spring-less & Friction-less when K=0 & C=0

Solution Parameters based on ODE Coefficients and Initial Conditions:

Mass-less when M=0

Solution Parameters based on ODE Coefficients and Initial Conditions:

Mass-less & Friction-less when M=0 & C=0

Solution Parameters based on ODE Coefficients and Initial Conditions:

Mass-less & Spring-less when M=0 and K=0

Solution Parameters based on ODE Coefficients and Initial Conditions:

Adaptive Time-stepping

Detecting changes in four ODE regimes adaptively and adjusting the piston motion timestep and/or fluid timestep is critical to accurate results.

If the “Mechanical Only” adaptive option is chosen, then is equal to . The piston-cylinder solver will shorten the piston motion ODE interval. A second (or sometimes more) piston motion solution intervals will be required to reach . Therefore, there will be piecewise piston motion solution in on the interval from . However, the linearly varying force assumption, which is usually due to time-varying pressures, still applies on the interval from . This means you a not resolving the fluid pressure changes at the shorter timestep.

If the “Coarse” adaptive option is chosen, then is set equal to , so the fluid timestep is shortened as well. This achieves greater accuracy.

The “Fine” adaptive option is like “Coarse”, expect with an addition constraint. If there is a spring, then there is a natural oscillation period based on the spring stiffness and piston-rod mass. The timestep will not be allowed to go higher than 1/N of the natural oscillation period. The default N is 12. You can change this input in the Advanced Options menu. This option should be used for Cylinder components with high coupling sensitivity to the fluid model – and where oscillations of the piston must be captured to correctly model the physical behavior of the system.

Finally, the “Off” option does not shorten either mechanical or fluid timestep. However, it does use a numerical trick to prevent spurious oscillations that occur when you require the piston to leave the end at the “wrong time”. It uses a smoothing option, described in detail in the following section, but only for the time-steps where the piston becomes unstuck.

Smooth Solution Option

The Cylinder Component provides an option to smooth the solution of oscillating pistons. Usually the fluid network will dampen the piston motion. Alternatively, you can specify viscous friction. However, in some cases, these two options may not be sufficient. If you do not know how much friction or fluid mass that is required for damping, but you know the piston motion should be smooth (non-oscillatory), then you can use the Smooth option. If the solution is in the underdamped mode, the smooth option will solve for piston motion including oscillations, and only at the last stage it zeros out the amplitude of the sine and cosine terms. For piston position, this is not usually important. However, for piston velocity, which directly affects cylinder port flow rates that are returned to the fluid solver, this can provide critical stabilization. This will aid both convergence robustness and solution accuracy. For all the other ODE regimes besides underdamped, the Smooth option does not take effect.

Mechanical & Fluid Coupling Method

To couple the piston motion with the fluid solver, the mass flow rate at the piston ports must be calculated. Here we define as the rate of mass entering the cylinder through the port.

The velocity term needs special treatment that is analogous to unwinding in CFD. This helps to achieve a stable and accurate solution.

In other words, Flow Simulator uses the exact solution for velocity to calculate mass flow rate when the piston is not near the end positions. However, when the piston becomes stuck or unstuck from the end position, Flow Simulator uses a conservative differencing formula. This makes sure the rest of fluid circuit experiences the correct amount of mass release or mass removal by the cylinder.

Cylinder Component Outputs

The following listing provides details about Cylinder Component output variables.

Note: Currently, Cylinder Piston Component results are only in ENG unit.

Name Description Units
MAIN SOLVER ITER CVG? Information whether the main solver converged or not and if converged at which iteration Flag
PISTON POSITION Position of the piston at each time step in
PISTON VELOCITY Velocity of the piston ft/s
PISTON ACCEL Acceleration of the piston during the motion ft/s2
TOTAL LOAD FORCE Total load force Lbf
TOTAL PRESSURE FORCE Total pressure force Lbf
TOTAL SPRING FORCE Total spring force (F_Spring + K_Spring*Distance) Lbf
VISCOUS DRAG FORCE Viscous drag force Lbf
COULOMB DRAG FORCE Coulomb drag force (coulomb_fric*pressure_force) Lbf
STICTION DRAG FORCE Stiction drag force Lbf
CYLINDER END-WALL FORCE Force on the cylinder end-wall Lbf
NOT ROD FORCE Rod force based on external pressure on rod Lbf
INNER PS (SIDE 1 & SIDE 2) Static Pressure inside psi
ACCUM (SIDE 1 & SIDE 2) Accumulated mass Lbm
ACCUM RATE (SIDE 1 & SIDE 2) Accumulation rate pps
LEAKAGE UPTAKE (SIDE 1 & SIDE 2) Leakage value pps
EXTERNAL UPTAKE (SIDE 1 & SIDE 2) External force pps
DYNVISC (SIDE 1 & SIDE 2) Dynamic viscosity of the fluid Lbm/ft/hr
PRESSURE FORCE (SIDE 1 & SIDE 2) Side 1’s pressure force Lbf
SPRING FORCE (SIDE 1 & SIDE 2) Side 1’s spring force Lbf