SPH Cells Distribution

The particles should be created in a hexagonal compact, face centered cubic or a cubic net packing.

The hexagonal compact net and face centered cubic are recommended for use in Radioss and give similar results. A face center cubic particle distribution can be created using the HyperMesh SPH panel. A HyperMesh Tcl macro to generate hexagonal compact net is available from the Altair Connect website.

Hexagonal Compact Net

A cubic centered faces net realizes a hexagonal compact distribution, which can be useful to build the net.


Figure 1. Local Views of the Hexagonal Compact Net


Figure 2. Perspective View of the Cubic Centered Faces Net

The hexagonal compact net distribution can be created in HyperMesh using a Tcl script available by searching in Altair Connect for the Hexagonal Compact Net Tcl script. When using this script, the pitch or distance between any particle and its closest neighbor is entered as h0.

The mass of the particle mp is defined in the property /PROP/SPH.

The SPH particle mass relates to the material density ρ and the pitch h0 of the hexagonal compact net. This particle mass can be calculated as:(1)
m p h 0 3 2 ρ
Since the space can be partitioned into polyhedras surrounding each particle of the net, each one with a volume:(2)
V p = h 0 3 2
Due to discretization differences in the volume, the mass can be more accurately represented by:(3)
m P = ρV n MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacaWGTbWaaSbaaSqaaiaadcfaaeqaaOGaeyypa0ZaaSaaaeaacqaH bpGCcaWGwbaabaGaamOBaaaaaaa@3F79@
Where,
V MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacaWGwbaaaa@39B3@
Total volume filled by the particles.
n MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacaWGUbaaaa@39CB@
Total number of particles distributed in the volume.

For hexagonal compact net, the recommended smoothing length h in /PROP/SPH is the pitch h0 which is the smallest distance between the particles. A smoothing length smaller than is this can only be used when there is no tension physical problems material. If the material does include tensile behavior, then a smoothing length larger than h0 can be used to increase stability but there will be an increase in the computational cost.

The following table shows the number of neighbors with when different smoothing lengths h in /PROP/SPH are used. The accuracy and computational cost of the simulation improves as the smooth length increases.
Table 1. Number of Neighbors in a Hexagonal Compact Net
Distance d Number of particles at distance d Number of particles within distance d
h0 12 12
2 h 0 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aadaGcaaqaaiaaikdaaSqabaGccaWGObWaaSbaaSqaaiaaicdaaeqa aaaa@3B8C@ 6 18
3 h 0 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aadaGcaaqaaiaaiodaaSqabaGccaWGObWaaSbaaSqaaiaaicdaaeqa aaaa@3B8D@ 24 42
2h0 12 54
5 h 0 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aadaGcaaqaaiaaiwdaaSqabaGccaWGObWaaSbaaSqaaiaaicdaaeqa aaaa@3B8F@ 24 78

Face Centered Cubic

The face centered cubic arranges the particles in groups of 14, forming the corners and the center of each face of a cube.


Figure 3.
Similar to hexagonal compact net, each particle has 12 neighbors and the mass of a particle is:(4)
m P = ρV n MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacaWGTbWaaSbaaSqaaiaadcfaaeqaaOGaeyypa0ZaaSaaaeaacqaH bpGCcaWGwbaabaGaamOBaaaaaaa@3F79@

For face centered cubic the recommended smoothing length h in /PROP/SPH is the pitch entered when created the sph mesh in HyperMesh. The pitch h0 is the smallest distance between the particles. A smoothing length smaller than h0 can only be used when there is no tension physical problems material. If the material does include tensile behavior, then a smoothing length larger than h0 can be used to increase stability but there will be an increase in the computational cost.

Simple Cubic Net

Let, c the side length of each elementary cube into the net.


Figure 4.
The mass of the particles mp should relate to the density of the material ρ and to the size c of the net, with respect to the following equation:(5)
m p = c 3 ρ
Table 2. Number of Neighbors in a Cubic Net
Distance d Number of particles at distance d Number of particles within distance d
c 6 6
2 c MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aadaGcaaqaaiaaikdaaSqabaGccaWGJbaaaa@3AA1@ 12 18
3 c MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aadaGcaaqaaiaaikdaaSqabaGccaWGJbaaaa@3AA1@ 8 26
2c 6 32
5 c MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aadaGcaaqaaiaaikdaaSqabaGccaWGJbaaaa@3AA1@ 24 56
6 c MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aadaGcaaqaaiaaikdaaSqabaGccaWGJbaaaa@3AA1@ 24 80
2 c MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aadaGcaaqaaiaaikdaaSqabaGccaWGJbaaaa@3AA1@ 12 92
3c 6 98

From experience using cubic net, a higher smoothing length compared to face centered cubic or hexagonal compact net is needed to solve the tension instability. This higher smoothing length increases the computational cost since more neighbor particles have to be included in the calculation for each particle.

For cubic net, a smoothing length h between 1.25c and 1.5c is recommended in /PROP/SPH.