STAT_NON_LINE >>> INCREMENT=_F(LIST_INST... Problem
- MGolbs
-
Topic Author
- Offline
- Platinum Member
-
14 years 3 months ago #5198
by MGolbs
Dem Überflüssigen nachlaufen, heißt das Wesentliche verpassen.
Jules Saliège
STAT_NON_LINE >>> INCREMENT=_F(LIST_INST... Problem was created 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
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
- kwou
-
- Offline
- Moderator
-
14 years 3 months ago #5199
by kwou
Interest: structural mechanics, solar energy (picture at 'my location' shows too little pv panels)
--
kind regards - kees
Replied by kwou on topic Re:STAT_NON_LINE >>> INCREMENT=_F(LIST_INST... Problem
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
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
- MGolbs
-
Topic Author
- Offline
- Platinum Member
-
14 years 3 months ago #5203
by MGolbs
Dem Überflüssigen nachlaufen, heißt das Wesentliche verpassen.
Jules Saliège
Replied by MGolbs on topic Re:STAT_NON_LINE >>> INCREMENT=_F(LIST_INST... Problem
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.
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
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.
!
!
! <A> <MODELISA4_8> !
! !
! !
! -> 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). !
! !
! -> Risque & 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 "message". !
! 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 ! !
!
!
!
!
! <A> <DISCRETISATION_4> !
! !
! !
! 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 ! !
!
!
!
!
! <A> <COMPOR1_47> !
! !
! !
! La modelisation choisie <C_PLAN> est incompatible avec le comportement <VMIS_ISOT_TRAC>. 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 ! !
!
!
!
!
! <EXCEPTION> <ALGORITH7_9> !
! !
! !
! 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 !
!
!
!
!
! <A> <UTILITAI4_35> !
! !
! !
! 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 ! !
!
!
<I> <FIN> FERMETURE DE LA BASE "GLOBALE" EFFECTUEE.
<I> <FIN> ARRET NORMAL DANS "FIN" PAR APPEL A "JEFINI".
********************************************************************************
* 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&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
!
!
! <EXCEPTION> <ALGORITH7_9> !
! !
! !
! C_PLAN méthode DEBORST et SIMO_MIEHE incompatibles !
! !
! !
! !
! !
!
!
<S> ERREUR UTILISATEUR RECUPEREE PAR LE SUPERVISEUR
<I> <FIN> FERMETURE DE LA BASE "GLOBALE" EFFECTUEE.
<I> <FIN> ARRET NORMAL DANS "FIN" PAR APPEL A "JEFINI".
<I> 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 : <F>_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
Moderators: catux
Time to create page: 0.136 seconds