Skip to content

[examples] Validation: beam extension/compression #6294

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
alxbilger wants to merge 1 commit into
base: master
Choose a base branch
from
Open

[examples] Validation: beam extension/compression #6294

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
41 changes: 41 additions & 0 deletions examples/Validation/beam_extension_compression/beam.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,41 @@
<Node name="Group">

<Node name="plugins">
<RequiredPlugin pluginName="Sofa.Component.Constraint.Projective"/> <!-- Needed to use components [FixedProjectiveConstraint] -->
<RequiredPlugin pluginName="Sofa.Component.Engine.Select"/> <!-- Needed to use components [BoxROI] -->
<RequiredPlugin pluginName="Sofa.Component.LinearSolver.Direct"/> <!-- Needed to use components [SparseLDLSolver] -->
<RequiredPlugin pluginName="Sofa.Component.LinearSolver.Iterative"/> <!-- Needed to use components [CGLinearSolver] -->
<RequiredPlugin pluginName="Sofa.Component.LinearSolver.Ordering"/> <!-- Needed to use components [NaturalOrderingMethod] -->
<RequiredPlugin pluginName="Sofa.Component.LinearSystem"/> <!-- Needed to use components [ConstantSparsityPatternSystem] -->
<RequiredPlugin pluginName="Sofa.Component.Mass"/> <!-- Needed to use components [MeshMatrixMass] -->
<RequiredPlugin pluginName="Sofa.Component.ODESolver.Backward"/> <!-- Needed to use components [EulerImplicitSolver] -->
<RequiredPlugin pluginName="Sofa.Component.SolidMechanics.FEM.Elastic"/>
<RequiredPlugin pluginName="Sofa.Component.StateContainer"/> <!-- Needed to use components [MechanicalObject] -->
<RequiredPlugin pluginName="Sofa.Component.Topology.Container.Dynamic"/> <!-- Needed to use components [TetrahedronSetGeometryAlgorithms,TetrahedronSetTopologyContainer,TetrahedronSetTopologyModifier] -->
<RequiredPlugin pluginName="Sofa.Component.Topology.Container.Grid"/> <!-- Needed to use components [RegularGridTopology] -->
<RequiredPlugin pluginName="Sofa.Component.Topology.Mapping"/> <!-- Needed to use components [Hexa2TetraTopologicalMapping] -->
<RequiredPlugin pluginName="Sofa.Component.Visual"/> <!-- Needed to use components [LineAxis,VisualGrid,VisualStyle] -->
<RequiredPlugin pluginName="Sofa.GL.Component.Rendering3D"/> <!-- Needed to use components [OglSceneFrame] -->
</Node>

<DefaultAnimationLoop/>
<VisualStyle displayFlags="hideBehaviorModels showForceFields" />

<VisualGrid size="0.1"/>
<LineAxis size="0.1"/>

<EulerImplicitSolver name="backward_Euler" rayleighStiffness="0.1" rayleighMass="0.1" />

<RegularGridTopology name="grid" min="-0.05 -0.01 -0.01" max="0.05 0.01 0.01" n="20 5 5"/>
<MechanicalObject template="Vec3" name="state" showObject="true"/>

<NodalMassDensity property="1100"/>
<FEMMass template="Vec3,Hexahedron"/>

<BoxROI template="Vec3" name="box_roi_left" box="-0.051 -0.011 -0.011 -0.049 0.011 0.011" drawBoxes="1" />
<LinearMovementProjectiveConstraint template="Vec3" keyTimes="0 1 3 4 1e12" movements="0 0 0 -0.02 0 0 0.02 0 0 0 0 0 0 0 0" indices="@box_roi_left.indices" />

<BoxROI template="Vec3" name="box_roi_right" box="0.049 -0.011 -0.011 0.051 0.011 0.011" drawBoxes="1" />
<LinearMovementProjectiveConstraint template="Vec3" keyTimes="0 1 3 4 1e12" movements="0 0 0 0.02 0 0 -0.02 0 0 0 0 0 0 0 0" indices="@box_roi_right.indices" />

</Node>
15 changes: 15 additions & 0 deletions ...sion_compression/hexahedron/corotational/parallel/assembled/CorotationalFEMForceField.scn
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,15 @@
<?xml version="1.0"?>
<Node name="root" dt="0.01" gravity="0 0 0">
<include href="../../../../beam.xml"/>

<ConstantSparsityPatternSystem template="CompressedRowSparseMatrix" name="A" checkIndices="false"/>
<NaturalOrderingMethod/>
<SparseLDLSolver name="linear_solver" template="CompressedRowSparseMatrix"/>

<CorotationalFEMForceField name="FEM" template="Vec3,Hexahedron"
youngModulus="2e6" poissonRatio="0.45" topology="@grid"
computeForceStrategy="parallel" computeForceDerivStrategy="parallel"/>
<VonMisesStress template="Vec3,Hexahedron" name="stress" topology="@grid" stressEvaluator="@FEM"
colorMap="green yellow orange red purple #221C35"/>

</Node>
13 changes: 13 additions & 0 deletions ...ion_compression/hexahedron/corotational/parallel/matrixfree/CorotationalFEMForceField.scn
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,13 @@
<?xml version="1.0"?>
<Node name="root" dt="0.01" gravity="0 0 0">
<include href="../../../../beam.xml"/>

<CGLinearSolver iterations="250" name="linear_solver" tolerance="1.0e-12" threshold="1.0e-12" />

<CorotationalFEMForceField name="FEM" template="Vec3,Hexahedron"
youngModulus="2e6" poissonRatio="0.45" topology="@grid"
computeForceStrategy="parallel" computeForceDerivStrategy="parallel"/>
<VonMisesStress template="Vec3,Hexahedron" name="stress" topology="@grid" stressEvaluator="@FEM"
colorMap="green yellow orange red purple #221C35"/>

</Node>
15 changes: 15 additions & 0 deletions ...on_compression/hexahedron/corotational/sequential/assembled/CorotationalFEMForceField.scn
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,15 @@
<?xml version="1.0"?>
<Node name="root" dt="0.01" gravity="0 0 0">
<include href="../../../../beam.xml"/>

<ConstantSparsityPatternSystem template="CompressedRowSparseMatrix" name="A" checkIndices="false"/>
<NaturalOrderingMethod/>
<SparseLDLSolver name="linear_solver" template="CompressedRowSparseMatrix"/>

<CorotationalFEMForceField name="FEM" template="Vec3,Hexahedron"
youngModulus="2e6" poissonRatio="0.45" topology="@grid"
computeForceStrategy="sequenced" computeForceDerivStrategy="sequenced"/>
<VonMisesStress template="Vec3,Hexahedron" name="stress" topology="@grid" stressEvaluator="@FEM"
colorMap="green yellow orange red purple #221C35"/>

</Node>
13 changes: 13 additions & 0 deletions ...n_compression/hexahedron/corotational/sequential/matrixfree/CorotationalFEMForceField.scn
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,13 @@
<?xml version="1.0"?>
<Node name="root" dt="0.01" gravity="0 0 0">
<include href="../../../../beam.xml"/>

<CGLinearSolver iterations="250" name="linear_solver" tolerance="1.0e-12" threshold="1.0e-12" />

<CorotationalFEMForceField name="FEM" template="Vec3,Hexahedron"
youngModulus="2e6" poissonRatio="0.45" topology="@grid"
computeForceStrategy="sequenced" computeForceDerivStrategy="sequenced"/>
<VonMisesStress template="Vec3,Hexahedron" name="stress" topology="@grid" stressEvaluator="@FEM"
colorMap="green yellow orange red purple #221C35"/>

</Node>
15 changes: 15 additions & 0 deletions ...nsion_compression/hexahedron/linear/parallel/assembled/LinearSmallStrainFEMForceField.scn
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,15 @@
<?xml version="1.0"?>
<Node name="root" dt="0.01" gravity="0 0 0">
<include href="../../../../beam.xml"/>

<ConstantSparsityPatternSystem template="CompressedRowSparseMatrix" name="A" checkIndices="false"/>
<NaturalOrderingMethod/>
<SparseLDLSolver name="linear_solver" template="CompressedRowSparseMatrix"/>

<LinearSmallStrainFEMForceField name="FEM" template="Vec3,Hexahedron"
youngModulus="2e6" poissonRatio="0.45" topology="@grid"
computeForceStrategy="parallel" computeForceDerivStrategy="parallel"/>
<VonMisesStress template="Vec3,Hexahedron" name="stress" topology="@grid" stressEvaluator="@FEM"
colorMap="green yellow orange red purple #221C35"/>

</Node>
13 changes: 13 additions & 0 deletions ...sion_compression/hexahedron/linear/parallel/matrixfree/LinearSmallStrainFEMForceField.scn
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,13 @@
<?xml version="1.0"?>
<Node name="root" dt="0.01" gravity="0 0 0">
<include href="../../../../beam.xml"/>

<CGLinearSolver iterations="250" name="linear_solver" tolerance="1.0e-12" threshold="1.0e-12" />

<LinearSmallStrainFEMForceField name="FEM" template="Vec3,Hexahedron"
youngModulus="2e6" poissonRatio="0.45" topology="@grid"
computeForceStrategy="parallel" computeForceDerivStrategy="parallel"/>
<VonMisesStress template="Vec3,Hexahedron" name="stress" topology="@grid" stressEvaluator="@FEM"
colorMap="green yellow orange red purple #221C35"/>

</Node>
15 changes: 15 additions & 0 deletions ...ion_compression/hexahedron/linear/sequential/assembled/LinearSmallStrainFEMForceField.scn
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,15 @@
<?xml version="1.0"?>
<Node name="root" dt="0.01" gravity="0 0 0">
<include href="../../../../beam.xml"/>

<ConstantSparsityPatternSystem template="CompressedRowSparseMatrix" name="A" checkIndices="false"/>
<NaturalOrderingMethod/>
<SparseLDLSolver name="linear_solver" template="CompressedRowSparseMatrix"/>

<LinearSmallStrainFEMForceField name="FEM" template="Vec3,Hexahedron"
youngModulus="2e6" poissonRatio="0.45" topology="@grid"
computeForceStrategy="sequenced" computeForceDerivStrategy="sequenced"/>
<VonMisesStress template="Vec3,Hexahedron" name="stress" topology="@grid" stressEvaluator="@FEM"
colorMap="green yellow orange red purple #221C35"/>

</Node>
13 changes: 13 additions & 0 deletions ...on_compression/hexahedron/linear/sequential/matrixfree/LinearSmallStrainFEMForceField.scn
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,13 @@
<?xml version="1.0"?>
<Node name="root" dt="0.01" gravity="0 0 0">
<include href="../../../../beam.xml"/>

<CGLinearSolver iterations="250" name="linear_solver" tolerance="1.0e-12" threshold="1.0e-12" />

<LinearSmallStrainFEMForceField name="FEM" template="Vec3,Hexahedron"
youngModulus="2e6" poissonRatio="0.45" topology="@grid"
computeForceStrategy="sequenced" computeForceDerivStrategy="sequenced"/>
<VonMisesStress template="Vec3,Hexahedron" name="stress" topology="@grid" stressEvaluator="@FEM"
colorMap="green yellow orange red purple #221C35"/>

</Node>
22 changes: 22 additions & 0 deletions ...ion_compression/tetrahedron/corotational/parallel/assembled/CorotationalFEMForceField.scn
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,22 @@
<?xml version="1.0"?>
<Node name="root" dt="0.01" gravity="0 0 0">
<include href="../../../../beam.xml"/>

<ConstantSparsityPatternSystem template="CompressedRowSparseMatrix" name="A" checkIndices="false"/>
<NaturalOrderingMethod/>
<SparseLDLSolver name="linear_solver" template="CompressedRowSparseMatrix"/>

<Node name="fem">
<TetrahedronSetTopologyContainer name="Tetra_topo"/>
<TetrahedronSetTopologyModifier name="Modifier" />
<TetrahedronSetGeometryAlgorithms template="Vec3" name="GeomAlgo" drawTetrahedra="false"/>
<Hexa2TetraTopologicalMapping input="@grid" output="@Tetra_topo" swapping="true"/>

<CorotationalFEMForceField name="FEM" template="Vec3,Tetrahedron"
youngModulus="2e6" poissonRatio="0.45" topology="@Tetra_topo"
computeForceStrategy="parallel" computeForceDerivStrategy="parallel"/>
<VonMisesStress template="Vec3,Tetrahedron" name="stress" topology="@Tetra_topo" stressEvaluator="@FEM"
colorMap="green yellow orange red purple #221C35"/>
</Node>

</Node>
20 changes: 20 additions & 0 deletions ...on_compression/tetrahedron/corotational/parallel/matrixfree/CorotationalFEMForceField.scn
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,20 @@
<?xml version="1.0"?>
<Node name="root" dt="0.01" gravity="0 0 0">
<include href="../../../../beam.xml"/>

<CGLinearSolver iterations="250" name="linear_solver" tolerance="1.0e-12" threshold="1.0e-12" />

