jeanpierreaubry Thanks Jean-Pierre for willing to look into the files. I am attaching the .geo file, the .comm file and the .export file. I am opening the .geo file in Gmsh and then 1D meshing the geometry and exporting the mesh in .unv format (ticking the option save groups of nodes ).
I am running standalone Code Aster and using Gmsh for meshing.
The .geo file:
" cl1 = 100;
Point(1) = {0, -1000, 0, cl1};
Point(2) = {0, -1000, 1000, cl1};
Point(3) = {0, -500, 1000, cl1};
Point(4) = {0, 0, 1000, cl1};
Line(1) = {1, 2};
Line(2) = {2, 3};
Line(3) = {3, 4};
Symmetry {0 , 1 , 0 , 0} {
Duplicata { Line{2 , 3 , 1 } ; }
}
Physical Curve("topbeam") = {2, 3, 5, 4};
Physical Curve("mast", 8) = {1, 6};
Physical Point("groundS", 9) = {1};
Physical Point("groundN", 10) = {13};
Physical Point("loadS", 11) = {3};
Physical Point("massN", 12) = {6};"
The .comm file:
"
DEBUT(LANG='EN');
mesh = LIRE_MAILLAGE(
FORMAT='IDEAS',
UNITE=3,
);
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,),
);
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',
),
),
);
steel = DEFI_MATERIAU(ELAS=F(E=210000, NU=0.3,RHO=8E+9),);
#U4.43.03
material = AFFE_MATERIAU(
MAILLAGE = mesh,
AFFE=F(GROUP_MA=('topbeam','mast'),MATER=steel,),);
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, ),
),
),
#in the next line we give to the discrete element
#the property of a point mass
#(CARA=’M_T_D_N’), and give it
#the value of 0.01 tonnes e.g. 10 kg
DISCRET = (
_F(
GROUP_MA='massN',
CARA='M_T_D_N',
VALE=(0.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',),
),
);
ground = AFFE_CHAR_MECA(
MODELE=model,
DDL_IMPO=_F(
GROUP_NO=('groundS','groundN'),
DX=0, DY=0, DZ=0,DRX=0,DRY=0,DRZ=0,
),
);
selfwght = AFFE_CHAR_MECA(
MODELE = model,
PESANTEUR = F(
GRAVITE=10000,
DIRECTION=(0,0,-1),
GROUP_MA=('topbeam', 'mast', 'massN'),
),
);
cc = AFFE_CHAR_MECA(
MODELE=model,
FORCE_NODALE=F(GROUP_NO=('loadS',),FZ=-135,)
);
cr = AFFE_CHAR_MECA(
MODELE=model,
#FORCE_POUTRE=F(GROUP_MA=(’topbeam’,),FZ=-0.1,),
FORCE_NODALE=F(GROUP_NO=('topbeam',),FZ=-0.1*2000/17,
#17 is the number of nodes in the group ’topbeam’
),
);
selfw_m = DEFI_FONCTION(
NOM_PARA='INST',
VALE=(2,0,3,1.35,5,1.35,6,0),
PROL_GAUCHE='CONSTANT',PROL_DROITE='CONSTANT',
);
cc_m = DEFI_FONCTION(
NOM_PARA='INST',
VALE=(3,0,4,1,5,1,6,0),
PROL_GAUCHE='CONSTANT',PROL_DROITE='CONSTANT',
);
cr_m = DEFI_FONCTION(
NOM_PARA='INST',
VALE=(4,0,5,1.5,6,1),
PROL_GAUCHE='CONSTANT',PROL_DROITE='CONSTANT',
);
liste = DEFI_LIST_REEL(DEBUT =2.0,
INTERVALLE=_F(JUSQU_A=6, PAS=1.0,)
,
);
stat = MECA_STATIQUE(
MODELE=model,
CHAM_MATER=material,
CARA_ELEM=elemcar,
#with the load, or boundary condition defined in EXCIT
#with the applied step function where needed
EXCIT=(
_F(CHARGE=ground,),
_F(CHARGE=selfwght,FONC_MULT=selfw_m,),
_F(CHARGE=cc,TYPE_CHARGE='FIXE_CSTE',FONC_MULT=cc_m),
_F(CHARGE=cr,TYPE_CHARGE='FIXE_CSTE',FONC_MULT=cr_m),
),
#the calculation is made for all instant in this list
LIST_INST=liste,
#we can give a title to this study
#TITRE=’my_title’
);
stat = CALC_CHAMP(
reuse=stat,RESULTAT=stat,
CONTRAINTE = ('SIEF_ELNO','SIPO_ELNO','SIPM_ELNO'),
FORCE=('REAC_NODA'),
);
masse = POST_ELEM(
RESULTAT=stat,
MODELE=model,
MASS_INER=_F(GROUP_MA=('topbeam','mast','massN',),),
);
IMPR_TABLE(TABLE=masse,
FORMAT_R='1PE12.3',)
sum_reac=POST_RELEVE_T(
ACTION=_F(
INTITULE='sum reactions',
GROUP_NO=('groundS','groundN'),
RESULTAT=stat,
NOM_CHAM='REAC_NODA',
TOUT_ORDRE='OUI',
RESULTANTE=('DX','DY','DZ'),
MOMENT=('DRX','DRY','DRZ'),
POINT=(0,0,0),
OPERATION='EXTRACTION',
),
);
IMPR_TABLE(TABLE=sum_reac,FORMAT_R='1PE12.3',)
#then in tabular format per group of node
reac1=POST_RELEVE_T(
ACTION=F(
INTITULE='reactionsS' ,
GROUP_NO=('groundS' , ) ,
RESULTAT=stat ,
NOM_CHAM='REAC_NODA' ,
TOUT_ORDRE='OUI' ,
RESULTANTE=('DX' ,'DY' ,'DZ' ) ,
OPERATION='EXTRACTION' ,
) ,
) ;
#very simple form
#IMPR_TABLE (TABLE=reac1,)
#or more detailed with field separator
IMPR_TABLE(
TABLE=reac1 ,
FORMAT='TABLEAU' ,
UNITE=8 , #this is also the default value
#whichever separator that suits the needs
SEPARATEUR='*' ,
FORMAT_R='1PE12.3' ,
TITRE='reaction_1' ,
INFO=2 ,
) ;
reac2=POST_RELEVE_T(
ACTION=F(
INTITULE='reactionsN',
GROUP_NO=('groundN' , ) ,
RESULTAT=stat ,
NOM_CHAM='REAC_NODA' ,
TOUT_ORDRE='OUI' ,
RESULTANTE=('DX' ,'DY' ,'DZ' ) ,
OPERATION='EXTRACTION' ,
) ,
) ;
IMPR_TABLE (TABLE=reac2 , FORMAT_R='1PE12.3' , )
#now we print the individual reactions
#in the .resu file in RESULTAT format
IMPR_RESU(
MODELE=model ,
FORMAT='RESULTAT' ,
RESU=_F(
NOM_CHAM='REAC_NODA' ,
GROUP_NO=('groundS' ,'groundN' , ) ,
RESULTAT=stat ,
FORMAT_R='1PE12.3' ,
) ,
) ;
IMPR_RESU(
MODELE=model ,
FORMAT='RESULTAT' ,
RESU=(
_F(
RESULTAT=stat ,
NOM_CHAM='DEPL' ,
NOM_CMP='DZ' ,
GROUP_MA=('topbeam' , ) ,
FORMAT_R='1PE12.3' ,
VALE_MAX='OUI' ,VALE_MIN='OUI' ,
) ,
_F(
RESULTAT=stat ,
NOM_CHAM='SIEF_ELNO',
NOM_CMP='N' ,
GROUP_MA=('mast' , ) ,
FORMAT_R='1PE12.3' ,
VALE_MAX='OUI' ,VALE_MIN='OUI' ,
) ,
_F(
RESULTAT=stat ,
NOM_CHAM='SIEF_ELNO' ,
NOM_CMP=('MFY' ,'MFZ' , ) ,
GROUP_MA=('topbeam' , ) ,
FORMAT_R='1PE12.3' ,
VALE_MAX='OUI' ,VALE_MIN='OUI' ,
) ,
) ,
) ;
IMPR_RESU(
MODELE=model ,
FORMAT='RESULTAT' ,
RESU=(
_F(
RESULTAT=stat ,
NOM_CHAM='SIPO_ELNO' ,
GROUP_MA=('topbeam' ,'mast' , ) ,
FORMAT_R='1PE12.3' ,
) ,
) ,
) ;
IMPR_RESU(
FORMAT='MED' ,
UNITE=80 ,
RESU=_F(
#following lines print on the named groups,
GROUP_MA=('topbeam' ,'mast' , ) ,
RESULTAT=stat ,
NOM_CHAM=(
'DEPL' ,
'SIEF_ELNO' ,
'SIPO_ELNO' ,
'SIPM_ELNO' ,
'REAC_NODA' ,
) ,
) ,
) ;
#U4.11.02
FIN ( ) ;"
and the .export file:
"P actions make_etude
P consbtc yes
P mem_aster 100
P memjob 4146176
P memory_limit 4049.0
P mode interactif
P mpi_nbcpu 1
P mpi_nbnoeud 1
P origine ASTK 2020.0
P rep_trav /tmp/simo-simo-VirtualBox-interactif_7499
P soumbtc yes
P time_limit 35996400.0
P tpsjob 599941
P version /home/.../14.6/share/aster
A memjeveux 506.125
A tpmax 35996400.0
F comm /home/.../frame1.comm D 1
F mmed /home/.../frame1.unv D 3
F mess /home/.../frame1.mess R 6
F resu /home/.../frame1.resu R 8
F rmed /home/.../frame1.rmed R 80
F mmed /home/.../frame1-v.med R 71
R base /home/.../frame1base.base RC 0"
