Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

add three-point-bending test #40

Merged
merged 1 commit into from
Jun 20, 2019
Merged

add three-point-bending test #40

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
214 changes: 214 additions & 0 deletions tests/threepoint_1.mpirun=2.output
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,214 @@
Running on 2 cores
Cells: 280

DoFs: 650 solid + 325 phase = 975

===================================================================
Parameters
==========
h (min): 0.231898
k: 1e-10
eps: 0.463796
G_c: 0.25
gamma penal: 0
Poisson nu: 0.3
E modulus: 1000
Lame mu: 8000
Lame lambda: 12000

Write solution 0


=======================================================================
Timestep 0: 0 (0.005) Cells: 280 DoFs: 975
-----------------------------------------------------------------------

It. #A.Set Residual Reduction LSrch #LinIts
0 1.210342e+02
1 170 1.152352e-01 9.520877e-04 0 2
2 74 1.871533e-13 1.624099e-12 0 1
3 0 1.997143e-03 1.067117e+10 10 1
4 0 2.134712e-13 1.068883e-10 0 1

No 0 time 0.005 bulk energy: 0.0110312 crack energy: 5.03491e-05 P11: 8.20409 PStress: 0.000508479
Write solution 1


=======================================================================
Timestep 1: 0.005 (0.005) Cells: 280 DoFs: 975
-----------------------------------------------------------------------

It. #A.Set Residual Reduction LSrch #LinIts
0 1.212139e+02
1 317 8.320615e-02 6.864406e-04 0 2
2 119 4.749894e-06 5.708586e-05 0 2
3 76 4.218550e-13 8.881356e-08 0 2
4 47 4.107307e-13 9.736300e-01 0 1
5 40 3.480733e-13 8.474489e-01 0 1
6 40 3.432909e-13 9.862606e-01 2 2

No 1 time 0.01 bulk energy: 0.0437771 crack energy: 0.000238815 P11: 15.8262 PStress: 0.00100935
Write solution 2


=======================================================================
Timestep 2: 0.01 (0.005) Cells: 280 DoFs: 975
-----------------------------------------------------------------------

It. #A.Set Residual Reduction LSrch #LinIts
0 1.210704e+02
1 155 5.723470e-02 4.727392e-04 0 2
2 21 7.368617e-07 1.287439e-05 0 2
3 6 5.437424e-13 7.379165e-07 0 2
4 2 4.744945e-13 8.726459e-01 0 1
5 1 3.639479e-05 7.670224e+07 10 1
6 1 5.066310e-13 1.392042e-08 0 1

No 2 time 0.015 bulk energy: 0.0968276 crack energy: 0.00114633 P11: 24.0067 PStress: 0.00149059
Write solution 3


=======================================================================
Timestep 3: 0.015 (0.005) Cells: 280 DoFs: 975
-----------------------------------------------------------------------

It. #A.Set Residual Reduction LSrch #LinIts
0 1.211129e+02
1 141 7.446003e-02 6.147987e-04 0 2
2 39 1.186125e-06 1.592968e-05 0 2
3 1 9.455601e-13 7.971843e-07 0 2
4 0 8.306558e-13 8.784801e-01 0 1
5 0 7.510008e-13 9.041060e-01 0 2

No 3 time 0.02 bulk energy: 0.167347 crack energy: 0.00382448 P11: 31.2408 PStress: 0.00193301
Write solution 4


=======================================================================
Timestep 4: 0.02 (0.005) Cells: 280 DoFs: 975
-----------------------------------------------------------------------

It. #A.Set Residual Reduction LSrch #LinIts
0 1.212247e+02
1 162 1.770591e-01 1.460585e-03 0 2
2 20 9.798712e-06 5.534149e-05 0 2
3 0 1.099841e-11 1.122434e-06 0 2
4 0 8.095483e-13 7.360594e-02 0 2

No 4 time 0.025 bulk energy: 0.250962 crack energy: 0.0101618 P11: 37.409 PStress: 0.00232195
Write solution 5


=======================================================================
Timestep 5: 0.025 (0.005) Cells: 280 DoFs: 975
-----------------------------------------------------------------------

It. #A.Set Residual Reduction LSrch #LinIts
0 1.213023e+02
1 151 3.402768e-01 2.805196e-03 0 2
2 23 4.677682e-05 1.374670e-04 0 2
3 0 8.463428e-11 1.809321e-06 0 2
4 0 1.002588e-12 1.184612e-02 0 2

No 5 time 0.03 bulk energy: 0.340422 crack energy: 0.0238955 P11: 42.5056 PStress: 0.00262118
Write solution 6


=======================================================================
Timestep 6: 0.03 (0.005) Cells: 280 DoFs: 975
-----------------------------------------------------------------------

It. #A.Set Residual Reduction LSrch #LinIts
0 1.213986e+02
1 155 5.764120e-01 4.748096e-03 0 2
2 22 1.565866e-03 2.716575e-03 0 2
3 0 3.173206e-09 2.026486e-06 0 2
4 0 1.473039e-12 4.642117e-04 0 2

