Skip to content
Snippets Groups Projects
Commit 6f6384b9 authored by Ludovic Noels's avatar Ludovic Noels
Browse files

reference value

parent 3c741e70
Branches
Tags
No related merge requests found
......@@ -80,7 +80,7 @@ if Solving_mode == 1:
ftime =1.0 # Final time
tol=1.e-5 # Relative tolerance for NR scheme
tolAbs = 1.e-6 # Absolute tolerance for NR scheme
nstepArch= 1 # Number of step between 2 archiving
nstepArch= 10 # Number of step between 2 archiving
nstepArchEnergy = 1 # Number of step between 2 energy computation and archiving
nstepArchForce = 1 # Number of step between 2 force archiving
MaxIter = 200 # Maximum number of iterations
......@@ -95,7 +95,7 @@ elif Solving_mode == 2:
ftime =1.0e-2 # Final time
tol=1.e-4 # Relative tolerance for NR scheme
tolAbs = 1.e-6 # Absolute tolerance for NR scheme
nstepArch= 1 # Number of step between 2 archiving
nstepArch= 10 # Number of step between 2 archiving
nstepArchEnergy = 1 # Number of step between 2 energy computation and archiving
nstepArchForce = 1 # Number of step between 2 force archiving
MaxIter = 200 # Maximum number of iterations
......@@ -232,10 +232,11 @@ if LoadingMode == 0:
mysolver.setPathFollowingIncrementAdaptation(True,4,0.3)
mysolver.setPathFollowingControlType(0)
mysolver.setPathFollowingCorrectionMethod(0)
mysolver.setPathFollowingArcLengthStep(1e-1)
mysolver.setBoundsOfPathFollowingArcLengthSteps(0,2.);
mysolver.setPathFollowingArcLengthStep(3e-2)
mysolver.setBoundsOfPathFollowingArcLengthSteps(0,3e-2);
mysolver.archivingForceOnPhysicalGroup("Face", 400, 1, nstepArchForce)
mysolver.archivingForceOnPhysicalGroup("Face", 200, 1, nstepArchForce)
# mysolver.pathFollowing(True,1)
# mysolver.setPathFollowingIncrementAdaptation(False,4)
# mysolver.setPathFollowingLocalSteps(9.0e-1,1.0e-3)
......@@ -321,7 +322,7 @@ mysolver.archivingNodeDisplacement(111,1, nstepArchForce)
mysolver.solve()
check = TestCheck()
check.equal(-6.645770e+05,mysolver.getArchivedForceOnPhysicalGroup("face", 400, 1),1.e-4)
check.equal(-2966.175372583341,mysolver.getArchivedForceOnPhysicalGroup("Face", 400, 1),1.e-4)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment