error using pilotage
- alessandro
- Topic Author
- Offline
- New Member
-
Less
More
- Posts: 10
- Thank you received: 0
8 years 11 months ago - 8 years 11 months ago #8098
by alessandro
error using pilotage was created by alessandro
Hi, I am using pilotage to follow the postbuckling behaviour of a shell.
The simulation does not start and I have the following error that I am not able to solve
<S> Exception utilisateur levee mais pas interceptee. !
! Les bases sont fermees. !
! Type de l'exception : error !
! !
! Pilotage. !
! La composante <DX> n'a pas été trouvée dans la numérotation. !
! Vérifier NOM_CMP dans le mot-clef PILOTAGE.
Any help would be very appreciated.
Best regards.
Here is my .mess file
The simulation does not start and I have the following error that I am not able to solve
<S> Exception utilisateur levee mais pas interceptee. !
! Les bases sont fermees. !
! Type de l'exception : error !
! !
! Pilotage. !
! La composante <DX> n'a pas été trouvée dans la numérotation. !
! Vérifier NOM_CMP dans le mot-clef PILOTAGE.
Any help would be very appreciated.
Best regards.
Here is my .mess file
# ------------------------------------------------------------------------------------------
# Impression du contenu du fichier de commandes �� ex��cuter :
# ------------------------------------------------------------------------------------------
#---------------------------------------------------------------------------------------
#
# UNITS: MASS LENTH TIME FORCE STRESS ENERGY
# Kg m s N Pa N-m
#---------------------------------------------------------------------------------------
######## Model Parameters ########
th=0.002;
E=7e10; #[Pa]
nu=0.3;
rho=2700.0; #[Kg/m^3]
imp_disp=-0.005; #[m]
P_int=12000; #[Pa]
P_ext=0.; #[Pa]
presss=P_int;
######## Simulation Parameters ########
stopTime=1 #[s] Time to interrupt the simulation
tf_0=0.1*stopTime #[s] Stop Time for first Calculation Run
tf_1=0.8*stopTime #[s] Stop Time for first Calculation Run
tf_2=1*stopTime #[s] Stop Time for first Calculation Run
maxdT = 1e-2*stopTime #[s] Max Time Step
archivT1= 0.1*stopTime #[s] Time increment to write results
precis=1e-6
tries=30
DEBUT();
#----------------------------------------------------
# (1) Reading of the mesh
#----------------------------------------------------
meshI=LIRE_MAILLAGE ( UNITE=20 ,FORMAT='MED' , INFO_MED=1,) ;
Qmesh=CREA_MAILLAGE(MAILLAGE=meshI,
LINE_QUAD=_F(TOUT='OUI',PREF_NOEUD='NNS'),);
mesh=CREA_MAILLAGE(MAILLAGE=Qmesh,
MODI_MAILLE=(
_F(TOUT='OUI',OPTION='QUAD8_9',),
),
);
#mesh=DEFI_GROUP (reuse =mesh,MAILLAGE=mesh ,CREA_GROUP_MA=( _F ( NOM='TOUT' , TOUT='OUI' , ) , ), CREA_GROUP_NO=_F(TOUT_GROUP_MA='OUI' ,) ,);
mesh=DEFI_GROUP (reuse =mesh,MAILLAGE=mesh ,
CREA_GROUP_NO=_F(TOUT_GROUP_MA='OUI' ,) ,);
mesh=MODI_MAILLAGE (
reuse =mesh ,
MAILLAGE=mesh ,
ORIE_NORM_COQUE=(_F (GROUP_MA=('surface',),),),
);
#---------------------------------------------------
# MODELISATION MECHANIC
# Element types
#---------------------------------------------------
Model=AFFE_MODELE(
MAILLAGE=mesh,
AFFE=(_F(
TOUT='OUI',
PHENOMENE='MECANIQUE',
MODELISATION='COQUE_3D', #COQUE_3D is for curved shell
),
),
);
#---------------------------------------------------
# MATERIAL AND MATERIAL FIELD
#---------------------------------------------------
Alu=DEFI_MATERIAU(ELAS=_F(E=E,NU=nu,RHO=rho,),);
Mat=AFFE_MATERIAU(MAILLAGE=mesh,AFFE=_F(
TOUT='OUI',
MATER=Alu,),);
#---------------------------------------------------
# ELEMENT PROPERTIES
#---------------------------------------------------
CARA=AFFE_CARA_ELEM(MODELE=Model,
COQUE=(_F(
GROUP_MA='surface',
EPAIS=th, #thickness of the shell
COQUE_NCOU=5, #number of nodes on the thickness
COEF_RIGI_DRZ =0.001,
),
),
);
#---------------------------------------------------
# MECHANICAL LOADS AND BOUNDARY CONDITIONS
#---------------------------------------------------
# fix support
BCnd=AFFE_CHAR_MECA (MODELE=Model ,
DDL_IMPO =(_F( GROUP_NO='equator',DX=0.0,DY=0.0,DZ=0.0, ),),
);
# indentation
displac=AFFE_CHAR_MECA(MODELE=Model ,DDL_IMPO =_F(GROUP_NO='top',DY=imp_disp ) ,);
# pressure
#appP=AFFE_CHAR_MECA (MODELE=Model ,PRES_REP =_F(GROUP_MA='in',PRES=presss, ) ,);
#---------------------------------------------------
# DEFINE CONTACT
#---------------------------------------------------
#cnt=DEFI_CONTACT(MODELE=Model,FORMULATION='CONTINUE',ALGO_RESO_GEOM='NEWTON',ZONE=(_F(GROUP_MA_MAIT='MASTER1',GROUP_MA_ESCL='SLAVE1',),),);
#---------------------------------------------------
# TIME CONTROL
#---------------------------------------------------
# List of time instant for the INCREMENTAL SOLUTION
L_INST1=DEFI_LIST_REEL(DEBUT=0.,INTERVALLE=_F(JUSQU_A=stopTime,NOMBRE=round(stopTime/maxdT,),),INFO=2,);
# List of time instant for ARCHIVAGE
saveT1=DEFI_LIST_REEL(DEBUT=0.,INTERVALLE=_F(JUSQU_A=stopTime,NOMBRE=round(stopTime/maxdT,),),);
# Variable Time Step
DEFLIST1=DEFI_LIST_INST(
DEFI_LIST=_F(METHODE='AUTO',LIST_INST=L_INST1,PAS_MINI=1e-10,),
ECHEC=_F(EVENEMENT='ERREUR',SUBD_PAS=4,),
ADAPTATION=_F(EVENEMENT='SEUIL',NB_INCR_SEUIL=3,NOM_PARA='NB_ITER_NEWTON',CRIT_COMP='LE',PCENT_AUGM=200.0,),
);
RampeP1=DEFI_FONCTION(NOM_PARA='INST',VALE=(0.0,0.0,stopTime,1.0,),INTERPOL='LIN',PROL_DROITE='CONSTANT',PROL_GAUCHE='CONSTANT',);
#---------------------------------------------------
# MECHANICAL SOLUTION, nonlinear
#---------------------------------------------------
result=STAT_NON_LINE(MODELE=Model ,
CHAM_MATER=Mat,
CARA_ELEM=CARA,
EXCIT=(
_F(CHARGE=BCnd,),
_F(#CHARGE=appP ,TYPE_CHARGE ='SUIV',
CHARGE=displac,
#FONC_MULT=RampeP1,
TYPE_CHARGE = 'FIXE_PILO'),
),
COMP_ELAS=_F( RELATION = 'ELAS',
DEFORMATION = 'GROT_GDEP',
),
INCREMENT=_F(LIST_INST=DEFLIST1,INST_FIN=tf_0,),
PILOTAGE=_F( TYPE = 'LONG_ARC',
SELECTION='ANGL_INCR_DEPL',
GROUP_NO = 'surface',
NOM_CMP = ('DX', 'DY', 'DZ',),
ETA_PILO_MAX=1.,
COEF_MULT = 120.),
NEWTON=_F( REAC_ITER = 1,MATRICE = 'TANGENTE'),
CONVERGENCE=_F(ITER_GLOB_MAXI = 30,
RESI_GLOB_RELA = 1.E-6),
)
# calc more fields
result=CALC_CHAMP(reuse =result,
MODELE=Model,
RESULTAT=result,
#CONTRAINTE='SIGM_NOEU',
#DEFORMATION='EPSI_NOEU',
FORCE=('FORC_NODA','REAC_NODA'),
#VARI_INTERNE=('VARI_NOEU',),
TOUT='OUI',
);
#----------------------------------------------------
# (8) Printing
#----------------------------------------------------
Tdisp=POST_RELEVE_T(
ACTION=(_F (
OPERATION='EXTRACTION' ,
INTITULE='displ_middle' ,
RESULTAT=result,
TOUT='OUI' ,
NOM_CHAM='DEPL' ,
TOUT_ORDRE='OUI' ,
GROUP_NO='top' ,
RESULTANTE = ( 'DZ' , ) ,
),
),
);
Tforce=POST_RELEVE_T(
ACTION=_F(RESULTAT=result,
GROUP_NO='top',
OPERATION='EXTRACTION',
NOM_CHAM='REAC_NODA',
RESULTANTE=('DX','DY','DZ'),
MOYE_NOEUD='OUI',
INTITULE='',
),
);
D=RECU_FONCTION(TABLE=Tdisp,
PARA_X='INST',
PARA_Y='DZ',
);
F=RECU_FONCTION(TABLE=Tforce,
PARA_X='INST',
PARA_Y='DZ',
);
IMPR_FONCTION(
FORMAT='TABLEAU',
UNITE=81,
COURBE=_F(FONC_X=D,
FONC_Y=F,),
COMMENTAIRE = '%',
COMM_PARA = '%',
);
asseX_0=POST_RELEVE_T(
ACTION=(_F (
GROUP_NO='asseX' ,
INTITULE='asseX' ,
RESULTAT=result,
NOM_CHAM='DEPL' ,
NOM_CMP=('DX' , 'DY' , 'DZ' , ) ,
TOUT_ORDRE='OUI' , #all time steps
OPERATION='EXTRACTION' ,
),
),
TITRE='asseX' ,
);
IMPR_TABLE (FORMAT='TABLEAU',
COMMENTAIRE = '%',
COMM_PARA = '%',
UNITE=82,
TABLE=asseX_0 ,);
asseZ_0=POST_RELEVE_T(
ACTION=(_F (
GROUP_NO='asseZ' ,
INTITULE='asseZ' ,
RESULTAT=result,
NOM_CHAM='DEPL' ,
NOM_CMP=('DX' , 'DY' , 'DZ' , ) ,
TOUT_ORDRE='OUI' , #all time steps
OPERATION='EXTRACTION' ,
),
),
TITRE='asseZ' ,
);
IMPR_TABLE (FORMAT='TABLEAU',
COMMENTAIRE = '%',
COMM_PARA = '%',
UNITE=83,
TABLE=asseZ_0 ,);
#resultI=PROJ_CHAMP(RESULTAT=result,MODELE_1=Model, # project fields of this model to MODELE_2=ModelI,); # field of model_2
#IMPR_RESU(FORMAT='MED',UNITE=80,RESU=_F(MAILLAGE=meshI,RESULTAT=resultI,NOM_CHAM=('DEPL'),),);
FIN();
# ------------------------------------------------------------------------------------------
-- CODE_ASTER -- VERSION : EXPLOITATION (stable) --
Version 11.7.0 modifi��e le 10/12/2014
r��vision 7a39a0fb6aea - branche 'v11'
Copyright EDF R&D 1991 - 2015
Ex��cution du : Wed May 20 06:40:14 2015
Nom de la machine : naxos
Architecture : 64bit
Type de processeur : x86_64
Syst��me d'exploitation : Linux CentOS 6.4 Final 2.6.32-358.23.2.el6.x86_64
Langue des messages : en (UTF-8)
Version de Python : 2.7.3
Version de NumPy : 1.8.1
Parall��lisme MPI : inactif
Parall��lisme OpenMP : actif
Nombre de processus utilis��s : 1
Version de la librairie HDF5 : 1.8.10
Version de la librairie MED : 3.0.8
Version de la librairie MUMPS : 4.10.0
Version de la librairie SCOTCH : 5.1.10
M��moire limite pour l'ex��cution : 15120.00 Mo
consomm��e par l'initialisation : 307.90 Mo
par les objets du jeu de commandes : 2.21 Mo
reste pour l'allocation dynamique : 14809.37 Mo
Taille limite des fichiers d'��change : 48.00 Go
# ------------------------------------------------------------------------------------------
# Commande No : 0001 Concept de type : -
# ------------------------------------------------------------------------------------------
DEBUT(INFO=1,
PAR_LOT='OUI',
IMPR_MACRO='NON',
MEMOIRE=_F(TAILLE_BLOC=800.0,
TAILLE_GROUP_ELEM=1000,),
RESERVE_CPU=_F(BORNE=900,),
MESURE_TEMPS=_F(MOYENNE='OUI',
NIVE_DETAIL=1,),
)
Valeur initiale du temps CPU maximum = 108900 secondes
Valeur du temps CPU maximum pass�� aux commandes = 108000 secondes
R��serve CPU pr��vue = 900 secondes
D��but de lecture...
R��ouverture de la base
Nom de la base : ELEMBASE
Cr����e avec la version : 11.07.00
Nombre d'enregistrements utilis��s : 24
Nombre d'enregistrements maximum : 512
Longueur d'enregistrement (octets) : 819200
Nombre d'identificateurs utilis��s : 123
Taille maximum du r��pertoire : 300
Pourcentage d'utilisation du r��pertoire : 41 %
Fermeture de la base
Nom de la base : ELEMBASE
Nombre d'enregistrements utilis��s : 24
Nombre d'enregistrements maximum : 512
Longueur d'enregistrement (octets) : 819200
Nombre total d'acc��s en lecture : 49
Volume des acc��s en lecture : 38.28 Mo.
Nombre total d'acc��s en ��criture : 3
Volume des acc��s en ��criture : 2.34 Mo.
Nombre d'identificateurs utilis��s : 123
Taille maximum du r��pertoire : 300
Pourcentage d'utilisation du r��pertoire : 41 %
Relecture des catalogues des ��l��ments faite.
Fin de lecture (dur��e 0.000000 s.)
# M��moire (Mo) : 338.41 / 333.93 / 27.33 / 14.95 (VmPeak / VmSize / Optimum / Minimum)
# Fin commande No : 0001 user+syst: 0.06s (syst: 0.04s, elaps: 0.05s)
# ------------------------------------------------------------------------------------------
# ------------------------------------------------------------------------------------------
# Commande No : 0002 Concept de type : maillage_sdaster
# ------------------------------------------------------------------------------------------
meshI = LIRE_MAILLAGE(FORMAT='MED',
UNITE=20,
INFO_MED=1,
INFO=1,
VERI_MAIL=_F(APLAT=1.E-3,
VERIF='OUI',),
)
====== VERIFICATION DU MAILLAGE ======
MAILLE POI1 M1 INCLUSE DANS UNE AUTRE
------------ MAILLAGE meshI - IMPRESSIONS NIVEAU 1 ------------
MED file generated by Gmsh
NOMBRE DE NOEUDS 765
NOMBRE DE MAILLES 861
POI1 1
SEG2 128
QUAD4 732
NOMBRE DE GROUPES DE MAILLES 5
top 1
asseX 32
equator 64
asseZ 32
surface 732
--------------------------------------------------------------------------------
# M��moire (Mo) : 338.41 / 334.32 / 27.33 / 14.95 (VmPeak / VmSize / Optimum / Minimum)
# Fin commande No : 0002 user+syst: 0.01s (syst: 0.00s, elaps: 0.01s)
# ------------------------------------------------------------------------------------------
# ------------------------------------------------------------------------------------------
# Commande No : 0003 Concept de type : maillage_sdaster
# ------------------------------------------------------------------------------------------
Qmesh = CREA_MAILLAGE(MAILLAGE=meshI,
LINE_QUAD=_F(PREF_NOEUD='NNS',
TOUT='OUI',
PREF_NUME=1,),
INFO=1,)
MOT CLE FACTEUR "LINE_QUAD", OCCURRENCE 1
TRANSFORMATION DE 128 MAILLES SEG2 EN SEG3
TRANSFORMATION DE 732 MAILLES QUAD4 EN QUAD8
------------ MAILLAGE Qmesh - IMPRESSIONS NIVEAU 1 ------------
ASTER 11.07.00 CONCEPT Qmesh CALCULE LE 20/05/2015 A 06:40:14 DE TYPE
MAILLAGE_SDASTER
NOMBRE DE NOEUDS 2261
NOMBRE DE MAILLES 861
POI1 1
SEG3 128
QUAD8 732
NOMBRE DE GROUPES DE MAILLES 5
top 1
asseX 32
equator 64
asseZ 32
surface 732
--------------------------------------------------------------------------------
# M��moire (Mo) : 338.41 / 334.43 / 27.33 / 14.95 (VmPeak / VmSize / Optimum / Minimum)
# Fin commande No : 0003 user+syst: 0.01s (syst: 0.00s, elaps: 0.01s)
# ------------------------------------------------------------------------------------------
# ------------------------------------------------------------------------------------------
# Commande No : 0004 Concept de type : maillage_sdaster
# ------------------------------------------------------------------------------------------
mesh = CREA_MAILLAGE(MAILLAGE=Qmesh,
MODI_MAILLE=_F(OPTION='QUAD8_9',
TOUT='OUI',
PREF_NOEUD='NS',
PREF_NUME=1,),
INFO=1,)
MOT CLE FACTEUR "MODI_MAILLE", OCCURRENCE 1
MODIFICATION DE 732 MAILLES QUAD8 EN QUAD9
------------ MAILLAGE mesh - IMPRESSIONS NIVEAU 1 ------------
ASTER 11.07.00 CONCEPT mesh CALCULE LE 20/05/2015 A 06:40:14 DE TYPE
MAILLAGE_SDASTER
NOMBRE DE NOEUDS 2993
NOMBRE DE MAILLES 861
POI1 1
SEG3 128
QUAD9 732
NOMBRE DE GROUPES DE MAILLES 5
top 1
asseX 32
equator 64
asseZ 32
surface 732
--------------------------------------------------------------------------------
# M��moire (Mo) : 338.41 / 334.69 / 27.33 / 14.95 (VmPeak / VmSize / Optimum / Minimum)
# Fin commande No : 0004 user+syst: 0.02s (syst: 0.00s, elaps: 0.02s)
# ------------------------------------------------------------------------------------------
# ------------------------------------------------------------------------------------------
# Commande No : 0005 Concept de type : maillage_sdaster
# ------------------------------------------------------------------------------------------
mesh = DEFI_GROUP(reuse=mesh,
MAILLAGE=mesh,
CREA_GROUP_NO=_F(TOUT_GROUP_MA='OUI',),
ALARME='OUI',)
NOMBRE DE GROUPES DE NOEUDS CREES : 5
=======================================
------------------------------------------------------
! NOM DU GROUPE ! NBRE DE NOEUDS DU !
! NOEUDS ! GROUPE_NO !
------------------------------------------------------
! top ! 1 !
! asseX ! 65 !
! equator ! 128 !
! asseZ ! 65 !
! surface ! 2993 !
------------------------------------------------------
# M��moire (Mo) : 338.41 / 334.71 / 27.33 / 14.95 (VmPeak / VmSize / Optimum / Minimum)
# Fin commande No : 0005 user+syst: 0.00s (syst: 0.00s, elaps: 0.00s)
# ------------------------------------------------------------------------------------------
# ------------------------------------------------------------------------------------------
# Commande No : 0006 Concept de type : maillage_sdaster
# ------------------------------------------------------------------------------------------
mesh = MODI_MAILLAGE(reuse=mesh,
MAILLAGE=mesh,
ORIE_NORM_COQUE=_F(GROUP_MA=('surface', ),),
INFO=1,)
TRAITEMENT DU GROUP_MA: surface DE 732 MAILLES
0 MAILLE(S) ONT ETE ORIENTEE(S)
# M��moire (Mo) : 338.41 / 334.71 / 27.33 / 14.95 (VmPeak / VmSize / Optimum / Minimum)
# Fin commande No : 0006 user+syst: 0.01s (syst: 0.00s, elaps: 0.02s)
# ------------------------------------------------------------------------------------------
# ------------------------------------------------------------------------------------------
# Commande No : 0007 Concept de type : modele_sdaster
# ------------------------------------------------------------------------------------------
Model = AFFE_MODELE(MAILLAGE=mesh,
AFFE=_F(PHENOMENE='MECANIQUE',
TOUT='OUI',
MODELISATION=('COQUE_3D', ),),
INFO=1,
VERI_JACOBIEN='OUI',
PARTITION=_F(PARALLELISME='GROUP_ELEM',),
)
SUR LES 861 MAILLES DU MAILLAGE mesh
ON A DEMANDE L'AFFECTATION DE 861
ON A PU EN AFFECTER 860
MODELISATION ELEMENT FINI TYPE MAILLE NOMBRE
COQUE_3D MEBOCQ3 SEG3 128
COQUE_3D MEC3QU9H QUAD9 732
#2 Calculs elementaires et assemblages CPU (USER+SYST/SYST/ELAPS): 0.00 0.00 0.00
# M��moire (Mo) : 341.75 / 341.71 / 27.33 / 21.06 (VmPeak / VmSize / Optimum / Minimum)
# Fin commande No : 0007 user+syst: 0.02s (syst: 0.01s, elaps: 0.01s)
# ------------------------------------------------------------------------------------------
# ------------------------------------------------------------------------------------------
# Commande No : 0008 Concept de type : mater_sdaster
# ------------------------------------------------------------------------------------------
Alu = DEFI_MATERIAU(ELAS=_F(E=7.E10,
RHO=2700.0,
NU=0.3,),
INFO=1,)
MATERIAU : Alu
RELATION DE COMPORTEMENT: ELAS
# M��moire (Mo) : 341.75 / 341.71 / 27.33 / 21.06 (VmPeak / VmSize / Optimum / Minimum)
# Fin commande No : 0008 user+syst: 0.00s (syst: 0.00s, elaps: 0.00s)
# ------------------------------------------------------------------------------------------
# ------------------------------------------------------------------------------------------
# Commande No : 0009 Concept de type : cham_mater
# ------------------------------------------------------------------------------------------
Mat = AFFE_MATERIAU(MAILLAGE=mesh,
AFFE=_F(MATER=(Alu, ),
TOUT='OUI',),
INFO=1,
)
# M��moire (Mo) : 341.75 / 341.71 / 27.33 / 21.06 (VmPeak / VmSize / Optimum / Minimum)
# Fin commande No : 0009 user+syst: 0.00s (syst: 0.00s, elaps: 0.00s)
# ------------------------------------------------------------------------------------------
# ------------------------------------------------------------------------------------------
# Commande No : 0010 Concept de type : cara_elem
# ------------------------------------------------------------------------------------------
CARA = AFFE_CARA_ELEM(MODELE=Model,
COQUE=_F(EPAIS=2.E-3,
GROUP_MA=('surface', ),
COQUE_NCOU=5,
COEF_RIGI_DRZ=1.E-3,
A_CIS=0.8333333,
MODI_METRIQUE='NON',),
INFO=1,)
LE MODELE Model CONTIENT :
860 ELEMENT(S) COQUE
#2 Calculs elementaires et assemblages CPU (USER+SYST/SYST/ELAPS): 0.00 0.00 0.00
# M��moire (Mo) : 341.82 / 341.72 / 27.33 / 21.10 (VmPeak / VmSize / Optimum / Minimum)
# Fin commande No : 0010 user+syst: 0.01s (syst: 0.00s, elaps: 0.01s)
# ------------------------------------------------------------------------------------------
# ------------------------------------------------------------------------------------------
# Commande No : 0011 Concept de type : char_meca
# ------------------------------------------------------------------------------------------
BCnd = AFFE_CHAR_MECA(MODELE=Model,
DDL_IMPO=_F(GROUP_NO=('equator', ),
DZ=0.0,
DX=0.0,
DY=0.0,),
INFO=1,
VERI_NORM='OUI',
)
# M��moire (Mo) : 342.09 / 341.95 / 27.33 / 21.10 (VmPeak / VmSize / Optimum / Minimum)
# Fin commande No : 0011 user+syst: 0.03s (syst: 0.00s, elaps: 0.03s)
# ------------------------------------------------------------------------------------------
# ------------------------------------------------------------------------------------------
# Commande No : 0012 Concept de type : char_meca
# ------------------------------------------------------------------------------------------
displac = AFFE_CHAR_MECA(MODELE=Model,
DDL_IMPO=_F(GROUP_NO=('top', ),
DY=-5.E-3,),
INFO=1,
VERI_NORM='OUI',
)
# M��moire (Mo) : 342.30 / 342.07 / 27.33 / 21.10 (VmPeak / VmSize / Optimum / Minimum)
# Fin commande No : 0012 user+syst: 0.00s (syst: 0.00s, elaps: 0.00s)
# ------------------------------------------------------------------------------------------
# ------------------------------------------------------------------------------------------
# Commande No : 0013 Concept de type : listr8_sdaster
# ------------------------------------------------------------------------------------------
L_INST1 = DEFI_LIST_REEL(INFO=2,
DEBUT=0.0,
INTERVALLE=_F(NOMBRE=100,
JUSQU_A=1,),
)
-------------------------------------------------------------------------------
ASTER 11.07.00 CONCEPT L_INST1 CALCULE LE 20/05/2015 A 06:40:14 DE TYPE
LISTR8_SDASTER
INTERVALLE DEBUT JUSQU_A PAR_PAS NOMBRE
1 0.00000E+00 1.00000E+00 1.00000E-02 100
IMPRESSION DE LA LISTE DE REELS
1 - 0.000000000E+00 1.000000000E-02 2.000000000E-02 3.000000000E-02 4.000000000E-02
6 - 5.000000000E-02 6.000000000E-02 7.000000000E-02 8.000000000E-02 9.000000000E-02
11 - 1.000000000E-01 1.100000000E-01 1.200000000E-01 1.300000000E-01 1.400000000E-01
16 - 1.500000000E-01 1.600000000E-01 1.700000000E-01 1.800000000E-01 1.900000000E-01
21 - 2.000000000E-01 2.100000000E-01 2.200000000E-01 2.300000000E-01 2.400000000E-01
26 - 2.500000000E-01 2.600000000E-01 2.700000000E-01 2.800000000E-01 2.900000000E-01
31 - 3.000000000E-01 3.100000000E-01 3.200000000E-01 3.300000000E-01 3.400000000E-01
36 - 3.500000000E-01 3.600000000E-01 3.700000000E-01 3.800000000E-01 3.900000000E-01
41 - 4.000000000E-01 4.100000000E-01 4.200000000E-01 4.300000000E-01 4.400000000E-01
46 - 4.500000000E-01 4.600000000E-01 4.700000000E-01 4.800000000E-01 4.900000000E-01
51 - 5.000000000E-01 5.100000000E-01 5.200000000E-01 5.300000000E-01 5.400000000E-01
56 - 5.500000000E-01 5.600000000E-01 5.700000000E-01 5.800000000E-01 5.900000000E-01
61 - 6.000000000E-01 6.100000000E-01 6.200000000E-01 6.300000000E-01 6.400000000E-01
66 - 6.500000000E-01 6.600000000E-01 6.700000000E-01 6.800000000E-01 6.900000000E-01
71 - 7.000000000E-01 7.100000000E-01 7.200000000E-01 7.300000000E-01 7.400000000E-01
76 - 7.500000000E-01 7.600000000E-01 7.700000000E-01 7.800000000E-01 7.900000000E-01
81 - 8.000000000E-01 8.100000000E-01 8.200000000E-01 8.300000000E-01 8.400000000E-01
86 - 8.500000000E-01 8.600000000E-01 8.700000000E-01 8.800000000E-01 8.900000000E-01
91 - 9.000000000E-01 9.100000000E-01 9.200000000E-01 9.300000000E-01 9.400000000E-01
96 - 9.500000000E-01 9.600000000E-01 9.700000000E-01 9.800000000E-01 9.900000000E-01
101 - 1.000000000E+00
# M��moire (Mo) : 342.30 / 342.07 / 27.33 / 21.10 (VmPeak / VmSize / Optimum / Minimum)
# Fin commande No : 0013 user+syst: 0.00s (syst: 0.00s, elaps: 0.00s)
# ------------------------------------------------------------------------------------------
# ------------------------------------------------------------------------------------------
# Commande No : 0014 Concept de type : listr8_sdaster
# ------------------------------------------------------------------------------------------
saveT1 = DEFI_LIST_REEL(DEBUT=0.0,
INTERVALLE=_F(NOMBRE=100,
JUSQU_A=1,),
INFO=1,)
# M��moire (Mo) : 342.30 / 342.07 / 27.33 / 21.10 (VmPeak / VmSize / Optimum / Minimum)
# Fin commande No : 0014 user+syst: 0.00s (syst: 0.00s, elaps: 0.00s)
# ------------------------------------------------------------------------------------------
# ------------------------------------------------------------------------------------------
# Commande No : 0015 Concept de type : list_inst
# ------------------------------------------------------------------------------------------
DEFLIST1 = DEFI_LIST_INST(DEFI_LIST=_F(METHODE='AUTO',
LIST_INST=L_INST1,
PAS_MINI=1.E-10,
NB_PAS_MAXI=1000000,),
ECHEC=_F(EVENEMENT='ERREUR',
SUBD_PAS=4,
SUBD_NIVEAU=3,
SUBD_PAS_MINI=0.0,
SUBD_METHODE='MANUEL',
ACTION='DECOUPE',),
ADAPTATION=_F(EVENEMENT='SEUIL',
NB_INCR_SEUIL=3,
CRIT_COMP='LE',
NOM_PARA='NB_ITER_NEWTON',
PCENT_AUGM=200.0,
MODE_CALCUL_TPLUS='FIXE',),
INFO=1,)
# M��moire (Mo) : 342.30 / 342.07 / 27.33 / 21.10 (VmPeak / VmSize / Optimum / Minimum)
# Fin commande No : 0015 user+syst: 0.01s (syst: 0.00s, elaps: 0.01s)
# ------------------------------------------------------------------------------------------
# ------------------------------------------------------------------------------------------
# Commande No : 0016 Concept de type : fonction_sdaster
# ------------------------------------------------------------------------------------------
RampeP1 = DEFI_FONCTION(PROL_DROITE='CONSTANT',
VALE=(0.0, 0.0, 1, 1.0),
NOM_PARA='INST',
PROL_GAUCHE='CONSTANT',
INTERPOL=('LIN', ),
INFO=1,
NOM_RESU='TOUTRESU',
VERIF='CROISSANT',)
# M��moire (Mo) : 342.30 / 342.07 / 27.33 / 21.10 (VmPeak / VmSize / Optimum / Minimum)
# Fin commande No : 0016 user+syst: 0.00s (syst: 0.00s, elaps: 0.00s)
# ------------------------------------------------------------------------------------------
# ------------------------------------------------------------------------------------------
# Commande No : 0017 Concept de type : evol_noli
# ------------------------------------------------------------------------------------------
result = STAT_NON_LINE(CARA_ELEM=CARA,
CHAM_MATER=Mat,
MODELE=Model,
EXCIT=(_F(CHARGE=BCnd,
TYPE_CHARGE='FIXE_CSTE',),
_F(TYPE_CHARGE='FIXE_PILO',
CHARGE=displac,),
),
PILOTAGE=_F(SELECTION='ANGL_INCR_DEPL',
NOM_CMP=('DX', 'DY', 'DZ'),
GROUP_NO=('surface', ),
COEF_MULT=120.0,
TYPE='LONG_ARC',
ETA_PILO_MAX=1.0,
EVOL_PARA='SANS',
PROJ_BORNES='OUI',),
INCREMENT=_F(LIST_INST=DEFLIST1,
INST_FIN=0.1,
PRECISION=1.E-06,),
CONVERGENCE=_F(RESI_GLOB_RELA=1.E-06,
ITER_GLOB_MAXI=30,
ARRET='OUI',
ITER_GLOB_ELAS=25,),
COMP_ELAS=_F(DEFORMATION='GROT_GDEP',
RELATION='ELAS',
RESI_INTE_RELA=1.E-06,
ITER_INTE_MAXI=10,),
NEWTON=_F(MATRICE='TANGENTE',
REAC_ITER=1,
REAC_INCR=1,
REAC_ITER_ELAS=0,),
METHODE='NEWTON',
ARCHIVAGE=_F(PRECISION=1.E-06,
CRITERE='RELATIF',),
SOLVEUR=_F(RENUM='METIS',
STOP_SINGULIER='OUI',
NPREC=8,
METHODE='MULT_FRONT',
SYME='NON',),
)
Mot-clef : COMP_ELAS - occurrence 1 : comportement ELAS - nombre de variables internes : 1
Noms des variables internes :
V1 : VIDE
Le syst��me lin��aire �� r��soudre contient 3763 noeuds dont:
- 2993 noeuds portant des degr��s de libert�� physiques
- 770 noeuds portant des degr��s de libert�� de Lagrange
Pour un total de 16532 ��quations.
La matrice est de taille 16532 ��quations.
Elle contient 746256 termes non nuls si elle est sym��trique et 1475980 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.540 %.
!---------------------------------------------------------------!
! <EXCEPTION> <MECANONLINE5_50> !
! !
! Pilotage. !
! La composante <DX> n'a pas ��t�� trouv��e dans la num��rotation. !
! V��rifier NOM_CMP dans le mot-clef PILOTAGE. !
!---------------------------------------------------------------!
Validation du concept 'result'.
# ------------------------------------------------------------------------------------------
# Commande No : 0018 Concept de type : -
# ------------------------------------------------------------------------------------------
FIN(RETASSAGE='NON',
INFO_RESU='OUI',
FORMAT_HDF='NON',
UNITE=6,)
!-------------------------------------------------------------------------!
! <I> <CATAMESS_89> !
! !
! 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 : !
! aucune alarme !
!-------------------------------------------------------------------------!
---------------------------------------------------------------------------
Concepts de la base: G
Nom Type Taille (Mo) Nombre Nombre de
d'objets segments
TOTAL 16.99 240 428
meshI MAILLAGE_SDASTER 0.11 35 53
Qmesh MAILLAGE_SDASTER 0.19 34 52
mesh MAILLAGE_SDASTER 0.25 36 67
Model MODELE_SDASTER 0.16 8 13
Alu MATER_SDASTER 0.00 4 4
Mat CHAM_MATER 0.00 5 10
CARA CARA_ELEM 0.03 8 13
BCnd CHAR_MECA 0.26 19 39
displac CHAR_MECA 0.12 19 39
L_INST1 LISTR8_SDASTER 0.00 5 5
saveT1 LISTR8_SDASTER 0.00 5 5
DEFLIST1 LIST_INST 0.00 9 9
RampeP1 FONCTION_SDASTER 0.00 3 3
result EVOL_NOLI 0.58 5 10
&FOZERO 0.00 2 2
&&_NUM_C 0.00 1 1
&CATA.AC 0.00 2 4
&CATA.CL 0.26 1 3
&CATA.GD 0.04 4 11
&CATA.ME 0.14 2 4
&CATA.OP 0.28 4 19
&CATA.PH 0.00 1 1
&CATA.PR 0.00 2 4
&CATA.TE 14.53 17 42
&CATA.TH 0.01 2 4
&CATA.TM 0.02 7 11
---------------------------------------------------------------------------
Fermeture de la base
Nom de la base : GLOBALE
Nombre d'enregistrements utilis��s : 29
Nombre d'enregistrements maximum : 62914
Longueur d'enregistrement (octets) : 819200
Nombre total d'acc��s en lecture : 7
Volume des acc��s en lecture : 5.47 Mo.
Nombre total d'acc��s en ��criture : 43
Volume des acc��s en ��criture : 33.59 Mo.
Nombre d'identificateurs utilis��s : 425
Taille maximum du r��pertoire : 2000
Pourcentage d'utilisation du r��pertoire : 21 %
Fermeture de la base
Nom de la base : VOLATILE
Nombre d'enregistrements utilis��s : 3
Nombre d'enregistrements maximum : 62914
Longueur d'enregistrement (octets) : 819200
Nombre total d'acc��s en lecture : 0
Volume des acc��s en lecture : 0.00 Mo.
Nombre total d'acc��s en ��criture : 4
Volume des acc��s en ��criture : 3.12 Mo.
Nombre d'identificateurs utilis��s : 626
Taille maximum du r��pertoire : 2000
Pourcentage d'utilisation du r��pertoire : 31 %
<FIN> Arr��t normal dans "FIN".
La m��moire demand��e au lancement est surestim��e, elle est de 15120.00 Mo.
Le pic m��moire utilis��e est de 363.84 Mo.
<I> FERMETURE DES BASES EFFECTUEE
STATISTIQUES CONCERNANT L'ALLOCATION DYNAMIQUE :
TAILLE CUMULEE MAXIMUM : 45 Mo.
TAILLE CUMULEE LIBEREE : 18 Mo.
NOMBRE TOTAL D'ALLOCATIONS : 28483
NOMBRE TOTAL DE LIBERATIONS : 28483
APPELS AU MECANISME DE LIBERATION : 0
TAILLE MEMOIRE CUMULEE RECUPEREE : 0 Mo.
VOLUME DES LECTURES : 0 Mo.
VOLUME DES ECRITURES : 0 Mo.
MEMOIRE JEVEUX MINIMALE REQUISE POUR L'EXECUTION : 28.22 Mo
- IMPOSE DE NOMBREUX ACCES DISQUE
- RALENTIT LA VITESSE D'EXECUTION
MEMOIRE JEVEUX OPTIMALE REQUISE POUR L'EXECUTION : 45.44 Mo
- LIMITE LES ACCES DISQUE
- AMELIORE LA VITESSE D'EXECUTION
MAXIMUM DE MEMOIRE UTILISEE PAR LE PROCESSUS : 363.84 Mo
- COMPREND LA MEMOIRE CONSOMMEE PAR JEVEUX,
LE SUPERVISEUR PYTHON, LES LIBRAIRIES EXTERNES
<I> FIN D'EXECUTION LE : ME-20-MAI -2015 06:40:14
<I> Informations sur les temps d'ex��cution
Temps cpu total 0.62 s
Temps cpu user total 0.47 s
Temps cpu systeme total 0.15 s
Temps cpu restant 107999.38 s
La signature de la base sauvegard��e est (�� l'adresse 0) :
2daa2b29a61a4f6b9721af8d7190d81759354fd2
********************************************************************************
* COMMAND : USER : SYSTEM : USER+SYS : ELAPSED *
********************************************************************************
* init (jdc) : 0.21 : 0.03 : 0.24 : 0.24 *
* . compile : 0.00 : 0.00 : 0.00 : 0.00 *
* . exec_compile : 0.05 : 0.01 : 0.06 : 0.06 *
* . report : 0.01 : 0.00 : 0.01 : 0.01 *
* . build : 0.00 : 0.00 : 0.00 : 0.00 *
* DEBUT : 0.02 : 0.04 : 0.06 : 0.05 *
* LIRE_MAILLAGE : 0.01 : 0.00 : 0.01 : 0.01 *
* CREA_MAILLAGE : 0.01 : 0.00 : 0.01 : 0.01 *
* CREA_MAILLAGE : 0.02 : 0.00 : 0.02 : 0.02 *
* DEFI_GROUP : 0.00 : 0.00 : 0.00 : 0.00 *
* MODI_MAILLAGE : 0.01 : 0.00 : 0.01 : 0.02 *
* AFFE_MODELE : 0.01 : 0.01 : 0.02 : 0.01 *
* DEFI_MATERIAU : 0.00 : 0.00 : 0.00 : 0.00 *
* AFFE_MATERIAU : 0.00 : 0.00 : 0.00 : 0.00 *
* AFFE_CARA_ELEM : 0.01 : 0.00 : 0.01 : 0.01 *
* AFFE_CHAR_MECA : 0.03 : 0.00 : 0.03 : 0.03 *
* AFFE_CHAR_MECA : 0.00 : 0.00 : 0.00 : 0.00 *
* DEFI_LIST_REEL : 0.00 : 0.00 : 0.00 : 0.00 *
* DEFI_LIST_REEL : 0.00 : 0.00 : 0.00 : 0.00 *
* DEFI_LIST_INST : 0.01 : 0.00 : 0.01 : 0.01 *
* DEFI_FONCTION : 0.00 : 0.00 : 0.00 : 0.00 *
* STAT_NON_LINE : 0.12 : 0.07 : 0.19 : 0.19 *
* FIN : 0.01 : 0.04 : 0.05 : 0.04 *
* . part Superviseur : 0.26 : 0.07 : 0.33 : 0.32 *
* . part Fortran : 0.21 : 0.08 : 0.29 : 0.30 *
********************************************************************************
* TOTAL_JOB : 0.47 : 0.15 : 0.62 : 0.62 *
********************************************************************************
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. !
! Type de l'exception : error !
! !
! Pilotage. !
! La composante <DX> n'a pas ��t�� trouv��e dans la num��rotation. !
! V��rifier NOM_CMP dans le mot-clef PILOTAGE. !
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
fin CR d'execution de JDC en MIXTE
>> JDC.py : FIN RAPPORT
EXECUTION_CODE_ASTER_EXIT_14987=1
Last edit: 8 years 11 months ago by alessandro.
Please Log in or Create an account to join the conversation.
- RichardS
- Offline
- Elite Member
-
8 years 11 months ago #8099
by RichardS
SimScale - Engineering Simulation in your browser!
Replied by RichardS on topic Re: error using pilotage
Hello,
I never really used PILOTAGE,
but my guess is that aster complains that you want to use pilotage for the DOF DX but it cannot find this DOF in the hassigned boundary condition.
Since you only impose a defined displacemeent in DY directrion, I would suggest that you have to remove DX and DZ fromin the PILOTAGE settings.
Best,
Richard
I never really used PILOTAGE,
but my guess is that aster complains that you want to use pilotage for the DOF DX but it cannot find this DOF in the hassigned boundary condition.
Since you only impose a defined displacemeent in DY directrion, I would suggest that you have to remove DX and DZ from
NOM_CMP = ('DX', 'DY', 'DZ',),Best,
Richard
SimScale - Engineering Simulation in your browser!
Please Log in or Create an account to join the conversation.
- alessandro
- Topic Author
- Offline
- New Member
-
Less
More
- Posts: 10
- Thank you received: 0
8 years 11 months ago - 8 years 11 months ago #8100
by alessandro
Replied by alessandro on topic Re: error using pilotage
Hi, thanks for the answer.
I tried that but the problem remains in the form
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! <S> Exception utilisateur levee mais pas interceptee. !
! Les bases sont fermees. !
! Type de l'exception : error !
! !
! Pilotage. !
! La composante <DY> n'a pas été trouvée dans la numérotation. !
! Vérifier NOM_CMP dans le mot-clef PILOTAGE. !
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
also the command with 'DX','DY','DZ' is used in other examples.
For instance I did the case SSNS501 and there it works fine, here is the .mess file for that case
[code]
#
# Impression du contenu du fichier de commandes �� ex��cuter :
#
# person_in_charge: patrick.massin at edf.fr
# DATE : 12/07/2000
# AUTEUR : F. LEBOUVIER (DELTACAD)
# TITRE : GRANDS DEPLACEMENTS D'UN PANNEAU CYLINDRIQUE
# ======================================================================
# COPYRIGHT (C) 1991 - 2012 EDF R&D WWW.CODE-ASTER.ORG
# THIS PROGRAM IS FREE SOFTWARE; YOU CAN REDISTRIBUTE IT AND/OR MODIFY
# IT UNDER THE TERMS OF THE GNU GENERAL PUBLIC LICENSE AS PUBLISHED BY
# THE FREE SOFTWARE FOUNDATION; EITHER VERSION 2 OF THE LICENSE, OR
# (AT YOUR OPTION) ANY LATER VERSION.
#
# THIS PROGRAM IS DISTRIBUTED IN THE HOPE THAT IT WILL BE USEFUL, BUT
# WITHOUT ANY WARRANTY; WITHOUT EVEN THE IMPLIED WARRANTY OF
# MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. SEE THE GNU
# GENERAL PUBLIC LICENSE FOR MORE DETAILS.
#
# YOU SHOULD HAVE RECEIVED A COPY OF THE GNU GENERAL PUBLIC LICENSE
# ALONG WITH THIS PROGRAM; IF NOT, WRITE TO EDF R&D CODE_ASTER,
# 1 AVENUE DU GENERAL DE GAULLE, 92141 CLAMART CEDEX, FRANCE.
# ======================================================================
#
# - CALCUL DE STABILITE D'UN PANNEAU CYLINDRIQUE SIMPLEMENT SUPPORTE
# SOUMIS A UN EFFORT CONCENTRE EN SON CENTRE. VALIDATION DANS LE
# DOMAINE QUASI_STATIQUE NON_LINEAIRE GEOMETRIQUE EN PRESENCE DE
# FORTES INSTABILITES
# - MODELISATION AVEC DES ELEMENTS 'COQUE_3D' (QUAD9)
#
#
DEBUT( )
#
MAIL=LIRE_MAILLAGE(FORMAT='MED', )
MAIL=DEFI_GROUP( reuse=MAIL,
MAILLAGE=MAIL,CREA_GROUP_NO=(
_F( GROUP_MA = 'TOUT_ELT', ),)
)
#
MAILLAGE=CREA_MAILLAGE( MAILLAGE=MAIL,
MODI_MAILLE=_F( TOUT = 'OUI',
OPTION = 'QUAD8_9')
)
#
MODELE=AFFE_MODELE( MAILLAGE=MAILLAGE,
AFFE=_F( TOUT = 'OUI',
PHENOMENE = 'MECANIQUE',
MODELISATION = 'COQUE_3D')
)
#
#
# CARACTERISTIQUES ELEMENTAIRES
#
#
CARA_ELT=AFFE_CARA_ELEM( MODELE=MODELE,
COQUE=_F( GROUP_MA = 'TOUT_ELT',
EPAIS = 0.00635,
# A_CIS : 2500.
COEF_RIGI_DRZ = 0.001)
)
#
#
# CARACTERISTIQUES MATERIAUX
#
#
MATERIAU=DEFI_MATERIAU( ELAS=_F( E = 3.10275E+09,
NU = 0.3)
)
#
CHAM_MAT=AFFE_MATERIAU( MAILLAGE=MAILLAGE,
AFFE=_F( TOUT = 'OUI', MATER = MATERIAU)
)
#
#
# CHARGEMENTS
#
#
COND_LIM=AFFE_CHAR_MECA( MODELE=MODELE,DDL_IMPO=(
_F( GROUP_NO = 'COTEAB',
DX = 0., DRY = 0., DRZ = 0.),
_F( GROUP_NO = 'COTECD',
DX = 0., DY = 0., DZ = 0.),
_F( GROUP_NO = 'COTEAD',
DY = 0., DRX = 0., DRZ = 0.))
)
#
CHARGE=AFFE_CHAR_MECA( MODELE=MODELE,
#FORCE_NODALE=_F( GROUP_NO = 'A',FZ = -500.)
DDL_IMPO =_F(GROUP_NO='A',DZ=-0.1 )
)
F_CHARGE=DEFI_FONCTION( NOM_PARA='INST',
PROL_DROITE='CONSTANT',
PROL_GAUCHE='EXCLU',
VALE=( 0., 0.,
1., 1.,
)
)
#
#
# RESOLUTION
#
#
LIS_INST=DEFI_LIST_REEL(DEBUT=0.0,INTERVALLE=(
_F( JUSQU_A = 10.0, NOMBRE = 100))
)
RESULTAT=STAT_NON_LINE(MODELE=MODELE,
CHAM_MATER=CHAM_MAT,
CARA_ELEM=CARA_ELT,
EXCIT=(
_F( CHARGE = COND_LIM),
_F( CHARGE = CHARGE,
TYPE_CHARGE = 'FIXE_PILO')
),
COMP_ELAS=_F( RELATION = 'ELAS',
DEFORMATION = 'GROT_GDEP',
),
INCREMENT=_F( LIST_INST = LIS_INST,),
PILOTAGE=_F( TYPE = 'LONG_ARC',
SELECTION='ANGL_INCR_DEPL',
GROUP_NO = 'TOUT_ELT',
NOM_CMP = ('DX', 'DY', 'DZ',),
ETA_PILO_MAX=1.,
COEF_MULT = 120.),
NEWTON=_F( REAC_ITER = 1,
MATRICE = 'TANGENTE'),
CONVERGENCE=_F(
ITER_GLOB_MAXI = 30,
RESI_GLOB_RELA = 1.E-6)
)
#
# calc more fields
RESULTAT=CALC_CHAMP(reuse =RESULTAT,
MODELE=MODELE,
RESULTAT=RESULTAT,
FORCE=('FORC_NODA','REAC_NODA'),
TOUT='OUI',
);
#
# IMPRESSION DES RESULTATS
#
DEPZ=RECU_FONCTION ( RESULTAT = RESULTAT,
NOM_CMP = 'DZ',
NOM_CHAM = 'DEPL',
GROUP_NO = 'A',
TOUT_ORDRE = 'OUI')
REAZ=RECU_FONCTION ( RESULTAT = RESULTAT,
NOM_CMP = 'DZ',
NOM_CHAM = 'REAC_NODA',
GROUP_NO = 'A',
TOUT_ORDRE = 'OUI')
FORZ=RECU_FONCTION ( RESULTAT = RESULTAT,
NOM_CMP = 'DZ',
NOM_CHAM = 'FORC_NODA',
GROUP_NO = 'A',
TOUT_ORDRE = 'OUI')
IMPR_FONCTION(FORMAT='TABLEAU',UNITE=82,COMMENTAIRE = '%',COMM_PARA = '%',
COURBE=_F( FONCTION = REAZ),);
IMPR_FONCTION(FORMAT='TABLEAU',UNITE=83,COMMENTAIRE = '%',COMM_PARA = '%',
COURBE=_F(FONC_X=DEPZ,FONC_Y=FORZ,),);
ETA=RECU_FONCTION (RESULTAT = RESULTAT,
NOM_PARA_RESU = 'ETA_PILOTAGE',
TOUT_ORDRE = 'OUI')
#IMPR_FONCTION (FORMAT='AGRAF', UNITE=27, UNITE_DIGR=28,COURBE=_F( FONCTION = ETA),TITRE='Eta en fonction du temps')
IMPR_FONCTION(FORMAT='TABLEAU',
UNITE=81,
COURBE=_F( FONCTION = ETA),
COMMENTAIRE = '%',
COMM_PARA = '%',
);
IMPR_RESU(FORMAT='MED',UNITE=80,RESU=_F(MAILLAGE=MAIL,RESULTAT=RESULTAT,NOM_CHAM=('DEPL','REAC_NODA','FORC_NODA'),),);
#
FIN()
#
#
#
-- CODE_ASTER -- VERSION : EXPLOITATION (stable) --
Version 11.7.0 modifi��e le 10/12/2014
r��vision 7a39a0fb6aea - branche 'v11'
Copyright EDF R&D 1991 - 2015
Ex��cution du : Tue May 19 07:57:03 2015
Nom de la machine : naxos
Architecture : 64bit
Type de processeur : x86_64
Syst��me d'exploitation : Linux CentOS 6.4 Final 2.6.32-358.23.2.el6.x86_64
Langue des messages : en (UTF-8)
Version de Python : 2.7.3
Version de NumPy : 1.8.1
Parall��lisme MPI : inactif
Parall��lisme OpenMP : actif
Nombre de processus utilis��s : 1
Version de la librairie HDF5 : 1.8.10
Version de la librairie MED : 3.0.8
Version de la librairie MUMPS : 4.10.0
Version de la librairie SCOTCH : 5.1.10
M��moire limite pour l'ex��cution : 15120.00 Mo
consomm��e par l'initialisation : 307.90 Mo
par les objets du jeu de commandes : 2.05 Mo
reste pour l'allocation dynamique : 14809.80 Mo
Taille limite des fichiers d'��change : 48.00 Go
#
# Commande No : 0001 Concept de type : -
#
DEBUT(INFO=1,
PAR_LOT='OUI',
IMPR_MACRO='NON',
MEMOIRE=_F(TAILLE_BLOC=800.0,
TAILLE_GROUP_ELEM=1000,),
RESERVE_CPU=_F(BORNE=900,),
MESURE_TEMPS=_F(MOYENNE='OUI',
NIVE_DETAIL=1,),
)
Valeur initiale du temps CPU maximum = 108900 secondes
Valeur du temps CPU maximum pass�� aux commandes = 108000 secondes
R��serve CPU pr��vue = 900 secondes
D��but de lecture...
R��ouverture de la base
Nom de la base : ELEMBASE
Cr����e avec la version : 11.07.00
Nombre d'enregistrements utilis��s : 24
Nombre d'enregistrements maximum : 512
Longueur d'enregistrement (octets) : 819200
Nombre d'identificateurs utilis��s : 123
Taille maximum du r��pertoire : 300
Pourcentage d'utilisation du r��pertoire : 41 %
Fermeture de la base
Nom de la base : ELEMBASE
Nombre d'enregistrements utilis��s : 24
Nombre d'enregistrements maximum : 512
Longueur d'enregistrement (octets) : 819200
Nombre total d'acc��s en lecture : 49
Volume des acc��s en lecture : 38.28 Mo.
Nombre total d'acc��s en ��criture : 3
Volume des acc��s en ��criture : 2.34 Mo.
Nombre d'identificateurs utilis��s : 123
Taille maximum du r��pertoire : 300
Pourcentage d'utilisation du r��pertoire : 41 %
Relecture des catalogues des ��l��ments faite.
Fin de lecture (dur��e 0.010000 s.)
# M��moire (Mo) : 337.97 / 333.50 / 27.33 / 14.95 (VmPeak / VmSize / Optimum / Minimum)
# Fin commande No : 0001 user+syst: 0.05s (syst: 0.03s, elaps: 0.06s)
#
#
# Commande No : 0002 Concept de type : maillage_sdaster
#
MAIL = LIRE_MAILLAGE(FORMAT='MED',
INFO_MED=1,
INFO=1,
UNITE=20,
VERI_MAIL=_F(APLAT=1.E-3,
VERIF='OUI',),
)
====== VERIFICATION DU MAILLAGE ======
MAILLAGE MAIL - IMPRESSIONS NIVEAU 1
CODE_ASTER - 11.01.04 - VERSION DE DEVELOPPEMENT - FRONTAL1 - L
NOMBRE DE NOEUDS 133
NOMBRE DE MAILLES 36
QUAD8 36
NOMBRE DE GROUPES DE NOEUDS 4
COTEAB 13
COTEAD 13
A 1
COTECD 13
NOMBRE DE GROUPES DE MAILLES 2
COUL_7 36
TOUT_ELT 36
# M��moire (Mo) : 337.97 / 334.63 / 27.33 / 14.95 (VmPeak / VmSize / Optimum / Minimum)
# Fin commande No : 0002 user+syst: 0.01s (syst: 0.00s, elaps: 0.01s)
#
#
# Commande No : 0003 Concept de type : maillage_sdaster
#
MAIL = DEFI_GROUP(reuse=MAIL,
MAILLAGE=MAIL,
CREA_GROUP_NO=_F(GROUP_MA=('TOUT_ELT', ),
CRIT_NOEUD='TOUS',),
ALARME='OUI',)
NOMBRE DE GROUPES DE NOEUDS CREES : 1
=======================================
! NOM DU GROUPE ! NBRE DE NOEUDS DU !
! NOEUDS ! GROUPE_NO !
! TOUT_ELT ! 133 !
# M��moire (Mo) : 337.97 / 334.63 / 27.33 / 14.95 (VmPeak / VmSize / Optimum / Minimum)
# Fin commande No : 0003 user+syst: 0.00s (syst: 0.00s, elaps: 0.00s)
#
#
# Commande No : 0004 Concept de type : maillage_sdaster
#
MAILLAGE = CREA_MAILLAGE(MAILLAGE=MAIL,
MODI_MAILLE=_F(OPTION='QUAD8_9',
TOUT='OUI',
PREF_NOEUD='NS',
PREF_NUME=1,),
INFO=1,)
MOT CLE FACTEUR "MODI_MAILLE", OCCURRENCE 1
MODIFICATION DE 36 MAILLES QUAD8 EN QUAD9
MAILLAGE MAILLAGE - IMPRESSIONS NIVEAU 1
ASTER 11.07.00 CONCEPT MAILLAGE CALCULE LE 19/05/2015 A 07:57:03 DE TYPE
MAILLAGE_SDASTER
NOMBRE DE NOEUDS 169
NOMBRE DE MAILLES 36
QUAD9 36
NOMBRE DE GROUPES DE NOEUDS 5
COTEAB 13
COTEAD 13
A 1
COTECD 13
TOUT_ELT 133
NOMBRE DE GROUPES DE MAILLES 2
COUL_7 36
TOUT_ELT 36
# M��moire (Mo) : 337.97 / 334.77 / 27.33 / 14.95 (VmPeak / VmSize / Optimum / Minimum)
# Fin commande No : 0004 user+syst: 0.01s (syst: 0.00s, elaps: 0.00s)
#
#
# Commande No : 0005 Concept de type : modele_sdaster
#
MODELE = AFFE_MODELE(MAILLAGE=MAILLAGE,
AFFE=_F(PHENOMENE='MECANIQUE',
TOUT='OUI',
MODELISATION=('COQUE_3D', ),),
INFO=1,
VERI_JACOBIEN='OUI',
PARTITION=_F(PARALLELISME='GROUP_ELEM',),
)
SUR LES 36 MAILLES DU MAILLAGE MAILLAGE
ON A DEMANDE L'AFFECTATION DE 36
ON A PU EN AFFECTER 36
MODELISATION ELEMENT FINI TYPE MAILLE NOMBRE
COQUE_3D MEC3QU9H QUAD9 36
#2 Calculs elementaires et assemblages CPU (USER+SYST/SYST/ELAPS): 0.00 0.00 0.00
# M��moire (Mo) : 341.52 / 341.50 / 27.33 / 20.91 (VmPeak / VmSize / Optimum / Minimum)
# Fin commande No : 0005 user+syst: 0.00s (syst: 0.00s, elaps: 0.01s)
#
#
# Commande No : 0006 Concept de type : cara_elem
#
CARA_ELT = AFFE_CARA_ELEM(MODELE=MODELE,
COQUE=_F(EPAIS=6.35E-3,
GROUP_MA=('TOUT_ELT', ),
COEF_RIGI_DRZ=1.E-3,
A_CIS=0.8333333,
MODI_METRIQUE='NON',
COQUE_NCOU=1,),
INFO=1,)
LE MODELE MODELE CONTIENT :
36 ELEMENT(S) COQUE
#2 Calculs elementaires et assemblages CPU (USER+SYST/SYST/ELAPS): 0.00 0.00 0.00
# M��moire (Mo) : 341.52 / 341.50 / 27.33 / 20.91 (VmPeak / VmSize / Optimum / Minimum)
# Fin commande No : 0006 user+syst: 0.01s (syst: 0.00s, elaps: 0.01s)
#
#
# Commande No : 0007 Concept de type : mater_sdaster
#
MATERIAU = DEFI_MATERIAU(ELAS=_F(E=3.102750000E9,
NU=0.3,),
INFO=1,)
MATERIAU : MATERIAU
RELATION DE COMPORTEMENT: ELAS
# M��moire (Mo) : 341.52 / 341.50 / 27.33 / 20.91 (VmPeak / VmSize / Optimum / Minimum)
# Fin commande No : 0007 user+syst: 0.00s (syst: 0.00s, elaps: 0.00s)
#
#
# Commande No : 0008 Concept de type : cham_mater
#
CHAM_MAT = AFFE_MATERIAU(MAILLAGE=MAILLAGE,
AFFE=_F(MATER=(MATERIAU, ),
TOUT='OUI',),
INFO=1,
)
# M��moire (Mo) : 341.52 / 341.50 / 27.33 / 20.91 (VmPeak / VmSize / Optimum / Minimum)
# Fin commande No : 0008 user+syst: 0.01s (syst: 0.00s, elaps: 0.01s)
#
#
# Commande No : 0009 Concept de type : char_meca
#
COND_LIM = AFFE_CHAR_MECA(MODELE=MODELE,
DDL_IMPO=(_F(DRY=0.0,
DRZ=0.0,
GROUP_NO=('COTEAB', ),
DX=0.0,),
_F(GROUP_NO=('COTECD', ),
DZ=0.0,
DX=0.0,
DY=0.0,),
_F(DRX=0.0,
DRZ=0.0,
GROUP_NO=('COTEAD', ),
DY=0.0,),
),
INFO=1,
VERI_NORM='OUI',
)
# M��moire (Mo) : 341.74 / 341.50 / 27.33 / 20.91 (VmPeak / VmSize / Optimum / Minimum)
# Fin commande No : 0009 user+syst: 0.01s (syst: 0.00s, elaps: 0.01s)
#
#
# Commande No : 0010 Concept de type : char_meca
#
CHARGE = AFFE_CHAR_MECA(MODELE=MODELE,
DDL_IMPO=_F(GROUP_NO=('A', ),
DZ=-0.1,),
INFO=1,
VERI_NORM='OUI',
)
# M��moire (Mo) : 341.75 / 341.50 / 27.33 / 20.91 (VmPeak / VmSize / Optimum / Minimum)
# Fin commande No : 0010 user+syst: 0.00s (syst: 0.00s, elaps: 0.00s)
#
#
# Commande No : 0011 Concept de type : fonction_sdaster
#
F_CHARGE = DEFI_FONCTION(PROL_DROITE='CONSTANT',
VALE=(0.0, 0.0, 1.0, 1.0),
NOM_PARA='INST',
PROL_GAUCHE='EXCLU',
INFO=1,
INTERPOL=('LIN', ),
NOM_RESU='TOUTRESU',
VERIF='CROISSANT',)
# M��moire (Mo) : 341.75 / 341.64 / 27.33 / 20.91 (VmPeak / VmSize / Optimum / Minimum)
# Fin commande No : 0011 user+syst: 0.00s (syst: 0.00s, elaps: 0.01s)
#
#
# Commande No : 0012 Concept de type : listr8_sdaster
#
LIS_INST = DEFI_LIST_REEL(DEBUT=0.0,
INTERVALLE=_F(NOMBRE=100,
JUSQU_A=10.0,),
INFO=1,)
# M��moire (Mo) : 341.75 / 341.64 / 27.33 / 20.91 (VmPeak / VmSize / Optimum / Minimum)
# Fin commande No : 0012 user+syst: 0.01s (syst: 0.00s, elaps: 0.00s)
#
#
# Commande No : 0013 Concept de type : evol_noli
#
RESULTAT = STAT_NON_LINE(CARA_ELEM=CARA_ELT,
CHAM_MATER=CHAM_MAT,
MODELE=MODELE,
EXCIT=(_F(CHARGE=COND_LIM,
TYPE_CHARGE='FIXE_CSTE',),
_F(TYPE_CHARGE='FIXE_PILO',
CHARGE=CHARGE,),
),
PILOTAGE=_F(SELECTION='ANGL_INCR_DEPL',
NOM_CMP=('DX', 'DY', 'DZ'),
GROUP_NO=('TOUT_ELT', ),
COEF_MULT=120.0,
TYPE='LONG_ARC',
ETA_PILO_MAX=1.0,
EVOL_PARA='SANS',
PROJ_BORNES='OUI',),
INCREMENT=_F(LIST_INST=LIS_INST,
PRECISION=1.E-06,),
CONVERGENCE=_F(RESI_GLOB_RELA=1.E-06,
ITER_GLOB_MAXI=30,
ARRET='OUI',
ITER_GLOB_ELAS=25,),
COMP_ELAS=_F(DEFORMATION='GROT_GDEP',
RELATION='ELAS',
RESI_INTE_RELA=1.E-06,
ITER_INTE_MAXI=10,),
NEWTON=_F(MATRICE='TANGENTE',
REAC_ITER=1,
REAC_INCR=1,
REAC_ITER_ELAS=0,),
METHODE='NEWTON',
ARCHIVAGE=_F(PRECISION=1.E-06,
CRITERE='RELATIF',),
SOLVEUR=_F(RENUM='METIS',
STOP_SINGULIER='OUI',
NPREC=8,
METHODE='MULT_FRONT',
SYME='NON',),
)
Mot-clef : COMP_ELAS - occurrence 1 : comportement ELAS - nombre de variables internes : 1
Noms des variables internes :
V1 : VIDE
Le syst��me lin��aire �� r��soudre contient 401 noeuds dont:
- 169 noeuds portant des degr��s de libert�� physiques
- 232 noeuds portant des degr��s de libert�� de Lagrange
Pour un total de 1138 ��quations.
La matrice est de taille 1138 ��quations.
Elle contient 39741 termes non nuls si elle est sym��trique et 78344 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 6.050 %.
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
Instant de calcul: 1.000000000000e-01 - Niveau de d��coupe: 1
| NEWTON | RESIDU | RESIDU | PILOTAGE | OPTION |
| ITERATION | RELATIF | ABSOLU | COEFFICIENT | ASSEMBLAGE |
| | RESI_GLOB_RELA | RESI_GLOB_MAXI | ETA | |
| 0 X | 2.56881E-01 X | 7.07797E+01 | 2.71212E-02 |TANGENTE |
| 1 X | 2.08575E-03 X | 4.24425E-01 | 2.64308E-02 |TANGENTE |
| 2 | 3.87930E-07 | 7.88830E-05 | 2.64233E-02 |TANGENTE |
Crit��re(s) de convergence atteint(s)
Le r��sidu de type <RESI_GLOB_RELA> vaut 3.879296069712e-07 au noeud et degr�� de libert�� <NO102 DY>
Temps CPU consomm�� dans ce pas de temps : 0.430 s
* Temps moyen par it��ration de Newton : 0.143 s ( 3 it��rations )
* Temps total int��gration comportement : 0.340 s ( 4 int��grations )
* Temps assemblage matrice : 0.000 s
* Temps total factorisation matrice : 0.020 s ( 3 factorisations )
* Temps construction second membre : 0.040 s
* Temps total r��solution K.U=F : 0.010 s ( 3 r��solutions )
* Temps autres op��rations : 0.020 s
M��moire (Mo) : 348.62 / 346.99 / 29.39 / 23.14 (VmPeak / VmSize / Optimum / Minimum)
Archivage des champs
Champ stock�� <DEPL> �� l'instant 1.000000000000e-01 pour le num��ro d'ordre 1
Champ stock�� <SIEF_ELGA> �� l'instant 1.000000000000e-01 pour le num��ro d'ordre 1
Champ stock�� <VARI_ELGA> �� l'instant 1.000000000000e-01 pour le num��ro d'ordre 1
Champ stock�� <COMPORTEMENT> �� l'instant 1.000000000000e-01 pour le num��ro d'ordre 1
...
EXECUTION_CODE_ASTER_EXIT_10486=0
[code]
I tried that but the problem remains in the form
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! <S> Exception utilisateur levee mais pas interceptee. !
! Les bases sont fermees. !
! Type de l'exception : error !
! !
! Pilotage. !
! La composante <DY> n'a pas été trouvée dans la numérotation. !
! Vérifier NOM_CMP dans le mot-clef PILOTAGE. !
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
also the command with 'DX','DY','DZ' is used in other examples.
For instance I did the case SSNS501 and there it works fine, here is the .mess file for that case
[code]
#
# Impression du contenu du fichier de commandes �� ex��cuter :
#
# person_in_charge: patrick.massin at edf.fr
# DATE : 12/07/2000
# AUTEUR : F. LEBOUVIER (DELTACAD)
# TITRE : GRANDS DEPLACEMENTS D'UN PANNEAU CYLINDRIQUE
# ======================================================================
# COPYRIGHT (C) 1991 - 2012 EDF R&D WWW.CODE-ASTER.ORG
# THIS PROGRAM IS FREE SOFTWARE; YOU CAN REDISTRIBUTE IT AND/OR MODIFY
# IT UNDER THE TERMS OF THE GNU GENERAL PUBLIC LICENSE AS PUBLISHED BY
# THE FREE SOFTWARE FOUNDATION; EITHER VERSION 2 OF THE LICENSE, OR
# (AT YOUR OPTION) ANY LATER VERSION.
#
# THIS PROGRAM IS DISTRIBUTED IN THE HOPE THAT IT WILL BE USEFUL, BUT
# WITHOUT ANY WARRANTY; WITHOUT EVEN THE IMPLIED WARRANTY OF
# MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. SEE THE GNU
# GENERAL PUBLIC LICENSE FOR MORE DETAILS.
#
# YOU SHOULD HAVE RECEIVED A COPY OF THE GNU GENERAL PUBLIC LICENSE
# ALONG WITH THIS PROGRAM; IF NOT, WRITE TO EDF R&D CODE_ASTER,
# 1 AVENUE DU GENERAL DE GAULLE, 92141 CLAMART CEDEX, FRANCE.
# ======================================================================
#
# - CALCUL DE STABILITE D'UN PANNEAU CYLINDRIQUE SIMPLEMENT SUPPORTE
# SOUMIS A UN EFFORT CONCENTRE EN SON CENTRE. VALIDATION DANS LE
# DOMAINE QUASI_STATIQUE NON_LINEAIRE GEOMETRIQUE EN PRESENCE DE
# FORTES INSTABILITES
# - MODELISATION AVEC DES ELEMENTS 'COQUE_3D' (QUAD9)
#
#
DEBUT( )
#
MAIL=LIRE_MAILLAGE(FORMAT='MED', )
MAIL=DEFI_GROUP( reuse=MAIL,
MAILLAGE=MAIL,CREA_GROUP_NO=(
_F( GROUP_MA = 'TOUT_ELT', ),)
)
#
MAILLAGE=CREA_MAILLAGE( MAILLAGE=MAIL,
MODI_MAILLE=_F( TOUT = 'OUI',
OPTION = 'QUAD8_9')
)
#
MODELE=AFFE_MODELE( MAILLAGE=MAILLAGE,
AFFE=_F( TOUT = 'OUI',
PHENOMENE = 'MECANIQUE',
MODELISATION = 'COQUE_3D')
)
#
#
# CARACTERISTIQUES ELEMENTAIRES
#
#
CARA_ELT=AFFE_CARA_ELEM( MODELE=MODELE,
COQUE=_F( GROUP_MA = 'TOUT_ELT',
EPAIS = 0.00635,
# A_CIS : 2500.
COEF_RIGI_DRZ = 0.001)
)
#
#
# CARACTERISTIQUES MATERIAUX
#
#
MATERIAU=DEFI_MATERIAU( ELAS=_F( E = 3.10275E+09,
NU = 0.3)
)
#
CHAM_MAT=AFFE_MATERIAU( MAILLAGE=MAILLAGE,
AFFE=_F( TOUT = 'OUI', MATER = MATERIAU)
)
#
#
# CHARGEMENTS
#
#
COND_LIM=AFFE_CHAR_MECA( MODELE=MODELE,DDL_IMPO=(
_F( GROUP_NO = 'COTEAB',
DX = 0., DRY = 0., DRZ = 0.),
_F( GROUP_NO = 'COTECD',
DX = 0., DY = 0., DZ = 0.),
_F( GROUP_NO = 'COTEAD',
DY = 0., DRX = 0., DRZ = 0.))
)
#
CHARGE=AFFE_CHAR_MECA( MODELE=MODELE,
#FORCE_NODALE=_F( GROUP_NO = 'A',FZ = -500.)
DDL_IMPO =_F(GROUP_NO='A',DZ=-0.1 )
)
F_CHARGE=DEFI_FONCTION( NOM_PARA='INST',
PROL_DROITE='CONSTANT',
PROL_GAUCHE='EXCLU',
VALE=( 0., 0.,
1., 1.,
)
)
#
#
# RESOLUTION
#
#
LIS_INST=DEFI_LIST_REEL(DEBUT=0.0,INTERVALLE=(
_F( JUSQU_A = 10.0, NOMBRE = 100))
)
RESULTAT=STAT_NON_LINE(MODELE=MODELE,
CHAM_MATER=CHAM_MAT,
CARA_ELEM=CARA_ELT,
EXCIT=(
_F( CHARGE = COND_LIM),
_F( CHARGE = CHARGE,
TYPE_CHARGE = 'FIXE_PILO')
),
COMP_ELAS=_F( RELATION = 'ELAS',
DEFORMATION = 'GROT_GDEP',
),
INCREMENT=_F( LIST_INST = LIS_INST,),
PILOTAGE=_F( TYPE = 'LONG_ARC',
SELECTION='ANGL_INCR_DEPL',
GROUP_NO = 'TOUT_ELT',
NOM_CMP = ('DX', 'DY', 'DZ',),
ETA_PILO_MAX=1.,
COEF_MULT = 120.),
NEWTON=_F( REAC_ITER = 1,
MATRICE = 'TANGENTE'),
CONVERGENCE=_F(
ITER_GLOB_MAXI = 30,
RESI_GLOB_RELA = 1.E-6)
)
#
# calc more fields
RESULTAT=CALC_CHAMP(reuse =RESULTAT,
MODELE=MODELE,
RESULTAT=RESULTAT,
FORCE=('FORC_NODA','REAC_NODA'),
TOUT='OUI',
);
#
# IMPRESSION DES RESULTATS
#
DEPZ=RECU_FONCTION ( RESULTAT = RESULTAT,
NOM_CMP = 'DZ',
NOM_CHAM = 'DEPL',
GROUP_NO = 'A',
TOUT_ORDRE = 'OUI')
REAZ=RECU_FONCTION ( RESULTAT = RESULTAT,
NOM_CMP = 'DZ',
NOM_CHAM = 'REAC_NODA',
GROUP_NO = 'A',
TOUT_ORDRE = 'OUI')
FORZ=RECU_FONCTION ( RESULTAT = RESULTAT,
NOM_CMP = 'DZ',
NOM_CHAM = 'FORC_NODA',
GROUP_NO = 'A',
TOUT_ORDRE = 'OUI')
IMPR_FONCTION(FORMAT='TABLEAU',UNITE=82,COMMENTAIRE = '%',COMM_PARA = '%',
COURBE=_F( FONCTION = REAZ),);
IMPR_FONCTION(FORMAT='TABLEAU',UNITE=83,COMMENTAIRE = '%',COMM_PARA = '%',
COURBE=_F(FONC_X=DEPZ,FONC_Y=FORZ,),);
ETA=RECU_FONCTION (RESULTAT = RESULTAT,
NOM_PARA_RESU = 'ETA_PILOTAGE',
TOUT_ORDRE = 'OUI')
#IMPR_FONCTION (FORMAT='AGRAF', UNITE=27, UNITE_DIGR=28,COURBE=_F( FONCTION = ETA),TITRE='Eta en fonction du temps')
IMPR_FONCTION(FORMAT='TABLEAU',
UNITE=81,
COURBE=_F( FONCTION = ETA),
COMMENTAIRE = '%',
COMM_PARA = '%',
);
IMPR_RESU(FORMAT='MED',UNITE=80,RESU=_F(MAILLAGE=MAIL,RESULTAT=RESULTAT,NOM_CHAM=('DEPL','REAC_NODA','FORC_NODA'),),);
#
FIN()
#
#
#
-- CODE_ASTER -- VERSION : EXPLOITATION (stable) --
Version 11.7.0 modifi��e le 10/12/2014
r��vision 7a39a0fb6aea - branche 'v11'
Copyright EDF R&D 1991 - 2015
Ex��cution du : Tue May 19 07:57:03 2015
Nom de la machine : naxos
Architecture : 64bit
Type de processeur : x86_64
Syst��me d'exploitation : Linux CentOS 6.4 Final 2.6.32-358.23.2.el6.x86_64
Langue des messages : en (UTF-8)
Version de Python : 2.7.3
Version de NumPy : 1.8.1
Parall��lisme MPI : inactif
Parall��lisme OpenMP : actif
Nombre de processus utilis��s : 1
Version de la librairie HDF5 : 1.8.10
Version de la librairie MED : 3.0.8
Version de la librairie MUMPS : 4.10.0
Version de la librairie SCOTCH : 5.1.10
M��moire limite pour l'ex��cution : 15120.00 Mo
consomm��e par l'initialisation : 307.90 Mo
par les objets du jeu de commandes : 2.05 Mo
reste pour l'allocation dynamique : 14809.80 Mo
Taille limite des fichiers d'��change : 48.00 Go
#
# Commande No : 0001 Concept de type : -
#
DEBUT(INFO=1,
PAR_LOT='OUI',
IMPR_MACRO='NON',
MEMOIRE=_F(TAILLE_BLOC=800.0,
TAILLE_GROUP_ELEM=1000,),
RESERVE_CPU=_F(BORNE=900,),
MESURE_TEMPS=_F(MOYENNE='OUI',
NIVE_DETAIL=1,),
)
Valeur initiale du temps CPU maximum = 108900 secondes
Valeur du temps CPU maximum pass�� aux commandes = 108000 secondes
R��serve CPU pr��vue = 900 secondes
D��but de lecture...
R��ouverture de la base
Nom de la base : ELEMBASE
Cr����e avec la version : 11.07.00
Nombre d'enregistrements utilis��s : 24
Nombre d'enregistrements maximum : 512
Longueur d'enregistrement (octets) : 819200
Nombre d'identificateurs utilis��s : 123
Taille maximum du r��pertoire : 300
Pourcentage d'utilisation du r��pertoire : 41 %
Fermeture de la base
Nom de la base : ELEMBASE
Nombre d'enregistrements utilis��s : 24
Nombre d'enregistrements maximum : 512
Longueur d'enregistrement (octets) : 819200
Nombre total d'acc��s en lecture : 49
Volume des acc��s en lecture : 38.28 Mo.
Nombre total d'acc��s en ��criture : 3
Volume des acc��s en ��criture : 2.34 Mo.
Nombre d'identificateurs utilis��s : 123
Taille maximum du r��pertoire : 300
Pourcentage d'utilisation du r��pertoire : 41 %
Relecture des catalogues des ��l��ments faite.
Fin de lecture (dur��e 0.010000 s.)
# M��moire (Mo) : 337.97 / 333.50 / 27.33 / 14.95 (VmPeak / VmSize / Optimum / Minimum)
# Fin commande No : 0001 user+syst: 0.05s (syst: 0.03s, elaps: 0.06s)
#
#
# Commande No : 0002 Concept de type : maillage_sdaster
#
MAIL = LIRE_MAILLAGE(FORMAT='MED',
INFO_MED=1,
INFO=1,
UNITE=20,
VERI_MAIL=_F(APLAT=1.E-3,
VERIF='OUI',),
)
====== VERIFICATION DU MAILLAGE ======
MAILLAGE MAIL - IMPRESSIONS NIVEAU 1
CODE_ASTER - 11.01.04 - VERSION DE DEVELOPPEMENT - FRONTAL1 - L
NOMBRE DE NOEUDS 133
NOMBRE DE MAILLES 36
QUAD8 36
NOMBRE DE GROUPES DE NOEUDS 4
COTEAB 13
COTEAD 13
A 1
COTECD 13
NOMBRE DE GROUPES DE MAILLES 2
COUL_7 36
TOUT_ELT 36
# M��moire (Mo) : 337.97 / 334.63 / 27.33 / 14.95 (VmPeak / VmSize / Optimum / Minimum)
# Fin commande No : 0002 user+syst: 0.01s (syst: 0.00s, elaps: 0.01s)
#
#
# Commande No : 0003 Concept de type : maillage_sdaster
#
MAIL = DEFI_GROUP(reuse=MAIL,
MAILLAGE=MAIL,
CREA_GROUP_NO=_F(GROUP_MA=('TOUT_ELT', ),
CRIT_NOEUD='TOUS',),
ALARME='OUI',)
NOMBRE DE GROUPES DE NOEUDS CREES : 1
=======================================
! NOM DU GROUPE ! NBRE DE NOEUDS DU !
! NOEUDS ! GROUPE_NO !
! TOUT_ELT ! 133 !
# M��moire (Mo) : 337.97 / 334.63 / 27.33 / 14.95 (VmPeak / VmSize / Optimum / Minimum)
# Fin commande No : 0003 user+syst: 0.00s (syst: 0.00s, elaps: 0.00s)
#
#
# Commande No : 0004 Concept de type : maillage_sdaster
#
MAILLAGE = CREA_MAILLAGE(MAILLAGE=MAIL,
MODI_MAILLE=_F(OPTION='QUAD8_9',
TOUT='OUI',
PREF_NOEUD='NS',
PREF_NUME=1,),
INFO=1,)
MOT CLE FACTEUR "MODI_MAILLE", OCCURRENCE 1
MODIFICATION DE 36 MAILLES QUAD8 EN QUAD9
MAILLAGE MAILLAGE - IMPRESSIONS NIVEAU 1
ASTER 11.07.00 CONCEPT MAILLAGE CALCULE LE 19/05/2015 A 07:57:03 DE TYPE
MAILLAGE_SDASTER
NOMBRE DE NOEUDS 169
NOMBRE DE MAILLES 36
QUAD9 36
NOMBRE DE GROUPES DE NOEUDS 5
COTEAB 13
COTEAD 13
A 1
COTECD 13
TOUT_ELT 133
NOMBRE DE GROUPES DE MAILLES 2
COUL_7 36
TOUT_ELT 36
# M��moire (Mo) : 337.97 / 334.77 / 27.33 / 14.95 (VmPeak / VmSize / Optimum / Minimum)
# Fin commande No : 0004 user+syst: 0.01s (syst: 0.00s, elaps: 0.00s)
#
#
# Commande No : 0005 Concept de type : modele_sdaster
#
MODELE = AFFE_MODELE(MAILLAGE=MAILLAGE,
AFFE=_F(PHENOMENE='MECANIQUE',
TOUT='OUI',
MODELISATION=('COQUE_3D', ),),
INFO=1,
VERI_JACOBIEN='OUI',
PARTITION=_F(PARALLELISME='GROUP_ELEM',),
)
SUR LES 36 MAILLES DU MAILLAGE MAILLAGE
ON A DEMANDE L'AFFECTATION DE 36
ON A PU EN AFFECTER 36
MODELISATION ELEMENT FINI TYPE MAILLE NOMBRE
COQUE_3D MEC3QU9H QUAD9 36
#2 Calculs elementaires et assemblages CPU (USER+SYST/SYST/ELAPS): 0.00 0.00 0.00
# M��moire (Mo) : 341.52 / 341.50 / 27.33 / 20.91 (VmPeak / VmSize / Optimum / Minimum)
# Fin commande No : 0005 user+syst: 0.00s (syst: 0.00s, elaps: 0.01s)
#
#
# Commande No : 0006 Concept de type : cara_elem
#
CARA_ELT = AFFE_CARA_ELEM(MODELE=MODELE,
COQUE=_F(EPAIS=6.35E-3,
GROUP_MA=('TOUT_ELT', ),
COEF_RIGI_DRZ=1.E-3,
A_CIS=0.8333333,
MODI_METRIQUE='NON',
COQUE_NCOU=1,),
INFO=1,)
LE MODELE MODELE CONTIENT :
36 ELEMENT(S) COQUE
#2 Calculs elementaires et assemblages CPU (USER+SYST/SYST/ELAPS): 0.00 0.00 0.00
# M��moire (Mo) : 341.52 / 341.50 / 27.33 / 20.91 (VmPeak / VmSize / Optimum / Minimum)
# Fin commande No : 0006 user+syst: 0.01s (syst: 0.00s, elaps: 0.01s)
#
#
# Commande No : 0007 Concept de type : mater_sdaster
#
MATERIAU = DEFI_MATERIAU(ELAS=_F(E=3.102750000E9,
NU=0.3,),
INFO=1,)
MATERIAU : MATERIAU
RELATION DE COMPORTEMENT: ELAS
# M��moire (Mo) : 341.52 / 341.50 / 27.33 / 20.91 (VmPeak / VmSize / Optimum / Minimum)
# Fin commande No : 0007 user+syst: 0.00s (syst: 0.00s, elaps: 0.00s)
#
#
# Commande No : 0008 Concept de type : cham_mater
#
CHAM_MAT = AFFE_MATERIAU(MAILLAGE=MAILLAGE,
AFFE=_F(MATER=(MATERIAU, ),
TOUT='OUI',),
INFO=1,
)
# M��moire (Mo) : 341.52 / 341.50 / 27.33 / 20.91 (VmPeak / VmSize / Optimum / Minimum)
# Fin commande No : 0008 user+syst: 0.01s (syst: 0.00s, elaps: 0.01s)
#
#
# Commande No : 0009 Concept de type : char_meca
#
COND_LIM = AFFE_CHAR_MECA(MODELE=MODELE,
DDL_IMPO=(_F(DRY=0.0,
DRZ=0.0,
GROUP_NO=('COTEAB', ),
DX=0.0,),
_F(GROUP_NO=('COTECD', ),
DZ=0.0,
DX=0.0,
DY=0.0,),
_F(DRX=0.0,
DRZ=0.0,
GROUP_NO=('COTEAD', ),
DY=0.0,),
),
INFO=1,
VERI_NORM='OUI',
)
# M��moire (Mo) : 341.74 / 341.50 / 27.33 / 20.91 (VmPeak / VmSize / Optimum / Minimum)
# Fin commande No : 0009 user+syst: 0.01s (syst: 0.00s, elaps: 0.01s)
#
#
# Commande No : 0010 Concept de type : char_meca
#
CHARGE = AFFE_CHAR_MECA(MODELE=MODELE,
DDL_IMPO=_F(GROUP_NO=('A', ),
DZ=-0.1,),
INFO=1,
VERI_NORM='OUI',
)
# M��moire (Mo) : 341.75 / 341.50 / 27.33 / 20.91 (VmPeak / VmSize / Optimum / Minimum)
# Fin commande No : 0010 user+syst: 0.00s (syst: 0.00s, elaps: 0.00s)
#
#
# Commande No : 0011 Concept de type : fonction_sdaster
#
F_CHARGE = DEFI_FONCTION(PROL_DROITE='CONSTANT',
VALE=(0.0, 0.0, 1.0, 1.0),
NOM_PARA='INST',
PROL_GAUCHE='EXCLU',
INFO=1,
INTERPOL=('LIN', ),
NOM_RESU='TOUTRESU',
VERIF='CROISSANT',)
# M��moire (Mo) : 341.75 / 341.64 / 27.33 / 20.91 (VmPeak / VmSize / Optimum / Minimum)
# Fin commande No : 0011 user+syst: 0.00s (syst: 0.00s, elaps: 0.01s)
#
#
# Commande No : 0012 Concept de type : listr8_sdaster
#
LIS_INST = DEFI_LIST_REEL(DEBUT=0.0,
INTERVALLE=_F(NOMBRE=100,
JUSQU_A=10.0,),
INFO=1,)
# M��moire (Mo) : 341.75 / 341.64 / 27.33 / 20.91 (VmPeak / VmSize / Optimum / Minimum)
# Fin commande No : 0012 user+syst: 0.01s (syst: 0.00s, elaps: 0.00s)
#
#
# Commande No : 0013 Concept de type : evol_noli
#
RESULTAT = STAT_NON_LINE(CARA_ELEM=CARA_ELT,
CHAM_MATER=CHAM_MAT,
MODELE=MODELE,
EXCIT=(_F(CHARGE=COND_LIM,
TYPE_CHARGE='FIXE_CSTE',),
_F(TYPE_CHARGE='FIXE_PILO',
CHARGE=CHARGE,),
),
PILOTAGE=_F(SELECTION='ANGL_INCR_DEPL',
NOM_CMP=('DX', 'DY', 'DZ'),
GROUP_NO=('TOUT_ELT', ),
COEF_MULT=120.0,
TYPE='LONG_ARC',
ETA_PILO_MAX=1.0,
EVOL_PARA='SANS',
PROJ_BORNES='OUI',),
INCREMENT=_F(LIST_INST=LIS_INST,
PRECISION=1.E-06,),
CONVERGENCE=_F(RESI_GLOB_RELA=1.E-06,
ITER_GLOB_MAXI=30,
ARRET='OUI',
ITER_GLOB_ELAS=25,),
COMP_ELAS=_F(DEFORMATION='GROT_GDEP',
RELATION='ELAS',
RESI_INTE_RELA=1.E-06,
ITER_INTE_MAXI=10,),
NEWTON=_F(MATRICE='TANGENTE',
REAC_ITER=1,
REAC_INCR=1,
REAC_ITER_ELAS=0,),
METHODE='NEWTON',
ARCHIVAGE=_F(PRECISION=1.E-06,
CRITERE='RELATIF',),
SOLVEUR=_F(RENUM='METIS',
STOP_SINGULIER='OUI',
NPREC=8,
METHODE='MULT_FRONT',
SYME='NON',),
)
Mot-clef : COMP_ELAS - occurrence 1 : comportement ELAS - nombre de variables internes : 1
Noms des variables internes :
V1 : VIDE
Le syst��me lin��aire �� r��soudre contient 401 noeuds dont:
- 169 noeuds portant des degr��s de libert�� physiques
- 232 noeuds portant des degr��s de libert�� de Lagrange
Pour un total de 1138 ��quations.
La matrice est de taille 1138 ��quations.
Elle contient 39741 termes non nuls si elle est sym��trique et 78344 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 6.050 %.
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
Instant de calcul: 1.000000000000e-01 - Niveau de d��coupe: 1
| NEWTON | RESIDU | RESIDU | PILOTAGE | OPTION |
| ITERATION | RELATIF | ABSOLU | COEFFICIENT | ASSEMBLAGE |
| | RESI_GLOB_RELA | RESI_GLOB_MAXI | ETA | |
| 0 X | 2.56881E-01 X | 7.07797E+01 | 2.71212E-02 |TANGENTE |
| 1 X | 2.08575E-03 X | 4.24425E-01 | 2.64308E-02 |TANGENTE |
| 2 | 3.87930E-07 | 7.88830E-05 | 2.64233E-02 |TANGENTE |
Crit��re(s) de convergence atteint(s)
Le r��sidu de type <RESI_GLOB_RELA> vaut 3.879296069712e-07 au noeud et degr�� de libert�� <NO102 DY>
Temps CPU consomm�� dans ce pas de temps : 0.430 s
* Temps moyen par it��ration de Newton : 0.143 s ( 3 it��rations )
* Temps total int��gration comportement : 0.340 s ( 4 int��grations )
* Temps assemblage matrice : 0.000 s
* Temps total factorisation matrice : 0.020 s ( 3 factorisations )
* Temps construction second membre : 0.040 s
* Temps total r��solution K.U=F : 0.010 s ( 3 r��solutions )
* Temps autres op��rations : 0.020 s
M��moire (Mo) : 348.62 / 346.99 / 29.39 / 23.14 (VmPeak / VmSize / Optimum / Minimum)
Archivage des champs
Champ stock�� <DEPL> �� l'instant 1.000000000000e-01 pour le num��ro d'ordre 1
Champ stock�� <SIEF_ELGA> �� l'instant 1.000000000000e-01 pour le num��ro d'ordre 1
Champ stock�� <VARI_ELGA> �� l'instant 1.000000000000e-01 pour le num��ro d'ordre 1
Champ stock�� <COMPORTEMENT> �� l'instant 1.000000000000e-01 pour le num��ro d'ordre 1
...
EXECUTION_CODE_ASTER_EXIT_10486=0
[code]
Last edit: 8 years 11 months ago by alessandro.
Please Log in or Create an account to join the conversation.
- RichardS
- Offline
- Elite Member
-
8 years 11 months ago - 8 years 11 months ago #8102
by RichardS
SimScale - Engineering Simulation in your browser!
Replied by RichardS on topic Re: error using pilotage
It seems that you still have included nodes to the PILOITAGE that have no DY.
For COQUE_3D shells you have 9 noded QUADs where the central node doesn't have ANY displacement DOFs.
try to exlude them from the pilotage group (by creating the node group before the "quad8_9 operation, just like in the example)
Best,
richard
For COQUE_3D shells you have 9 noded QUADs where the central node doesn't have ANY displacement DOFs.
try to exlude them from the pilotage group (by creating the node group before the "quad8_9 operation, just like in the example)
Best,
richard
SimScale - Engineering Simulation in your browser!
Last edit: 8 years 11 months ago by RichardS.
Please Log in or Create an account to join the conversation.
- alessandro
- Topic Author
- Offline
- New Member
-
Less
More
- Posts: 10
- Thank you received: 0
8 years 11 months ago #8103
by alessandro
Replied by alessandro on topic Re: error using pilotage
It worked!
Thanks for the help
Thanks for the help
Please Log in or Create an account to join the conversation.
Moderators: catux
Time to create page: 0.202 seconds