Sw. Attaway et al., COUPLING OF SMOOTH PARTICLE HYDRODYNAMICS WITH THE FINITE-ELEMENT METHOD, Nuclear Engineering and Design, 150(2-3), 1994, pp. 199-205
A gridless technique called smooth particle hydrodynamics (SPH) has be
en coupled with the transient dynamics finite element code PRONTO. In
this paper, a new weighted residual derivation for the SPH method will
be presented, and the methods used to embed SPH within PRONTO will be
outlined. Example SPH PRONTO calculations will also be presented. One
major difficulty associated with the Lagrangian finite element method
is modeling materials with no shear strength; for example, gases, flu
ids and explosive biproducts. Typically, these materials can be modele
d for only a short time with a Lagrangian finite element code. Large d
istortions cause tangling of the mesh, which will eventually lead to n
umerical difficulties, such as negative element area or ''bow tie'' el
ements. Remeshing will allow the problem to continue for a short while
, but the large distortions can prevent a complete analysis. SPH is a
gridless Lagrangian technique. Requiring no mesh, SPH has the potentia
l to model material fracture, large shear flows and penetration. SPH c
omputes the strain rate and the stress divergence based on the nearest
neighbors of a particle, which are determined using an efficient part
icle-sorting technique. Embedding the SPH method within PRONTO allows
part of the problem to be modeled with quadrilateral finite elements,
while other parts are modeled with the gridless SPH method. SPH elemen
ts are coupled to the quadrilateral elements through a contact-like al
gorithm.