Code_Aster contact analysis - train rail

Hi,

Thanks again for the help, I have been getting to grips with google translate and been using the glossary a bit on the code aster website...

I have made the first few changes and attach the latest version of the study. I will attempt to find some time tonight to make further progress, I am using the version of code_aster in a flat CAELinux install, I believe 10.5, although 11 seems to be installed also, SalomeMeca seems to select 10.5 by default...

Kev

Hi kees,

Sorry there has been little going on with this, I have been away for a week. Since getting back I have done quite a bit, I have re-written the command file (eficas) and currently have what I think* is a nearly working study.

Code_aster exits with what looks like relative success, but no results seem to materialise, which I think is due to the mesh errors.

any help that could be given with this would be very gratefully received, if it is considered it would be worthwhile, I will write the final into a tutorial???

I have attached the files...

Many Thanks,

Kev

Hay Kees,

I have tried to again, however increased the initial clearance between the wheel and the rail to 0.005m, I think the MODI_MAILLAGE command now exectutes ok, I now have error in the STAT_NON_LINE calculation.

Below is an extract from the .mess file, for some reason I am unable to upload files from this computer... I will upload all the files from another later.

Thanks for the help!

Kev

Extract from new_case_test.mess
#

---

# Commande No : 0011 Concept de type : evol_noli
#

---

resu = STAT_NON_LINE(CHAM_MATER=Appl_mat,
MODELE=model,
CONTACT=cont,
EXCIT=(_F(CHARGE=boun_con,
TYPE_CHARGE='FIXE_CSTE',),
_F(CHARGE=load,
TYPE_CHARGE='FIXE_CSTE',),
),
COMP_INCR=_F(RELATION='ELAS',
ITER_INTE_PAS=0,
RESI_CPLAN_RELA=9.9999999999999995E-07,
PARM_THETA=1.0,
ITER_CPLAN_MAXI=1,
RESI_INTE_RELA=9.9999999999999995E-07,
DEFORMATION='PETIT',
PARM_ALPHA=1.0,
ITER_INTE_MAXI=20,),
INCREMENT=_F(LIST_INST=list,
PRECISION=9.9999999999999995E-07,),
NEWTON=_F(REAC_INCR=1,
MATRICE='TANGENTE',
REAC_ITER=0,
REAC_ITER_ELAS=0,),
CONVERGENCE=_F(ARRET='OUI',
ITER_GLOB_MAXI=10,
ITER_GLOB_ELAS=25,),
METHODE='NEWTON',
ARCHIVAGE=_F(PRECISION=9.9999999999999995E-07,
CRITERE='RELATIF',),
SOLVEUR=_F(RENUM='METIS',
STOP_SINGULIER='OUI',
NPREC=8,
METHODE='MULT_FRONT',
SYME='NON',),
)

La mémoire consommée actuellement hors JEVEUX est de 260.79 Mo.
La limite de l'allocation dynamique JEVEUX est fixée à 501.21 Mo.

Mot-clef : COMP_INCR - occurrence 1 : comportement ELAS - nombre de variables internes : 1
Noms des variables internes :
V1 : VIDE

Le système linéaire à résoudre contient 5194 noeuds dont:
- 5194 noeuds portant des degrés de liberté physiques
- 0 noeuds portant des degrés de liberté de Lagrange
Pour un total de 15582 équations.

La matrice est de taille 15582 équations.
Elle contient 276162 termes non nuls si elle est symétrique et 536742 termes non nuls si elle n'est pas symétrique (le nombre de
termes non nuls est susceptible de varier si l'on utilise le contact en formulation continue ou la méthode XFEM avec contact).
Soit un taux de remplissage de 0.221 %.

!

---

!
! <A> <CONTACT_1> !
! !
! Contact méthodes discrètes. !
! -> Les méthodes de contact discrètes supposent la symétrie de la matrice obtenue après assemblage. !
! Si votre modélisation produit une matrice non-symétrique, on force donc sa symétrie pour résoudre !
! le contact. !
! -> Risque & Conseil : !
! Ce changement peut conduire à des difficultés de convergence dans le processus de Newton mais en !
! aucun cas il ne produit des résultats faux. !
! !
! Si la matrice de rigidité de votre structure est symétrique, vous pouvez ignorer ce qui précède. !
! Enfin, il est possible de supprimer l'affichage de cette alarme en renseignant SYME='OUI' !
! sous le mot-clé facteur SOLVEUR. !
! !
! !
! Ceci est une alarme. Si vous ne comprenez pas le sens de cette !
! alarme, vous pouvez obtenir des résultats inattendus ! !
!

---

!


Nombre total de noeuds esclaves pour le contact: 32

Il n'y a pas d'état initial défini. On prend un état initial nul.
Le champ <DEPL> est initialisé a zéro
Le champ <SIEF_ELGA> est initialisé a zéro
Le champ <VARI_ELGA> est initialisé a zéro

