by J. Beaurain, EDF R&D / AMA
To study the durability of nuclear facilities and to represent the irreversible damage of materials, softening constitutive laws are generally used. For that kind of laws, the uniqueness of the solution is lost so that multiple solutions which respect the first order mechanical evolution of the structures could emerge.
These solutions can however be found using alternative solvers or continuation methods. Taking into account that possible existence of multiple solutions, the stability of numerical solutions is studied by means of the second order mechanical evolution to assure their physical observance.
The criterion we choose to define the stability is the positivity of the second derivative of the total energy considering only evolutions which increase damage, to respect the unilateral constraint. This leads to the minimization of a quadratic quantity, subjected to inequality constraints of positivity.
Several algorithms of inequality constraints minimization exist in the literature but their efficiency is limited in the case of important number of degrees of freedom. We overcome that limitation in Code_Aster by using the power method with projection combined with the Sorensen model reduction based on sub-space decomposition which contains the most relevant eigenvectors. To use that feature, users have to declared the option DDL_STAB into the analysis operator CRIT_STAB of STAT_NON_LINE.
So that it is now possible to:
This work has resulted in three scientific papers in national and international conferences :