Skip to content
Snippets Groups Projects
Commit e2bd920b authored by Ludovic Noels's avatar Ludovic Noels
Browse files
parents a6c614da b6325467
No related branches found
No related tags found
No related merge requests found
SetFactory('OpenCASCADE'); SetFactory('OpenCASCADE');
Mesh.SecondOrderLinear = 1; Mesh.SecondOrderLinear = 1;
Mesh.CharacteristicLengthMin =0.25; Mesh.CharacteristicLengthMin =0.525;
Mesh.CharacteristicLengthMax =0.5; Mesh.CharacteristicLengthMax =0.5;
Box(1) = {0,0,0,10.0,10.0,5.0}; Box(1) = {0,0,0,10.0,10.0,5.0};
Cylinder(2) = {0,0,0,0,10.0,0,0.5,2*Pi}; Cylinder(2) = {0,0,0,0,10.0,0,0.5,2*Pi};
......
This diff is collapsed.
...@@ -153,10 +153,10 @@ else: ...@@ -153,10 +153,10 @@ else:
# solver # solver
sol = 2 # Gmm=0 (default) Taucs=1 PETsc=2 sol = 2 # Gmm=0 (default) Taucs=1 PETsc=2
soltype =1 # StaticLinear=0 (default) StaticNonLinear=1 soltype =1 # StaticLinear=0 (default) StaticNonLinear=1
nstep = 2000 # number of step (used only if soltype=1) nstep = 10000 # number of step (used only if soltype=1)
ftime=0.5 ftime=0.5
tol=1.e-4 # relative tolerance for NR scheme (used only if soltype=1) tol=1.e-4 # relative tolerance for NR scheme (used only if soltype=1)
nstepArch=20 # Number of step between 2 archiving (used only if soltype=1) nstepArch=100 # Number of step between 2 archiving (used only if soltype=1)
# creation of Solver # creation of Solver
mysolver = nonLinearMechSolver(1000) mysolver = nonLinearMechSolver(1000)
...@@ -172,7 +172,7 @@ mysolver.snlData(nstep,ftime,tol,tol/100.) ...@@ -172,7 +172,7 @@ mysolver.snlData(nstep,ftime,tol,tol/100.)
mysolver.stepBetweenArchiving(nstepArch) mysolver.stepBetweenArchiving(nstepArch)
mysolver.options("-ksp_type preonly -pc_type lu -pc_factor_mat_solver_type petsc") mysolver.options("-ksp_type preonly -pc_type lu -pc_factor_mat_solver_type petsc")
mysolver.snlManageTimeStep(150,5,2.,500) mysolver.snlManageTimeStep(150,2,2.,500)
stiffModif = BFGSStiffnessModificationMonitoring(20, False) # stiffness will be re-computed after 20 iterations in each step stiffModif = BFGSStiffnessModificationMonitoring(20, False) # stiffness will be re-computed after 20 iterations in each step
mysolver.stiffnessModification(stiffModif) mysolver.stiffnessModification(stiffModif)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment