Welcome to the forums. Please post in English or French.

You are not logged in.

#1 Re: Code_Aster usage » Code_Aster used in thesis » 2021-12-27 13:44:47

I found it funny. Analysis of Potato with code-aster.
Title- FEM Analysis of effects of mechanical impact parameters on fruit characteristics
cigrjournal.org/index.php/Ejounral/article/view/3391/2196

#2 Code_Aster installation » Salome meca windows software throws error while running » 2021-12-17 09:54:22

nirmaljoshi
Replies: 1

I want to use windows version of code-aster. I downloaded it and after extracting the zip, I launched it.

It launches well. No problem.

Geometry and Mesh module works well. I created the geometry and mesh. No problem.

I launched Aster module and created command file. No problem.

But when I run the analysis, there is always some error such as

EXECUTION_CODE_ASTER_EXIT_XXXX=0
where XXXX changes every time I try to execute the analysis.

The .comm file is pretty simple and surely it has no error. What is going on? I need help.

Log file is copied below. MED and COMM file is attached, just in case.

I must mention that the windows pc I am trying to run is under LAN and controlled by IT guys of my office.

















Active code page: 65001


--------------------------------------------------------------------------------

Code_Aster execution


<INFO> prepare environment in C:\Users\JOSHI~1.NIR\AppData\Local\Temp\joshi_nirmalraj-PC0133-interactif_5552


--------------------------------------------------------------------------------

Copying datas


copying ...\joshi.nirmalraj-PC0133-Fri-17-173448.778\RunCase_3_Stage_1.comm...

                                                                        copying ...\joshi.nirmalraj-PC0133-Fri-17-173448.778\mesh.med...        <INFO> Parameters : memory 2048 MB - time limit 900 s


--------------------------------------------------------------------------------

Content of C:\Users\JOSHI~1.NIR\AppData\Local\Temp\joshi_nirmalraj-PC0133-interactif_5552 before execution


b''


--------------------------------------------------------------------------------

Code_Aster run


<INFO> Command line 1 :
<INFO> C:\sm-2019-w64-1.2\PRODUCTS\code_aster\14.4\bin\aster C:\sm-2019-w64-1.2\PRODUCTS\code_aster\14.4\lib\aster\Execution\E_SUPERV.py -commandes fort.1 --num_job=5552 --mode=interactif --rep_outils=C:\sm-2019-w64-1.2\PRODUCTS\code_aster\outils --rep_mat=C:\sm-2019-w64-1.2\PRODUCTS\code_aster\14.4\share\aster\materiau --rep_dex=C:\sm-2019-w64-1.2\PRODUCTS\code_aster\14.4\share\aster\datg --numthreads=1 --suivi_batch --tpmax=900.0 --memjeveux=256.0
# ------------------------------------------------------------------------------------------
EXECUTION_CODE_ASTER_EXIT_5552=0
<INFO> Code_Aster run ended, diagnostic : <F>_ABNORMAL_ABORT


--------------------------------------------------------------------------------

Content of C:\Users\JOSHI~1.NIR\AppData\Local\Temp\joshi_nirmalraj-PC0133-interactif_5552 after execution


b''


--------------------------------------------------------------------------------

Size of bases




--------------------------------------------------------------------------------

Copying results




<A>_COPYFILE       no such file or directory: fort.3


copying ...\fort.6...                                                   

<F>_ABNORMAL_ABORT Code_Aster run ended






---------------------------------------------------------------------------------

                                            cpu     system    cpu+sys    elapsed

---------------------------------------------------------------------------------

   Preparation of environment              0.00       0.00       0.00       0.00

   Copying datas                           0.00       0.00       0.00       0.06

   Code_Aster run                          0.00       0.00       0.00       0.75

   Copying results                         0.00       0.00       0.00       0.04

---------------------------------------------------------------------------------

   Total                                   0.00       0.06       0.06       1.21

---------------------------------------------------------------------------------


as_run 2018.1


------------------------------------------------------------

--- DIAGNOSTIC JOB : <F>_ABNORMAL_ABORT

------------------------------------------------------------




EXIT_CODE=4

#3 Code_Aster usage » How to generate zero thickness volume element in salome » 2021-05-10 08:25:35

nirmaljoshi
Replies: 0

Hello, I need some help.
How to generate zero thickness volume element in salome for interface element. An example is attached from the validation folder. See mesh goup G_DEG

#4 Re: Code_Aster usage » Code_Aster used in thesis » 2021-04-15 14:12:42

I used CA to analyse a dam's deflection.
"Time-Dependent Deformation of a Concrete Arch Dam in Thailand - Numerical Study on Effect of Alkali Silica Reaction on Deflection of Arch"
doi.org/10.3151/jact.19.181

#5 Re: Code_Aster usage » suction between two surface » 2020-09-23 07:00:24

You can use DEFI_FONCTION to define the force that depends on x,y,or z coordinate and apply via AFFE_CHAR_MECA_FO.

#6 Re: Code_Aster usage » Memory question related to IMPR_TABLE » 2020-09-23 06:56:08

seems you are appending the output to the same file by IMPR_TABLE

#8 Re: Code_Aster usage » Correct model for plastic material » 2020-09-23 06:51:46

Without your command file, it is difficult to know what you are upto, but ELAS_xxx refers to elastic analysis. You should perform non-linear analysis with proper material behaviour (or Comportment in CA language) to simulate the plasticity.

#9 Code_Aster usage » Which material model supports thermo mechanical analysis of concrete? » 2020-09-23 06:43:56

nirmaljoshi
Replies: 0

I checked thermo-mechanical analysis for concrete. In the model, I performed linear thermal analysis with increasing temperature. This temperature field was then used for expansion of concrete and do the non-linear mechanical analysis. I tried few models. I want to confirm the followings if I am correct. I am asking because the documentation is unclear.
1. MAZAR model can incorporate the thermal expansion
2. ENDO_ISOT_BETON can incorporate the thermal expansion
3. ENDO_ORTH_BETON Does NOT incorporate thermal expansion

Attached the file too for you to play with. It is small cube of 15 cm, fixed at one end.

#10 Re: Code_Aster usage » ERROR in documentation » 2020-05-18 06:19:14

V7.31.121

1. Dimension of mesh file provided is smaller by factor of 10.
2. rectify section 1.2 :: sigma_y=6E6 Pa; ET=−6E5 Pa
3. rectify section 3.2 :: number of nodes=128

#11 Code_Aster usage » SCHEMA_THM in STAT_NON_LINE gives error » 2020-05-18 05:55:10

nirmaljoshi
Replies: 1

The command below is taken from validation file "wtnv121a.comm". The issue is in the line "SCHEMA_THM=_F(PARM_THETA=0.57,)".

---------
REPONSE=STAT_NON_LINE( MODELE     = MO,
                       CHAM_MATER = CHMAT,
                       EXCIT      = ( _F( CHARGE = CLIM),
                                      _F( CHARGE = CHARGT, FONC_MULT=CAP)),
                       SCHEMA_THM=_F(PARM_THETA=0.57,),
                       COMPORTEMENT  = _F( RELATION     = 'KIT_HHM',
                                        DEFORMATION  = 'PETIT',
                                        RELATION_KIT = ( 'ENDO_ISOT_BETON',
                                                         'LIQU_GAZ',
                                                         'HYDR_ENDO',      ),),
                      ETAT_INIT   = _F( SIGM      = SIGINIT,
                                        DEPL      = DEPINIT),
                      INCREMENT   = _F( INST_FIN  = 1.0,
                                        LIST_INST = L_INST,),
                      NEWTON      = _F( MATRICE   = 'TANGENTE',
                                        REAC_ITER = 1     ),
                      CONVERGENCE = _F( RESI_GLOB_RELA = 1.E-6,
                                        ITER_GLOB_MAXI = 20),
                      SOLVEUR=_F(METHODE='MULT_FRONT',RENUM='MDA'),                                   
                     )
--------
The program works fine if I remove this line. But when this command is inserted, it gives error saying
      !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
      ! Unknown key words: SCHEMA_THM !
      !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

U4.51.3 section "3.7 Keyword SCHEMA_THM" shows that it can be used. So why there is this error?

#12 Re: Code_Aster usage » mechanical non linear analysis with time varying load » 2020-04-28 16:48:59

thanks for the reply. May be i could not make it clear in my question. The solution you proposed gives the load that varies with time. This is fine. This has already been defined in my problem. My question was : now at each of these time steps, i need to perform non linear analysis with for some material behaviour(say Mazer). In other words, i need to perform incremental loading again (from zero to load value) at each time step to capture material nonlinearity. How to do it?
Or am i having some mis conception?

#13 Code_Aster usage » mechanical non linear analysis with time varying load » 2020-04-28 10:45:57

nirmaljoshi
Replies: 5

Hello everyone,
I request for help on following.

I have a mechanical model where the load depends on time (Time1, see below). How to perform non-linear analysis such that, the load in specific time instant (say Time1=5) is used as the basic load and incremental loading is done for this load (say in the increment of 0.1). Tentatively shown in the false code below:

Time1=0 to N step 1
   someLoad=functionOf(Time1)
   Time2=0 to 1 step 0.1
       non_linear_analysis(someLoad in steping of Time2, someComportment)
   loop Time2
loop Time1

Or should i be doing non-linear analysis for each step separately?

