Posts by massimobrivio

massimobrivio · Code_Aster usage

Good morning,

I'm using DKT elements and from the documentation i thought that transverse shear was not considered. But when i request the results i see that the values of Qx and Qy are not 0 (in my case they are even big).
Moreover the energy 'CISAILLE' (extracted with ENEL_ELNO) which i assume is the transverse shear, is non zero.
Someone knows what is going on? I probably misunderstand something.

Massimo

massimobrivio · Salome-Meca usage

I'm working on a semicircular geometry and orthotropic material properties using shell DST elements. I need to set the reference direction for each element in the radial direction so that the main orthotropic direction is pointing radially. Is there a way of quickly doing this from the Salome interface?
I attached a picture of the geometry for clarity.
Thanks,
Massimo

massimobrivio · Salome-Meca usage

Good Afternoon,

I was wondering if there is a way of getting the field of displacement calulated by MECA_STATIQUE as a field of type CHAM_ELEM. I need to get the elements' names in a table together with the results of the displacement on every node. For fields calculated in CALC_CHAMP i can do this in POST_RELEVE_T by selecting MOYE_NOEUD = 'NON', but this option only works with fields of type CHAM_ELEM as the documentation U4.81.21 says.
Do you know a way?

Best,
Massimo

massimobrivio · Code_Aster development

Good Morning,

I would need to extract the stiffness matrix relative to every element of the FE model. I had a look to documentation U1.03.02 in particular the section "Access method to an unspecified structure of data". But onestly i don't really understand how this works.
I also ad a look to the test case sdll123c but again is really not well explained.
What are the basics to access these objects. Which steps should i follow? How do i access the fortnan routines? How do i know which routine does what?
I assume that the routine where the total stiffness matrix for the structure is assembled also contains the elemental stiffness matrixes that i need. How to get them?

Regards,
Massimo

massimobrivio · Code_Aster usage

Thanks.
E-mail sent.

Massimo

massimobrivio · Code_Aster usage

hi Mario! Thanks a lot for the indications.
And in paraview what does it happen? Do you know?
Can i contact you directly somehow?

Massimo

massimobrivio · Code_Aster usage

Sorry for resuming an old post.
How did you get through this issue. I'm facing the exact same issue (i'm not running from flesheur tho).
I cannot figure it out.

Massimo

massimobrivio · Salome-Meca usage

Thanks Johannes,

I will have a look soon and leave a feedback.

Massimo

massimobrivio · Code_Aster usage

Thanks Mario,

I solved. I read in another post that aster change the element ID with the respect of the ones created in the mesh file. But i checked and it's not the case in salome 2019. So i can directly read them out from the impr_table generated file.

Thanks for the help,
Massimo

massimobrivio · Code_Aster usage

Good Afternoon,

I'm running the command as_run inside a python loop because i want to do some optimization over a structure. The first iteration works fine but at the second one i get the following error:

2020/10/29-15:58:47.429 [Mess@/home/massimo/salome_meca/V2019.0.3_universal/tools/Code_aster_frontend-20190/lib/python3.6/site-packages/asrun/run.py:639]
Warning or error raised :
<A>_COPY_DATA '/home/massimo/_simulazioni/_13_notched_element_ext/notched_laminate.comm' overwrites 'fort.1.1'
Traceback:
File "<string>", line 1, in <module>
File "/home/massimo/salome_meca/V2019.0.3_universal/tools/Code_aster_frontend-20190/lib/python3.6/site-packages/asrun/main.py", line 109, in main
start()
File "/home/massimo/salome_meca/V2019.0.3_universal/tools/Code_aster_frontend-20190/lib/python3.6/site-packages/asrun/main.py", line 92, in start
meth(run, *args)
File "/home/massimo/salome_meca/V2019.0.3_universal/tools/Code_aster_frontend-20190/lib/python3.6/site-packages/asrun/execute.py", line 434, in RunAster
copyfiles(run, 'DATA', prof)
File "/home/massimo/salome_meca/V2019.0.3_universal/tools/Code_aster_frontend-20190/lib/python3.6/site-packages/asrun/execution.py", line 532, in copyfiles
icomm = copyfileD(run, df, icomm, ncomm)
File "/home/massimo/salome_meca/V2019.0.3_universal/tools/Code_aster_frontend-20190/lib/python3.6/site-packages/asrun/execution.py", line 556, in copyfileD
run.Mess(ufmt(_("'%s' overwrites '%s'"), df['path'], dest), '<A>_COPY_DATA')
File "/home/massimo/salome_meca/V2019.0.3_universal/tools/Code_aster_frontend-20190/lib/python3.6/site-packages/asrun/run.py", line 639, in Mess
print_traceback=True)
File "/home/massimo/salome_meca/V2019.0.3_universal/tools/Code_aster_frontend-20190/lib/python3.6/site-packages/asrun/run.py", line 817, in DBG
stack_id=stack_id, all=all, prefix=prefix, *args)
File "/home/massimo/salome_meca/V2019.0.3_universal/tools/Code_aster_frontend-20190/lib/python3.6/site-packages/asrun/run.py", line 773, in _printDBG
stack = traceback.format_stack(limit=10)