<Node name="fem">
<TetrahedronSetTopologyContainer name="Tetra_topo"/>
<TetrahedronSetTopologyModifier name="Modifier" />
<TetrahedronSetGeometryAlgorithms template="Vec3" name="GeomAlgo" drawTetrahedra="false"/>
<Hexa2TetraTopologicalMapping input="@grid" output="@Tetra_topo" swapping="true"/>

<CorotationalFEMForceField name="FEM" template="Vec3,Tetrahedron"
youngModulus="2e6" poissonRatio="0.45" topology="@Tetra_topo"
computeForceStrategy="parallel" computeForceDerivStrategy="parallel"/>
<VonMisesStress template="Vec3,Tetrahedron" name="stress" topology="@Tetra_topo" stressEvaluator="@FEM"
colorMap="green yellow orange red purple #221C35"/>
</Node>

</Node>
22 changes: 22 additions & 0 deletions ...n_compression/tetrahedron/corotational/sequential/assembled/CorotationalFEMForceField.scn
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,22 @@
<?xml version="1.0"?>
<Node name="root" dt="0.01" gravity="0 0 0">
<include href="../../../../beam.xml"/>

<ConstantSparsityPatternSystem template="CompressedRowSparseMatrix" name="A" checkIndices="false"/>
<NaturalOrderingMethod/>
<SparseLDLSolver name="linear_solver" template="CompressedRowSparseMatrix"/>

<Node name="fem">
<TetrahedronSetTopologyContainer name="Tetra_topo"/>
<TetrahedronSetTopologyModifier name="Modifier" />
<TetrahedronSetGeometryAlgorithms template="Vec3" name="GeomAlgo" drawTetrahedra="false"/>
<Hexa2TetraTopologicalMapping input="@grid" output="@Tetra_topo" swapping="true"/>

<CorotationalFEMForceField name="FEM" template="Vec3,Tetrahedron"
youngModulus="2e6" poissonRatio="0.45" topology="@Tetra_topo"
computeForceStrategy="sequenced" computeForceDerivStrategy="sequenced"/>
<VonMisesStress template="Vec3,Tetrahedron" name="stress" topology="@Tetra_topo" stressEvaluator="@FEM"
colorMap="green yellow orange red purple #221C35"/>
</Node>

</Node>
20 changes: 20 additions & 0 deletions ..._compression/tetrahedron/corotational/sequential/matrixfree/CorotationalFEMForceField.scn
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,20 @@
<?xml version="1.0"?>
<Node name="root" dt="0.01" gravity="0 0 0">
<include href="../../../../beam.xml"/>

<CGLinearSolver iterations="250" name="linear_solver" tolerance="1.0e-12" threshold="1.0e-12" />

<Node name="fem">
<TetrahedronSetTopologyContainer name="Tetra_topo"/>
<TetrahedronSetTopologyModifier name="Modifier" />
<TetrahedronSetGeometryAlgorithms template="Vec3" name="GeomAlgo" drawTetrahedra="false"/>
<Hexa2TetraTopologicalMapping input="@grid" output="@Tetra_topo" swapping="true"/>

<CorotationalFEMForceField name="FEM" template="Vec3,Tetrahedron"
youngModulus="2e6" poissonRatio="0.45" topology="@Tetra_topo"
computeForceStrategy="sequenced" computeForceDerivStrategy="sequenced"/>
<VonMisesStress template="Vec3,Tetrahedron" name="stress" topology="@Tetra_topo" stressEvaluator="@FEM"
colorMap="green yellow orange red purple #221C35"/>
</Node>

</Node>
22 changes: 22 additions & 0 deletions ...sion_compression/tetrahedron/linear/parallel/assembled/LinearSmallStrainFEMForceField.scn
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,22 @@
<?xml version="1.0"?>
<Node name="root" dt="0.01" gravity="0 0 0">
<include href="../../../../beam.xml"/>

<ConstantSparsityPatternSystem template="CompressedRowSparseMatrix" name="A" checkIndices="false"/>
<NaturalOrderingMethod/>
<SparseLDLSolver name="linear_solver" template="CompressedRowSparseMatrix"/>

<Node name="fem">
<TetrahedronSetTopologyContainer name="Tetra_topo"/>
<TetrahedronSetTopologyModifier name="Modifier" />
<TetrahedronSetGeometryAlgorithms template="Vec3" name="GeomAlgo" drawTetrahedra="false"/>
<Hexa2TetraTopologicalMapping input="@grid" output="@Tetra_topo" swapping="true"/>