#14 Re: Code_Aster usage » Using CREA_CHAMP to assign stress values to individual Gauss points » 2020-04-24 06:58:41

You do not have to worry about element indices. CREA_CHAMP will do it for you. Basically every formula application needs three steps:
a. apply formula (OPERATION=AFFE)
b. evaluate it (OPERATION=EVAL) (again you may need secondary CREA_CHAMP to put values inside the function of (a))
c. transform the evaluated values to gauss point or nodal points etc (OPERATION=ASSE)

transformation to node and gauss points can be done with OPERATION=DISC (discritize)

Again, please refer to zzzz130a, it shows all the possibilities.

#15 Re: Code_Aster usage » Using CREA_CHAMP to assign stress values to individual Gauss points » 2020-04-23 15:30:51

I am not sure about the python. But way i know is either of following
a)You need to either read it from external file.
b) do all calculation in code-aster using the Crea_cham using formulas; this allows to put all values inside the code aster.

See test file zzzz130a.

#16 Re: Code_Aster usage » reading and applying field (of strain) from table to model » 2020-04-18 14:30:22

I revised the code with new table format (shown below),but it still gives error saying "the node ?????? does not form part of the mesh". I checkd the mesh "mesh0", there exists the node "N1". What am i doing wrong?

table format
---
MESH NOEUD SIXX SIYY SIZZ
mesh0 N1 -2.62931E-06  3.24443E-06 -9.16076E-06
...
------

#17 Code_Aster usage » reading and applying field (of strain) from table to model » 2020-04-18 05:34:59

nirmaljoshi
Replies: 1

Hi,
I have a table of EPSI_NOEU as follows :
NOEUD    EPXX    EPYY    EPZZ    EPXY    EPXZ    EPYZ
N1    7.200E-06    7.200E-06    7.200E-06    0.000E+00    0.000E+00    0.000E+00
N2    6.310E-24    6.310E-24    6.310E-24    0.000E+00    0.000E+00    0.000E+00
N3    6.310E-24    6.310E-24    6.310E-24    0.000E+00    0.000E+00    0.000E+00

I could read it using
A_jEPSI = LIRE_TABLE(
  FORMAT='TABLEAU',
  UNITE=3
)

I want to apply it to the mesh as a field.
jepsi = CREA_CHAMP(
  INFO=2,
  MAILLAGE=mesh0,
  MODELE=M_Mod,
  OPERATION='EXTR',
  TABLE=A_jEPSI,
  TYPE_CHAM='NOEU_TEMP_R'
)
At this point, it gives error (shown below). I do not understand what is the exact format to import this field? Or am i doing something wrong (U4.72.04, section 3.9.7.4)?
------------------------------------------------------------------------
User error in command CREA_CHAMP / EXTR / COUNT:           !
   !  In the array A_jEPSI, to create a field of the type NOEU, !
   !  certain parameters                                        !
   ! are compulsory and others are prohibited:                  !
   !                                                            !
   !  NOEU:                                                     !
   !  compulsory: Prohibited                                    !
   !  NODE: NET, POINT, SOUS_POINT                              !
   !  ELGA:                                                     !
   !  compulsory: NET,                                          !
   ! prohibited                                                 !
   !  POINT: NODE                                               !
   !  ELNO:                                                     !
   !  compulsory: NET, NODE (or POINT)                          !
   !  ELEM:                                                     !
   !  compulsory: NET                                           !
   !  prohibited: NODE, POINT 
------------------------------

#18 Code_Aster usage » Creating field of time for multisteps » 2020-04-16 12:50:24

nirmaljoshi
Replies: 0

Hello everyone,
I have an application where the strain depends on temperature and time (e.g. field=INST*TEMP).  There is an example in  U4.72.04 -Section 3.6.3 with almost same criteria (see below, notice CHINSTG-time step field).

