0% found this document useful (0 votes)

44 views

Loading Path For A Structure: Code - Aster, Salome-Meca Course Material

1) The document discusses errors that can occur when numerically integrating elastoplastic constitutive laws with implicit schemes. It presents examples of loading paths applied to a material point to illustrate the issues. 2) One source of error is approximating the plastic strain increment, which can be mitigated by refining the time discretization when loading is not perfectly radial or there are load reversals within a time step. 3) An indicator called RADI_T is proposed to measure how non-radial the loading path is, to help determine when time discretization needs to be refined to reduce integration errors.

Uploaded by

Stefano Milani

Available Formats

Download as PDF, TXT or read online on Scribd

0% found this document useful (0 votes)

44 views

Loading Path For A Structure: Code - Aster, Salome-Meca Course Material

Uploaded by

Stefano Milani

Available Formats

Download as PDF, TXT or read online on Scribd

You are on page 1/ 15

Loading path for a structure

Code_Aster, Salome-Meca course material

GNU FDL licence (http://www.gnu.org/copyleft/fdl.html)
Problem to be solved

The chosen model is suited to loading path ?

What is the error in the integration of elastoplastic constitutive

law ?

How to avoid these errors?

2 - Code_Aster and Salome-Meca course material GNU FDL Licence

Problem to be solved

Know the limitations of theoretical models used

Traction loading + 90 ° phase shifted torsion

3 - Code_Aster and Salome-Meca course material GNU FDL Licence

A (very) simple problem

Material point subjected to a loading path


Elastoplasticity with linear isotropic hardening
y  ET
 Y
ET E = 195000 MPa
1 2
Y  = 0.3
E = 195000 MPa
= 0.3
Y = 181 MPa
1 mm E
EY T= 181
= 1930
MPa MPa
E 
3 4 x E T = 1930 MPa
1 mm


Loading
D D D
D A (MPa) (MPa)
d D A 151.2 93.1
D B
D D
B 257.2 33.1
0 
C C 259.3 0

4 - Code_Aster and Salome-Meca course material GNU FDL Licence

Time discretization error

3 methods of resolution
Analytical resolution -> reference values for  xy
Model A : coarse time discretization
Model B : fine time discretization
(between two time steps, the cumulated plastic strain is constrained not to increase no
more than 0.1 10-2)

Relative error of XY

A
Time step Model A Model B

1 s (point A) 0,003 % 0,003%

2 s (point B) 4,3% 1,1 % B C

5 - Code_Aster and Salome-Meca course material GNU FDL Licence

Implicit discretization

Continuous Discrete

(t )  (t ) (ti )  i  (ti )  i

      i  
       
 i 
Constitutive law
 
   ( , ) Choice of the  i 1   ( i 1 ,  i 1 ,  i ,  i )
scheme

1)   A : (   ) p
R  R( p) 1)    A : (    p  ) R   R( p  )
3  ~  p 3    p 
Isotropic
2)   
p
p   2)  
hardening +
2  eq t 2 t  eq t t
implicit scheme
3) f  0   0 f 0 3)  f   0   0 f 0

Implicit order 1 scheme

6 - Code_Aster and Salome-Meca course material GNU FDL Licence

Error

1
t  dp 
p  t  dt   d   OK
dt

t  d t  d
p
2 
  p
t  dt  t  n( )dt   n( )d
dt dt

 
    

n( )  n' ( ).(   ) d  n( )  n' ( ).
2
2
 

Error!

7 - Code_Aster and Salome-Meca course material GNU FDL Licence

Error

3 : if   0 f  0 from t - to t+ : OK, except if : load + undload in 1

single time step.

Loading path calculated


Actual loading path

8 - Code_Aster and Salome-Meca course material GNU FDL Licence

Error

 2
Main source of error :  n' ( ).
2

Exact integration
Permanent load during the loading step
The normal to the criterion is not reversing ("radial" load )

Refine the time step outside these two cases!

But how to detect?

9 - Code_Aster and Salome-Meca course material GNU FDL Licence

Evaluation of the error

Relative error on integration :

2 
 p n ' ( ) n ' ( )
  2  2
 p n( )  1

n  1
  n( )
 2 2

10 - Code_Aster and Salome-Meca course material GNU FDL Licence

An indicator of loss of "radiality"

1  sin q The error is measured

  n  sin 
2 2 2 by  and q
One define DERA_ELGA(CALC_CHAMP)
2 
RADI_T = sin
2
~
n RADI_V = 1 cosq 
 
Equals 0 when the loading is radial
q
 n
~

1
(  ) eq ( ) eq


One define RESI_RADI_RELA

(COMP_INCR)
RADI_T
VMIS_xx
11 - Code_Aster and Salome-Meca course material GNU FDL Licence
Sufficient conditions of radiality

"Radiality" = when the threshold normal is not reversing

 ( x, t )   (t )  0 ( x)
Radiality conditions
External loads are proportional
Initial state undeformed and not hardened
Von Mises plasticity with isotropic hardening R( p)  Cp
n

Advantages
Exact local integration
Easy global convergence