Skip to content

TriangularFEMForceFieldOptim

Corotational Triangular finite elements.

Vec3d

Templates:

  • Vec3d

Target: Sofa.Component.SolidMechanics.FEM.Elastic

namespace: sofa::component::solidmechanics::fem::elastic

parents:

  • BaseLinearElasticityFEMForceField

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
rayleighStiffness Rayleigh damping - stiffness matrix coefficient 0
poissonRatio FEM Poisson Ratio in Hooke's law [0,0.5[ 0.45
youngModulus FEM Young's Modulus in Hooke's law 5000
triangleInfo Internal triangle data (persistent)
triangleState Internal triangle data (time-dependent)
damping Ratio damping/stiffness 0
restScale Scale factor applied to rest positions (to simulate pre-stretched materials) 1
computePrincipalStress Compute principal stress for each triangle 0
stressMaxValue Max stress value computed over the triangulation 0
Visualization
showStressVector Flag activating rendering of stress directions within each triangle 0
showStressThreshold Threshold value to render only stress vectors higher to this threshold 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
mechanicalStates List of mechanical states to which this component is associated BaseMechanicalState
mstate MechanicalState used by this component MechanicalState<Vec3d>
topology Link to a topology BaseMeshTopology

Examples

TriangularFEMForceFieldOptim.scn

<!-- Mechanical TriangularFEMForceFieldOptim Example -->
<Node name="root" dt="0.05" gravity="0 10 10" showBoundingTree="0">
    <RequiredPlugin pluginName="Sofa.Component.Constraint.Projective"/> <!-- Needed to use components [FixedProjectiveConstraint] -->
    <RequiredPlugin pluginName="Sofa.Component.IO.Mesh"/> <!-- Needed to use components [MeshGmshLoader] -->
    <RequiredPlugin pluginName="Sofa.Component.LinearSolver.Iterative"/> <!-- Needed to use components [CGLinearSolver] -->
    <RequiredPlugin pluginName="Sofa.Component.Mapping.Linear"/> <!-- Needed to use components [IdentityMapping] -->
    <RequiredPlugin pluginName="Sofa.Component.Mass"/> <!-- Needed to use components [DiagonalMass] -->
    <RequiredPlugin pluginName="Sofa.Component.IntegrationScheme.Backward"/> <!-- Needed to use components [EulerImplicitIntegrationScheme] -->
    <RequiredPlugin pluginName="Sofa.Component.SolidMechanics.FEM.Elastic"/> <!-- Needed to use components [TriangularFEMForceFieldOptim] -->
    <RequiredPlugin pluginName="Sofa.Component.StateContainer"/> <!-- Needed to use components [MechanicalObject] -->
    <RequiredPlugin pluginName="Sofa.Component.Topology.Container.Dynamic"/> <!-- Needed to use components [TriangleSetGeometryAlgorithms TriangleSetTopologyContainer TriangleSetTopologyModifier] -->
    <RequiredPlugin pluginName="Sofa.Component.Visual"/> <!-- Needed to use components [VisualStyle] -->
    <RequiredPlugin pluginName="Sofa.GL.Component.Rendering3D"/> <!-- Needed to use components [OglModel] -->

    <VisualStyle displayFlags="showVisual showBehaviorModels showForceFields showWireframe" />
    <DefaultAnimationLoop/>
    <!-- Activate this loader to use a square mesh with only two triangles (useful to debug base equations) -->
    <!--<MeshGmshLoader name="loaderSquare" triangles="0 1 3  1 2 3" position="0 0 0  1 0 0  1 1 0  0 1 0" />-->
    <!-- Activate this loader to load a square mesh with many triangles -->
    <MeshGmshLoader filename="mesh/square3.msh" name="loaderSquare" />
    <Node name="SquareGravity1">
        <EulerImplicitIntegrationScheme name="odesolver1" printLog="0"  rayleighStiffness="0.1" rayleighMass="0.1" />
        <CGLinearSolver printLog="0" iterations="25" name="linearsolver1" tolerance="1.0e-12" threshold="1.0e-12" />
        <TriangleSetTopologyContainer name="Container" src="@../loaderSquare" />
        <MechanicalObject name="DOFs" src="@../loaderSquare" scale="100"  />
        <TriangleSetTopologyModifier name="Modifier" />
        <TriangleSetGeometryAlgorithms name="GeomAlgo" template="Vec3" />
        <DiagonalMass massDensity="0.005" />
        <FixedProjectiveConstraint indices="0 1 2" />
        <TriangularFEMForceFieldOptim name="FEM" youngModulus="600" poissonRatio="0.3" method="large" printLog="1"/>
        <Node name="VisuA">
            <OglModel name="Visual" color="yellow" />
            <IdentityMapping name="visualMapping" input="@../DOFs" output="@Visual" />
        </Node>
    </Node>
    <!-- Activate this version to compare computed stiffness matrix (addKToMatrix) with addDForce -->
    <!--<Node name="SquareGravityTestMatrixConstruction">
        <EulerImplicitIntegrationScheme name="odesolver2" printLog="0" />
        <CGLinearSolver template="SparseMatrix" verbose="0" printLog="1" iterations="25" name="linearsolver2" tolerance="1.0e-12" threshold="1.0e-12" />
        <TriangleSetTopologyContainer name="Container" src="@../loaderSquare" />
        <MechanicalObject name="DOFs" src="@../loaderSquare" scale="100"  />
        <TriangleSetTopologyModifier name="Modifier" />
        <TriangleSetGeometryAlgorithms name="GeomAlgo" template="Vec3" />
        <DiagonalMass massDensity="0.005" />
        <FixedProjectiveConstraint indices="0 1 2" />
        <TriangularFEMForceFieldOptim name="FEM" youngModulus="600" poissonRatio="0.3" method="large" printLog="1"/>
        <Node name="VisuA">
            <OglModel name="Visual" color="yellow" />
            <IdentityMapping name="visualMapping" input="@../DOFs" output="@Visual" />
        </Node>
    </Node>-->