RHO=1000.
G=10.
KP=3.
SIZZ = FORMULA (NOM_PARA = ‘Z’, VALE = ‘RHO*G*Z’)
SIXX = FORMULA (NOM_PARA = (‘Z’, ‘INST’), VALE = ‘KP*SIZZ (Z) +3.*INST’)
# Assignment of the functions:
# -----------------------------------------------
SIG1=CRÉA_CHAMP (OPERATION=' AFFE', TYPE_CHAM=' ELGA_NEUT_F',
MODELE=MO, PROL_ZERO=' OUI',
AFFE=_F (TOUT=' OUI', NOM_CMP= (‘X1’, ‘X2’),
VALE_F= (‘SIXX’, ‘SIZZ’)))
# Evaluation of the functions:
# -----------------------------------------------
SIG2= CREA_CHAMP (OPERATION=' EVAL', TYPE_CHAM=' ELGA_NEUT_R',
MODELE=MO, CHAM_F=SIG1, CHAM_PARA= (CHGEOMG, CHINSTG) )
# transformation of the field of NEUT_R into SIEF_R:
# -----------------------------------------------
SIG3=CRÉA_CHAMP (OPERATION=' ASSE', TYPE_CHAM=' ELGA_SIEF_R',
MODELE=MO, PROL_ZERO=' OUI',
ASSE=_F (TOUT=' OUI', CHAM_GD=SIG2,
NOM_CMP= (‘X1’ , ‘X2’),

The example runs smoothly.

My problem is , I have to do same thing for  timesteps t=0 to 100 @step size=1. Is there way in code_Aster alone to assign (calculate) the time filed for all steps in a single command?

#19 Re: Code_Aster usage » (SOLVED) Force Nodale in DKT and COQUE_3D » 2020-04-16 01:50:21

Could you please mark it solved in the subject line too (edit the first post). Thanks!

#20 Re: Code_Aster usage » [solved]applying time and space varying function as boundary condition » 2020-04-16 01:44:33

Thanks for the reply.

P_func_Z = FORMULE(NOM_PARA=('INST','Z',),
                   VALE='INST * Z * 0.1');

loadWatF = AFFE_CHAR_MECA_F(MODELE=model,
                            PRES_REP=_F(GROUP_MA=('Group_loadWater260', ),
                                        PRES=P_func_Z))

With the code shown above, the pressure "P_func_Z" will be applied thoughout the "Group_loadWater260" from Z=0 to Z=HeightOfDam. This is not the case. When water level rises (with time), the pressure should be applied upto the water level at that instant only, any other part should have pressure zero.

As  I mentioned earlier, at every time step z1, z2 and p1, p2 will change (see code below); so, the forumla mentioned above will not be correct.

P_func_Z = DEFI_FONCTION(NOM_PARA='Z',
                         VALE=(z1, p1, z2, p2))

#21 Re: Code_Aster usage » [solved]applying time and space varying function as boundary condition » 2020-04-15 10:42:02

That means some python will come in scene. Ok, I will test it.

For now, I did work around by putting the function values in text file:
----------
Time
0.00
1.00
2.00
...

Z    P2    P3    P4    P5    P6    P7    P8    P9    P10    P11
3.7E+01    1.5E+06    2.9E+06    4.4E+06    5.9E+06    7.4E+06    8.8E+06    1.0E+07    1.2E+07    1.3E+07    1.5E+07
5.2E+01    0.0E+00    1.5E+06    2.9E+06    4.4E+06    5.9E+06    7.4E+06    8.8E+06    1.0E+07    1.2E+07    1.3E+07
6.7E+01    0.0E+00    0.0E+00    1.5E+06    2.9E+06    4.4E+06    5.9E+06    7.4E+06    8.8E+06    1.0E+07    1.2E+07
8.2E+01    0.0E+00    0.0E+00    0.0E+00    1.5E+06    2.9E+06    4.4E+06    5.9E+06    7.4E+06    8.8E+06    1.0E+07
...
-------------------


Then I read it using LIRE_FONCTION:
-----
P_func_Z = LIRE_FONCTION(
  DEFI_FONCTION=(_F(
      INDIC_RESU=(2, 2)
    ), _F(
      INDIC_RESU=(2, 3)
    ), _F(
      INDIC_RESU=(2, 4)
    ), _F(
      INDIC_RESU=(2, 5)
    ), _F(
      INDIC_RESU=(2, 6)
    ), _F(
      INDIC_RESU=(2, 7)
    ), _F(
      INDIC_RESU=(2, 8)
    ), _F(
      INDIC_RESU=(2, 9)
    ), _F(
      INDIC_RESU=(2, 10)
    ), _F(
      INDIC_RESU=(2, 11)
    )),
  FORMAT='LIBRE',
  INDIC_ABSCISSE=(2, 1),
  INDIC_PARA=(1, 1),
  NOM_PARA='INST',
  NOM_PARA_FONC='Z',
  TYPE='NAPPE',
  UNITE=3
)
------

I hope it works. I will test your method and post again.

#22 Re: Code_Aster usage » Modal Analysis: How to print the base? » 2020-04-15 02:56:13

If you are using Parvis for post processing, you can have two results separately in pipe-line browser. Then you can compare the results.

#23 Re: Code_Aster usage » (SOLVED) Force Nodale in DKT and COQUE_3D » 2020-04-15 02:53:41

Posting your file could be more helpful to solve it.
Just stating "too high deflection" is meaningless unless we know other parameters.

Generally, you may need to refine mesh. Force_nodale should work fine.

#25 Re: Code_Aster usage » [solved] how to apply load in curved section normal to each element » 2020-04-14 07:31:36

Here is the verification with reaction forces. For displacement, needs finer mesh.