[Seria]

How can we use Carving on Grid Topology dear @Hugo ?

[Seria]

Dear @Hugo
Actually based on the below code my topology is triangles topology and i dont need mapping tetrahedron to triangle.

<?xml version="1.0"?>
<Node 	 name="1"  gravity="0 -10 0"  dt="0.01"  >
	<VisualStyle name="visualStyle1"  displayFlags="showVisual showBehavior showCollisionModels hideBoundingCollisionModels hideMapping hideOptions" />
	<Gravity name="gravity0" />
	<RequiredPlugin name="requiredPlugin3"  printLog="1"  pluginName="SPHFluidSurfaceMapping SPHFluidForceField" />
	<DefaultCollisionGroupManager name="defaultCollisionGroupManager4" />
	<BruteForceDetection name="bruteForceDetection3" />
	<MinProximityIntersection name="minProximityIntersection11"  alarmDistance="0.1" />
	<DefaultPipeline name="defaultPipeline2" />
	<DefaultContactManager name="Response"  response="default" />
	<EulerImplicitSolver name="eulerImplicitSolver5"  rayleighStiffness="0.1"  rayleighMass="0.1" />
	<CGLinearSolver template="GraphScattered" name="cGLinearSolver6" />
	<CarvingManager name="carvingManager5"  listening="1"  modelTool="Carv/tSphereModel6"  modelSurface="2/child/tTriangleModel5"  active="1" />
	<SimedixKeyboard name="sk"  listening="1"  Step="0.6"  tool="18 0 0 0 0 0 1" />
	<Node 	 name="center"  >
		<MechanicalObject template="Rigid3d" name="m"  position="0 0 0 0 0 0 1"  velocity="0 0 0 0 0 0"  force="0 0 0 0 0 0"  externalForce="0 0 0 0 0 0"  size="1" />
	</Node>
	<Node 	 name="ref_Point"  >
		<MechanicalObject template="Rigid3d" name="m"  position="@sk.tool"  velocity="0 0 0 0 0 0"  force="0 0 0 0 0 0"  externalForce="0 0 0 0 0 0"  translation="0 0 0"  size="1" />
		<PartialFixedConstraint template="Rigid3d" name="partialFixedConstraint23"  indices="0"  fixedDirections="1 1 1 1 1 1" />
	</Node>
	<Node 	 name="instrument"  >
		<MechanicalObject template="Rigid3d" name="m"  position="0 0 0 0 0 0 1"  velocity="0 0 0 0 0 0"  force="0 0 0 0 0 0"  externalForce="0 0 0 0 0 0"  translation="20 0 0"  size="1" />
		<Node 	 name="Vis"  gravity="0 0 0"  >
			<MeshObjLoader name="mesh"  filename="/home/masoud/Desktop/SiMedix_Scenes (copy)/data/mesh/phaco_col_3.obj"  translation="0 0 0" />
			<VisualStyle name="visualStyle8"  displayFlags="showVisual hideBehavior hideCollision hideMapping hideOptions" />
			<MeshTopology name="meshTopology57"  position="@mesh.position"  triangles="@mesh.triangles" />
			<MechanicalObject template="Vec3d" name="m"  position="0 0 0"  velocity="0 0 0"  force="0 0 0"  externalForce="0 0 0"  size="1" />
			<RigidMapping template="Rigid3d,Vec3d" name="mapp"  input="@../"  output="@./" />
			<Node 	 name="vis"  >
				<OglModel template="ExtVec3f" name="oglModel18"  primitiveType="DEFAULT"  blendEquation="GL_FUNC_ADD"  sfactor="GL_SRC_ALPHA"  dfactor="GL_ONE_MINUS_SRC_ALPHA" />
				<IdentityMapping template="Vec3d,ExtVec3f" name="identityMap1"  input="@../"  output="@./" />
			</Node>
		</Node>
		<JointSpringForceField template="Rigid3d" name="jointSpringFF16"  spring="BEGIN_SPRING  0 0  END_SPRING
