QuadBendingFEMForceField
Bending Quad finite elements
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 |
| quadInfo | Internal quad data | |
| vertexInfo | Internal node data | |
| edgeInfo | Internal edge data | |
| method | large: large displacements, small: small displacements | small |
| poissonRatio | Poisson ratio in Hooke's law (vector) | 0.45 |
| youngModulus | Young modulus in Hooke's law (vector) | 1000 |
| thickness | Thickness of the elements | 1 |
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/QuadBendingFEMForceField.scn
<!-- Mechanical QuadbendingFEMForceField Example -->
<Node name="root" gravity="-1 0 0" dt="0.02">
<RequiredPlugin name="Sofa.Component.Collision.Detection.Algorithm"/> <!-- Needed to use components [BVHNarrowPhase BruteForceBroadPhase CollisionPipeline] -->
<RequiredPlugin name="Sofa.Component.Collision.Detection.Intersection"/> <!-- Needed to use components [MinProximityIntersection] -->
<RequiredPlugin name="Sofa.Component.Collision.Response.Contact"/> <!-- Needed to use components [CollisionResponse] -->
<RequiredPlugin name="Sofa.Component.Constraint.Projective"/> <!-- Needed to use components [FixedProjectiveConstraint] -->
<RequiredPlugin name="Sofa.Component.IO.Mesh"/> <!-- Needed to use components [MeshVTKLoader] -->
<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 [QuadBendingFEMForceField] -->
<RequiredPlugin name="Sofa.Component.StateContainer"/> <!-- Needed to use components [MechanicalObject] -->
<RequiredPlugin name="Sofa.Component.Topology.Container.Dynamic"/> <!-- Needed to use components [QuadSetGeometryAlgorithms QuadSetTopologyContainer QuadSetTopologyModifier] -->
<RequiredPlugin name="Sofa.Component.Visual"/> <!-- Needed to use components [VisualStyle] -->
<RequiredPlugin name="Sofa.GL.Component.Rendering3D"/> <!-- Needed to use components [OglModel] -->
<DefaultAnimationLoop/>
<VisualStyle displayFlags="showVisual showForceFields showBehaviorModels" />
<CollisionPipeline name="CollisionPipeline" verbose="0" />
<BruteForceBroadPhase/>
<BVHNarrowPhase/>
<CollisionResponse name="collision response" response="PenalityContactForceField" />
<MinProximityIntersection name="proximity" alarmDistance="0.5" contactDistance="0.33"/>
<Node name="skin" gravity="0 0 -9.81">
<EulerImplicitSolver name="cg_odesolver" />
<CGLinearSolver iterations="25" name="linear solver" tolerance="1.0e-9" threshold="1.0e-9" />
<MeshVTKLoader filename="mesh/test_quad.vtk" name="loader" />
<MechanicalObject template="Vec3" src="@loader" name="DOFs" />
<QuadSetTopologyContainer src="@loader" name="topology" />
<QuadSetTopologyModifier name="Modifier" />
<QuadSetGeometryAlgorithms name="GeomAlgo" template="Vec3" />
<DiagonalMass massDensity="0.5" name="mass" />
<QuadBendingFEMForceField name="FEM" youngModulus="2000" poissonRatio="0.49" thickness="0.5" method="small" />
<FixedProjectiveConstraint indices="0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58" />
<Node name="Visual">
<OglModel name="Visual" color="yellow" />
<IdentityMapping template="Vec3,Vec3" name="visualMapping" input="@../DOFs" output="@Visual" />
</Node>
</Node>
</Node>
def createScene(rootNode):
root = rootNode.addChild('root', gravity="-1 0 0", dt="0.02")
root.addObject('RequiredPlugin', name="Sofa.Component.Collision.Detection.Algorithm")
root.addObject('RequiredPlugin', name="Sofa.Component.Collision.Detection.Intersection")
root.addObject('RequiredPlugin', name="Sofa.Component.Collision.Response.Contact")
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('DefaultAnimationLoop')
root.addObject('VisualStyle', displayFlags="showVisual showForceFields showBehaviorModels")
root.addObject('CollisionPipeline', name="CollisionPipeline", verbose="0")
root.addObject('BruteForceBroadPhase')
root.addObject('BVHNarrowPhase')
root.addObject('CollisionResponse', name="collision response", response="PenalityContactForceField")
root.addObject('MinProximityIntersection', name="proximity", alarmDistance="0.5", contactDistance="0.33")
skin = root.addChild('skin', gravity="0 0 -9.81")
skin.addObject('EulerImplicitSolver', name="cg_odesolver")
skin.addObject('CGLinearSolver', iterations="25", name="linear solver", tolerance="1.0e-9", threshold="1.0e-9")
skin.addObject('MeshVTKLoader', filename="mesh/test_quad.vtk", name="loader")
skin.addObject('MechanicalObject', template="Vec3", src="@loader", name="DOFs")
skin.addObject('QuadSetTopologyContainer', src="@loader", name="topology")
skin.addObject('QuadSetTopologyModifier', name="Modifier")
skin.addObject('QuadSetGeometryAlgorithms', name="GeomAlgo", template="Vec3")
skin.addObject('DiagonalMass', massDensity="0.5", name="mass")
skin.addObject('QuadBendingFEMForceField', name="FEM", youngModulus="2000", poissonRatio="0.49", thickness="0.5", method="small")
skin.addObject('FixedProjectiveConstraint', indices="0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58")
Visual = skin.addChild('Visual')
Visual.addObject('OglModel', name="Visual", color="yellow")
Visual.addObject('IdentityMapping', template="Vec3,Vec3", name="visualMapping", input="@../DOFs", output="@Visual")