Numerical Simulation Study on Flow and Fracture of Granular Materials
Algorithm 4
Beam element generation algorithm.
(1)
Obtain the detailed contact information of tetrahedron according to the algorithm in Section 2.2.
(2)
Cycle for all tetrahedrons (for (i = 0; i < num; i++) (num is the number of tetrahedrons)).
(a) Cycle all tetrahedrons in contact with tetrahedron j (for (j = 0; j < neighborNum; j++) (neighborNum is the number of tetrahedrons in contact with tetrahedron i)).
If there is no beam between tetrahedron i and tetrahedron j.
① Get the centroid of tetrahedron i and neighbor unit j.
② The starting points and ending points of the beam are set to be the mass centers of two tetrahedrons, respectively, and the length of the beam is the distance between the centers of mass.
③ Set the cross-sectional area of the beam as the contact area.
④ Girder is endowed with mechanical properties parameters.
