Skip to content

TetrahedralTensorMassForceField

Linear Elastic Tetrahedral Mesh Supports GPU-side computations using CUDA

Templates:

  • Vec3d

Target: Sofa.Component.SolidMechanics.TensorMass

namespace: sofa::component::solidmechanics::tensormass

parents:

  • ForceField

Data:

Name Description Default value
name object name unnamed
printLog if true, emits extra messages at runtime. 0
tags list of the subsets the objet 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
isCompliance Consider the component as a compliance, else as a stiffness 0
rayleighStiffness Rayleigh damping - stiffness matrix coefficient 0
poissonRatio Poisson ratio in Hooke's law 0.3
youngModulus Young's modulus in Hooke's law 1000
edgeInfo Internal edge data

Links:

Name Description
context Graph Node containing this object (or BaseContext::getDefault() if no graph is used)
slaves Sub-objects used internally by this object
master nullptr for regular objects, or master object for which this object is one sub-objects
mechanicalStates List of mechanical states to which this component is associated
mstate MechanicalState used by this component
topology link to the topology container

Examples

Component/SolidMechanics/TensorMass/TetrahedralTensorMassForceField.scn

<Node name="root">
    <Node name="plugins">
        <RequiredPlugin name="Sofa.Component.Constraint.Projective"/> <!-- Needed to use components [FixedProjectiveConstraint] -->
        <RequiredPlugin name="Sofa.Component.Engine.Select"/> <!-- Needed to use components [BoxROI] -->
        <RequiredPlugin name="Sofa.Component.LinearSolver.Direct"/> <!-- Needed to use components [EigenSimplicialLDLT] -->
        <RequiredPlugin name="Sofa.Component.Mapping.Linear"/> <!-- Needed to use components [IdentityMapping] -->
        <RequiredPlugin name="Sofa.Component.Mass"/> <!-- Needed to use components [UniformMass] -->
        <RequiredPlugin name="Sofa.Component.ODESolver.Backward"/> <!-- Needed to use components [EulerImplicitSolver] -->
        <RequiredPlugin name="Sofa.Component.SolidMechanics.TensorMass"/> <!-- Needed to use components [TetrahedralTensorMassForceField] -->
        <RequiredPlugin name="Sofa.Component.StateContainer"/> <!-- Needed to use components [MechanicalObject] -->
        <RequiredPlugin name="Sofa.Component.Topology.Container.Dynamic"/> <!-- Needed to use components [QuadSetGeometryAlgorithms QuadSetTopologyContainer QuadSetTopologyModifier TetrahedronSetGeometryAlgorithms TetrahedronSetTopologyContainer TetrahedronSetTopologyModifier] -->
        <RequiredPlugin name="Sofa.Component.Topology.Container.Grid"/> <!-- Needed to use components [RegularGridTopology] -->
        <RequiredPlugin name="Sofa.Component.Topology.Mapping"/> <!-- Needed to use components [Hexa2QuadTopologicalMapping Hexa2TetraTopologicalMapping] -->
        <RequiredPlugin name="Sofa.Component.Visual"/> <!-- Needed to use components [VisualStyle] -->
        <RequiredPlugin name="Sofa.GL.Component.Rendering3D"/> <!-- Needed to use components [OglModel] -->
    </Node>

    <VisualStyle displayFlags="showBehaviorModels showForceFields" />

    <DefaultAnimationLoop name="animationLoop"/>
    <DefaultVisualManagerLoop name="visualLoop"/>

    <EulerImplicitSolver name="odesolver" rayleighStiffness="0.1" rayleighMass="0.1" />
    <EigenSimplicialLDLT template="CompressedRowSparseMatrixMat3x3"/>
    <MechanicalObject name="DoFs" />
    <UniformMass name="mass" totalMass="320" />
    <RegularGridTopology name="grid" nx="4" ny="4" nz="20" xmin="-9" xmax="-6" ymin="0" ymax="3" zmin="0" zmax="19" />
    <BoxROI name="box" box="-10 -1 -0.0001  -5 4 0.0001"/>
    <FixedProjectiveConstraint indices="@box.indices" />

    <TetrahedronSetTopologyContainer name="Tetra_topo"/>
    <TetrahedronSetTopologyModifier name="Modifier" />
    <TetrahedronSetGeometryAlgorithms template="Vec3" name="GeomAlgo" />
    <Hexa2TetraTopologicalMapping input="@grid" output="@Tetra_topo" />
    <TetrahedralTensorMassForceField name="deformable" youngModulus="100000" poissonRatio="0.4" />

    <Node name="quads">
        <QuadSetTopologyContainer  name="Container" />
        <QuadSetTopologyModifier   name="Modifier" />
        <QuadSetGeometryAlgorithms name="GeomAlgo" template="Vec3" />
        <Hexa2QuadTopologicalMapping input="@../grid" output="@Container" />
        <Node name="Visu">
            <OglModel name="Visual" color="yellow" />
            <IdentityMapping input="@../../DoFs" output="@Visual" />
        </Node>
    </Node>
