TUTEF03-Axi-Symmetric analysis command file

January 21, 2014 No Comments »
TUTEF03-Axi-Symmetric analysis command file

Let us look at how to generate a command file for an Axi-symmetric Pressure Vessel analysis.

This Vertical Pressure Vessel has Internal Radius of 1000mm and a thickness of 10mm.  It has a hemispherical head on top and for the bottom part is symmetrical, so only top half of the Pressure Vessel is modeled.  Note for Axi-symmetric analysis, the axis of revolution should be Y axis.

Bottom part of shell is fixed in Y direction for symmetry and a pressure of 1MPa is applied on the inside wall of the Pressure Vessel.

Symmetric Mesh group is given a name of “FixDY” and the segment on which Pressure is applied is given a name of “Press”

The mesh is in MED format

So lets start generating command file for this analysis

Note: All screenshots are from a Windows 7 machine, using a different Operating system will not change the command file.

Step1:

Start Efficient (Click Here to see how to do that)

On the first tab, keep everything as default shown in screenshot below

TutEF02-01

Step2:

Click on Tab “Analysis Type” where we will select “Mechanical – Axi-Symmetric” in the left-top drop down box (highlighted yellow in screenshot below) and leave everything else as it is.

TutEF03-02

Step3:

Click on Tab “Material” where we will add material to the analysis.

Enter “CS” for the Material Name

Keep default Youngs Modulus of 2.1e5 which is in MPa, Poisson’s Ratio of 0.3 and Density of 7.8e-9 (tonnes/mm^3)

Enter “All” for Is assigned to, to specify that entire Beam has this material.

Click on “Add” and the material will be added to the command file

TutEF02-03

 

Step4:

For Axi-Symmetric analysis we do not need to add Element Type so skip the Tab “Element Type”

Step5:

Click on Tab “Boundary Condition(s)” where we will add information about the Symmetry Boundary Condition to the analysis.

Select “D.O.F (DDL) onMesh Group” for Type of Boundary Cond. as we are adding a boundary condition on Group of Mesh

Enter “SymmDY” for Boundary Condition Name.  This is required as this name will be used in the command file from now on.  This can be the same as the Mesh Group name defined in the mesh, but it does not need to be the same.  Make sure the length of the Boundary Condition Name is less than or equal to 8 (Limitation of Fortran).

Enter “FixDY” for Is assigned to, which is the name of the Mesh Group name defined in the mesh.  This has to be exactly the same name (case sensitive) as given in the mesh.

Enter “0” (Zero) for DY and leave everything else Blank as we want to make this node as a Symmetric boundary condition.

Click on “Add” and the boundary condition will be added to the command file.

TutEF03-05

Step6:

Click on Tab “Load(s)” where we will add information about the Pressure to be applied to the analysis.

Enter “Pressure” for Load Name,  This is required as this name will be used in the command file from now on.  This can be the same as the Node Group name defined in the mesh, but it does not need to be the same.  Make sure the length of the Load Name is less than or equal to 8 (Limitation of Fortran).

Select “Pressure” for Load Type as we will be applying Pressure on the internal wall of the Pressure Vessel.

Enter “Press” for Is Assigned to, which is the name of the Mesh Group name defined in the mesh.  This has to be exactly the same name (case sensitive) as given in the mesh.

Enter “1” for Amount of Pressure as we want to apply 1MPa Pressure.

Click on “Add” and the pressure will be added to the command file.

TutEF03-06

Step7:

Click on Tab “Output” where we will add information about the types of results that we want in the analysis.

Efficient is intelligent to know that the type of Analysis is “Axi-Symmetric” and so it showed only those output Options that are related to Axi-Symmetric analysis.

If you don’t ADD anything in this Tab, only Deflection will be saved in MED file

As we want all options that are available for this Analysis we will click “Add ALL”

TutEF03-07

Step8:

Now the only thing left is the Click “Save .comm File”, Efficient will ask for a location where you want to save this file, select a folder where you want the file to be saved and give it a name, Here I have given it a name “AxiSymmetry.comm”.  (Remember to add .comm at the end of the File name as by default this is not written by Efficient.

TutEF03-08

Click on “Save” and the file will be saved.  Once Efficient has saved the file, it will show message as follows

TutEF02-09

Hurray, you have a working command file within few seconds.

Your comm file should look like below

###############################
#File created by Efficient Software version 0.1.0
#Version of Code_Aster is 11.x
#Author of Efficient is Dhramit Thakore
#http://engineering.moonish.biz
###############################
#U4.11.01
#@Eff@#StartCont#DEBUT
DEBUT();
#U4.21.01
#@Eff@#MeshType#MED
mesh=LIRE_MAILLAGE(FORMAT='MED',);
#U4.23.04
mesh=MODI_MAILLAGE(reuse=mesh, MAILLAGE=mesh, ORIE_PEAU_2D=(_F(GROUP_MA='Press',), ),);
#U4.41.01
#@Eff@#AnalysisType#Mechanical - Axi-Symmetric
model=AFFE_MODELE(MAILLAGE=mesh, AFFE=_F(TOUT='OUI', PHENOMENE='MECANIQUE', MODELISATION='AXIS',),);
#U4.43.01
#@Eff@#MaterialList#CS;2.1e5;0.3;7.8e-9;All
CS=DEFI_MATERIAU(ELAS=_F(E=2.1e5, NU=0.3, RHO=7.8e-9,),);
#U4.43.03
material=AFFE_MATERIAU(MAILLAGE=mesh, AFFE=(_F(TOUT='OUI', MATER=CS,),),);
#U4.44.01
#@Eff@#BCList#D.O.F (DDL) on Mesh Group;SymmDY;FixDY;NA;0;NA;NA;NA;NA;NA;NA
SymmDY=AFFE_CHAR_MECA(MODELE=model, DDL_IMPO=_F(GROUP_MA='FixDY',DY=0,),);
#U4.44.01
#@Eff@#LoadList#Pressure;Pressure;Press;1;NA;NA;NA;NA;NA
Pressure=AFFE_CHAR_MECA(MODELE=model, PRES_REP=(_F(GROUP_MA='Press', PRES=1,),),);
result=MECA_STATIQUE(MODELE=model, CHAM_MATER=material, EXCIT=(_F(CHARGE=SymmDY,),_F(CHARGE=Pressure,),),);
#U4.81.04
result=CALC_CHAMP(reuse=result, RESULTAT=result, CONTRAINTE=('SIGM_ELNO','SIGM_NOEU','SIEF_ELNO',), CRITERES=('SIEQ_ELNO',), FORCE=('REAC_NODA',),);
#U4.91.01
#@Eff@#OutputMEDList#SIGM_ELNO
#@Eff@#OutputMEDList#SIGM_NOEU
#@Eff@#OutputMEDList#SIEF_ELNO
#@Eff@#OutputMEDList#SIEQ_ELNO
#@Eff@#OutputMEDList#REAC_NODA
IMPR_RESU(FORMAT='MED', UNITE=80, RESU=_F(MAILLAGE=mesh, RESULTAT=result, NOM_CHAM=('DEPL', 'SIGM_ELNO', 'SIGM_NOEU', 'SIEF_ELNO', 'SIEQ_ELNO', 'REAC_NODA',),),);
#U4.11.02
FIN();

Use this command file in one of your analysis and check results.

Related Posts

Comments are closed.