Skip to content

refactor: Phase-field and damage solid model cleanup, refactor #4096

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

Open
frankfeifan wants to merge 12 commits into
base: develop
Choose a base branch
from
+1,673 −1,187
Open

refactor: Phase-field and damage solid model cleanup, refactor #4096

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
12 commits
Select commit Hold shift + click to select a range
0cb5390
PhaseFieldFractureSolver cleanup
frankfeifan Jul 16, 2026
baf40f8
PhaseFieldDamageFEM cleanup, fix typo, condense, rename
frankfeifan Jul 16, 2026
24f37b0
added a check for nucleation model to see if the dissipaiton option i…
frankfeifan Jul 16, 2026
9b33c2c
move localDissipationoption to constitutive model and added FractureM…
frankfeifan Jul 18, 2026
f2818e7
updated xmls
frankfeifan Jul 18, 2026
7dad7b3
renamed sed to crackDrivingForce
frankfeifan Jul 19, 2026
0188a02
cleaned up PhaseFieldPoromechanicsSolver
frankfeifan Jul 19, 2026
f492395
apply Damage default field values after array allocation
frankfeifan Jul 19, 2026
50c014c
refactored and added spectral/vol-dev split test decks for all fractu…
frankfeifan Jul 22, 2026
488e537
fix save history of crack driving force
frankfeifan Jul 29, 2026
801b32b
enabled multi-region run with certain defined regions have no damage dof
frankfeifan Aug 4, 2026
41ea53a
Add phase-field benchmark
frankfeifan Aug 6, 2026
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
28 changes: 26 additions & 2 deletions inputFiles/phaseField/PhaseFieldFracture.ats
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -28,21 +28,45 @@ decks = [
check_step=10,
restartcheck_params=RestartcheckParameters(**restartcheck_params)),
TestDeck(
name="PhaseFieldFracture_SpectralSplit",
name="PhaseFieldFracture_SpectralSplit_Cohesive",
description=
'Testing the spectral split of the stress tensor in PF Fracture',
partitions=((1, 1, 1), (2, 2, 1)),
restart_step=15,
check_step=30,
restartcheck_params=RestartcheckParameters(**restartcheck_params)),
TestDeck(
name="PhaseFieldFracture_VolDevSplit",
name="PhaseFieldFracture_SpectralSplit_Brittle",
description=
'Testing the spectral split with the brittle degradation function in PF Fracture',
partitions=((1, 1, 1), (2, 2, 1)),
restart_step=15,
check_step=30,
restartcheck_params=RestartcheckParameters(**restartcheck_params)),
TestDeck(
name="PhaseFieldFracture_VolDevSplit_Brittle",
description=
'Testing the volumetric-deviatoric split of the stress tensor in PF Fracture',
partitions=((1, 1, 1), (2, 2, 1)),
restart_step=15,
check_step=30,
restartcheck_params=RestartcheckParameters(**restartcheck_params)),
TestDeck(
name="PhaseFieldFracture_VolDevSplit_Cohesive",
description=
'Testing the volumetric-deviatoric split with the cohesive degradation function in PF Fracture',
partitions=((1, 1, 1), (2, 2, 1)),
restart_step=15,
check_step=30,
restartcheck_params=RestartcheckParameters(**restartcheck_params)),
TestDeck(
name="PhaseFieldFracture_MultiRegion",
description=
'Testing PF Fracture on one region of a multi-region problem, with mechanics only elsewhere',
partitions=((1, 1, 1), (2, 2, 1)),
restart_step=10,
check_step=20,
restartcheck_params=RestartcheckParameters(**restartcheck_params)),
TestDeck(name="PhaseFieldFracture_Nucleation_smoke",
description='Testing the nucleation model for PF Fracture',
partitions=((1, 1, 1), (2, 2, 1)),
Expand Down
9 changes: 0 additions & 9 deletions inputFiles/phaseField/PhaseFieldFracture_CohesiveModel.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -32,7 +32,6 @@
discretization="FE1"
timeIntegrationOption="SteadyState"
fieldName="Damage"
localDissipation="Linear"
logLevel="1"
targetRegions="{ Region1 }">
<NonlinearSolverParameters
Expand Down Expand Up @@ -119,14 +118,6 @@
objectPath="nodeManager"
scale="0"
setNames="{xneg, xpos, yneg, ypos, zneg, zpos}"/> -->
<FieldSpecification
name="Gc"
initialCondition="1"
setNames="{ all }"
objectPath="ElementRegions/Region1"
fieldName="shale_criticalFractureEnergy"
scale="2.7"/>

<FieldSpecification
name="xpos"
objectPath="nodeManager"
Expand Down
13 changes: 2 additions & 11 deletions inputFiles/phaseField/PhaseFieldFracture_DamageAndLoad.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -35,8 +35,6 @@
discretization="FE1"
timeIntegrationOption="SteadyState"
fieldName="Damage"

localDissipation="Quadratic"
logLevel="1"
targetRegions="{ Region1 }">
<NonlinearSolverParameters
Expand Down Expand Up @@ -114,18 +112,11 @@
defaultShearModulus="3.00e5"
lengthScale="0.2"
defaultCriticalFractureEnergy="2.7"
criticalStrainEnergy="0"/>
criticalStrainEnergy="0"
localDissipationOption="Quadratic"/>
</Constitutive>

<FieldSpecifications>
<FieldSpecification
name="Gc"
initialCondition="1"
setNames="{ all }"
objectPath="ElementRegions/Region1"
fieldName="shale_criticalFractureEnergy"
scale="2.7"/>

<FieldSpecification
name="xtraction"
fieldName="totalDisplacement"
Expand Down
10 changes: 0 additions & 10 deletions inputFiles/phaseField/PhaseFieldFracture_DamageOnly.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -35,8 +35,6 @@
discretization="FE1"
timeIntegrationOption="SteadyState"
fieldName="Damage"

localDissipation="Linear"
logLevel="1"
targetRegions="{ Region1 }">
<NonlinearSolverParameters
Expand Down Expand Up @@ -120,14 +118,6 @@
</Constitutive>

<FieldSpecifications>
<FieldSpecification
name="Gc"
initialCondition="1"
setNames="{ all }"
objectPath="ElementRegions/Region1"
fieldName="shale_criticalFractureEnergy"
scale="2.7"/>

<FieldSpecification
name="xboundary"
fieldName="totalDisplacement"
Expand Down
76 changes: 3 additions & 73 deletions ...eField/PhaseFieldFracture_VolDevSplit.xml → ...FieldFracture_DecompositionSplit_base.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -7,11 +7,11 @@
name="PhaseFieldSolve"
solidSolverName="SolidMechSolve"
damageSolverName="DamageSolve"
logLevel="2"
logLevel="1"
targetRegions="{ Region1 }">
<NonlinearSolverParameters
couplingType="Sequential"
lineSearchAction="None"
lineSearchAction="None"
newtonMaxIter="100"/>
</PhaseFieldFracture>

Expand All @@ -33,8 +33,6 @@
discretization="FE1"
timeIntegrationOption="SteadyState"
fieldName="Damage"

localDissipation="Quadratic"
logLevel="1"
targetRegions="{ Region1 }">
<NonlinearSolverParameters
Expand Down Expand Up @@ -73,37 +71,6 @@
xMax="{ 1.01, 1.01, 0.11 }"/>
</Geometry>

<Events
maxTime="30.0">
<!-- <SoloEvent name="siloOutput0"
target="/Outputs/SiloOutput"/> -->
<!-- This event is applied every cycle, and overrides the
solver time-step request -->
<PeriodicEvent
name="solverApplications"
forceDt="1.0"
target="/Solvers/PhaseFieldSolve"/>

<!-- This event is applied every 1.0s. The targetExactTimestep
flag allows this event to request a dt modification to match an
integer multiple of the timeFrequency. -->
<PeriodicEvent
name="vtkOutput"
timeFrequency="1.0"
target="/Outputs/Mode2VolDev_VTKOutput"/>

<PeriodicEvent
name="siloOutputs"
timeFrequency="1.0"
target="/Outputs/Mode2VolDev_SiloOutput"/>

<PeriodicEvent
name="restarts"
timeFrequency="15.0"
targetExactTimestep="0"
target="/Outputs/restartOutput"/>
</Events>

<NumericalMethods>
<FiniteElements>
<FiniteElementSpace
Expand All @@ -119,26 +86,7 @@
materialList="{ shale }"/>
</ElementRegions>

<Constitutive>
<DamageVolDevElasticIsotropic
name="shale"
defaultDensity="2700"
defaultBulkModulus="1.750e5"
defaultShearModulus="8.077e4"
lengthScale="0.25"
defaultCriticalFractureEnergy="54.0"
criticalStrainEnergy="14.88"/>
</Constitutive>

<FieldSpecifications>
<FieldSpecification
name="Gc"
initialCondition="1"
setNames="{ all }"
objectPath="ElementRegions/Region1"
fieldName="shale_criticalFractureEnergy"
scale="54.0"/>

<FieldSpecification
name="shearLoad"
fieldName="totalDisplacement"
Expand Down Expand Up @@ -178,22 +126,4 @@
values="{ -1, 100 }"/>
</Functions>

<Outputs>
<Silo
name="Mode2VolDev_SiloOutput"/>

<VTK
name="Mode2VolDev_VTKOutput"
plotFileRoot="plot"/>

<Restart
name="restartOutput"/>
</Outputs>

<!-- <Geometry> -->
<!-- <Box name="bottom" xMin="{-0.51, -0.51, -0.01" xMax="0.51, -0.49, +0.11"/> -->
<!-- <Box name="top" xMin="{-0.51, 0.49, -0.01" xMax="0.51, 0.51, +0.11"/> -->
<!-- <Box name="left" xMin="{-0.51, -0.51, -0.01" xMax="-0.49, 0.51, +0.11"/> -->
<!-- <Box name="right" xMin="0.49, -0.51, -0.01" xMax="0.51, 0.51, +0.11"/> -->
<!-- </Geometry> -->
</Problem>
</Problem>
Loading