</Node>

def createScene(rootNode):

    root = rootNode.addChild('root')

    plugins = root.addChild('plugins')
    plugins.addObject('RequiredPlugin', name="Sofa.Component.Constraint.Projective")
    plugins.addObject('RequiredPlugin', name="Sofa.Component.Engine.Select")
    plugins.addObject('RequiredPlugin', name="Sofa.Component.LinearSolver.Direct")
    plugins.addObject('RequiredPlugin', name="Sofa.Component.Mapping.Linear")
    plugins.addObject('RequiredPlugin', name="Sofa.Component.Mass")
    plugins.addObject('RequiredPlugin', name="Sofa.Component.ODESolver.Backward")
    plugins.addObject('RequiredPlugin', name="Sofa.Component.SolidMechanics.TensorMass")
    plugins.addObject('RequiredPlugin', name="Sofa.Component.StateContainer")
    plugins.addObject('RequiredPlugin', name="Sofa.Component.Topology.Container.Dynamic")
    plugins.addObject('RequiredPlugin', name="Sofa.Component.Topology.Container.Grid")
    plugins.addObject('RequiredPlugin', name="Sofa.Component.Topology.Mapping")
    plugins.addObject('RequiredPlugin', name="Sofa.Component.Visual")
    plugins.addObject('RequiredPlugin', name="Sofa.GL.Component.Rendering3D")
    root.addObject('VisualStyle', displayFlags="showBehaviorModels showForceFields")
    root.addObject('DefaultAnimationLoop', name="animationLoop")
    root.addObject('DefaultVisualManagerLoop', name="visualLoop")
    root.addObject('EulerImplicitSolver', name="odesolver", rayleighStiffness="0.1", rayleighMass="0.1")
    root.addObject('EigenSimplicialLDLT', template="CompressedRowSparseMatrixMat3x3")
    root.addObject('MechanicalObject', name="DoFs")
    root.addObject('UniformMass', name="mass", totalMass="320")
    root.addObject('RegularGridTopology', name="grid", nx="4", ny="4", nz="20", xmin="-9", xmax="-6", ymin="0", ymax="3", zmin="0", zmax="19")
    root.addObject('BoxROI', name="box", box="-10 -1 -0.0001  -5 4 0.0001")
    root.addObject('FixedProjectiveConstraint', indices="@box.indices")
    root.addObject('TetrahedronSetTopologyContainer', name="Tetra_topo")
    root.addObject('TetrahedronSetTopologyModifier', name="Modifier")
    root.addObject('TetrahedronSetGeometryAlgorithms', template="Vec3", name="GeomAlgo")
    root.addObject('Hexa2TetraTopologicalMapping', input="@grid", output="@Tetra_topo")
    root.addObject('TetrahedralTensorMassForceField', name="deformable", youngModulus="100000", poissonRatio="0.4")

    quads = root.addChild('quads')
    quads.addObject('QuadSetTopologyContainer', name="Container")
    quads.addObject('QuadSetTopologyModifier', name="Modifier")
    quads.addObject('QuadSetGeometryAlgorithms', name="GeomAlgo", template="Vec3")
    quads.addObject('Hexa2QuadTopologicalMapping', input="@../grid", output="@Container")

    Visu = quads.addChild('Visu')
    Visu.addObject('OglModel', name="Visual", color="yellow")
    Visu.addObject('IdentityMapping', input="@../../DoFs", output="@Visual")