Add HyperreducedFixedWeakConstraint

bakpaul · bakpaul · commit f2e8974dde96 · 2025-04-03T15:35:44.000+02:00
diff --git a/CMakeLists.txt b/CMakeLists.txt
@@ -27,7 +27,8 @@ set(HEADER_FILES
     ${MOR_SRC}/component/forcefield/HyperReducedHexahedronFEMForceField.h
     ${MOR_SRC}/component/forcefield/HyperReducedHexahedronFEMForceField.inl
     ${MOR_SRC}/component/forcefield/HyperReducedRestShapeSpringsForceField.h
-    ${MOR_SRC}/component/forcefield/HyperReducedRestShapeSpringsForceField.inl
+    ${MOR_SRC}/component/forcefield/HyperReducedFixedWeakConstraint.h
+    ${MOR_SRC}/component/forcefield/HyperReducedFixedWeakConstraint.inl
     ${MOR_SRC}/component/forcefield/HyperReducedTetrahedralCorotationalFEMForceField.h
     ${MOR_SRC}/component/forcefield/HyperReducedTetrahedralCorotationalFEMForceField.inl
     ${MOR_SRC}/component/forcefield/HyperReducedTetrahedronFEMForceField.h
@@ -58,6 +59,7 @@ set(SOURCE_FILES
     ${MOR_SRC}/component/forcefield/HyperReducedHelper.cpp
     ${MOR_SRC}/component/forcefield/HyperReducedHexahedronFEMForceField.cpp
     ${MOR_SRC}/component/forcefield/HyperReducedRestShapeSpringsForceField.cpp
+    ${MOR_SRC}/component/forcefield/HyperReducedFixedWeakConstraint.cpp
     ${MOR_SRC}/component/forcefield/HyperReducedTetrahedralCorotationalFEMForceField.cpp
     ${MOR_SRC}/component/forcefield/HyperReducedTetrahedronFEMForceField.cpp
     ${MOR_SRC}/component/forcefield/HyperReducedTetrahedronHyperelasticityFEMForceField.cpp
diff --git a/examples/organs/liver/reduced/liver/reduced_liverFine.py b/examples/organs/liver/reduced/liver/reduced_liverFine.py
@@ -77,7 +77,7 @@ def Reduced_test(
     liver.addObject('BoxROI' , name= 'boxROIactuation' , orientedBox= newBox([[-5.0, 0.5, -0.5], [-5.0, 0.0, -0.5], [-4.0, 0.0, -0.5]] , [0.0, 0.0, 0.0],translation,rotation,[0, 0, 0.5],scale) + multiply(scale[2],[1.0]).tolist(),drawBoxes=True)
     liver.addObject('UniformMass' , totalMass = 0.3)
     liver.addObject('HyperReducedTetrahedronFEMForceField' , poissonRatio = '0.3', youngModulus = '5000', name = 'reducedFF_liver_0', nbModes = nbrOfModes, performECSW = hyperReduction, modesPath = path + r'/data/modes.txt', RIDPath = path + r'/data/reducedFF_liver_0_RID.txt', weightsPath = path + r'/data/reducedFF_liver_0_weight.txt')
-    liver.addObject('HyperReducedRestShapeSpringsForceField' , indices = '@ROI1.indices', stiffness = '1e8', name = 'reducedFF_liver_1', nbModes = nbrOfModes, performECSW = hyperReduction, modesPath = path + r'/data/modes.txt', RIDPath = path + r'/data/reducedFF_liver_1_RID.txt', weightsPath = path + r'/data/reducedFF_liver_1_weight.txt')
+    liver.addObject('HyperReducedFixedWeakConstraint' , indices = '@ROI1.indices', stiffness = '1e8', name = 'reducedFF_liver_1', nbModes = nbrOfModes, performECSW = hyperReduction, modesPath = path + r'/data/modes.txt', RIDPath = path + r'/data/reducedFF_liver_1_RID.txt', weightsPath = path + r'/data/reducedFF_liver_1_weight.txt')
     liver.addObject('HyperReducedRestShapeSpringsForceField' , externalIndices = [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10], indices = '@boxROIactuation.indices', name = 'reducedFF_liver_2', stiffness = '1e8', external_rest_shape = '@actuator/actuatorState', nbModes = nbrOfModes, performECSW = hyperReduction, modesPath = path + r'/data/modes.txt', RIDPath = path + r'/data/reducedFF_liver_2_RID.txt', weightsPath = path + r'/data/reducedFF_liver_2_weight.txt')
     liver.addObject('ModelOrderReductionMapping' , input = '@../MechanicalObject', modesPath = path + r'/data/modes.txt', output = '@./MechanicalObject')
 
diff --git a/examples/organs/liver/reduced/liver_hyperElastic/reduced_liverFineHyperElastic.py b/examples/organs/liver/reduced/liver_hyperElastic/reduced_liverFineHyperElastic.py
@@ -77,7 +77,7 @@ def Reduced_test(
     liver.addObject('BoxROI' , name= 'boxROIactuation' , orientedBox= newBox([[-5.0, 0.5, -0.5], [-5.0, 0.0, -0.5], [-4.0, 0.0, -0.5]] , [0.0, 0.0, 0.0],translation,rotation,[0, 0, 0.5],scale) + multiply(scale[2],[1.0]).tolist(),drawBoxes=True)
     liver.addObject('UniformMass' , totalMass = 0.3)
     liver.addObject('HyperReducedTetrahedronHyperelasticityFEMForceField' , materialName = 'NeoHookean', ParameterSet = '1923.076923076923 4166.666666666666', AnisotropyDirections = '', name = 'reducedFF_liver_0', nbModes = nbrOfModes, performECSW = hyperReduction, modesPath = path + r'/data/modes.txt', RIDPath = path + r'/data/reducedFF_liver_0_RID.txt', weightsPath = path + r'/data/reducedFF_liver_0_weight.txt')
-    liver.addObject('HyperReducedRestShapeSpringsForceField' , indices = '@ROI1.indices', stiffness = '1e8', name = 'reducedFF_liver_1', nbModes = nbrOfModes, performECSW = hyperReduction, modesPath = path + r'/data/modes.txt', RIDPath = path + r'/data/reducedFF_liver_1_RID.txt', weightsPath = path + r'/data/reducedFF_liver_1_weight.txt')
+    liver.addObject('HyperReducedFixedWeakConstraint' , indices = '@ROI1.indices', stiffness = '1e8', name = 'reducedFF_liver_1', nbModes = nbrOfModes, performECSW = hyperReduction, modesPath = path + r'/data/modes.txt', RIDPath = path + r'/data/reducedFF_liver_1_RID.txt', weightsPath = path + r'/data/reducedFF_liver_1_weight.txt')
     liver.addObject('HyperReducedRestShapeSpringsForceField' , externalIndices = [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10], indices = '@boxROIactuation.indices', name = 'reducedFF_liver_2', stiffness = '1e8', external_rest_shape = '@actuator/actuatorState', nbModes = nbrOfModes, performECSW = hyperReduction, modesPath = path + r'/data/modes.txt', RIDPath = path + r'/data/reducedFF_liver_2_RID.txt', weightsPath = path + r'/data/reducedFF_liver_2_weight.txt')
     liver.addObject('ModelOrderReductionMapping' , input = '@../MechanicalObject', modesPath = path + r'/data/modes.txt', output = '@./MechanicalObject')
 
diff --git a/examples/others/hexaBeam/reduced/reduced_hexaBeam.py b/examples/others/hexaBeam/reduced/reduced_hexaBeam.py
@@ -73,7 +73,7 @@ def Reduced_test(
     M1.addObject('MechanicalObject')
     M1.addObject('UniformMass' , totalMass = '0.1')
     M1.addObject('RegularGridTopology' , nx = '4', ny = '4', nz = '20', xmin = '-9', xmax = '-6', ymin = '0', ymax = '3', zmin = '0', zmax = '19')
-    M1.addObject('HyperReducedRestShapeSpringsForceField' , indices = '0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15', stiffness = 1000000000000.0, name = 'reducedFF_M1_0', nbModes = nbrOfModes, performECSW = hyperReduction, modesPath = path + r'/data/modes.txt', RIDPath = path + r'/data/reducedFF_M1_0_RID.txt', weightsPath = path + r'/data/reducedFF_M1_0_weight.txt')
+    M1.addObject('HyperReducedFixedWeakConstraint' , indices = '0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15', stiffness = 1000000000000.0, name = 'reducedFF_M1_0', nbModes = nbrOfModes, performECSW = hyperReduction, modesPath = path + r'/data/modes.txt', RIDPath = path + r'/data/reducedFF_M1_0_RID.txt', weightsPath = path + r'/data/reducedFF_M1_0_weight.txt')
     M1.addObject('HyperReducedHexahedronFEMForceField' , name = 'reducedFF_M1_1', youngModulus = '4000', poissonRatio = '0.3', method = 'large', nbModes = nbrOfModes, performECSW = hyperReduction, modesPath = path + r'/data/modes.txt', RIDPath = path + r'/data/reducedFF_M1_1_RID.txt', weightsPath = path + r'/data/reducedFF_M1_1_weight.txt')
     M1.addObject('ModelOrderReductionMapping' , input = '@../MechanicalObject', modesPath = path + r'/data/modes.txt', output = '@./MechanicalObject')
 
diff --git a/examples/softRobots/diamond/reduced/diamond/reduced_diamondRobot.py b/examples/softRobots/diamond/reduced/diamond/reduced_diamondRobot.py
@@ -87,7 +87,7 @@ def Reduced_test(
 
     FixedBox = modelNode.addChild('FixedBox')
     FixedBox.addObject('BoxROI' , name= 'BoxROI' , orientedBox= newBox([[-15, 15, -40], [-15, -15, -40], [15, -15, -40]] , [0.0, 0.0, 35],translation,rotation,[0, 0, 25.0],scale) + multiply(scale[2],[50]).tolist(),drawBoxes=True)
-    FixedBox.addObject('HyperReducedRestShapeSpringsForceField' , indices = '@BoxROI.indices', stiffness = 1000000000000.0, name = 'reducedFF_FixedBox_1', nbModes = nbrOfModes, performECSW = hyperReduction, modesPath = path + r'/data/modes.txt', RIDPath = path + r'/data/reducedFF_FixedBox_1_RID.txt', weightsPath = path + r'/data/reducedFF_FixedBox_1_weight.txt')
+    FixedBox.addObject('HyperReducedFixedWeakConstraint' , indices = '@BoxROI.indices', stiffness = 1000000000000.0, name = 'reducedFF_FixedBox_1', nbModes = nbrOfModes, performECSW = hyperReduction, modesPath = path + r'/data/modes.txt', RIDPath = path + r'/data/reducedFF_FixedBox_1_RID.txt', weightsPath = path + r'/data/reducedFF_FixedBox_1_weight.txt')
 
 
     north = modelNode.addChild('north')
diff --git a/examples/softRobots/diamond/reduced/diamond_hyperElastic/reduced_diamondRobotHyperElastic.py b/examples/softRobots/diamond/reduced/diamond_hyperElastic/reduced_diamondRobotHyperElastic.py
@@ -81,7 +81,7 @@ def Reduced_test(
 
     FixedBox = modelNode.addChild('FixedBox')
     FixedBox.addObject('BoxROI' , name= 'BoxROI' , orientedBox= newBox([[-15, 15, -40], [-15, -15, -40], [15, -15, -40]] , [0.0, 0.0, 35],translation,rotation,[0, 0, 25.0],scale) + multiply(scale[2],[50]).tolist(),drawBoxes=True)
-    FixedBox.addObject('HyperReducedRestShapeSpringsForceField' , indices = '@BoxROI.indices', stiffness = 1000000000000.0, name = 'reducedFF_FixedBox_1', nbModes = nbrOfModes, performECSW = hyperReduction, modesPath = path + r'/data/modes.txt', RIDPath = path + r'/data/reducedFF_FixedBox_1_RID.txt', weightsPath = path + r'/data/reducedFF_FixedBox_1_weight.txt')
+    FixedBox.addObject('HyperReducedFixedWeakConstraint' , indices = '@BoxROI.indices', stiffness = 1000000000000.0, name = 'reducedFF_FixedBox_1', nbModes = nbrOfModes, performECSW = hyperReduction, modesPath = path + r'/data/modes.txt', RIDPath = path + r'/data/reducedFF_FixedBox_1_RID.txt', weightsPath = path + r'/data/reducedFF_FixedBox_1_weight.txt')
 
 
     north = modelNode.addChild('north')
diff --git a/examples/softRobots/finger/reduced/reduced_finger.py b/examples/softRobots/finger/reduced/reduced_finger.py
@@ -77,7 +77,7 @@ def Reduced_test(
     finger.addObject('UniformMass' , totalMass = 0.05)
     finger.addObject('HyperReducedTetrahedronFEMForceField' , poissonRatio = '0.45', youngModulus = '600', name = 'reducedFF_finger_0', nbModes = nbrOfModes, performECSW = hyperReduction, modesPath = path + r'/data/modes.txt', RIDPath = path + r'/data/reducedFF_finger_0_RID.txt', weightsPath = path + r'/data/reducedFF_finger_0_weight.txt')
     finger.addObject('BoxROI' , name= 'ROI1' , orientedBox= newBox([[-15, 10, 0], [-15, 0, 0], [5, 0, 0]] , [0.0, 0.0, 0.0],translation,rotation,[0, 0, 7.5],scale) + multiply(scale[2],[15]).tolist(),drawBoxes=True)
-    finger.addObject('HyperReducedRestShapeSpringsForceField' , indices = '@ROI1.indices', stiffness = 1000000000000.0, name = 'reducedFF_finger_1', nbModes = nbrOfModes, performECSW = hyperReduction, modesPath = path + r'/data/modes.txt', RIDPath = path + r'/data/reducedFF_finger_1_RID.txt', weightsPath = path + r'/data/reducedFF_finger_1_weight.txt')
+    finger.addObject('HyperReducedFixedWeakConstraint' , indices = '@ROI1.indices', stiffness = 1000000000000.0, name = 'reducedFF_finger_1', nbModes = nbrOfModes, performECSW = hyperReduction, modesPath = path + r'/data/modes.txt', RIDPath = path + r'/data/reducedFF_finger_1_RID.txt', weightsPath = path + r'/data/reducedFF_finger_1_weight.txt')
     finger.addObject('ModelOrderReductionMapping' , input = '@../MechanicalObject', modesPath = path + r'/data/modes.txt', output = '@./MechanicalObject')
 
 
diff --git a/examples/softRobots/sofiaLeg/reduced/reduced_sofiaLeg.py b/examples/softRobots/sofiaLeg/reduced/reduced_sofiaLeg.py
@@ -76,7 +76,7 @@ def Reduced_SofiaLeg(
     SofiaLeg.addObject('UniformMass' , totalMass = 0.01)
     SofiaLeg.addObject('HyperReducedTetrahedronFEMForceField' , youngModulus = 300, poissonRatio = 0.45, name = 'reducedFF_SofiaLeg_0', nbModes = nbrOfModes, performECSW = hyperReduction, modesPath = path + r'/data/modes.txt', RIDPath = path + r'/data/reducedFF_SofiaLeg_0_RID.txt', weightsPath = path + r'/data/reducedFF_SofiaLeg_0_weight.txt')
     SofiaLeg.addObject('BoxROI' , name = 'boxROITop', orientedBox = [[-12.0, 53.0, 0.0], [12.0, 53.0, 0.0], [12.0, 64.0, 0.0], 16.0], drawBoxes = True)
-    SofiaLeg.addObject('HyperReducedRestShapeSpringsForceField' , name = 'reducedFF_SofiaLeg_1', indices = '@boxROITop.indices', stiffness = '1e8', nbModes = nbrOfModes, performECSW = hyperReduction, modesPath = path + r'/data/modes.txt', RIDPath = path + r'/data/reducedFF_SofiaLeg_1_RID.txt', weightsPath = path + r'/data/reducedFF_SofiaLeg_1_weight.txt')
+    SofiaLeg.addObject('HyperReducedFixedWeakConstraint' , name = 'reducedFF_SofiaLeg_1', indices = '@boxROITop.indices', stiffness = '1e8', nbModes = nbrOfModes, performECSW = hyperReduction, modesPath = path + r'/data/modes.txt', RIDPath = path + r'/data/reducedFF_SofiaLeg_1_RID.txt', weightsPath = path + r'/data/reducedFF_SofiaLeg_1_weight.txt')
     SofiaLeg.addObject('BoxROI' , name = 'boxROICollision', orientedBox = [[-25.0, -41.0, -7.0], [25.0, -42.0, -7.0], [25.0, -39.0, -7.0], 2.0, [-25.0, -42.0, 7.0], [25.0, -42.0, 7.0], [25.0, -39.0, 7.0], 2.0], drawPoints = '0', computeEdges = '0', computeTriangles = '0', computeTetrahedra = '0', computeHexahedra = '0', computeQuad = '0', drawSize = 5, drawBoxes = True)
     SofiaLeg.addObject('BoxROI' , name = 'boxROIMiddle', orientedBox = [[-2.5, -8.5, 0.0], [2.5, -8.5, 0.0], [2.5, -3.5, 0.0], 18.0], drawBoxes = True)
     SofiaLeg.addObject('HyperReducedRestShapeSpringsForceField' , externalIndices = [0, 1, 2], indices = '@boxROIMiddle.indices', name = 'reducedFF_SofiaLeg_2', stiffness = '1e8', external_rest_shape = '@../../SofiaLeg_actuator/actuatorState', nbModes = nbrOfModes, performECSW = hyperReduction, modesPath = path + r'/data/modes.txt', RIDPath = path + r'/data/reducedFF_SofiaLeg_2_RID.txt', weightsPath = path + r'/data/reducedFF_SofiaLeg_2_weight.txt')
diff --git a/python/mor/utility/sceneCreation.py b/python/mor/utility/sceneCreation.py
@@ -34,7 +34,8 @@
                             'HyperReducedHexahedronFEMForceField':'hexahedra',
                             'HyperReducedTetrahedronFEMForceField':'tetrahedra',
                             'HyperReducedTriangleFEMForceField':'triangles',
-                            'HyperReducedRestShapeSpringsForceField':'indices'
+                            'HyperReducedRestShapeSpringsForceField':'indices',
+                            'FixedWeakConstraint':'indices'
                         }
 
 import Sofa
diff --git a/python/mor/wrapper/replaceAndSave.py b/python/mor/wrapper/replaceAndSave.py
@@ -37,7 +37,8 @@
                             'HexahedronFEMForceField':('HyperReducedHexahedronFEMForceField','hexahedra'),
                             'TetrahedronFEMForceField':('HyperReducedTetrahedronFEMForceField','tetrahedra'),
                             'TriangleFEMForceField':('HyperReducedTriangleFEMForceField','triangles'),
-                            'RestShapeSpringsForceField':('HyperReducedRestShapeSpringsForceField','indices')
+                            'RestShapeSpringsForceField':('HyperReducedRestShapeSpringsForceField','indices'),
+                            'FixedWeakConstraint':('HyperReducedFixedWeakConstraint','indices')
 
                         }
 
diff --git a/src/ModelOrderReduction/component/forcefield/HyperReducedFixedWeakConstraint.cpp b/src/ModelOrderReduction/component/forcefield/HyperReducedFixedWeakConstraint.cpp
@@ -0,0 +1,38 @@
+/******************************************************************************
+*            Model Order Reduction plugin for SOFA                            *
+*                         version 1.0                                         *
+*                       Copyright © Inria                                     *
+*                       All rights reserved                                   *
+*                       2018                                                  *
+*                                                                             *
+* This software is under the GNU General Public License v2 (GPLv2)            *
+*            https://www.gnu.org/licenses/licenses.en.html                    *
+*                                                                             *
+*                                                                             *
+*                                                                             *
+* Authors: Olivier Goury, Felix Vanneste                                      *
+*                                                                             *
+* Contact information: https://project.inria.fr/modelorderreduction/contact   *
+******************************************************************************/
+#define SOFA_COMPONENT_FORCEFIELD_HYPERREDUCEDFIXEDWEAKCONSTRAINT_CPP
+
+#include <ModelOrderReduction/component/forcefield/HyperReducedFixedWeakConstraint.inl>
+
+#include <sofa/helper/visual/DrawTool.h>
+#include <sofa/core/ObjectFactory.h>
+
+namespace sofa::component::solidmechanics::spring
+{
+
+using namespace sofa::defaulttype;
+
+
+SOFA_DECL_CLASS(HyperReducedFixedWeakConstraint)
+
+
+int HyperReducedFixedWeakConstraintClass = core::RegisterObject("Simple elastic springs applied to given degrees of freedom between their current and rest shape position")
+        .add< HyperReducedFixedWeakConstraint<Vec3Types> >()
+        ;
+
+template class SOFA_MODELORDERREDUCTION_API HyperReducedFixedWeakConstraint<Vec3Types>;
+} // namespace sofa::component::solidmechanics::spring
diff --git a/src/ModelOrderReduction/component/forcefield/HyperReducedFixedWeakConstraint.h b/src/ModelOrderReduction/component/forcefield/HyperReducedFixedWeakConstraint.h
diff --git a/src/ModelOrderReduction/component/forcefield/HyperReducedFixedWeakConstraint.inl b/src/ModelOrderReduction/component/forcefield/HyperReducedFixedWeakConstraint.inl
diff --git a/src/ModelOrderReduction/component/forcefield/HyperReducedRestShapeSpringsForceField.cpp b/src/ModelOrderReduction/component/forcefield/HyperReducedRestShapeSpringsForceField.cpp
diff --git a/src/ModelOrderReduction/component/forcefield/HyperReducedRestShapeSpringsForceField.h b/src/ModelOrderReduction/component/forcefield/HyperReducedRestShapeSpringsForceField.h

Original file line number	Diff line number	Diff line change
`@@ -34,7 +34,8 @@`
`34`	`34`	`'HyperReducedHexahedronFEMForceField':'hexahedra',`
`35`	`35`	`'HyperReducedTetrahedronFEMForceField':'tetrahedra',`
`36`	`36`	`'HyperReducedTriangleFEMForceField':'triangles',`
`37`		`- 'HyperReducedRestShapeSpringsForceField':'indices'`
	`37`	`+ 'HyperReducedRestShapeSpringsForceField':'indices',`
	`38`	`+ 'FixedWeakConstraint':'indices'`
`38`	`39`	`}`
`39`	`40`
`40`	`41`	`import Sofa`
Original file line number	Diff line number	Diff line change
`@@ -37,7 +37,8 @@`
`37`	`37`	`'HexahedronFEMForceField':('HyperReducedHexahedronFEMForceField','hexahedra'),`
`38`	`38`	`'TetrahedronFEMForceField':('HyperReducedTetrahedronFEMForceField','tetrahedra'),`
`39`	`39`	`'TriangleFEMForceField':('HyperReducedTriangleFEMForceField','triangles'),`
`40`		`- 'RestShapeSpringsForceField':('HyperReducedRestShapeSpringsForceField','indices')`
	`40`	`+ 'RestShapeSpringsForceField':('HyperReducedRestShapeSpringsForceField','indices'),`
	`41`	`+ 'FixedWeakConstraint':('HyperReducedFixedWeakConstraint','indices')`
`41`	`42`
`42`	`43`	`}`
`43`	`44`