no such file or directory: fort.80

More
8 years 2 months ago #8563 by F
Dear all,
I try to run my first tutorial with Salome+Code Aster.
Unfortunatly the simulation does not start and I get following error:

<A> COPYFILE: No such file or directory fort.80

Where can lie the error?

Here the .comm file

U4.11.02
DEBUT();

#U4.21.01
mesh=LIRE_MAILLAGE(UNITE=20,
FORMAT='MED',
INFO=1,);

#U4.22.01
mesh=DEFI_GROUP(reuse =mesh,
MAILLAGE=mesh,
CREA_GROUP_MA=_F(NOM='TOUT',
TOUT='OUI',),
CREA_GROUP_NO=_F(TOUT_GROUP_MA='OUI',),);

IMPR_RESU(FORMAT='MED',
UNITE=71,
RESU=_F(MAILLAGE=mesh,),);

#U4.41.01
model=AFFE_MODELE(MAILLAGE=mesh,
AFFE=(_F(GROUP_MA=('topbeam','mast',),
PHENOMENE='MECANIQUE',
MODELISATION='POU_D_T',),
_F(GROUP_MA='massN',
PHENOMENE='MECANIQUE',
MODELISATION='DIS_T',),),);

#F MATERIALI ##########################################################################################

#U4.43.01
steel=DEFI_MATERIAU(ELAS=_F(E=210000.0,
NU=0.3,
RHO=8e-09,),);

#U4.43.03
material=AFFE_MATERIAU(MAILLAGE=mesh,
AFFE=_F(GROUP_MA=('topbeam','mast',),
MATER=steel,),);

#U4.42.01
elemcar=AFFE_CARA_ELEM (
MODELE=model ,
POUTRE = (
#the vertical members are rectangular section
#(40x20 mm) with a thickness of 1.5 mm
_F (
GROUP_MA = ( 'mast' , ) ,
SECTION='RECTANGLE' ,
CARA = ( 'HY' , 'HZ' , 'EPY' , 'EPZ', ) ,
VALE = ( 40 , 20 , 1.5 , 1.5, ) ,
),
#same with the horizontal bar
_F (
GROUP_MA = ( 'topbeam' , ) , SECTION='RECTANGLE' ,
CARA = ( 'HY' , 'HZ' , 'EPY' , 'EPZ', ) , VALE = ( 40 , 20 , 1.5 , 1.5, ) ,
),
),
DISCRET = (
_F ( GROUP_MA='massN' , CARA='M_T_D_N' , VALE = ( .01 ) , ) ,
#following block set stiffness of point element 'massN'
#to null stiffness
#although this is not necessary,
#commenting this block would raise a warning,
#unimportant in this case
_F (
GROUP_MA = ( 'massN' , ) ,
CARA='K_T_D_N' ,
VALE = ( 0 , 0 , 0 , ) ,
REPERE='GLOBAL' ,
),
),
);


#F BOUNDARY CONDITIONS ################################################################################

#F incastro ai nodi di terra
#U4.44.01
ground=AFFE_CHAR_MECA(MODELE=model,
DDL_IMPO=_F(GROUP_NO=('groundS','groundN',),
DX=0.0,
DY=0.0,
DZ=0.0,
DRX=0.0,
DRY=0.0,
DRZ=0.0,),);

#F LOADS ##############################################################################################

#F gravità (PESANTEUR)
selfwght=AFFE_CHAR_MECA(MODELE=model,
PESANTEUR=_F(GROUP_MA=('topbeam','mast','massN',),
GRAVITE=10000.0,
DIRECTION=(0.0,0.0,-1.0,),),);

#F forza in un punto della struttura ('loadS')
cc=AFFE_CHAR_MECA(MODELE=model,
FORCE_NODALE=_F(GROUP_NO='loadS',
FZ=-135.0,),);

#F forza distribuita ('topbeam') (ASTER non accetta due forze distribuite (gravità più questa)
# quindi la azssegno nodo per nodo)
cr=AFFE_CHAR_MECA(MODELE=model,
FORCE_NODALE=_F(GROUP_NO='topbeam',
FZ=-11.7647058824,),);

#F stepping di applicazione delle forze
selfw_m=DEFI_FONCTION(
NOM_PARA='INST',
VALE=(2.0 ,0.0 ,
3.0 ,1.35 ,
5.0 ,1.35 ,
6.0 ,0.0 ,),
PROL_DROITE='CONSTANT',
PROL_GAUCHE='CONSTANT',);

cc_m=DEFI_FONCTION(
NOM_PARA='INST',
VALE=(3.0 ,0.0 ,
4.0 ,1.0 ,
5.0 ,1.0 ,
6.0 ,0.0 ,),
PROL_DROITE='CONSTANT',
PROL_GAUCHE='CONSTANT',);

cr_m=DEFI_FONCTION(
NOM_PARA='INST',
VALE=(4.0 ,0.0 ,
5.0 ,1.5 ,
6.0 ,1.0 ,),
PROL_DROITE='CONSTANT',
PROL_GAUCHE='CONSTANT',);

# SOLUZIONE ###########################################################################################

#F stepping the soluzione
liste=DEFI_LIST_REEL(DEBUT=2.0,
INTERVALLE=_F(JUSQU_A=6.0,
PAS=1.0,),);