No 6 time 0.035 bulk energy: 0.421071 crack energy: 0.0537582 P11: 45.9175 PStress: 0.00275701
Write solution 7


=======================================================================
Timestep 7: 0.035 (0.005) Cells: 280 DoFs: 975
-----------------------------------------------------------------------

It. #A.Set Residual Reduction LSrch #LinIts
0 1.216195e+02
1 156 8.740304e-01 7.186600e-03 0 2
2 15 1.697626e-03 1.942296e-03 0 2
3 0 2.446063e-08 1.440873e-05 0 2
4 0 1.996862e-12 8.163578e-05 0 2

DoFs: 770 solid + 385 phase = 1155
MESH CHANGED!


=======================================================================
Timestep 7: 0.035 (0.005) Cells: 328 DoFs: 1155
-----------------------------------------------------------------------

It. #A.Set Residual Reduction LSrch #LinIts
0 1.218326e+02
1 177 1.522263e+00 1.249471e-02 0 2
2 177 2.183220e-02 1.434194e-02 0 2
3 133 9.675654e-06 4.431827e-04 0 2
4 110 4.623936e-12 4.778939e-07 0 2
5 90 1.465700e-12 3.169810e-01 0 1
6 73 1.379436e-12 9.411452e-01 0 1
7 58 2.860842e-05 2.073921e+07 10 1
8 58 1.486521e-12 5.196097e-08 0 1

No 7 time 0.04 bulk energy: 0.445505 crack energy: 0.121808 P11: 45.7943 PStress: 0.00243689
Write solution 8


=======================================================================
Timestep 8: 0.04 (0.005) Cells: 328 DoFs: 1155
-----------------------------------------------------------------------

It. #A.Set Residual Reduction LSrch #LinIts
0 1.217683e+02
1 200 1.304519e+01 1.071312e-01 0 2
2 281 1.856793e+00 1.423355e-01 0 2
3 235 1.213703e+00 6.536556e-01 0 2
4 207 1.654739e-01 1.363380e-01 0 2
5 190 9.703670e-04 5.864169e-03 0 2
6 180 1.736472e-07 1.789500e-04 0 2
7 178 1.454823e-12 8.378043e-06 0 2
8 178 1.442770e-12 9.917147e-01 0 2

DoFs: 898 solid + 449 phase = 1347
MESH CHANGED!


=======================================================================
Timestep 8: 0.04 (0.005) Cells: 388 DoFs: 1347
-----------------------------------------------------------------------

It. #A.Set Residual Reduction LSrch #LinIts
0 1.217445e+02
1 228 1.287882e+01 1.057857e-01 0 2
2 285 1.867728e+00 1.450232e-01 0 2
3 233 1.208468e+00 6.470255e-01 0 2
4 202 1.726925e-01 1.429021e-01 0 2
5 187 1.072328e-03 6.209466e-03 0 2
6 178 2.311829e-07 2.155896e-04 0 2
7 178 1.666521e-12 7.208669e-06 0 2

No 8 time 0.045 bulk energy: 0.31501 crack energy: 0.253812 P11: 37.5968 PStress: 0.00088921
Write solution 9

Finishing time step loop: 0.35967


+---------------------------------------------+------------+------------+
+---------------------------------+-----------+------------+------------+
+---------------------------------+-----------+------------+------------+



+---------------------------------------------+------------+------------+
+---------------------------------+-----------+------------+------------+
+---------------------------------+-----------+------------+------------+

52 changes: 52 additions & 0 deletions tests/threepoint_1.prm
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,52 @@
# Three point bending test from
# Miehe, Hofacker, Welschinger, 2010
# --------------------------------------------------------

subsection Global parameters
set Global pre-refinement steps = 0
set Local pre-refinement steps = 0
set Adaptive refinement cycles = 1

set Max No of timesteps = 8
set Timestep size = 5.0e-3
set Timestep size to switch to = 1.0e-3
set Switch timestep after steps = 35


set outer solver = active set
set test case = three point bending
# phase field,
# phase field three point top
# mix
set ref strategy = phase field

set value phase field for refinement = 0.5
set Output filename = solution_tp_
end

subsection Problem dependent parameters
set K reg = 1.0e-10
set Eps reg = 2.0*h
set Gamma penalization = 0.0

# Material and problem-rhs parameters
set Pressure = 0.0
set Fracture toughness G_c = 0.25
set Density solid = 1.0
set Poisson ratio nu = 0.3
set E modulus = 1e+3
set Lame mu = 8.0e+3
set Lame lambda = 12.0e+3
end

subsection Solver parameters
set Use Direct Inner Solver = false
set Newton lower bound = 1.0e-6
set Newton maximum steps = 30
set Upper Newton rho = 0.99
set Line search maximum steps = 10
set Line search damping = 0.6
set Decompose stress in rhs = 1.0
set Decompose stress in matrix = 1.0

end