Attention, en cas d'erreur (contact, loi de comportement, pilotage, ...), le pas de temps
ne sera pas redécoupé.

Archivage de l'état initial

Archivage des champs
Champ stocké <DEPL> à l'instant 0.000000000000e+00 pour le numéro d'ordre 0
Champ stocké <SIEF_ELGA> à l'instant 0.000000000000e+00 pour le numéro d'ordre 0
Champ stocké <VARI_ELGA> à l'instant 0.000000000000e+00 pour le numéro d'ordre 0
Champ stocké <COMPORTEMENT> à l'instant 0.000000000000e+00 pour le numéro d'ordre 0
Champ stocké <VALE_CONT> à l'instant 0.000000000000e+00 pour le numéro d'ordre 0

---

Instant de calcul: 5.000000000000e-02 - Niveau de découpe: 1

---

| CONTACT | NEWTON | RESIDU | RESIDU | OPTION | CONTACT | CONTACT |
| BCL. GEOM. | ITERATION | RELATIF | ABSOLU | ASSEMBLAGE | DISCRET | CRITERE |
| ITERATION | | RESI_GLOB_RELA | RESI_GLOB_MAXI | | NB. ITER | VALEUR |

---

!

---

!
! <EXCEPTION> <FACTOR_10> !
! !
! Problème : la matrice n'est pas factorisable : !
! Lors de la factorisation de la matrice, on a rencontré un problème !
! (pivot nul ou presque nul) Ã la ligne 6954 !
! Le degré de liberté associé à cette ligne est de type : (A) !
! Le nombre de décimales "perdues" pour cette ligne est : 15. !
! Ce nombre de décimales perdues est lié au degré de singularité de la matrice. Plus il est grand, !
! plus le système est singulier. Quand il dépasse 8, on considère que l'on a perdu 50% de la précision !
! sur les nombres réels (qui ont 15 décimales environ). !
! !
! Les différents types du degré de liberté sont : !
! * (A) : Degré liberté physique associé au noeud : N2318 et à la composante : DZ. !
! * ( : Degré de liberté de Lagrange associé au blocage : (XXXX). !
! * (C) : Degré de liberté de Lagrange associé à une relation linéaire entre plusieurs degrés de liberté. !
! La liste des noeuds impliqués dans cette relation linéaire est imprimée ci-dessus. !
! * (D) : Degré de liberté d'un système "généralisé". !
! Nom du noeud : XXXX !
! Composante : XXXX !
! Information complémentaire (éventuellement) : XXXX !
! !
! Conventions : !
! Ce message peut être émis dans des situations différentes pour lesquelles on ne connaît !
! pas toutes les informations imprimées ci-dessus. !
! On adopte alors les conventions suivantes : !
! * Si le numéro de la ligne est -999 : !
! Soit la matrice est vraiment singulière et la factorisation n'a pu aller jusqu'au bout. !
! Soit on ne sait pas attribuer la singularité de la matrice à une ligne de la matrice. !
! * Si le nombre de décimales perdues est -999 : !
! On ne sait pas déterminer la perte de décimales sur la ligne incriminée, !
! !
! Risques et conseils : !
! * Si la ligne correspond a un degré de liberté physique (A), il s'agit probablement d'un mouvement !
! de corps rigide mal bloqué. !
! Vérifiez les conditions aux limites. !
! Si vous faites du contact, il ne faut pas que la structure ne "tienne" que par le contact. !
! Vérifiez également les caractéristiques matériaux (module d'Young, ...). !
! !
! * Si la ligne correspond a un degré de liberté de Lagrange ( ou (C), il s'agit sans doute d'une condition !
! limite redondante. !
! En particulier, il se peut que la relation linéaire surabondante provienne des conditions de contact. !
! Peut-être devriez vous exclure certains noeuds des conditions de contact !
! (mots clés SANS_NOEUD et SANS_GROUP_NO). !
! !
! * Si le solveur utilisé est LDLT ou MULT_FRONT, vous pouvez utiliser le solveur MUMPS !
! car celui-ci est le seul à pouvoir factoriser les matrices qui ne sont pas définies positives. !
! !
! * Il se peut aussi que ce phénomène soit tout à fait normal avec X-FEM si la fissure passe !
! très près d'un noeud. !
! Si le nombre de décimales perdues n'est pas trop grand (max 10 décimales), !
! vous pouvez relancer le calcul en augmentant le nombre de décimales perdues autorisé : !
! mot-clé NPREC du mot clé facteur SOLVEUR. !
! Sinon, contactez l'équipe de développement. !
!

---

!


Validation du concept 'resu'.

Files attached, again thanks for the help! if you do get a chance to try and run the simulation that would be great!!!

Kev