Applying mechanical states of IdentityMultiMapping to a model that is encompassed by these mechanical states

[zbounik]

Hi,
There is an object (sphere) that is partly surrounded by deformable cube A. The rest of this sphere is surrounded by deformable cube B. The part of the sphere that is inside cube A should move according to it, and the other part should move according to cube B. I used IdentityMultiMapping to concatanate mechanical states of A and B, and then barycentricMapping to map these DOFs to the sphere, but this is not giving the desirable result. I am wondering if using IdentityMultiMapping is not a good solution for mapping the cube deformation to a part of the sphere, or I am doing something wrong here.

I appreciate any hint in this regard.
Below you can find the scene:

<VisualStyle name="visualStyle1"  displayFlags="showVisual showBehaviorModels" />

<!--my change in visualStyle2-->
<VisualStyle name="visualStyle2"  displayFlags="showBehaviorModels" />
<DefaultVisualManagerLoop name="defaultVisualManagerLoop1" />
<FreeMotionAnimationLoop name="FreeMotionAnimationLoop" />
<DefaultPipeline name="CollisionPipeline" />
<MinProximityIntersection name="Proximity"  alarmDistance="0.5"  contactDistance="1" />
<DefaultContactManager name="Response"  response="default" />
<LCPConstraintSolver name="lCPConstraintSolver1"  build_lcp="0"  tolerance="0.0001"  maxIt="1000"  group="0" />
<BruteForceBroadPhase/>
<BVHNarrowPhase/>
<EulerImplicitSolver name="eulerImplicitSolver3" />
<CGLinearSolver iterations="25" name="linear solver" tolerance="1.0e-9" threshold="1.0e-9" />


<Node 	 name="A"  >
	<MeshObjLoader name="meshLoader_0"  filename="mesh/cube.obj" translation="-1 0 0"/>
	<OglModel template="Vec3d" name="Visual" src="@meshLoader_0" position="@meshLoader_0.position"  edges="@meshLoader_0.edges"  triangles="@meshLoader_0.triangles"  quads="@meshLoader_0.quads"  material="Default Diffuse 1 0.2 1 1 0.5 Ambient 1 0.04 0.2 0.2 0.5 Specular 0 0.2 1 1 0.5 Emissive 0 0.2 1 1 0.5 Shininess 0 45 "  primitiveType="DEFAULT"  blendEquation="GL_FUNC_ADD"  sfactor="GL_SRC_ALPHA"  dfactor="GL_ONE_MINUS_SRC_ALPHA" />
	<MechanicalObject name="MechA" src="@meshLoader_0"/>
	<UniformMass totalMass="60"  name="uniformMass" />
	<TriangleSetTopologyContainer name="triangleSetTopologyContainer10" src="@meshLoader_0"/>
	<TriangleSetTopologyModifier name="triangleSetTopologyModifier13" />
	<TriangleSetGeometryAlgorithms template="Vec3d" name="triangleSetGeometryAlgorithms19" />
	<TriangularFEMForceField template="Vec3d" name="triangularFEMFF21" poissonRatio="0.49" youngModulus="10000" method="large"/>
	<BarycentricMapping input="@MechA" output="@Visual" name="visual mapping" />
	<UncoupledConstraintCorrection />
    <TriangleCollisionModel />
    <FixedConstraint indices="3 6 1" />

</Node>

<Node 	 name="B"  >
	<MeshObjLoader name="meshLoader_0"  filename="mesh/cube.obj" translation="1 0 0"/>
	<OglModel template="Vec3d" name="Visual" src="@meshLoader_0" position="@meshLoader_0.position"  edges="@meshLoader_0.edges"  triangles="@meshLoader_0.triangles"  quads="@meshLoader_0.quads"  material="Default Diffuse 1 1 1 0.2 0.5 Ambient 1 0.2 0.2 0.04 0.5 Specular 0 1 1 0.2 0.5 Emissive 0 1 1 0.2 0.5 Shininess 0 45 "  primitiveType="DEFAULT"  blendEquation="GL_FUNC_ADD"  sfactor="GL_SRC_ALPHA"  dfactor="GL_ONE_MINUS_SRC_ALPHA" />
	<MechanicalObject name="MechB" src="@meshLoader_0" />
	<UniformMass totalMass="60"  name="uniformMass" />
	<TriangleSetTopologyContainer name="triangleSetTopologyContainer10" src="@meshLoader_0"/>
	<TriangleSetTopologyModifier name="triangleSetTopologyModifier13" />
	<TriangleSetGeometryAlgorithms template="Vec3d" name="triangleSetGeometryAlgorithms19" />
	<TriangularFEMForceField template="Vec3d" name="triangularFEMFF21" poissonRatio="0.49" youngModulus="10000" method="large"/>
	<BarycentricMapping input="@MechB" output="@Visual" name="visual mapping" />
	<UncoupledConstraintCorrection />
    <TriangleCollisionModel />
    <FixedConstraint indices="0 1 4" />


</Node>


<Node name="concatenation">
    <MechanicalObject template="Vec3d" name="dofall" showObject="1"/>
    <IdentityMultiMapping template="Vec3d,Vec3d" input="@../A/MechA @../B/MechB" output="@dofall" />
	<TriangleSetTopologyContainer name="triangleSetTopologyContainer10"/>
	<TriangleSetTopologyModifier name="triangleSetTopologyModifier13" />
	<TriangleSetGeometryAlgorithms template="Vec3d" name="triangleSetGeometryAlgorithms19" />

</Node>

<Node name="sphere"  >
	<MeshObjLoader name="meshLoader_0"  filename="mesh/sphere.obj"/>
	<OglModel template="Vec3d" name="Visual"  position="@meshLoader_0.position"  edges="@meshLoader_0.edges"  triangles="@meshLoader_0.triangles"  quads="@meshLoader_0.quads"  material="Default Diffuse 1 1 1 1 1 Ambient 1 0.2 0.2 0.2 1 Specular 0 1 1 1 1 Emissive 0 1 1 1 1 Shininess 0 45 "  primitiveType="DEFAULT"  blendEquation="GL_FUNC_ADD"  sfactor="GL_SRC_ALPHA"  dfactor="GL_ONE_MINUS_SRC_ALPHA" />
	<TriangleSetTopologyContainer name="triangleSetTopologyContainer10" src="@meshLoader_0"/>
	<TriangleSetGeometryAlgorithms template="Vec3d" name="triangleSetGeometryAlgorithms19" />
	<MechanicalObject name="MechT" src="@meshLoader_0" scale="0.5"/>
	<BarycentricMapping input="@../concatenation/dofall" output="@Visual" />

</Node>

[zbounik]

I used mergeMesh along with IdentityMultiMapping to also concatenate the meshes of the two cubes, and then used BarycentricMapping for mapping to the sphere. This solved the problem.