Skip to content
Snippets Groups Projects
Commit ac625ed1 authored by Vinayak Gholap's avatar Vinayak Gholap
Browse files

Update test strong coupling SMP heating shape recovery

parent 913a8e9f
No related branches found
No related tags found
1 merge request!419Updates in strong coupling tests (change cpvac, Kvac to get right delatT in...
......@@ -73,8 +73,8 @@ wm0 = 5.
xi = 0. # rate parameter for Tt
alphag0 = 0 #0.8e-3*EAM/Eg #12.0e-5
alphar0 = 8.e-4 #15.0e-5
alphaAD0= -0.015 #-0.0#-0.024
alphar0 = 0. #8.e-4 #15.0e-5
alphaAD0= 0. #-0.015 #-0.0#-0.024
#crsytal contraction for low zg only
crystallizationRate = piecewiseScalarFunction() #surface should be 1
crystallizationRate.put(0.,20.)
......@@ -89,7 +89,7 @@ TaylorQuineyG=TaylorQuineyR=TaylorQuineyAM= 1. # fraction of plastic energy con
cfG = 1.0
cfR = 0.0
alphaAM0= 1.e-5#alphar0
alphaAM0= 0. #1.e-5#alphar0
cAM = 2020.*rho
KxAM=KyAM=KzAM = Kxr
......@@ -227,7 +227,7 @@ lawSMPPheno= PhenomenologicalSMPDG3DMaterialLaw(lawnum+7, rho, alpha, beta, gamm
TaylorQuineyG,TaylorQuineyR,TaylorQuineyAM, zAM)
lawSMPPheno.setFunctionCrystallizationVolumeRate(crystallizationRate)
lawSMPPheno.setStrainOrder(-1)
lawSMPPheno.setStrainOrder(5)
lawSMPPheno.setJ2IsotropicHardeningTractionG(_j2IHGt)
lawSMPPheno.setJ2IsotropicHardeningCompressionG(_j2IHGc)
lawSMPPheno.setKinematicHardeningG(hardenkG)
......@@ -339,7 +339,7 @@ meshfile="undeformedSMPFullDomain.msh" # name of mesh file
# solver
sol = 2 # Gmm=0 (default) Taucs=1 PETsc=2
soltype = 1 # StaticLinear=0 (default) StaticNonLinear=1, Explicit=2, Multi=3, Implicit=4, Eigen=5
nstep = 900 # 8*EMncycles # number of step (used only if soltype=1)
nstep = 600 # 8*EMncycles # number of step (used only if soltype=1)
ftime = 200. # EMncycles/freq; # according to characteristic time
tol=1.e-4 # relative tolerance for NR scheme (used only if soltype=1)
nstepArch=50 #(int)(nstep/100) # Number of step between 2 archiving (used only if soltype=1)
......@@ -461,7 +461,8 @@ mysolver.addMaterialLaw(lawvac)
mysolver.Scheme(soltype)
mysolver.Solver(sol)
mysolver.snlData(nstep,ftime,tol,tol/100.)#forsmp
mysolver.options("-ksp_type preonly -pc_type lu -pc_factor_mat_solver_type petsc")
mysolver.options("-ksp_type preonly -pc_type lu")
#mysolver.options("-pc_factor_mat_solver_type petsc")
mysolver.stepBetweenArchiving(nstepArch)
mysolver.snlManageTimeStep(50, 3, 2, 10)
......@@ -760,18 +761,21 @@ check.equal(0.0,mysolver.getArchivedNodalValue(SmpRefPoint,4,mysolver.displaceme
# ===[CONTINUE SIMULATION]=========================================================================
# ===Reset boundary conditions=====================================================================
mysolver.clearNumStepSeries() # used to clear setNumStepTimeInterval option settings
LastForce= mysolver.getArchivedForceOnPhysicalGroup("Face", SmpRef, 1)
#LastForce= mysolver.getArchivedForceOnPhysicalGroup("Face", SmpRef, 1)
mysolver.resetBoundaryConditions()
# creation of Solver
nstep2=1600;
nstepArch2= 1; #(int)(nstep2/100)
nstepArch2= 5; #(int)(nstep2/100)
#soltype2=4; #StaticLinear=0, StaticNonLinear=1, Explicit=2, Multi=3, Implicit=4, Eigen=5
soltype2=1;
sol2 = 2 # Gmm=0 (default) Taucs=1 PETsc=2
ftime2 = 1.e-3 #100. # s
mysolver.Scheme(soltype2)
mysolver.Solver(sol)
mysolver.Solver(sol2)
mysolver.snlData(nstep2,ftime2,tol,tol/100.)
mysolver.options("-ksp_type preonly -pc_type lu")
#mysolver.options("-pc_factor_mat_solver_type petsc")
#mysolver.implicitSpectralRadius(0.01)
mysolver.snlManageTimeStep(50, 3, 2, 10)
mysolver.stepBetweenArchiving(nstepArch2)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment