Changeset 19394

Timestamp:

06/10/15 00:12:52 (10 years ago)

Author:

cborstad

Message:

CHG: check stabilization when setting up constraints

File:

• issm/trunk-jpl/src/c/analyses/DamageEvolutionAnalysis.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/analyses/DamageEvolutionAnalysis.cpp

@@ -11 +11 @@
         iomodel->Constant(&finiteelement,DamageElementinterpEnum);
 
-        IoModelToConstraintsx(constraints,iomodel,DamageSpcdamageEnum,DamageEvolutionAnalysisEnum,finiteelement);
-
+        /*Fetch parameters: */
+        int stabilization;
+        iomodel->Constant(&stabilization,DamageStabilizationEnum);
+
+        /*Do not add constraints in DG,  they are weakly imposed*/
+        if(stabilization!=3){
+                IoModelToConstraintsx(constraints,iomodel,DamageSpcdamageEnum,DamageEvolutionAnalysisEnum,finiteelement);
+        }
+
+        /*FCT, constraints are imposed using penalties*/
+        if(stabilization==4){
+                constraints->ActivatePenaltyMethod(DamageEvolutionAnalysisEnum);
+        }
 }/*}}}*/
 void DamageEvolutionAnalysis::CreateLoads(Loads* loads, IoModel* iomodel){/*{{{*/