#U4.51.01

stat=MECA_STATIQUE(MODELE=model,
CHAM_MATER=material,
CARA_ELEM=elemcar,
EXCIT=(_F(CHARGE=ground,),
_F(CHARGE=selfwght,
FONC_MULT=selfw_m,),
_F(CHARGE=cc,
FONC_MULT=cc_m,
TYPE_CHARGE='FIXE',),
_F(CHARGE=cr,
FONC_MULT=cr_m,
TYPE_CHARGE='FIXE',),),
LIST_INST=liste,);

# RISULTATI ###########################################################################################
#U4.80.01

stat=CALC_CHAMP(reuse =stat,
RESULTAT=stat,
CONTRAINTE=('SIEF_ELNO','SIPO_ELNO','SIPM_ELNO',),
FORCE='REAC_NODA',);

#F calcola le masse e le restituisce in una tabella

masse=POST_ELEM(MASS_INER=_F(GROUP_MA=('topbeam','mast','massN',),),
MODELE=model,
RESULTAT=stat,
TITRE='masse',);

#U4.91.03

IMPR_TABLE(TABLE=masse,
FORMAT_R='1PE12.3',);

#F calcola le reazioni vincolari e le restituisce in una tabella
#U4.81.21

sum_reac=POST_RELEVE_T(ACTION=_F(OPERATION='EXTRACTION',
INTITULE='sum reactions',
RESULTAT=stat,
NOM_CHAM='REAC_NODA',
TOUT_ORDRE='OUI',
GROUP_NO=('groundS','groundN',),
RESULTANTE=('DX','DY','DZ',),
MOMENT=('DRX','DRY','DRZ',),
POINT=(0.0,0.0,0.0,),),);

IMPR_TABLE(TABLE=sum_reac,
FORMAT_R='1PE12.3',);

#F calcola le reazioni vincolari e le restituisce in una tabella
#then in tabular format per group of node

reac1=POST_RELEVE_T(ACTION=_F(OPERATION='EXTRACTION',
INTITULE='reactionsS',
RESULTAT=stat,
NOM_CHAM='REAC_NODA',
TOUT_ORDRE='OUI',
GROUP_NO='groundS',
RESULTANTE=('DX','DY','DZ',),),);

#very simple form
#IMPR_TABLE (TABLE=reac1,)
#or more detailed with field separator

IMPR_TABLE(TABLE=reac1,
FORMAT='TABLEAU',
UNITE=8,
FORMAT_R='1PE12.3',
SEPARATEUR=' * ',
TITRE='reaction_1',
INFO=2,);

reac2=POST_RELEVE_T(ACTION=_F(OPERATION='EXTRACTION',
INTITULE='reactionsN',
RESULTAT=stat,
NOM_CHAM='REAC_NODA',
TOUT_ORDRE='OUI',
GROUP_NO='groundN',
RESULTANTE=('DX','DY','DZ',),),);

IMPR_TABLE(TABLE=reac2,
FORMAT_R='1PE12.3',);

#now we print the individual reactions
#in the .resu file in RESULTAT format
#U4.91.01

IMPR_RESU(FORMAT='RESULTAT',
MODELE=model,
RESU=_F(RESULTAT=stat,
NOM_CHAM='REAC_NODA',
GROUP_NO=('groundS','groundN',),
FORMAT_R='1PE12.3',),);

IMPR_RESU(FORMAT='RESULTAT',
MODELE=model,
RESU=(_F(RESULTAT=stat,
NOM_CHAM='DEPL',
NOM_CMP='DZ',
GROUP_MA='topbeam',
VALE_MAX='OUI',
VALE_MIN='OUI',
FORMAT_R='1PE12.3',),
_F(RESULTAT=stat,
NOM_CHAM='SIEF_ELNO',
NOM_CMP='N',
GROUP_MA='mast',
VALE_MAX='OUI',
VALE_MIN='OUI',
FORMAT_R='1PE12.3',),
_F(RESULTAT=stat,
NOM_CHAM='SIEF_ELNO',
NOM_CMP=('MFY','MFZ',),
GROUP_MA='topbeam',
VALE_MAX='OUI',
VALE_MIN='OUI',
FORMAT_R='1PE12.3',),),);

IMPR_RESU(FORMAT='RESULTAT',
MODELE=model,
RESU=_F(RESULTAT=stat,
NOM_CHAM='SIPO_ELNO',
GROUP_MA=('topbeam','mast',),
FORMAT_R='1PE12.3',),);

# stampa risultati per via grafica
#U7.05.21
IMPR_RESU(FORMAT='MED',
UNITE=80,
RESU=_F(RESULTAT=stat,
NOM_CHAM=('DEPL','SIEF_ELNO','SIPO_ELNO','SIPM_ELNO','REAC_NODA',),
GROUP_MA=('topbeam','mast',),),);

#U4.11.02
FIN();
#CHECKSUM:6dd3acf7402fc55ad65faab30775c6fd -:FIN CHECKSUM[/color]

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

More
8 years 2 months ago #8564 by kwou
Replied by kwou on topic Re: no such file or directory: fort.80
hoi Felice

<A> COPYFILE: No such file or directory fort.80

probably means that the result file (coupled to fort.80) has not been created.
This may be caused by anything, so the mess file is more accurate.


groeten

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.

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