Skip to content

FastTriangularBendingSprings

Springs added to a triangular mesh to prevent bending.

Vec3d

Templates:

  • Vec3d

Target: Sofa.Component.SolidMechanics.Spring

namespace: sofa::component::solidmechanics::spring

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
rayleighStiffness Rayleigh damping - stiffness matrix coefficient 0
bendingStiffness Bending stiffness of the material 1
minDistValidity Distance under which a spring is not valid 1e-06
edgeInfo Internal edge data
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 BaseObject
master nullptr for regular objects, or master object for which this object is one sub-objects BaseObject
mechanicalStates List of mechanical states to which this component is associated BaseMechanicalState
mstate MechanicalState used by this component MechanicalState<Vec3d>
topology link to the topology container BaseMeshTopology

Examples

FastTriangularBendingSprings.scn

<?xml version="1.0" ?>
<Node name="root" dt="0.01" gravity="0 0 -1">
    <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.IO.Mesh"/> <!-- Needed to use components [MeshOBJLoader] -->
    <RequiredPlugin name="Sofa.Component.LinearSolver.Iterative"/> <!-- Needed to use components [CGLinearSolver] -->
    <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.FEM.Elastic"/> <!-- Needed to use components [TriangularFEMForceFieldOptim] -->
    <RequiredPlugin name="Sofa.Component.SolidMechanics.Spring"/> <!-- Needed to use components [FastTriangularBendingSprings] -->
    <RequiredPlugin name="Sofa.Component.StateContainer"/> <!-- Needed to use components [MechanicalObject] -->
    <RequiredPlugin name="Sofa.Component.Topology.Container.Constant"/> <!-- Needed to use components [MeshTopology] -->
    <RequiredPlugin name="Sofa.Component.Topology.Container.Dynamic"/> <!-- Needed to use components [EdgeSetGeometryAlgorithms] -->
    <RequiredPlugin name="Sofa.Component.Visual"/> <!-- Needed to use components [VisualStyle] -->

    <DefaultAnimationLoop/>
    <VisualStyle displayFlags="showBehavior hideCollision hideVisual " />
    <EulerImplicitSolver  rayleighStiffness="0.1" rayleighMass="0.1" />
    <CGLinearSolver iterations="25" tolerance="1e-5" threshold="1e-5"/>
    <Node name="Thin shell">
        <MeshOBJLoader name="loader" filename="mesh/triangleGrid_10_10.obj" />
        <MeshTopology src="@loader" />
        <EdgeSetGeometryAlgorithms />
        <MechanicalObject name="defoDOF" template="Vec3"  src="@loader" showObject="1"/>
        <BoxROI name="box1" box="-0.5 -0.5 -0.5  100.5 0.005 0.005  " />
        <FixedProjectiveConstraint indices="@box1.indices"/>
        <TriangularFEMForceFieldOptim name="FEM1" youngModulus="20000" poissonRatio="0.3" method="large" />
        <FastTriangularBendingSprings bendingStiffness="10000" />
        <UniformMass totalMass="2500" printLog="0" />
    </Node>
</Node>

def createScene(root_node):

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

   root.addObject('RequiredPlugin', name="Sofa.Component.Constraint.Projective")
   root.addObject('RequiredPlugin', name="Sofa.Component.Engine.Select")
   root.addObject('RequiredPlugin', name="Sofa.Component.IO.Mesh")
   root.addObject('RequiredPlugin', name="Sofa.Component.LinearSolver.Iterative")
   root.addObject('RequiredPlugin', name="Sofa.Component.Mass")
   root.addObject('RequiredPlugin', name="Sofa.Component.ODESolver.Backward")
   root.addObject('RequiredPlugin', name="Sofa.Component.SolidMechanics.FEM.Elastic")
   root.addObject('RequiredPlugin', name="Sofa.Component.SolidMechanics.Spring")
   root.addObject('RequiredPlugin', name="Sofa.Component.StateContainer")
   root.addObject('RequiredPlugin', name="Sofa.Component.Topology.Container.Constant")
   root.addObject('RequiredPlugin', name="Sofa.Component.Topology.Container.Dynamic")
   root.addObject('RequiredPlugin', name="Sofa.Component.Visual")
   root.addObject('DefaultAnimationLoop', )
   root.addObject('VisualStyle', displayFlags="showBehavior hideCollision hideVisual ")
   root.addObject('EulerImplicitSolver', rayleighStiffness="0.1", rayleighMass="0.1")
   root.addObject('CGLinearSolver', iterations="25", tolerance="1e-5", threshold="1e-5")

   thin_shell = root.addChild('Thin shell')

   thin_shell.addObject('MeshOBJLoader', name="loader", filename="mesh/triangleGrid_10_10.obj")
   thin_shell.addObject('MeshTopology', src="@loader")
   thin_shell.addObject('EdgeSetGeometryAlgorithms', )
   thin_shell.addObject('MechanicalObject', name="defoDOF", template="Vec3", src="@loader", showObject="1")
   thin_shell.addObject('BoxROI', name="box1", box="-0.5 -0.5 -0.5  100.5 0.005 0.005  ")
   thin_shell.addObject('FixedProjectiveConstraint', indices="@box1.indices")
   thin_shell.addObject('TriangularFEMForceFieldOptim', name="FEM1", youngModulus="20000", poissonRatio="0.3", method="large")
   thin_shell.addObject('FastTriangularBendingSprings', bendingStiffness="10000")
   thin_shell.addObject('UniformMass', totalMass="2500", printLog="0")