# Creep calculation and results

**URL:** <https://calculix.discourse.group/t/creep-calculation-and-results/868>\
**Category:** Uncategorized\
**Created:** [November 28, 2021, 12:14am UTC](https://calculix.discourse.group/t/creep-calculation-and-results/868 "2021-11-28T00:14:44Z")\
**Posts on this page:** 5\
**Page:** 2

<div class="post-metadata">

**Author:** ![ZFDong](https://avatars.discourse-cdn.com/v4/letter/z/838e76/32.png) [@ZFDong](https://calculix.discourse.group/u/ZFDong)\
**Post date:** [December 2, 2021, 10:48pm UTC](https://calculix.discourse.group/t/creep-calculation-and-results/868/22 "2021-12-02T22:48:46Z")

</div>

Calc\_em:

You can try this one, it can be solved by Abaqus, but even run a few step in Calculix, the results are not correct at all. Thanks.

\*HEADING  
BASIC VERIFICATION OF CREEP (CREEP TEST)  
\*NODE  
1,1.,10.,1.  
2,2.,10.,1.  
3,2.,10.,0.  
4,1.,10.,0.  
11,1.,5.,1.  
12,2.,5.,1.  
13,2.,5.,0.  
14,1.,5.,0.  
21,1.,0.,1.  
22,2.,0.,1.  
23,2.,0.,0.  
24,1.,0.,0.  
\*NSET,NSET=BOT  
21,22,23,24  
\*NSET,NSET=TOP  
1,2,3,4  
\*NSET,NSET=LEFT  
21,24  
\*ELEMENT,TYPE=C3D4  
1,21,22,23,12  
2,21,12,23,14  
3,21,14,23,24  
4,21,12,14,11  
5,14,12,23,13  
6,11,12,14,1  
7,1,12,14,3  
8,14,3,12,13  
9,12,3,1,2  
10,1,3,14,4  
\*ELSET,ELSET=DATTIM  
1,2,3,4,5,6,7,8,9,10  
\*ELSET,ELSET=CREEPERS  
DATTIM  
\*SOLID SECTION,ELSET=DATTIM,MATERIAL=A1  
\*MATERIAL,NAME= A1  
\*ELASTIC  
20E6,.3  
\*PLASTIC  
0.,0.  
0.,10.  
\*CREEP,LAW=NORTON  
2.5E-27,5.,-0.2  
\*\*2.5E-27,5.,0.0  
\*BOUNDARY  
BOT,1,1  
BOT,2,2  
BOT,3,3  
\*\*LEFT,1

* * *

\*STEP  
\*\* STEP 1 ----- ELASTIC LOADING  
\*STATIC  
1.0E-7,1.0E-7  
\*DLOAD  
1,P1,20000.  
2,P1,20000.  
\*EL FILE,FREQUENCY=20  
SDV,  
S,  
CEEQ  
\*NODE FILE, FREQUENCY=20  
U  
\*END STEP  
\*STEP,INC=2500000  
\*\* STEP 2 ----- CREEP TEST  
\*VISCO, CETOL=5.E-3  
1.e-8,100.,1.E-18  
\*\*PRINT,FREQUENCY=5,RESIDUAL=NO  
\*EL PRINT,ELSET=CREEPERS,FREQUENCY=20  
SDV,S  
CEEQ,  
\*EL FILE,ELSET=CREEPERS,FREQUENCY=20  
SDV,S,CEEQ,  
\*NODE FILE,FREQUENCY=20  
RF,U  
\*END STEP

---

<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:** [December 2, 2021, 11:22pm UTC](https://calculix.discourse.group/t/creep-calculation-and-results/868/23 "2021-12-02T23:22:39Z")

</div>

> [@ZFDong](#):
>
> \*PLASTIC  
> 0.,0.  
> 0.,10.

This definition won’t pass in Abaqus. And I would avoid such coarse meshes consisting of C3D4 elements for these tests.

To actually confirm the incorrectness of results obtained with CalculiX’s in this case, it is necessary to achieve convergence in both solvers (CalculiX and the reference one - Abaqus) and test different element types with a sufficiently dense mesh. I would expect C3D8R elements to be the most accurate but it’s also worth checking out other types. Here you can find my results of plastic bending test with different elements: [Elasto-plastic bending - elements, accuracy](https://calculix.discourse.group/t/elasto-plastic-bending-elements-accuracy/700)

---

<div class="post-metadata">

**Author:** ![ZFDong](https://avatars.discourse-cdn.com/v4/letter/z/838e76/32.png) [@ZFDong](https://calculix.discourse.group/u/ZFDong)\
**Post date:** [December 3, 2021, 12:57am UTC](https://calculix.discourse.group/t/creep-calculation-and-results/868/24 "2021-12-03T00:57:48Z")

</div>

Calc\_em:  
Thanks.  
I have checked part of CCX source code for using the Norton law’s contribution to the creep analysis.  
I only can find one location to use this, that is, umat\_aniso\_creep.f uses this formulation. But this function is for user material, ANISO\_CREEP material. I wonder how the CCX counts the Norton’s Law in the code for the default creep material.  
I have tried many ways to get the creep results. The results from CCX are all not reasonable.

---

<div class="post-metadata">

**Author:** ![demirdag](https://avatars.discourse-cdn.com/v4/letter/d/8edcca/32.png) [@demirdag](https://calculix.discourse.group/u/demirdag)\
**Post date:** [April 4, 2022, 8:47am UTC](https://calculix.discourse.group/t/creep-calculation-and-results/868/25 "2022-04-04T08:47:43Z")

</div>

Is there any other solver than Pardiso that working gut with Calculix? Or is it possible to get that plot with gnuplot from Calculix?

---

<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:** [September 8, 2022, 7:11am UTC](https://calculix.discourse.group/t/creep-calculation-and-results/868/26 "2022-09-08T07:11:30Z")

</div>

Hi,

I have start to search about Creep and I have recently found a very interesting explanation of why convergence could be difficult with Calculix in some cases. The origin could be the Creep LAW=NORTON formulation itself. (Explained in the first two minutes of the vid).

According to the manual the stress is not normalized in Calculix, and it could lead to a tiny A parameter which further leads to rounding/convergence issues. Reviewing the example shown in the previous posts of ZFDong , I can see he is using an A parameter of A=2.5E-27. That seems a bad starting point for any numerical analisys. Maybe other unit system for the model would help to reduce that A value.  
¿What do you think? ¿Or maybe directly normalizing the NORTON LAW and changing A by the initial Strain rate parameter in Calculix?

[![](https://global.discourse-cdn.com/free1/uploads/calculix/original/2X/e/eb0edff477e4022b1f57627dcb45a67c98d7e6f0.jpeg "Abaqus Creep Models: Which One Should You Use?") ](https://www.youtube.com/watch?v=W6bIDjkOX9o)

![screenshot.38](https://global.discourse-cdn.com/free1/uploads/calculix/original/1X/3144bc8528922578e00c8fa8163262433ef59615.png)

Regards,

[Previous page](https://calculix.discourse.group/t/creep-calculation-and-results/868.md?page=1)
