pfacet engines

Hi,

I am trying to simulate the triaxial test with a flexible membrane in a cylindrical shape.

Background information:
Spherical particle.
Top and bottom plates use facets.
Lateral membrane use pfacet.

################################################
Right now, I am confused about the engines.
The following is an example of the engines.
#################################################
 ###################################################################
###################################################################
O.engines=[
 ForceResetter(),
 InsertionSortCollider([
  ####### for sphere and facet wall
  Bo1_Sphere_Aabb(), Bo1_Facet_Aabb(),
  #### for pfacet membrane and gridconnection
  Bo1_PFacet_Aabb(), Bo1_GridConnection_Aabb()]),
 InteractionLoop(
  [
  ## for sphere contact with sphere and facet wall
  Ig2_Sphere_Sphere_ScGeom6D(),
  Ig2_Facet_Sphere_ScGeom(),
  ## for contacts between: sphere/pfacet, node/node, connection/connection, pfacet/pfacet and so on.
  Ig2_GridNode_GridNode_GridNodeGeom6D(),
  Ig2_GridConnection_GridConnection_GridCoGridCoGeom(),
  Ig2_Sphere_PFacet_ScGridCoGeom(),
  Ig2_Facet_Sphere_ScGeom6D(),
  Ig2_Sphere_GridConnection_ScGridCoGeom(),
  Ig2_PFacet_PFacet_ScGeom(),
  Ig2_GridConnection_PFacet_ScGeom(),
  ],
  [
  Ip2_FrictMat_FrictMat_FrictPhys(),
  Ip2_CohFrictMat_CohFrictMat_CohFrictPhys(label="cohesiveIp")
  ],
  [
  Law2_ScGeom_FrictPhys_CundallStrack(),
  Law2_ScGeom6D_CohFrictPhys_CohesionMoment(useIncrementalForm=True,always_use_moment_law=False,label='cohesiveLaw'),
  Law2_GridCoGridCoGeom_FrictPhys_CundallStrack(),
  Law2_ScGridCoGeom_FrictPhys_CundallStrack(),
  ],label="iloop"),
 GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=25,timestepSafetyCoefficient=0.8),
 NewtonIntegrator(gravity=(0,0,0),damping=0.1,label='newton')]
###########################################################
O.engines[2].lawDispatcher.functors[1].always_use_moment_law=False
O.engines[2].physDispatcher.functors[1].setCohesionOnNewContacts=False
#######################################################################
#######################################################################

Questions:

1: Ig2_Facet_Sphere_ScGeom() and Ig2_Facet_Sphere_ScGeom6D().

Do I need to include both of them?
one is for sphere-facet, another one is for sphere-facet (facet from pfacet, correct me if I am wrong).
because as mentioned in [1], this is for the virtual sphere, so I need to add Ig2_Facet_Sphere_ScGeom() for the sphere which contacts the top and bottom facets.

2: For the sphere and PFacet contact, Do I need to add all the 7 Ig2 functors?
########################################################
  Ig2_GridNode_GridNode_GridNodeGeom6D(),
  Ig2_GridConnection_GridConnection_GridCoGridCoGeom(),
  Ig2_Sphere_PFacet_ScGridCoGeom(),
  Ig2_Facet_Sphere_ScGeom6D(),
  Ig2_Sphere_GridConnection_ScGridCoGeom(),
  Ig2_PFacet_PFacet_ScGeom(),
  Ig2_GridConnection_PFacet_ScGeom(),
#######################################################
Because from [2] (the author also simulates the triaxial test with pfacet), the code does not include the following Ig2 factors.
  Ig2_Sphere_GridConnection_ScGridCoGeom(),
  Ig2_PFacet_PFacet_ScGeom(),
  Ig2_GridConnection_PFacet_ScGeom(),
#######################################################

3: what is the difference between Ig2_Facet_Sphere_ScGeom6D() and Ig2_Sphere_PFacet_ScGridCoGeom()?
Because the pfacet-sphere contact also uses the virtual sphere to detect the contact with real spheres.

thank you!

references:
[1]https://yade-dev.gitlab.io/trunk/yade.wrapper.html#yade.wrapper.Ig2_Facet_Sphere_ScGeom6D
[2]https://www.proquest.com/openview/057ef8332b9c247ae1b67a35db05587e/1?pq-origsite=gscholar&cbl=18750&diss=y