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

change parameters

parent c9ab9c20
No related branches found
No related tags found
No related merge requests found
...@@ -33,7 +33,7 @@ sol = 2 # Gmm=0 (default) Taucs=1 PETsc=2 ...@@ -33,7 +33,7 @@ sol = 2 # Gmm=0 (default) Taucs=1 PETsc=2
soltype = 4 #3 # StaticLinear=0 (default) StaticNonLinear=1 soltype = 4 #3 # StaticLinear=0 (default) StaticNonLinear=1
nstep = 6000 # number of step (used only if soltype=1) nstep = 6000 # number of step (used only if soltype=1)
ftime =3.e-2 # Final time (used only if soltype=1) ftime =3.e-2 # Final time (used only if soltype=1)
tol=1.e-6 # relative tolerance for NR scheme (used only if soltype=1) tol=1.e-5 # relative tolerance for NR scheme (used only if soltype=1)
nstepArch=20 # Number of step between 2 archiving (used only if soltype=1) nstepArch=20 # Number of step between 2 archiving (used only if soltype=1)
fullDg =1 #O = CG, 1 = DG fullDg =1 #O = CG, 1 = DG
space1 = 0 # function space (Lagrange=0) space1 = 0 # function space (Lagrange=0)
...@@ -73,7 +73,7 @@ mysolver.addMaterialLaw(lawcoh1) ...@@ -73,7 +73,7 @@ mysolver.addMaterialLaw(lawcoh1)
mysolver.addMaterialLaw(lawfrac1) mysolver.addMaterialLaw(lawfrac1)
mysolver.Scheme(soltype) mysolver.Scheme(soltype)
mysolver.Solver(sol) mysolver.Solver(sol)
mysolver.snlData(nstep,ftime,tol) mysolver.snlData(nstep,ftime,tol,tol/10.)
mysolver.snlManageTimeStep(30,5,20,10) mysolver.snlManageTimeStep(30,5,20,10)
#mysolver.lineSearch(bool(1)) #mysolver.lineSearch(bool(1))
mysolver.stepBetweenArchiving(nstepArch) mysolver.stepBetweenArchiving(nstepArch)
...@@ -128,8 +128,8 @@ mysolver.archivingNodeDisplacement(57,1,10) ...@@ -128,8 +128,8 @@ mysolver.archivingNodeDisplacement(57,1,10)
mysolver.solve() mysolver.solve()
check = TestCheck() check = TestCheck()
check.equal(-8.431679e+01,mysolver.getArchivedForceOnPhysicalGroup("Face", 55, 1),5.e-2) check.equal(-4.139801e+02,mysolver.getArchivedForceOnPhysicalGroup("Face", 55, 1),5.e-2)
check.equal(9.2e-04,mysolver.getArchivedNodalValue(57,1,mysolver.displacement),1e-6) check.equal(8.999786e-04,mysolver.getArchivedNodalValue(57,1,mysolver.displacement),1e-6)
try: try:
import linecache import linecache
...@@ -139,5 +139,5 @@ except: ...@@ -139,5 +139,5 @@ except:
import os import os
os._exit(1) os._exit(1)
else: else:
check.equal(8.659006e-01,float(linesDam.split(';')[1])) check.equal(9.179512e-01,float(linesDam.split(';')[1]))
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment