STAT_NON_LINE >>> INCREMENT=_F(LIST_INST... Problem

  • MGolbs
  • MGolbs's Avatar Topic Author
  • Offline
  • Platinum Member
  • Platinum Member
More
13 years 3 months ago #5198 by MGolbs
Hallo,

habe ein Tria3 DKT Modell und will nichtlinear rechnen. Habe als Basis:

www.caelinux.org/wiki/index.php/Contrib:KeesWouters/plasticity

versucht zu verwenden. Leider bekomme ich schon im Eficas unter
... INCREMENT=_F(LIST_INST=time,.. rote Markierungen angezeigt.

STAT_NON_LINE(MODELE=pmode,
CHAM_MATER=matprops,
EXCIT=(_F(CHARGE=LoadFix,),
_F(CHARGE=LoadPres,FONC_MULT=ramp,),),
COMP_INCR=_F(RELATION='VMIS_ISOT_TRAC',
DEFORMATION='SIMO_MIEHE',
TOUT='OUI',),
INCREMENT=_F(LIST_INST=time,
SUBD_PAS=10,
SUBD_PAS_MINI=0.005,),
NEWTON=_F(REAC_INCR=1,
MATRICE='TANGENTE',
REAC_ITER=1,),
CONVERGENCE=_F(ITER_GLOB_MAXI=20,),
ARCHIVAGE=_F(PAS_ARCH=1,),);


Woran könnte das liegen?



Gruß und Dank Markus<br /><br />Post edited by: MGolbs, at: 2011/01/17 14:44

Dem Überflüssigen nachlaufen, heißt das Wesentliche verpassen.
Jules Saliège

Please Log in or Create an account to join the conversation.

More
13 years 3 months ago #5199 by kwou
Hoi Markus

warscheinlich unterschiedliche Versionen von C-Aster
Ich werde mal nachgucken (oder Klaus macht das ;-))
Sehe zB www.code-aster.org/forum2/viewtopic.php?id=14769

Bis bald - Kees


Bin wieder da:
Bitte aendern oder zufuegen in der comm Datei:

[code:1]
....
deflist = DEFI_LIST_INST(DEFI_LIST=_F(METHODE ='AUTO',
LIST_INST = time,
PAS_MINI = 0.0005),)
....

Presul=STAT_NON_LINE(MODELE=pmode,
CHAM_MATER=matprops,
EXCIT=(_F(CHARGE=LoadFix,),
_F(CHARGE=LoadPres,FONC_MULT=ramp,),),
COMP_INCR=_F(RELATION='VMIS_ISOT_TRAC',
DEFORMATION='SIMO_MIEHE',
TOUT='OUI',),
INCREMENT=_F(LIST_INST= deflist,), #time,
# INCREMENT=_F(LIST_INST=time,
# SUBD_PAS=10,
# SUBD_PAS_MINI=0.0005,),
NEWTON=_F(REAC_INCR=1,
MATRICE='TANGENTE',
REAC_ITER=1,),
CONVERGENCE=_F(ITER_GLOB_MAXI=20,),
ARCHIVAGE=_F(PAS_ARCH=1,),);
[/code:1]

Ich hoffe das funktioniert dan alles wieder. Ich werde es nachgucken und in der Wiki aendern.

Siehe auch:
www.caelinux.org/wiki/index.php/Contrib:...plasticity/solidbeam

Auf wiederhoeren - kees<br /><br />Post edited by: Kees Wouters, at: 2011/01/17 21:10

Interest: structural mechanics, solar energy (picture at 'my location' shows too little pv panels)

--
kind regards - kees

Please Log in or Create an account to join the conversation.

  • MGolbs
  • MGolbs's Avatar Topic Author
  • Offline
  • Platinum Member
  • Platinum Member
More
13 years 3 months ago #5203 by MGolbs
Hallo,

danke Kees. Werde den *.comm Code einfügen und noch mal den Solverlauf starten.

Gruß und Dank Markus

Hallo,

Solver bringt weiter Fehlermeldung, tippe aber darauf, dass auch meine Definitionen noch nicht stimmen.

!


!
! &lt;A&gt; &lt;MODELISA4_8&gt; !
! !
! !
! -&gt; Phase de vérification du maillage : présence de mailles doubles (ou triples, ...) !
! Les mailles multiples sont des mailles de noms différents qui ont la meme connectivité !
! (elles s'appuient sur les memes noeuds). !
! !
! -&gt; Risque &amp; Conseil : !
! Le risque est de modéliser 2 fois (ou plus) l'espace. On peut par exemple avoir !
! un modèle 2 fois trop lourd ou 2 fois trop rigide. !
! Remarque : les mailles concernées sont imprimées dans le fichier &quot;message&quot;. !
! Sur ce maillage, il est imprudent d'affecter des quantités avec le mot clé TOUT='OUI'. !
! !
! !
! !
! Ceci est une alarme. Si vous ne comprenez pas le sens de cette !
! alarme, vous pouvez obtenir des résultats inattendus ! !
!
!


!
!
! &lt;A&gt; &lt;DISCRETISATION_4&gt; !
! !
! !
! Avec une gestion automatique de la liste d'instants, la notion de NIVEAU de subdivisions est ignorée. !
! !
! !
! !
! Ceci est une alarme. Si vous ne comprenez pas le sens de cette !
! alarme, vous pouvez obtenir des résultats inattendus ! !
!
!


!
!
! &lt;A&gt; &lt;COMPOR1_47&gt; !
! !
! !
! La modelisation choisie &lt;C_PLAN&gt; est incompatible avec le comportement &lt;VMIS_ISOT_TRAC&gt;. Pour modeliser !
! des contraintes planes (ou des coques) avec ce comportement, on utilise ALGO_C_PLAN='DEBORST' !
! !
! !
! !
! Ceci est une alarme. Si vous ne comprenez pas le sens de cette !
! alarme, vous pouvez obtenir des résultats inattendus ! !
!
!


!
!
! &lt;EXCEPTION&gt; &lt;ALGORITH7_9&gt; !
! !
! !
! C_PLAN méthode DEBORST et SIMO_MIEHE incompatibles !
! !
! !
! !
! !
!
!


!
!
! 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 : !
! COMPOR1_47 émise 1 fois !
! DISCRETISATION_4 émise 1 fois !
! MODELISA4_8 émise 1 fois !
!
!


!
!
! &lt;A&gt; &lt;UTILITAI4_35&gt; !
! !
! !
! pas de numero d'ordre pour Solution !
! !
! !
! !
! Ceci est une alarme. Si vous ne comprenez pas le sens de cette !
! alarme, vous pouvez obtenir des résultats inattendus ! !
!
!

&lt;I&gt; &lt;FIN&gt; FERMETURE DE LA BASE &quot;GLOBALE&quot; EFFECTUEE.
&lt;I&gt; &lt;FIN&gt; ARRET NORMAL DANS &quot;FIN&quot; PAR APPEL A &quot;JEFINI&quot;.

********************************************************************************
* COMMAND : USER : SYSTEM : USER+SYS : ELAPSED *
********************************************************************************
* init (jdc) : 0.78 : 1.35 : 2.13 : 2.15 *
* . compile : 0.00 : 0.00 : 0.00 : 0.00 *
* . exec_compile : 0.18 : 0.16 : 0.34 : 0.35 *
* . report : 0.02 : 0.00 : 0.02 : 0.02 *
* . build : 0.00 : 0.00 : 0.00 : 0.00 *
* DEBUT : 0.02 : 0.49 : 0.51 : 0.53 *
* LIRE_MAILLAGE : 0.24 : 0.06 : 0.30 : 0.31 *
* MODI_MAILLAGE : 10.17 : 0.20 : 10.37 : 10.87 *
* DEFI_FONCTION : 0.02 : 0.00 : 0.02 : 0.01 *
* DEFI_MATERIAU : 0.00 : 0.00 : 0.00 : 0.01 *
* AFFE_MATERIAU : 0.02 : 0.00 : 0.02 : 0.01 *
* AFFE_MODELE : 0.05 : 0.00 : 0.05 : 0.33 *
* AFFE_CARA_ELEM : 0.29 : 0.01 : 0.30 : 0.30 *
* AFFE_CHAR_MECA : 0.09 : 0.01 : 0.10 : 0.10 *
* DEFI_FONCTION : 0.00 : 0.00 : 0.00 : 0.00 *
* AFFE_CHAR_MECA : 11.85 : 0.03 : 11.88 : 12.12 *
* DEFI_LIST_REEL : 0.00 : 0.00 : 0.00 : 0.00 *
* DEFI_LIST_INST : 0.02 : 0.00 : 0.02 : 0.02 *
* STAT_NON_LINE : 6.19 : 6.18 : 12.37 : 13.36 *
* FIN : 0.36 : 0.61 : 0.97 : 1.69 *
* . part Superviseur : 0.81 : 1.84 : 2.65 : 2.69 *
* . part Fortran : 28.94 : 6.49 : 35.43 : 37.44 *
********************************************************************************
* TOTAL_JOB : 29.75 : 8.33 : 38.08 : 40.13 *
********************************************************************************


Content of ERROR file

-- CODE_ASTER -- VERSION DE DEVELOPPEMENT STABILISEE 10.01 --

VERSION PARALLELE COMPILEE AVEC MPI

COPYRIGHT EDF-R&amp;D 1991 - 2009

EXECUTION DU : MA-18-JANV-2011 11:31:41

PLATE-FORME : HAMM-CAE-LINUX

NB MAX PROC OpenMP: 1

SYSTEME : LINUX

CPU : X86_64


!
!
! &lt;EXCEPTION&gt; &lt;ALGORITH7_9&gt; !
! !
! !
! C_PLAN méthode DEBORST et SIMO_MIEHE incompatibles !
! !
! !
! !
! !
!
!

&lt;S&gt; ERREUR UTILISATEUR RECUPEREE PAR LE SUPERVISEUR
&lt;I&gt; &lt;FIN&gt; FERMETURE DE LA BASE &quot;GLOBALE&quot; EFFECTUEE.
&lt;I&gt; &lt;FIN&gt; ARRET NORMAL DANS &quot;FIN&quot; PAR APPEL A &quot;JEFINI&quot;.
&lt;I&gt; FERMETURE DES BASES EFFECTUEE


cpu system cpu+sys elapsed
Preparation of environment 0.04 0.18 0.22 0.28
Copying datas 0.03 0.03 0.06 0.09
Code_Aster run 60.27 18.25 78.52 42.64
Code_Aster run - rank #0 29.75 8.33 38.08 0.00
Copying results 0.09 0.56 0.65 0.62
Total 60.53 19.24 79.77 43.90
(*) cpu and system times are not correctly counted using mpirun.

as_run 1.8.0

--- DIAGNOSTIC JOB : &lt;F&gt;_COPY_ERROR


EXIT_CODE=4



Hier noch meine *.comm

[code:1]
DEBUT();

MAIL=LIRE_MAILLAGE(UNITE=20,
FORMAT='MED',
NOM_MED='ausgleic',
INFO=1,);

MAIL=MODI_MAILLAGE(reuse =MAIL,
MAILLAGE=MAIL,
ORIE_NORM_COQUE=_F(GROUP_MA='EleTria3',),);

Traction=DEFI_FONCTION(NOM_PARA='EPSI',VALE=(0.003143,220.0,
0.005199,290.0,
0.010668,330.0,
0.022036,365.0,
0.041728,405.0,
0.072294,440.0,
0.116391,470.0,
0.176765,510.0,
0.256243,545.0,
0.357725,585.0,
0.484174,620.0,
0.638615,650.0,
0.824129,660.0,
),INTERPOL='LIN',PROL_DROITE='LINEAIRE',PROL_GAUCHE='EXCLU',);

Stahl=DEFI_MATERIAU(ELAS=_F(E=210000.0,
NU=0.3,),
TRACTION=_F(SIGM=Traction,),);

MatMod=AFFE_MATERIAU(MAILLAGE=MAIL,
AFFE=_F(TOUT='OUI',
MATER=Stahl,),);

EleMod=AFFE_MODELE(MAILLAGE=MAIL,
INFO=2,
AFFE=_F(GROUP_MA='EleTria3',
PHENOMENE='MECANIQUE',
MODELISATION='DKT',),);

EleProp=AFFE_CARA_ELEM(MODELE=EleMod,
COQUE=_F(GROUP_MA='EleTria3',
EPAIS=5.0,
VECTEUR=(12.0,19.0,27.0,),),);

Ground=AFFE_CHAR_MECA(VERI_NORM='OUI',
MODELE=EleMod,
DDL_IMPO=_F(GROUP_MA='LagerEle',
DX=0.0,
DY=0.0,
DZ=0.0,
DRX=0.0,
DRY=0.0,
DRZ=0.0,),);

tsteps = 0.1;

ramp=DEFI_FONCTION(NOM_PARA='INST',VALE=(0.00,0.00,
1.00,1.00,
),INFO=2,TITRE='ramp',);

Pilot=AFFE_CHAR_MECA(VERI_NORM='OUI',
MODELE=EleMod,
PRES_REP=_F(GROUP_MA='EleTria3',
PRES=-10.0,),);

Time=DEFI_LIST_REEL(DEBUT=0,
INTERVALLE=_F(JUSQU_A=1,
NOMBRE=1,),);

deflist = DEFI_LIST_INST(DEFI_LIST=_F(METHODE ='AUTO',
LIST_INST = Time,
PAS_MINI = 0.0005),)

Solution=STAT_NON_LINE(MODELE=EleMod,
CHAM_MATER=MatMod,
CARA_ELEM=EleProp,
EXCIT=(_F(CHARGE=Ground,),
_F(CHARGE=Pilot,FONC_MULT=ramp,),),
COMP_INCR=_F(RELATION='VMIS_ISOT_TRAC',
DEFORMATION='SIMO_MIEHE',
TOUT='OUI',),
INCREMENT=_F(LIST_INST= deflist,), #Time,
# INCREMENT=_F(LIST_INST=Time,
# SUBD_PAS=10,
# SUBD_PAS_MINI=0.0005,),
NEWTON=_F(REAC_INCR=1,
MATRICE='TANGENTE',
REAC_ITER=1,),
CONVERGENCE=_F(ITER_GLOB_MAXI=20,),
ARCHIVAGE=_F(PAS_ARCH=1,),);

Solution=CALC_ELEM(reuse =Solution,
MODELE=EleMod,
RESULTAT=Solution,
OPTION=('SIEF_ELNO_ELGA','EPSG_ELNO_DEPL','EPSI_ELNO_DEPL',),);

IMPR_RESU(MODELE=EleMod,
FORMAT='MED',
RESU=_F(MAILLAGE=MAIL,
RESULTAT=Solution,
INFO_MAILLAGE='OUI',
PARTIE='REEL',
TOUT_CHAM='OUI',
TOUT_CMP='OUI',
TOUT='OUI',),
INFO=2,);

FIN(FORMAT_HDF='OUI',);
[/code:1]

Gruß und Dank Markus<br /><br />Post edited by: MGolbs, at: 2011/01/18 12:35

Dem Überflüssigen nachlaufen, heißt das Wesentliche verpassen.
Jules Saliège

Please Log in or Create an account to join the conversation.

Moderators: catux
Time to create page: 0.126 seconds
Powered by Kunena Forum