-
Notifications
You must be signed in to change notification settings - Fork 311
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
[SofaMiscFEM] ADD: new scene for FastTetrahedralCorotationalFEM, and …
…rename the old one as _validation as it compares to TetrahedralCororationalFEM.
- Loading branch information
Showing
2 changed files
with
62 additions
and
32 deletions.
There are no files selected for viewing
53 changes: 21 additions & 32 deletions
53
examples/Components/forcefield/FastTetrahedronCorotationalForceField.scn
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -1,41 +1,30 @@ | ||
<?xml version="1.0" ?> | ||
<Node name="root" dt="2" showBoundingTree="0" gravity="0 0 0"> | ||
<Node name="root" dt="0.05" showBoundingTree="0" gravity="0 0 0"> | ||
<VisualStyle displayFlags="showBehaviorModels showForceFields showCollisionModels showVisual" /> | ||
<CollisionPipeline verbose="0" name="CollisionPipeline" /> | ||
<BruteForceDetection name="N2" /> | ||
<CollisionResponse response="default" name="collision response" /> | ||
<DiscreteIntersection /> | ||
<Node name="MeshTopology"> | ||
<TetrahedronSetTopologyContainer name="Container1" fileTopology="mesh/cylinder.msh" /> | ||
<TetrahedronSetGeometryAlgorithms name="GeomAlgo" /> | ||
<MechanicalObject name="dofs" /> | ||
<Node name="FastTetrahedronCorotationalForceField"> | ||
<CGLinearSolver iterations="500" name="linear solver" tolerance="1.0e-9" threshold="1.0e-9" /> | ||
<StaticSolver /> | ||
<TetrahedronSetTopologyContainer name="Container2" /> | ||
<TetrahedronSetGeometryAlgorithms drawEdges="1" drawColorEdges="0 1 0" name="GeomAlgo" /> | ||
<SimpleTesselatedTetraTopologicalMapping input="@Container1" output="@Container2" /> | ||
<MechanicalObject translation="0 0 0" /> | ||
<MeshMatrixMass name="mass" lumping="1" printMass="0" massDensity="1" /> | ||
<BoxROI box="-1.2 -1.2 -0.01 1.2 1.2 0.01" drawBoxes="1" name="fixedPlane2" /> | ||
<FixedConstraint indices="@fixedPlane2.indices" /> | ||
<BoxROI box="-1.2 -1.2 0.99 1.2 1.2 1.01" drawBoxes="0" name="pressurePlane2" /> | ||
<ConstantForceField points="@pressurePlane2.indices" totalForce="0.05 0 -0.2" /> | ||
<FastTetrahedralCorotationalForceField poissonRatio="0.45" youngModulus="30" method="large" /> | ||
</Node> | ||
<Node name="TetrahedronCorotationalForceField"> | ||
<CGLinearSolver iterations="500" name="linear solver" tolerance="1.0e-9" threshold="1.0e-9" /> | ||
<StaticSolver /> | ||
<TetrahedronSetTopologyContainer name="Container3" /> | ||
<TetrahedronSetGeometryAlgorithms drawEdges="1" drawColorEdges="1 0 0" name="GeomAlgo" /> | ||
<SimpleTesselatedTetraTopologicalMapping input="@Container1" output="@Container3" /> | ||
<MechanicalObject translation="0 0 0" /> | ||
<MeshMatrixMass name="mass" lumping="1" printMass="0" massDensity="1" /> | ||
<BoxROI box="-1.2 -1.2 -0.01 1.2 1.2 0.01" drawBoxes="1" name="fixedPlane" /> | ||
<FixedConstraint indices="@fixedPlane.indices" /> | ||
<BoxROI box="-1.2 -1.2 0.99 1.2 1.2 1.01" drawBoxes="0" name="pressurePlane" /> | ||
<ConstantForceField points="@pressurePlane.indices" totalForce="0.05 0 -0.2" /> | ||
<TetrahedralCorotationalFEMForceField poissonRatio="0.45" youngModulus="30" method="large" /> | ||
<Node name="TT_FAST"> | ||
<EulerImplicit name="cg_odesolver" printLog="false" rayleighStiffness="0.1" rayleighMass="0.1" /> | ||
<CGLinearSolver iterations="25" name="linear solver" tolerance="1.0e-9" threshold="1.0e-9" /> | ||
<MeshGmshLoader name="loader" filename="mesh/cylinder_8x30x6.msh" /> | ||
<MechanicalObject src="@loader" name="Volume" translation="1 0 0"/> | ||
<include href="Objects/TetrahedronSetTopology.xml" src="@loader" tags=" " /> | ||
<DiagonalMass massDensity="0.5" /> | ||
<FixedPlaneConstraint direction="0 0 1" dmin="-0.1" dmax="0.1" /> | ||
<FixedConstraint indices="0" /> | ||
<FastTetrahedralCorotationalForceField poissonRatio="0.3" youngModulus="3600" method="large" /> | ||
|
||
<Node name="T_FAST"> | ||
<include href="Objects/TriangleSetTopology.xml" src="@" tags=" " /> | ||
<Tetra2TriangleTopologicalMapping input="@../Container" output="@Container" /> | ||
<TrianglePressureForceField normal="0 0 1" dmin="0.9" dmax="1.1" pressure="0.4 0 0" /> | ||
<TriangleSet /> | ||
<Node name="Visu"> | ||
<OglModel name="Visual" color="blue" /> | ||
<IdentityMapping input="@../../Volume" output="@Visual" /> | ||
</Node> | ||
</Node> | ||
</Node> | ||
</Node> |
41 changes: 41 additions & 0 deletions
41
examples/Components/forcefield/FastTetrahedronCorotationalForceField_validation.scn
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -0,0 +1,41 @@ | ||
<?xml version="1.0" ?> | ||
<Node name="root" dt="2" showBoundingTree="0" gravity="0 0 0"> | ||
<VisualStyle displayFlags="showBehaviorModels showForceFields showCollisionModels showVisual" /> | ||
<CollisionPipeline verbose="0" name="CollisionPipeline" /> | ||
<BruteForceDetection name="N2" /> | ||
<CollisionResponse response="default" name="collision response" /> | ||
<DiscreteIntersection /> | ||
<Node name="MeshTopology"> | ||
<TetrahedronSetTopologyContainer name="Container1" fileTopology="mesh/cylinder.msh" /> | ||
<TetrahedronSetGeometryAlgorithms name="GeomAlgo" /> | ||
<MechanicalObject name="dofs" /> | ||
<Node name="FastTetrahedronCorotationalForceField"> | ||
<CGLinearSolver iterations="500" name="linear solver" tolerance="1.0e-9" threshold="1.0e-9" /> | ||
<StaticSolver /> | ||
<TetrahedronSetTopologyContainer name="Container2" /> | ||
<TetrahedronSetGeometryAlgorithms drawEdges="1" drawColorEdges="0 1 0" name="GeomAlgo" /> | ||
<SimpleTesselatedTetraTopologicalMapping input="@Container1" output="@Container2" /> | ||
<MechanicalObject translation="0 0 0" /> | ||
<MeshMatrixMass name="mass" lumping="1" printMass="0" massDensity="1" /> | ||
<BoxROI box="-1.2 -1.2 -0.01 1.2 1.2 0.01" drawBoxes="1" name="fixedPlane2" /> | ||
<FixedConstraint indices="@fixedPlane2.indices" /> | ||
<BoxROI box="-1.2 -1.2 0.99 1.2 1.2 1.01" drawBoxes="0" name="pressurePlane2" /> | ||
<ConstantForceField points="@pressurePlane2.indices" totalForce="0.05 0 -0.2" /> | ||
<FastTetrahedralCorotationalForceField poissonRatio="0.45" youngModulus="30" method="large" /> | ||
</Node> | ||
<Node name="TetrahedronCorotationalForceField"> | ||
<CGLinearSolver iterations="500" name="linear solver" tolerance="1.0e-9" threshold="1.0e-9" /> | ||
<StaticSolver /> | ||
<TetrahedronSetTopologyContainer name="Container3" /> | ||
<TetrahedronSetGeometryAlgorithms drawEdges="1" drawColorEdges="1 0 0" name="GeomAlgo" /> | ||
<SimpleTesselatedTetraTopologicalMapping input="@Container1" output="@Container3" /> | ||
<MechanicalObject translation="0 0 0" /> | ||
<MeshMatrixMass name="mass" lumping="1" printMass="0" massDensity="1" /> | ||
<BoxROI box="-1.2 -1.2 -0.01 1.2 1.2 0.01" drawBoxes="1" name="fixedPlane" /> | ||
<FixedConstraint indices="@fixedPlane.indices" /> | ||
<BoxROI box="-1.2 -1.2 0.99 1.2 1.2 1.01" drawBoxes="0" name="pressurePlane" /> | ||
<ConstantForceField points="@pressurePlane.indices" totalForce="0.05 0 -0.2" /> | ||
<TetrahedralCorotationalFEMForceField poissonRatio="0.45" youngModulus="30" method="large" /> | ||
</Node> | ||
</Node> | ||
</Node> |