Add peridynamics module #11562

hchen139 · 2018-05-21T17:15:04Z

closes #11561

hchen139 · 2018-05-21T17:25:50Z

@bwspenc Here is the PR.

moosebuild · 2018-05-21T17:55:22Z

Job Documentation on d603440 wanted to post the following:

View the site here

This comment will be updated on new commits.

stale · 2018-07-17T16:58:07Z

This pull request has been automatically marked as stale because it has not had recent activity. It will be closed if no further activity occurs. Thank you for your contributions.

permcody · 2018-07-17T19:13:19Z

@hchen139 - Please fix this up. We've had discussions about this and will add it to the modules when it's updated. Thanks for your patience!

dschwen · 2018-08-20T14:24:53Z

Who will be working on this now that @hchen139 is gone?

hchen139 · 2018-08-24T14:16:50Z

@dschwen I will continue to work on getting it merged into MOOSE.

bwspenc · 2018-09-21T23:18:57Z

I'll get to this early next week, stale bot.

naveen-pr · 2018-09-24T21:17:35Z

modules/peridynamics/src/mesh/FileMeshPD.C

+    Real Y = (_pdnode[i].coord)(1);
+    Real Z = (_pdnode[i].coord)(2);
+    Real dis = std::sqrt(X * X + Y * Y + Z * Z);
+    if (dis < 0.001)


@hchen139 Instead of hard coding in a length parameter like 0.001 here, you think it would be better to make it a function of horizon (1/3 or 1/2 horizon or something like that) ?

You are right, but we need to pre calculate a quantity such as average horizon size for the irregular mesh. Sometimes this may be problematic since a global average horizon size may not necessarily presenting the local horizon size accurately for some irregular meshes.

I think we can calculate the minimum spacing between two material points for a irregular discretization, and use this value to assign boundary ids. Will do this.

naveen-pr · 2018-09-25T03:31:45Z

modules/peridynamics/include/kernels/MechanicsBaseNOSPD.h

+  virtual RankTwoTensor computeDSDU(unsigned int component, unsigned int nd);
+
+  ///@{ Material point based material properties
+  const MaterialProperty<Real> & _multi;


@hchen139 Not able to figure out what 'multi' does (?)

It is a "factor" for the force state and its value is calculated in material class. Its value is calculated as: weight function * v_i * v_j * bond_associated_horizon_volume / horizon_volume.

Ok got it, so it is a factor for the irregular discretization I'm guessing ?

This factor is not specific to irregular discretization. But the expression, weight function * v_i * v_j * bond_associated_horizon_volume / horizon_volume, is specific to bond-associated correspondence formulation. For conventional correspondence formulation, it is "weight function * v_i * v_j".

naveen-pr · 2018-09-26T20:53:18Z

modules/peridynamics/src/kernels/FiniteStrainMechanicsNOSPD.C

+
+  std::vector<RankTwoTensor> nodal_force(_nnodes);
+  for (unsigned int nd = 0; nd < _nnodes; ++nd)
+    nodal_force[nd] = (_dgrad[nd].det() * _stress[nd] * _dgrad[nd].inverse().transpose()) *


@hchen139 Referring to your other comment, where does the 'multi' here for the finite strain mechanics get calculated ?

We have it on the line 55.

naveen-pr · 2018-09-27T17:28:59Z

modules/peridynamics/src/kernels/FiniteStrainMechanicsNOSPD.C

+  for (unsigned int nd = 0; nd < _nnodes; ++nd)
+    nodal_force[nd] = (_dgrad[nd].det() * _stress[nd] * _dgrad[nd].inverse().transpose()) *
+                      _shape[nd].inverse() * _multi[nd];
+


// Residual = - (T * xi - T' * -xi )
// = - (T * xi + T' * xi )
// = - (T + T') * xi
// T, T' are nodal 'forces' (but are rank two and operate on _origin_vec_ij)

bwspenc · 2018-10-10T21:27:30Z

I'm not sure what's up with all the merge commits. Do 'git pull --rebase upstream devel', and then push -f this, and you should just have just your commits.

