# Von Mises vs PEEQ

**URL:** <https://calculix.discourse.group/t/von-mises-vs-peeq/1855>\
**Category:** Uncategorized\
**Created:** [October 4, 2023, 3:04am UTC](https://calculix.discourse.group/t/von-mises-vs-peeq/1855 "2023-10-04T03:04:57Z")\
**Posts on this page:** 11\
**Page:** 2

<div class="post-metadata">

**Author:** ![Calc\_em](https://avatars.discourse-cdn.com/v4/letter/c/43a26b/32.png) [@Calc\_em](https://calculix.discourse.group/u/Calc_em)\
**Post date:** [October 5, 2023, 5:52pm UTC](https://calculix.discourse.group/t/von-mises-vs-peeq/1855/24 "2023-10-05T17:52:25Z")

</div>

You have to set “Anyone with the link” in the access settings for that file or we won’t be able to download it.

---

<div class="post-metadata">

**Author:** ![Juseche](https://avatars.discourse-cdn.com/v4/letter/j/258eb7/32.png) [@Juseche](https://calculix.discourse.group/u/Juseche)\
**Post date:** [October 5, 2023, 9:09pm UTC](https://calculix.discourse.group/t/von-mises-vs-peeq/1855/25 "2023-10-05T21:09:09Z")

</div>

Here,

> **[pipe\_dent.inp](https://wetransfer.com/downloads/4f7e2553306f2a41fba8f570764b4dd220231005210816/885e94)**
>
> 1 file sent via WeTransfer, the simplest way to send your files around the world

---

<div class="post-metadata">

**Author:** ![Juseche](https://avatars.discourse-cdn.com/v4/letter/j/258eb7/32.png) [@Juseche](https://calculix.discourse.group/u/Juseche)\
**Post date:** [October 5, 2023, 9:12pm UTC](https://calculix.discourse.group/t/von-mises-vs-peeq/1855/26 "2023-10-05T21:12:24Z")

</div>

Here, please use this link (“Anywone with the link”):

> **[pipe\_dent.inp](https://drive.google.com/file/d/1H5VWG3qsJgxHL-CqpPWWebrssWRYeLrX/view?usp=sharing)**
>
> Google Drive file.

---

<div class="post-metadata">

**Author:** ![Disla](https://yyz2.discourse-cdn.com/free1/user_avatar/calculix.discourse.group/disla/32/3548_2.png) [@Disla](https://calculix.discourse.group/u/Disla)\
**Post date:** [October 6, 2023, 7:06am UTC](https://calculix.discourse.group/t/von-mises-vs-peeq/1855/27 "2023-10-06T07:06:24Z")

</div>

HI,

I have torture three bricks of API-579-2021. One under Torsion, one under bending and other under a combo of torsion and bending. All three with the same mesh density and C3D20R.

 ![screenshot.16](https://global.discourse-cdn.com/free1/uploads/calculix/original/2X/7/7cf54b8301140259f74037243d3e2bb5697f1866.jpeg)  
 ![imagen](https://global.discourse-cdn.com/free1/uploads/calculix/original/2X/7/7f316c5ac2e75e1596c00fcf05329922c5ec4529.png)

I have later represented all the final VM Stress State versus PEEQ Value for the whole set of nodes.  
I would expect some discrepancy with respect to the ideal Axial, Biaxial and Planar Stress States (which has shown to behave as expected) due to interpolation but the results are worse than I expected.  
I’m jet not sure how to read it as some of the PEEQ values are negative and I can’t make them to approach to the curve with refinement.

These are the results which in some way seems to confirm Juseche discrepancy when plasticity and Bending/Torsion is present.

 ![screenshot.17](https://global.discourse-cdn.com/free1/uploads/calculix/original/2X/d/de12883b8c1f2222c68891b2803848cf23ada2fc.jpeg)  
 ![screenshot.18](https://global.discourse-cdn.com/free1/uploads/calculix/original/2X/c/ca522f5e50afa8238c001cb81852fd8bb9a08abc.jpeg)

---

<div class="post-metadata">

**Author:** ![Juseche](https://avatars.discourse-cdn.com/v4/letter/j/258eb7/32.png) [@Juseche](https://calculix.discourse.group/u/Juseche)\
**Post date:** [October 6, 2023, 11:26am UTC](https://calculix.discourse.group/t/von-mises-vs-peeq/1855/28 "2023-10-06T11:26:36Z")

</div>

I think the problem comes from elements with boundary conditions.

Try running my model and plot only points at the center of the dent.

Jack

---

<div class="post-metadata">

**Author:** ![Disla](https://yyz2.discourse-cdn.com/free1/user_avatar/calculix.discourse.group/disla/32/3548_2.png) [@Disla](https://calculix.discourse.group/u/Disla)\
**Post date:** [October 6, 2023, 12:46pm UTC](https://calculix.discourse.group/t/von-mises-vs-peeq/1855/29 "2023-10-06T12:46:02Z")

</div>

Hi, Juseche,

Nice finding.!!!. It confirms my believe that it is always worth to go all the way and clear up all doubts.

I have selected nodes away from BC and It cleans many of the Stress Strain States below the curve but I still have some questions.

 ![screenshot.23](https://global.discourse-cdn.com/free1/uploads/calculix/original/2X/7/7c673b371dc23379dbbdd0ec0468bf1288ff195a.jpeg)

-¿How can PEEQ be negative? It is the most surprising to me.

-¿How can VM end up above the curve?. Error due to extrapolating values from the integration points or ¿ Maybe some stress components do not induce plastic Strain?.

-¿How can VM end up below the curve?. Error due to extrapolating values from the integration points or ¿VM evolution at some point has start to decrease due to the load path?

 ![imagen](https://global.discourse-cdn.com/free1/uploads/calculix/original/2X/1/1dd808e320068e738fe271fff278fee7b90eb60e.png)

---

<div class="post-metadata">

**Author:** ![Juseche](https://avatars.discourse-cdn.com/v4/letter/j/258eb7/32.png) [@Juseche](https://calculix.discourse.group/u/Juseche)\
**Post date:** [October 6, 2023, 3:13pm UTC](https://calculix.discourse.group/t/von-mises-vs-peeq/1855/30 "2023-10-06T15:13:31Z")

</div>

I don’t know why PEEQ is negative and VM ends up above/below the curve… it’s weird!

However, I found that my results went better with C3D8R for SVM vs PEEQ. Please use this element in your model and check it. I think, results will go better (far from BC).

Jack

---

<div class="post-metadata">

**Author:** ![Disla](https://yyz2.discourse-cdn.com/free1/user_avatar/calculix.discourse.group/disla/32/3548_2.png) [@Disla](https://calculix.discourse.group/u/Disla)\
**Post date:** [October 6, 2023, 7:11pm UTC](https://calculix.discourse.group/t/von-mises-vs-peeq/1855/31 "2023-10-06T19:11:55Z")

</div>

I have extract PEEQ and VM (computed from the stress components) at the integration points to remove any uncertainty from the interpolation process.  
Result is much clean. 😀  
I was not aware how much noise the interpolation process was introducing into the results.  
I’m postprocesing PEEQ and VM from the integration points in the future.  
There is still an small amount of points below the curve which I need to examine carefully (Not only the BC elements). ¿Any idea why some points end up there?

> [@Juseche](#):
>
> my results went better

The fact that data points end up on top of the curve doesn’t mean the result is correct.  
I would add at leats 3-4 layers through thickness to your model. Only one is not enought.

 ![imagen](https://global.discourse-cdn.com/free1/uploads/calculix/original/2X/3/3015d2b06bec99099729217a1c309320b9732888.png)  
 ![imagen](https://global.discourse-cdn.com/free1/uploads/calculix/original/2X/9/995b76b8d93152d07f6c93f464054f27d69e1f2d.png)

---

<div class="post-metadata">

**Author:** ![Juseche](https://avatars.discourse-cdn.com/v4/letter/j/258eb7/32.png) [@Juseche](https://calculix.discourse.group/u/Juseche)\
**Post date:** [October 7, 2023, 4:48pm UTC](https://calculix.discourse.group/t/von-mises-vs-peeq/1855/32 "2023-10-07T16:48:08Z")

</div>

Very interesting!

I agree with use more elements along thickness. But elements turns bad shaped (small relation thk/long). Could this affects the results?

By the way, how to plotting not- interpolated-to-nodes element results in Calculix graphics?

Thanks!

---

<div class="post-metadata">

**Author:** ![Disla](https://yyz2.discourse-cdn.com/free1/user_avatar/calculix.discourse.group/disla/32/3548_2.png) [@Disla](https://calculix.discourse.group/u/Disla)\
**Post date:** [October 8, 2023, 7:39am UTC](https://calculix.discourse.group/t/von-mises-vs-peeq/1855/33 "2023-10-08T07:39:45Z")

</div>

> [@Juseche](#):
>
> Could this affects the results?

I think so.

> [@Juseche](#):
>
> how to plotting

Check this post.

> **[Stresses at the integration points in ParaView](https://prepomax.discourse.group/t/stresses-at-the-integration-points-in-paraview/797)**
>
> One of the limitations of CalculiX is that it can only output averaged nodal stresses in .frd files. If you need them at the integration points (where they are originally calculated and most accurate) you can use .dat file output but that’s not...

---

<div class="post-metadata">

**Author:** ![Juseche](https://avatars.discourse-cdn.com/v4/letter/j/258eb7/32.png) [@Juseche](https://calculix.discourse.group/u/Juseche)\
**Post date:** [October 19, 2023, 9:07pm UTC](https://calculix.discourse.group/t/von-mises-vs-peeq/1855/34 "2023-10-19T21:07:23Z")

</div>

Thank a lot for your help!

Jairo

[Previous page](https://calculix.discourse.group/t/von-mises-vs-peeq/1855.md?page=1)
