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

Damage Plasticity Model #394

Merged
merged 7 commits into from
Aug 1, 2024
Merged

Damage Plasticity Model #394

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
7 commits
Select commit Hold shift + click to select a range
3eecc9b
Major Changes : Constants (ft0,fc0,at), Functions(f,fbar,df,dfbar), WIP
vprithiv Oct 16, 2023
1013d4a
changes as per f and fbar, and derivatives (beta)
vprithiv Nov 14, 2023
e508432
Refined declarations and definitions
vprithiv Jan 8, 2024
752d3a3
Changes in at calc, fbar derivative, gamma calc and derivative, dpote…
vprithiv May 16, 2024
6e8b1a2
Final cleanup and gold files update ref #184
vprithiv Jul 3, 2024
fadd3b6
Addressed review comments - Minor changes ref #184
vprithiv Jul 30, 2024
6e45724
Minor updates to address remaining review comments ref #184
bwspenc Jul 31, 2024
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
125 changes: 61 additions & 64 deletions doc/content/source/materials/DamagePlasticityStressUpdate.md
View file
Open in desktop

Large diffs are not rendered by default.

41 changes: 19 additions & 22 deletions include/materials/DamagePlasticityStressUpdate.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -53,6 +53,10 @@ class DamagePlasticityStressUpdate : public MultiParameterPlasticityStressUpdate
const Real _fc;
///fracture energy in compression (user parameter)
const Real _FEc;
/// Maximum stress in tension without damage
const Real _ft0;
/// Maximum stress in compression without damage
const Real _fc0;

///@{
/** The following variables are intermediate and are calculated based on the user parameters given
Expand All @@ -61,14 +65,9 @@ class DamagePlasticityStressUpdate : public MultiParameterPlasticityStressUpdate
const Real _ac;
const Real _zt;
const Real _zc;
const Real _dPhit;
const Real _dPhic;
const Real _sqrtPhit_max;
const Real _sqrtPhic_max;
const Real _dt_bt;
const Real _dc_bc;
const Real _ft0;
const Real _fc0;

///@}

/// Intermediate variable calculated using user parameter tip_smoother
Expand Down Expand Up @@ -111,33 +110,32 @@ class DamagePlasticityStressUpdate : public MultiParameterPlasticityStressUpdate
///damaged maximum principal stress
MaterialProperty<Real> & _sigma2;
/**
* Obtain the tensile strength
* Obtain the undamaged strength
* @param intnl (Array containing damage states in tension and compression, respectively)
* @return value of ft (tensile strength)
*/
Real ft(const std::vector<Real> & intnl) const;

Real fbar(const Real & f0, const Real & a, const Real & exponent, const Real & kappa) const;
// Real ftbar(const std::vector<Real> & intnl) const;
/**
* Obtain the partial derivative of the tensile strength to the damage state
* Obtain the partial derivative of the undamaged tensile strength to the damage state
* @param intnl (Array containing damage states in tension and compression, respectively)
* @return value of dft (partial derivative of the tensile strength to the damage state)
*/
Real dft(const std::vector<Real> & intnl) const;
Real
dfbar_dkappa(const Real & f0, const Real & a, const Real & exponent, const Real & kappa) const;

/**
* Obtain the conpressive strength
* @param intnl (Array containing damage states in tension and compression, respectively)
* @return value of fc (conpressive strength)
*/
Real fc(const std::vector<Real> & intnl) const;
// * Obtain the damaged tensile strength
// * @param intnl (Array containing damage states in tension and compression, respectively)
// * @return value of ft (tensile strength)
// */
Real f(const Real & f0, const Real & a, const Real & kappa) const;

/**
* Obtain the partial derivative of the compressive strength to the damage state
* Obtain the partial derivative of the undamaged strength to the damage state
* @param intnl (Array containing damage states in tension and compression, respectively)
* @return value of dfc
* @return value of dft (partial derivative of the tensile strength to the damage state)
*/
Real dfc(const std::vector<Real> & intnl) const;

Real df_dkappa(const Real & f0, const Real & a, const Real & kappa) const;
/**
* beta is a dimensionless constant, which is a component of the yield function
* It is defined in terms of tensile strength, compressive strength, and another
Expand Down Expand Up @@ -233,7 +231,6 @@ class DamagePlasticityStressUpdate : public MultiParameterPlasticityStressUpdate
* @return dr_dstress (dflowpotential_dstress)
*/
virtual void dflowPotential_dstress(const std::vector<Real> & stress_params,
const std::vector<Real> & intnl,
std::vector<std::vector<Real>> & dr_dstress) const;
/**
* This function calculates the derivative of the flow potential with the damage states
Expand Down
Loading