<LinearSmallStrainFEMForceField name="FEM" template="Vec3,Tetrahedron"
youngModulus="2e6" poissonRatio="0.45" topology="@Tetra_topo"
computeForceStrategy="parallel" computeForceDerivStrategy="parallel"/>
<VonMisesStress template="Vec3,Tetrahedron" name="stress" topology="@Tetra_topo" stressEvaluator="@FEM"
colorMap="green yellow orange red purple #221C35"/>
</Node>

</Node>
20 changes: 20 additions & 0 deletions ...ion_compression/tetrahedron/linear/parallel/matrixfree/LinearSmallStrainFEMForceField.scn
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,20 @@
<?xml version="1.0"?>
<Node name="root" dt="0.01" gravity="0 0 0">
<include href="../../../../beam.xml"/>

<CGLinearSolver iterations="250" name="linear_solver" tolerance="1.0e-12" threshold="1.0e-12" />

<Node name="fem">
<TetrahedronSetTopologyContainer name="Tetra_topo"/>
<TetrahedronSetTopologyModifier name="Modifier" />
<TetrahedronSetGeometryAlgorithms template="Vec3" name="GeomAlgo" drawTetrahedra="false"/>
<Hexa2TetraTopologicalMapping input="@grid" output="@Tetra_topo" swapping="true"/>

<LinearSmallStrainFEMForceField name="FEM" template="Vec3,Tetrahedron"
youngModulus="2e6" poissonRatio="0.45" topology="@Tetra_topo"
computeForceStrategy="parallel" computeForceDerivStrategy="parallel"/>
<VonMisesStress template="Vec3,Tetrahedron" name="stress" topology="@Tetra_topo" stressEvaluator="@FEM"
colorMap="green yellow orange red purple #221C35"/>
</Node>

</Node>
22 changes: 22 additions & 0 deletions ...on_compression/tetrahedron/linear/sequential/assembled/LinearSmallStrainFEMForceField.scn
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,22 @@
<?xml version="1.0"?>
<Node name="root" dt="0.01" gravity="0 0 0">
<include href="../../../../beam.xml"/>

<ConstantSparsityPatternSystem template="CompressedRowSparseMatrix" name="A" checkIndices="false"/>
<NaturalOrderingMethod/>
<SparseLDLSolver name="linear_solver" template="CompressedRowSparseMatrix"/>

<Node name="fem">
<TetrahedronSetTopologyContainer name="Tetra_topo"/>
<TetrahedronSetTopologyModifier name="Modifier" />
<TetrahedronSetGeometryAlgorithms template="Vec3" name="GeomAlgo" drawTetrahedra="false"/>
<Hexa2TetraTopologicalMapping input="@grid" output="@Tetra_topo" swapping="true"/>

<LinearSmallStrainFEMForceField name="FEM" template="Vec3,Tetrahedron"
youngModulus="2e6" poissonRatio="0.45" topology="@Tetra_topo"
computeForceStrategy="sequenced" computeForceDerivStrategy="sequenced"/>
<VonMisesStress template="Vec3,Tetrahedron" name="stress" topology="@Tetra_topo" stressEvaluator="@FEM"
colorMap="green yellow orange red purple #221C35"/>
</Node>

</Node>
20 changes: 20 additions & 0 deletions ...n_compression/tetrahedron/linear/sequential/matrixfree/LinearSmallStrainFEMForceField.scn
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,20 @@
<?xml version="1.0"?>
<Node name="root" dt="0.01" gravity="0 0 0">
<include href="../../../../beam.xml"/>

<CGLinearSolver iterations="250" name="linear_solver" tolerance="1.0e-12" threshold="1.0e-12" />

<Node name="fem">
<TetrahedronSetTopologyContainer name="Tetra_topo"/>
<TetrahedronSetTopologyModifier name="Modifier" />
<TetrahedronSetGeometryAlgorithms template="Vec3" name="GeomAlgo" drawTetrahedra="false"/>
<Hexa2TetraTopologicalMapping input="@grid" output="@Tetra_topo" swapping="true"/>

<LinearSmallStrainFEMForceField name="FEM" template="Vec3,Tetrahedron"
youngModulus="2e6" poissonRatio="0.45" topology="@Tetra_topo"
computeForceStrategy="sequenced" computeForceDerivStrategy="sequenced"/>
<VonMisesStress template="Vec3,Tetrahedron" name="stress" topology="@Tetra_topo" stressEvaluator="@FEM"
colorMap="green yellow orange red purple #221C35"/>
</Node>

</Node>
Loading