</Node>

def createScene(root_node):

   root = root_node.addChild('root', dt="0.05", gravity="0 10 10", showBoundingTree="0")

   root.addObject('RequiredPlugin', pluginName="Sofa.Component.Constraint.Projective")
   root.addObject('RequiredPlugin', pluginName="Sofa.Component.IO.Mesh")
   root.addObject('RequiredPlugin', pluginName="Sofa.Component.LinearSolver.Iterative")
   root.addObject('RequiredPlugin', pluginName="Sofa.Component.Mapping.Linear")
   root.addObject('RequiredPlugin', pluginName="Sofa.Component.Mass")
   root.addObject('RequiredPlugin', pluginName="Sofa.Component.IntegrationScheme.Backward")
   root.addObject('RequiredPlugin', pluginName="Sofa.Component.SolidMechanics.FEM.Elastic")
   root.addObject('RequiredPlugin', pluginName="Sofa.Component.StateContainer")
   root.addObject('RequiredPlugin', pluginName="Sofa.Component.Topology.Container.Dynamic")
   root.addObject('RequiredPlugin', pluginName="Sofa.Component.Visual")
   root.addObject('RequiredPlugin', pluginName="Sofa.GL.Component.Rendering3D")
   root.addObject('VisualStyle', displayFlags="showVisual showBehaviorModels showForceFields showWireframe")
   root.addObject('DefaultAnimationLoop', )
   root.addObject('MeshGmshLoader', filename="mesh/square3.msh", name="loaderSquare")

   square_gravity1 = root.addChild('SquareGravity1')

   square_gravity1.addObject('EulerImplicitIntegrationScheme', name="odesolver1", printLog="0", rayleighStiffness="0.1", rayleighMass="0.1")
   square_gravity1.addObject('CGLinearSolver', printLog="0", iterations="25", name="linearsolver1", tolerance="1.0e-12", threshold="1.0e-12")
   square_gravity1.addObject('TriangleSetTopologyContainer', name="Container", src="@../loaderSquare")
   square_gravity1.addObject('MechanicalObject', name="DOFs", src="@../loaderSquare", scale="100")
   square_gravity1.addObject('TriangleSetTopologyModifier', name="Modifier")
   square_gravity1.addObject('TriangleSetGeometryAlgorithms', name="GeomAlgo", template="Vec3")
   square_gravity1.addObject('DiagonalMass', massDensity="0.005")
   square_gravity1.addObject('FixedProjectiveConstraint', indices="0 1 2")
   square_gravity1.addObject('TriangularFEMForceFieldOptim', name="FEM", youngModulus="600", poissonRatio="0.3", method="large", printLog="1")

   visu_a = SquareGravity1.addChild('VisuA')

   visu_a.addObject('OglModel', name="Visual", color="yellow")
   visu_a.addObject('IdentityMapping', name="visualMapping", input="@../DOFs", output="@Visual")