15.4. Example Model Using Hydrostatic Fluid Elements

The following example shows an axisymmetric model of a cylinder-piston arrangement. The contained fluid in the cylinder is modeled with hydrostatic fluid elements, shown as triangular elements in the figure below.

Figure 15.2: Meshed Model

Meshed Model

Due to the gap between the cylinder and piston, a hydrostatic fluid element with no underlying solid elements is added manually. The node numbering on this hydrostatic fluid element is reversed to achieve a negative fluid volume that balances the extra fluid volume created by hydrostatic fluid elements attached to the cylinder above the piston. The figure below shows a close-up view of the node numbering in the area where the hydrostatic fluid element with negative volume and no underlying solid element is defined.

Figure 15.3: Close-up View of Meshed Model

Close-up View of Meshed Model

The following sequence of commands generates this model, including the hydrostatic fluid elements, material models, and boundary conditions.

/batch,list
/title, Example input data for defining hydrostatic fluid elements
/prep7

!
! Define element types
!

et,1,182	!plane182
keyopt,1,3,1	!axisymmetric option
et,2,241	!hsfld241 
keyopt,2,3,1
!
! Define material type
!

! Material properties for Steel
mp,ex,1,2.0e5	!Young's Modulus in N/mm^2
mp,nuxy,1,0.3	!Poisson's ratio
mp,dens,1,7.7e-9!Density in ton/mm^3

! Material properties for hydraulic fluid
tb,fluid,2,,,,liquid
tbdata,1,2000.0	!Bulk Modulus in N/mm^2
tbdata,2,6.4e-4	!Coefficient of thermal expansion
tbdata,3,8.42e-5!Density in ton/mm^3

!
! Define nodes, keypoints, areas
!

! Nodes
n,1,0.0,50.0	!pressure node

! Keypoints
k, 1, 0.0,  0.0
k, 2,25.0,  0.0
k, 3,25.0,100.0
k, 4,10.0,100.0
k, 5,10.0, 95.0
k, 6,20.0, 95.0
k, 7,20.0,  5.0
k, 8, 0.0,  5.0
k, 9, 0.0, 80.0
k,10,19.0, 80.0
k,11,19.0, 85.0
k,12, 2.5, 85.0
k,13, 2.5,185.0
k,14, 0.0,185.0

! Areas
a,1,2,3,4,5,6,7,8
a,9,10,11,12,13,14

!
! Define mesh
!

! Planar elements
type,1
mat,1
esize,2.5
amesh,all

! Hydrostatic fluid elements--automatic generation
type,2
mat,2
lsel,s,line,,6,9
nsll,s,1
esurf,1
allsel,all

! Hydrostatic fluid elements--E command
e,61,162,,1	!element with negative volume--nodes 61 and 162 swapped

!
! Define degree-of-freedom constraints (boundary conditions)
!

! Symmetry boundary condition on pressure node
d,1,ux,0.0

! Symmetry boundary condition on planar mesh
lsel,s,line,,8,14,6
nsll,s,1
d,all,ux,0.0
allsel,all

! Fixed boundary condition on planar mesh
lsel,s,line,,1
nsll,s,1
d,all,all,0.0
allsel,all

finish