"  object1="@/ref_Point/m"  object2="@./" />
	</Node>
	<Node 	 name="Orient"  >
		<MechanicalObject template="Rigid3d" name="m"  position="-4 0 0 0 0 0 1"  velocity="0 0 0 0 0 0"  force="0 0 0 0 0 0"  externalForce="0 0 0 0 0 0"  size="1" />
		<RigidRigidMapping template="Rigid3d,Rigid3d" name="rigidRigidMap1"  index="0"  input="@/instrument/m"  output="@./" />
	</Node>
	<Node 	 name="instrutip"  >
		<MechanicalObject template="Rigid3d" name="m"  position="-2 0 0 0 0 0 1"  velocity="0 0 0 0 0 0"  force="0 0 0 0 0 0"  externalForce="0 0 0 0 0 0"  size="1" />
		<RigidRigidMapping template="Rigid3d,Rigid3d" name="rigidRigidMap1"  index="0"  input="@/instrument/m"  output="@./" />
	</Node>
	<Node 	 name="Carv"  >
		<MechanicalObject template="Rigid3d" name="m"  position="-2 0 0 0 0 0 1"  velocity="0 0 0 0 0 0"  force="0 0 0 0 0 0"  externalForce="0 0 0 0 0 0"  size="1" />
		<TSphereModel template="Rigid3d" name="tSphereModel6"  contactStiffness="0"  color="1 0 0 1"  group="1"  radius="1.1" />
		<RigidRigidMapping template="Rigid3d,Rigid3d" name="rigidRigidMap1"  index="0"  input="@/instrument/m"  output="@./" />
	</Node>
	<Node 	 name="2"  >
		<LightManager name="lightManager2"  listening="1" />
		<DirectionalLight name="directionalLight1"  color="1 1 1 1"  modelViewMatrix="0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0"  projectionMatrix="0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0" />
		<RegularGridTopology name="g"  n="5 40 5"  min="-1.5 -3 -1.5"  max="0 12 0"  p0="-1.5 -3 -1.5" />
		<MechanicalObject template="Vec3d" name="MModel"  position="0 0 0"  velocity="0 0 0"  force="0 0 0"  externalForce="0 0 0"  size="1" />
		<SpatialGridContainer template="Vec3d" name="spatialGridContainer1"  listening="1"  cellWidth="0.75"  autoUpdate="1"  sortPoints="1" />
		<SPHFluidForceField template="Vec3d" name="sPHFluidFF1"  printLog="0"  radius="0.745"  pressure="100"  density="15"  viscosity="40"  surfaceTension="-1000"  kernelType="1"  viscosityType="2" />
		<UniformVelocityDampingForceField template="Vec3d" name="uniformVelocityDampingFF1"  dampingCoefficient="12" />
		<UniformMass template="Vec3d" name="uniformMass8"  filename="unused"  handleTopoChange="1" />
		<SimedixVaccumSphereForceField template="Vec3d" name="simedixVaccumSphereFF4"  center="@/instrutip/m.position"  radius="3"  apply_to_cone="1"  orientation="@/Orient/m.position"  angle="3"  stiffness="-20"  damping="0" />
		<Node 	 name="node7"  >
			<OglModel template="ExtVec3f" name="oglModel0"  primitiveType="DEFAULT"  blendEquation="GL_FUNC_ADD"  sfactor="GL_SRC_ALPHA"  dfactor="GL_ONE_MINUS_SRC_ALPHA" />
			<SPHFluidSurfaceMapping template="Vec3d,ExtVec3f" name="sPHFluidSurfaceMap2"  radius="0.745"  input="@../"  output="@./" />
		</Node>
		<Node 	 name="child"  >
			<TriangleSetTopologyAlgorithms template="Vec3d" name="triangleSetTopologyAlgorithms8" />
			<TriangleSetGeometryAlgorithms template="Vec3d" name="triangleSetGeometryAlgorithms7" />
			<TriangleSetTopologyModifier name="triangleSetTopologyModifier6" />
			<TriangleSetTopologyContainer name="triangleSetTopologyContainer5"  position="@MModel.position"  edges="@g.edges"  checkConnexity="1"  triangles="@g.triangles" />
			<TTriangleModel template="Vec3d" name="tTriangleModel5"  color="1 0 0 1"  group="1" />
		</Node>
		<InteractionEllipsoidForceField template="Vec3d,Rigid3d" name="ellips"  center="0 0 0"  vradius="10 10 10"  stiffness="-10000"  damping="1"  color="0.909804 0.427451 0.427451 1"  object2_forces="0"  object1="@./"  object2="@center/m" />
	</Node>
