TriangularFEMForceFieldOptim
Corotational Triangular finite elements Supports GPU-side computations using CUDA
Templates:
Vec3d
Target: Sofa.Component.SolidMechanics.FEM.Elastic
namespace: sofa::component::solidmechanics::fem::elastic
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 |
| triangleInfo | Internal triangle data (persistent) | |
| triangleState | Internal triangle data (time-dependent) | |
| poissonRatio | Poisson ratio in Hooke's law | 0.3 |
| youngModulus | Young modulus in Hooke's law | 1000 |
| 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 |
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/FEM/TriangularFEMForceFieldOptim.scn
<!-- Mechanical TriangularFEMForceFieldOptim Example -->
<Node name="root" dt="0.05" gravity="0 10 10" showBoundingTree="0">
<RequiredPlugin name="Sofa.Component.Constraint.Projective"/> <!-- Needed to use components [FixedProjectiveConstraint] -->
<RequiredPlugin name="Sofa.Component.IO.Mesh"/> <!-- Needed to use components [MeshGmshLoader] -->
<RequiredPlugin name="Sofa.Component.LinearSolver.Iterative"/> <!-- Needed to use components [CGLinearSolver] -->
<RequiredPlugin name="Sofa.Component.Mapping.Linear"/> <!-- Needed to use components [IdentityMapping] -->
<RequiredPlugin name="Sofa.Component.Mass"/> <!-- Needed to use components [DiagonalMass] -->
<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.StateContainer"/> <!-- Needed to use components [MechanicalObject] -->
<RequiredPlugin name="Sofa.Component.Topology.Container.Dynamic"/> <!-- Needed to use components [TriangleSetGeometryAlgorithms TriangleSetTopologyContainer TriangleSetTopologyModifier] -->
<RequiredPlugin name="Sofa.Component.Visual"/> <!-- Needed to use components [VisualStyle] -->
<RequiredPlugin name="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">
<EulerImplicitSolver name="odesolver1" printLog="0" rayleighStiffness="0.1" rayleighMass="0.1" />
<CGLinearSolver verbose="0" printLog="0" iterations="25" name="linearsolver1" tolerance="1.0e-9" threshold="1.0e-9" />
<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">
<EulerImplicitSolver name="odesolver2" printLog="0" />
<CGLinearSolver template="SparseMatrix" verbose="0" printLog="1" iterations="25" name="linearsolver2" tolerance="1.0e-9" threshold="1.0e-9" />
<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(rootNode):
root = rootNode.addChild('root', dt="0.05", gravity="0 10 10", showBoundingTree="0")
root.addObject('RequiredPlugin', name="Sofa.Component.Constraint.Projective")
root.addObject('RequiredPlugin', name="Sofa.Component.IO.Mesh")
root.addObject('RequiredPlugin', name="Sofa.Component.LinearSolver.Iterative")
root.addObject('RequiredPlugin', name="Sofa.Component.Mapping.Linear")
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.StateContainer")
root.addObject('RequiredPlugin', name="Sofa.Component.Topology.Container.Dynamic")
root.addObject('RequiredPlugin', name="Sofa.Component.Visual")
root.addObject('RequiredPlugin', name="Sofa.GL.Component.Rendering3D")
root.addObject('VisualStyle', displayFlags="showVisual showBehaviorModels showForceFields showWireframe")
root.addObject('DefaultAnimationLoop')
root.addObject('MeshGmshLoader', filename="mesh/square3.msh", name="loaderSquare")
SquareGravity1 = root.addChild('SquareGravity1')
SquareGravity1.addObject('EulerImplicitSolver', name="odesolver1", printLog="0", rayleighStiffness="0.1", rayleighMass="0.1")
SquareGravity1.addObject('CGLinearSolver', verbose="0", printLog="0", iterations="25", name="linearsolver1", tolerance="1.0e-9", threshold="1.0e-9")
SquareGravity1.addObject('TriangleSetTopologyContainer', name="Container", src="@../loaderSquare")
SquareGravity1.addObject('MechanicalObject', name="DOFs", src="@../loaderSquare", scale="100")
SquareGravity1.addObject('TriangleSetTopologyModifier', name="Modifier")
SquareGravity1.addObject('TriangleSetGeometryAlgorithms', name="GeomAlgo", template="Vec3")
SquareGravity1.addObject('DiagonalMass', massDensity="0.005")
SquareGravity1.addObject('FixedProjectiveConstraint', indices="0 1 2")
SquareGravity1.addObject('TriangularFEMForceFieldOptim', name="FEM", youngModulus="600", poissonRatio="0.3", method="large", printLog="1")
VisuA = SquareGravity1.addChild('VisuA')
VisuA.addObject('OglModel', name="Visual", color="yellow")
VisuA.addObject('IdentityMapping', name="visualMapping", input="@../DOFs", output="@Visual")