<A>_COPY_DATA '/home/massimo/_simulazioni/_13_notched_element_ext/notched_laminate.comm' overwrites 'fort.1.1'

.
.
.
.
.

<F>_COPY_ERROR error during copying fort.1 to fort.1
message : cp: 'fort.1' and 'fort.1' are the same file

.
.
.

as_run 2018.1

------------------------------------------------------------
--- DIAGNOSTIC JOB : <F>_COPY_ERROR
------------------------------------------------------------

EXIT_CODE=4

The command i use in the python loop is: os.system(PATH_to_salome + 'salome shell as_run -- ' + exportfile.export). As i mentioned this command works for single jobs but not in a loop, it seems.
I guess that it fails to overwrite the files of the old run (1st iteration) and to write the new files from the new iteration.

How can i solve the error?

Have a good day,
Massimo

massimobrivio · Code_Aster usage

Good Evening,

Does someone know how i can recover the element IDs of elements inside a group?
I define some groups into the mesh file and then i assign properties to each of these groups in asterstudy. I would like to change those properties based on the results (on each group) that i obtain from the analysis. I can't manage to do that because in the table that i create as a result (POST_RELEVE_T) i have the results listed for each node (and also element) of the mesh but i can't identify which element is part of the groups i defined in the mesh.

Regards,
Massimo

massimobrivio · Code_Aster usage

Hi, sameer!

Right now since is for my thesis i can't really do it. But at the end of the thesis when i have more clear ideas i can definetly pass on what i have learnt doing a tutorial or so.

Massimo

massimobrivio · Code_Aster usage

Hello Johannes,

Unfortunately i cannot find anymore the command EXTR_MATR() in salome_meca 2019.0.3 GUI.
Can i still use it in the command file?

Massimo

massimobrivio · Salome-Meca usage

Good morning,

thanks for the reply. I'm working with orthotropy and layered shell elements (DEFI_COMPOSITE).
Two things:
1) i wrote by mistake AFFE_MATERIAU() but orientation is done in AFFE_CARA_ELEM()
2) by now i was doing orientation by mean of the keyword VECTEUR or ANGLE_REP in COQUE = _F(VECTEUR ())

