# Higer stress from results than in elasto-plastic curve

**URL:** <https://prepomax.discourse.group/t/higer-stress-from-results-than-in-elasto-plastic-curve/1115>\
**Category:** General Questions\
**Tags:** material-definition\
**Created:** [September 25, 2023, 3:25pm UTC](https://prepomax.discourse.group/t/higer-stress-from-results-than-in-elasto-plastic-curve/1115 "2023-09-25T15:25:27Z")\
**Posts on this page:** 13\
**Page:** 1

<div class="post-metadata">

**Author:** ![lorebergo](https://avatars.discourse-cdn.com/v4/letter/l/ecc23a/32.png) [@lorebergo](https://prepomax.discourse.group/u/lorebergo)\
**Post date:** [September 25, 2023, 3:25pm UTC](https://prepomax.discourse.group/t/higer-stress-from-results-than-in-elasto-plastic-curve/1115/1 "2023-09-25T15:25:27Z")

</div>

Dear all, I am using elasto-plastic material curve with these properties:

 ![PrePoMax_nbmoFLjofc](https://global.discourse-cdn.com/free1/uploads/prepomax/original/2X/8/82ebbac8828aa177ea2c61afe7610b9f3200ed3a.png)  
Maximum stress is roughly 111 MPa. When I solve the model (nonlinear analysis, nnlin geom ON), I get maximum stress values (both Von Mises and principal maximum) higer then the highest value in the stress-strain curve (this affect a wide part of the model, it isn’t just a few elements). This sound strange to me, expecially because the material has a perfect plastic behaviour from a certain strain. Here the two images of the model:  
**Principal stress**  
 ![GBCMyveLQL](https://global.discourse-cdn.com/free1/uploads/prepomax/original/2X/4/4864ac5c78863c715715a9003a776383d2539e83.png)  
**Von Mises**  
 ![PrePoMax_tcbLDXTHLG](https://global.discourse-cdn.com/free1/uploads/prepomax/original/2X/2/2b282e1679868adcb7907ade06dad77c061d4dd3.png)

Thanks in advance for any help!

---

<div class="post-metadata">

**Author:** ![FEAnalyst](https://yyz2.discourse-cdn.com/free1/user_avatar/prepomax.discourse.group/feanalyst/32/688_2.png) [@FEAnalyst](https://prepomax.discourse.group/u/FEAnalyst)\
**Post date:** [September 25, 2023, 4:02pm UTC](https://prepomax.discourse.group/t/higer-stress-from-results-than-in-elasto-plastic-curve/1115/2 "2023-09-25T16:02:38Z")

</div>

In the majority of cases, this is caused by the extrapolation of stresses from the integration points to nodes. You would have to check the stresses directly at the integration points (it can be done in the .dat file or with a script: [Stresses at the integration points in ParaView - #2 by mkraska](https://prepomax.discourse.group/t/stresses-at-the-integration-points-in-paraview/797/2)).

---

<div class="post-metadata">

**Author:** ![synt](https://yyz2.discourse-cdn.com/free1/user_avatar/prepomax.discourse.group/synt/32/4658_2.png) [@synt](https://prepomax.discourse.group/u/synt)\
**Post date:** [September 25, 2023, 4:13pm UTC](https://prepomax.discourse.group/t/higer-stress-from-results-than-in-elasto-plastic-curve/1115/3 "2023-09-25T16:13:59Z")

</div>

> [@lorebergo](#):
>
> This sound strange to me, expecially because the material has a perfect plastic behaviour from a certain strain.

when magnify is high, it can be used as sign to do mesh convergence study. refinement is required in these areas. still discrepancy will exist due to extrapolation.

---

<div class="post-metadata">

**Author:** ![lorebergo](https://avatars.discourse-cdn.com/v4/letter/l/ecc23a/32.png) [@lorebergo](https://prepomax.discourse.group/u/lorebergo)\
**Post date:** [September 26, 2023, 7:26pm UTC](https://prepomax.discourse.group/t/higer-stress-from-results-than-in-elasto-plastic-curve/1115/4 "2023-09-26T19:26:22Z")

</div>

@synt thank you for your advice: as you suggested I adopted a finer mesh in the area of interest, going from an average el. size of 0.5mm down to 0.1mm: results didn’t change, I still have values around 140MPa instead of the maximum value in the curve, 111 MPa.  
Moreover, I tried with another material (same model, same loads, just different material) less ductile than the first I used, again adoptin gelasto-plastic behaviour such as this:

 ![image](https://global.discourse-cdn.com/free1/uploads/prepomax/original/2X/9/9b3a25e8184ff6bd021e9fe8b0f8b71713c977d0.png)  
Young’s modulus is roughly 3 times higher than the initial material.

@FEAnalyst I see your point and I would agree if for one reason: maximum stress in stress-strain curve is lower than the value I should find in any of the integration points. Interpolating integration points to the nodes, should give me at least not accurate stress, but I don’t expect this huge difference (110 vs 140MPa in wide areas of the model is roughly 25% of difference!!!)

---

<div class="post-metadata">

**Author:** ![ANYS](https://yyz2.discourse-cdn.com/free1/user_avatar/prepomax.discourse.group/anys/32/7737_2.png) [@ANYS](https://prepomax.discourse.group/u/ANYS)\
**Post date:** [September 26, 2023, 8:42pm UTC](https://prepomax.discourse.group/t/higer-stress-from-results-than-in-elasto-plastic-curve/1115/5 "2023-09-26T20:42:28Z")

</div>

> [@lorebergo](#):
>
> (both Von Mises and principal maximum)

Hi Lorebergo,

¿What stress are you comparing with?. Not clear to me as your curve says “DATA”.

---

<div class="post-metadata">

**Author:** ![synt](https://yyz2.discourse-cdn.com/free1/user_avatar/prepomax.discourse.group/synt/32/4658_2.png) [@synt](https://prepomax.discourse.group/u/synt)\
**Post date:** [September 26, 2023, 8:45pm UTC](https://prepomax.discourse.group/t/higer-stress-from-results-than-in-elasto-plastic-curve/1115/6 "2023-09-26T20:45:11Z")

</div>

> [@lorebergo](#):
>
> I adopted a finer mesh in the area of interest, going from an average el. size of 0.5mm down to 0.1mm: results didn’t change, I still have values around 140MPa instead of the maximum value in the curve, 111 MPa.

is maximum plastic strain located at the same nodes? if not, probably cause of over-constraint in model and simplification, thus unrealistic results.

---

<div class="post-metadata">

**Author:** ![FEAnalyst](https://yyz2.discourse-cdn.com/free1/user_avatar/prepomax.discourse.group/feanalyst/32/688_2.png) [@FEAnalyst](https://prepomax.discourse.group/u/FEAnalyst)\
**Post date:** [September 26, 2023, 8:52pm UTC](https://prepomax.discourse.group/t/higer-stress-from-results-than-in-elasto-plastic-curve/1115/7 "2023-09-26T20:52:18Z")

</div>

Check this thread as well, it’s an interesting discussion about a similar problem: [Load variation in non-linear material problem](https://prepomax.discourse.group/t/load-variation-in-non-linear-material-problem/417)

---

<div class="post-metadata">

**Author:** ![ANYS](https://yyz2.discourse-cdn.com/free1/user_avatar/prepomax.discourse.group/anys/32/7737_2.png) [@ANYS](https://prepomax.discourse.group/u/ANYS)\
**Post date:** [September 26, 2023, 8:52pm UTC](https://prepomax.discourse.group/t/higer-stress-from-results-than-in-elasto-plastic-curve/1115/8 "2023-09-26T20:52:44Z")

</div>

> [@lorebergo](#):
>
> nnlin geom ON

Maybe your plastic card is not right ?¿? Make sure you wrote it correctly and there are not warnings on the monitor like:

\*WARNING reading \*STEP: parameter not recognized:  
NNLINGEOM=ON  
\*WARNING reading \*STEP. Card image:  
\*STEP,NNLINGEOM=ON

It should say  
NLGEOM=ON

Sometimes, if the error is not critical the analysis can end up but with unexpected result.

---

<div class="post-metadata">

**Author:** ![FEAnalyst](https://yyz2.discourse-cdn.com/free1/user_avatar/prepomax.discourse.group/feanalyst/32/688_2.png) [@FEAnalyst](https://prepomax.discourse.group/u/FEAnalyst)\
**Post date:** [September 26, 2023, 8:57pm UTC](https://prepomax.discourse.group/t/higer-stress-from-results-than-in-elasto-plastic-curve/1115/9 "2023-09-26T20:57:19Z")

</div>

> [@ANYS](#):
>
> your curve says “DATA”

It’s because temperature is plotted in the same graph so there are different quantities on the Y axis. But temperature can be ignored here so it’s just yield stress - plastic strain plot of the input data for \*PLASTIC.

---

<div class="post-metadata">

**Author:** ![synt](https://yyz2.discourse-cdn.com/free1/user_avatar/prepomax.discourse.group/synt/32/4658_2.png) [@synt](https://prepomax.discourse.group/u/synt)\
**Post date:** [September 26, 2023, 10:03pm UTC](https://prepomax.discourse.group/t/higer-stress-from-results-than-in-elasto-plastic-curve/1115/10 "2023-09-26T22:03:28Z")

</div>

> [@ANYS](#):
>
> \*WARNING reading \*STEP: parameter not recognized:  
> NNLINGEOM=ON  
> \*WARNING reading \*STEP. Card image:  
> \*STEP,NNLINGEOM=ON

it seems not possible to do typo, since these feature is available in menus.of PrePoMax

also, being restricted to modify accidentally by the user in Edit CalculiX Keyword.

any reason to do manually typing or modify/overwrites even possible?

 ![2023-09-27 04_55_18-Edit Step_ Step-1](https://global.discourse-cdn.com/free1/uploads/prepomax/original/2X/7/75584780f648f670f2336e9abadc868ae67fc366.png)

\*edited  
when typo as guessed, there will be no results shown and a warning message appears.

 ![2023-09-27 07_11_43-Monitor](https://global.discourse-cdn.com/free1/uploads/prepomax/original/2X/e/e43fdbb77491cefb768876e942b19d49e0a455b7.png)

---

<div class="post-metadata">

**Author:** ![Matej](https://avatars.discourse-cdn.com/v4/letter/m/a698b9/32.png) [@Matej](https://prepomax.discourse.group/u/Matej)\
**Post date:** [September 27, 2023, 6:54am UTC](https://prepomax.discourse.group/t/higer-stress-from-results-than-in-elasto-plastic-curve/1115/11 "2023-09-27T06:54:27Z")

</div>

If you use the Field Output (on one of the most loaded elements), you will get the stress results in the integration points. Check this stress to see its real value.

---

<div class="post-metadata">

**Author:** ![FEAnalyst](https://yyz2.discourse-cdn.com/free1/user_avatar/prepomax.discourse.group/feanalyst/32/688_2.png) [@FEAnalyst](https://prepomax.discourse.group/u/FEAnalyst)\
**Post date:** [September 27, 2023, 7:02am UTC](https://prepomax.discourse.group/t/higer-stress-from-results-than-in-elasto-plastic-curve/1115/12 "2023-09-27T07:02:57Z")

</div>

> [@Matej](#):
>
> If you use the Field Output (on one of the most loaded elements)

Don’t you mean history output ?

---

<div class="post-metadata">

**Author:** ![Matej](https://avatars.discourse-cdn.com/v4/letter/m/a698b9/32.png) [@Matej](https://prepomax.discourse.group/u/Matej)\
**Post date:** [September 27, 2023, 7:44am UTC](https://prepomax.discourse.group/t/higer-stress-from-results-than-in-elasto-plastic-curve/1115/13 "2023-09-27T07:44:03Z")

</div>

Yes, I am sorry, the History output.
