Skip to content

Hexa2PrismTopologicalMapping

Topological mapping where HexahedronSetTopology is converted to PrismSetTopology

Target: Sofa.Component.Topology.Mapping

namespace: sofa::component::topology::mapping

parents:

  • TopologicalMapping

Data

Name Description Default value
name object name unnamed
printLog if true, emits extra messages at runtime. 0
tags list of the subsets the object belongs to
bbox this object bounding box
componentState The state of the component among (Dirty, Valid, Undefined, Loading, Invalid). Undefined
listening if true, handle the events, otherwise ignore the events 0
Name Description Destination type name
context Graph Node containing this object (or BaseContext::getDefault() if no graph is used) BaseContext
slaves Sub-objects used internally by this object BaseComponent
master nullptr for regular objects, or master object for which this object is one sub-objects BaseComponent
input Input topology to map BaseMeshTopology
output Output topology to map BaseMeshTopology

Examples

Hexa2PrismTopologicalMapping.scn

<?xml version="1.0"?>
<Node name="root" dt="0.01" gravity="0 -9.81 0">

    <Node name="plugins">
        <RequiredPlugin pluginName="Sofa.Component.Constraint.Projective"/> <!-- Needed to use components [FixedProjectiveConstraint] -->
        <RequiredPlugin pluginName="Sofa.Component.Engine.Select"/> <!-- Needed to use components [BoxROI] -->
        <RequiredPlugin pluginName="Sofa.Component.LinearSolver.Direct"/> <!-- Needed to use components [SparseLDLSolver] -->
        <RequiredPlugin pluginName="Sofa.Component.Mass"/> <!-- Needed to use components [MeshMatrixMass] -->
        <RequiredPlugin pluginName="Sofa.Component.SolidMechanics.FEM.Elastic"/>
        <RequiredPlugin pluginName="Sofa.Component.StateContainer"/> <!-- Needed to use components [MechanicalObject] -->
        <RequiredPlugin pluginName="Sofa.Component.Topology.Container.Constant"/> <!-- Needed to use components [MeshTopology] -->
        <RequiredPlugin pluginName="Sofa.Component.Topology.Container.Grid"/> <!-- Needed to use components [RegularGridTopology] -->
        <RequiredPlugin pluginName="Sofa.Component.Topology.Mapping"/> <!-- Needed to use components [Hexa2PrismTopologicalMapping] -->
        <RequiredPlugin pluginName="Sofa.Component.Visual"/> <!-- Needed to use components [VisualMesh] -->
    </Node>

    <DefaultAnimationLoop/>
    <VisualStyle displayFlags="showBehaviorModels showForceFields" />

    <VisualGrid size="0.1"/>
    <LineAxis size="0.1"/>
    <OglSceneFrame/>

    <EulerImplicitSolver name="backward Euler" rayleighStiffness="0.1" rayleighMass="0.1" />
    <SparseLDLSolver/>

    <RegularGridTopology name="grid" min="-0.01 -0.01 0" max="0.01 0.01 0.2" n="5 5 30"/>
    <MechanicalObject template="Vec3" name="state" showObject="true"/>

    <MeshMatrixMass massDensity="1100"/>

    <Node name="prisms">
        <MeshTopology name="prism_topology"/>
        <Hexa2PrismTopologicalMapping input="@grid" output="@prism_topology" />

        <PrismCorotationalFEMForceField name="FEM" youngModulus="2e6" poissonRatio="0.45" topology="@prism_topology"
                                        rotationMethod="polar" computeForceStrategy="sequenced" computeForceDerivStrategy="sequenced"/>

        <VisualMesh position="@../state.position" topology="@prism_topology" enable="true"/>
    </Node>

    <BoxROI template="Vec3" name="box_roi" box="-0.011 -0.011 -0.0001   0.011 0.011 0.0001" drawBoxes="1" />
    <FixedProjectiveConstraint template="Vec3" indices="@box_roi.indices" />
</Node>

def createScene(root_node):

   root = root_node.addChild('root', dt="0.01", gravity="0 -9.81 0")

   plugins = root.addChild('plugins')

   plugins.addObject('RequiredPlugin', pluginName="Sofa.Component.Constraint.Projective")
   plugins.addObject('RequiredPlugin', pluginName="Sofa.Component.Engine.Select")
   plugins.addObject('RequiredPlugin', pluginName="Sofa.Component.LinearSolver.Direct")
   plugins.addObject('RequiredPlugin', pluginName="Sofa.Component.Mass")
   plugins.addObject('RequiredPlugin', pluginName="Sofa.Component.SolidMechanics.FEM.Elastic")
   plugins.addObject('RequiredPlugin', pluginName="Sofa.Component.StateContainer")
   plugins.addObject('RequiredPlugin', pluginName="Sofa.Component.Topology.Container.Constant")
   plugins.addObject('RequiredPlugin', pluginName="Sofa.Component.Topology.Container.Grid")
   plugins.addObject('RequiredPlugin', pluginName="Sofa.Component.Topology.Mapping")
   plugins.addObject('RequiredPlugin', pluginName="Sofa.Component.Visual")

   root.addObject('DefaultAnimationLoop', )
   root.addObject('VisualStyle', displayFlags="showBehaviorModels showForceFields")
   root.addObject('VisualGrid', size="0.1")
   root.addObject('LineAxis', size="0.1")
   root.addObject('OglSceneFrame', )
   root.addObject('EulerImplicitSolver', name="backward Euler", rayleighStiffness="0.1", rayleighMass="0.1")
   root.addObject('SparseLDLSolver', )
   root.addObject('RegularGridTopology', name="grid", min="-0.01 -0.01 0", max="0.01 0.01 0.2", n="5 5 30")
   root.addObject('MechanicalObject', template="Vec3", name="state", showObject="true")
   root.addObject('MeshMatrixMass', massDensity="1100")

   prisms = root.addChild('prisms')

   prisms.addObject('MeshTopology', name="prism_topology")
   prisms.addObject('Hexa2PrismTopologicalMapping', input="@grid", output="@prism_topology")
   prisms.addObject('PrismCorotationalFEMForceField', name="FEM", youngModulus="2e6", poissonRatio="0.45", topology="@prism_topology", rotationMethod="polar", computeForceStrategy="sequenced", computeForceDerivStrategy="sequenced")
   prisms.addObject('VisualMesh', position="@../state.position", topology="@prism_topology", enable="true")

   root.addObject('BoxROI', template="Vec3", name="box_roi", box="-0.011 -0.011 -0.0001   0.011 0.011 0.0001", drawBoxes="1")
   root.addObject('FixedProjectiveConstraint', template="Vec3", indices="@box_roi.indices")