Is the the any difference with the MASSIF keyword you said? Since in the documentation i read that also MASSIF = _F() is used to "assign to elements 3D or 2D local axes (which can be for example used to define directions
of orthotropism (cf. DEFI_MATERIAU [U4.43.01], DEFI_COMPOR [U4.43.06]."

Moreover, in the documentation you provided u1.03.02 there is a reference to EXT_MATR() but the command is not present anymore in salome meca V2019.0.3. Is there an equivalent?

Have nice day,
Massimo

massimobrivio · Code_Aster usage

Hello,

I was wondering wether is possible to extract in python the stiffness matrix of each element (element stiffness matrixes that then are used to assemble the global stiffness matrix). I read the command CALC_MATR_ELEM() documentation and it calculate what i want, the only thing is that i would like to print it to a file or convert it into a python object. It seems there was such a command EXTR_MATRIX() or IMPR_MATRIX() but both of them are not available anymore in my salome version 2019.0.3.

How can i do it?

Massimo

massimobrivio · Code_Aster usage

How? i can't manage to find anything to click to mark it as solved.

massimobrivio · Code_Aster usage

Ok I solved. Thanks Fernando for the help. I managed to run a salome meca simulation in YACS and integrate it with other python files. Copy pasting the OpenTurns function in a inline python function block in YACS does the trick, as simple as that.

massimobrivio · Salome-Meca usage

Good morning,

I'm have to assign different material orientation to each single element of a mesh. the orientation is calculated in python by mean of an algorithm.
I'm not sure how this could be possible.
something like

fieldmat[j] = AFFE_MATERIAU(...
VECTEUR='vector[j]',
GROUP_MA='element[j]')

where "element" contains a list of all the elements' IDs of the mesh and "vector" a list of directions.
My question is on how to get the list "element".

Thanks in advance,
Massimo

massimobrivio · Code_Aster usage

Thanks Ioannis, i solved by trying as you suggested. The stresses are always calculated in the user reference mark. But you can contruct a field with the given rotations in CALC_CHAMP().

Thanks
massimo

massimobrivio · Code_Aster usage

Good Morning,

I read on documentation U2.02.01 and U2.01.05 that stresses in elements of type coque are calculated in the user reference mark defined in AFFE_CARA_ELEM specifying VECTEUR or ANGL_REP. But when deinfing a layered shell with DEFI_COMPOSITE, i also specify the angle of orientation of each layer. When extracting the stresses at the layers with EXTR_COQUE in POST_CHAMP the extracted stresses are expressed in the reference mark defined in AFFE_CARA_ELEM or they are already considering the orientation of the layer in consideration?

Have a good day,
Massimo

massimobrivio · Code_Aster usage

Hi Fernando,

many thanks for the detailed reply. I just watched the tutorial and it seems this strategy of using OT to drive CA can be what i'm searching for. My aim would be to wrap a simulation into a custo-maid optimization algorithm. I saw that into OT optimization module can be used with several algorthms already available. Would in case be possible to include an equivalent python script as the one you copy pasted into the optimization routine that i created in python?

Massimo

massimobrivio · Code_Aster usage

Thanks a lot Irvise. This is the strategy i was following. But I was wondering wheter is possible to use POST_RELEVE_T and output its results in my python environment (spyder), without having to export the table in a .dat file and then import it in spyder for the calculations. This sounds to me a bit unefficient whether would be nice that the values extracted through POST_RELEVE_T would already be stored into a python variable in my environement after the simulation ends.

max

massimobrivio · Code_Aster usage

Good Morning community,

I'm running a simulation and extracting some results into a table. My purpose is to be able to read this table into a python main file which define some parameters, run the simulation and then get the results ready to analyze. I would like to use something like Spyder or Jupyter notebook, since i find it very convenient for post processing and manipulating data.

By now i'm able to set parameters and run the simulation from the "main file" in spyder, but i cannot manage to create, in the main file, a python variable containing the results of the table i created.
I used the command EXTRACT_TABLE() at the end of the command file but, how do i access the variable created from my main file in spyder?

Attached here the command file.

Best,
Massimo

massimobrivio · Introduce yourself / Présentez vous

Hello everyone, happy to be now part of the code_aster community. I'm a Master student in mechanical engineering, writing now his master thesis using code_aster.

Massimo

#1 Code_Aster usage » line loads QX, QY when using DKT elements » 2021-02-23 11:50:16

#2 Salome-Meca usage » radial projections for local reference systems (element orientation) » 2021-01-21 15:51:18

#3 Salome-Meca usage » DEPL field as CHAM_ELEM » 2020-12-09 15:36:15

#4 Code_Aster development » Extracting elemental stiffness matrixes » 2020-11-03 12:03:25

#5 Re: Code_Aster usage » [SOLVED ]Get elements in group » 2020-10-30 12:03:12

#6 Re: Code_Aster usage » [SOLVED ]Get elements in group » 2020-10-30 10:43:09

#7 Re: Code_Aster usage » copy error » 2020-10-29 22:19:36

#8 Re: Salome-Meca usage » Apply different properties to each elements » 2020-10-29 16:52:07

#9 Re: Code_Aster usage » [SOLVED ]Get elements in group » 2020-10-29 16:49:50

#10 Code_Aster usage » COPY_ERROR when executing in python loop » 2020-10-29 16:46:35

#11 Code_Aster usage » [SOLVED ]Get elements in group » 2020-10-26 17:27:40

#12 Re: Code_Aster usage » [SOLVED] Code aster and python post processing » 2020-10-23 13:48:20

#13 Re: Code_Aster usage » element stiffness matrix extraction » 2020-10-22 14:08:22

#14 Re: Salome-Meca usage » Apply different properties to each elements » 2020-10-22 11:28:24

#15 Code_Aster usage » element stiffness matrix extraction » 2020-10-21 16:44:06

#16 Re: Code_Aster usage » [SOLVED] Code aster and python post processing » 2020-10-21 10:27:36

#17 Re: Code_Aster usage » [SOLVED] Code aster and python post processing » 2020-10-20 16:02:05

#18 Salome-Meca usage » Apply different properties to each elements » 2020-10-20 10:02:32

#19 Re: Code_Aster usage » [SOLVED] CALC_CHAMP and POST_CHAMP in COQUE » 2020-10-20 09:37:29

#20 Code_Aster usage » [SOLVED] CALC_CHAMP and POST_CHAMP in COQUE » 2020-10-14 12:40:38

#21 Re: Code_Aster usage » [SOLVED] Code aster and python post processing » 2020-10-07 11:59:18

#22 Re: Code_Aster usage » [SOLVED] Code aster and python post processing » 2020-10-06 10:20:37

#23 Code_Aster usage » [SOLVED] Code aster and python post processing » 2020-09-30 10:39:41

#24 Introduce yourself / Présentez vous » Greetings from ETH Zurich » 2020-09-30 08:42:15

Board footer