</Node>

So the above example that you mentioned cant work on this type.
Thank you in advance for your responses dear @Hugo .

[Seria]

Thank you for your response. i will check it and make you informed about the result on the project.

[Seria]

Thank you @Hugo for your patient.
This make me happy working with you.

Honestly, i wanna use Carving on a SPHFluidForceField. but carving work only with triangle collision. so i need to mapp grid to triangle.

here is my example of SPHFluid:

<?xml version="1.0"?>
<Node 	 name="1"  gravity="0 -10 0"  dt="0.01"  >
	<VisualStyle name="visualStyle1"  displayFlags="showBehaviorModels showForceFields showCollisionModels" />
	<RequiredPlugin name="requiredPlugin3"  printLog="1"  pluginName="SPHFluidSurfaceMapping SPHFluidForceField" />
	<Node 	 name="2"  >
		<LightManager name="lightManager2"  listening="1" />
		<DirectionalLight name="directionalLight1"  color="1 1 1 1"  modelViewMatrix="0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0"  projectionMatrix="0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0" />
		<CGLinearSolver template="GraphScattered" name="cGLinearSolver7"  iterations="12"  tolerance="1e-10"  threshold="1e-10" />
		<EulerImplicitSolver name="eulerImplicitSolver6"  rayleighStiffness="0.1"  rayleighMass="0.1" />
		<MechanicalObject template="Vec3d" name="MModel"  position="0 0 0"  velocity="0 0 0"  force="0 0 0"  externalForce="0 0 0"  size="1" />
		<RegularGridTopology name="regularGridTopology1"  n="5 40 5"  min="-1.5 -3 -1.5"  max="0 12 0"  p0="-1.5 -3 -1.5" />
		<SpatialGridContainer template="Vec3d" name="spatialGridContainer1"  listening="1"  cellWidth="0.75"  autoUpdate="1"  sortPoints="1" />
		<SPHFluidForceField template="Vec3d" name="sPHFluidFF1"  printLog="0"  radius="0.745"  pressure="100"  density="15"  viscosity="40"  surfaceTension="-1000"  kernelType="1"  viscosityType="2" />
		<UniformVelocityDampingForceField template="Vec3d" name="uniformVelocityDampingFF1"  dampingCoefficient="12" />
		<UniformMass template="Vec3d" name="uniformMass8"  filename="unused"  handleTopoChange="1" />
		<PlaneForceField template="Vec3d" name="planeFF1"  normal="1 0 0"  d="-4"  color="0 0.498039 0.2 1" />
		<PlaneForceField template="Vec3d" name="planeFF2"  normal="-1 0 0"  d="-4" />
		<PlaneForceField template="Vec3d" name="planeFF3"  normal="0.5 1 0.1"  d="-4"  draw="1"  color="0 0.498039 0.2 1" />
		<PlaneForceField template="Vec3d" name="planeFF4"  normal="0 0 1"  d="-4" />
		<PlaneForceField template="Vec3d" name="planeFF5"  normal="0 0 -1"  d="-4" />
		<Node 	 name="node7"  >
			<OglModel template="ExtVec3f" name="oglModel0"  primitiveType="DEFAULT"  blendEquation="GL_FUNC_ADD"  sfactor="GL_SRC_ALPHA"  dfactor="GL_ONE_MINUS_SRC_ALPHA" />
			<SPHFluidSurfaceMapping template="Vec3d,ExtVec3f" name="sPHFluidSurfaceMap2"  radius="0.745"  input="@../"  output="@./" />
		</Node>
	</Node>
</Node>

[Seria]

@Hugo
So good Dear Hugo. you really help us with SOFA. I really enjoy getting involved with problems. This plugin helped me learn working with plugins.
Thank you.

[Seria]

Hi dear @Hugo
Actually is solved my problem with different approach. However you suggestion was good. but the developer of this component should work better on his plugin because this approach cause unstable situation to the simulation,

Thank you.

[Author]

Posts