2D contact problem

I'm new using aster (Salome-meca2009) and I'm trying to solve a simple 2D contact problem, I try with:
The JML tutorial (Nº5)
The press tutorial
The advanced contact example
PinAssembly

the problem is a square with half a circle on top, the bottom of the square is fixed, as in JML tutorial



Until now I been able to solve many problems but in all tutorials I always get stuck with:

Arret : absence de convergence avec le nombre d'itérations requis.

The important parts of the .mess file are (or at least I think is important):


[code:1]
--- NOMBRE TOTAL DE NOEUDS : 19079 DONT : 648 NOEUDS "LAGRANGE"
--- NOMBRE TOTAL D'EQUATIONS : 37510
--- TAILLE DU PROFIL MORSE DE LA TRIANGULAIRE SUPERIEURE (FORMAT SCR): 440625
--- DONC LA TAILLE DE LA MATRICE EST:
--- EN SYMETRIQUE NNZ= 440625
--- EN NON SYMETRIQUE NNZ= 843740

--- NOMBRE TOTAL DE NOEUDS ESCLAVES : 127

!

---

!
! <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 ! !
!

---

!
[/code:1]



[code:1]
SUBDIVISION DU PAS DE TEMPS EN 20 SOUS PAS

Méthode de subdivision : Pendant la nouvelle subdivision
Le pas minimal de la subdivision est atteint.
Pas de Temps actuel = <7.812500E-11>
Pas de Temps minimum imposé = <1.000000E-10>
Niveau de subdivision = <6>
Méthode de subdivision = <SIMPLE>
Conseil :
Diminuer SUBD_PAS_MINI. Il est également possible d'ajuster SUBD_NIVEAU pour
indiquer le nombre successif de subdivision d'un pas de temps.
Si les 2 mots clefs SUBD_NIVEAU et SUBD_PAS_MINI sont utilisés la subdivision
s'arrete dès que l'un des 2 critères est vérifié.

!

---

!
! <EXCEPTION> <MECANONLINE_83> !
! !
! !
! Arret : absence de convergence avec le nombre d'itérations requis. !
! !
! !
! !
! !
!

---

!


Validation du concept 'solution'.

!

---

!
! Liste des alarmes émises lors de l'exécution du calcul. !
! !
! Les alarmes que vous avez choisies d'ignorer sont précédées de (*). !
! Nombre d'occurrences pour chacune des alarmes : !
! CONTACT_1 émise 1 fois !
!

---

!
[/code:1]

[code:1]
JDC.py : ERREUR A L'EXECUTION - INTERRUPTION
>> JDC.py : DEBUT RAPPORT
CR d'execution de JDC en MIXTE
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! <S> Exception utilisateur levee mais pas interceptee. !
! Les bases sont fermees. !
! NonConvergenceError : <EXCEPTION LEVEE> MECANONLINE_83 !
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
fin CR d'execution de JDC en MIXTE

>> JDC.py : FIN RAPPORT
EXECUTION_CODE_ASTER_EXIT_4081=1

<F>_ABNORMAL_ABORT exit code = 1
[/code:1]

I try to change the values for SUBD_NIVEAU and SUBD_PAS_MINI, but I think that this is not the problem (1e-10 and 100 should be more than enought I think)

I would aprecciate it if any body can help I been stuck with this for some days and although I have learn somethings I run out of patience.

I attach the files if somebody have some time to look at them, the template for the comm is from The advanced contact example with some modifications

NOTE: the mesh is very fine because I try everything

Post edited by: Danilo Estay, at: 2010/03/16 20:56

EDIT: I have problems uploading the files...<br /><br />Post edited by: Danilo Estay, at: 2010/03/16 20:59

Here are the files attached