Good afternoon
Is possible simulate on calculix Creep with temperature on time line?
Thank you
Yes, you can set the dependency of creep properties on temperature. Check the *CREEP keyword in the documentation. Here’s an example:
*CREEP
1.E-10,5.,0.,100.
2.E-10,5.,0.,200.
The last item on the data line is temperature (100-200).
1 Like
Thank you calc_em
Is possible by 2 step simulate residual stress + Creep with thermal flux for time very long?
It should be possible with the *Coupled temperature-displacement step. Give it a try.
I have try to create 2 step:
step-1 displacement to deform the piece.
step-2 remove displacement apply surface temperature (500°C)
this generate an error
*ERROR in CalculiX: in nonlinear calculations
energy output must be selected in the first step
where i wrong?
**
** Heading +++++++++++++++++++++++++++++++++++++++++++++++++
**
*Heading
Hash: PpYJYxo1, Date: 09/06/2022, Unit system: MM_TON_S_C
**
** Nodes +++++++++++++++++++++++++++++++++++++++++++++++++++
**
*Node
... hidden data ...
**
** Elements ++++++++++++++++++++++++++++++++++++++++++++++++
**
*Element, Type=C3D10, Elset=Solid_part-1
... hidden data ...
**
** Node sets +++++++++++++++++++++++++++++++++++++++++++++++
**
*Nset, Nset=Internal-1_fix
... hidden data ...
*Nset, Nset=Internal-1_load
... hidden data ...
*Nset, Nset=Internal_Selection-1_Temperature-1
... hidden data ...
*Nset, Nset=Internal_Selection-1_Concentrated_flux-1
... hidden data ...
*Nset, Nset=Internal_Selection-2_Temperature-1
... hidden data ...
**
** Element sets ++++++++++++++++++++++++++++++++++++++++++++
**
*Elset, Elset=Internal-1_fix_S4
... hidden data ...
*Elset, Elset=Internal-1_fix_S1
... hidden data ...
*Elset, Elset=Internal-1_fix_S3
... hidden data ...
*Elset, Elset=Internal-1_fix_S2
... hidden data ...
*Elset, Elset=Internal-1_load_S4
... hidden data ...
*Elset, Elset=Internal-1_load_S2
... hidden data ...
*Elset, Elset=Internal-1_load_S3
... hidden data ...
*Elset, Elset=Internal-1_load_S1
... hidden data ...
**
** Surfaces ++++++++++++++++++++++++++++++++++++++++++++++++
**
*Surface, Name=fix, Type=Element
Internal-1_fix_S4, S4
Internal-1_fix_S1, S1
Internal-1_fix_S3, S3
Internal-1_fix_S2, S2
*Surface, Name=load, Type=Element
Internal-1_load_S4, S4
Internal-1_load_S2, S2
Internal-1_load_S3, S3
Internal-1_load_S1, S1
**
** Physical constants ++++++++++++++++++++++++++++++++++++++
**
**
** Materials +++++++++++++++++++++++++++++++++++++++++++++++
**
*Material, Name=S355_mio
*Density
7.8E-09
*Elastic
210000, 0.28
*Expansion, Zero=20
1.1E-05
*Conductivity
14
*Specific heat
440000000
*Plastic
355, 0.002
450, 0.0025
550, 0.006
600, 0.015
**
** Sections ++++++++++++++++++++++++++++++++++++++++++++++++
**
*Solid section, Elset=Solid_part-1, Material=S355_mio
**
** Pre-tension sections ++++++++++++++++++++++++++++++++++++
**
**
** Constraints +++++++++++++++++++++++++++++++++++++++++++++
**
**
** Surface interactions ++++++++++++++++++++++++++++++++++++
**
**
** Contact pairs +++++++++++++++++++++++++++++++++++++++++++
**
**
** Amplitudes ++++++++++++++++++++++++++++++++++++++++++++++
**
*Amplitude, Name=Tabular-1
0, 0, 0.5, 1, 1, 0
*Amplitude, Name=Tabular-2
0, 0, 0.5, 1, 1, 0
*Amplitude, Name=Tabular-3
0, 500, 500, 500
**
** Initial conditions ++++++++++++++++++++++++++++++++++++++
**
** Name: Temperature-1
*Initial conditions, Type=Temperature
Internal_Selection-1_Temperature-1, 20
**
** Steps +++++++++++++++++++++++++++++++++++++++++++++++++++
**
**
** Step-1 ++++++++++++++++++++++++++++++++++++++++++++++++++
**
*Step, Nlgeom, Inc=100
*Static
0.001, 1, 1E-05, 0.02
**
** Boundary conditions +++++++++++++++++++++++++++++++++++++
**
*Boundary, op=New
** Name: Fixed-1
*Boundary
Internal-1_fix, 1, 6, 0
** Name: Displacement_rotation-1
*Boundary
Internal-1_load, 2, 2, 0.3
**
** Loads +++++++++++++++++++++++++++++++++++++++++++++++++++
**
*Cload, op=New
*Dload, op=New
** Name: Surface_traction-1: Deactivated
**
** Defined fields ++++++++++++++++++++++++++++++++++++++++++
**
**
** History outputs +++++++++++++++++++++++++++++++++++++++++
**
*Node print, Nset=Internal-1_fix, Totals=Yes
RF
*Node print, Nset=Internal-1_load, Totals=Yes
RF
*Node print, Nset=Internal-1_fix, Totals=Only
RF
**
** Field outputs +++++++++++++++++++++++++++++++++++++++++++
**
*Node file
RF, U
*El file
S, E, PEEQ
**
** End step ++++++++++++++++++++++++++++++++++++++++++++++++
**
*End step
**
** Step-2 ++++++++++++++++++++++++++++++++++++++++++++++++++
**
*Step, Inc=100
*Coupled temperature-displacement, Direct
1, 500
**
** Boundary conditions +++++++++++++++++++++++++++++++++++++
**
*Boundary, op=New
** Name: Fixed-1
*Boundary
Internal-1_fix, 1, 6, 0
** Name: Displacement_rotation-1: Deactivated
** Name: Temperature-1
*Boundary
Internal_Selection-2_Temperature-1, 11, 11, 500
**
** Loads +++++++++++++++++++++++++++++++++++++++++++++++++++
**
*Cload, op=New
*Dload, op=New
*Cflux, op=New
** Name: Surface_traction-1: Deactivated
** Name: Concentrated_flux-1: Deactivated
**
** Defined fields ++++++++++++++++++++++++++++++++++++++++++
**
**
** History outputs +++++++++++++++++++++++++++++++++++++++++
**
**
** Field outputs +++++++++++++++++++++++++++++++++++++++++++
**
*Node file
RF, U, NT, RFL
*El file
PEEQ
**
** End step ++++++++++++++++++++++++++++++++++++++++++++++++
**
*End step
Try adding ENER to *EL FILE in the first step.
1 Like
thank you is better, but now there is an error:
*ERROR: increment size smaller than minimum
How I can repair?
**
** Heading +++++++++++++++++++++++++++++++++++++++++++++++++
**
*Heading
Hash: PpYJYxo1, Date: 09/07/2022, Unit system: MM_TON_S_C
**
** Nodes +++++++++++++++++++++++++++++++++++++++++++++++++++
**
*Node
... hidden data ...
**
** Elements ++++++++++++++++++++++++++++++++++++++++++++++++
**
*Element, Type=C3D10, Elset=Solid_part-1
... hidden data ...
**
** Node sets +++++++++++++++++++++++++++++++++++++++++++++++
**
*Nset, Nset=Internal-1_fix
... hidden data ...
*Nset, Nset=Internal-1_load
... hidden data ...
*Nset, Nset=Internal_Selection-1_Temperature-1
... hidden data ...
*Nset, Nset=Internal_Selection-1_Concentrated_flux-1
... hidden data ...
*Nset, Nset=Internal_Selection-2_Temperature-1
... hidden data ...
*Nset, Nset=Internal-1_Internal_Selection-1_Convective_film-1
... hidden data ...
**
** Element sets ++++++++++++++++++++++++++++++++++++++++++++
**
*Elset, Elset=Internal-1_fix_S4
... hidden data ...
*Elset, Elset=Internal-1_fix_S1
... hidden data ...
*Elset, Elset=Internal-1_fix_S3
... hidden data ...
*Elset, Elset=Internal-1_fix_S2
... hidden data ...
*Elset, Elset=Internal-1_load_S4
... hidden data ...
*Elset, Elset=Internal-1_load_S2
... hidden data ...
*Elset, Elset=Internal-1_load_S3
... hidden data ...
*Elset, Elset=Internal-1_load_S1
... hidden data ...
*Elset, Elset=Internal-1_Internal_Selection-1_Convective_film-1_S1
... hidden data ...
*Elset, Elset=Internal-1_Internal_Selection-1_Convective_film-1_S4
... hidden data ...
*Elset, Elset=Internal-1_Internal_Selection-1_Convective_film-1_S3
... hidden data ...
*Elset, Elset=Internal-1_Internal_Selection-1_Convective_film-1_S2
... hidden data ...
**
** Surfaces ++++++++++++++++++++++++++++++++++++++++++++++++
**
*Surface, Name=fix, Type=Element
Internal-1_fix_S4, S4
Internal-1_fix_S1, S1
Internal-1_fix_S3, S3
Internal-1_fix_S2, S2
*Surface, Name=load, Type=Element
Internal-1_load_S4, S4
Internal-1_load_S2, S2
Internal-1_load_S3, S3
Internal-1_load_S1, S1
*Surface, Name=Internal_Selection-1_Convective_film-1, Type=Element
Internal-1_Internal_Selection-1_Convective_film-1_S1, S1
Internal-1_Internal_Selection-1_Convective_film-1_S4, S4
Internal-1_Internal_Selection-1_Convective_film-1_S3, S3
Internal-1_Internal_Selection-1_Convective_film-1_S2, S2
**
** Physical constants ++++++++++++++++++++++++++++++++++++++
**
**
** Materials +++++++++++++++++++++++++++++++++++++++++++++++
**
*Material, Name=S355_mio
*Density
7.8E-09
*Elastic
210000, 0.28
*Expansion, Zero=20
1.1E-05
*Conductivity
14
*Specific heat
440000000
*Plastic
355, 0.002
450, 0.0025
550, 0.006
600, 0.015
**
** Sections ++++++++++++++++++++++++++++++++++++++++++++++++
**
*Solid section, Elset=Solid_part-1, Material=S355_mio
**
** Pre-tension sections ++++++++++++++++++++++++++++++++++++
**
**
** Constraints +++++++++++++++++++++++++++++++++++++++++++++
**
**
** Surface interactions ++++++++++++++++++++++++++++++++++++
**
**
** Contact pairs +++++++++++++++++++++++++++++++++++++++++++
**
**
** Amplitudes ++++++++++++++++++++++++++++++++++++++++++++++
**
*Amplitude, Name=Tabular-1
0, 0, 0.5, 1, 1, 0
*Amplitude, Name=Tabular-2
0, 0, 0.5, 1, 1, 0
*Amplitude, Name=Tabular-3
0, 500, 500, 500
**
** Initial conditions ++++++++++++++++++++++++++++++++++++++
**
** Name: Temperature-1
*Initial conditions, Type=Temperature
Internal_Selection-1_Temperature-1, 20
**
** Steps +++++++++++++++++++++++++++++++++++++++++++++++++++
**
**
** Step-1 ++++++++++++++++++++++++++++++++++++++++++++++++++
**
*Step, Nlgeom, Inc=100
*Static
0.001, 1, 1E-05, 0.02
**
** Boundary conditions +++++++++++++++++++++++++++++++++++++
**
*Boundary, op=New
** Name: Fixed-1
*Boundary
Internal-1_fix, 1, 6, 0
** Name: Displacement_rotation-1
*Boundary
Internal-1_load, 2, 2, 0.3
**
** Loads +++++++++++++++++++++++++++++++++++++++++++++++++++
**
*Cload, op=New
*Dload, op=New
** Name: Surface_traction-1: Deactivated
**
** Defined fields ++++++++++++++++++++++++++++++++++++++++++
**
**
** History outputs +++++++++++++++++++++++++++++++++++++++++
**
*Node print, Nset=Internal-1_fix, Totals=Yes
RF
*Node print, Nset=Internal-1_load, Totals=Yes
RF
*Node print, Nset=Internal-1_fix, Totals=Only
RF
**
** Field outputs +++++++++++++++++++++++++++++++++++++++++++
**
*Node file
RF, U
*El file
S, E, PEEQ, ENER
**
** End step ++++++++++++++++++++++++++++++++++++++++++++++++
**
*End step
**
** Step-2 ++++++++++++++++++++++++++++++++++++++++++++++++++
**
*Step, Inc=100
*Coupled temperature-displacement
1, 500, 1E-05, 1E+30
**
** Boundary conditions +++++++++++++++++++++++++++++++++++++
**
*Boundary, op=New
** Name: Fixed-1
*Boundary
Internal-1_fix, 1, 6, 0
** Name: Displacement_rotation-1: Deactivated
** Name: Temperature-1: Deactivated
**
** Loads +++++++++++++++++++++++++++++++++++++++++++++++++++
**
*Cload, op=New
*Dload, op=New
*Cflux, op=New
*Film, op=New
** Name: Surface_traction-1: Deactivated
** Name: Concentrated_flux-1: Deactivated
** Name: Convective_film-1
*Film
Internal-1_Internal_Selection-1_Convective_film-1_S1, F1, 120, 10
Internal-1_Internal_Selection-1_Convective_film-1_S4, F4, 120, 10
Internal-1_Internal_Selection-1_Convective_film-1_S3, F3, 120, 10
Internal-1_Internal_Selection-1_Convective_film-1_S2, F2, 120, 10
**
** Defined fields ++++++++++++++++++++++++++++++++++++++++++
**
**
** History outputs +++++++++++++++++++++++++++++++++++++++++
**
**
** Field outputs +++++++++++++++++++++++++++++++++++++++++++
**
*Node file
RF, U, NT, RFL
*El file
PEEQ
**
** End step ++++++++++++++++++++++++++++++++++++++++++++++++
**
*End step
That’s a typical convergence issue. Try simplifying the model (removing features one by one) to see what causes the problem. It can be Nlgeom, too large load or incorrect boundary conditions (also thermal ones) among the others.
Creep is not included yet intentionally, right ?
1 Like
I’m try to semplify model and BC.
Yes the Creep is included intentionally, why Is wrong?
Thank you
Now I have try to modify BC but whit this error:
*ERROR: max. # of increments reached
Is better use *visco configuration or is wrong for creep simulation?
Here the code:
**
** Heading +++++++++++++++++++++++++++++++++++++++++++++++++
**
*Heading
Hash: PpYJYxo1, Date: 09/14/2022, Unit system: MM_TON_S_C
**
** Nodes +++++++++++++++++++++++++++++++++++++++++++++++++++
**
*Node
... hidden data ...
**
** Elements ++++++++++++++++++++++++++++++++++++++++++++++++
**
*Element, Type=C3D10, Elset=Solid_part-1
... hidden data ...
**
** Node sets +++++++++++++++++++++++++++++++++++++++++++++++
**
*Nset, Nset=Internal_Selection-1_Temperature-1
... hidden data ...
*Nset, Nset=Internal_Selection-1_Concentrated_flux-1
... hidden data ...
*Nset, Nset=Internal-1_fix
... hidden data ...
*Nset, Nset=Internal-1_load
... hidden data ...
*Nset, Nset=Internal-1_Internal_Selection-1_Convective_film-1
... hidden data ...
*Nset, Nset=Internal_Selection-2_Temperature-1
... hidden data ...
**
** Element sets ++++++++++++++++++++++++++++++++++++++++++++
**
*Elset, Elset=Internal-1_fix_S3
... hidden data ...
*Elset, Elset=Internal-1_fix_S2
... hidden data ...
*Elset, Elset=Internal-1_fix_S4
... hidden data ...
*Elset, Elset=Internal-1_load_S2
... hidden data ...
*Elset, Elset=Internal-1_load_S4
... hidden data ...
*Elset, Elset=Internal-1_load_S3
... hidden data ...
*Elset, Elset=Internal-1_load_S1
... hidden data ...
*Elset, Elset=Internal-1_Internal_Selection-1_Convective_film-1_S4
... hidden data ...
*Elset, Elset=Internal-1_Internal_Selection-1_Convective_film-1_S2
... hidden data ...
*Elset, Elset=Internal-1_Internal_Selection-1_Convective_film-1_S1
... hidden data ...
*Elset, Elset=Internal-1_Internal_Selection-1_Convective_film-1_S3
... hidden data ...
**
** Surfaces ++++++++++++++++++++++++++++++++++++++++++++++++
**
*Surface, Name=fix, Type=Element
Internal-1_fix_S3, S3
Internal-1_fix_S2, S2
Internal-1_fix_S4, S4
*Surface, Name=load, Type=Element
Internal-1_load_S2, S2
Internal-1_load_S4, S4
Internal-1_load_S3, S3
Internal-1_load_S1, S1
*Surface, Name=Internal_Selection-1_Convective_film-1, Type=Element
Internal-1_Internal_Selection-1_Convective_film-1_S4, S4
Internal-1_Internal_Selection-1_Convective_film-1_S2, S2
Internal-1_Internal_Selection-1_Convective_film-1_S1, S1
Internal-1_Internal_Selection-1_Convective_film-1_S3, S3
**
** Physical constants ++++++++++++++++++++++++++++++++++++++
**
**
** Materials +++++++++++++++++++++++++++++++++++++++++++++++
**
*Material, Name=S355_mio
*Density
7.8E-09
*Elastic
210000, 0.28
*Expansion, Zero=20
1.1E-05
*Conductivity
14
*Specific heat
440000000
*Plastic
355, 0.002
450, 0.0025
550, 0.006
600, 0.015
**
** Sections ++++++++++++++++++++++++++++++++++++++++++++++++
**
*Solid section, Elset=Solid_part-1, Material=S355_mio
**
** Pre-tension sections ++++++++++++++++++++++++++++++++++++
**
**
** Constraints +++++++++++++++++++++++++++++++++++++++++++++
**
**
** Surface interactions ++++++++++++++++++++++++++++++++++++
**
**
** Contact pairs +++++++++++++++++++++++++++++++++++++++++++
**
**
** Amplitudes ++++++++++++++++++++++++++++++++++++++++++++++
**
*Amplitude, Name=Tabular-1
0, 0, 0.5, 1, 1, 0
*Amplitude, Name=Tabular-2
0, 0, 0.5, 1, 1, 0
*Amplitude, Name=Tabular-3
0, 500, 500, 500
**
** Initial conditions ++++++++++++++++++++++++++++++++++++++
**
** Name: Temperature-1
*Initial conditions, Type=Temperature
Internal_Selection-1_Temperature-1, 20
**
** Steps +++++++++++++++++++++++++++++++++++++++++++++++++++
**
**
** Step-1 ++++++++++++++++++++++++++++++++++++++++++++++++++
**
*Step, Inc=100
*Static
0.001, 1, 1E-05, 0.02
**
** Boundary conditions +++++++++++++++++++++++++++++++++++++
**
*Boundary, op=New
** Name: Fixed-1
*Boundary
Internal-1_fix, 1, 6, 0
** Name: Displacement_rotation-1
*Boundary
Internal-1_load, 2, 2, 0.3
**
** Loads +++++++++++++++++++++++++++++++++++++++++++++++++++
**
*Cload, op=New
*Dload, op=New
** Name: Surface_traction-1: Deactivated
**
** Defined fields ++++++++++++++++++++++++++++++++++++++++++
**
**
** History outputs +++++++++++++++++++++++++++++++++++++++++
**
*Node print, Nset=Internal-1_fix, Totals=Yes
RF
*Node print, Nset=Internal-1_load, Totals=Yes
RF
*Node print, Nset=Internal-1_fix, Totals=Only
RF
**
** Field outputs +++++++++++++++++++++++++++++++++++++++++++
**
*Node file
RF, U
*El file
S, E, PEEQ, ENER
**
** End step ++++++++++++++++++++++++++++++++++++++++++++++++
**
*End step
**
** Step-2 ++++++++++++++++++++++++++++++++++++++++++++++++++
**
*Step, Inc=200
*Coupled temperature-displacement
1E-10, 20, 1E-10, 1E+30
**
** Boundary conditions +++++++++++++++++++++++++++++++++++++
**
*Boundary, op=New
** Name: Fixed-1
*Boundary
Internal-1_fix, 1, 6, 0
** Name: Displacement_rotation-1: Deactivated
** Name: Temperature-1
*Boundary
Internal_Selection-2_Temperature-1, 11, 11, 400
**
** Loads +++++++++++++++++++++++++++++++++++++++++++++++++++
**
*Cload, op=New
*Dload, op=New
*Cflux, op=New
*Film, op=New
** Name: Surface_traction-1: Deactivated
** Name: Concentrated_flux-1: Deactivated
** Name: Convective_film-1: Deactivated
**
** Defined fields ++++++++++++++++++++++++++++++++++++++++++
**
**
** History outputs +++++++++++++++++++++++++++++++++++++++++
**
**
** Field outputs +++++++++++++++++++++++++++++++++++++++++++
**
*Node file
RF, U, NT, RFL
*El file
PEEQ
**
** End step ++++++++++++++++++++++++++++++++++++++++++++++++
**
*End step
Visco step is meant specifically for analyses involving creep but it can’t account for heat transfer.
Try increasing the maximum allowed number of increments. However, it might be hard to reach convergence anyway: Creep calculation and results