hchen139 · 2018-10-11T14:11:46Z

I'm not sure what's up with all the merge commits. Do 'git pull --rebase upstream devel', and then push -f this, and you should just have just your commits.

@bwspenc I thought this was what I have done for my previous commit. Anyway, redid that and this one looks right now.

tophmatthews · 2018-10-22T22:24:54Z

@bwspenc @hchen139 @naveen-pr How is this PR going? Think it can be merged sometime soon? @naveen-pr has some enhancements he's been working on he'd like to merge after it's all in.

hchen139 · 2018-10-24T13:30:24Z

@tophmatthews We don't have exact time when this will be merged. It all depends on the review process. If someone besides Ben can help to review the code, it will definitely shorten the time.

@naveen-pr For now, I think you can submit your changes to my repo and I will push your changes into the main PR.

naveen-pr · 2018-11-30T17:37:03Z

modules/peridynamics/src/materials/SmallStrainNOSPD.C

+SmallStrainNOSPD::computeQpTotalStrain()
+{
+  // the green-lagrange strain tensor
+  _total_strain[_qp] = 0.5 * (_deformation_gradient[_qp].transpose() * _deformation_gradient[_qp] -


Why use the green-lagrange strain for small strain ?

Since it is small strain, it should doesn't matter which measure you choose. And the green-lagrange strain was use here.

True, but if by some chance you use small strain NOSPD but with large deformation, won't this calculate large strains but do the rest the of mechanics with the small strain assumption ?

This is possible. But we have the FiniteStrainNOSPD for large/finite strain case.

bwspenc · 2018-12-12T15:49:01Z

@hchen139 Would you mind rebasing, resolving the comments, and force-pushing this branch to resolve the conflicts?

…daholab#11561)

…cs userobject for multi-threading (idaholab#11561)

moosebuild · 2018-12-12T21:30:55Z

Job Precheck on d603440 wanted to post the following:

Warning: This PR changes repo size by 6.04 MiB.

stale · 2019-01-12T05:02:30Z

This pull request has been automatically marked as stale because it has not had recent activity. It will be closed if no further activity occurs. Thank you for your contributions.

stale · 2019-01-13T09:51:02Z

Closing due to 30 days of inactivity. See http://mooseframework.org/moose/framework_development/patch_to_code.html

hchen139 force-pushed the peridynamics branch from ca67985 to 713eaa9 Compare May 21, 2018 17:21

hchen139 changed the title ~~Add peridynamics module (#11561)~~ Add peridynamics module May 21, 2018

hchen139 force-pushed the peridynamics branch 4 times, most recently from 08d31d8 to 58e7252 Compare May 24, 2018 18:56

hchen139 force-pushed the peridynamics branch 2 times, most recently from ae1a507 to a459f42 Compare June 6, 2018 22:18

moosebuild added the PR: Failed but allowed label Jun 6, 2018

stale bot added the wontfix label Jul 17, 2018

stale bot removed the wontfix label Jul 17, 2018

hchen139 force-pushed the peridynamics branch from a459f42 to fc1cc02 Compare July 31, 2018 14:43

moosebuild removed the PR: Failed but allowed label Jul 31, 2018

hchen139 force-pushed the peridynamics branch from fc1cc02 to 9743370 Compare August 5, 2018 05:43

hchen139 force-pushed the peridynamics branch 3 times, most recently from dacf2ec to e0d5644 Compare August 25, 2018 23:35

moosebuild added the PR: Failed but allowed label Aug 26, 2018

moosebuild removed the PR: Failed but allowed label Sep 6, 2018

naveen-pr reviewed Sep 24, 2018

View reviewed changes

naveen-pr reviewed Sep 25, 2018

View reviewed changes

naveen-pr reviewed Sep 26, 2018

View reviewed changes

naveen-pr reviewed Sep 27, 2018

View reviewed changes

hchen139 force-pushed the peridynamics branch from a26fa0e to 080afa1 Compare October 9, 2018 14:37

moosebuild added PR: Failed but allowed and removed PR: Failed but allowed labels Oct 9, 2018

moosebuild added the PR: Failed but allowed label Oct 10, 2018

hchen139 force-pushed the peridynamics branch from 0650d54 to 517d39b Compare October 11, 2018 14:01

moosebuild removed the PR: Failed but allowed label Oct 11, 2018

hchen139 force-pushed the peridynamics branch from 517d39b to 449d9a2 Compare October 11, 2018 15:05

moosebuild added the PR: Failed but allowed label Oct 11, 2018

rwcarlsen changed the base branch from devel to next November 7, 2018 19:30

moosebuild removed the PR: Failed but allowed label Nov 7, 2018

naveen-pr reviewed Nov 30, 2018

View reviewed changes

hchen139 and others added 4 commits December 12, 2018 11:15

Add peridynamics module (idaholab#11561)

96af544

First round of review changes ref idaholab#11561

d5d0343

Updated all test files to meet SQA requirement and PeridynamicsApp (i…

0a9dcd7

…daholab#11561)

Updated PeridynamicsTestApp and ModulesApp, fixed issue of peridynami…

0813b6b

…cs userobject for multi-threading (idaholab#11561)

hchen139 force-pushed the peridynamics branch from 74fb671 to 2b19bd1 Compare December 12, 2018 16:20

Fix test failure in PR (idaholab#11561)

d603440

hchen139 force-pushed the peridynamics branch from 2b19bd1 to d603440 Compare December 12, 2018 21:28

moosebuild added the PR: Failed but allowed label Dec 12, 2018

stale bot added the stale PRs that have reached or exceeded 90 days with no activity label Jan 12, 2019

stale bot closed this Jan 13, 2019

tophmatthews mentioned this pull request Jun 13, 2019

Frictionless and frictional contact with new mortar #13508

Merged

Provide feedback

Saved searches

Use saved searches to filter your results more quickly

Add peridynamics module #11562

Add peridynamics module #11562

hchen139 commented May 21, 2018

hchen139 commented May 21, 2018

moosebuild commented May 21, 2018 •

edited

Loading

stale bot commented Jul 17, 2018

permcody commented Jul 17, 2018

dschwen commented Aug 20, 2018

hchen139 commented Aug 24, 2018

bwspenc commented Sep 21, 2018

naveen-pr Sep 24, 2018

hchen139 Sep 25, 2018

hchen139 Sep 27, 2018

naveen-pr Sep 25, 2018

hchen139 Sep 25, 2018

naveen-pr Sep 26, 2018

hchen139 Sep 26, 2018

naveen-pr Sep 26, 2018

hchen139 Sep 26, 2018

naveen-pr Sep 27, 2018 •

edited

Loading

bwspenc commented Oct 10, 2018

hchen139 commented Oct 11, 2018

tophmatthews commented Oct 22, 2018

hchen139 commented Oct 24, 2018

naveen-pr Nov 30, 2018

hchen139 Nov 30, 2018

naveen-pr Nov 30, 2018

hchen139 Nov 30, 2018

bwspenc commented Dec 12, 2018

moosebuild commented Dec 12, 2018

stale bot commented Jan 12, 2019

stale bot commented Jan 13, 2019

Add peridynamics module #11562

Add peridynamics module #11562

Conversation

hchen139 commented May 21, 2018

hchen139 commented May 21, 2018

moosebuild commented May 21, 2018 • edited Loading

stale bot commented Jul 17, 2018

permcody commented Jul 17, 2018

dschwen commented Aug 20, 2018

hchen139 commented Aug 24, 2018

bwspenc commented Sep 21, 2018

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

naveen-pr Sep 27, 2018 • edited Loading

Choose a reason for hiding this comment

bwspenc commented Oct 10, 2018

hchen139 commented Oct 11, 2018

tophmatthews commented Oct 22, 2018

hchen139 commented Oct 24, 2018

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

bwspenc commented Dec 12, 2018

moosebuild commented Dec 12, 2018

stale bot commented Jan 12, 2019

stale bot commented Jan 13, 2019

moosebuild commented May 21, 2018 •

edited

Loading

naveen-pr Sep 27, 2018 •

edited

Loading